BmnDch.h

Go to the documentation of this file.

// -------------------------------------------------------------------------
// -----                         BmnDch header file                    -----
// -------------------------------------------------------------------------
 
#ifndef BMNDCH_H
#define BMNDCH_H
 
#include "BmnDchPoint.h"
#include "FairDetector.h"
#include "TClonesArray.h"
#include "TGeoMedium.h"
#include "TLorentzVector.h"
#include "TString.h"
#include "TVector3.h"
 
#include <map>
 
class BmnDch : public FairDetector
{
  public:
    // *@param name    detector name
    // *@param active  sensitivity flag
    BmnDch(const char* name, Bool_t active);
 
    BmnDch();
    virtual ~BmnDch();
 
    // Defines the action to be taken when a step is inside the
    // active volume. Creates BmnDchPoints 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 DCH geometry
    virtual void ConstructGeometry();
    // Construct the geometry from an ASCII geometry file
    virtual void ConstructAsciiGeometry();
    // Construct the geometry from a GDML geometry file
    virtual void ConstructGDMLGeometry(TGeoMatrix*);
    void ExpandNodeForGdml(TGeoNode* node);
    map<TString, TGeoMedium*> fFixedMedia;   // List of media "repaired" after importing GMDL
 
    // 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);
 
    // static Int_t GetWheel(Int_t uid){ return ((uid-1) >>3); };  //lsp  [0-1] == [inner,outer]
    // static Int_t GetProj(Int_t uid){ return ((uid-1) & 6)>>1; }; //lsp [0-3] == [x,y,u,v]
    // static Int_t GetGasGap(Int_t uid){ return ((uid-1) & 0x0001); }; //lsp [0-1] == [inner,outer]
 
  private:
    // Track information to be stored until the track leaves the active volume.
    Int_t fTrackID;        
    Int_t fVolumeID;       
    Int_t fPlaneNumber;    
    TVector3 fPos;         
    TVector3 fPosLocal;    
    TLorentzVector fMom;   
    Double32_t fTime;      
    Double32_t fLength;    
    Double32_t fELoss;     
    Int_t fIsPrimary;      
    Double_t fCharge;      
    Double_t fRadius;      
    Int_t fPdgId;          
    Int_t fwheel;
    TVector3 fPosIn;    
    TVector3 fPosOut;   
 
    TVector3 fPosInTmp;   
    Int_t fTrackIDTmp;    
    Int_t fwheelTmp;
 
    Int_t fPosIndex;                  
    TClonesArray* fPointCollection;   
 
    int DistAndPoints(TVector3 p3, TVector3 p4, TVector3& pa, TVector3& pb);
    TVector3 GlobalToLocal(TVector3& global);
    TVector3 LocalToGlobal(TVector3& local);
 
    // Adds a BmnDchPoint to the HitCollection
    BmnDchPoint* AddHit(Int_t trackID,
                        Int_t detID,
                        TVector3 pos,
                        Double_t radius,
                        TVector3 mom,
                        Double_t time,
                        Double_t length,
                        Double_t eLoss,
                        Int_t isPrimary,
                        Double_t charge,
                        Int_t fPdgId,
                        TVector3 trackPos);
 
    // Resets the private members for the track parameters
    void ResetParameters();
    BmnDch(const BmnDch&) = delete;
    BmnDch operator=(const BmnDch&) = delete;
 
    ClassDef(BmnDch, 1)
};
 
//------------------------------------------------------------------------------------------------------------------------
inline void BmnDch::ResetParameters()
{
    fTrackID = -1;
    fVolumeID = fwheel = 0;
    fPos.SetXYZ(0.0, 0.0, 0.0);
    fMom.SetXYZM(0.0, 0.0, 0.0, 0.0);
    fTime = fLength = fELoss = 0;
    fPosIndex = 0;
 
    fPosInTmp.SetXYZ(0.0, 0.0, 0.0);
    fTrackIDTmp = -1;
    fwheelTmp = -1;
};
 
#endif