1 #ifndef ALIFIDUCIALCUT_H
2 #define ALIFIDUCIALCUT_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
7 //_________________________________________________________________________
8 // Class for track/cluster acceptance selection
9 // Selection in Central barrel, EMCAL and PHOS
11 // Several selection regions possible for the different
14 //*-- Author: Gustavo Conesa (INFN-LNF)
16 // --- ROOT system ---
21 class TLorentzVector ;
23 //--- AliRoot system ---
26 class AliFiducialCut : public TObject {
30 AliFiducialCut() ; // ctor
31 AliFiducialCut(const AliFiducialCut & g) ; // cpy ctor
32 AliFiducialCut & operator = (const AliFiducialCut & g) ;//cpy assignment
33 virtual ~AliFiducialCut() ;//virtual dtor
35 void InitParameters();
37 void Print(const Option_t * opt)const;
38 Bool_t CheckFiducialRegion(const TLorentzVector lv, const TArrayF* minphi, const TArrayF* maxphi, const TArrayF* mineta, const TArrayF* maxeta) const ;
39 Bool_t IsInFiducialCut (const TLorentzVector lv, const TString det) const ;
41 void DoCTSFiducialCut(Bool_t b) {fCTSFiducialCut = b; }
42 void DoEMCALFiducialCut(Bool_t b) {fEMCALFiducialCut = b; }
43 void DoPHOSFiducialCut(Bool_t b) {fPHOSFiducialCut = b; }
45 Bool_t GetCTSFiducialCutStatus() const {return fCTSFiducialCut ; }
46 Bool_t GetEMCALFiducialCut() const {return fEMCALFiducialCut ; }
47 Bool_t GetPHOSFiducialCutStatus() const {return fPHOSFiducialCut ; }
49 void SetSimpleCTSFiducialCut(const Float_t abseta, const Float_t phimin, const Float_t phimax) ;
50 void SetSimpleEMCALFiducialCut(const Float_t abseta, const Float_t phimin, const Float_t phimax) ;
51 void SetSimplePHOSFiducialCut(const Float_t abseta, const Float_t phimin, const Float_t phimax) ;
53 void AddCTSFidCutMaxEtaArray(TArrayF & array)
54 { fCTSFidCutMaxEta = new TArrayF(array) ; }
55 TArrayF * GetCTSFidCutMaxEtaArray() const {return fCTSFidCutMaxEta;}
57 void AddCTSFidCutMaxPhiArray(TArrayF & array)
58 { fCTSFidCutMaxPhi = new TArrayF(array) ; }
59 TArrayF * GetCTSFidCutMaxPhiArray() const {return fCTSFidCutMaxPhi;}
61 void AddCTSFidCutMinEtaArray(TArrayF & array)
62 { fCTSFidCutMinEta = new TArrayF(array) ; }
63 TArrayF * GetCTSFidCutMinEtaArray() const {return fCTSFidCutMinEta;}
65 void AddCTSFidCutMinPhiArray(TArrayF & array)
66 { fCTSFidCutMinPhi = new TArrayF(array) ; }
67 TArrayF * GetCTSFidCutMinPhiArray() const {return fCTSFidCutMinPhi;}
69 void AddEMCALFidCutMaxEtaArray(TArrayF & array)
70 { fEMCALFidCutMaxEta = new TArrayF(array) ; }
71 TArrayF * GetEMCALFidCutMaxEtaArray() const {return fEMCALFidCutMaxEta;}
73 void AddEMCALFidCutMaxPhiArray(TArrayF & array)
74 { fEMCALFidCutMaxPhi = new TArrayF(array) ; }
75 TArrayF * GetEMCALFidCutMaxPhiArray() const {return fEMCALFidCutMaxPhi;}
77 void AddEMCALFidCutMinEtaArray(TArrayF & array)
78 { fEMCALFidCutMinEta = new TArrayF(array) ; }
79 TArrayF * GetEMCALFidCutMinEtaArray() const {return fEMCALFidCutMinEta;}
81 void AddEMCALFidCutMinPhiArray(TArrayF & array)
82 { fEMCALFidCutMinPhi = new TArrayF(array) ; }
83 TArrayF * GetEMCALFidCutMinPhiArray() const {return fEMCALFidCutMinPhi;}
85 void AddPHOSFidCutMaxEtaArray(TArrayF & array)
86 { fPHOSFidCutMaxEta = new TArrayF(array) ; }
87 TArrayF * GetPHOSFidCutMaxEtaArray() const {return fPHOSFidCutMaxEta;}
89 void AddPHOSFidCutMaxPhiArray(TArrayF & array)
90 { fPHOSFidCutMaxPhi = new TArrayF(array) ; }
91 TArrayF * GetPHOSFidCutMaxPhiArray() const {return fPHOSFidCutMaxPhi;}
93 void AddPHOSFidCutMinEtaArray(TArrayF & array)
94 { fPHOSFidCutMinEta = new TArrayF(array) ; }
95 TArrayF * GetPHOSFidCutMinEtaArray() const {return fPHOSFidCutMinEta;}
97 void AddPHOSFidCutMinPhiArray(TArrayF & array)
98 { fPHOSFidCutMinPhi = new TArrayF(array) ; }
99 TArrayF * GetPHOSFidCutMinPhiArray() const {return fPHOSFidCutMinPhi;}
105 //Detector acceptance cuts
106 Bool_t fEMCALFiducialCut ; // Apply fiducial cuts to EMCAL clusters
107 Bool_t fPHOSFiducialCut ;// Apply fiducial cuts to PHOS clusters
108 Bool_t fCTSFiducialCut ;//Apply fiducial cuts to CTS tracks
110 TArrayF * fCTSFidCutMinEta ; //Take particles in CTS with eta > fCTSFidCutMinEta
111 TArrayF * fCTSFidCutMinPhi ; //Take particles in CTS with phi > fCTSFidCutMinPhi
112 TArrayF * fCTSFidCutMaxEta ; //Take particles in CTS with eta < fCTSFidCutMaxEta
113 TArrayF * fCTSFidCutMaxPhi ; //Take particles in CTS with phi > fCTSFidCutMaxPhi
115 TArrayF * fEMCALFidCutMinEta ; //Take particles in EMCAL with eta > fEMCALFidCutMinEta
116 TArrayF * fEMCALFidCutMinPhi ; //Take particles in EMCAL with phi > fEMCALFidCutMinPhi
117 TArrayF * fEMCALFidCutMaxEta ; //Take particles in EMCAL with eta < fEMCALFidCutMaxEta
118 TArrayF * fEMCALFidCutMaxPhi ; //Take particles in EMCAL with phi > fEMCALFidCutMaxPhi
120 TArrayF * fPHOSFidCutMinEta ; //Take particles in PHOS with eta > fPHOSFidCutMinEta
121 TArrayF * fPHOSFidCutMinPhi ; //Take particles in PHOS with phi > fPHOSFidCutMinPhi
122 TArrayF * fPHOSFidCutMaxEta ; //Take particles in PHOS with eta < fPHOSFidCutMaxEta
123 TArrayF * fPHOSFidCutMaxPhi ; //Take particles in PHOS with phi > fPHOSFidCutMaxPhi
125 ClassDef(AliFiducialCut,1)
129 #endif //ALIFIDUCIALCUT_H