BmnSiBT.h

Go to the documentation of this file.

#ifndef BMNSiBT_H
#define BMNSiBT_H
 
#include "BmnSiBTPoint.h"
#include "FairDetector.h"
#include "TClonesArray.h"
#include "TGeoMedium.h"
#include "TVector3.h"
 
#include <map>
 
class BmnSiBT : public FairDetector
{
  public:
    // *@param name    detector name
    // *@param active  sensitivity flag
    BmnSiBT(const char* name, Bool_t active);
 
    BmnSiBT();
    virtual ~BmnSiBT();
 
    // Defines the action to be taken when a step is inside the
    // active volume. Creates BmnSiBTPoints 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 SiBT 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 BmnSiBTPoint to the HitCollection
    BmnSiBTPoint* 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();
    BmnSiBT(const BmnSiBT&) = delete;
    BmnSiBT operator=(const BmnSiBT&) = delete;
 
    ClassDef(BmnSiBT, 1)
};
 
//------------------------------------------------------------------------------
inline void BmnSiBT::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 /* BMNSiBT_H */