]>
Commit | Line | Data |
---|---|---|
1 | //Class to extract data to do ITS+TPC global Spectra | |
2 | //Autor Marek Chojnacki | |
3 | //emali Marek.Chojnacki@cern.ch | |
4 | //Used on 2009 data | |
5 | //last line of comments | |
6 | ||
7 | ||
8 | #include "AliAnalysisChargedHadronSpectraITSTruncatedMeanTask.h" | |
9 | #include "AliAnalysisManager.h" | |
10 | ||
11 | #include "AliESDEvent.h" | |
12 | ||
13 | ||
14 | //#include "AliESDtrack.h" | |
15 | ||
16 | #include "Riostream.h" | |
17 | #include "AliInputEventHandler.h" | |
18 | #include "AliStack.h" | |
19 | //#include "AliMCEventHandler.h" | |
20 | #include "AliMCEvent.h" | |
21 | ||
22 | #include "TGraph.h" | |
23 | #include "TH1F.h" | |
24 | #include "TH2F.h" | |
25 | #include "TMath.h" | |
26 | #include "TH3F.h" | |
27 | //#include "TMCProcess.h" | |
28 | #include "AliVEvent.h" | |
29 | ||
30 | #include "AliESDtrackCuts.h" | |
31 | //#include "AliESDpidCuts.h" | |
32 | //#include "AliESDpid.h" | |
33 | #include "AliCentrality.h" | |
34 | #include "AliESDUtils.h" | |
35 | #include "AliMultiplicity.h" | |
36 | ||
37 | class AliMCEventHandler; | |
38 | class Riostream; | |
39 | ||
40 | using namespace std; | |
41 | ClassImp(AliAnalysisChargedHadronSpectraITSTruncatedMeanTask) | |
42 | ||
43 | //________________________________________________________________________ | |
44 | AliAnalysisChargedHadronSpectraITSTruncatedMeanTask::AliAnalysisChargedHadronSpectraITSTruncatedMeanTask(const char *name) | |
45 | :AliAnalysisTaskSE(name),fESD(0),fCuts(0),fCutsMul(0),fMC(0), | |
46 | fLowMultiplicity(-1),fUpMultiplicity(-1),fLowCentrality(-10.0),fUpCentrality(-10.0),fSPD(0),fUsePilerejection(0), | |
47 | fYCut(100.0),fsigmacut(3.0),fnsigmaxy(7.0),fnsigmaz(5.0),fchargeCut(0.0), | |
48 | fCorrectSDD(0),fCorrectSSD(0),fHIsettings(0),fdovertexrescuts(0), | |
49 | fK0weight(0),flambdaweight(0),fAntilambdaweight(0), | |
50 | fHistStats(0),fHistZVertexBeforeCut(0),fHistZVertexAfterCut(0),fHistXYVertexBeforeCut(0),fHistXYVertexAfterCut(0), | |
51 | fHistPhiPtBeforeCuts(0),fHistPhiPtAfterCuts(0),fHistEtaPtBeforeCuts(0),fHistEtaPtAfterCuts(0),fHistDCABeforeCuts(0),fHistDCAAfterCuts(0), | |
52 | fHistPminusTPCinPAfterCuts(0),fHistPminusTPCinPglobalAfterCuts(0), | |
53 | fHistMydEPpositive(0),fHistMydETPCinPpositive(0),fHistMydETPCinPglobalpositive(0), | |
54 | fHistMydEPnegative(0),fHistMydETPCinPnegative(0),fHistMydETPCinPglobalnegative(0), | |
55 | fHistL3dEP(0),fHistL4dEP(0),fHistL5dEP(0),fHistL6dEP(0),fHistL3dETPCinP(0), | |
56 | fHistL4dETPCinP(0),fHistL5dETPCinP(0),fHistL6dETPCinP(0),fHistwhichhasmin(0),fHistMysignalminusESD(0), | |
57 | fHistminsignalifPionP(0),fHistminsignalifKaonP(0),fHistminsignalifProtonP(0),fHistminsignalifAntiPionP(0),fHistminsignalifAntiKaonP(0),fHistminsignalifAntiProtonP(0), | |
58 | fDCAXYZforcleanPions(0),fDCAXYZforcleanAntiPions(0),fDCAXYZforcleanProtons(0),fDCAXYZforcleanAntiProtons(0), | |
59 | fDCAXYZOpenforcleanPions(0),fDCAXYZOpenforcleanAntiPions(0),fDCAXYZOpenforcleanProtons(0),fDCAXYZOpenforcleanAntiProtons(0), | |
60 | fHistNtrackwithstandardcuts(0),fHistNtrackwithITSPIDcuts(0), | |
61 | fHistSignalinTPCKaonforstandardcuts(0),fHistSignalinTPCKaonforITSPIDcuts(0),fHistSignalinTPCAntiKaonforstandardcuts(0),fHistSignalinTPCAntiKaonforITSPIDcuts(0), | |
62 | fHistSignalinTPCProtonforstandardcuts(0),fHistSignalinTPCProtonforITSPIDcuts(0),fHistSignalinTPCAntiProtonforstandardcuts(0),fHistSignalinTPCAntiProtonforITSPIDcuts(0), | |
63 | fHistStandartMul(0),fHistMytrackMul(0),fHistStandartMulvSPD2(0), | |
64 | fHistminsignalifPionPPrimary(0),fHistminsignalifKaonPPrimary(0),fHistminsignalifProtonPPrimary(0),fHistminsignalifProtonPPrimaryfake(0), | |
65 | fHistminsignalifAntiPionPPrimary(0),fHistminsignalifAntiKaonPPrimary(0),fHistminsignalifAntiProtonPPrimary(0),fHistminsignalifAntiProtonPPrimaryfake(0), | |
66 | fHistminsignalifPionPSecondary(0),fHistminsignalifKaonPSecondary(0), | |
67 | fHistminsignalifProtonPSecondaryWD(0),fHistminsignalifProtonPSecondaryHI(0),fHistminsignalifProtonPSecondaryRest(0), | |
68 | fHistminsignalifProtonPSecondaryWDfake(0),fHistminsignalifProtonPSecondaryHIfake(0), | |
69 | fHistminsignalifAntiPionPSecondary(0),fHistminsignalifAntiKaonPSecondary(0), | |
70 | fHistminsignalifAntiProtonPSecondaryWD(0),fHistminsignalifAntiProtonPSecondaryHI(0), fHistminsignalifAntiProtonPSecondaryRest(0), | |
71 | fHistminsignalifAntiProtonPSecondaryWDfake(0),fHistminsignalifAntiProtonPSecondaryHIfake(0), | |
72 | fHistminsignalifMuEPositiveP(0),fHistminsignalifMuENegativeP(0), | |
73 | fHistminsignalifPionPrimaryfake(0),fHistminsignalifKaonPrimaryfake(0),fHistminsignalifAntiPionPrimaryfake(0),fHistminsignalifAntiKaonPrimaryfake(0), | |
74 | fHistminsignalifPionSecondaryfake(0),fHistminsignalifKaonSecondaryfake(0),fHistminsignalifAntiPionSecondaryfake(0),fHistminsignalifAntiKaonSecondaryfake(0), | |
75 | fHistminsignalifPionPMCPrimary(0),fHistminsignalifKaonPMCPrimary(0),fHistminsignalifProtonPMCPrimary(0), | |
76 | fHistminsignalifAntiPionPMCPrimary(0),fHistminsignalifAntiKaonPMCPrimary(0),fHistminsignalifAntiProtonPMCPrimary(0), | |
77 | fHistminsignalifPionPMCPrimaryBeforeEventCuts(0),fHistminsignalifKaonPMCPrimaryBeforeEventCuts(0),fHistminsignalifProtonPMCPrimaryBeforeEventCuts(0), | |
78 | fHistminsignalifAntiPionPMCPrimaryBeforeEventCuts(0),fHistminsignalifAntiKaonPMCPrimaryBeforeEventCuts(0),fHistminsignalifAntiProtonPMCPrimaryBeforeEventCuts(0), | |
79 | fHistminsignalifPionPMCPrimaryBeforeEventCutswithgoodZvertex(0),fHistminsignalifKaonPMCPrimaryBeforeEventCutswithgoodZvertex(0),fHistminsignalifProtonPMCPrimaryBeforeEventCutswithgoodZvertex(0), | |
80 | fHistminsignalifAntiPionPMCPrimaryBeforeEventCutswithgoodZvertex(0),fHistminsignalifAntiKaonPMCPrimaryBeforeEventCutswithgoodZvertex(0),fHistminsignalifAntiProtonPMCPrimaryBeforeEventCutswithgoodZvertex(0), | |
81 | fHistminsignalifPionPMCPrimaryAfterEventCutsBeforeVertexZ(0),fHistminsignalifKaonPMCPrimaryAfterEventCutsBeforeVertexZ(0),fHistminsignalifProtonPMCPrimaryAfterEventCutsBeforeVertexZ(0), | |
82 | fHistminsignalifAntiPionPMCPrimaryAfterEventCutsBeforeVertexZ(0),fHistminsignalifAntiKaonPMCPrimaryAfterEventCutsBeforeVertexZ(0),fHistminsignalifAntiProtonPMCPrimaryAfterEventCutsBeforeVertexZ(0), | |
83 | fDCAXYZforcleanPionsMCPrimary(0),fDCAXYZforcleanAntiPionsMCPrimary(0),fDCAXYZforcleanProtonsMCPrimary(0),fDCAXYZforcleanAntiProtonsMCPrimary(0), | |
84 | fDCAXYZforcleanPionsWD(0),fDCAXYZforcleanAntiPionsWD(0),fDCAXYZforcleanProtonsWD(0), fDCAXYZforcleanAntiProtonsWD(0),fDCAXYZforcleanPionsHI(0),fDCAXYZforcleanAntiPionsHI(0), | |
85 | fDCAXYZforcleanProtonsHI(0),fDCAXYZforcleanAntiProtonsHI(0),fDCAXYZforcleanPionsMEPrimary(0),fDCAXYZforcleanAntiPionsMEPrimary(0),fDCAXYZforcleanPionsMESecondary(0),fDCAXYZforcleanAntiPionsMESecondary(0),fDCAXYZforcleanPionsR(0),fDCAXYZforcleanAntiPionsR(0),fDCAXYZforcleanProtonsR(0),fDCAXYZforcleanAntiProtonsR(0), | |
86 | fDCAXYZOpenforcleanPionsMCPrimary(0),fDCAXYZOpenforcleanAntiPionsMCPrimary(0),fDCAXYZOpenforcleanProtonsMCPrimary(0),fDCAXYZOpenforcleanAntiProtonsMCPrimary(0), | |
87 | fDCAXYZOpenforcleanPionsWD(0),fDCAXYZOpenforcleanAntiPionsWD(0),fDCAXYZOpenforcleanProtonsWD(0), fDCAXYZOpenforcleanAntiProtonsWD(0),fDCAXYZOpenforcleanPionsHI(0),fDCAXYZOpenforcleanAntiPionsHI(0), | |
88 | fDCAXYZOpenforcleanProtonsHI(0),fDCAXYZOpenforcleanAntiProtonsHI(0),fDCAXYZOpenforcleanPionsMEPrimary(0),fDCAXYZOpenforcleanAntiPionsMEPrimary(0),fDCAXYZOpenforcleanPionsMESecondary(0),fDCAXYZOpenforcleanAntiPionsMESecondary(0),fDCAXYZOpenforcleanPionsR(0),fDCAXYZOpenforcleanAntiPionsR(0),fDCAXYZOpenforcleanProtonsR(0),fDCAXYZOpenforcleanAntiProtonsR(0), | |
89 | fElectronsource(0),fAntiElectronsource(0), | |
90 | fMuonsource(0),fAntiMuonsource(0), | |
91 | fPionNTPCClusters(0),fAntiPionNTPCClusters(0),fKaonNTPCClusters(0),fAntiKaonNTPCClusters(0),fProtonNTPCClusters(0),fAntiProtonNTPCClusters(0), | |
92 | fPionchi2(0),fAntiPionchi2(0),fKaonchi2(0),fAntiKaonchi2(0),fProtonchi2(0),fAntiProtonchi2(0), | |
93 | fTracksCutmonitoring(0),fParticlesCutmonitoring(0),fVertexshift(0),fPtESDminusPtMCvPtESDafterallcuts(0),fPtESDminusPtMCvPtESDafterTPCcuts(0),fMulESDMulMCVz(0), | |
94 | fTPCPIDCUT(0), fESDpid(0),fPrimaryElectronsMother(0), | |
95 | flist(0) | |
96 | { | |
97 | //Constructor | |
98 | fESDpid=new AliESDpid(); | |
99 | // fESDpid->GetTPCResponse().SetBetheBlochParameters(0.0283086,2.63394e+01,5.04114e-11,2.12543e+00,4.88663e+00); | |
100 | fESDpid->GetTPCResponse().SetBetheBlochParameters(1.28949/50., 2.74095e+01, TMath::Exp(-3.21763e+01), 2.44026, 6.58800); | |
101 | ||
102 | ||
103 | ||
104 | fCutsMul=new AliESDtrackCuts("Mul","Mul"); | |
105 | fCutsMul->SetMinNClustersTPC(70); | |
106 | fCutsMul->SetMaxChi2PerClusterTPC(4); | |
107 | fCutsMul->SetAcceptKinkDaughters(kFALSE); | |
108 | fCutsMul->SetRequireTPCRefit(kTRUE); | |
109 | // ITS | |
110 | fCutsMul->SetRequireITSRefit(kTRUE); | |
111 | fCutsMul->SetClusterRequirementITS(AliESDtrackCuts::kSPD, | |
112 | AliESDtrackCuts::kAny); | |
113 | fCutsMul->SetMaxDCAToVertexXYPtDep("0.0182+0.0350/pt^1.01"); | |
114 | ||
115 | fCutsMul->SetMaxDCAToVertexZ(2); | |
116 | fCutsMul->SetDCAToVertex2D(kFALSE); | |
117 | fCutsMul->SetRequireSigmaToVertex(kFALSE); | |
118 | ||
119 | fCutsMul->SetEtaRange(-0.8,+0.8); | |
120 | fCutsMul->SetPtRange(0.15, 1e10); | |
121 | ||
122 | fdcaxypar[0]=0.0050; | |
123 | fdcaxypar[1]=0.0060; | |
124 | fdcaxypar[2]=0.9; | |
125 | ||
126 | fdcazpar[0]=0.0146; | |
127 | fdcazpar[1]=0.0070; | |
128 | fdcazpar[2]=1.114758; | |
129 | fdcazpar[3]=0.0216; | |
130 | ||
131 | ||
132 | ||
133 | ||
134 | Printf("end of AliAnalysisChargedHadronSpectraITSTruncatedMeanTask"); | |
135 | DefineOutput(1, TList::Class()); | |
136 | } | |
137 | ||
138 | ||
139 | //________________________________________________________________________ | |
140 | void AliAnalysisChargedHadronSpectraITSTruncatedMeanTask::UserCreateOutputObjects() | |
141 | { | |
142 | //UserCreateOutputObject | |
143 | Printf("AliAnalysisChargedHadronSpectraITSTruncatedMeanTask UserCreateOutputObjects"); | |
144 | flist=new TList(); | |
145 | flist->SetOwner(); | |
146 | Float_t ptmax=2.0; | |
147 | Float_t etamax=1.0; | |
148 | ||
149 | Int_t netabins=100; | |
150 | Int_t nptbins=40; | |
151 | Double_t dcamax=3.7; | |
152 | const Int_t ndec=2; | |
153 | Int_t startvalue=-1; | |
154 | const Int_t npredec=50; | |
155 | Double_t tabx[ndec*npredec+1]; | |
156 | for (Int_t i=0;i<ndec;i++) | |
157 | { | |
158 | for (Int_t j=0;j<npredec;j++) | |
159 | { | |
160 | tabx[npredec*i+j]=TMath::Power(10,((Double_t)i)+((Double_t)startvalue)+((Double_t)j)/((Double_t)npredec)); | |
161 | } | |
162 | } | |
163 | tabx[ndec*npredec]=TMath::Power(10,ndec+startvalue); | |
164 | ||
165 | const Int_t ny=600; | |
166 | const Double_t jump=1.5; | |
167 | const Double_t starty=0.0; | |
168 | ||
169 | ||
170 | Int_t kPtBins=30; | |
171 | Double_t binsPtDummy[kPtBins+1]; | |
172 | binsPtDummy[0]=0.0; | |
173 | for(int i=1;i<=kPtBins+1;i++) | |
174 | { | |
175 | if(binsPtDummy[i-1]+0.05<1.01) | |
176 | binsPtDummy[i]=binsPtDummy[i-1]+0.05; | |
177 | else | |
178 | binsPtDummy[i]=binsPtDummy[i-1]+0.1; | |
179 | } | |
180 | //{0., 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0}; | |
181 | ||
182 | ||
183 | fHistStats=new TH1F("StatsHist","StatsHist",10,-0.5,9.5); | |
184 | fHistStats->GetXaxis()->SetBinLabel(1,"Nevents"); | |
185 | fHistStats->GetXaxis()->SetBinLabel(2,"Nevents_physics"); | |
186 | fHistStats->GetXaxis()->SetBinLabel(3,"Nevents_physics_with_good_SPD_vertex"); | |
187 | fHistStats->GetXaxis()->SetBinLabel(4,"Nevents_physics_with_good_SPD_vertex_and_|z|<10.0"); | |
188 | fHistStats->GetXaxis()->SetBinLabel(5,"N_tracks_with_3_or_4_in_SSD_SDD"); | |
189 | fHistStats->GetXaxis()->SetBinLabel(6,"N_tracks_with_3_or_4_with_goodq_in_SSD_SDD"); | |
190 | fHistStats->GetXaxis()->SetBinLabel(7,"e_in_pi"); | |
191 | fHistStats->GetXaxis()->SetBinLabel(8,"mu_in_pi"); | |
192 | fHistStats->GetXaxis()->SetBinLabel(9,"MC_event"); | |
193 | fHistStats->GetXaxis()->SetBinLabel(10,"MC_event_with_z<10.0"); | |
194 | flist->Add(fHistStats); | |
195 | ||
196 | fHistZVertexBeforeCut=new TH1F("HistZVertexBeforeCut","ZVertex;z[cm];N_{counts}",400,-20,20); | |
197 | flist->Add(fHistZVertexBeforeCut); | |
198 | fHistZVertexAfterCut=new TH1F("HistZVertexAfterCut","ZVertex;z[cm];N_{counts}",400,-20,20); | |
199 | flist->Add(fHistZVertexAfterCut); | |
200 | fHistXYVertexBeforeCut=new TH2F("HistXYVertexBeforeCut","XYVertex;x[cm];y[cm];N_{conuts}",100,-0.4,0.4,100,-0.4,0.4); | |
201 | flist->Add(fHistXYVertexBeforeCut); | |
202 | fHistXYVertexAfterCut=new TH2F("HistXYVertexAfterCut","XYVertex;x[cm];y[cm];N_{conuts}",100,-0.4,0.4,100,-0.4,0.4); | |
203 | flist->Add(fHistXYVertexAfterCut); | |
204 | ||
205 | fHistPhiPtBeforeCuts=new TH2F("HistPhiPtBeforeCuts",";#phi;pt[GeV/c]",70,0,2.0*TMath::Pi(),nptbins,-1.0*ptmax,ptmax); | |
206 | flist->Add(fHistPhiPtBeforeCuts); | |
207 | ||
208 | fHistPhiPtAfterCuts=new TH2F("HistPhiPtAfterCuts",";#phi;pt[GeV/c]",70,0,2.0*TMath::Pi(),nptbins,-1.0*ptmax,ptmax); | |
209 | flist->Add(fHistPhiPtAfterCuts); | |
210 | ||
211 | fHistEtaPtBeforeCuts=new TH2F("HistEtaPtBeforeCuts",";#eta;pt[GeV/c]",netabins,-1.0*etamax,etamax,nptbins,-1.0*ptmax,ptmax);//eta pt | |
212 | flist->Add(fHistEtaPtBeforeCuts); | |
213 | ||
214 | fHistEtaPtAfterCuts=new TH2F("HistEtaPtAfterCuts",";#eta;pt[GeV/c]",netabins,-1.0*etamax,etamax,nptbins,-1.0*ptmax,ptmax);//eta pt | |
215 | flist->Add(fHistEtaPtAfterCuts); | |
216 | ||
217 | fHistDCABeforeCuts=new TH2F("HistDCABeforeCut",";dcaxy[cm];dcaz[cm]",100,-1.0*dcamax,dcamax,100,-1.0*dcamax,1.0*dcamax); | |
218 | flist->Add(fHistDCABeforeCuts); | |
219 | ||
220 | fHistDCAAfterCuts=new TH2F("HistDCAAfterCut",";dcaxy[cm];dcaz[cm]",100,-1.0*dcamax,dcamax,100,-1.0*dcamax,1.0*dcamax); | |
221 | flist->Add(fHistDCAAfterCuts); | |
222 | ||
223 | ||
224 | fHistPminusTPCinPAfterCuts= new TH2F("HistPminusTPCinPVPTPCinAfterCuts",";P-PTPCin [GeV/c];PTPCin",100,-0.5,0.5,ndec*npredec,tabx); | |
225 | flist->Add(fHistPminusTPCinPAfterCuts); | |
226 | ||
227 | fHistPminusTPCinPglobalAfterCuts= new TH2F("HistPminusTPCinPVPTPCinglobalAfterCuts",";P-PTPCinglobal [GeV/c];PTPCin",100,-0.5,0.5,ndec*npredec,tabx); | |
228 | flist->Add(fHistPminusTPCinPglobalAfterCuts); | |
229 | ||
230 | fHistMydEPpositive=new TH2F("HistMydEPpositive",";P[GeV/c];dE[in 300#mum]",ndec*npredec,tabx,ny,starty,ny*jump+starty); | |
231 | flist->Add(fHistMydEPpositive); | |
232 | ||
233 | fHistMydETPCinPpositive=new TH2F("HistMydETPCinPpositive",";P[GeV/c];dE[in 300#mum]",ndec*npredec,tabx,ny,starty,ny*jump+starty); | |
234 | flist->Add(fHistMydETPCinPpositive); | |
235 | ||
236 | fHistMydETPCinPglobalpositive=new TH2F("HistMydETPCinPglobalpositive",";P[GeV/c];dE[in 300#mum]",ndec*npredec,tabx,ny,starty,ny*jump+starty); | |
237 | flist->Add(fHistMydETPCinPglobalpositive); | |
238 | ||
239 | fHistMydEPnegative=new TH2F("HistMydEPnegative",";P[GeV/c];dE[in 300#mum]",ndec*npredec,tabx,ny,starty,ny*jump+starty); | |
240 | flist->Add(fHistMydEPnegative); | |
241 | ||
242 | fHistMydETPCinPnegative=new TH2F("HistMydETPCinPnegative",";P[GeV/c];dE[in 300#mum]",ndec*npredec,tabx,ny,starty,ny*jump+starty); | |
243 | flist->Add(fHistMydETPCinPnegative); | |
244 | ||
245 | fHistMydETPCinPglobalnegative=new TH2F("HistMydETPCinPglobalnegative",";P[GeV/c];dE[in 300#mum]",ndec*npredec,tabx,ny,starty,ny*jump+starty); | |
246 | flist->Add(fHistMydETPCinPglobalnegative); | |
247 | ||
248 | ||
249 | fHistL3dEP=new TH2F("HistL3dEP",";P[GeV/c];dE[in 300#mum]",ndec*npredec,tabx,ny,starty,ny*jump+starty); | |
250 | flist->Add(fHistL3dEP); | |
251 | ||
252 | fHistL4dEP=new TH2F("HistL4dEP",";P[GeV/c];dE[in 300#mum]",ndec*npredec,tabx,ny,starty,ny*jump+starty); | |
253 | flist->Add(fHistL4dEP); | |
254 | ||
255 | fHistL5dEP=new TH2F("HistL5dEP",";P[GeV/c];dE[in 300#mum]",ndec*npredec,tabx,ny,starty,ny*jump+starty); | |
256 | flist->Add(fHistL5dEP); | |
257 | ||
258 | fHistL6dEP=new TH2F("HistL6dEP",";P[GeV/c];dE[in 300#mum]",ndec*npredec,tabx,ny,starty,ny*jump+starty); | |
259 | flist->Add(fHistL6dEP); | |
260 | ||
261 | fHistL3dETPCinP=new TH2F("HistL3dETPCinP",";P[GeV/c];dE[in 300#mum]",ndec*npredec,tabx,ny,starty,ny*jump+starty); | |
262 | flist->Add(fHistL3dETPCinP); | |
263 | ||
264 | fHistL4dETPCinP=new TH2F("HistL4dETPCinP",";P[GeV/c];dE[in 300#mum]",ndec*npredec,tabx,ny,starty,ny*jump+starty); | |
265 | flist->Add(fHistL4dETPCinP); | |
266 | ||
267 | fHistL5dETPCinP=new TH2F("HistL5dETPCinP",";P[GeV/c];dE[in 300#mum]",ndec*npredec,tabx,ny,starty,ny*jump+starty); | |
268 | flist->Add(fHistL5dETPCinP); | |
269 | ||
270 | fHistL6dETPCinP=new TH2F("HistL6dETPCinP",";P[GeV/c];dE[in 300#mum]",ndec*npredec,tabx,ny,starty,ny*jump+starty); | |
271 | flist->Add(fHistL6dETPCinP); | |
272 | ||
273 | ||
274 | ||
275 | fHistwhichhasmin=new TH2F("Histwhichhasmin","Histwhichhasmin;L;Q",4,-0.5,3.5,100,0,1000); | |
276 | fHistwhichhasmin->GetXaxis()->SetBinLabel(1,"SDD1"); | |
277 | fHistwhichhasmin->GetXaxis()->SetBinLabel(2,"SDD2"); | |
278 | fHistwhichhasmin->GetXaxis()->SetBinLabel(3,"SSD1"); | |
279 | fHistwhichhasmin->GetXaxis()->SetBinLabel(4,"SSD2"); | |
280 | flist->Add(fHistwhichhasmin); | |
281 | fHistMysignalminusESD=new TH1F("HistMysignalminus","HistMysignalminus;my-ESD;N",100,-0.2,0.2); | |
282 | flist->Add(fHistMysignalminusESD); | |
283 | ||
284 | ||
285 | fHistminsignalifPionP=new TH2F("HistminsignalifPionP",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
286 | flist->Add(fHistminsignalifPionP); | |
287 | fHistminsignalifKaonP=new TH2F("HistminsignalifKaonP",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
288 | flist->Add(fHistminsignalifKaonP); | |
289 | fHistminsignalifProtonP=new TH2F("HistminsignalifProtonP",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
290 | flist->Add(fHistminsignalifProtonP); | |
291 | ||
292 | ||
293 | fHistminsignalifAntiPionP=new TH2F("HistminsignalifAntiPionP",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
294 | flist->Add(fHistminsignalifAntiPionP); | |
295 | fHistminsignalifAntiKaonP=new TH2F("HistminsignalifAntiKaonP",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
296 | flist->Add(fHistminsignalifAntiKaonP); | |
297 | fHistminsignalifAntiProtonP=new TH2F("HistminsignalifAntiProtonP",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
298 | flist->Add(fHistminsignalifAntiProtonP); | |
299 | ||
300 | Int_t kDCABins=20+2+40; | |
301 | Double_t binsDCADummy[62+1]={-3.0,-2.7,-2.4,-2.1,-1.8,-1.5,-1.2,-0.9,-0.6,-0.3,-0.25,-0.2,-0.19,-0.18,-0.17,-0.16,-0.15,-0.14,-0.13,-0.12,-0.11,-0.10,-0.09,-0.08,-0.07,-0.06,-0.05,-0.04,-0.03,-0.02,-0.01,0.0,0.01,0.02,0.03,0.04,0.05,0.06,0.07,0.08,0.09,0.1,0.11,0.12,0.13,0.14,0.15,0.16,0.17,0.18,0.19,0.2,0.25,0.3,0.6,0.9,1.2,1.5,1.8,2.1,2.4,2.7,3.0}; | |
302 | ||
303 | ||
304 | ||
305 | ||
306 | fDCAXYZforcleanPions=new TH3F("fDCAXYZforcleanPions",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
307 | flist->Add(fDCAXYZforcleanPions); | |
308 | fDCAXYZforcleanAntiPions=new TH3F("fDCAXYZforcleanAntiPions",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
309 | flist->Add(fDCAXYZforcleanAntiPions); | |
310 | fDCAXYZforcleanProtons=new TH3F("fDCAXYZforcleanProtons",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
311 | flist->Add(fDCAXYZforcleanProtons); | |
312 | fDCAXYZforcleanAntiProtons=new TH3F("fDCAXYZforcleanAntiProtons",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
313 | flist->Add(fDCAXYZforcleanAntiProtons); | |
314 | ||
315 | ||
316 | fDCAXYZOpenforcleanPions=new TH3F("fDCAXYZOpenforcleanPions",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
317 | flist->Add(fDCAXYZOpenforcleanPions); | |
318 | fDCAXYZOpenforcleanAntiPions=new TH3F("fDCAXYZOpenforcleanAntiPions",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
319 | flist->Add(fDCAXYZOpenforcleanAntiPions); | |
320 | fDCAXYZOpenforcleanProtons=new TH3F("fDCAXYZOpenforcleanProtons",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
321 | flist->Add(fDCAXYZOpenforcleanProtons); | |
322 | fDCAXYZOpenforcleanAntiProtons=new TH3F("fDCAXYZOpenforcleanAntiProtons",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
323 | flist->Add(fDCAXYZOpenforcleanAntiProtons); | |
324 | ||
325 | fHistNtrackwithstandardcuts=new TH2F ("fHistNtrackwithstandardcuts",";Pt[GeV/c];type;counts",kPtBins,binsPtDummy,3,-0.5,2.5); | |
326 | flist->Add(fHistNtrackwithstandardcuts); | |
327 | fHistNtrackwithITSPIDcuts=new TH2F ("fHistNtrackwithITSPIDcuts",";Pt[GeV/c];type;counts",kPtBins,binsPtDummy,3,-0.5,2.5); | |
328 | flist->Add(fHistNtrackwithITSPIDcuts); | |
329 | ||
330 | ||
331 | fHistSignalinTPCKaonforstandardcuts= new TH2F("fHistSignalinTPCKaonforstandardcuts",";Pt[GeV/c];signal",kPtBins,binsPtDummy,100,-1.0,1.0) ; | |
332 | flist->Add(fHistSignalinTPCKaonforstandardcuts); | |
333 | fHistSignalinTPCKaonforITSPIDcuts= new TH2F("fHistSignalinTPCKaonforITSPIDcuts",";Pt[GeV/c];signal",kPtBins,binsPtDummy,100,-1.0,1.0) ; | |
334 | flist->Add(fHistSignalinTPCKaonforITSPIDcuts); | |
335 | ||
336 | fHistSignalinTPCAntiKaonforstandardcuts= new TH2F("fHistSignalinTPCAntiKaonforstandardcuts",";Pt[GeV/c];signal",kPtBins,binsPtDummy,100,-1.0,1.0) ; | |
337 | flist->Add(fHistSignalinTPCAntiKaonforstandardcuts); | |
338 | fHistSignalinTPCAntiKaonforITSPIDcuts= new TH2F("fHistSignalinTPCAntiKaonforITSPIDcuts",";Pt[GeV/c];signal",kPtBins,binsPtDummy,100,-1.0,1.0) ; | |
339 | flist->Add(fHistSignalinTPCAntiKaonforITSPIDcuts); | |
340 | ||
341 | fHistSignalinTPCProtonforstandardcuts= new TH2F("fHistSignalinTPCProtonforstandardcuts",";Pt[GeV/c];signal",kPtBins,binsPtDummy,100,-1.0,1.0) ; | |
342 | flist->Add(fHistSignalinTPCProtonforstandardcuts); | |
343 | fHistSignalinTPCProtonforITSPIDcuts= new TH2F("fHistSignalinTPCProtonforITSPIDcuts",";Pt[GeV/c];signal",kPtBins,binsPtDummy,100,-1.0,1.0) ; | |
344 | flist->Add(fHistSignalinTPCProtonforITSPIDcuts); | |
345 | ||
346 | ||
347 | fHistSignalinTPCAntiProtonforstandardcuts= new TH2F("fHistSignalinTPCAntiProtonforstandardcuts",";Pt[GeV/c];signal",kPtBins,binsPtDummy,100,-1.0,1.0) ; | |
348 | flist->Add(fHistSignalinTPCAntiProtonforstandardcuts); | |
349 | fHistSignalinTPCAntiProtonforITSPIDcuts= new TH2F("fHistSignalinTPCAntiProtonforITSPIDcuts",";Pt[GeV/c];signal",kPtBins,binsPtDummy,100,-1.0,1.0) ; | |
350 | flist->Add(fHistSignalinTPCAntiProtonforITSPIDcuts); | |
351 | ||
352 | fPionNTPCClusters=new TH2F("fPionNTPCClusters","fPionNTPCClusters;Pt [GeV/c];NClusters",kPtBins,binsPtDummy,20,60,160); | |
353 | flist->Add(fPionNTPCClusters); | |
354 | fAntiPionNTPCClusters=new TH2F("fAntiPionNTPCClusters","fAntiPionNTPCClusters;Pt [GeV/c];NClusters",kPtBins,binsPtDummy,20,60,160); | |
355 | flist->Add(fAntiPionNTPCClusters); | |
356 | fKaonNTPCClusters=new TH2F("fKaonNTPCClusters","fKaonNTPCClusters;Pt [GeV/c];NClusters",kPtBins,binsPtDummy,20,60,160); | |
357 | flist->Add(fKaonNTPCClusters); | |
358 | fAntiKaonNTPCClusters=new TH2F("fAntiKaonNTPCClusters","fAntiKaonNTPCClusters;Pt [GeV/c];NClusters",kPtBins,binsPtDummy,20,60,160); | |
359 | flist->Add(fAntiKaonNTPCClusters); | |
360 | fProtonNTPCClusters=new TH2F("fProtonNTPCClusters","fProtonNTPCClusters;Pt [GeV/c];NClusters",kPtBins,binsPtDummy,20,60,160); | |
361 | flist->Add(fProtonNTPCClusters); | |
362 | fAntiProtonNTPCClusters=new TH2F("fAntiProtonNTPCClusters","fAntiProtonNTPCClusters;Pt [GeV/c];NClusters",kPtBins,binsPtDummy,20,60,160); | |
363 | flist->Add(fAntiProtonNTPCClusters); | |
364 | ||
365 | ||
366 | fPionchi2=new TH2F("fPionchi2","fPionchi2;Pt [GeV/c];NClusters",kPtBins,binsPtDummy,20,0,6); | |
367 | flist->Add(fPionchi2); | |
368 | fAntiPionchi2=new TH2F("fAntiPionchi2","fAntiPionchi2;Pt [GeV/c];NClusters",kPtBins,binsPtDummy,20,0,6); | |
369 | flist->Add(fAntiPionchi2); | |
370 | fKaonchi2=new TH2F("fKaonchi2","fKaonchi2;Pt [GeV/c];NClusters",kPtBins,binsPtDummy,20,0,6); | |
371 | flist->Add(fKaonchi2); | |
372 | fAntiKaonchi2=new TH2F("fAntiKaonchi2","fAntiKaonchi2;Pt [GeV/c];NClusters",kPtBins,binsPtDummy,20,0,6); | |
373 | flist->Add(fAntiKaonchi2); | |
374 | fProtonchi2=new TH2F("fProtonchi2","fProtonchi2;Pt [GeV/c];NClusters",kPtBins,binsPtDummy,20,0,6); | |
375 | flist->Add(fProtonchi2); | |
376 | fAntiProtonchi2=new TH2F("fAntiProtonchi2","fAntiProtonchi2;Pt [GeV/c];NClusters",kPtBins,binsPtDummy,20,0,6); | |
377 | flist->Add(fAntiProtonchi2); | |
378 | ||
379 | if(fHIsettings) | |
380 | { | |
381 | fHistStandartMul=new TH1F("fHistStandartMul",";Ntracks;counts",300,0,3000); | |
382 | flist->Add(fHistStandartMul); | |
383 | fHistMytrackMul=new TH1F("fHistMytrackMul",";Ntracks;counts",300,0,3000); | |
384 | flist->Add(fHistMytrackMul); | |
385 | fHistStandartMulvSPD2=new TH2F("fHistStandartMulvSPD2",";Ntracks;nSPD2;counts",300,0,3000,300,0,3000); | |
386 | flist->Add(fHistStandartMulvSPD2); | |
387 | } | |
388 | else | |
389 | { | |
390 | fHistStandartMul=new TH1F("fHistStandartMul",";Ntracks;counts",300,0,300); | |
391 | flist->Add(fHistStandartMul); | |
392 | fHistMytrackMul=new TH1F("fHistMytrackMul",";Ntracks;counts",300,0,300); | |
393 | flist->Add(fHistMytrackMul); | |
394 | fHistStandartMulvSPD2=new TH2F("fHistStandartMulvSPD2",";Ntracks;nSPD2;counts",300,0,300,300,0,300); | |
395 | flist->Add(fHistStandartMulvSPD2); | |
396 | } | |
397 | fTracksCutmonitoring=new TH2F("fTracksCutmonitoring",";cut;pt[GeV/c];N_{entries}",4,0.5,4.5,kPtBins,binsPtDummy); | |
398 | fTracksCutmonitoring->GetXaxis()->SetBinLabel(1,"TPCin"); | |
399 | fTracksCutmonitoring->GetXaxis()->SetBinLabel(2,"standard"); | |
400 | fTracksCutmonitoring->GetXaxis()->SetBinLabel(3,"ITSpid"); | |
401 | fTracksCutmonitoring->GetXaxis()->SetBinLabel(4,"DCA"); | |
402 | flist->Add(fTracksCutmonitoring); | |
403 | ||
404 | ||
405 | ||
406 | ||
407 | if(!fMC) | |
408 | { | |
409 | ||
410 | Printf("end of CreateOutputObjects no MC"); | |
411 | PostData(1,flist); | |
412 | return; | |
413 | } | |
414 | ||
415 | ||
416 | ||
417 | fHistminsignalifPionPPrimary=new TH2F("HistminsignalifPionPPrimary",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
418 | flist->Add(fHistminsignalifPionPPrimary); | |
419 | fHistminsignalifKaonPPrimary=new TH2F("HistminsignalifKaonPPrimary",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
420 | flist->Add(fHistminsignalifKaonPPrimary); | |
421 | fHistminsignalifProtonPPrimary=new TH2F("HistminsignalifProtonPPrimary",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
422 | flist->Add(fHistminsignalifProtonPPrimary); | |
423 | fHistminsignalifProtonPPrimaryfake=new TH2F("HistminsignalifProtonPPrimaryfake",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
424 | flist->Add(fHistminsignalifProtonPPrimaryfake); | |
425 | ||
426 | fHistminsignalifAntiPionPPrimary=new TH2F("HistminsignalifAntiPionPPrimary",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
427 | flist->Add(fHistminsignalifAntiPionPPrimary); | |
428 | fHistminsignalifAntiKaonPPrimary=new TH2F("HistminsignalifAntiKaonPPrimary",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
429 | flist->Add(fHistminsignalifAntiKaonPPrimary); | |
430 | fHistminsignalifAntiProtonPPrimary=new TH2F("HistminsignalifAntiProtonPPrimary",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
431 | flist->Add(fHistminsignalifAntiProtonPPrimary); | |
432 | fHistminsignalifAntiProtonPPrimaryfake=new TH2F("HistminsignalifAntiProtonPPrimaryfake",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
433 | flist->Add(fHistminsignalifAntiProtonPPrimaryfake); | |
434 | ||
435 | fHistminsignalifPionPSecondary=new TH2F("HistminsignalifPionPSecondary",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
436 | flist->Add(fHistminsignalifPionPSecondary); | |
437 | fHistminsignalifKaonPSecondary=new TH2F("HistminsignalifKaonPSecondary",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
438 | flist->Add(fHistminsignalifKaonPSecondary); | |
439 | fHistminsignalifProtonPSecondaryWD=new TH2F("HistminsignalifProtonPSecondaryWD",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
440 | flist->Add(fHistminsignalifProtonPSecondaryWD); | |
441 | fHistminsignalifProtonPSecondaryHI=new TH2F("HistminsignalifProtonPSecondaryHI",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
442 | flist->Add(fHistminsignalifProtonPSecondaryHI); | |
443 | fHistminsignalifProtonPSecondaryRest=new TH2F("HistminsignalifProtonPSecondaryRest",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
444 | flist->Add(fHistminsignalifProtonPSecondaryRest); | |
445 | fHistminsignalifProtonPSecondaryWDfake=new TH2F("HistminsignalifProtonPSecondaryWDfake",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
446 | flist->Add(fHistminsignalifProtonPSecondaryWDfake); | |
447 | fHistminsignalifProtonPSecondaryHIfake=new TH2F("HistminsignalifProtonPSecondaryHIfake",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
448 | flist->Add(fHistminsignalifProtonPSecondaryHIfake); | |
449 | ||
450 | ||
451 | fHistminsignalifAntiPionPSecondary=new TH2F("HistminsignalifAntiPionPSecondary",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
452 | flist->Add(fHistminsignalifAntiPionPSecondary); | |
453 | fHistminsignalifAntiKaonPSecondary=new TH2F("HistminsignalifAntiKaonPSecondary",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
454 | flist->Add(fHistminsignalifAntiKaonPSecondary); | |
455 | fHistminsignalifAntiProtonPSecondaryWD=new TH2F("HistminsignalifAntiProtonPSecondaryWD",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
456 | flist->Add(fHistminsignalifAntiProtonPSecondaryWD); | |
457 | fHistminsignalifAntiProtonPSecondaryHI=new TH2F("HistminsignalifAntiProtonPSecondaryHI",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
458 | flist->Add(fHistminsignalifAntiProtonPSecondaryHI); | |
459 | fHistminsignalifAntiProtonPSecondaryRest=new TH2F("HistminsignalifAntiProtonPSecondaryRest",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
460 | flist->Add(fHistminsignalifAntiProtonPSecondaryRest); | |
461 | fHistminsignalifAntiProtonPSecondaryWDfake=new TH2F("HistminsignalifAntiProtonPSecondaryWDfake",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
462 | flist->Add(fHistminsignalifAntiProtonPSecondaryWDfake); | |
463 | fHistminsignalifAntiProtonPSecondaryHIfake=new TH2F("HistminsignalifAntiProtonPSecondaryHIfake",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
464 | flist->Add(fHistminsignalifAntiProtonPSecondaryHIfake); | |
465 | ||
466 | fHistminsignalifMuEPositiveP=new TH2F("HistminsignalifMuEPositiveP",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
467 | flist->Add(fHistminsignalifMuEPositiveP); | |
468 | fHistminsignalifMuENegativeP=new TH2F("HistminsignalifMuENegativeP",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
469 | flist->Add(fHistminsignalifMuENegativeP); | |
470 | ||
471 | ||
472 | fHistminsignalifPionPrimaryfake=new TH2F("HistminsignalifPionPrimaryfake",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
473 | flist->Add(fHistminsignalifPionPrimaryfake); | |
474 | fHistminsignalifKaonPrimaryfake=new TH2F("HistminsignalifKaonPrimaryfake",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
475 | flist->Add(fHistminsignalifKaonPrimaryfake); | |
476 | ||
477 | fHistminsignalifAntiPionPrimaryfake=new TH2F("HistminsignalifAntiPionPrimaryfake",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
478 | flist->Add(fHistminsignalifAntiPionPrimaryfake); | |
479 | fHistminsignalifAntiKaonPrimaryfake=new TH2F("HistminsignalifAntiKaonPrimaryfake",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
480 | flist->Add(fHistminsignalifAntiKaonPrimaryfake); | |
481 | ||
482 | ||
483 | fHistminsignalifPionSecondaryfake=new TH2F("HistminsignalifPionSecondaryfake",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
484 | flist->Add(fHistminsignalifPionSecondaryfake); | |
485 | fHistminsignalifKaonSecondaryfake=new TH2F("HistminsignalifKaonSecondaryfake",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
486 | flist->Add(fHistminsignalifKaonSecondaryfake); | |
487 | ||
488 | fHistminsignalifAntiPionSecondaryfake=new TH2F("HistminsignalifAntiPionSecondaryfake",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
489 | flist->Add(fHistminsignalifAntiPionSecondaryfake); | |
490 | fHistminsignalifAntiKaonSecondaryfake=new TH2F("HistminsignalifAntiKaonSecondaryfake",";Pt[GeV/c];log(dE_real)-log(dE_fit)",kPtBins,binsPtDummy,ny,-4,4); | |
491 | flist->Add(fHistminsignalifAntiKaonSecondaryfake); | |
492 | ||
493 | ||
494 | fHistminsignalifPionPMCPrimary=new TH1F("HistminsignalifPionPMCPrimary",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
495 | flist->Add(fHistminsignalifPionPMCPrimary); | |
496 | fHistminsignalifKaonPMCPrimary=new TH1F("HistminsignalifKaonPMCPrimary",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
497 | flist->Add(fHistminsignalifKaonPMCPrimary); | |
498 | fHistminsignalifProtonPMCPrimary=new TH1F("HistminsignalifProtonPMCPrimary",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
499 | flist->Add(fHistminsignalifProtonPMCPrimary); | |
500 | ||
501 | fHistminsignalifAntiPionPMCPrimary=new TH1F("HistminsignalifAntiPionPMCPrimary",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
502 | flist->Add(fHistminsignalifAntiPionPMCPrimary); | |
503 | fHistminsignalifAntiKaonPMCPrimary=new TH1F("HistminsignalifAntiKaonPMCPrimary",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
504 | flist->Add(fHistminsignalifAntiKaonPMCPrimary); | |
505 | fHistminsignalifAntiProtonPMCPrimary=new TH1F("HistminsignalifAntiProtonPMCPrimary",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
506 | flist->Add(fHistminsignalifAntiProtonPMCPrimary); | |
507 | ||
508 | ||
509 | fHistminsignalifPionPMCPrimaryBeforeEventCuts=new TH1F("HistminsignalifPionPMCPrimaryBeforeEventCuts",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
510 | flist->Add(fHistminsignalifPionPMCPrimaryBeforeEventCuts); | |
511 | fHistminsignalifKaonPMCPrimaryBeforeEventCuts=new TH1F("HistminsignalifKaonPMCPrimaryBeforeEventCuts",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
512 | flist->Add(fHistminsignalifKaonPMCPrimaryBeforeEventCuts); | |
513 | fHistminsignalifProtonPMCPrimaryBeforeEventCuts=new TH1F("HistminsignalifProtonPMCPrimaryBeforeEventCuts",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
514 | flist->Add(fHistminsignalifProtonPMCPrimaryBeforeEventCuts); | |
515 | ||
516 | fHistminsignalifAntiPionPMCPrimaryBeforeEventCuts=new TH1F("HistminsignalifAntiPionPMCPrimaryBeforeEventCuts",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
517 | flist->Add(fHistminsignalifAntiPionPMCPrimaryBeforeEventCuts); | |
518 | fHistminsignalifAntiKaonPMCPrimaryBeforeEventCuts=new TH1F("HistminsignalifAntiKaonPMCPrimaryBeforeEventCuts",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
519 | flist->Add(fHistminsignalifAntiKaonPMCPrimaryBeforeEventCuts); | |
520 | fHistminsignalifAntiProtonPMCPrimaryBeforeEventCuts=new TH1F("HistminsignalifAntiProtonPMCPrimaryBeforeEventCuts",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
521 | flist->Add(fHistminsignalifAntiProtonPMCPrimaryBeforeEventCuts); | |
522 | ||
523 | fHistminsignalifPionPMCPrimaryBeforeEventCutswithgoodZvertex=new TH1F("HistminsignalifPionPMCPrimaryBeforeEventCutswithgoodZvertex",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
524 | flist->Add(fHistminsignalifPionPMCPrimaryBeforeEventCutswithgoodZvertex); | |
525 | fHistminsignalifKaonPMCPrimaryBeforeEventCutswithgoodZvertex=new TH1F("HistminsignalifKaonPMCPrimaryBeforeEventCutswithgoodZvertex",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
526 | flist->Add(fHistminsignalifKaonPMCPrimaryBeforeEventCutswithgoodZvertex); | |
527 | fHistminsignalifProtonPMCPrimaryBeforeEventCutswithgoodZvertex=new TH1F("HistminsignalifProtonPMCPrimaryBeforeEventCutswithgoodZvertex",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
528 | flist->Add(fHistminsignalifProtonPMCPrimaryBeforeEventCutswithgoodZvertex); | |
529 | ||
530 | fHistminsignalifAntiPionPMCPrimaryBeforeEventCutswithgoodZvertex=new TH1F("HistminsignalifAntiPionPMCPrimaryBeforeEventCutswithgoodZvertex",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
531 | flist->Add(fHistminsignalifAntiPionPMCPrimaryBeforeEventCutswithgoodZvertex); | |
532 | fHistminsignalifAntiKaonPMCPrimaryBeforeEventCutswithgoodZvertex=new TH1F("HistminsignalifAntiKaonPMCPrimaryBeforeEventCutswithgoodZvertex",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
533 | flist->Add(fHistminsignalifAntiKaonPMCPrimaryBeforeEventCutswithgoodZvertex); | |
534 | fHistminsignalifAntiProtonPMCPrimaryBeforeEventCutswithgoodZvertex=new TH1F("HistminsignalifAntiProtonPMCPrimaryBeforeEventCutswithgoodZvertex",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
535 | flist->Add(fHistminsignalifAntiProtonPMCPrimaryBeforeEventCutswithgoodZvertex); | |
536 | ||
537 | fHistminsignalifPionPMCPrimaryAfterEventCutsBeforeVertexZ=new TH1F("HistminsignalifPionPMCPrimaryAfterEventCutsBeforeVertexZ",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
538 | flist->Add(fHistminsignalifPionPMCPrimaryAfterEventCutsBeforeVertexZ); | |
539 | fHistminsignalifKaonPMCPrimaryAfterEventCutsBeforeVertexZ=new TH1F("HistminsignalifKaonPMCPrimaryAfterEventCutsBeforeVertexZ",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
540 | flist->Add(fHistminsignalifKaonPMCPrimaryAfterEventCutsBeforeVertexZ); | |
541 | fHistminsignalifProtonPMCPrimaryAfterEventCutsBeforeVertexZ=new TH1F("HistminsignalifProtonPMCPrimaryAfterEventCutsBeforeVertexZ",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
542 | flist->Add(fHistminsignalifProtonPMCPrimaryAfterEventCutsBeforeVertexZ); | |
543 | ||
544 | fHistminsignalifAntiPionPMCPrimaryAfterEventCutsBeforeVertexZ=new TH1F("HistminsignalifAntiPionPMCPrimaryAfterEventCutsBeforeVertexZ",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
545 | flist->Add(fHistminsignalifAntiPionPMCPrimaryAfterEventCutsBeforeVertexZ); | |
546 | fHistminsignalifAntiKaonPMCPrimaryAfterEventCutsBeforeVertexZ=new TH1F("HistminsignalifAntiKaonPMCPrimaryAfterEventCutsBeforeVertexZ",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
547 | flist->Add(fHistminsignalifAntiKaonPMCPrimaryAfterEventCutsBeforeVertexZ); | |
548 | fHistminsignalifAntiProtonPMCPrimaryAfterEventCutsBeforeVertexZ=new TH1F("HistminsignalifAntiProtonPMCPrimaryAfterEventCutsBeforeVertexZ",";Pt[GeV/c]",kPtBins,binsPtDummy); | |
549 | flist->Add(fHistminsignalifAntiProtonPMCPrimaryAfterEventCutsBeforeVertexZ); | |
550 | ||
551 | ||
552 | ||
553 | fDCAXYZforcleanPionsMCPrimary=new TH3F("fDCAXYZforcleanPionsMCPrimary",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
554 | flist->Add(fDCAXYZforcleanPionsMCPrimary); | |
555 | fDCAXYZforcleanAntiPionsMCPrimary=new TH3F("fDCAXYZforcleanAntiPionsMCPrimary",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
556 | flist->Add(fDCAXYZforcleanAntiPionsMCPrimary); | |
557 | fDCAXYZforcleanProtonsMCPrimary=new TH3F("fDCAXYZforcleanProtonsMCPrimary",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
558 | flist->Add(fDCAXYZforcleanProtonsMCPrimary); | |
559 | fDCAXYZforcleanAntiProtonsMCPrimary=new TH3F("fDCAXYZforcleanAntiProtonsMCPrimary",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
560 | flist->Add(fDCAXYZforcleanAntiProtonsMCPrimary); | |
561 | ||
562 | fDCAXYZforcleanPionsWD=new TH3F("fDCAXYZforcleanPionsWD",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
563 | flist->Add(fDCAXYZforcleanPionsWD); | |
564 | fDCAXYZforcleanAntiPionsWD=new TH3F("fDCAXYZforcleanAntiPionsWD",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
565 | flist->Add(fDCAXYZforcleanAntiPionsWD); | |
566 | //Secondrary Protons weak deacy | |
567 | ||
568 | fDCAXYZforcleanProtonsWD=new TH3F("fDCAXYZforcleanProtonsWD",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
569 | flist->Add(fDCAXYZforcleanProtonsWD); | |
570 | fDCAXYZforcleanAntiProtonsWD=new TH3F("fDCAXYZforcleanAntiProtonsWD",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
571 | flist->Add(fDCAXYZforcleanAntiProtonsWD); | |
572 | ||
573 | fDCAXYZforcleanPionsHI=new TH3F("fDCAXYZforcleanPionsHI",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
574 | flist->Add(fDCAXYZforcleanPionsHI); | |
575 | fDCAXYZforcleanAntiPionsHI=new TH3F("fDCAXYZforcleanAntiPionsHI",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
576 | flist->Add(fDCAXYZforcleanAntiPionsHI); | |
577 | //Secondrary Protons Hadronic | |
578 | fDCAXYZforcleanProtonsHI=new TH3F("fDCAXYZforcleanProtonsHI",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
579 | flist->Add(fDCAXYZforcleanProtonsHI); | |
580 | fDCAXYZforcleanAntiProtonsHI=new TH3F("fDCAXYZforcleanAntiProtonsHI",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
581 | flist->Add(fDCAXYZforcleanAntiProtonsHI); | |
582 | //Secondrary Pions mu el | |
583 | fDCAXYZforcleanPionsMEPrimary=new TH3F("fDCAXYZforcleanPionsMEPrimary",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
584 | flist->Add(fDCAXYZforcleanPionsMEPrimary); | |
585 | fDCAXYZforcleanAntiPionsMEPrimary=new TH3F("fDCAXYZforcleanAntiPionsMEPrimary",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
586 | flist->Add(fDCAXYZforcleanAntiPionsMEPrimary); | |
587 | fDCAXYZforcleanPionsMESecondary=new TH3F("fDCAXYZforcleanPionsMESecondary",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
588 | flist->Add(fDCAXYZforcleanPionsMESecondary); | |
589 | fDCAXYZforcleanAntiPionsMESecondary=new TH3F("fDCAXYZforcleanAntiPionsMESecondary",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
590 | flist->Add(fDCAXYZforcleanAntiPionsMESecondary); | |
591 | ||
592 | fDCAXYZforcleanPionsR=new TH3F("fDCAXYZforcleanPionsR",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
593 | flist->Add(fDCAXYZforcleanPionsR); | |
594 | fDCAXYZforcleanAntiPionsR=new TH3F("fDCAXYZforcleanAntiPionsR",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
595 | flist->Add(fDCAXYZforcleanAntiPionsR); | |
596 | //Secondrary Protons Hadronic | |
597 | fDCAXYZforcleanProtonsR=new TH3F("fDCAXYZforcleanProtonsR",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
598 | flist->Add(fDCAXYZforcleanProtonsR); | |
599 | fDCAXYZforcleanAntiProtonsR=new TH3F("fDCAXYZforcleanAntiProtonsR",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
600 | flist->Add(fDCAXYZforcleanAntiProtonsR); | |
601 | ||
602 | ||
603 | ||
604 | fDCAXYZOpenforcleanPionsMCPrimary=new TH3F("fDCAXYZOpenforcleanPionsMCPrimary",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
605 | flist->Add(fDCAXYZOpenforcleanPionsMCPrimary); | |
606 | fDCAXYZOpenforcleanAntiPionsMCPrimary=new TH3F("fDCAXYZOpenforcleanAntiPionsMCPrimary",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
607 | flist->Add(fDCAXYZOpenforcleanAntiPionsMCPrimary); | |
608 | fDCAXYZOpenforcleanProtonsMCPrimary=new TH3F("fDCAXYZOpenforcleanProtonsMCPrimary",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
609 | flist->Add(fDCAXYZOpenforcleanProtonsMCPrimary); | |
610 | fDCAXYZOpenforcleanAntiProtonsMCPrimary=new TH3F("fDCAXYZOpenforcleanAntiProtonsMCPrimary",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
611 | flist->Add(fDCAXYZOpenforcleanAntiProtonsMCPrimary); | |
612 | ||
613 | fDCAXYZOpenforcleanPionsWD=new TH3F("fDCAXYZOpenforcleanPionsWD",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
614 | flist->Add(fDCAXYZOpenforcleanPionsWD); | |
615 | fDCAXYZOpenforcleanAntiPionsWD=new TH3F("fDCAXYZOpenforcleanAntiPionsWD",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
616 | flist->Add(fDCAXYZOpenforcleanAntiPionsWD); | |
617 | //Secondrary Protons weak deacy | |
618 | ||
619 | fDCAXYZOpenforcleanProtonsWD=new TH3F("fDCAXYZOpenforcleanProtonsWD",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
620 | flist->Add(fDCAXYZOpenforcleanProtonsWD); | |
621 | fDCAXYZOpenforcleanAntiProtonsWD=new TH3F("fDCAXYZOpenforcleanAntiProtonsWD",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
622 | flist->Add(fDCAXYZOpenforcleanAntiProtonsWD); | |
623 | ||
624 | fDCAXYZOpenforcleanPionsHI=new TH3F("fDCAXYZOpenforcleanPionsHI",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
625 | flist->Add(fDCAXYZOpenforcleanPionsHI); | |
626 | fDCAXYZOpenforcleanAntiPionsHI=new TH3F("fDCAXYZOpenforcleanAntiPionsHI",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
627 | flist->Add(fDCAXYZOpenforcleanAntiPionsHI); | |
628 | //Secondrary Protons Hadronic | |
629 | fDCAXYZOpenforcleanProtonsHI=new TH3F("fDCAXYZOpenforcleanProtonsHI",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
630 | flist->Add(fDCAXYZOpenforcleanProtonsHI); | |
631 | fDCAXYZOpenforcleanAntiProtonsHI=new TH3F("fDCAXYZOpenforcleanAntiProtonsHI",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
632 | flist->Add(fDCAXYZOpenforcleanAntiProtonsHI); | |
633 | //Secondrary Pions mu el | |
634 | ||
635 | fDCAXYZOpenforcleanPionsMEPrimary=new TH3F("fDCAXYZOpenforcleanPionsMEPrimary",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
636 | flist->Add(fDCAXYZOpenforcleanPionsMEPrimary); | |
637 | fDCAXYZOpenforcleanAntiPionsMEPrimary=new TH3F("fDCAXYZOpenforcleanAntiPionsMEPrimary",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
638 | flist->Add(fDCAXYZOpenforcleanAntiPionsMEPrimary); | |
639 | fDCAXYZOpenforcleanPionsMESecondary=new TH3F("fDCAXYZOpenforcleanPionsMESecondary",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
640 | flist->Add(fDCAXYZOpenforcleanPionsMESecondary); | |
641 | fDCAXYZOpenforcleanAntiPionsMESecondary=new TH3F("fDCAXYZOpenforcleanAntiPionsMESecondary",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
642 | flist->Add(fDCAXYZOpenforcleanAntiPionsMESecondary); | |
643 | ||
644 | fDCAXYZOpenforcleanPionsR=new TH3F("fDCAXYZOpenforcleanPionsR",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
645 | flist->Add(fDCAXYZOpenforcleanPionsR); | |
646 | fDCAXYZOpenforcleanAntiPionsR=new TH3F("fDCAXYZOpenforcleanAntiPionsR",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
647 | flist->Add(fDCAXYZOpenforcleanAntiPionsR); | |
648 | //Secondrary Protons Hadronic | |
649 | fDCAXYZOpenforcleanProtonsR=new TH3F("fDCAXYZOpenforcleanProtonsR",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
650 | flist->Add(fDCAXYZOpenforcleanProtonsR); | |
651 | fDCAXYZOpenforcleanAntiProtonsR=new TH3F("fDCAXYZOpenforcleanAntiProtonsR",";Pt[GeV/c];dcaxy[cm];dcaz[cm]",kPtBins,binsPtDummy,kDCABins,binsDCADummy,kDCABins,binsDCADummy); | |
652 | flist->Add(fDCAXYZOpenforcleanAntiProtonsR); | |
653 | ||
654 | ||
655 | ||
656 | fElectronsource=new TH2F("fElectronsource","electrons;Pt [GeV/c];Process",kPtBins,binsPtDummy, kMaxMCProcess,0,kMaxMCProcess); | |
657 | flist->Add(fElectronsource); | |
658 | fAntiElectronsource=new TH2F("fAntiElectronsource","positrons;Pt [GeV/c];Process",kPtBins,binsPtDummy, kMaxMCProcess,0,kMaxMCProcess); | |
659 | flist->Add(fAntiElectronsource); | |
660 | fMuonsource=new TH2F("fMuonsource","electrons;Pt [GeV/c];Process",kPtBins,binsPtDummy, kMaxMCProcess,0,kMaxMCProcess); | |
661 | flist->Add(fMuonsource); | |
662 | fAntiMuonsource=new TH2F("fAntiMuonsource","positrons;Pt [GeV/c];Process",kPtBins,binsPtDummy, kMaxMCProcess,0,kMaxMCProcess); | |
663 | flist->Add(fAntiMuonsource); | |
664 | ||
665 | fPrimaryElectronsMother=new TH1F("fPrimaryElectronsMother",";pdg code",4990,10.5,5000.5); | |
666 | flist->Add(fPrimaryElectronsMother); | |
667 | ||
668 | Double_t type[13]={-0.5,0.5,1.5,2.5,3.5,4.5,5.5,6.5,7.5,8.5,9.5,10.5,11.5}; | |
669 | Double_t cutlevel[10]={0.5,1.5,2.5,3.5,4.5,5.5,6.5,7.5,8.5,9.5}; | |
670 | ||
671 | fParticlesCutmonitoring=new TH3F("fParticlesCutmonitoring",";particle;cut;Pt [GeV/c]",12,type,9,cutlevel,kPtBins,binsPtDummy); | |
672 | ||
673 | fParticlesCutmonitoring->GetXaxis()->SetBinLabel(1,"pion"); | |
674 | fParticlesCutmonitoring->GetXaxis()->SetBinLabel(2,"kaon"); | |
675 | fParticlesCutmonitoring->GetXaxis()->SetBinLabel(3,"proton"); | |
676 | fParticlesCutmonitoring->GetXaxis()->SetBinLabel(4,"antipion"); | |
677 | fParticlesCutmonitoring->GetXaxis()->SetBinLabel(5,"antikaon"); | |
678 | fParticlesCutmonitoring->GetXaxis()->SetBinLabel(6,"antiproton"); | |
679 | fParticlesCutmonitoring->GetXaxis()->SetBinLabel(7,"pionfake"); | |
680 | fParticlesCutmonitoring->GetXaxis()->SetBinLabel(8,"kaonfake"); | |
681 | fParticlesCutmonitoring->GetXaxis()->SetBinLabel(9,"protonfake"); | |
682 | fParticlesCutmonitoring->GetXaxis()->SetBinLabel(10,"antipionfake"); | |
683 | fParticlesCutmonitoring->GetXaxis()->SetBinLabel(11,"antikaonfake"); | |
684 | fParticlesCutmonitoring->GetXaxis()->SetBinLabel(12,"antiprotonfake"); | |
685 | ||
686 | fParticlesCutmonitoring->GetYaxis()->SetBinLabel(1,"TPCin"); | |
687 | fParticlesCutmonitoring->GetYaxis()->SetBinLabel(2,"TPCrefit"); | |
688 | fParticlesCutmonitoring->GetYaxis()->SetBinLabel(3,"nTPCclu"); | |
689 | fParticlesCutmonitoring->GetYaxis()->SetBinLabel(4,"chi2"); | |
690 | fParticlesCutmonitoring->GetYaxis()->SetBinLabel(5,"ITSrefit"); | |
691 | fParticlesCutmonitoring->GetYaxis()->SetBinLabel(6,"SPDany"); | |
692 | fParticlesCutmonitoring->GetYaxis()->SetBinLabel(7,"standard"); | |
693 | fParticlesCutmonitoring->GetYaxis()->SetBinLabel(8,"ITSpid"); | |
694 | fParticlesCutmonitoring->GetYaxis()->SetBinLabel(9,"DCA"); | |
695 | flist->Add(fParticlesCutmonitoring); | |
696 | ||
697 | fVertexshift=new TH3F("fVertexshift",";#delta_{x};#delta_{y};#delta_{z}",50,-0.06,0.06,50,-0.06,0.06,50,-2,2); | |
698 | flist->Add(fVertexshift); | |
699 | ||
700 | Double_t deltapttpc[41]; | |
701 | Double_t deltaptall[41]; | |
702 | for(int i=0;i<41;i++) | |
703 | { | |
704 | deltapttpc[i]=-0.8+i*(1.6/40); | |
705 | deltaptall[i]=-0.2+i*(0.4/40); | |
706 | } | |
707 | fPtESDminusPtMCvPtESDafterallcuts= new TH3F("fPtESDminusPtMCvPtESDafterallcuts",";#delta_{PtESD-PtMC};PtESD;type",40,deltaptall,kPtBins,binsPtDummy,12,type); | |
708 | flist->Add(fPtESDminusPtMCvPtESDafterallcuts); | |
709 | fPtESDminusPtMCvPtESDafterallcuts->GetZaxis()->SetBinLabel(1,"pion"); | |
710 | fPtESDminusPtMCvPtESDafterallcuts->GetZaxis()->SetBinLabel(2,"kaon"); | |
711 | fPtESDminusPtMCvPtESDafterallcuts->GetZaxis()->SetBinLabel(3,"proton"); | |
712 | fPtESDminusPtMCvPtESDafterallcuts->GetZaxis()->SetBinLabel(4,"antipion"); | |
713 | fPtESDminusPtMCvPtESDafterallcuts->GetZaxis()->SetBinLabel(5,"antikaon"); | |
714 | fPtESDminusPtMCvPtESDafterallcuts->GetZaxis()->SetBinLabel(6,"antiproton"); | |
715 | fPtESDminusPtMCvPtESDafterallcuts->GetZaxis()->SetBinLabel(7,"pionfake"); | |
716 | fPtESDminusPtMCvPtESDafterallcuts->GetZaxis()->SetBinLabel(8,"kaonfake"); | |
717 | fPtESDminusPtMCvPtESDafterallcuts->GetZaxis()->SetBinLabel(9,"protonfake"); | |
718 | fPtESDminusPtMCvPtESDafterallcuts->GetZaxis()->SetBinLabel(10,"antipionfake"); | |
719 | fPtESDminusPtMCvPtESDafterallcuts->GetZaxis()->SetBinLabel(11,"antikaonfake"); | |
720 | fPtESDminusPtMCvPtESDafterallcuts->GetZaxis()->SetBinLabel(12,"antiprotonfake"); | |
721 | ||
722 | fPtESDminusPtMCvPtESDafterTPCcuts= new TH3F("fPtESDminusPtMCvPtESDafterTPCcuts",";#delta_{PtESD-PtMC};PtESD;type",40,deltapttpc,kPtBins,binsPtDummy,12,type); | |
723 | flist->Add(fPtESDminusPtMCvPtESDafterTPCcuts); | |
724 | fPtESDminusPtMCvPtESDafterTPCcuts->GetZaxis()->SetBinLabel(1,"pion"); | |
725 | fPtESDminusPtMCvPtESDafterTPCcuts->GetZaxis()->SetBinLabel(2,"kaon"); | |
726 | fPtESDminusPtMCvPtESDafterTPCcuts->GetZaxis()->SetBinLabel(3,"proton"); | |
727 | fPtESDminusPtMCvPtESDafterTPCcuts->GetZaxis()->SetBinLabel(4,"antipion"); | |
728 | fPtESDminusPtMCvPtESDafterTPCcuts->GetZaxis()->SetBinLabel(5,"antikaon"); | |
729 | fPtESDminusPtMCvPtESDafterTPCcuts->GetZaxis()->SetBinLabel(6,"antiproton"); | |
730 | fPtESDminusPtMCvPtESDafterTPCcuts->GetZaxis()->SetBinLabel(7,"pionfake"); | |
731 | fPtESDminusPtMCvPtESDafterTPCcuts->GetZaxis()->SetBinLabel(8,"kaonfake"); | |
732 | fPtESDminusPtMCvPtESDafterTPCcuts->GetZaxis()->SetBinLabel(9,"protonfake"); | |
733 | fPtESDminusPtMCvPtESDafterTPCcuts->GetZaxis()->SetBinLabel(10,"antipionfake"); | |
734 | fPtESDminusPtMCvPtESDafterTPCcuts->GetZaxis()->SetBinLabel(11,"antikaonfake"); | |
735 | fPtESDminusPtMCvPtESDafterTPCcuts->GetZaxis()->SetBinLabel(12,"antiprotonfake"); | |
736 | ||
737 | ||
738 | fMulESDMulMCVz=new TH3F("fMulESDMulMCVz",";NtracksESD;NparticlesMC;Vrt_z ",50,0,50,100,0,100,20,-10,10); | |
739 | flist->Add(fMulESDMulMCVz); | |
740 | PostData(1, flist); | |
741 | Printf("end of CreateOutputObjects with MC"); | |
742 | } | |
743 | //________________________________________________________________________ | |
744 | void AliAnalysisChargedHadronSpectraITSTruncatedMeanTask::LocalInit() | |
745 | { | |
746 | //LocalInit | |
747 | Printf("end of LocalInit"); | |
748 | } | |
749 | ||
750 | //________________________________________________________________________ | |
751 | void AliAnalysisChargedHadronSpectraITSTruncatedMeanTask::UserExec(Option_t *) | |
752 | { | |
753 | //UserExec | |
754 | Bool_t isphysevent=0; | |
755 | Bool_t isgoodvertex=0; | |
756 | Bool_t isvxerteinZ=0; | |
757 | fESD = dynamic_cast<AliESDEvent*> (InputEvent()); | |
758 | if (!fESD) | |
759 | { | |
760 | Printf("ERROR: fESD not available"); | |
761 | return; | |
762 | } | |
763 | ||
764 | Float_t refmultiplicity=fCutsMul->CountAcceptedTracks(fESD); | |
765 | if(!fSPD) | |
766 | { | |
767 | if(fLowMultiplicity>-1) | |
768 | { | |
769 | if(refmultiplicity<fLowMultiplicity) | |
770 | return; | |
771 | } | |
772 | if(fUpMultiplicity>-1) | |
773 | { | |
774 | if(refmultiplicity>fUpMultiplicity) | |
775 | return; | |
776 | } | |
777 | } | |
778 | AliStack* stack=0x0; | |
779 | Double_t mcXvertex=0.0; | |
780 | Double_t mcYvertex=0.0; | |
781 | Double_t mcZvertex=0.0; | |
782 | ||
783 | if(fMC) | |
784 | { | |
785 | AliMCEvent* mcEvent = (AliMCEvent*) MCEvent(); | |
786 | //Printf("MC particles: %d", mcEvent->GetNumberOfTracks()); | |
787 | stack = mcEvent->Stack(); | |
788 | mcXvertex=mcEvent->GetPrimaryVertex()->GetX(); | |
789 | mcYvertex=mcEvent->GetPrimaryVertex()->GetY(); | |
790 | mcZvertex=mcEvent->GetPrimaryVertex()->GetZ(); | |
791 | } | |
792 | ||
793 | ||
794 | fHistStats->Fill(0); | |
795 | //Event selection | |
796 | //if( ((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()==0) | |
797 | UInt_t isSelected = 0; | |
798 | if(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))) | |
799 | isSelected=((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected(); | |
800 | //Printf("Mask_selection %u %u", isSelected,AliVEvent::kMB); | |
801 | if(!(isSelected&AliVEvent::kMB)) | |
802 | { | |
803 | isphysevent=0; | |
804 | // Printf("No phys event.........\n"); | |
805 | } | |
806 | else | |
807 | { | |
808 | isphysevent=1; | |
809 | if(!fHIsettings) | |
810 | fHistStats->Fill(1); | |
811 | } | |
812 | if(isphysevent&&fHIsettings) | |
813 | { | |
814 | AliCentrality *centrality = fESD->GetCentrality(); | |
815 | if(!(fLowCentrality<0.0)&&fUpCentrality>0.0) | |
816 | { | |
817 | if(!centrality->IsEventInCentralityClass(fLowCentrality,fUpCentrality,"V0M")) | |
818 | return; | |
819 | else | |
820 | fHistStats->Fill(1); | |
821 | } | |
822 | } | |
823 | //Printf("MC 1"); | |
824 | //Good vertex | |
825 | const AliESDVertex *vertex = 0x0; | |
826 | if(isphysevent) | |
827 | { | |
828 | vertex = fESD->GetPrimaryVertexTracks(); | |
829 | if(vertex->GetNContributors()<1) | |
830 | { | |
831 | // SPD vertex | |
832 | vertex = fESD->GetPrimaryVertexSPD(); | |
833 | if(vertex->GetNContributors()<1) | |
834 | { | |
835 | //Printf("No good Vertex.........\n"); | |
836 | isgoodvertex=0; | |
837 | } | |
838 | else | |
839 | { | |
840 | isgoodvertex=1; | |
841 | //fHistStats->Fill(2); | |
842 | fHistZVertexBeforeCut->Fill(vertex ->GetZ()); | |
843 | fHistXYVertexBeforeCut->Fill(vertex ->GetX(),vertex ->GetY()); | |
844 | } | |
845 | } | |
846 | else | |
847 | { | |
848 | isgoodvertex=1; | |
849 | //fHistStats->Fill(2); | |
850 | fHistZVertexBeforeCut->Fill(vertex ->GetZ()); | |
851 | fHistXYVertexBeforeCut->Fill(vertex ->GetX(),vertex ->GetY()); | |
852 | } | |
853 | if(isgoodvertex&&fUsePilerejection) | |
854 | { | |
855 | if(fESD->IsPileupFromSPDInMultBins()) | |
856 | isgoodvertex=0; | |
857 | } | |
858 | if(isgoodvertex) | |
859 | { | |
860 | if(TMath::Abs(vertex ->GetZ())>10.0) | |
861 | { | |
862 | //Printf("No good Z of Vertex.........\n"); | |
863 | isvxerteinZ=0; | |
864 | } | |
865 | else | |
866 | isvxerteinZ=1; | |
867 | } | |
868 | } | |
869 | ||
870 | ||
871 | if(fdovertexrescuts&&fMC) | |
872 | { | |
873 | if(vertex) | |
874 | { | |
875 | cout<<TMath::Abs(vertex->GetX()-mcXvertex)<<" "<<TMath::Abs(vertex->GetY()-mcYvertex)<<" "<<TMath::Abs(vertex->GetZ()-mcZvertex)<<endl; | |
876 | if(TMath::Abs(vertex->GetX()-mcXvertex)>0.015||TMath::Abs(vertex->GetY()-mcYvertex)>0.015||TMath::Abs(vertex->GetZ()-mcZvertex)>0.15) | |
877 | isvxerteinZ=0; | |
878 | } | |
879 | } | |
880 | Float_t spdCorr=-1.0; | |
881 | if(isgoodvertex) | |
882 | { | |
883 | const AliMultiplicity *mult = fESD->GetMultiplicity(); | |
884 | Float_t nClusters[6]={0.0,0.0,0.0,0.0,0.0,0.0}; | |
885 | for(Int_t ilay=0; ilay<6; ilay++) | |
886 | { | |
887 | nClusters[ilay] = (Float_t)mult->GetNumberOfITSClusters(ilay); | |
888 | } | |
889 | spdCorr = AliESDUtils::GetCorrSPD2(nClusters[1],vertex->GetZ()); | |
890 | if(fSPD) | |
891 | { | |
892 | if(fLowMultiplicity>-1) | |
893 | { | |
894 | if(((Int_t)spdCorr)<fLowMultiplicity) | |
895 | { | |
896 | PostData(1, flist); | |
897 | return; | |
898 | } | |
899 | } | |
900 | if(fUpMultiplicity>-1) | |
901 | { | |
902 | if(((Int_t)spdCorr)>fUpMultiplicity) | |
903 | { | |
904 | PostData(1, flist); | |
905 | return; | |
906 | } | |
907 | } | |
908 | } | |
909 | fHistStats->Fill(2); | |
910 | } | |
911 | ||
912 | //Printf("MC 2"); | |
913 | Int_t fMCmult=0; | |
914 | if(stack&&fMC)//Looping over MC information of all events | |
915 | { | |
916 | Float_t minpt=0.0; | |
917 | Float_t maxpt=0.0; | |
918 | Float_t mineta=0.0; | |
919 | Float_t maxeta=0.0; | |
920 | //Printf("MC 12"); | |
921 | fCutsMul->GetPtRange(minpt,maxpt); | |
922 | fCutsMul->GetEtaRange(mineta,maxeta); | |
923 | fHistStats->Fill(8); | |
924 | if(TMath::Abs(mcZvertex)<10.0) | |
925 | fHistStats->Fill(9); | |
926 | for (int imc=0;imc<stack->GetNtrack();imc++) | |
927 | { | |
928 | if(!(stack->IsPhysicalPrimary(imc))) | |
929 | continue; | |
930 | TParticle *particleMC = stack->Particle(imc); | |
931 | if(!particleMC) | |
932 | continue; | |
933 | Int_t pdgcodeMC = particleMC->GetPdgCode(); | |
934 | if(!(pdgcodeMC==211||pdgcodeMC==-211||pdgcodeMC==321||pdgcodeMC==-321||pdgcodeMC==2212||pdgcodeMC==-2212)) | |
935 | continue; | |
936 | if(particleMC->Pt()>minpt&&particleMC->Pt()<maxpt&&particleMC->Eta()>mineta&&particleMC->Eta()<maxeta) | |
937 | fMCmult++; | |
938 | if (TMath::Abs(particleMC->Y())>fYCut) | |
939 | continue; | |
940 | if (particleMC->Pt()>2.0) | |
941 | continue; | |
942 | //Printf("%d aa",imc); | |
943 | // Printf("MC 22"); | |
944 | if(pdgcodeMC==211) | |
945 | fHistminsignalifPionPMCPrimaryBeforeEventCuts->Fill(particleMC->Pt()); | |
946 | if(pdgcodeMC==-211) | |
947 | fHistminsignalifAntiPionPMCPrimaryBeforeEventCuts->Fill(particleMC->Pt()); | |
948 | if(pdgcodeMC==321) | |
949 | fHistminsignalifKaonPMCPrimaryBeforeEventCuts->Fill(particleMC->Pt()); | |
950 | if(pdgcodeMC==-321) | |
951 | fHistminsignalifAntiKaonPMCPrimaryBeforeEventCuts->Fill(particleMC->Pt()); | |
952 | if(pdgcodeMC==2212) | |
953 | fHistminsignalifProtonPMCPrimaryBeforeEventCuts->Fill(particleMC->Pt()); | |
954 | if(pdgcodeMC==-2212) | |
955 | fHistminsignalifAntiProtonPMCPrimaryBeforeEventCuts->Fill(particleMC->Pt()); | |
956 | ||
957 | if(TMath::Abs(mcZvertex)<10.0) | |
958 | { | |
959 | if(pdgcodeMC==211) | |
960 | fHistminsignalifPionPMCPrimaryBeforeEventCutswithgoodZvertex->Fill(particleMC->Pt()); | |
961 | if(pdgcodeMC==-211) | |
962 | fHistminsignalifAntiPionPMCPrimaryBeforeEventCutswithgoodZvertex->Fill(particleMC->Pt()); | |
963 | if(pdgcodeMC==321) | |
964 | fHistminsignalifKaonPMCPrimaryBeforeEventCutswithgoodZvertex->Fill(particleMC->Pt()); | |
965 | if(pdgcodeMC==-321) | |
966 | fHistminsignalifAntiKaonPMCPrimaryBeforeEventCutswithgoodZvertex->Fill(particleMC->Pt()); | |
967 | if(pdgcodeMC==2212) | |
968 | fHistminsignalifProtonPMCPrimaryBeforeEventCutswithgoodZvertex->Fill(particleMC->Pt()); | |
969 | if(pdgcodeMC==-2212) | |
970 | fHistminsignalifAntiProtonPMCPrimaryBeforeEventCutswithgoodZvertex->Fill(particleMC->Pt()); | |
971 | } | |
972 | if(!isphysevent) | |
973 | continue; | |
974 | if(!isgoodvertex) | |
975 | continue; | |
976 | if(pdgcodeMC==211) | |
977 | fHistminsignalifPionPMCPrimaryAfterEventCutsBeforeVertexZ->Fill(particleMC->Pt()); | |
978 | if(pdgcodeMC==-211) | |
979 | fHistminsignalifAntiPionPMCPrimaryAfterEventCutsBeforeVertexZ->Fill(particleMC->Pt()); | |
980 | if(pdgcodeMC==321) | |
981 | fHistminsignalifKaonPMCPrimaryAfterEventCutsBeforeVertexZ->Fill(particleMC->Pt()); | |
982 | if(pdgcodeMC==-321) | |
983 | fHistminsignalifAntiKaonPMCPrimaryAfterEventCutsBeforeVertexZ->Fill(particleMC->Pt()); | |
984 | if(pdgcodeMC==2212) | |
985 | fHistminsignalifProtonPMCPrimaryAfterEventCutsBeforeVertexZ->Fill(particleMC->Pt()); | |
986 | if(pdgcodeMC==-2212) | |
987 | fHistminsignalifAntiProtonPMCPrimaryAfterEventCutsBeforeVertexZ->Fill(particleMC->Pt()); | |
988 | if(!isvxerteinZ) | |
989 | continue; | |
990 | if(pdgcodeMC==211) | |
991 | { | |
992 | fHistminsignalifPionPMCPrimary->Fill(particleMC->Pt()); | |
993 | } | |
994 | if(pdgcodeMC==-211) | |
995 | { | |
996 | fHistminsignalifAntiPionPMCPrimary->Fill(particleMC->Pt()); | |
997 | } | |
998 | if(pdgcodeMC==321) | |
999 | { | |
1000 | fHistminsignalifKaonPMCPrimary->Fill(particleMC->Pt()); | |
1001 | } | |
1002 | if(pdgcodeMC==-321) | |
1003 | { | |
1004 | fHistminsignalifAntiKaonPMCPrimary->Fill(particleMC->Pt()); | |
1005 | } | |
1006 | if(pdgcodeMC==2212) | |
1007 | { | |
1008 | fHistminsignalifProtonPMCPrimary->Fill(particleMC->Pt()); | |
1009 | } | |
1010 | if(pdgcodeMC==-2212) | |
1011 | { | |
1012 | fHistminsignalifAntiProtonPMCPrimary->Fill(particleMC->Pt()); | |
1013 | } | |
1014 | ||
1015 | } | |
1016 | } | |
1017 | else if(fMC) | |
1018 | return; | |
1019 | else | |
1020 | Printf("Data mode \n"); | |
1021 | ||
1022 | if(!(isphysevent&&isgoodvertex&&isvxerteinZ)) | |
1023 | { | |
1024 | //Printf("No Good event.........\n"); | |
1025 | PostData(1, flist); | |
1026 | //Printf("end of Exec"); | |
1027 | return; | |
1028 | } | |
1029 | ||
1030 | ||
1031 | ||
1032 | fHistStandartMulvSPD2->Fill(refmultiplicity,spdCorr); | |
1033 | fHistStats->Fill(3); | |
1034 | ||
1035 | fHistZVertexAfterCut->Fill(vertex ->GetZ()); | |
1036 | fHistXYVertexAfterCut->Fill(vertex ->GetX(),vertex ->GetY()); | |
1037 | ||
1038 | if(fMC) | |
1039 | { | |
1040 | fVertexshift->Fill(vertex->GetX()-mcXvertex,vertex->GetY()-mcYvertex,vertex->GetZ()-mcZvertex); | |
1041 | fMulESDMulMCVz->Fill(refmultiplicity,fMCmult,vertex->GetZ()); | |
1042 | } | |
1043 | if(fCuts==0) | |
1044 | { | |
1045 | //Printf("No CUTS Defined.........\n"); | |
1046 | PostData(1, flist); | |
1047 | //Printf("end of Exec"); | |
1048 | return; | |
1049 | } | |
1050 | ||
1051 | //Printf("There are %d tracks in this event", fESD->GetNumberOfTracks()); | |
1052 | Int_t nTracks=fESD->GetNumberOfTracks(); | |
1053 | ||
1054 | Int_t mynumberoftracks=0; | |
1055 | AliESDtrack *trackESD=0; | |
1056 | ||
1057 | const Float_t pionmass=AliPID::ParticleMass(2); | |
1058 | const Float_t kaonmass=AliPID::ParticleMass(3); | |
1059 | const Float_t protonmass=AliPID::ParticleMass(4); | |
1060 | ||
1061 | for(int tr1=0;tr1<nTracks;tr1++) | |
1062 | { | |
1063 | ||
1064 | trackESD=fESD->GetTrack(tr1); | |
1065 | //fHistStats->Fill(2); | |
1066 | ||
1067 | Double_t pt=trackESD->Pt()*trackESD->GetSign(); | |
1068 | Double_t p=trackESD->P(); | |
1069 | Double_t eta=trackESD->Eta(); | |
1070 | Double_t phi=trackESD->Phi(); | |
1071 | Float_t dcaxy=0.0; | |
1072 | Float_t dcaz=0.0; | |
1073 | trackESD->GetImpactParameters(dcaxy,dcaz); | |
1074 | Double_t pz=trackESD->Pz(); | |
1075 | UShort_t nTPCclusters=trackESD->GetTPCNcls(); | |
1076 | Float_t chi2=trackESD->GetTPCchi2(); | |
1077 | if(nTPCclusters>0) | |
1078 | chi2=chi2/((Float_t)nTPCclusters); | |
1079 | else | |
1080 | chi2=-1.0; | |
1081 | if(!trackESD->IsOn(AliESDtrack::kTPCin)) | |
1082 | continue; | |
1083 | //Y assumtion | |
1084 | Float_t yforpion=0.5*TMath::Log((TMath::Sqrt(pionmass*pionmass+p*p)+pz)/(TMath::Sqrt(pionmass*pionmass+p*p)-pz)); | |
1085 | Float_t yforkaon=0.5*TMath::Log((TMath::Sqrt(kaonmass*kaonmass+p*p)+pz)/(TMath::Sqrt(kaonmass*kaonmass+p*p)-pz)); | |
1086 | Float_t yforproton=0.5*TMath::Log((TMath::Sqrt(protonmass*protonmass+p*p)+pz)/(TMath::Sqrt(protonmass*protonmass+p*p)-pz)); | |
1087 | ||
1088 | if(TMath::Abs(yforpion)>fYCut&&TMath::Abs(yforkaon)>fYCut&&TMath::Abs(yforproton)>fYCut) //go trought one y cut | |
1089 | continue; | |
1090 | Int_t label=-1; | |
1091 | if(fMC) | |
1092 | label=trackESD->GetLabel(); | |
1093 | //if(label<0) | |
1094 | // Printf("label %d %f %f %f %f %d %f %f\n",label,p,pt,eta,chi2,nTPCclusters,dcaxy,dcaz); | |
1095 | Int_t pdgcode=0; | |
1096 | Int_t primary=0; | |
1097 | Double_t chargeMC=1.0; | |
1098 | Float_t etaMC=10.0; | |
1099 | Float_t ptMC=10.0; | |
1100 | Int_t uniqueID=-1; | |
1101 | Int_t pdgcodefake=0; | |
1102 | Int_t primaryfake=0; | |
1103 | TParticle *particle2=0x0; | |
1104 | if(label>=0&&stack&&fMC) | |
1105 | { | |
1106 | primary=stack->IsPhysicalPrimary(TMath::Abs(label)); | |
1107 | particle2 = stack->Particle(TMath::Abs(label)); | |
1108 | pdgcode=particle2->GetPdgCode(); | |
1109 | chargeMC=particle2->GetPDG(0)->Charge()/3.0; | |
1110 | etaMC=particle2->Eta(); | |
1111 | ptMC=particle2->Pt(); | |
1112 | uniqueID=particle2->GetUniqueID(); | |
1113 | } | |
1114 | if(label<0&&stack&&fMC) | |
1115 | { | |
1116 | primaryfake=stack->IsPhysicalPrimary(TMath::Abs(label)); | |
1117 | particle2 = stack->Particle(TMath::Abs(label)); | |
1118 | pdgcodefake=particle2->GetPdgCode(); | |
1119 | uniqueID=particle2->GetUniqueID(); | |
1120 | ||
1121 | } | |
1122 | ||
1123 | Int_t typeParticle=-10; | |
1124 | if((primaryfake||primary)) | |
1125 | { | |
1126 | ||
1127 | if((pdgcodefake==211||pdgcode==211)&&TMath::Abs(yforpion)<fYCut) | |
1128 | typeParticle=0; | |
1129 | if((pdgcodefake==321||pdgcode==321)&&TMath::Abs(yforkaon)<fYCut) | |
1130 | typeParticle=1; | |
1131 | if((pdgcodefake==2212||pdgcode==2212)&&TMath::Abs(yforproton)<fYCut) | |
1132 | typeParticle=2; | |
1133 | if((pdgcodefake==-211||pdgcode==-211)&&TMath::Abs(yforpion)<fYCut) | |
1134 | typeParticle=3; | |
1135 | if((pdgcodefake==-321||pdgcode==-321)&&TMath::Abs(yforkaon)<fYCut) | |
1136 | typeParticle=4; | |
1137 | if((pdgcodefake==-2212||pdgcode==-2212)&&TMath::Abs(yforproton)<fYCut) | |
1138 | typeParticle=5; | |
1139 | ||
1140 | if(primaryfake) | |
1141 | typeParticle+=6; | |
1142 | } | |
1143 | ||
1144 | fTracksCutmonitoring->Fill(1,TMath::Abs(pt)); | |
1145 | if(fMC) | |
1146 | { | |
1147 | fParticlesCutmonitoring->Fill(typeParticle,1,TMath::Abs(pt)); | |
1148 | if(trackESD->IsOn(AliESDtrack::kTPCrefit)) | |
1149 | { | |
1150 | fParticlesCutmonitoring->Fill(typeParticle,2,TMath::Abs(pt)); | |
1151 | if(nTPCclusters>70) | |
1152 | { | |
1153 | fParticlesCutmonitoring->Fill(typeParticle,3,TMath::Abs(pt)); | |
1154 | if(chi2<4.0) | |
1155 | { | |
1156 | fParticlesCutmonitoring->Fill(typeParticle,4,TMath::Abs(pt)); | |
1157 | fPtESDminusPtMCvPtESDafterTPCcuts->Fill(TMath::Abs(pt)-particle2->Pt(),TMath::Abs(pt),typeParticle); | |
1158 | if(trackESD->IsOn(AliESDtrack::kITSrefit)) | |
1159 | { | |
1160 | fParticlesCutmonitoring->Fill(typeParticle,5,TMath::Abs(pt)); | |
1161 | if(trackESD->HasPointOnITSLayer(0)||trackESD->HasPointOnITSLayer(1)) | |
1162 | { | |
1163 | fParticlesCutmonitoring->Fill(typeParticle,6,TMath::Abs(pt)); | |
1164 | } | |
1165 | } | |
1166 | } | |
1167 | ||
1168 | } | |
1169 | } | |
1170 | } | |
1171 | ||
1172 | ||
1173 | fHistPhiPtBeforeCuts->Fill(phi,pt);//phi pt | |
1174 | fHistEtaPtBeforeCuts->Fill(eta,pt); | |
1175 | fHistDCABeforeCuts->Fill(dcaxy,dcaz); | |
1176 | ||
1177 | ||
1178 | ||
1179 | ||
1180 | //standart cuts | |
1181 | if(fCuts->AcceptTrack(trackESD)==kFALSE) | |
1182 | continue; | |
1183 | fTracksCutmonitoring->Fill(2,TMath::Abs(pt)); | |
1184 | if(fMC) | |
1185 | fParticlesCutmonitoring->Fill(typeParticle,7,TMath::Abs(pt)); | |
1186 | //Tpc pid cut for debug | |
1187 | Double_t pinTPC=trackESD->GetTPCInnerParam()->GetP();//momentum in primary vertex taken from TPC tracking | |
1188 | Double_t pinTPCglobal=trackESD->GetInnerParam()->GetP();//momentum at the inner wall of the TPC taken from global tracking | |
1189 | Float_t sigKaon= fESDpid->GetTPCResponse().GetExpectedSignal(pinTPCglobal, AliPID::kKaon); | |
1190 | Float_t sigProton= fESDpid->GetTPCResponse().GetExpectedSignal(pinTPCglobal, AliPID::kProton); | |
1191 | Double_t tpcSignal =trackESD ->GetTPCsignal(); | |
1192 | ||
1193 | if(fTPCPIDCUT) | |
1194 | { | |
1195 | if(fTPCPIDCUT->AcceptTrack(trackESD,fESD)==kFALSE) | |
1196 | continue; | |
1197 | } | |
1198 | ||
1199 | Bool_t cutDCA=SelectOnImpPar(trackESD); | |
1200 | ||
1201 | //fHistStats->Fill(4); | |
1202 | Double_t tmpQESD[4]={-1.0,-1.0,-1.0,-1.0}; | |
1203 | trackESD->GetITSdEdxSamples(tmpQESD); | |
1204 | if(fCorrectSDD&&fMC) | |
1205 | CorrectSDD(tmpQESD); | |
1206 | if(fCorrectSSD&&fMC) | |
1207 | CorrectSSD(tmpQESD); | |
1208 | Int_t nSSDSDD=0; | |
1209 | Int_t nSSDSDDall=0; | |
1210 | ||
1211 | ||
1212 | ||
1213 | if(TMath::Abs(yforpion)<fYCut&&cutDCA) | |
1214 | fHistNtrackwithstandardcuts->Fill(TMath::Abs(pt),0); | |
1215 | if(TMath::Abs(yforkaon)<fYCut&&cutDCA) | |
1216 | { | |
1217 | fHistNtrackwithstandardcuts->Fill(TMath::Abs(pt),1); | |
1218 | if(pt>0.0) | |
1219 | fHistSignalinTPCKaonforstandardcuts->Fill(pt,(tpcSignal-sigKaon)/sigKaon); | |
1220 | else | |
1221 | fHistSignalinTPCAntiKaonforstandardcuts->Fill(TMath::Abs(pt),(tpcSignal-sigKaon)/sigKaon); | |
1222 | } | |
1223 | if(TMath::Abs(yforproton)<fYCut&&cutDCA) | |
1224 | { | |
1225 | fHistNtrackwithstandardcuts->Fill(TMath::Abs(pt),2); | |
1226 | if(pt>0.0) | |
1227 | fHistSignalinTPCProtonforstandardcuts->Fill(pt,(tpcSignal-sigProton)/sigProton); | |
1228 | else | |
1229 | fHistSignalinTPCAntiProtonforstandardcuts->Fill(TMath::Abs(pt),(tpcSignal-sigProton)/sigProton); | |
1230 | } | |
1231 | ||
1232 | ||
1233 | for (int iL=0;iL<4;iL++) | |
1234 | { | |
1235 | if(tmpQESD[iL]>fchargeCut) | |
1236 | nSSDSDD++; | |
1237 | if(tmpQESD[iL]>0.0) | |
1238 | nSSDSDDall++; | |
1239 | } | |
1240 | if(nSSDSDDall>=3) | |
1241 | fHistStats->Fill(4); | |
1242 | //ITS PId cut | |
1243 | if(nSSDSDD<3) | |
1244 | { | |
1245 | // cout<<"BAD "<<stack->Particle(TMath::Abs(trackESD->GetLabel()))->GetPdgCode()<<endl; | |
1246 | continue; | |
1247 | } | |
1248 | fTracksCutmonitoring->Fill(3,TMath::Abs(pt)); | |
1249 | if(fMC) | |
1250 | fParticlesCutmonitoring->Fill(typeParticle,8,TMath::Abs(pt)); | |
1251 | if(cutDCA) | |
1252 | { | |
1253 | fTracksCutmonitoring->Fill(4,TMath::Abs(pt)); | |
1254 | if(fMC) | |
1255 | { | |
1256 | fParticlesCutmonitoring->Fill(typeParticle,9,TMath::Abs(pt)); | |
1257 | fPtESDminusPtMCvPtESDafterallcuts->Fill(TMath::Abs(pt)-particle2->Pt(),TMath::Abs(pt),typeParticle); | |
1258 | } | |
1259 | } | |
1260 | if(TMath::Abs(yforpion)<fYCut&&cutDCA) | |
1261 | { | |
1262 | fHistNtrackwithITSPIDcuts->Fill(TMath::Abs(pt),0); | |
1263 | if(pt>0.0) | |
1264 | { | |
1265 | fPionNTPCClusters->Fill(pt,nTPCclusters); | |
1266 | fPionchi2->Fill(pt,chi2); | |
1267 | } | |
1268 | else | |
1269 | { | |
1270 | fAntiPionNTPCClusters->Fill(TMath::Abs(pt),nTPCclusters); | |
1271 | fAntiPionchi2->Fill(TMath::Abs(pt),chi2); | |
1272 | } | |
1273 | } | |
1274 | if(TMath::Abs(yforkaon)<fYCut&&cutDCA) | |
1275 | { | |
1276 | fHistNtrackwithITSPIDcuts->Fill(TMath::Abs(pt),1); | |
1277 | if(pt>0.0) | |
1278 | fHistSignalinTPCKaonforITSPIDcuts->Fill(pt,(tpcSignal-sigKaon)/sigKaon); | |
1279 | else | |
1280 | fHistSignalinTPCAntiKaonforITSPIDcuts->Fill(TMath::Abs(pt),(tpcSignal-sigKaon)/sigKaon); | |
1281 | } | |
1282 | if(TMath::Abs(yforproton)<fYCut&&cutDCA) | |
1283 | { | |
1284 | fHistNtrackwithITSPIDcuts->Fill(TMath::Abs(pt),2); | |
1285 | if(pt>0.0) | |
1286 | fHistSignalinTPCProtonforITSPIDcuts->Fill(pt,(tpcSignal-sigProton)/sigProton); | |
1287 | else | |
1288 | fHistSignalinTPCAntiProtonforITSPIDcuts->Fill(TMath::Abs(pt),(tpcSignal-sigProton)/sigProton); | |
1289 | } | |
1290 | fHistStats->Fill(5); | |
1291 | fHistPhiPtAfterCuts->Fill(phi,pt); | |
1292 | fHistEtaPtAfterCuts->Fill(eta,pt); | |
1293 | fHistDCAAfterCuts->Fill(dcaxy,dcaz); | |
1294 | fHistPminusTPCinPAfterCuts->Fill(pinTPC-p,pinTPC); | |
1295 | fHistPminusTPCinPglobalAfterCuts->Fill(pinTPC-p,pinTPCglobal); | |
1296 | if(tmpQESD[0]>fchargeCut) | |
1297 | { | |
1298 | fHistL3dEP->Fill(p,tmpQESD[0]); | |
1299 | fHistL3dETPCinP->Fill(pinTPC,tmpQESD[0]); | |
1300 | } | |
1301 | if(tmpQESD[1]>fchargeCut) | |
1302 | { | |
1303 | fHistL4dEP->Fill(p,tmpQESD[1]); | |
1304 | fHistL4dETPCinP->Fill(pinTPC,tmpQESD[1]); | |
1305 | } | |
1306 | if(tmpQESD[2]>fchargeCut) | |
1307 | { | |
1308 | fHistL5dEP->Fill(p,tmpQESD[2]); | |
1309 | fHistL5dETPCinP->Fill(pinTPC,tmpQESD[2]); | |
1310 | } | |
1311 | if(tmpQESD[3]>fchargeCut) | |
1312 | { | |
1313 | fHistL6dEP->Fill(p,tmpQESD[3]); | |
1314 | fHistL6dETPCinP->Fill(pinTPC,tmpQESD[3]); | |
1315 | } | |
1316 | Float_t myITSsignal=0.0; | |
1317 | Float_t minITSsignal=0.0; | |
1318 | Int_t whichLmin=-1; | |
1319 | Int_t nosignaL=-1; | |
1320 | if(nSSDSDD==3) | |
1321 | { | |
1322 | Double_t tmp2QESD[3]={-1.0,-1.0,-1.0}; | |
1323 | Int_t iLnotZero=0; | |
1324 | for (int iL=0;iL<4;iL++) | |
1325 | { | |
1326 | if(tmpQESD[iL]>fchargeCut) | |
1327 | { | |
1328 | tmp2QESD[iLnotZero]=tmpQESD[iL]; | |
1329 | iLnotZero++; | |
1330 | } | |
1331 | else | |
1332 | nosignaL=iL; | |
1333 | } | |
1334 | whichLmin=TMath::LocMin(3,tmp2QESD); | |
1335 | if(nosignaL>-1&&nosignaL<=whichLmin) | |
1336 | whichLmin++; | |
1337 | minITSsignal=TMath::MinElement(3,tmp2QESD); | |
1338 | myITSsignal=MyITSsignalusing3points(tmp2QESD); | |
1339 | } | |
1340 | if(nSSDSDD==4) | |
1341 | { | |
1342 | myITSsignal=MyITSsignalusing4points(tmpQESD); | |
1343 | whichLmin=TMath::LocMin(4,tmpQESD); | |
1344 | minITSsignal=TMath::MinElement(4,tmpQESD); | |
1345 | } | |
1346 | ||
1347 | if(whichLmin==0) | |
1348 | fHistwhichhasmin->Fill(0.0,tmpQESD[0]); | |
1349 | if(whichLmin==1) | |
1350 | fHistwhichhasmin->Fill(1.0,tmpQESD[1]); | |
1351 | if(whichLmin==2) | |
1352 | fHistwhichhasmin->Fill(2.0,tmpQESD[2]); | |
1353 | if(whichLmin==3) | |
1354 | fHistwhichhasmin->Fill(3.0,tmpQESD[3]); | |
1355 | if(pt>0.0) | |
1356 | { | |
1357 | fHistMydEPpositive->Fill(p,myITSsignal); | |
1358 | fHistMydETPCinPglobalpositive->Fill(pinTPCglobal,myITSsignal); | |
1359 | fHistMydETPCinPpositive->Fill(pinTPC,myITSsignal); | |
1360 | } | |
1361 | else | |
1362 | { | |
1363 | fHistMydEPnegative->Fill(p,myITSsignal); | |
1364 | fHistMydETPCinPglobalnegative->Fill(pinTPCglobal,myITSsignal); | |
1365 | fHistMydETPCinPnegative->Fill(pinTPC,myITSsignal); | |
1366 | } | |
1367 | Float_t signaltouse=myITSsignal; | |
1368 | ||
1369 | ||
1370 | Float_t itspidsignalforpions=TMath::Log(signaltouse)-TMath::Log(fESDpid->GetITSResponse().Bethe(p,pionmass,kFALSE)); | |
1371 | Float_t itspidsignalforkaons=TMath::Log(signaltouse)-TMath::Log(fESDpid->GetITSResponse().Bethe(p,kaonmass,kFALSE)); | |
1372 | Float_t itspidsignalforprotons=TMath::Log(signaltouse)-TMath::Log(fESDpid->GetITSResponse().Bethe(p,protonmass,kFALSE)); | |
1373 | if(cutDCA) | |
1374 | { | |
1375 | mynumberoftracks++; | |
1376 | if(nSSDSDD==4) | |
1377 | fHistMysignalminusESD->Fill((signaltouse-trackESD->GetITSsignal())/signaltouse); | |
1378 | } | |
1379 | //Printf("Select on clean \n"); | |
1380 | if(TMath::Abs(yforpion)<=fYCut) | |
1381 | { | |
1382 | Float_t weight=1.0; | |
1383 | if(fMC) | |
1384 | weight=GetWeight(label,stack); | |
1385 | if(pt>0.0) | |
1386 | { | |
1387 | if(cutDCA) | |
1388 | { | |
1389 | if(fMC) | |
1390 | fHistminsignalifPionP->Fill(pt,itspidsignalforpions,weight); | |
1391 | else | |
1392 | fHistminsignalifPionP->Fill(pt,itspidsignalforpions); | |
1393 | if(itspidsignalforpions>-0.5&&itspidsignalforpions<0.2) //select on clean | |
1394 | { | |
1395 | fDCAXYZforcleanPions->Fill(pt,dcaxy,dcaz); | |
1396 | if(fMC) | |
1397 | { | |
1398 | if(primary&&pdgcode==211) | |
1399 | fDCAXYZforcleanPionsMCPrimary->Fill(pt,dcaxy,dcaz); | |
1400 | else if(!primary&&pdgcode==211&&uniqueID==kPDecay) | |
1401 | fDCAXYZforcleanPionsWD->Fill(pt,dcaxy,dcaz); | |
1402 | else if(!primary&&pdgcode==211&&uniqueID==kPHadronic) | |
1403 | fDCAXYZforcleanPionsHI->Fill(pt,dcaxy,dcaz); | |
1404 | else if(primary&&(pdgcode==-11||pdgcode==-13)) | |
1405 | fDCAXYZforcleanPionsMEPrimary->Fill(pt,dcaxy,dcaz); | |
1406 | else if(!primary&&(pdgcode==-11||pdgcode==-13)) | |
1407 | fDCAXYZforcleanPionsMESecondary->Fill(pt,dcaxy,dcaz); | |
1408 | else | |
1409 | fDCAXYZforcleanPionsR->Fill(pt,dcaxy,dcaz); | |
1410 | ||
1411 | } | |
1412 | } //select on clean | |
1413 | if(primary) | |
1414 | { | |
1415 | if(pdgcode==211) | |
1416 | fHistminsignalifPionPPrimary->Fill(pt,itspidsignalforpions); | |
1417 | } | |
1418 | else | |
1419 | { | |
1420 | if(pdgcode==211) | |
1421 | fHistminsignalifPionPSecondary->Fill(pt,itspidsignalforpions,weight); | |
1422 | //cout<<pdgcode<<" "<< uniqueID<<" "<<kPDecay<<" "<<kPHadronic<<endl; | |
1423 | } | |
1424 | if(pdgcode==-11||pdgcode==-13) | |
1425 | { | |
1426 | fHistminsignalifMuEPositiveP->Fill(pt,itspidsignalforpions); | |
1427 | if(!primary) | |
1428 | { | |
1429 | if(pdgcode==-11) | |
1430 | { | |
1431 | fHistStats->Fill(6); | |
1432 | fAntiElectronsource->Fill(pt,uniqueID); | |
1433 | } | |
1434 | else if(pdgcode==-13) | |
1435 | { | |
1436 | fHistStats->Fill(7); | |
1437 | fAntiMuonsource->Fill(pt,uniqueID); | |
1438 | } | |
1439 | } | |
1440 | else if(primary&&pdgcode==-11) | |
1441 | { | |
1442 | Printf("%d Mom",particle2->GetFirstMother()); | |
1443 | if(particle2->GetFirstMother()>-1) | |
1444 | fPrimaryElectronsMother->Fill(TMath::Abs(stack->Particle(particle2->GetFirstMother())->GetPdgCode())); | |
1445 | else | |
1446 | fPrimaryElectronsMother->Fill(-1); | |
1447 | fAntiElectronsource->Fill(pt,0); | |
1448 | } | |
1449 | else if(primary&&pdgcode==-13) | |
1450 | fAntiMuonsource->Fill(pt,0); | |
1451 | } | |
1452 | if(pdgcodefake==211) | |
1453 | { | |
1454 | if(primaryfake) | |
1455 | fHistminsignalifPionPrimaryfake->Fill(pt,itspidsignalforpions); | |
1456 | else | |
1457 | fHistminsignalifPionSecondaryfake->Fill(pt,itspidsignalforpions); | |
1458 | } | |
1459 | } | |
1460 | if(itspidsignalforpions>-0.5&&itspidsignalforpions<0.2) //select on clean | |
1461 | { | |
1462 | fDCAXYZOpenforcleanPions->Fill(pt,dcaxy,dcaz); | |
1463 | if(fMC) | |
1464 | { | |
1465 | if(primary&&pdgcode==211) | |
1466 | fDCAXYZOpenforcleanPionsMCPrimary->Fill(pt,dcaxy,dcaz); | |
1467 | else if(!primary&&pdgcode==211&&uniqueID==kPDecay) | |
1468 | fDCAXYZOpenforcleanPionsWD->Fill(pt,dcaxy,dcaz); | |
1469 | else if(!primary&&pdgcode==211&&uniqueID==kPHadronic) | |
1470 | fDCAXYZOpenforcleanPionsHI->Fill(pt,dcaxy,dcaz); | |
1471 | else if(primary&&(pdgcode==-11||pdgcode==-13)) | |
1472 | fDCAXYZOpenforcleanPionsMEPrimary->Fill(pt,dcaxy,dcaz); | |
1473 | else if(!primary&&(pdgcode==-11||pdgcode==-13)) | |
1474 | fDCAXYZOpenforcleanPionsMESecondary->Fill(pt,dcaxy,dcaz); | |
1475 | else | |
1476 | fDCAXYZOpenforcleanPionsR->Fill(pt,dcaxy,dcaz); | |
1477 | } | |
1478 | } | |
1479 | } | |
1480 | else | |
1481 | { | |
1482 | if(cutDCA) | |
1483 | { | |
1484 | if(fMC) | |
1485 | fHistminsignalifAntiPionP->Fill(TMath::Abs(pt),itspidsignalforpions,weight); | |
1486 | else | |
1487 | fHistminsignalifAntiPionP->Fill(TMath::Abs(pt),itspidsignalforpions); | |
1488 | if(itspidsignalforpions>-0.5&&itspidsignalforpions<0.2)//select on clean | |
1489 | { | |
1490 | fDCAXYZforcleanAntiPions->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1491 | if(fMC) | |
1492 | { | |
1493 | if(primary&&pdgcode==-211) | |
1494 | fDCAXYZforcleanAntiPionsMCPrimary->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1495 | else if(!primary&&pdgcode==-211&&uniqueID==kPDecay) | |
1496 | fDCAXYZforcleanAntiPionsWD->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1497 | else if(!primary&&pdgcode==-211&&uniqueID==kPHadronic) | |
1498 | fDCAXYZforcleanAntiPionsHI->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1499 | else if(primary&&(pdgcode==11||pdgcode==13)) | |
1500 | fDCAXYZforcleanAntiPionsMEPrimary->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1501 | else if(!primary&&(pdgcode==11||pdgcode==13)) | |
1502 | fDCAXYZforcleanAntiPionsMESecondary->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1503 | else | |
1504 | fDCAXYZforcleanAntiPionsR->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1505 | } | |
1506 | } //select on clean | |
1507 | if(primary) | |
1508 | { | |
1509 | if(pdgcode==-211) | |
1510 | fHistminsignalifAntiPionPPrimary->Fill(TMath::Abs(pt),itspidsignalforpions); | |
1511 | } | |
1512 | else | |
1513 | { | |
1514 | if(pdgcode==-211) | |
1515 | fHistminsignalifAntiPionPSecondary->Fill(TMath::Abs(pt),itspidsignalforpions,weight); | |
1516 | //cout<<pdgcode<<" "<< uniqueID<<" "<<kPDecay<<" "<<kPHadronic<<endl; | |
1517 | } | |
1518 | if(pdgcode==11||pdgcode==13) | |
1519 | { | |
1520 | fHistminsignalifMuENegativeP->Fill(TMath::Abs(pt),itspidsignalforpions); | |
1521 | if(!primary) | |
1522 | { | |
1523 | if(pdgcode==11) | |
1524 | { | |
1525 | fHistStats->Fill(6); | |
1526 | fElectronsource->Fill(TMath::Abs(pt),uniqueID); | |
1527 | } | |
1528 | else if(pdgcode==13) | |
1529 | { | |
1530 | fHistStats->Fill(7); | |
1531 | fMuonsource->Fill(TMath::Abs(pt),uniqueID); | |
1532 | } | |
1533 | } | |
1534 | else if(primary&&pdgcode==11) | |
1535 | { | |
1536 | Printf("%d Mom",particle2->GetFirstMother()); | |
1537 | if(particle2->GetFirstMother()>-1) | |
1538 | fPrimaryElectronsMother->Fill(TMath::Abs(stack->Particle(particle2->GetFirstMother())->GetPdgCode())); | |
1539 | else | |
1540 | fPrimaryElectronsMother->Fill(-1); | |
1541 | fElectronsource->Fill(TMath::Abs(pt),0); | |
1542 | } | |
1543 | else if(primary&&pdgcode==13) | |
1544 | fMuonsource->Fill(TMath::Abs(pt),0); | |
1545 | } | |
1546 | if(pdgcodefake==-211) | |
1547 | { | |
1548 | if(primaryfake) | |
1549 | fHistminsignalifAntiPionPrimaryfake->Fill(TMath::Abs(pt),itspidsignalforpions); | |
1550 | else | |
1551 | fHistminsignalifAntiPionSecondaryfake->Fill(TMath::Abs(pt),itspidsignalforpions); | |
1552 | } | |
1553 | } | |
1554 | if(itspidsignalforpions>-0.5&&itspidsignalforpions<0.2)//select on clean | |
1555 | { | |
1556 | fDCAXYZOpenforcleanAntiPions->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1557 | if(fMC) | |
1558 | { | |
1559 | if(primary&&pdgcode==-211) | |
1560 | fDCAXYZOpenforcleanAntiPionsMCPrimary->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1561 | else if(!primary&&pdgcode==-211&&uniqueID==kPDecay) | |
1562 | fDCAXYZOpenforcleanAntiPionsWD->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1563 | else if(!primary&&pdgcode==-211&&uniqueID==kPHadronic) | |
1564 | fDCAXYZOpenforcleanAntiPionsHI->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1565 | else if(primary&&(pdgcode==11||pdgcode==13)) | |
1566 | fDCAXYZOpenforcleanAntiPionsMEPrimary->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1567 | else if(!primary&&(pdgcode==11||pdgcode==13)) | |
1568 | fDCAXYZOpenforcleanAntiPionsMESecondary->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1569 | else | |
1570 | fDCAXYZOpenforcleanAntiPionsR->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1571 | } | |
1572 | } //sel | |
1573 | } | |
1574 | } | |
1575 | if(TMath::Abs(yforkaon)<=fYCut&&cutDCA) | |
1576 | { | |
1577 | if(pt>0.0) | |
1578 | { | |
1579 | fHistminsignalifKaonP->Fill(pt,itspidsignalforkaons); | |
1580 | if((itspidsignalforkaons>-0.2)&&(itspidsignalforkaons<0.2)) | |
1581 | { | |
1582 | fKaonNTPCClusters->Fill(pt,nTPCclusters); | |
1583 | fKaonchi2->Fill(pt,chi2); | |
1584 | } | |
1585 | if(primary) | |
1586 | { | |
1587 | if(pdgcode==321) | |
1588 | fHistminsignalifKaonPPrimary->Fill(pt,itspidsignalforkaons); | |
1589 | } | |
1590 | else | |
1591 | { | |
1592 | if(pdgcode==321) | |
1593 | fHistminsignalifKaonPSecondary->Fill(pt,itspidsignalforkaons); | |
1594 | } | |
1595 | if(pdgcodefake==321) | |
1596 | { | |
1597 | if(primaryfake) | |
1598 | fHistminsignalifKaonPrimaryfake->Fill(pt,itspidsignalforkaons); | |
1599 | else | |
1600 | fHistminsignalifKaonSecondaryfake->Fill(pt,itspidsignalforkaons); | |
1601 | } | |
1602 | ||
1603 | ||
1604 | } | |
1605 | else | |
1606 | { | |
1607 | fHistminsignalifAntiKaonP->Fill(TMath::Abs(pt),itspidsignalforkaons); | |
1608 | if((itspidsignalforkaons>-0.2)&&(itspidsignalforkaons<0.2)) | |
1609 | { | |
1610 | fAntiKaonNTPCClusters->Fill(TMath::Abs(pt),nTPCclusters); | |
1611 | fAntiKaonchi2->Fill(TMath::Abs(pt),chi2); | |
1612 | } | |
1613 | if(primary) | |
1614 | { | |
1615 | if(pdgcode==-321) | |
1616 | fHistminsignalifAntiKaonPPrimary->Fill(TMath::Abs(pt),itspidsignalforkaons); | |
1617 | } | |
1618 | else | |
1619 | { | |
1620 | if(pdgcode==-321) | |
1621 | fHistminsignalifAntiKaonPSecondary->Fill(TMath::Abs(pt),itspidsignalforkaons); | |
1622 | } | |
1623 | if(pdgcodefake==-321) | |
1624 | { | |
1625 | if(primaryfake) | |
1626 | fHistminsignalifAntiKaonPrimaryfake->Fill(TMath::Abs(pt),itspidsignalforkaons); | |
1627 | else | |
1628 | fHistminsignalifAntiKaonSecondaryfake->Fill(TMath::Abs(pt),itspidsignalforkaons); | |
1629 | ||
1630 | } | |
1631 | } | |
1632 | } | |
1633 | if(TMath::Abs(yforproton)<=fYCut) | |
1634 | { | |
1635 | Float_t weight=1.0; | |
1636 | if(fMC) | |
1637 | weight=GetWeight(label,stack); | |
1638 | if(pt>0.0) | |
1639 | { | |
1640 | if(cutDCA) | |
1641 | { | |
1642 | if(fMC) | |
1643 | fHistminsignalifProtonP->Fill(pt,itspidsignalforprotons,weight); | |
1644 | else | |
1645 | fHistminsignalifProtonP->Fill(pt,itspidsignalforprotons); | |
1646 | ||
1647 | //if(((itspidsignalforprotons))>(TMath::Abs(pt)<0.45?-0.2:0.0))&&(itspidsignalforprotons))<0.2)//select on | |
1648 | // clean | |
1649 | if((itspidsignalforprotons>-0.2)&&(itspidsignalforprotons<0.5)) | |
1650 | { | |
1651 | fDCAXYZforcleanProtons->Fill(pt,dcaxy,dcaz); | |
1652 | fProtonNTPCClusters->Fill(pt,nTPCclusters); | |
1653 | fProtonchi2->Fill(pt,chi2); | |
1654 | if(fMC) | |
1655 | { | |
1656 | if(primary&&pdgcode==2212) | |
1657 | fDCAXYZforcleanProtonsMCPrimary->Fill(pt,dcaxy,dcaz); | |
1658 | else if(!primary&&!primaryfake&&(pdgcode==2212||pdgcodefake==2212)&&uniqueID==kPDecay) | |
1659 | fDCAXYZforcleanProtonsWD->Fill(pt,dcaxy,dcaz); | |
1660 | else if(!primary&&!primaryfake&&(pdgcode==2212||pdgcodefake==2212)&&uniqueID==kPHadronic) | |
1661 | fDCAXYZforcleanProtonsHI->Fill(pt,dcaxy,dcaz); | |
1662 | else | |
1663 | fDCAXYZforcleanProtonsR->Fill(pt,dcaxy,dcaz); | |
1664 | } | |
1665 | }//select on clean | |
1666 | if(primary) | |
1667 | { | |
1668 | if(pdgcode==2212) | |
1669 | fHistminsignalifProtonPPrimary->Fill(pt,itspidsignalforprotons); | |
1670 | } | |
1671 | else if(primaryfake) | |
1672 | { | |
1673 | if(pdgcodefake==2212) | |
1674 | fHistminsignalifProtonPPrimaryfake->Fill(pt,itspidsignalforprotons); | |
1675 | } | |
1676 | else | |
1677 | { | |
1678 | if(pdgcode==2212&&uniqueID==kPDecay) | |
1679 | fHistminsignalifProtonPSecondaryWD->Fill(pt,itspidsignalforprotons,weight); | |
1680 | else if(pdgcode==2212&&uniqueID==kPHadronic) | |
1681 | fHistminsignalifProtonPSecondaryHI->Fill(pt,itspidsignalforprotons); | |
1682 | else if(pdgcodefake==2212&&uniqueID==kPDecay) | |
1683 | fHistminsignalifProtonPSecondaryWDfake->Fill(pt,itspidsignalforprotons,weight); | |
1684 | else if(pdgcodefake==2212&&uniqueID==kPHadronic) | |
1685 | fHistminsignalifProtonPSecondaryHIfake->Fill(pt,itspidsignalforprotons); | |
1686 | else if(fMC) | |
1687 | fHistminsignalifProtonPSecondaryRest->Fill(pt,itspidsignalforprotons); | |
1688 | } | |
1689 | } | |
1690 | ||
1691 | if((itspidsignalforprotons>-0.2)&&(itspidsignalforprotons<0.5)) | |
1692 | { | |
1693 | fDCAXYZOpenforcleanProtons->Fill(pt,dcaxy,dcaz); | |
1694 | if(fMC) | |
1695 | { | |
1696 | if(primary&&pdgcode==2212) | |
1697 | fDCAXYZOpenforcleanProtonsMCPrimary->Fill(pt,dcaxy,dcaz); | |
1698 | else if(!primary&&!primaryfake&&(pdgcode==2212||pdgcodefake==2212)&&uniqueID==kPDecay) | |
1699 | fDCAXYZOpenforcleanProtonsWD->Fill(pt,dcaxy,dcaz); | |
1700 | else if(!primary&&!primaryfake&&(pdgcode==2212||pdgcodefake==2212)&&uniqueID==kPHadronic) | |
1701 | fDCAXYZOpenforcleanProtonsHI->Fill(pt,dcaxy,dcaz); | |
1702 | else | |
1703 | fDCAXYZOpenforcleanProtonsR->Fill(pt,dcaxy,dcaz); | |
1704 | } | |
1705 | }//select on clean | |
1706 | ||
1707 | } | |
1708 | else | |
1709 | { | |
1710 | if(cutDCA) | |
1711 | { | |
1712 | if(fMC) | |
1713 | fHistminsignalifAntiProtonP->Fill(TMath::Abs(pt),itspidsignalforprotons,weight); | |
1714 | else | |
1715 | fHistminsignalifAntiProtonP->Fill(TMath::Abs(pt),itspidsignalforprotons); | |
1716 | if((itspidsignalforprotons>-0.2)&&(itspidsignalforprotons<0.5)) | |
1717 | {//select on clean | |
1718 | fDCAXYZforcleanAntiProtons->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1719 | ||
1720 | if(fMC) | |
1721 | { | |
1722 | if(primary&&pdgcode==-2212) | |
1723 | fDCAXYZforcleanAntiProtonsMCPrimary->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1724 | else if(!primary&&!primaryfake&&(pdgcode==-2212||pdgcodefake==-2212)&&uniqueID==kPDecay) | |
1725 | fDCAXYZforcleanAntiProtonsWD->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1726 | else if(!primary&&!primaryfake&&(pdgcode==-2212||pdgcodefake==-2212)&&uniqueID==kPHadronic) | |
1727 | fDCAXYZforcleanAntiProtonsHI->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1728 | else if(fMC) | |
1729 | fDCAXYZforcleanAntiProtonsR->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1730 | } | |
1731 | }//select on clean | |
1732 | if(primary) | |
1733 | { | |
1734 | if(pdgcode==-2212) | |
1735 | fHistminsignalifAntiProtonPPrimary->Fill(TMath::Abs(pt),itspidsignalforprotons); | |
1736 | } | |
1737 | else if(primaryfake) | |
1738 | { | |
1739 | if(pdgcodefake==-2212) | |
1740 | fHistminsignalifAntiProtonPPrimaryfake->Fill(TMath::Abs(pt),itspidsignalforprotons); | |
1741 | } | |
1742 | else | |
1743 | { | |
1744 | if(pdgcode==-2212&&uniqueID==kPDecay) | |
1745 | fHistminsignalifAntiProtonPSecondaryWD->Fill(TMath::Abs(pt),itspidsignalforprotons,weight); | |
1746 | else if(pdgcode==-2212&&uniqueID==kPHadronic) | |
1747 | fHistminsignalifAntiProtonPSecondaryHI->Fill(TMath::Abs(pt),itspidsignalforprotons); | |
1748 | else if(pdgcodefake==-2212&&uniqueID==kPDecay) | |
1749 | fHistminsignalifAntiProtonPSecondaryWDfake->Fill(TMath::Abs(pt),itspidsignalforprotons,weight); | |
1750 | else if(pdgcodefake==-2212&&uniqueID==kPHadronic) | |
1751 | fHistminsignalifAntiProtonPSecondaryHIfake->Fill(TMath::Abs(pt),itspidsignalforprotons); | |
1752 | else if(fMC) | |
1753 | fHistminsignalifAntiProtonPSecondaryRest->Fill(TMath::Abs(pt),itspidsignalforprotons); | |
1754 | } | |
1755 | } | |
1756 | if((itspidsignalforprotons>-0.2)&&(itspidsignalforprotons<0.5)) | |
1757 | {//select on clean | |
1758 | fDCAXYZOpenforcleanAntiProtons->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1759 | fAntiProtonNTPCClusters->Fill(TMath::Abs(pt),nTPCclusters); | |
1760 | fAntiProtonchi2->Fill(TMath::Abs(pt),chi2); | |
1761 | if(fMC) | |
1762 | { | |
1763 | if(primary&&pdgcode==-2212) | |
1764 | fDCAXYZOpenforcleanAntiProtonsMCPrimary->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1765 | else if(!primary&&!primaryfake&&(pdgcode==-2212||pdgcodefake==-2212)&&uniqueID==kPDecay) | |
1766 | fDCAXYZOpenforcleanAntiProtonsWD->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1767 | else if(!primary&&!primaryfake&&(pdgcode==-2212||pdgcodefake==-2212)&&uniqueID==kPHadronic) | |
1768 | fDCAXYZOpenforcleanAntiProtonsHI->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1769 | else if(fMC) | |
1770 | fDCAXYZOpenforcleanAntiProtonsR->Fill(TMath::Abs(pt),dcaxy,dcaz); | |
1771 | } | |
1772 | } | |
1773 | } | |
1774 | } | |
1775 | } | |
1776 | fHistStandartMul->Fill(refmultiplicity); | |
1777 | fHistMytrackMul->Fill(mynumberoftracks); | |
1778 | ||
1779 | ||
1780 | // Post output data. | |
1781 | Printf("Done..........\n"); | |
1782 | PostData(1, flist); | |
1783 | //Printf("....................Done!\n"); | |
1784 | //Printf("end of Exec"); | |
1785 | } | |
1786 | ||
1787 | //________________________________________________________________________ | |
1788 | void AliAnalysisChargedHadronSpectraITSTruncatedMeanTask::Terminate(Option_t *) | |
1789 | { | |
1790 | //Terminate | |
1791 | if(fCuts) | |
1792 | fCuts->Dump(); | |
1793 | Printf("YCut=%f",fYCut); | |
1794 | Printf("nsigma=%f",fsigmacut); | |
1795 | Printf("DCA cut xy sigma =%f zsigma=%f", fnsigmaxy, fnsigmaz); | |
1796 | Printf("ChargeCut=%f ", fchargeCut); | |
1797 | Printf("DCAxy parameters %f %f %f",fdcaxypar[0],fdcaxypar[1],fdcaxypar[2]); | |
1798 | Printf("DCAz parameters %f %f %f %f",fdcazpar[0],fdcazpar[1],fdcazpar[2],fdcazpar[3]); | |
1799 | if(fTPCPIDCUT) | |
1800 | fTPCPIDCUT->Dump(); | |
1801 | if(fMC) | |
1802 | Printf("MC On\n"); | |
1803 | if(fCorrectSDD) | |
1804 | Printf("correct SDD On\n"); | |
1805 | if(fCorrectSSD) | |
1806 | Printf("correct SSD On\n"); | |
1807 | if(fK0weight) | |
1808 | { | |
1809 | Printf("weigth for pions"); | |
1810 | fK0weight->Print("All"); | |
1811 | } | |
1812 | if(flambdaweight) | |
1813 | { | |
1814 | Printf("weigth for protons"); | |
1815 | flambdaweight->Print("All"); | |
1816 | } | |
1817 | if(fAntilambdaweight) | |
1818 | { | |
1819 | Printf("weigth for antiprotons"); | |
1820 | fAntilambdaweight->Print("All"); | |
1821 | } | |
1822 | Printf("Mul low %d Mul up %d",fLowMultiplicity, fUpMultiplicity); | |
1823 | Printf("cent low %f cent up %f",fLowCentrality,fUpCentrality); | |
1824 | if(fdovertexrescuts) | |
1825 | Printf("Veretx resolution cut"); | |
1826 | Printf("end of Terminate"); | |
1827 | } | |
1828 | //___________________________________________________ | |
1829 | Float_t AliAnalysisChargedHadronSpectraITSTruncatedMeanTask::MyITSsignalusing4points(Double_t* tmpQESD) const | |
1830 | { | |
1831 | //dE signal in case of 4 points in ITS | |
1832 | Int_t indexes[4]={-1,-1,-1,-1}; | |
1833 | TMath::Sort(4,tmpQESD,indexes,0); | |
1834 | return 0.5*(tmpQESD[indexes[0]]+tmpQESD[indexes[1]]); | |
1835 | } | |
1836 | //________________________________________________________ | |
1837 | Float_t AliAnalysisChargedHadronSpectraITSTruncatedMeanTask::MyITSsignalusing3points( Double_t* tmpQESD) const | |
1838 | { | |
1839 | //dE signal in case of 3 points in ITS | |
1840 | Int_t indexes[3]={-1,-1,-1}; | |
1841 | TMath::Sort(3,tmpQESD,indexes,0); | |
1842 | //cout<<tmpQESD[indexes[0]]<<" "<<tmpQESD[indexes[1]]<<" "<<tmpQESD[indexes[2]]<<endl; | |
1843 | return (tmpQESD[indexes[0]]+tmpQESD[indexes[1]]*0.5)/1.5; | |
1844 | } | |
1845 | //____________________________________________________________________________________________________ | |
1846 | void AliAnalysisChargedHadronSpectraITSTruncatedMeanTask::SetFunctionParam( Double_t * const par) | |
1847 | { | |
1848 | fESDpid->GetITSResponse().SetBetheBlochParamsITSTPC(par); | |
1849 | } | |
1850 | //_____________________________________________________________________________________________________ | |
1851 | void AliAnalysisChargedHadronSpectraITSTruncatedMeanTask::CorrectSDD(Double_t *tmpQESD) const | |
1852 | { | |
1853 | //correction of SDD signal | |
1854 | if(tmpQESD[0]>0.0) | |
1855 | { | |
1856 | tmpQESD[0]=tmpQESD[0]*3.34/5.43; | |
1857 | if(tmpQESD[0]<30.0) | |
1858 | tmpQESD[0]=-1.0; | |
1859 | } | |
1860 | if(tmpQESD[1]>0.0) | |
1861 | { | |
1862 | tmpQESD[1]=tmpQESD[1]*3.34/5.43; | |
1863 | if(tmpQESD[1]<30.0) | |
1864 | tmpQESD[1]=-1.0; | |
1865 | } | |
1866 | } | |
1867 | //_____________________________________________________________________________________________________ | |
1868 | void AliAnalysisChargedHadronSpectraITSTruncatedMeanTask::CorrectSSD(Double_t *tmpQESD) const | |
1869 | { | |
1870 | //Correction of SSD signal | |
1871 | tmpQESD[2]=(85.0/77.0)*tmpQESD[2]; | |
1872 | tmpQESD[3]=(85.0/77.0)*tmpQESD[3]; | |
1873 | } | |
1874 | //_______________________________________________________________________________________________________ | |
1875 | Bool_t AliAnalysisChargedHadronSpectraITSTruncatedMeanTask::SelectOnImpPar(AliESDtrack* t) const | |
1876 | { | |
1877 | // | |
1878 | // cut on transverse impact parameter | |
1879 | // | |
1880 | Float_t d0z0[2],covd0z0[3]; | |
1881 | t->GetImpactParameters(d0z0,covd0z0); | |
1882 | Float_t sigma= fdcaxypar[0]+fdcaxypar[1]/TMath::Power(t->Pt(),fdcaxypar[2]); | |
1883 | Float_t d0max = fnsigmaxy*sigma; | |
1884 | // | |
1885 | Float_t sigmaZ = fdcazpar[0]+fdcazpar[1]/TMath::Power(t->Pt(),fdcazpar[2]); | |
1886 | if (t->Pt() > 1) | |
1887 | sigmaZ = fdcazpar[3]; | |
1888 | Float_t d0maxZ = fnsigmaz*sigmaZ; | |
1889 | // | |
1890 | if(TMath::Abs(d0z0[0]) < d0max && TMath::Abs(d0z0[1]) < d0maxZ) //error | |
1891 | return kTRUE; | |
1892 | return kFALSE; | |
1893 | } | |
1894 | //__________________________________________________________________________________________________ | |
1895 | Float_t AliAnalysisChargedHadronSpectraITSTruncatedMeanTask::GetWeight(Int_t t,AliStack* const stack) const | |
1896 | { | |
1897 | ||
1898 | //Get weigth for pions protons and antiprotons | |
1899 | if(stack->IsPhysicalPrimary(TMath::Abs(t))) | |
1900 | return 1.0; | |
1901 | TParticle *particleMC = stack->Particle(TMath::Abs(t)); | |
1902 | Int_t pdgcodeMC = particleMC->GetPdgCode(); | |
1903 | if(TMath::Abs(pdgcodeMC)!=211&&TMath::Abs(pdgcodeMC)!=2212) | |
1904 | return 1.0; | |
1905 | if(!stack->IsPhysicalPrimary(TMath::Abs(particleMC->GetFirstMother()))) | |
1906 | return 1.0; | |
1907 | TParticle *particleMother=stack->Particle(TMath::Abs(particleMC->GetFirstMother())); | |
1908 | Int_t pdgcodeMother = particleMother->GetPdgCode(); | |
1909 | Float_t motherpt=particleMother-> Pt(); | |
1910 | if(TMath::Abs(pdgcodeMC)==211&&pdgcodeMother==310&&fK0weight) | |
1911 | return fK0weight->Eval(motherpt); | |
1912 | else if (pdgcodeMother==3122&&flambdaweight) | |
1913 | return flambdaweight->Eval(motherpt); | |
1914 | else if(pdgcodeMother==-3122&&fAntilambdaweight) | |
1915 | return fAntilambdaweight->Eval(motherpt); | |
1916 | return 1.0; | |
1917 | } | |
1918 | //________________________________________________________________________________________________ | |
1919 | void AliAnalysisChargedHadronSpectraITSTruncatedMeanTask::SetDCA2010() | |
1920 | { | |
1921 | //setting the DCA for 2010 | |
1922 | fdcaxypar[0]=0.0026; | |
1923 | fdcaxypar[1]=0.005; | |
1924 | fdcaxypar[2]=1.01; | |
1925 | ||
1926 | fdcazpar[0]=1000000.0; | |
1927 | fdcazpar[1]=0.0; | |
1928 | fdcazpar[2]=1.0; | |
1929 | fdcazpar[3]=1000000.0; | |
1930 | } | |
1931 | //______________________________________________________________________________________________________________ | |
1932 | void AliAnalysisChargedHadronSpectraITSTruncatedMeanTask::SetCentralityCut(Float_t low, Float_t up) | |
1933 | { | |
1934 | //centrality cut setter | |
1935 | if((up>low)&&(!(low<0.0))&&(!(up>100.0))) | |
1936 | { | |
1937 | SetHImode(); | |
1938 | fLowCentrality=low; | |
1939 | fUpCentrality=up; | |
1940 | } | |
1941 | } | |
1942 | //_____________________________________________________________________________________________________________ | |
1943 | void AliAnalysisChargedHadronSpectraITSTruncatedMeanTask::SetMultiplicityCut(Int_t low, Int_t up) | |
1944 | { | |
1945 | //mulyiplicty cut setter | |
1946 | if((!(up>low))&&low>=0&&up>=0) | |
1947 | { | |
1948 | fLowMultiplicity=-1; | |
1949 | fUpMultiplicity=-1; | |
1950 | } | |
1951 | else | |
1952 | { | |
1953 | fLowMultiplicity=low; | |
1954 | fUpMultiplicity=up; | |
1955 | } | |
1956 | } |