]>
Commit | Line | Data |
---|---|---|
10eaad41 | 1 | #ifndef ALIANALYSISKINKESDMC_H |
2 | #define ALIANALYSISKINKESDMC_H | |
3 | ||
4 | /* See cxx source for full Copyright notice */ | |
5 | ||
6 | //----------------------------------------------------------------- | |
7 | // AliAnalysisKinkESDMC class | |
8 | // This task is an example of an analysis task | |
9 | // for kink topology Study | |
10 | // Authors: Martha Spyropoulou-Stassinaki | |
11 | // and members of the Greek group at the | |
12 | // Physics Department of Athens University | |
13 | // mspyrop@phys.uoa.gr | |
14 | //----------------------------------------------------------------- | |
15 | ||
a0271296 | 16 | class AliPIDResponse; |
17 | class AliESDEvent; | |
be1a7181 | 18 | class AliESDVertex; |
894840ad | 19 | class AliESDtrack; |
894840ad | 20 | class TF1; |
10eaad41 | 21 | class TH1F; |
22 | class TH2F; | |
23 | class TH1D; | |
894840ad | 24 | class TH2D; |
92adf4f6 | 25 | class TList; |
a0271296 | 26 | class AliESDtrackCuts; |
27 | class AliPhysicsSelection; | |
85f91829 | 28 | |
894840ad | 29 | #include "AliAnalysisTaskSE.h" |
30 | ||
10eaad41 | 31 | class AliAnalysisKinkESDMC : public AliAnalysisTaskSE { |
32 | public: | |
92adf4f6 | 33 | AliAnalysisKinkESDMC(const char *name = "AliAnalysisKinkESDMC"); |
10eaad41 | 34 | virtual ~AliAnalysisKinkESDMC() {} |
92adf4f6 | 35 | |
36 | virtual void UserCreateOutputObjects(); | |
37 | virtual void UserExec(Option_t *option); | |
10eaad41 | 38 | virtual void Terminate(Option_t *); |
92adf4f6 | 39 | |
10eaad41 | 40 | const AliESDVertex *GetEventVertex(const AliESDEvent* esd) const; |
a0271296 | 41 | // const AliESDVertex *GetEventVertex(AliESDEvent* esd) ; |
42 | void SetMulCut(Int_t low, Int_t up){fLowMulcut=low;fUpMulcut=up;} | |
43 | ||
44 | void SetKinkRadius(Float_t lRadiusKLow, Float_t lRadiusKUp) { fKinkRadLow=lRadiusKLow; fKinkRadUp=lRadiusKUp;} | |
45 | ||
10eaad41 | 46 | |
5220a87d | 47 | void SetNClusterCut(Int_t lowCluster){fLowCluster=lowCluster;} |
48 | void SetQtCut(Float_t lowQt){fLowQt=lowQt;} | |
2acbcadf | 49 | void SetYKRange(Float_t RapidityK){fRapiK=RapidityK;} |
5220a87d | 50 | |
10eaad41 | 51 | private: |
a0271296 | 52 | TH1F *fHistPtESD; //!Pt spectrum of all ESD inside eta, Pt cuts |
53 | TH1F *fHistPt; //!Pt spectrum of all ESD tracks | |
54 | TH1F *fHistQtAll; //!Qt spectrum of all kinks | |
55 | TH1F *fHistQt1; //!Qt spectrum of Kaon selected sample | |
56 | TH1F *fHistQt2; //!Qt spectrum in Qt region of kaons | |
57 | TH1F *fHistPtKaon; //!Pt Kaon spectrum of clean sample | |
58 | TH1F *fHistPtKPDG; //!Pt Kaon spectrum , confirmed by PDG,inside kaon Qt region | |
59 | TH1F *fHistEta; //!Eta spectrum of all kinks | |
60 | TH1F *fHistEtaK; //!Eta spectrum of kaons selected by kink topology | |
61 | TH1F *fptKMC; //!Pt Kaon spectrum MC, inside eta and pt cuts | |
62 | TH1F *fMultiplMC; //!charge multipl MC | |
63 | TH1F *fESDMult; //!ESD charged mult | |
64 | TH1F *frad; //!radius of kinks, MC , inside the eta nad Pt cuts | |
65 | TH1F *fradMC; //!radius of kinks, MC , inside the eta nad Pt cuts | |
66 | TH1F *fKinkKaon; //!Pt of PDG Kaons inside the selcted ones by the KInk topology | |
67 | TH1F *fKinkKaonBg; //!Pt of the BG inside the kink-Kaon identified spectrum | |
68 | TH1F *fM1kaon; //!inv mass of kink-tracks taken as kaons decaying to mu + neutrino | |
69 | TH1F *fgenPtEtR; //!MC Pt spectrum of kaons decaying to muon+neutrino and pi +pi, inside eta,Pt,Rad cuts | |
70 | TH1F *fPtKink; //!Pt spectrum of all kinks from track bank | |
71 | TH2F *fcodeH; //!PDG code(mother) vrs PDG dcode(daughter) of kinks with Qt <0.12 (fake) | |
72 | TH2F *fdcodeH; //!inks, code vrs dcode of BG,if mother code is 321 and daughter code > | |
73 | TH2F *fAngMomK; //! Decay angle vrs Mother Mom for pdg kaons | |
74 | TH2F *fAngMomPi; //! Decay angle vrs Mother Mom for pdg pions | |
75 | TH2F *fAngMomKC; //!Decay angle vrs Mother Mom for pdg kaons, inside the selected sample | |
76 | TH1F *fMultESDK; //!ESD charged mult | |
77 | TH1F *fMultMCK; //!MC K charged mult | |
78 | TH2F *fSignPtNcl;//!signPt vrs number of clusters in TPC for kaons from kink sele sample | |
79 | TH2F *fSignPtEta;//!signPt vrs Eta in TPC for kaons from kink sele sample | |
80 | TH2F *fSignPtEtaMC;//!signPt vrs Eta in TPC for kaons from kink sele sample | |
81 | TH1F *fSignPtMC;//!signPt in TPC for kaons | |
82 | TH2F *fEtaNcl;//!Eta vrs Nclu in TPC for kaons from kink sele sample | |
83 | TH1F *fSignPt;//!signPt in TPC for kaons from kink sele sample | |
84 | TH2F *fChi2NclTPC;//!chi2 vrs TPC Nclusters for kaons from kink sele sample | |
85 | TH1F *fRatChi2Ncl;//! Ratio chi2/ Ncl TPC for kaons from kink sele sample | |
86 | TH2F *fRadiusNcl;//! Radis f kink vetex for kaons from kink sele sample | |
87 | TH2F *fTPCSgnlP;//! TPC de/dx signal for kaons from kink sele sample | |
88 | TH2F *fTPCSgnlPa;//! TPC de/dx signal for kink sample | |
89 | TH1F *fSignPtGen;//!signPt in TPC for kaonsgenerated | |
90 | TH1D *fRpr;//! Radius of VTX at Y , X plane | |
91 | TH1D *fZpr;//!Z distrio of main vertex | |
92 | TH1D *fdcatoVxXY;//! dca to Vertex XY distr | |
93 | TH1F *fMCEtaKaon;//!MC eta for kaons | |
94 | TH2F *fZvXv;//! two dime of Z vrs X of vtx main | |
95 | TH2F *fZvYv;//! two dime of Z vrs Y of vtx main | |
96 | TH2F *fXvYv;//! two dime of X vrs Y of main tracks vtx main | |
97 | TH1F *fPtPrKink;//! pt of Primary PDG kaons inside the selected ones by the kink topology | |
98 | TH1F *fgenPtEtRP;//!MC Pt spectrum of kaons decaying to muon+neutrino and pi +pi, inside eta,Pt,Rad cuts | |
99 | TH1F *fgenPtEtRN;//!MC Pt spectrum of kaons decaying to muon+neutrino and pi +pi, inside eta,Pt,Rad cuts | |
100 | TH1F *fkinkKaonP;//!MC Pt spectrum of kaons decaying to muon+neutrino and pi +pi, inside eta,Pt,Rad cuts | |
101 | TH1F *fkinkKaonN;//!MC Pt spectrum of kaons decaying to muon+neutrino and pi +pi, inside eta,Pt,Rad cuts | |
102 | TH1F *frapidESDK;//!MC ESD K rapidity distr. | |
103 | TH1F *frapidKMC;//!MC MC K rapidity dis | |
104 | TH1F *fPtKPlMC;//!MC MC K rapidity dis | |
105 | TH1F *fPtKMnMC;//!MC MC K rapidity dis | |
106 | TH1F *fHistPtKaoP;//!MC MC K rapidity dis | |
107 | TH1F *fHistPtKaoN;//!MC MC K rapidity dis | |
108 | TH1F *fHiPtKPDGP;//!MC MC K rapidity dis | |
109 | TH1F *fHiPtKPDGN;//!MC MC K rapidity dis | |
110 | TH1F *fKinKBGP;//!MC MC K rapidity dis | |
111 | TH1F *fKinKBGN;//!MC MC K rapidity dis | |
112 | TH1F *fQtKMu;//!MC MC K Qt K to mu | |
113 | TH1F *fQtKPi;//!MC MC K Qt K to mu | |
114 | TH1F *fQtKEl;//!MC MC K Qt K to mu | |
115 | TH1F *fFakepipi;//!MC Fake pipi | |
116 | TH1F *fFakeKPi;//!MC Fake Kpi | |
117 | TH1F *fDCAkink;//!MC dcs kink | |
118 | TH1F *fDCAkinkBG;//!MC dcs kink | |
119 | TH2F *fPosiKink;//!MC position kink | |
120 | TH2F *fPosiKinkK;//!MC position kink | |
121 | TH2F *fPosiKinKXZ;//!MC position kink | |
122 | TH2F *fPosiKinKYZ;//!MC position kink | |
123 | TH2F *fPosiKinKBgZY;//!MC position kink | |
124 | TH2F *fcode2;//!PDG code(mother) vrs PDG dcode(daughter) of kinks with Qt <0.12 (fake) | |
125 | TH2F *fcode4;//!PDG code(mother) vrs PDG dcode(daughter) of kinks with Qt <0.12 (fake) | |
126 | TH2F *fZkinkZDau;//! z-position of kink z position of daughter bg ) | |
127 | TH1F *fQtKMuMC;//!MC MC K Qt K to mu | |
128 | TH1F *fQtKElMC;//!MC MC K Qt K to mu | |
129 | TH1F *fQtKPiMC;//!MC MC K Qt K to mu | |
130 | TH1F *fQtK3PiP;//!EDS K Qt K to 3Pi | |
131 | TH1F *fQtK3PiM;//!EDS K Qt K to 3Pi | |
132 | TH2F *fmaxAngMomKmu; //!Decay angle vrs Mother Mom for pdg kaons, inside the selected sample | |
133 | TH2F *fPosiKinKBgZX;//!MC position kink | |
134 | TH2F *fPosiKinKBgXY;//!MC position kink | |
135 | TH1F *fMinvPi;//!MC R life time | |
136 | TH1F *fMinvKa;//!MC R life time | |
137 | TH1F *fMinvPr;//!MC R life time | |
138 | TH2F *fTPCSgnlPtpc;//Kink mother moment vrs TPC signal | |
139 | TH2F *fTPCMomNSgnl;//kink mother TPC momentum vrs nsigmas of dEdx | |
140 | TH2F *fMothKinkMomSgnl;//kink mother TPC momentum vrs nsigmas of dEdx | |
141 | TH1F *fNSigmTPC;//kink mother TPC momentum vrs nsigmas of dEdx | |
142 | TH2F *fTPCSgnlKinkDau;//Kink mother moment vrs TPC signal | |
143 | TH2F *fcodeDau1;//!PDG code(mother) vrs PDG dcode(daughter) of kinks daughters | |
144 | TH2F *fcodeDau2;//!PDG code(mother) vrs PDG dcode(daughter) of kinks daughters , Bg study | |
145 | TH2F *fMothKinkMomSgnlD;//kink mother TPC momentum vrs nsigmas of dEdx | |
146 | TH1F *fInvMassMuNuAll;//kinks, Inv Mass all kinks MuNu | |
147 | TH2F *fInvMassMuNuPt;//kinks,Invariant Mass MuNu vs Pt | |
89db1464 | 148 | TH2F *fInvMassMuNuPtAll;//kinks,Invariant Mass MuNu vs Pt 22/8/2013 |
a0271296 | 149 | TH1F *fRatioCrossedRows; //ratio crossed rows |
150 | TH1F *fRatioCrossedRowsKink; //ratio crossed rows for kinks | |
151 | TH2F *fRadiusPt;//kinks, Radius vs Pt | |
152 | TH2F *fRadiusPtcln;//kinks, Radius vs Pt for clean kaons | |
5220a87d | 153 | TH2F *fRadiusPtKaon;//kinks, Radius vs Pt for clean kaons |
154 | TH2F *fRadiusPtPion;//kinks, Radius vs Pt for clean kaons | |
155 | TH2F *fRadiusPtFake;//kinks, Radius vs Pt for clean kaons | |
a0271296 | 156 | TH1F *fPtCut1; //K Pt spectrum of all kinks from track bank, K0 bins |
157 | TH1F *fPtCut2; //K Pt spectrum of all kinks from track bank, K0 bins | |
158 | TH1F *fPtCut3; //K Pt spectrum of all kinks from track bank, K0 bins | |
159 | TH2F *fAngMomKKinks;//kinks, Angle vs Momentum for K-kinks | |
160 | TH1F *flengthMCK;//!MC R life time | |
161 | TH1F *flifetiMCK;//!MC R life time | |
162 | TH1F *flifetim2;//!MC R life time | |
163 | TH1F *fLHelESDK;//!MC R life time | |
164 | TH1F *flifeInt;//!MC R life time | |
165 | TH1F *flifeYuri;//!MC R life time | |
166 | TH1F *flenYuri;//!MC R life time | |
167 | TH1F *flenTrRef;//!MC R life time | |
168 | TH1F *flifeSmall;//!MC R life time | |
169 | TH1F *flifetime;//!MC R life time | |
170 | TH1F *flifTiESDK;//!MC R life time | |
171 | TH1F *flifeKink;//!MC R life time | |
172 | TH1F *flenHelx;//!MC R life time | |
173 | TH3F *fradPtRapMC;//!MC R life time | |
174 | TH3F *fradPtRapDC;//!MC R life time | |
175 | TH3F *fradPtRapESD;//!MC R life time | |
176 | TH2F *fRadNclcln;//!MC R life time | |
177 | ||
178 | ||
179 | ||
180 | TF1 *f1; | |
181 | TF1 *f2; | |
182 | TList *fListOfHistos; //! list of histos | |
183 | ||
184 | Int_t fLowMulcut; // | |
185 | Int_t fUpMulcut; | |
60f5d945 | 186 | Float_t fKinkRadUp; |
187 | Float_t fKinkRadLow; | |
5220a87d | 188 | Int_t fLowCluster; |
189 | Float_t fLowQt; | |
2acbcadf | 190 | Float_t fRapiK; |
5220a87d | 191 | |
a0271296 | 192 | AliESDtrackCuts* fCutsMul; |
193 | ||
194 | AliESDtrackCuts* fMaxDCAtoVtxCut; | |
195 | AliPIDResponse *fPIDResponse; //! PID response object | |
10eaad41 | 196 | |
197 | AliAnalysisKinkESDMC(const AliAnalysisKinkESDMC&); // not implemented | |
198 | AliAnalysisKinkESDMC& operator=(const AliAnalysisKinkESDMC&); // not implemented | |
199 | ||
200 | ClassDef(AliAnalysisKinkESDMC, 1); // example of analysis | |
201 | }; | |
202 | ||
203 | #endif |