AliPhysics  5403132 (5403132)
AliConversionPhotonBase.h
Go to the documentation of this file.
1 #ifndef ALICONVERSIONPHOTONBASE_H
2 #define ALICONVERSIONPHOTONBASE_H
3 
4 #include "TMath.h"
5 #include "TParticle.h"
6 #include "AliLog.h"
7 #include "TObject.h"
8 #include "AliMCEvent.h"
9 #include "AliESDEvent.h"
10 #include "AliKFParticle.h"
11 #include "TParticle.h"
12 #include <vector>
13 #include "AliESDpid.h"
14 #include "TF1.h"
15 #include "TRandom3.h"
16 #include "AliPID.h"
17 #include "AliESDtrack.h"
18 #include "AliKFVertex.h"
19 #include "AliMCEventHandler.h"
20 #include "AliESDtrackCuts.h"
21 #include "AliGenCocktailEventHeader.h"
22 #include "TList.h"
23 
24 
25 using namespace std;
26 
28 
29  public:
30 
31  //Constructors
33 
34  //Copy Constructor
36  //assignment operator
37  AliConversionPhotonBase & operator = (const AliConversionPhotonBase & g);
38 
39  //Destructor
40  virtual ~AliConversionPhotonBase();
41 
43  void SetTag( Bool_t tagged ) { fTagged = tagged; }
44  Bool_t IsTagged(){return fTagged;}
45 
46  //Get the Chi2 of particle
47  void SetChi2perNDF(Float_t chi2) {fChi2perNDF = chi2;}
48  Float_t GetChi2perNDF() const {return fChi2perNDF;}
49 
50 
52  void SetLabelPositive(Int_t label){fLabel[0] = label;}
53  void SetLabelNegative(Int_t label){fLabel[1] = label;}
54  void SetTrackLabels(Int_t label1, Int_t label2){fLabel[0] = label1; fLabel[1] = label2;}
55  Int_t GetTrackLabelPositive() const{return fLabel[0];}
56  Int_t GetTrackLabelNegative() const {return fLabel[1];}
57  Int_t GetTrackLabel(Int_t i) const {return fLabel[i];}
58  virtual Int_t GetLabel(Int_t i) const { return GetTrackLabel(i); }
59  virtual Int_t GetLabel1() const { return GetTrackLabelPositive(); }
60  virtual Int_t GetLabel2() const { return GetTrackLabelNegative(); }
61 
62  // MC Label
63 
64  void SetMCLabel(Int_t* label){fMCLabel[0]=label[0];fMCLabel[1]=label[1];}
65  void SetMCLabelPositive(Int_t label){fMCLabel[0]=label;}
66  void SetMCLabelNegative(Int_t label){fMCLabel[1]=label;}
67  Int_t GetMCLabel(Int_t i) const{return fMCLabel[i];}
68  Int_t GetMCLabelPositive() const{return fMCLabel[0];}
69  Int_t GetMCLabelNegative() const{return fMCLabel[1];}
70  Int_t GetMCParticleLabel(AliMCEvent *mcEvent);
71 
72  // GetMCParticle
73 
74  Bool_t IsTruePhoton(AliMCEvent *mcEvent);
75  TParticle *GetMCParticle(AliMCEvent *mcEvent);
76  TParticle *GetPositiveMCDaughter(AliMCEvent *mcEvent){return GetMCDaughter(mcEvent,0);}
77  TParticle *GetNegativeMCDaughter(AliMCEvent *mcEvent){return GetMCDaughter(mcEvent,1);}
78  TParticle *GetMCDaughter(AliMCEvent *mcEvent,Int_t label);
79 
80  // V0Index
81  Int_t GetV0Index() const {return fV0Index;}
82  void SetV0Index(Int_t index) {fV0Index=index;}
83 
84  // leadingCellID - stored in fV0Index for GammaCalo analysis
85  Int_t GetLeadingCellID() const {return fV0Index;}
86  void SetLeadingCellID(Int_t index) {fV0Index=index;}
87 
88  // Conversion Point
89  void SetConversionPoint(Double_t convpoint[3]){fConversionPoint[0]=convpoint[0];fConversionPoint[1]=convpoint[1];fConversionPoint[2]=convpoint[2];}
90  void GetConversionPoint(Double_t convpoint[3]){convpoint[0]=fConversionPoint[0];convpoint[1]=fConversionPoint[1];convpoint[2]=fConversionPoint[2];}
91  Double_t GetConversionRadius() const {return TMath::Sqrt(fConversionPoint[0]*fConversionPoint[0]+fConversionPoint[1]*fConversionPoint[1]);}
92  Double_t GetConversionX() const {return fConversionPoint[0];}
93  Double_t GetConversionY() const {return fConversionPoint[1];}
94  Double_t GetConversionZ() const {return fConversionPoint[2];}
95  void GetDistanceOfClossetApproachToPrimVtx(const AliVVertex* primVertex, Float_t * dca);
96  void DeterminePhotonQuality(AliVTrack* negTrack, AliVTrack* posTrack);
97  UChar_t GetPhotonQuality() const {return fQuality;}
98  // Armenteros Qt Alpha
99  void GetArmenterosQtAlpha(Double_t qtalpha[2]){qtalpha[0]=fArmenteros[0];qtalpha[1]=fArmenteros[1];}
100  Double_t GetArmenterosQt() const {return fArmenteros[0];}
101  Double_t GetArmenterosAlpha() const {return fArmenteros[1];}
102 
103  // virtual functions to be implemented in KF/AOD classes
104 
105  virtual Double_t GetPhotonMass() const = 0;
106  virtual Double_t GetPhotonPt()const = 0;
107  virtual Double_t GetPhotonP() const = 0;
108  virtual Double_t GetPhotonEta() const = 0;
109  virtual Double_t GetPhotonPhi() const =0;
110 // virtual Double_t GetPhotonTheta() const =0;
111  virtual Double_t GetPx() const = 0;
112  virtual Double_t GetPy() const = 0;
113  virtual Double_t GetPz() const = 0;
114 
115 
116  Float_t GetMass() const { return fIMass; }
117  void SetMass( Float_t mass) { fIMass = mass; }
118 
119  Float_t GetPsiPair() const {return fPsiPair;}
120  void SetPsiPair(Float_t PsiPair){fPsiPair=PsiPair;}
121 
122 
123  protected:
124 
125  Int_t fLabel[2]; // Electron/Positron Track Label
126  Int_t fV0Index; // Index of the V0
127  Int_t fMCLabel[2]; // Electron/Positron MC Label
128  Float_t fChi2perNDF; // Chi2perNDF
129  Double_t fArmenteros[2]; // Armenteros Paramters
130  Double_t fConversionPoint[3]; // Conversion Point
131  Bool_t fTagged; // Is it tagged as decay pion (only for gammas)
132  Float_t fIMass; // Invariant Mass of dilepton pair
133  Float_t fPsiPair; // Psi Pair Value
134  UChar_t fQuality; //Photon Quality:
135  //0: garbage,
136  //1: both tracks TPC only,
137  //2: 1 track TPC only
138  //3: both tracks more than 1 ITS cluster
139  ClassDef(AliConversionPhotonBase,5);
140 };
141 
142 
143 #endif
144 
145 
146 
void GetArmenterosQtAlpha(Double_t qtalpha[2])
Double_t GetArmenterosAlpha() const
double Double_t
Definition: External.C:58
void SetConversionPoint(Double_t convpoint[3])
void SetMCLabelNegative(Int_t label)
Double_t mass
TParticle * GetPositiveMCDaughter(AliMCEvent *mcEvent)
virtual Int_t GetLabel(Int_t i) const
void GetConversionPoint(Double_t convpoint[3])
TParticle * GetNegativeMCDaughter(AliMCEvent *mcEvent)
int Int_t
Definition: External.C:63
Int_t GetMCLabel(Int_t i) const
float Float_t
Definition: External.C:68
void SetMCLabelPositive(Int_t label)
Int_t GetTrackLabel(Int_t i) const
void SetChi2perNDF(Float_t chi2)
void SetLeadingCellID(Int_t index)
virtual Int_t GetLabel2() const
void SetTrackLabels(Int_t label1, Int_t label2)
void SetLabelPositive(Int_t label)
Track labels.
void SetLabelNegative(Int_t label)
void SetTag(Bool_t tagged)
Set the tag for decay meson.
void SetPsiPair(Float_t PsiPair)
bool Bool_t
Definition: External.C:53
virtual Int_t GetLabel1() const
Double_t GetConversionRadius() const