3 void MakeDataReport(const char* outputFile="JpsiDataReport.pdf",
4 const char* histos="jpsi_HistosSE.root",
5 const char* cf="jpsi_CF.root")
8 // Make a pdf file with the efficiency report
13 AliDielectronCFdraw d(cf);
14 d.SetRangeUser("PairType",1,1);
15 d.SetRangeUser("Y",-.89,.9,"0");
18 TFile f("jpsi_HistosSE.root");
20 AliDielectronHistos h;
21 TIter nextHists((TList*)f.Get("Dielectron_Histos"));
23 TPaveText pt(.02,.6,.98,.8);
24 TText *t1=pt.AddText("");
25 TText *t2=pt.AddText("");
27 TCanvas *c1=new TCanvas;
32 // Invariant mass plots
40 t1->SetTitle("QA summary plots for");
42 while ( (list=(THashList*)nextHists()) ){
43 h.SetHistogramList(*list);
44 t2->SetTitle(list->GetName());
56 const Int_t NCont=255;
58 TStyle *st = new TStyle("mystyle","mystyle");
59 gROOT->GetStyle("Plain")->Copy((*st));
65 st->SetNumberContours(NCont);
67 st->SetOptStat("erm");
69 st->SetGridColor(kGray+1);
70 st->SetPadGridX(kTRUE);
71 st->SetPadGridY(kTRUE);
72 st->SetPadTickX(kTRUE);
73 st->SetPadTickY(kTRUE);
76 const Int_t NRGBs = 5;
77 Double_t stops[NRGBs] = { 0.00, 0.34, 0.61, 0.84, 1.00 };
78 Double_t red[NRGBs] = { 0.00, 0.00, 0.87, 1.00, 0.51 };
79 Double_t green[NRGBs] = { 0.00, 0.81, 1.00, 0.20, 0.00 };
80 Double_t blue[NRGBs] = { 0.51, 1.00, 0.12, 0.00, 0.00 };
82 TColor::CreateGradientColorTable(NRGBs, stops, red, green, blue, NCont);
87 TFile f("jpsi_debug.root");
88 // if (!f.IsOpen()) return;
90 TTree *t=(TTree*)f.Get("Pair");
95 gStyle->SetOptStat(0);
97 TLegend *leg=new TLegend(0.59,.79,.99,.99);
100 l.SetLineColor(kGreen-5);
103 leg->SetFillColor(10);
108 t->SetLineColor(kBlack);
109 t->Draw("M>>hAll(200,-.01,3.99)","","histe");
110 TH1 *hAll=(TH1*)gROOT->FindObject("hAll");
111 hAll->SetMinimum(0.1);
112 hAll->SetTitle(";M [GeV]; yield");
113 leg->AddEntry(hAll,"|n#sigma e|<2 + pt>0.3 GeV","l");
115 l.DrawLine(3.097,1,3.097,1e4);
117 t->SetLineColor(kOrange-5);
118 t->Draw("M>>hC11(200,-.01,3.99)","abs(Leg1_ImpactParXY)<.004&&abs(Leg2_ImpactParXY)<.004","histesame");
119 hAll=(TH1*)gROOT->FindObject("hC11");
120 leg->AddEntry(hAll,"|n#sigma e|<2 + pt>0.3 GeV + |dXY|<40#mum","l");
122 TCut d1_1="abs(Leg1_TPC_nSigma_Electrons)<1";
123 TCut d2_1="abs(Leg2_TPC_nSigma_Electrons)<1";
126 t->SetLineColor(kRed);
127 t->Draw("M>>hC1(200,-.01,3.99)","Leg2_Pt>1&&Leg1_Pt>1","histesame");
128 hAll=(TH1*)gROOT->FindObject("hC1");
129 leg->AddEntry(hAll,"|n#sigma e|<2 + pt>1 GeV","l");
132 t->SetLineColor(kGreen);
133 t->Draw("M>>hC2(200,-.01,3.99)",d_1+"Leg2_Pt>1&&Leg1_Pt>1","histesame");
134 hAll=(TH1*)gROOT->FindObject("hC2");
135 leg->AddEntry(hAll,"|n#sigma e|<1 + pt>1 GeV","l");
137 t->SetLineColor(kMagenta);
138 t->Draw("M>>hC3(200,-.01,3.99)","Leg1_Pt>2&&Leg2_Pt>2","histesame");
139 hAll=(TH1*)gROOT->FindObject("hC3");
140 leg->AddEntry(hAll,"|n#sigma e|<2 + pt>2 GeV","l");
142 t->SetLineColor(kBlue);
143 t->Draw("M>>hC4(200,-.01,3.99)","Leg1_Pt>3&&Leg2_Pt>3","histesame");
144 hAll=(TH1*)gROOT->FindObject("hC4");
145 leg->AddEntry(hAll,"|n#sigma e|<2 + pt>3 GeV","l");
151 Double_t alephParameters[5];
153 alephParameters[0] = 2.15898e+00/50.;
154 alephParameters[1] = 1.75295e+01;
155 alephParameters[2] = 3.40030e-09;
156 alephParameters[3] = 1.96178e+00;
157 alephParameters[4] = 3.91720e+00;
160 TF1 *foProton = new TF1("foProton", "50*AliExternalTrackParam::BetheBlochAleph(x/0.93827,[0],[1],[2],[3],[4])",0.05,20);
161 TF1 *foPion = new TF1("foPion", "50*AliExternalTrackParam::BetheBlochAleph(x/0.13957,[0],[1],[2],[3],[4])",0.05,20);
162 TF1 *foElec = new TF1("foElec", "50*AliExternalTrackParam::BetheBlochAleph(x/0.000511,[0],[1],[2],[3],[4])",0.05,20);
163 TF1 *foKaon = new TF1("foKaon", "50*AliExternalTrackParam::BetheBlochAleph(x/0.493677,[0],[1],[2],[3],[4])",0.05,20);
164 TF1 *foMuon = new TF1("foMuon", "50*AliExternalTrackParam::BetheBlochAleph(x/0.105658,[0],[1],[2],[3],[4])",0.05,20);
166 foProton->SetParameters(alephParameters);
167 foPion->SetParameters(alephParameters);
168 foElec->SetParameters(alephParameters);
169 foKaon->SetParameters(alephParameters);
170 foMuon->SetParameters(alephParameters);
172 foProton->SetLineColor(color);
173 foPion->SetLineColor(color);
174 foElec->SetLineColor(color);
175 foKaon->SetLineColor(color);
176 foMuon->SetLineColor(color);
179 foProton->SetLineWidth(lineWidth);
180 foPion->SetLineWidth(lineWidth);
181 foElec->SetLineWidth(lineWidth);
182 foKaon->SetLineWidth(lineWidth);
183 foMuon->SetLineWidth(lineWidth);
186 foProton->SetNpx(200);
192 foProton->Draw("same");
193 foPion->Draw("same");
194 foElec->Draw("same");
195 foKaon->Draw("same");
196 foMuon->Draw("same");
204 alephParameters[0] = 0.0283086;
205 alephParameters[1] = 2.63394e+01;
206 alephParameters[2] = 5.04114e-11;
207 alephParameters[3] = 2.12543e+00;
208 alephParameters[4] = 4.88663e+00;
211 alephParameters[0] = 0.0283086/0.97;
212 //alephParameters[0] = 0.0283086;
213 alephParameters[1] = 2.63394e+01;
214 alephParameters[2] = 5.04114e-11;
215 alephParameters[3] = 2.12543e+00;
216 alephParameters[4] = 4.88663e+00;
219 TF1 *foDataProton = new TF1("foDataProton", "50*AliExternalTrackParam::BetheBlochAleph(x/0.93827,[0],[1],[2],[3],[4])",0.05,20);
220 TF1 *foDataProtonP = new TF1("foDataProtonP",Form( "50*AliExternalTrackParam::BetheBlochAleph(x/0.93827,[0],[1],[2],[3],[4])*(1+%f)",3*res),0.05,20);
221 TF1 *foDataProtonM = new TF1("foDataProtonM", Form("50*AliExternalTrackParam::BetheBlochAleph(x/0.93827,[0],[1],[2],[3],[4])*(1-%f)",res),0.05,20);
223 TF1 *foDataPion = new TF1("foDataPion", "50*AliExternalTrackParam::BetheBlochAleph(x/0.13957,[0],[1],[2],[3],[4])",0.05,20);
224 TF1 *foDataPionP = new TF1("foDataPionP",Form( "50*AliExternalTrackParam::BetheBlochAleph(x/0.13957,[0],[1],[2],[3],[4])*(1+%f)",res),0.05,20);
225 TF1 *foDataPionM = new TF1("foDataPionM", Form("50*AliExternalTrackParam::BetheBlochAleph(x/0.93827,[0],[1],[2],[3],[4])*(1-%f)",res),0.05,20);
227 TF1 *foDataElec = new TF1("foDataElec", "50*AliExternalTrackParam::BetheBlochAleph(x/0.000511,[0],[1],[2],[3],[4])",0.05,20);
228 TF1 *foDataElecP = new TF1("foDataElecP",Form( "50*AliExternalTrackParam::BetheBlochAleph(x/0.93827,[0],[1],[2],[3],[4])*(1+%f)",res),0.05,20);
229 TF1 *foDataElecM = new TF1("foDataElecM", Form("50*AliExternalTrackParam::BetheBlochAleph(x/0.93827,[0],[1],[2],[3],[4])*(1-%f)",res),0.05,20);
231 TF1 *foDataKaon = new TF1("foDataKaon", "50*AliExternalTrackParam::BetheBlochAleph(x/0.493677,[0],[1],[2],[3],[4])",0.05,20);
232 TF1 *foDataKaonP = new TF1("foDataKaonP",Form( "50*AliExternalTrackParam::BetheBlochAleph(x/0.93827,[0],[1],[2],[3],[4])*(1+%f)",res),0.05,20);
233 TF1 *foDataKaonM = new TF1("foDataKaonM", Form("50*AliExternalTrackParam::BetheBlochAleph(x/0.93827,[0],[1],[2],[3],[4])*(1-%f)",res),0.05,20);
235 TF1 *foDataMuon = new TF1("foDataMuon", "50*AliExternalTrackParam::BetheBlochAleph(x/0.105658,[0],[1],[2],[3],[4])",0.05,20);
236 TF1 *foDataMuonP = new TF1("foDataMuonP",Form( "50*AliExternalTrackParam::BetheBlochAleph(x/0.93827,[0],[1],[2],[3],[4])*(1+%f)",res),0.05,20);
237 TF1 *foDataMuonM = new TF1("foDataMuonM", Form("50*AliExternalTrackParam::BetheBlochAleph(x/0.93827,[0],[1],[2],[3],[4])*(1-%f)",res),0.05,20);
240 foDataProton->SetParameters(alephParameters);
241 foDataProtonP->SetParameters(alephParameters);
242 foDataProtonM->SetParameters(alephParameters);
243 foDataPion->SetParameters(alephParameters);
244 foDataPionP->SetParameters(alephParameters);
245 foDataPionM->SetParameters(alephParameters);
246 foDataElec->SetParameters(alephParameters);
247 foDataElecP->SetParameters(alephParameters);
248 foDataElecM->SetParameters(alephParameters);
249 foDataKaon->SetParameters(alephParameters);
250 foDataKaonP->SetParameters(alephParameters);
251 foDataKaonM->SetParameters(alephParameters);
252 foDataMuon->SetParameters(alephParameters);
253 foDataMuonP->SetParameters(alephParameters);
254 foDataMuonM->SetParameters(alephParameters);
256 foDataProton->SetLineColor(color);
257 foDataProtonP->SetLineColor(color-4);
258 foDataProtonM->SetLineColor(color-4);
259 foDataPion->SetLineColor(color);
260 foDataPionP->SetLineColor(color-4);
261 foDataPionM->SetLineColor(color-4);
262 foDataElec->SetLineColor(color);
263 foDataElecP->SetLineColor(color-4);
264 foDataElecM->SetLineColor(color-4);
265 foDataKaon->SetLineColor(color);
266 foDataKaonP->SetLineColor(color-4);
267 foDataKaonM->SetLineColor(color-4);
268 foDataMuon->SetLineColor(color);
269 foDataMuonP->SetLineColor(color-4);
270 foDataMuonM->SetLineColor(color-4);
273 foDataProton->SetLineWidth(lineWidth);
274 foDataProtonP->SetLineWidth(lineWidth);
275 foDataProtonM->SetLineWidth(lineWidth);
276 foDataPion->SetLineWidth(lineWidth);
277 foDataPionP->SetLineWidth(lineWidth);
278 foDataPionM->SetLineWidth(lineWidth);
279 foDataElec->SetLineWidth(lineWidth);
280 foDataElecP->SetLineWidth(lineWidth);
281 foDataElecM->SetLineWidth(lineWidth);
282 foDataKaon->SetLineWidth(lineWidth);
283 foDataKaonP->SetLineWidth(lineWidth);
284 foDataKaonM->SetLineWidth(lineWidth);
285 foDataMuon->SetLineWidth(lineWidth);
286 foDataMuonP->SetLineWidth(lineWidth);
287 foDataMuonM->SetLineWidth(lineWidth);
290 foDataProtonP->SetLineStyle(2);
291 foDataProtonM->SetLineStyle(2);
292 foDataPionP->SetLineStyle(2);
293 foDataPionM->SetLineStyle(2);
294 foDataElecP->SetLineStyle(2);
295 foDataElecM->SetLineStyle(2);
296 foDataKaonP->SetLineStyle(2);
297 foDataKaonM->SetLineStyle(2);
298 foDataMuonP->SetLineStyle(2);
299 foDataMuonM->SetLineStyle(2);
302 foDataProton->SetNpx(200);
303 foDataProtonP->SetNpx(200);
304 foDataProtonM->SetNpx(200);
305 foDataPion->SetNpx(200);
306 foDataPionP->SetNpx(200);
307 foDataPionM->SetNpx(200);
308 foDataElec->SetNpx(200);
309 foDataKaon->SetNpx(200);
310 foDataMuon->SetNpx(200);
312 foDataProton->Draw("same");
313 foDataProtonP->Draw("same");
314 foDataProtonM->Draw("same");
315 foDataPion->Draw("same");
316 foDataElec->Draw("same");
317 foDataKaon->Draw("same");
318 foDataMuon->Draw("same");
326 Int_t baseColors[5]={kRed, kGreen+1, kAzure-4, kMagenta, kCyan+1};
327 Int_t sigmaColorOffset=1;
329 Int_t baseColors[5]={kRed, kGreen+1, kAzure-4, kMagenta, kCyan+1};
330 Int_t sigmaColorOffset=0;
331 Int_t baseColors[5]={kRed, kRed, kRed, kRed, kRed};
333 Double_t sigmas[5]={3,3,3,3,3};
336 for (Int_t i=0; i<AliPID::kSPECIES; ++i) masses[i]=AliPID::ParticleMass(i);
339 Double_t alephParameters[5];
341 alephParameters[0] = 0.0283086;
342 alephParameters[1] = 2.63394e+01;
343 alephParameters[2] = 5.04114e-11;
344 alephParameters[3] = 2.12543e+00;
345 alephParameters[4] = 4.88663e+00;
346 alephParameters[0] = 1.25202/50.; //was 1.79571/55.;
347 alephParameters[1] = 2.74992e+01; //was 22.0028;
348 alephParameters[2] = TMath::Exp(-3.31517e+01); //was1.55354e-11;
349 alephParameters[3] = 2.46246; //was 2.39804;
350 alephParameters[4] = 6.78938; //was 5.1209;
361 for (Int_t i=0; i<5; ++i){
362 fBethe[i] = new TF1(Form("fBethe%d",i), Form("%f*AliExternalTrackParam::BetheBlochAleph(x/%f,[0],[1],[2],[3],[4])",mip,masses[i]),0.05,20);
363 fBetheP[i] = new TF1(Form("fBethe%d",i), Form("%f*AliExternalTrackParam::BetheBlochAleph(x/%f,[0],[1],[2],[3],[4])*(1+%f*%f)",mip,masses[i],res,sigmas[i]),0.05,20);
364 fBetheM[i] = new TF1(Form("fBethe%d",i), Form("%f*AliExternalTrackParam::BetheBlochAleph(x/%f,[0],[1],[2],[3],[4])*(1-%f*%f)",mip,masses[i],res,sigmas[i]),0.05,20);
366 fBethe[i]->SetParameters(alephParameters);
367 fBetheP[i]->SetParameters(alephParameters);
368 fBetheM[i]->SetParameters(alephParameters);
370 fBethe[i]->SetLineColor(baseColors[i]);
371 fBetheP[i]->SetLineColor(baseColors[i]-sigmaColorOffset);
372 fBetheM[i]->SetLineColor(baseColors[i]-sigmaColorOffset);
374 fBethe[i]->SetLineWidth(lineWidth);
375 fBetheP[i]->SetLineWidth(lineWidth);
376 fBetheM[i]->SetLineWidth(lineWidth);
378 fBetheP[i]->SetLineStyle(2);
379 fBetheM[i]->SetLineStyle(2);
381 fBethe[i]->SetNpx(200);
382 fBetheP[i]->SetNpx(200);
383 fBetheM[i]->SetNpx(200);
386 for (Int_t i=0; i<5; ++i){
387 fBethe[i]->Draw("same");
388 // fBetheP[i]->Draw("same");
389 // fBetheM[i]->Draw("same");
398 Double_t resolution=0.052;
400 TF1 *ffPio=new TF1(Form("fBethe%d",AliPID::kPion), Form("(%f*%f+(AliExternalTrackParam::BetheBlochAleph(x/%f,[0],[1],[2],[3],[4])-AliExternalTrackParam::BetheBlochAleph(x/%f,[0],[1],[2],[3],[4])))/%f", nSigma,resolution, AliPID::ParticleMass(AliPID::kPion), AliPID::ParticleMass(AliPID::kElectron), resolution), 0.05,200.);
401 ffPio->SetParameters(0.0283086/0.97,2.63394e+01,5.04114e-11,2.12543e+00,4.88663e+00);
404 TF1 f("fP","-8*exp(-0.6*x)",0,40);
414 c->SetAlias("cut","PairType==1")
417 c->SetAlias("nCls","Leg1_NclsTPC>90&&Leg2_NclsTPC>90");
422 c->SetAlias("cutE","abs(Leg1_TPC_nSigma_Electrons)<3&&abs(Leg2_TPC_nSigma_Electrons)<3");
423 // c->SetAlias("cutE","Leg1_TPC_nSigma_Electrons>-1 && Leg2_TPC_nSigma_Electrons>-1");
424 c->SetAlias("cutPi","abs(Leg1_TPC_nSigma_Pions)>3&&abs(Leg2_TPC_nSigma_Pions)>3");
425 c->SetAlias("cutP","(Leg1_TPC_nSigma_Protons)>3&&(Leg2_TPC_nSigma_Protons)>3");
426 c->SetAlias("pidSig","cutE&&cutPi&&cutP");
429 // c->SetAlias("eleParam","Leg1_TPC_nSigma_Electrons<5&&Leg2_TPC_nSigma_Electrons<5&&Leg1_TPC_nSigma_Electrons>-2.65*exp(-0.6757*Leg1_P_InnerParam)&&Leg2_TPC_nSigma_Electrons>-8*exp(-0.6*Leg2_P_InnerParam)");
430 c->SetAlias("eleParam","Leg1_TPC_nSigma_Electrons<5&&Leg2_TPC_nSigma_Electrons<5&&Leg1_TPC_nSigma_Electrons>-6*exp(-0.6*Leg1_P_InnerParam)&&Leg2_TPC_nSigma_Electrons>-6*exp(-0.6*Leg2_P_InnerParam)");
431 c->SetAlias("pidParam","eleParam&&cutP");
435 c->SetAlias("LegEta","abs(Leg1_Eta)<0.9&&abs(Leg2_Eta<0.9)");
436 c->SetAlias("LegNcl","Leg1_NclsTPC>90&&Leg2_NclsTPC>90");
437 c->SetAlias("Rap","abs(Y)<0.9");
438 c->SetAlias("QA","LegNcl&&LegEta&&Rap");
439 c->SetAlias("spdFirst","(Leg1_ITS_clusterMap&1)==1 && (Leg2_ITS_clusterMap&1)==1");
440 c->SetAlias("LegNclDiffIter1","abs(Leg1_NclsTPC-Leg1_NclsTPCiter1)<10&&abs(Leg2_NclsTPC-Leg2_NclsTPCiter1)<10")
441 c->SetAlias("LegNclPID","(Leg1_NclsTPC-Leg1_TPCsignalN)<20&&(Leg2_NclsTPC-Leg2_TPCsignalN)<20")
443 c->SetAlias("cut","PairType==1&&QA&&pidSig")
445 c->SetMarkerStyle(20);
446 c->SetMarkerSize(.8);
447 c->SetMarkerColor(kBlack);
448 c->SetLineColor(kBlack);
450 // c->SetAlias("nCls","Leg1_NclsTPC>90&&Leg2_NclsTPC>90");
451 c->Draw("M>>hM(50,2,4)","cut","e");
453 c->SetMarkerColor(kBlue);
454 c->SetLineColor(kBlue);
455 // c->SetAlias("cut","PairType==1&&nCls&&pidParam&&LegEta&&Rap")
456 c->SetAlias("cut","PairType==1&&QA&&pidSig&&LegNclPID")
457 // c->SetAlias("nCls","Leg1_NclsTPC>140&&Leg2_NclsTPC>140");
458 c->Draw("M>>hM2(50,2,4)","cut","esame");
461 c->SetMarkerColor(kGreen);
462 c->SetLineColor(kGreen);
463 c->SetAlias("cut","PairType==1&&QA&&pidSig")
464 c->SetAlias("LegNcl","Leg1_NclsTPC>120&&Leg2_NclsTPC>120");
465 // c->SetAlias("nCls","Leg1_NclsTPC>150&&Leg2_NclsTPC>150");
466 c->Draw("M>>hM3(50,2,4)","cut","esame");
471 c->SetAlias("LegNclDiffIter1","(Leg1_NclsTPC-Leg1_NclsTPCiter1)>-1&&(Leg2_NclsTPC-Leg2_NclsTPCiter1)>-1")
472 c->SetAlias("cut","PairType==1&&LegNclDiffIter1")
474 AliDielectronHistos h("h","h");
475 h.AddClass("TPCsignal");
476 h.UserHistogram("TPCsignal","sigTPC","TPC signal;P [GeV];TPC signal [arb. Units]",400,.3,40,400,0.,200.,0,0,kTRUE,kFALSE)
477 h.GetHistogram("TPCsignal","sigTPC")->SetDirectory(gDirectory)
479 h.UserHistogram("TPCsignal","nSigE","TPC n #sigma Electrons;P [GeV];TPC n #sigma Electrons",200,.3,40.,100,-10.,10.,0,0,kTRUE,kFALSE)
480 h.GetHistogram("TPCsignal","nSigE")->SetDirectory(gDirectory)
481 h.UserHistogram("TPCsignal","nSigMu","TPC n #sigma Muons;P [GeV];TPC n #sigma Muons",400,.3,40.,500,-10.,10.,0,0,kTRUE,kFALSE)
482 h.GetHistogram("TPCsignal","nSigMu")->SetDirectory(gDirectory)
483 h.UserHistogram("TPCsignal","nSigPi","TPC n #sigma Pions;P [GeV];TPC n #sigma Pions",400,.3,40.,500,-10.,10.,0,0,kTRUE,kFALSE)
484 h.GetHistogram("TPCsignal","nSigPi")->SetDirectory(gDirectory)
485 h.UserHistogram("TPCsignal","nSigK","TPC n #sigma Kaons;P [GeV];TPC n #sigma Kaons",400,.3,40,500,-10,10,0,0,kTRUE,kFALSE)
486 h.GetHistogram("TPCsignal","nSigK")->SetDirectory(gDirectory)
487 h.UserHistogram("TPCsignal","nSigP","TPC n #sigma Protons;P [GeV];TPC n #sigma Protons",400,.3,40.,500,-10,10.,0,0,kTRUE,kFALSE)
488 h.GetHistogram("TPCsignal","nSigP")->SetDirectory(gDirectory)
490 h.UserHistogram("TPCsignal","nSigDiffP","ncls-nclsXX;P [GeV];ncls-nclsXX",400,.3,40.,200,-40,160.,0,0,kTRUE,kFALSE)
491 h.GetHistogram("TPCsignal","nSigDiffP")->SetDirectory(gDirectory)
494 c->Draw("Leg1_TPC_signal:Leg1_P_InnerParam>>sigTPC","cut","colz")
495 c->Draw("Leg2_TPC_signal:Leg2_P_InnerParam>>+sigTPC","cut","colz")
498 c->Draw("Leg1_TPC_nSigma_Electrons:Leg1_P_InnerParam>>nSigE","cut","colz")
499 c->Draw("Leg2_TPC_nSigma_Electrons:Leg2_P_InnerParam>>+nSigE","cut","colz")
501 c->Draw("Leg1_TPC_nSigma_Muos:Leg1_P_InnerParam>>nSigMu","cut","goff")
502 c->Draw("Leg2_TPC_nSigma_Muons:Leg1_P_InnerParam>>nSigMu","cut","goff")
503 c->Draw("Leg2_TPC_nSigma_Muons:Leg2_P_InnerParam>>+nSigMu","cut","colz")
505 c->Draw("Leg1_TPC_nSigma_Pions:Leg1_P_InnerParam>>nSigPi","cut","goff")
506 c->Draw("Leg2_TPC_nSigma_Pions:Leg2_P_InnerParam>>+nSigPi","cut","colz")
508 c->Draw("Leg1_TPC_nSigma_Kaons:Leg1_P_InnerParam>>nSigK","cut","goff")
509 c->Draw("Leg2_TPC_nSigma_Kaons:Leg2_P_InnerParam>>+nSigK","cut","colz")
511 c->Draw("Leg1_TPC_nSigma_Protons+.5:Leg1_P_InnerParam>>nSigP","cut","goff")
512 c->Draw("Leg2_TPC_nSigma_Protons+.5:Leg2_P_InnerParam>>+nSigP","cut","colz")
516 c->Draw("Leg1_TOF_nSigma_Electrons:Leg1_P_InnerParam>>nSigE","cut","colz")
517 c->Draw("Leg2_TOF_nSigma_Electrons:Leg2_P_InnerParam>>+nSigE","cut","colz")
520 c->Draw("Leg1_TOF_nSigma_Protons:Leg1_P_InnerParam>>nSigP","cut","goff")
521 c->Draw("Leg2_TOF_nSigma_Protons:Leg2_P_InnerParam>>+nSigP","cut","colz")
524 Pair->SetScanField(0)
525 Pair->Scan("EventInFile:File.GetString()","","colsize=1 col=3.d:100.s")
528 AliDielectronSignalFunc sig;
533 c->Draw("Leg1_NclsTPC>>hNcls(160,-0.5,159.5)","cut","goff");
534 c->Draw("Leg2_NclsTPC>>+hNcls","cut","goff");
535 hNclsPID->SetLineColor(kBlack)
537 c->Draw("Leg1_TPCsignalN>>hNclsPID(160,-0.5,159.5)","cut","goff");
538 c->Draw("Leg2_TPCsignalN>>+hNclsPID","cut","goff");
539 hNclsPID->SetLineColor(kBlue)
541 c->Draw("Leg1_NclsTPCiter1>>hNclsIter1(160,-0.5,159.5)","cut","goff");
542 c->Draw("Leg2_NclsTPCiter1>>+hNclsIter1","cut","goff");
543 hNclsIter1->SetLineColor(kGreen)
546 hNclsPID->Draw("same");
547 hNclsIter1->Draw("same");
552 c->Draw("Leg1_TPCsignalN:Leg1_NclsTPC>>hNclsPIDNcls(160,-0.5,159.5,160,-0.5,159.5)","cut","colz");
553 c->Draw("Leg2_TPCsignalN:Leg2_NclsTPC>>+hNclsPIDNcls","cut","colz");
556 c->Draw("Leg1_NclsTPC-Leg1_TPCsignalN:Leg1_P_InnerParam>>nSigDiffP","cut","colz");
557 c->Draw("Leg2_NclsTPC-Leg2_TPCsignalN:Leg1_P_InnerParam>>+nSigDiffP","cut","colz");
560 c->Draw("Leg1_NclsTPC-Leg1_NclsTPCiter1:Leg1_P_InnerParam>>nSigDiffP","cut","colz");
561 c->Draw("Leg2_NclsTPC-Leg2_NclsTPCiter1:Leg1_P_InnerParam>>+nSigDiffP","cut","colz");
571 c->Draw("Leg1_NclsTPC-Leg1_TPCsignalN:Leg1_TPC_nSigma_Electrons:Leg1_P_InnerParam>>hXX(100,0,10,20,-4,4)","cut","profcolz")
572 c->Draw("Leg1_NclsTPC-Leg1_TPCsignalN:Leg2_TPC_nSigma_Electrons:Leg2_P_InnerParam>>+hXX","cut","profcolz")
583 c->SetAlias("cutE","abs(Leg1_TPC_nSigma_Electrons)<3&&abs(Leg2_TPC_nSigma_Electrons)<3");
584 c->SetAlias("cutPi","abs(Leg1_TPC_nSigma_Pions)>5&&abs(Leg2_TPC_nSigma_Pions)>5");
585 // c->SetAlias("cutPi","Leg1_TPC_signal>65&&Leg2_TPC_signal>65");
586 c->SetAlias("cutP","(Leg1_TPC_nSigma_Protons)>5&&(Leg2_TPC_nSigma_Protons)>5");
587 c->SetAlias("cutK","abs(Leg1_TPC_nSigma_Kaons)>2&&abs(Leg2_TPC_nSigma_Kaons)>2");
588 c->SetAlias("pid","cutE&&cutPi&&cutP");
589 c->SetAlias("etaLeg","abs(Leg1_Eta)<.9&&abs(Leg2_Eta)<.9");
590 c->SetAlias("rap","abs(Y)<.9");
593 c->SetAlias("cutAdd","PairType==1&&abs(Leg1_ImpactParXY)<.02&&abs(Leg2_ImpactParXY)<.02&&Leg2_Pt>1.")
594 c->Draw("M>>hM(50,2,4)","PairType==1&&pid","e");
599 sig.Draw("samestat");
605 c->SetAlias("cutAdd","PairType==1&&abs(Leg1_ImpactParXY)<.03&&abs(Leg2_ImpactParXY)<.03&&Leg2_Pt>1")
606 c->SetAlias("cut","cutAdd&&pid&etaLeg&&rap")
608 c->Draw("M>>hM(50,1.99,3.99)","cut","e");
610 sig.Draw("samestat");
613 c->SetAlias("cut","PairType==1")
614 c->Draw("Leg1_TPC_signal:Leg1_P_InnerParam>>sigTPC","cut","colz")
615 c->Draw("Leg2_TPC_signal:Leg2_P_InnerParam>>+sigTPC","cut","colz")
618 c->SetAlias("cutE","abs(Leg1_TPC_nSigma_Electrons)<2&&abs(Leg2_TPC_nSigma_Electrons)<2");
619 c->SetAlias("cutPi","abs(Leg1_TPC_nSigma_Pions)>2&&abs(Leg2_TPC_nSigma_Pions)>2");
620 c->SetAlias("cutPi2","Leg1_TPC_signal>65&&Leg2_TPC_signal>65");
621 c->SetAlias("cutPi3","abs(Leg1_TPC_nSigma_Pions)>2.5&&abs(Leg2_TPC_nSigma_Pions)>2.5");
622 c->SetAlias("cutP","abs(Leg1_TPC_nSigma_Protons)>2&&abs(Leg2_TPC_nSigma_Protons)>2");
623 c->SetAlias("cutP2","abs(Leg1_TPC_nSigma_Protons)>1.5&&abs(Leg2_TPC_nSigma_Protons)>1.5");
624 c->SetAlias("cutK","abs(Leg1_TPC_nSigma_Kaons)>2&&abs(Leg2_TPC_nSigma_Kaons)>2");
625 c->SetAlias("cutdXY","abs(Leg1_ImpactParXY)<.02&&abs(Leg2_ImpactParXY)<.02")
626 c->SetAlias("cutPt","Leg2_Pt>1.")
627 c->SetAlias("cutPt2","Leg2_Pt>1.2")
630 c->SetMarkerSize(0.7);
632 c->SetMarkerStyle(20);
633 c->SetLineColor(kRed);
634 c->SetMarkerColor(kRed);
635 c->Draw("M>>hM(100,2,4)","cutPi&&cutPt","e");
637 c->SetMarkerStyle(21);
638 c->SetLineColor(kRed-1);
639 c->SetMarkerColor(kRed-2);
640 c->Draw("M>>hM2(100,2,4)","cutPi2&&cutPt","esame");
642 c->SetMarkerStyle(22);
643 c->SetLineColor(kRed-2);
644 c->SetMarkerColor(kRed-2);
645 c->Draw("M>>hM3(100,2,4)","cutPi3&&cutPt","esame");
648 c->SetMarkerStyle(20);
649 c->SetLineColor(kBlue);
650 c->SetMarkerColor(kBlue);
651 c->Draw("M>>hM4(100,2,4)","cutPi&&cutPt&&cutP","esame");
653 c->SetMarkerStyle(21);
654 c->SetLineColor(kBlue-1);
655 c->SetMarkerColor(kBlue-1);
656 c->Draw("M>>hM5(100,2,4)","cutPi2&&cutPt&&cutP","esame");
658 c->SetMarkerStyle(22);
659 c->SetLineColor(kBlue-2);
660 c->SetMarkerColor(kBlue-2);
661 c->Draw("M>>hM6(100,2,4)","cutPi3&&cutPt&&cutP","esame");
665 c->SetMarkerStyle(20);
666 c->SetLineColor(kGreen);
667 c->SetMarkerColor(kGreen);
668 c->Draw("M>>hM7(100,2,4)","cutPi&&cutPt&&cutP2","esame");
670 c->SetMarkerStyle(21);
671 c->SetLineColor(kGreen-1);
672 c->SetMarkerColor(kGreen-1);
673 c->Draw("M>>hM8(100,2,4)","cutPi2&&cutPt&&cutP2","esame");
675 c->SetMarkerStyle(22);
676 c->SetLineColor(kGreen-2);
677 c->SetMarkerColor(kGreen-2);
678 c->Draw("M>>hM9(100,2,4)","cutPi3&&cutPt&&cutP2","esame");
683 c->SetMarkerStyle(20);
684 c->SetLineColor(kMagentha);
685 c->SetMarkerColor(kMagentha);
686 c->Draw("M>>hM7(100,2,4)","cutPi&&cutPt&&cutP2","esame");
688 c->SetMarkerStyle(21);
689 c->SetLineColor(kMagentha-1);
690 c->SetMarkerColor(kMagentha-1);
691 c->Draw("M>>hM8(100,2,4)","cutPi2&&cutPt&&cutP2","esame");
693 c->SetMarkerStyle(22);
694 c->SetLineColor(kMagentha-2);
695 c->SetMarkerColor(kMagentha-2);
696 c->Draw("M>>hM9(100,2,4)","cutPi3&&cutPt&&cutP2","esame");
699 c->SetLineColor(kBlack);
700 c->SetMarkerColor(kBlue);
701 c->Draw("M>>hM4(100,2,4)","cutE&&cutPi&&cutK&&cutP&&cutdXY&&cutPt","esame");
707 c->SetAlias("cutE","Leg1_TPC_nSigma_Electrons>-1.5&&Leg2_TPC_nSigma_Electrons>-1.5");
708 // c->SetAlias("cutE","Leg1_TPC_signal>60&&Leg2_TPC_signal>60");
709 c->SetAlias("cutP","abs(Leg1_TPC_nSigma_Protons)>3&&abs(Leg2_TPC_nSigma_Protons)>3")
711 c->SetAlias("cutAdd","PairType==1&&abs(Leg1_ImpactParXY)<.03&&abs(Leg2_ImpactParXY)<.03&&Leg2_Pt>0")
712 c->SetAlias("cut","Leg2_Pt>1&&cutE&&cutP")
714 c->Draw("M>>hM(50,1.99,3.99)","cut","e");
716 c->SetAlias("cutAdd","PairType==1&&abs(Leg1_ImpactParXY)<.03&&abs(Leg2_ImpactParXY)<.03&&Leg2_Pt>.8")
717 c->Draw("M>>hM2(50,1.99,3.99)","cut","esame");
719 c->SetAlias("cutAdd","PairType==1&&abs(Leg1_ImpactParXY)<.03&&abs(Leg2_ImpactParXY)<.03&&Leg2_Pt>1")
720 c->Draw("M>>hM3(50,1.99,3.99)","cut","esame");
722 c->SetAlias("cutAdd","PairType==1&&abs(Leg1_ImpactParXY)<.03&&abs(Leg2_ImpactParXY)<.03&&Leg2_Pt>1.2")
723 c->Draw("M>>hM4(50,1.99,3.99)","cut","esame");
726 c->Draw("Leg1_TPC_signal:Leg1_P_InnerParam>>sigTPC","cut","goff")
727 c->Draw("Leg2_TPC_signal:Leg2_P_InnerParam>>+sigTPC","cut","goff")
728 c1->Modified();c1->Update()
731 c->SetAlias("cutE","Leg1_TPC_nSigma_Electrons>-1.5&&Leg2_TPC_nSigma_Electrons>-1.5");
732 c->SetAlias("cut","Leg2_P_InnerParam>1.5&&cutE")
733 c->SetMarkerStyle(21);
734 c->Draw("M>>hM(50,1.99,3.99)","cut","e");
736 c->SetAlias("cutE","Leg1_TPC_nSigma_Electrons>-1.5+.8&&Leg2_TPC_nSigma_Electrons>-1.5+.8");
737 c->SetAlias("cut","Leg2_P_InnerParam>1.5&&cutE")
738 c->SetMarkerStyle(20);
739 c->SetMarkerColor(kRed);
740 c->Draw("M>>hM3(50,1.99,3.99)","cut","esame");
745 //=================== pass 1 pass 2 comparison =======================
748 //------ cuts ------------
751 c->SetAlias("cutPro","Leg1_TPC_signal>50*exp(-.5*(Leg1_P_InnerParam-2))&&Leg2_TPC_signal>50*exp(-.5*(Leg2_P_InnerParam-2))&&Leg1_TPC_signal<85&&Leg2_TPC_signal<85");
753 c->SetAlias("Pcut","Leg1_P_InnerParam>1.3&&Leg2_P_InnerParam>1.3")
754 c->SetAlias("Ptcut","Leg1_Pt>1&&Leg2_Pt>1")
755 c->SetAlias("TOFcut","abs(Leg1_TOF_nSigma_Electrons)<3&&abs(Leg2_TOF_nSigma_Electrons)<3");
756 c->SetAlias("TOFcut2","(Leg1_P_InnerParam<1.3&&abs(Leg1_TOF_nSigma_Electrons)<3||Leg1_P_InnerParam>=1.3)&&(Leg2_P_InnerParam<1.3&&abs(Leg2_TOF_nSigma_Electrons)<3||Leg2_P_InnerParam>=1.3)");
757 c->SetAlias("TPCcut","(Leg1_TPC_signal>70*(1-exp(-1*(Leg1_P_InnerParam+2))))&&(Leg2_TPC_signal>70*(1-exp(-1*(Leg2_P_InnerParam+2))))")
758 c->SetAlias("NClcut","Leg1_NclsTPC>90&&Leg2_NclsTPC>90");
760 c->SetAlias("eleParam","Leg1_TPC_nSigma_Electrons<5&&Leg2_TPC_nSigma_Electrons<5&&Leg1_TPC_nSigma_Electrons>-2.65*exp(-0.6757*Leg1_P_InnerParam)&&Leg2_TPC_nSigma_Electrons>-2.65*exp(-0.6757*Leg2_P_InnerParam)")
761 c->SetAlias("cut","PairType==1&&eleParam&&Run<127719")
762 c->SetAlias("cut","1==1")
763 c->SetAlias("cut","NClcut")
765 c->SetAlias("cut","TOFcut&&TPCcut&&NClcut")
767 c->SetAlias("cut","TOFcut2&&TPCcut&&NClcut")
769 c->SetAlias("cut","cutPro&&TPCcut&&NClcut")
771 c->SetAlias("cut","Pcut&&TPCcut&&NClcut")
773 c->SetAlias("cut","Ptcut&&TPCcut&&NClcut")
777 //------------ plots --------------
780 c->SetAlias("cut","1==1")
783 c->SetLineColor(kBlack);
784 c->Draw("M>>hM(50,1.99,3.99)","cut","e");
785 hM->SetTitle(";Inv. Mass [GeV]; Pair (e+e-) / 40GeV")
788 c1->Print("pics/M_noCut.png");
793 c->SetAlias("cut","1==1")
796 gStyle->SetOptStat(0);
797 c->Draw("Leg1_NclsTPC:Leg1_TPC_signal:Leg1_P_InnerParam>>test(1000,0,40,200,60,100)","cut","profcolz")
798 c->Draw("Leg2_NclsTPC:Leg2_TPC_signal:Leg2_P_InnerParam>>+test","cut","profcolz")
800 test->SetTitle("mean TPC number of clusters;P_{TPC} [GeV]; TPC signal [arb. units]")
803 c1->Print("pics/TPCnCl_P.png");
808 c->SetAlias("cut","NClcut")
809 c->Draw("Leg1_NclsTPC:Leg1_TPC_signal:Leg1_P_InnerParam>>test(1000,0,40,200,60,100)","cut","profcolz")
810 c->Draw("Leg2_NclsTPC:Leg2_TPC_signal:Leg2_P_InnerParam>>+test","cut","profcolz")
812 test->SetTitle("mean TPC number of clusters;P_{TPC} [GeV]; TPC signal [arb. units]")
815 c1->Print("pics/TPCnCl_P_cutNcl.png");
819 //tpc signal + signal cut
824 h.GetHistogram("TPCsignal","sigTPC")->GetYaxis()->SetRangeUser(60,100);
825 c->SetAlias("cut","NClcut")
826 c->Draw("Leg1_TPC_signal:Leg1_P_InnerParam>>sigTPC","cut","colz")
827 c->Draw("Leg2_TPC_signal:Leg2_P_InnerParam>>+sigTPC","cut","colz")
828 TF1 f("f1","[0]*(1-exp(-[1]*(x-[2])))",0.3,40);
829 f.SetParameters(70,1,-2);
833 c1->Print("pics/TPCsignal_P_cutNcl.png");
840 c->SetAlias("cut","1==1")
841 c->SetLineColor(kBlack);
842 c->SetMarkerColor(kBlack);
843 c->Draw("M>>hM(50,1.99,3.99)","cut","e");
844 hM->SetTitle(";Inv. Mass [GeV]; Pair (e+e-) / 40GeV")
846 c->SetAlias("cut","NClcut")
847 c->SetLineColor(kRed);
848 c->SetMarkerColor(kRed);
849 c->Draw("M>>hM2(50,1.99,3.99)","cut","esame");
852 c1->Print("pics/M_nClCut.png");
856 //tpc signal: ncl + tpc cut
861 c->SetAlias("cut","TPCcut&&NClcut")
862 c->Draw("Leg1_TPC_signal:Leg1_P_InnerParam>>sigTPC","cut","colz")
863 c->Draw("Leg2_TPC_signal:Leg2_P_InnerParam>>+sigTPC","cut","colz")
866 c1->Print("pics/TPCsignal_P_cutNcl_tpc.png");
873 c->SetAlias("cut","1==1")
874 c->SetLineColor(kBlack);
875 c->SetMarkerColor(kBlack);
876 c->Draw("M>>hM(50,1.99,3.99)","cut","e");
877 hM->SetTitle(";Inv. Mass [GeV]; Pair (e+e-) / 40GeV")
879 c->SetAlias("cut","NClcut")
880 c->SetLineColor(kRed);
881 c->SetMarkerColor(kRed);
882 c->Draw("M>>hM2(50,1.99,3.99)","cut","esame");
883 c->SetAlias("cut","TPCcut&&NClcut")
884 c->SetLineColor(kGreen);
885 c->SetMarkerColor(kGreen);
886 c->Draw("M>>hM3(50,1.99,3.99)","cut","esame");
889 c1->Print("pics/M_nClCut_tpc.png");
893 //TPC signal: ncl + tpc + tof cut
898 c->SetAlias("cut","TOFcut2&&TPCcut&&NClcut")
899 c->Draw("Leg1_TPC_signal:Leg1_P_InnerParam>>sigTPC","cut","colz")
900 c->Draw("Leg2_TPC_signal:Leg2_P_InnerParam>>+sigTPC","cut","colz")
903 c1->Print("pics/TPCsignal_P_cutNcl_tpc.png");
910 c->SetAlias("cut","1==1")
911 c->SetAlias("cut","TPCcut&&NClcut")
912 c->SetLineColor(kGreen);
913 c->Draw("M>>hM(50,1.99,3.99)","cut","e");
914 hM->SetTitle(";Inv. Mass [GeV]; Pair (e+e-) / 40GeV")
916 c->SetAlias("cut","TOFcut2&&TPCcut&&NClcut")
917 c->SetLineColor(kBlue);
918 c->Draw("M>>hM2(50,1.99,3.99)","cut","esame");
921 c1->Print("pics/M_nClCut_tpc.png");
924 //Inv Mass: different cuts
930 c->SetAlias("cut","Ptcut&&TPCcut&&NClcut")
931 c->SetLineColor(kMagenta);
932 c->SetMarkerColor(kMagenta);
933 c->SetMarkerStyle(22);
934 c->Draw("M>>hM(50,1.99,3.99)","cut","e");
936 c->SetAlias("cut","Pcut&&TPCcut&&NClcut")
937 c->SetLineColor(kCyan+1);
938 c->SetMarkerColor(kCyan+1);
939 c->SetMarkerStyle(21);
940 c->Draw("M>>hM2(50,1.99,3.99)","cut","esame");
942 c->SetAlias("cut","TOFcut2&&TPCcut&&NClcut")
943 c->SetMarkerStyle(20);
944 c->SetLineColor(kBlue);
945 c->SetMarkerColor(kBlue);
946 c->Draw("M>>hM3(50,1.99,3.99)","cut","esame");
950 c1->Print("pics/M_nClCut_tpc_tof.png");
955 c->SetAlias("NClcut","Leg1_NclsTPC>90&&Leg2_NclsTPC>90");
956 c->SetAlias("Ptcut","Leg1_Pt>1&&Leg2_Pt>1")
957 c->SetAlias("PairT","PairType==1");
958 c->SetAlias("cut","NClcut&&Ptcut&&PairT")
961 c->SetAlias("cutE","abs(Leg1_TPC_nSigma_Electrons)<3&&abs(Leg2_TPC_nSigma_Electrons)<3");
962 c->SetAlias("cutPi","abs(Leg1_TPC_nSigma_Pions)>3&&abs(Leg2_TPC_nSigma_Pions)>3");
963 c->SetAlias("cutP","(Leg1_TPC_nSigma_Protons)>3&&(Leg2_TPC_nSigma_Protons)>3");
964 c->SetAlias("cut","NClcut&&Ptcut&&PairT&&cutE&&cutPi&&cutP")
966 c->SetAlias("eta","abs(Leg1_Eta)<0.88&&abs(Leg1_Eta)<0.88");
967 c->SetAlias("rap","abs(Y)<0.88");
968 c->SetAlias("spdFirst","(Leg1_ITS_clusterMap&1)==1&&(Leg2_ITS_clusterMap&1)==1")
969 c->SetAlias("cut","NClcut&&Ptcut&&PairT&&cutE&&cutPi&&cutP&&eta&&rap&&spdFirst")
970 c->Draw("M>>hM(50,1.98,1.98+50*.04)","cut","e")
973 c->SetAlias("cutProL1",Form("Leg1_TPC_nSigma_Electrons>(%f*%f+(AliExternalTrackParam::BetheBlochAleph(Leg1_P_InnerParam/%f,0.0283086/0.97,2.63394e+01,5.04114e-11,2.12543e+00,4.88663e+00)-AliExternalTrackParam::BetheBlochAleph(Leg1_P_InnerParam/%f,0.0283086/0.97,2.63394e+01,5.04114e-11,2.12543e+00,4.88663e+00)))/%f",nSigma,resolution, AliPID::ParticleMass(AliPID::kProton), AliPID::ParticleMass(AliPID::kElectron), resolution))
975 c->SetAlias("cutProL2",Form("Leg2_TPC_nSigma_Electrons>(%f*%f+(AliExternalTrackParam::BetheBlochAleph(Leg2_P_InnerParam/%f,0.0283086/0.97,2.63394e+01,5.04114e-11,2.12543e+00,4.88663e+00)-AliExternalTrackParam::BetheBlochAleph(Leg2_P_InnerParam/%f,0.0283086/0.97,2.63394e+01,5.04114e-11,2.12543e+00,4.88663e+00)))/%f",nSigma,resolution, AliPID::ParticleMass(AliPID::kProton), AliPID::ParticleMass(AliPID::kElectron), resolution))
977 c->SetAlias("cutPioL1",Form("Leg1_TPC_nSigma_Electrons>(%f*%f+(AliExternalTrackParam::BetheBlochAleph(Leg1_P_InnerParam/%f,0.0283086/0.97,2.63394e+01,5.04114e-11,2.12543e+00,4.88663e+00)-AliExternalTrackParam::BetheBlochAleph(Leg1_P_InnerParam/%f,0.0283086/0.97,2.63394e+01,5.04114e-11,2.12543e+00,4.88663e+00)))/%f",nSigma,resolution, AliPID::ParticleMass(AliPID::kPion), AliPID::ParticleMass(AliPID::kElectron), resolution))
979 c->SetAlias("cutPioL2",Form("Leg2_TPC_nSigma_Electrons>(%f*%f+(AliExternalTrackParam::BetheBlochAleph(Leg2_P_InnerParam/%f,0.0283086/0.97,2.63394e+01,5.04114e-11,2.12543e+00,4.88663e+00)-AliExternalTrackParam::BetheBlochAleph(Leg2_P_InnerParam/%f,0.0283086/0.97,2.63394e+01,5.04114e-11,2.12543e+00,4.88663e+00)))/%f",nSigma,resolution, AliPID::ParticleMass(AliPID::kPion), AliPID::ParticleMass(AliPID::kElectron), resolution))
981 c->SetAlias("cutPro","cutProL1&&cutProL2");
982 c->SetAlias("cutPio","cutPioL1&&cutPioL2");