PartCorr split in 2 Base and Dep; coding violations corrected; PHOS geometry can...
[u/mrichter/AliRoot.git] / PWG4 / PartCorrBase / AliFidutialCut.h
1 #ifndef ALIFIDUTIALCUT_H
2 #define ALIFIDUTIALCUT_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4  * See cxx source for full Copyright notice     */
5 /* $Id:  $ */
6
7 //_________________________________________________________________________
8 // Class for track/cluster acceptance selection
9 // Selection in Central barrel, EMCAL and PHOS
10 //
11 //*-- Author: Gustavo Conesa (INFN-LNF)
12
13 // --- ROOT system ---
14 #include <TObject.h> 
15 #include <TArrayF.h> 
16
17 class TString ;
18 class TLorentzVector ;
19
20 //--- AliRoot system ---
21
22 class AliLog ;
23
24 class AliFidutialCut : public TObject {
25
26 public: 
27
28   AliFidutialCut() ; // ctor
29   AliFidutialCut(const AliFidutialCut & g) ; // cpy ctor
30   AliFidutialCut & operator = (const AliFidutialCut & g) ;//cpy assignment
31   virtual ~AliFidutialCut() ;//virtual dtor
32
33   void InitParameters();
34   
35   void Print(const Option_t * opt)const;
36   
37   Bool_t IsInFidutialCut(TLorentzVector l, TString det) const ;
38   
39   void DoCTSFidutialCut(Bool_t b) {fCTSFidutialCut = b; }
40   void DoEMCALFidutialCut(Bool_t b) {fEMCALFidutialCut = b; }
41   void DoPHOSFidutialCut(Bool_t b) {fPHOSFidutialCut = b; }
42
43   Bool_t GetCTSFidutialCutStatus() const {return fCTSFidutialCut ; }
44   Bool_t GetEMCALFidutialCut() const {return fEMCALFidutialCut ; }
45   Bool_t GetPHOSFidutialCutStatus() const {return fPHOSFidutialCut ; }
46
47   void SetSimpleCTSFidutialCut(const Float_t abseta, const Float_t phimin, const Float_t phimax) ;
48   void SetSimpleEMCALFidutialCut(const Float_t abseta, const Float_t phimin, const Float_t phimax) ;
49   void SetSimplePHOSFidutialCut(const Float_t abseta, const Float_t phimin, const Float_t phimax) ;
50
51   void AddCTSFidCutMaxEtaArray(TArrayF & array)  
52   {  fCTSFidCutMaxEta  = new TArrayF(array) ; } 
53   TArrayF * GetCTSFidCutMaxEtaArray() const   {return   fCTSFidCutMaxEta;}
54   
55   void AddCTSFidCutMaxPhiArray(TArrayF & array)  
56   {  fCTSFidCutMaxPhi  = new TArrayF(array) ; }
57   TArrayF * GetCTSFidCutMaxPhiArray() const   {return   fCTSFidCutMaxPhi;}
58   
59   void AddCTSFidCutMinEtaArray(TArrayF & array)  
60   {  fCTSFidCutMinEta  = new TArrayF(array) ; } 
61   TArrayF * GetCTSFidCutMinEtaArray() const   {return   fCTSFidCutMinEta;}
62  
63   void AddCTSFidCutMinPhiArray(TArrayF & array)  
64   {  fCTSFidCutMinPhi  = new TArrayF(array) ; }
65   TArrayF * GetCTSFidCutMinPhiArray() const   {return   fCTSFidCutMinPhi;}
66   
67   void AddEMCALFidCutMaxEtaArray(TArrayF & array)  
68   {  fEMCALFidCutMaxEta  = new TArrayF(array) ; } 
69   TArrayF * GetEMCALFidCutMaxEtaArray() const   {return   fEMCALFidCutMaxEta;}
70   
71   void AddEMCALFidCutMaxPhiArray(TArrayF & array)  
72   {  fEMCALFidCutMaxPhi  = new TArrayF(array) ; }
73   TArrayF * GetEMCALFidCutMaxPhiArray() const   {return   fEMCALFidCutMaxPhi;}
74  
75   void AddEMCALFidCutMinEtaArray(TArrayF & array)  
76   {  fEMCALFidCutMinEta  = new TArrayF(array) ; } 
77   TArrayF * GetEMCALFidCutMinEtaArray() const   {return   fEMCALFidCutMinEta;}
78   
79   void AddEMCALFidCutMinPhiArray(TArrayF & array)  
80   {  fEMCALFidCutMinPhi  = new TArrayF(array) ; }
81   TArrayF * GetEMCALFidCutMinPhiArray() const   {return   fEMCALFidCutMinPhi;}
82   
83   void AddPHOSFidCutMaxEtaArray(TArrayF & array)  
84   {  fPHOSFidCutMaxEta  = new TArrayF(array) ; } 
85   TArrayF * GetPHOSFidCutMaxEtaArray() const   {return   fPHOSFidCutMaxEta;}
86   
87    void AddPHOSFidCutMaxPhiArray(TArrayF & array)  
88    {  fPHOSFidCutMaxPhi  = new TArrayF(array) ; }
89    TArrayF * GetPHOSFidCutMaxPhiArray() const   {return   fPHOSFidCutMaxPhi;}
90    
91    void AddPHOSFidCutMinEtaArray(TArrayF & array)  
92    {  fPHOSFidCutMinEta  = new TArrayF(array) ; } 
93    TArrayF * GetPHOSFidCutMinEtaArray() const   {return   fPHOSFidCutMinEta;}
94  
95    void AddPHOSFidCutMinPhiArray(TArrayF & array)  
96    {  fPHOSFidCutMinPhi  = new TArrayF(array) ; }
97    TArrayF * GetPHOSFidCutMinPhiArray() const   {return   fPHOSFidCutMinPhi;}
98    
99  
100    
101  protected:
102    
103    //Detector acceptance cuts
104    Bool_t     fEMCALFidutialCut ; // Apply fidutial cuts to EMCAL clusters
105   Bool_t     fPHOSFidutialCut ;// Apply fidutial cuts to PHOS clusters
106   Bool_t     fCTSFidutialCut ;//Apply fidutial cuts to  CTS tracks
107   
108   TArrayF * fCTSFidCutMinEta ; //Take particles in CTS with eta > fCTSFidCutMinEta
109   TArrayF * fCTSFidCutMinPhi ; //Take particles in CTS with phi > fCTSFidCutMinPhi
110   TArrayF * fCTSFidCutMaxEta ; //Take particles in CTS with eta < fCTSFidCutMaxEta
111   TArrayF * fCTSFidCutMaxPhi ; //Take particles in CTS with phi > fCTSFidCutMaxPhi
112   
113   TArrayF * fEMCALFidCutMinEta ; //Take particles in EMCAL with eta > fEMCALFidCutMinEta
114   TArrayF * fEMCALFidCutMinPhi ; //Take particles in EMCAL with phi > fEMCALFidCutMinPhi
115   TArrayF * fEMCALFidCutMaxEta ; //Take particles in EMCAL with eta < fEMCALFidCutMaxEta
116   TArrayF * fEMCALFidCutMaxPhi ; //Take particles in EMCAL with phi > fEMCALFidCutMaxPhi
117   
118   TArrayF * fPHOSFidCutMinEta ; //Take particles in PHOS with eta > fPHOSFidCutMinEta
119   TArrayF * fPHOSFidCutMinPhi ; //Take particles in PHOS with phi > fPHOSFidCutMinPhi
120   TArrayF * fPHOSFidCutMaxEta ; //Take particles in PHOS with eta < fPHOSFidCutMaxEta
121   TArrayF * fPHOSFidCutMaxPhi ; //Take particles in PHOS with phi > fPHOSFidCutMaxPhi
122
123   
124   ClassDef(AliFidutialCut,1)
125 } ;
126
127
128 #endif //ALIFIDUTIALCUT_H
129
130
131