BmnSiProf.h

Go to the documentation of this file.

#ifndef BMNSiProf_H
#define BMNSiProf_H
 
#include "BmnSiProfPoint.h"
#include "FairDetector.h"
#include "TClonesArray.h"
#include "TGeoMedium.h"
#include "TLorentzVector.h"
#include "TVector3.h"
 
#include <map>
 
class BmnSiProf : public FairDetector
{
 
  public:
    // *@param name    detector name
    // *@param active  sensitivity flag
    BmnSiProf(const char* name, Bool_t active);
 
    BmnSiProf();
    virtual ~BmnSiProf();
 
    // Defines the action to be taken when a step is inside the
    // active volume. Creates BmnSiProfPoints and adds them to the collection.
    // @param vol Pointer to the active volume
    virtual Bool_t ProcessHits(FairVolume* vol = 0);
 
    // If verbosity level is set, print hit collection at the
    // end of the event and resets it afterwards.
    virtual void EndOfEvent();
 
    // Registers the hit collection in the ROOT manager.
    virtual void Register();
 
    // Accessor to the hit collection
    virtual TClonesArray* GetCollection(Int_t iColl) const;
 
    // Screen output of hit collection.
    virtual void Print(Option_t*) const;
 
    // Clears the hit collection
    virtual void Reset();
 
    // *@param cl1     Origin
    // *@param cl2     Target
    // *@param offset  Index offset
    virtual void CopyClones(TClonesArray* cl1, TClonesArray* cl2, Int_t offset);
 
    // Constructs the SiProf geometry
    virtual void ConstructGeometry();
 
    // Construct the geometry from a GDML geometry file
    virtual void ConstructGDMLGeometry(TGeoMatrix*);
 
    void ExpandNodeForGdml(TGeoNode* node);
 
    // Check whether a volume is sensitive.
    // The decision is based on the volume name. Only used in case
    // of GDML and ROOT geometry.
    // @param name    Volume name
    // @value         kTRUE if volume is sensitive, else kFALSE
    virtual Bool_t CheckIfSensitive(std::string name);
 
  private:
    // Track information to be stored until the track leaves the active volume.
    Int_t fTrackID;       //!  track index
    Int_t fVolumeID;      //!  volume id
    TVector3 fPosIn;      //!  position (in)
    TVector3 fPosOut;     //!  position (out)
    TVector3 fMomIn;      //!  momentum (in)
    TVector3 fMomOut;     //!  momentum (out)
    Double32_t fTime;     //!  time
    Double32_t fLength;   //!  length
    Double32_t fELoss;    //!  energy loss
    Int_t fIsPrimary;     //!  is track primary?
    Double_t fCharge;     //!  track charge
    Int_t fPdgId;         //!  pdg id of particle
 
    Int_t fPosIndex;                  //!
    TClonesArray* fPointCollection;   //! Hit collection
 
    std::map<TString, TGeoMedium*> fFixedMedia;   // List of media "repaired" after importing GMDL
 
    // Adds a BmnSiProfPoint to the HitCollection
    BmnSiProfPoint* AddHit(Int_t trackID,
                           Int_t detID,
                           TVector3 posIn,
                           TVector3 posOut,
                           TVector3 momIn,
                           TVector3 momOut,
                           Double_t time,
                           Double_t length,
                           Double_t eLoss,
                           Int_t isPrimary,
                           Double_t charge,
                           Int_t pdgId,
                           Int_t stationNum,
                           Int_t moduleNum);
 
    // Resets the private members for the track parameters
    void ResetParameters();
    BmnSiProf(const BmnSiProf&) = delete;
    BmnSiProf operator=(const BmnSiProf&) = delete;
 
    ClassDef(BmnSiProf, 1)
};
 
//------------------------------------------------------------------------------
inline void BmnSiProf::ResetParameters()
{
    fTrackID = -1;
    fVolumeID = 0;
    fPosIn.SetXYZ(0.0, 0.0, 0.0);
    fPosOut.SetXYZ(0.0, 0.0, 0.0);
    fMomIn.SetXYZ(0.0, 0.0, 0.0);
    fMomOut.SetXYZ(0.0, 0.0, 0.0);
    fTime = fLength = fELoss = 0.0;
    fPosIndex = 0;
};
//------------------------------------------------------------------------------
 
#endif /* BMNSiProf_H */