In This Package:

GtMuoneratorTool.h

Go to the documentation of this file.

#ifndef GTMUONERATORTOOL_H
#define GTMUONERATORTOOL_H

#include "GaudiKernel/RndmGenerators.h"
#include "GaudiAlg/GaudiTool.h"

#include "GenTools/IHepMCEventMutator.h"

#include "CLHEP/Vector/ThreeVector.h"

#include <string>
#include <vector>

class TH1F;

class GtMuoneratorTool : public GaudiTool,
                         virtual public IHepMCEventMutator
{
public: 
    GtMuoneratorTool(const std::string& type,
                     const std::string& name, 
                     const IInterface* parent);
    virtual ~GtMuoneratorTool();

    // GaudiTool interface
    virtual StatusCode initialize();
    virtual StatusCode finalize();

    // HepMCEventMutator interface
    virtual StatusCode mutate(HepMC::GenEvent& event);
    
private:

    // Hold properties
    std::string m_whichSite;
    std::string m_muonFileName;
    std::string m_ratioFileName;
    std::string m_volume;
    bool m_rotation;

    // just data
    TH1F* m_pmratio;
    std::vector<double> m_muonE, m_muonPhi, m_muonTheta;
    double m_Sx, m_Sy, m_Sz, m_Smax;
    CLHEP::Hep3Vector m_detectorDim;    //the water cherenkov detector dimension.
    CLHEP::Hep3Vector m_detectorCenter; // detector center.
    double m_dphi;

    IRndmGenSvc* m_rgs;
    Rndm::Numbers m_uni;

    StatusCode ReadMuons();
    double GetMuPlusMinusRatio(double momInGeV);
    int GetPIDFromMomentum(double momInGeV);
};


#endif  // GTMUONERATORTOOL_H

---