]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWG3/dielectron/macros/MakeDataReport.C
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Add AliEMCALRecoUtils pointer to remove bad or exotic clusters for cluster histograms...
[u/mrichter/AliRoot.git] / PWG3 / dielectron / macros / MakeDataReport.C
CommitLineData
8df8e382 1void SetupStyle();
2
3void MakeDataReport(const char* outputFile="JpsiDataReport.pdf",
4 const char* histos="jpsi_HistosSE.root",
5 const char* cf="jpsi_CF.root")
6{
7 //
8 // Make a pdf file with the efficiency report
9 //
10
11 SetupStyle();
12
13 AliDielectronCFdraw d(cf);
14 d.SetRangeUser("PairType",1,1);
15 d.SetRangeUser("Y",-.89,.9,"0");
16
17
18 TFile f("jpsi_HistosSE.root");
19
20 AliDielectronHistos h;
21 TIter nextHists((TList*)f.Get("Dielectron_Histos"));
22
23 TPaveText pt(.02,.6,.98,.8);
24 TText *t1=pt.AddText("");
25 TText *t2=pt.AddText("");
26
27 TCanvas *c1=new TCanvas;
28
29 TPDF p(outputFile);
30
31 //
32 // Invariant mass plots
33 //
34
35
36 //
37 // Make QA info
38 //
39
40 t1->SetTitle("QA summary plots for");
41 THashList *list=0x0;
42 while ( (list=(THashList*)nextHists()) ){
43 h.SetHistogramList(*list);
44 t2->SetTitle(list->GetName());
45 pt.Draw();
46 c1->Update();
47 h.Draw();
48 c1->Clear();
49 }
50 p.Close();
51 delete c1;
52}
53
54void SetupStyle()
55{
56 const Int_t NCont=255;
57
58 TStyle *st = new TStyle("mystyle","mystyle");
59 gROOT->GetStyle("Plain")->Copy((*st));
60 st->SetTitleX(0.1);
61 st->SetTitleW(0.8);
62 st->SetTitleH(0.08);
63 st->SetStatX(.9);
64 st->SetStatY(.9);
65 st->SetNumberContours(NCont);
66 st->SetPalette(1,0);
67 st->SetOptStat("erm");
68 st->SetOptFit(0);
69 st->SetGridColor(kGray+1);
70 st->SetPadGridX(kTRUE);
71 st->SetPadGridY(kTRUE);
72 st->SetPadTickX(kTRUE);
73 st->SetPadTickY(kTRUE);
74 st->cd();
75
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 };
81
82 TColor::CreateGradientColorTable(NRGBs, stops, red, green, blue, NCont);
83
84}
85
86void DrawUnbinned(){
87 TFile f("jpsi_debug.root");
88// if (!f.IsOpen()) return;
89
90 TTree *t=(TTree*)f.Get("Pair");
91// if (!t) return;
92
93 TCanvas c1;
94 gPad->SetLogy();
95 gStyle->SetOptStat(0);
96
97 TLegend *leg=new TLegend(0.59,.79,.99,.99);
98 TLine l;
99
100 l.SetLineColor(kGreen-5);
101 l.SetLineWidth(2);
102 l.SetLineStyle(2);
103 leg->SetFillColor(10);
104
105 leg->Clear();
106
107
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");
114
115 l.DrawLine(3.097,1,3.097,1e4);
116
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");
121
122 TCut d1_1="abs(Leg1_TPC_nSigma_Electrons)<1";
123 TCut d2_1="abs(Leg2_TPC_nSigma_Electrons)<1";
124 TCut d_1=d1_1+d2_1;
125
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");
130
131
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");
136
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");
141
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");
146
147 leg->Draw();
148}
149
150/*
151 Double_t alephParameters[5];
152 // simulation
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;
158 Color_t color=kRed;
159
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);
165 //
166 foProton->SetParameters(alephParameters);
167 foPion->SetParameters(alephParameters);
168 foElec->SetParameters(alephParameters);
169 foKaon->SetParameters(alephParameters);
170 foMuon->SetParameters(alephParameters);
171 //
172 foProton->SetLineColor(color);
173 foPion->SetLineColor(color);
174 foElec->SetLineColor(color);
175 foKaon->SetLineColor(color);
176 foMuon->SetLineColor(color);
177 //
178 Int_t lineWidth=1;
179 foProton->SetLineWidth(lineWidth);
180 foPion->SetLineWidth(lineWidth);
181 foElec->SetLineWidth(lineWidth);
182 foKaon->SetLineWidth(lineWidth);
183 foMuon->SetLineWidth(lineWidth);
184
185 //
186 foProton->SetNpx(200);
187 foPion->SetNpx(200);
188 foElec->SetNpx(200);
189 foKaon->SetNpx(200);
190 foMuon->SetNpx(200);
191 //
192 foProton->Draw("same");
193 foPion->Draw("same");
194 foElec->Draw("same");
195 foKaon->Draw("same");
196 foMuon->Draw("same");
197
198
199
200
201
202 // data
554e40f8 203 Double_t res=5.2e-2;
8df8e382 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;
209 Color_t color=kRed;
210
554e40f8 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;
8df8e382 217
218
219 TF1 *foDataProton = new TF1("foDataProton", "50*AliExternalTrackParam::BetheBlochAleph(x/0.93827,[0],[1],[2],[3],[4])",0.05,20);
554e40f8 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);
8df8e382 221 TF1 *foDataProtonM = new TF1("foDataProtonM", Form("50*AliExternalTrackParam::BetheBlochAleph(x/0.93827,[0],[1],[2],[3],[4])*(1-%f)",res),0.05,20);
222
223 TF1 *foDataPion = new TF1("foDataPion", "50*AliExternalTrackParam::BetheBlochAleph(x/0.13957,[0],[1],[2],[3],[4])",0.05,20);
554e40f8 224 TF1 *foDataPionP = new TF1("foDataPionP",Form( "50*AliExternalTrackParam::BetheBlochAleph(x/0.13957,[0],[1],[2],[3],[4])*(1+%f)",res),0.05,20);
8df8e382 225 TF1 *foDataPionM = new TF1("foDataPionM", Form("50*AliExternalTrackParam::BetheBlochAleph(x/0.93827,[0],[1],[2],[3],[4])*(1-%f)",res),0.05,20);
226
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);
230
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);
234
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);
238
239 //
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);
255 //
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);
271 //
272 Int_t lineWidth=1;
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);
288
289 //
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);
300
301 //
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);
311 //
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");
319
320
321
322
323
324{
325
326 Int_t baseColors[5]={kRed, kGreen+1, kAzure-4, kMagenta, kCyan+1};
327 Int_t sigmaColorOffset=1;
328
329Int_t baseColors[5]={kRed, kGreen+1, kAzure-4, kMagenta, kCyan+1};
330 Int_t sigmaColorOffset=0;
2a14a7b1 331Int_t baseColors[5]={kRed, kRed, kRed, kRed, kRed};
8df8e382 332
48609e3d 333 Double_t sigmas[5]={3,3,3,3,3};
8df8e382 334 Double_t masses[5];
335
336 for (Int_t i=0; i<AliPID::kSPECIES; ++i) masses[i]=AliPID::ParticleMass(i);
337
338 Double_t res=7.e-2;
339 Double_t alephParameters[5];
340
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;
ffbede40 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;
351
2a14a7b1 352 Double_t mip=50;
8df8e382 353
354 Color_t color=kRed;
355 Int_t lineWidth=2;
356
357TF1 *fBethe[5];
358TF1 *fBetheP[5];
359TF1 *fBetheM[5];
360
361for (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);
365
366 fBethe[i]->SetParameters(alephParameters);
367 fBetheP[i]->SetParameters(alephParameters);
368 fBetheM[i]->SetParameters(alephParameters);
369
370 fBethe[i]->SetLineColor(baseColors[i]);
371 fBetheP[i]->SetLineColor(baseColors[i]-sigmaColorOffset);
372 fBetheM[i]->SetLineColor(baseColors[i]-sigmaColorOffset);
373
374 fBethe[i]->SetLineWidth(lineWidth);
375 fBetheP[i]->SetLineWidth(lineWidth);
376 fBetheM[i]->SetLineWidth(lineWidth);
377
378 fBetheP[i]->SetLineStyle(2);
379 fBetheM[i]->SetLineStyle(2);
380
381 fBethe[i]->SetNpx(200);
382 fBetheP[i]->SetNpx(200);
383 fBetheM[i]->SetNpx(200);
384}
385
386for (Int_t i=0; i<5; ++i){
387 fBethe[i]->Draw("same");
48609e3d 388// fBetheP[i]->Draw("same");
389// fBetheM[i]->Draw("same");
8df8e382 390}
391}
392
393*/
394
8df8e382 395
554e40f8 396/*
8df8e382 397
554e40f8 398 Double_t resolution=0.052;
399 Double_t nSigma=3.;
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);
402fPPio->Draw("same");
8df8e382 403
554e40f8 404 TF1 f("fP","-8*exp(-0.6*x)",0,40);
405f.Draw("same")
8df8e382 406
8df8e382 407
554e40f8 408*/
8df8e382 409
410
554e40f8 411/*
fb7d2d99 412//
413// Default
414//
415c->SetAlias("cutE","abs(Leg1_TPC_nSigma_Electrons)<3&&abs(Leg2_TPC_nSigma_Electrons)<3");
ba15fdfb 416c->SetAlias("cutPi","abs(Leg1_TPC_nSigma_Pions)>3.5&&abs(Leg2_TPC_nSigma_Pions)>3.5");
417c->SetAlias("cutP","(Leg1_TPC_nSigma_Protons)>3.5&&(Leg2_TPC_nSigma_Protons)>3.5");
fb7d2d99 418c->SetAlias("pidSig","cutE&&cutPi&&cutP");
419
420c->SetAlias("LegEta","abs(Leg1_Eta)<0.9&&abs(Leg2_Eta)<0.9");
ba15fdfb 421c->SetAlias("LegNcl","Leg1_NclsTPC>70&&Leg2_NclsTPC>70");
fb7d2d99 422c->SetAlias("LegPt","Leg1_Pt>1&&Leg2_Pt>1");
423c->SetAlias("Rap","abs(Y)<0.9");
424c->SetAlias("QA","LegNcl&&LegEta&&Rap&&LegPt");
425
426c->SetAlias("spdFirst","(Leg1_ITS_clusterMap&1)==1 && (Leg2_ITS_clusterMap&1)==1");
427
428c->SetAlias("cut","PairType==1&&QA&&pidSig");
429
430c->Draw("M>>hM(125,0.,5.)","cut","e");
8df8e382 431
8df8e382 432
fb7d2d99 433
434
435//---------------------
436
554e40f8 437c->SetAlias("cut","PairType==1")
8df8e382 438
8df8e382 439
554e40f8 440c->SetAlias("nCls","Leg1_NclsTPC>90&&Leg2_NclsTPC>90");
8df8e382 441
fb7d2d99 442
443c->SetAlias("cutE","Leg1_TPC_nSigma_Electrons>-1 && Leg2_TPC_nSigma_Electrons>-1");
444c->SetAlias("LegPt","Leg1_Pt>1&&Leg2_Pt>1");
445
446
447
554e40f8 448//--------PID
45b2b1b8 449//-Param dEdx
450c->SetAlias("cutPipardEdx","Leg1_TPC_signal>75-20*exp(-.7*Leg1_P_InnerParam)&&Leg2_TPC_signal>75-20*exp(-.7*Leg2_P_InnerParam)")
451 c->SetAlias("TOFe1r","(((Leg1_TrackStatus&32768)==32768)&&abs(Leg1_TOF_nSigma_Electrons)<3)");
452 c->SetAlias("TOFe2r","(((Leg2_TrackStatus&32768)==32768)&&abs(Leg2_TOF_nSigma_Electrons)<3)");
453
454c->SetAlias("cutPspecial","(abs(Leg1_TPC_nSigma_Protons)>3||(abs(Leg1_TPC_nSigma_Protons)<=3&&TOFe1r))&&(abs(Leg2_TPC_nSigma_Protons)>3||(abs(Leg2_TPC_nSigma_Protons)<=3&&TOFe2r))")
8df8e382 455
554e40f8 456//-- nsigma
457c->SetAlias("cutE","abs(Leg1_TPC_nSigma_Electrons)<3&&abs(Leg2_TPC_nSigma_Electrons)<3");
ffbede40 458// c->SetAlias("cutE","Leg1_TPC_nSigma_Electrons>-1 && Leg2_TPC_nSigma_Electrons>-1");
554e40f8 459c->SetAlias("cutPi","abs(Leg1_TPC_nSigma_Pions)>3&&abs(Leg2_TPC_nSigma_Pions)>3");
fb7d2d99 460c->SetAlias("cutP","(Leg1_TPC_nSigma_Protons)>3.&&(Leg2_TPC_nSigma_Protons)>3.3");
554e40f8 461c->SetAlias("pidSig","cutE&&cutPi&&cutP");
554e40f8 462//-- Pi param
45b2b1b8 463// c->SetAlias("eleParam","Leg1_TPC_nSigma_Electrons<5&&Leg2_TPC_nSigma_Electrons<5&&Leg1_TPC_nSigma_Electrons>-2.65*exp(-0.9*Leg1_P_InnerParam)&&Leg2_TPC_nSigma_Electrons>-8*exp(-0.6*Leg2_P_InnerParam)");
464c->SetAlias("eleParam","Leg1_TPC_nSigma_Electrons<5&&Leg2_TPC_nSigma_Electrons<5&&Leg1_TPC_nSigma_Electrons>-3.7*exp(-0.9*Leg1_P_InnerParam)-0.1&&Leg2_TPC_nSigma_Electrons>-3.7*exp(-0.6*Leg2_P_InnerParam)-0.1");
554e40f8 465c->SetAlias("pidParam","eleParam&&cutP");
8df8e382 466
467
468
554e40f8 469c->SetAlias("LegEta","abs(Leg1_Eta)<0.9&&abs(Leg2_Eta<0.9)");
ffbede40 470c->SetAlias("LegNcl","Leg1_NclsTPC>90&&Leg2_NclsTPC>90");
45b2b1b8 471c->SetAlias("LegPt","Leg1_Pt>1&&Leg2_Pt>1");
554e40f8 472c->SetAlias("Rap","abs(Y)<0.9");
473c->SetAlias("QA","LegNcl&&LegEta&&Rap");
474c->SetAlias("spdFirst","(Leg1_ITS_clusterMap&1)==1 && (Leg2_ITS_clusterMap&1)==1");
ffbede40 475c->SetAlias("LegNclDiffIter1","abs(Leg1_NclsTPC-Leg1_NclsTPCiter1)<10&&abs(Leg2_NclsTPC-Leg2_NclsTPCiter1)<10")
476c->SetAlias("LegNclPID","(Leg1_NclsTPC-Leg1_TPCsignalN)<20&&(Leg2_NclsTPC-Leg2_TPCsignalN)<20")
8df8e382 477
fb7d2d99 478c->SetAlias("LegNcl","Leg1_NFclsTPCrobust&&Leg2_NFclsTPCrobust");
479c->SetAlias("cut","PairType==1&&QA&&pidSig&&LegPt");
480
45b2b1b8 481c->SetAlias("cut","PairType==1&&QA&&cutPipardEdx&&cutPspecial")
482c->SetAlias("cut","PairType==1&&QA&&pidSig&&LegPt")
8df8e382 483
fb7d2d99 484c->SetMarkerStyle(22);
554e40f8 485c->SetMarkerSize(.8);
fb7d2d99 486
487//-------- nsigma
488
554e40f8 489c->SetMarkerColor(kBlack);
8df8e382 490c->SetLineColor(kBlack);
8df8e382 491
fb7d2d99 492c->SetAlias("cutPi","abs(Leg1_TPC_nSigma_Pions)>3&&abs(Leg2_TPC_nSigma_Pions)>3");
493c->SetAlias("cutP","(Leg1_TPC_nSigma_Protons)>3&&(Leg2_TPC_nSigma_Protons)>3");
494
495c->Draw("M>>hM(301,-.01,6.01)","cut","e");
496c->GetHistogram()->GetXaxis()->SetRangeUser(2,4);
497
498//--------
8df8e382 499
554e40f8 500c->SetMarkerColor(kBlue);
8df8e382 501c->SetLineColor(kBlue);
fb7d2d99 502
1201a1a9 503c->SetAlias("cutPi","abs(Leg1_TPC_nSigma_Pions)>3&&abs(Leg2_TPC_nSigma_Pions)>3");
504c->SetAlias("cutP","(Leg1_TPC_nSigma_Protons)>3&&(Leg2_TPC_nSigma_Protons)>3");
8df8e382 505
fb7d2d99 506c->Draw("M>>hM2(301,-.01,6.01)","cut","esame");
507
508//--------
ffbede40 509
554e40f8 510c->SetMarkerColor(kGreen);
8df8e382 511c->SetLineColor(kGreen);
fb7d2d99 512
513c->SetAlias("cutPi","abs(Leg1_TPC_nSigma_Pions)>4&&abs(Leg2_TPC_nSigma_Pions)>4");
514c->SetAlias("cutP","(Leg1_TPC_nSigma_Protons)>3.5&&(Leg2_TPC_nSigma_Protons)>3.5");
515
516c->Draw("M>>hM3(301,-.01,6.01)","cut","esame");
517
518//--------
519
520c->SetMarkerColor(kMagenta);
521c->SetLineColor(kMagenta);
522
45b2b1b8 523c->SetAlias("cutPi","abs(Leg1_TPC_nSigma_Pions)>3.5&&abs(Leg2_TPC_nSigma_Pions)>3.5");
fb7d2d99 524c->SetAlias("cutP","(Leg1_TPC_nSigma_Protons)>4&&(Leg2_TPC_nSigma_Protons)>4");
525
526c->Draw("M>>hM4(301,-.01,6.01)","cut","esame");
527
528
529
530
531
532//
533//-------- rapidity
534//
535
536c->SetMarkerColor(kBlack);
537c->SetLineColor(kBlack);
538
539c->Draw("M>>hM(301,-.01,6.01)","cut&&Y<=0","e");
540
541//--------
542
543c->SetMarkerColor(kBlue);
544c->SetLineColor(kBlue);
545
546c->Draw("M>>hM2(301,-.01,6.01)","cut&&Y>0","esame");
547
548
549
550
551
552c->SetAlias("cutE","Leg1_TPC_nSigma_Electrons>-1 && Leg2_TPC_nSigma_Electrons>-1");
553c->SetAlias("LegPt","Leg1_Pt>1.2&&Leg2_Pt>1.2");
554
555c->SetAlias("cut","PairType==1&&QA&&pidSig");
556
557
558//-------- binning
559
560c->SetMarkerColor(kBlack);
561c->SetLineColor(kBlack);
562
1201a1a9 563c->SetAlias("cutPi","abs(Leg1_TPC_nSigma_Pions)>3&&abs(Leg2_TPC_nSigma_Pions)>3");
564c->SetAlias("cutP","(Leg1_TPC_nSigma_Protons)>3&&(Leg2_TPC_nSigma_Protons)>3");
fb7d2d99 565
1201a1a9 566c->Draw("M>>hM(601,-.015,6.005)","cut","e");
fb7d2d99 567
568TGraphErrors gr0;
569TH1 *h=c->GetHistogram();
570
571gr0->SetLineColor(h->GetLineColor());
572gr0->SetMarkerColor(h->GetMarkerColor());
573for (Int_t i=0;i<h->GetNbinsX();++i){
574 gr0.SetPoint(i,h->GetXaxis()->GetBinCenter(i+1),h->GetBinContent(i+1));
575// gr0.SetPointError(i,h->GetXaxis()->GetBinWidth(i+1)/2,h->GetBinError(i+1));
576 gr0.SetPointError(i,0,h->GetBinError(i+1));
577}
578
579//--------
580
581c->SetMarkerColor(kBlue);
582c->SetLineColor(kBlue);
583
1201a1a9 584c->Draw("M>>hM2(601,-.01,6.01)","cut","egoff");
fb7d2d99 585
586TGraphErrors gr1;
587TH1 *h=c->GetHistogram();
588
589gr1->SetLineColor(h->GetLineColor());
590gr1->SetMarkerColor(h->GetMarkerColor());
591for (Int_t i=0;i<h->GetNbinsX();++i){
592 gr1.SetPoint(i,h->GetXaxis()->GetBinCenter(i+1),h->GetBinContent(i+1));
593// gr1.SetPointError(i,h->GetXaxis()->GetBinWidth(i+1)/2,h->GetBinError(i+1));
594 gr1.SetPointError(i,0,h->GetBinError(i+1));
595}
596
597//--------
598
599c->SetMarkerColor(kGreen);
600c->SetLineColor(kGreen);
601
1201a1a9 602c->Draw("M>>hM3(601,-.005,6.015)","cut","egoff");
fb7d2d99 603c->GetHistogram()->GetXaxis()->SetRangeUser(2.,4);
604
605TGraphErrors gr2;
606TH1 *h=c->GetHistogram();
607
608gr2->SetLineColor(h->GetLineColor());
609gr2->SetMarkerColor(h->GetMarkerColor());
610for (Int_t i=0;i<h->GetNbinsX();++i){
611 gr2.SetPoint(i,h->GetXaxis()->GetBinCenter(i+1),h->GetBinContent(i+1));
612// gr2.SetPointError(i,h->GetXaxis()->GetBinWidth(i+1)/2,h->GetBinError(i+1));
613 gr2.SetPointError(i,0,h->GetBinError(i+1));
614}
615
616//--------
617
618c->SetMarkerColor(kMagenta);
619c->SetLineColor(kMagenta);
620
1201a1a9 621c->Draw("M>>hM4(601,-.0,6.02)","cut","egoff");
fb7d2d99 622
623TGraphErrors gr3;
624TH1 *h=c->GetHistogram();
625
626gr3->SetLineColor(h->GetLineColor());
627gr3->SetMarkerColor(h->GetMarkerColor());
628for (Int_t i=0;i<h->GetNbinsX();++i){
629 gr3.SetPoint(i,h->GetXaxis()->GetBinCenter(i+1),h->GetBinContent(i+1));
630// gr3.SetPointError(i,h->GetXaxis()->GetBinWidth(i+1)/2,h->GetBinError(i+1));
631 gr3.SetPointError(i,0,h->GetBinError(i+1));
632}
633
634gr0->Draw("ap");
635gr0->GetHistogram()->GetXaxis()->SetRangeUser(2.,4);
636gr1->Draw("p");
637gr2->Draw("p");
638gr3->Draw("p");
8df8e382 639
640
641
642
ffbede40 643c->SetAlias("LegNclDiffIter1","(Leg1_NclsTPC-Leg1_NclsTPCiter1)>-1&&(Leg2_NclsTPC-Leg2_NclsTPCiter1)>-1")
644c->SetAlias("cut","PairType==1&&LegNclDiffIter1")
8df8e382 645// histos
646AliDielectronHistos h("h","h");
647h.AddClass("TPCsignal");
648h.UserHistogram("TPCsignal","sigTPC","TPC signal;P [GeV];TPC signal [arb. Units]",400,.3,40,400,0.,200.,0,0,kTRUE,kFALSE)
649h.GetHistogram("TPCsignal","sigTPC")->SetDirectory(gDirectory)
650
554e40f8 651h.UserHistogram("TPCsignal","nSigE","TPC n #sigma Electrons;P [GeV];TPC n #sigma Electrons",200,.3,40.,100,-10.,10.,0,0,kTRUE,kFALSE)
8df8e382 652h.GetHistogram("TPCsignal","nSigE")->SetDirectory(gDirectory)
554e40f8 653h.UserHistogram("TPCsignal","nSigMu","TPC n #sigma Muons;P [GeV];TPC n #sigma Muons",400,.3,40.,500,-10.,10.,0,0,kTRUE,kFALSE)
8df8e382 654h.GetHistogram("TPCsignal","nSigMu")->SetDirectory(gDirectory)
554e40f8 655h.UserHistogram("TPCsignal","nSigPi","TPC n #sigma Pions;P [GeV];TPC n #sigma Pions",400,.3,40.,500,-10.,10.,0,0,kTRUE,kFALSE)
8df8e382 656h.GetHistogram("TPCsignal","nSigPi")->SetDirectory(gDirectory)
554e40f8 657h.UserHistogram("TPCsignal","nSigK","TPC n #sigma Kaons;P [GeV];TPC n #sigma Kaons",400,.3,40,500,-10,10,0,0,kTRUE,kFALSE)
8df8e382 658h.GetHistogram("TPCsignal","nSigK")->SetDirectory(gDirectory)
554e40f8 659h.UserHistogram("TPCsignal","nSigP","TPC n #sigma Protons;P [GeV];TPC n #sigma Protons",400,.3,40.,500,-10,10.,0,0,kTRUE,kFALSE)
8df8e382 660h.GetHistogram("TPCsignal","nSigP")->SetDirectory(gDirectory)
661
ffbede40 662h.UserHistogram("TPCsignal","nSigDiffP","ncls-nclsXX;P [GeV];ncls-nclsXX",400,.3,40.,200,-40,160.,0,0,kTRUE,kFALSE)
663h.GetHistogram("TPCsignal","nSigDiffP")->SetDirectory(gDirectory)
664
8df8e382 665
666c->Draw("Leg1_TPC_signal:Leg1_P_InnerParam>>sigTPC","cut","colz")
667c->Draw("Leg2_TPC_signal:Leg2_P_InnerParam>>+sigTPC","cut","colz")
668
669
670c->Draw("Leg1_TPC_nSigma_Electrons:Leg1_P_InnerParam>>nSigE","cut","colz")
671c->Draw("Leg2_TPC_nSigma_Electrons:Leg2_P_InnerParam>>+nSigE","cut","colz")
672
673c->Draw("Leg1_TPC_nSigma_Muos:Leg1_P_InnerParam>>nSigMu","cut","goff")
674c->Draw("Leg2_TPC_nSigma_Muons:Leg1_P_InnerParam>>nSigMu","cut","goff")
675c->Draw("Leg2_TPC_nSigma_Muons:Leg2_P_InnerParam>>+nSigMu","cut","colz")
676
677c->Draw("Leg1_TPC_nSigma_Pions:Leg1_P_InnerParam>>nSigPi","cut","goff")
678c->Draw("Leg2_TPC_nSigma_Pions:Leg2_P_InnerParam>>+nSigPi","cut","colz")
679
680c->Draw("Leg1_TPC_nSigma_Kaons:Leg1_P_InnerParam>>nSigK","cut","goff")
681c->Draw("Leg2_TPC_nSigma_Kaons:Leg2_P_InnerParam>>+nSigK","cut","colz")
682
683c->Draw("Leg1_TPC_nSigma_Protons+.5:Leg1_P_InnerParam>>nSigP","cut","goff")
684c->Draw("Leg2_TPC_nSigma_Protons+.5:Leg2_P_InnerParam>>+nSigP","cut","colz")
685
686
687
61d106d3 688c->Draw("Leg1_TOF_nSigma_Electrons:Leg1_P_InnerParam>>nSigE","cut","colz")
689c->Draw("Leg2_TOF_nSigma_Electrons:Leg2_P_InnerParam>>+nSigE","cut","colz")
690
8df8e382 691
692c->Draw("Leg1_TOF_nSigma_Protons:Leg1_P_InnerParam>>nSigP","cut","goff")
693c->Draw("Leg2_TOF_nSigma_Protons:Leg2_P_InnerParam>>+nSigP","cut","colz")
694
695
696Pair->SetScanField(0)
697Pair->Scan("EventInFile:File.GetString()","","colsize=1 col=3.d:100.s")
698
699
700AliDielectronSignalFunc sig;
701sig.SetDefaults(1);
702
554e40f8 703
ffbede40 704//----------
705c->Draw("Leg1_NclsTPC>>hNcls(160,-0.5,159.5)","cut","goff");
706c->Draw("Leg2_NclsTPC>>+hNcls","cut","goff");
707hNclsPID->SetLineColor(kBlack)
708
709c->Draw("Leg1_TPCsignalN>>hNclsPID(160,-0.5,159.5)","cut","goff");
710c->Draw("Leg2_TPCsignalN>>+hNclsPID","cut","goff");
711hNclsPID->SetLineColor(kBlue)
712
713c->Draw("Leg1_NclsTPCiter1>>hNclsIter1(160,-0.5,159.5)","cut","goff");
714c->Draw("Leg2_NclsTPCiter1>>+hNclsIter1","cut","goff");
715hNclsIter1->SetLineColor(kGreen)
716
717hNcls->Draw();
718hNclsPID->Draw("same");
719hNclsIter1->Draw("same");
720//-----------
721
722
723
724c->Draw("Leg1_TPCsignalN:Leg1_NclsTPC>>hNclsPIDNcls(160,-0.5,159.5,160,-0.5,159.5)","cut","colz");
725c->Draw("Leg2_TPCsignalN:Leg2_NclsTPC>>+hNclsPIDNcls","cut","colz");
726
727
728c->Draw("Leg1_NclsTPC-Leg1_TPCsignalN:Leg1_P_InnerParam>>nSigDiffP","cut","colz");
729c->Draw("Leg2_NclsTPC-Leg2_TPCsignalN:Leg1_P_InnerParam>>+nSigDiffP","cut","colz");
730
731
732c->Draw("Leg1_NclsTPC-Leg1_NclsTPCiter1:Leg1_P_InnerParam>>nSigDiffP","cut","colz");
733c->Draw("Leg2_NclsTPC-Leg2_NclsTPCiter1:Leg1_P_InnerParam>>+nSigDiffP","cut","colz");
734
735
736
737
738
739
740
741
742
743c->Draw("Leg1_NclsTPC-Leg1_TPCsignalN:Leg1_TPC_nSigma_Electrons:Leg1_P_InnerParam>>hXX(100,0,10,20,-4,4)","cut","profcolz")
744c->Draw("Leg1_NclsTPC-Leg1_TPCsignalN:Leg2_TPC_nSigma_Electrons:Leg2_P_InnerParam>>+hXX","cut","profcolz")
745
746
747
748
749
750
751
752
554e40f8 753
8df8e382 754//WooJins cuts:
755c->SetAlias("cutE","abs(Leg1_TPC_nSigma_Electrons)<3&&abs(Leg2_TPC_nSigma_Electrons)<3");
554e40f8 756c->SetAlias("cutPi","abs(Leg1_TPC_nSigma_Pions)>5&&abs(Leg2_TPC_nSigma_Pions)>5");
8df8e382 757// c->SetAlias("cutPi","Leg1_TPC_signal>65&&Leg2_TPC_signal>65");
554e40f8 758c->SetAlias("cutP","(Leg1_TPC_nSigma_Protons)>5&&(Leg2_TPC_nSigma_Protons)>5");
8df8e382 759c->SetAlias("cutK","abs(Leg1_TPC_nSigma_Kaons)>2&&abs(Leg2_TPC_nSigma_Kaons)>2");
554e40f8 760c->SetAlias("pid","cutE&&cutPi&&cutP");
761c->SetAlias("etaLeg","abs(Leg1_Eta)<.9&&abs(Leg2_Eta)<.9");
762c->SetAlias("rap","abs(Y)<.9");
8df8e382 763
764
765c->SetAlias("cutAdd","PairType==1&&abs(Leg1_ImpactParXY)<.02&&abs(Leg2_ImpactParXY)<.02&&Leg2_Pt>1.")
554e40f8 766c->Draw("M>>hM(50,2,4)","PairType==1&&pid","e");
8df8e382 767h.Rebin();
768h.Rebin();
769h.Rebin();
770sig.Process(hM);
771sig.Draw("samestat");
772
773
774
775
776//test
8df8e382 777c->SetAlias("cutAdd","PairType==1&&abs(Leg1_ImpactParXY)<.03&&abs(Leg2_ImpactParXY)<.03&&Leg2_Pt>1")
554e40f8 778c->SetAlias("cut","cutAdd&&pid&etaLeg&&rap")
8df8e382 779
780c->Draw("M>>hM(50,1.99,3.99)","cut","e");
781sig.Process(hM);
782sig.Draw("samestat");
783
784
785c->SetAlias("cut","PairType==1")
786c->Draw("Leg1_TPC_signal:Leg1_P_InnerParam>>sigTPC","cut","colz")
787c->Draw("Leg2_TPC_signal:Leg2_P_InnerParam>>+sigTPC","cut","colz")
788
789////
790c->SetAlias("cutE","abs(Leg1_TPC_nSigma_Electrons)<2&&abs(Leg2_TPC_nSigma_Electrons)<2");
791c->SetAlias("cutPi","abs(Leg1_TPC_nSigma_Pions)>2&&abs(Leg2_TPC_nSigma_Pions)>2");
792c->SetAlias("cutPi2","Leg1_TPC_signal>65&&Leg2_TPC_signal>65");
793c->SetAlias("cutPi3","abs(Leg1_TPC_nSigma_Pions)>2.5&&abs(Leg2_TPC_nSigma_Pions)>2.5");
794c->SetAlias("cutP","abs(Leg1_TPC_nSigma_Protons)>2&&abs(Leg2_TPC_nSigma_Protons)>2");
795c->SetAlias("cutP2","abs(Leg1_TPC_nSigma_Protons)>1.5&&abs(Leg2_TPC_nSigma_Protons)>1.5");
796c->SetAlias("cutK","abs(Leg1_TPC_nSigma_Kaons)>2&&abs(Leg2_TPC_nSigma_Kaons)>2");
797c->SetAlias("cutdXY","abs(Leg1_ImpactParXY)<.02&&abs(Leg2_ImpactParXY)<.02")
798c->SetAlias("cutPt","Leg2_Pt>1.")
799c->SetAlias("cutPt2","Leg2_Pt>1.2")
800
801//----
802c->SetMarkerSize(0.7);
803
804c->SetMarkerStyle(20);
805c->SetLineColor(kRed);
806c->SetMarkerColor(kRed);
807c->Draw("M>>hM(100,2,4)","cutPi&&cutPt","e");
808
809c->SetMarkerStyle(21);
810c->SetLineColor(kRed-1);
811c->SetMarkerColor(kRed-2);
812c->Draw("M>>hM2(100,2,4)","cutPi2&&cutPt","esame");
813
814c->SetMarkerStyle(22);
815c->SetLineColor(kRed-2);
816c->SetMarkerColor(kRed-2);
817c->Draw("M>>hM3(100,2,4)","cutPi3&&cutPt","esame");
818
819//----
820c->SetMarkerStyle(20);
821c->SetLineColor(kBlue);
822c->SetMarkerColor(kBlue);
823c->Draw("M>>hM4(100,2,4)","cutPi&&cutPt&&cutP","esame");
824
825c->SetMarkerStyle(21);
826c->SetLineColor(kBlue-1);
827c->SetMarkerColor(kBlue-1);
828c->Draw("M>>hM5(100,2,4)","cutPi2&&cutPt&&cutP","esame");
829
830c->SetMarkerStyle(22);
831c->SetLineColor(kBlue-2);
832c->SetMarkerColor(kBlue-2);
833c->Draw("M>>hM6(100,2,4)","cutPi3&&cutPt&&cutP","esame");
834
835//----
836
837c->SetMarkerStyle(20);
838c->SetLineColor(kGreen);
839c->SetMarkerColor(kGreen);
840c->Draw("M>>hM7(100,2,4)","cutPi&&cutPt&&cutP2","esame");
841
842c->SetMarkerStyle(21);
843c->SetLineColor(kGreen-1);
844c->SetMarkerColor(kGreen-1);
845c->Draw("M>>hM8(100,2,4)","cutPi2&&cutPt&&cutP2","esame");
846
847c->SetMarkerStyle(22);
848c->SetLineColor(kGreen-2);
849c->SetMarkerColor(kGreen-2);
850c->Draw("M>>hM9(100,2,4)","cutPi3&&cutPt&&cutP2","esame");
851
852//----
853
854
855c->SetMarkerStyle(20);
856c->SetLineColor(kMagentha);
857c->SetMarkerColor(kMagentha);
858c->Draw("M>>hM7(100,2,4)","cutPi&&cutPt&&cutP2","esame");
859
860c->SetMarkerStyle(21);
861c->SetLineColor(kMagentha-1);
862c->SetMarkerColor(kMagentha-1);
863c->Draw("M>>hM8(100,2,4)","cutPi2&&cutPt&&cutP2","esame");
864
865c->SetMarkerStyle(22);
866c->SetLineColor(kMagentha-2);
867c->SetMarkerColor(kMagentha-2);
868c->Draw("M>>hM9(100,2,4)","cutPi3&&cutPt&&cutP2","esame");
869
870
871c->SetLineColor(kBlack);
872c->SetMarkerColor(kBlue);
873c->Draw("M>>hM4(100,2,4)","cutE&&cutPi&&cutK&&cutP&&cutdXY&&cutPt","esame");
874
875
876
877
878//
879c->SetAlias("cutE","Leg1_TPC_nSigma_Electrons>-1.5&&Leg2_TPC_nSigma_Electrons>-1.5");
880// c->SetAlias("cutE","Leg1_TPC_signal>60&&Leg2_TPC_signal>60");
881c->SetAlias("cutP","abs(Leg1_TPC_nSigma_Protons)>3&&abs(Leg2_TPC_nSigma_Protons)>3")
882
883c->SetAlias("cutAdd","PairType==1&&abs(Leg1_ImpactParXY)<.03&&abs(Leg2_ImpactParXY)<.03&&Leg2_Pt>0")
884c->SetAlias("cut","Leg2_Pt>1&&cutE&&cutP")
885
886c->Draw("M>>hM(50,1.99,3.99)","cut","e");
887
888c->SetAlias("cutAdd","PairType==1&&abs(Leg1_ImpactParXY)<.03&&abs(Leg2_ImpactParXY)<.03&&Leg2_Pt>.8")
889c->Draw("M>>hM2(50,1.99,3.99)","cut","esame");
890
891c->SetAlias("cutAdd","PairType==1&&abs(Leg1_ImpactParXY)<.03&&abs(Leg2_ImpactParXY)<.03&&Leg2_Pt>1")
892c->Draw("M>>hM3(50,1.99,3.99)","cut","esame");
893
894c->SetAlias("cutAdd","PairType==1&&abs(Leg1_ImpactParXY)<.03&&abs(Leg2_ImpactParXY)<.03&&Leg2_Pt>1.2")
895c->Draw("M>>hM4(50,1.99,3.99)","cut","esame");
896
897
898c->Draw("Leg1_TPC_signal:Leg1_P_InnerParam>>sigTPC","cut","goff")
899c->Draw("Leg2_TPC_signal:Leg2_P_InnerParam>>+sigTPC","cut","goff")
900c1->Modified();c1->Update()
901
902
903c->SetAlias("cutE","Leg1_TPC_nSigma_Electrons>-1.5&&Leg2_TPC_nSigma_Electrons>-1.5");
904c->SetAlias("cut","Leg2_P_InnerParam>1.5&&cutE")
905c->SetMarkerStyle(21);
906c->Draw("M>>hM(50,1.99,3.99)","cut","e");
907
908c->SetAlias("cutE","Leg1_TPC_nSigma_Electrons>-1.5+.8&&Leg2_TPC_nSigma_Electrons>-1.5+.8");
909c->SetAlias("cut","Leg2_P_InnerParam>1.5&&cutE")
910c->SetMarkerStyle(20);
911c->SetMarkerColor(kRed);
61d106d3 912c->Draw("M>>hM3(50,1.99,3.99)","cut","esame");
913
914
915
916
917//=================== pass 1 pass 2 comparison =======================
918
919
920//------ cuts ------------
921
922
923c->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");
924
925c->SetAlias("Pcut","Leg1_P_InnerParam>1.3&&Leg2_P_InnerParam>1.3")
926c->SetAlias("Ptcut","Leg1_Pt>1&&Leg2_Pt>1")
927c->SetAlias("TOFcut","abs(Leg1_TOF_nSigma_Electrons)<3&&abs(Leg2_TOF_nSigma_Electrons)<3");
928c->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)");
929c->SetAlias("TPCcut","(Leg1_TPC_signal>70*(1-exp(-1*(Leg1_P_InnerParam+2))))&&(Leg2_TPC_signal>70*(1-exp(-1*(Leg2_P_InnerParam+2))))")
2a14a7b1 930c->SetAlias("NClcut","Leg1_NclsTPC>90&&Leg2_NclsTPC>90");
61d106d3 931
932c->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)")
164bfb53 933c->SetAlias("cut","PairType==1&&eleParam&&Run<127719")
61d106d3 934c->SetAlias("cut","1==1")
935c->SetAlias("cut","NClcut")
936
937c->SetAlias("cut","TOFcut&&TPCcut&&NClcut")
938
939c->SetAlias("cut","TOFcut2&&TPCcut&&NClcut")
940
941c->SetAlias("cut","cutPro&&TPCcut&&NClcut")
942
943c->SetAlias("cut","Pcut&&TPCcut&&NClcut")
944
945c->SetAlias("cut","Ptcut&&TPCcut&&NClcut")
946
947
948
949//------------ plots --------------
950
951//no cut
952c->SetAlias("cut","1==1")
953c1->SetLogx(0)
954c1->SetLogz(0)
955c->SetLineColor(kBlack);
956c->Draw("M>>hM(50,1.99,3.99)","cut","e");
957hM->SetTitle(";Inv. Mass [GeV]; Pair (e+e-) / 40GeV")
958c1->Modified()
959c1->Update();
960c1->Print("pics/M_noCut.png");
961
962//=============
963//ncl: No cut
964//=============
965c->SetAlias("cut","1==1")
966c1->SetLogx(1)
967c1->SetLogz(0)
968gStyle->SetOptStat(0);
969c->Draw("Leg1_NclsTPC:Leg1_TPC_signal:Leg1_P_InnerParam>>test(1000,0,40,200,60,100)","cut","profcolz")
970c->Draw("Leg2_NclsTPC:Leg2_TPC_signal:Leg2_P_InnerParam>>+test","cut","profcolz")
971test->SetMinimum(80)
972test->SetTitle("mean TPC number of clusters;P_{TPC} [GeV]; TPC signal [arb. units]")
973c1->Modified()
974c1->Update();
975c1->Print("pics/TPCnCl_P.png");
976
977//=============
978//TPCsignal: ncl cut
979//=============
980c->SetAlias("cut","NClcut")
981c->Draw("Leg1_NclsTPC:Leg1_TPC_signal:Leg1_P_InnerParam>>test(1000,0,40,200,60,100)","cut","profcolz")
982c->Draw("Leg2_NclsTPC:Leg2_TPC_signal:Leg2_P_InnerParam>>+test","cut","profcolz")
983test->SetMinimum(80)
984test->SetTitle("mean TPC number of clusters;P_{TPC} [GeV]; TPC signal [arb. units]")
985c1->Modified()
986c1->Update();
987c1->Print("pics/TPCnCl_P_cutNcl.png");
988
989
990//=============
991//tpc signal + signal cut
992//=============
993c1->SetLogx(1)
994c1->SetLogy(0)
995c1->SetLogz(1)
996h.GetHistogram("TPCsignal","sigTPC")->GetYaxis()->SetRangeUser(60,100);
997c->SetAlias("cut","NClcut")
998c->Draw("Leg1_TPC_signal:Leg1_P_InnerParam>>sigTPC","cut","colz")
999c->Draw("Leg2_TPC_signal:Leg2_P_InnerParam>>+sigTPC","cut","colz")
1000TF1 f("f1","[0]*(1-exp(-[1]*(x-[2])))",0.3,40);
1001f.SetParameters(70,1,-2);
1002f.Draw("same");
1003c1->Modified();
1004c1->Update();
1005c1->Print("pics/TPCsignal_P_cutNcl.png");
1006
1007//------- Mass
1008
1009c1->SetLogx(0)
1010c1->SetLogy(1)
1011c1->SetLogz(0)
1012c->SetAlias("cut","1==1")
1013c->SetLineColor(kBlack);
1014c->SetMarkerColor(kBlack);
1015c->Draw("M>>hM(50,1.99,3.99)","cut","e");
1016hM->SetTitle(";Inv. Mass [GeV]; Pair (e+e-) / 40GeV")
1017hM->SetMinimum(5e2);
1018c->SetAlias("cut","NClcut")
1019c->SetLineColor(kRed);
1020c->SetMarkerColor(kRed);
8df8e382 1021c->Draw("M>>hM2(50,1.99,3.99)","cut","esame");
61d106d3 1022c1->Modified();
1023c1->Update();
1024c1->Print("pics/M_nClCut.png");
1025
1026
1027//==========
1028//tpc signal: ncl + tpc cut
1029//==========
1030c1->SetLogx(1)
1031c1->SetLogy(0)
1032c1->SetLogz(1)
1033c->SetAlias("cut","TPCcut&&NClcut")
1034c->Draw("Leg1_TPC_signal:Leg1_P_InnerParam>>sigTPC","cut","colz")
1035c->Draw("Leg2_TPC_signal:Leg2_P_InnerParam>>+sigTPC","cut","colz")
1036c1->Modified();
1037c1->Update();
1038c1->Print("pics/TPCsignal_P_cutNcl_tpc.png");
1039
1040/--- Mass
1041
1042c1->SetLogx(0)
1043c1->SetLogy(1)
1044c1->SetLogz(0)
1045c->SetAlias("cut","1==1")
1046c->SetLineColor(kBlack);
1047c->SetMarkerColor(kBlack);
1048c->Draw("M>>hM(50,1.99,3.99)","cut","e");
1049hM->SetTitle(";Inv. Mass [GeV]; Pair (e+e-) / 40GeV")
1050hM->SetMinimum(5);
1051c->SetAlias("cut","NClcut")
1052c->SetLineColor(kRed);
1053c->SetMarkerColor(kRed);
1054c->Draw("M>>hM2(50,1.99,3.99)","cut","esame");
1055c->SetAlias("cut","TPCcut&&NClcut")
1056c->SetLineColor(kGreen);
1057c->SetMarkerColor(kGreen);
1058c->Draw("M>>hM3(50,1.99,3.99)","cut","esame");
1059c1->Modified();
1060c1->Update();
1061c1->Print("pics/M_nClCut_tpc.png");
1062
1063
1064//========
1065//TPC signal: ncl + tpc + tof cut
1066//=======
1067c1->SetLogx(1)
1068c1->SetLogy(0)
1069c1->SetLogz(1)
1070c->SetAlias("cut","TOFcut2&&TPCcut&&NClcut")
1071c->Draw("Leg1_TPC_signal:Leg1_P_InnerParam>>sigTPC","cut","colz")
1072c->Draw("Leg2_TPC_signal:Leg2_P_InnerParam>>+sigTPC","cut","colz")
1073c1->Modified();
1074c1->Update();
1075c1->Print("pics/TPCsignal_P_cutNcl_tpc.png");
1076
1077//--- Mass
1078
1079c1->SetLogx(0)
1080c1->SetLogy(0)
1081c1->SetLogz(0)
1082c->SetAlias("cut","1==1")
1083c->SetAlias("cut","TPCcut&&NClcut")
1084c->SetLineColor(kGreen);
1085c->Draw("M>>hM(50,1.99,3.99)","cut","e");
1086hM->SetTitle(";Inv. Mass [GeV]; Pair (e+e-) / 40GeV")
1087hM->SetMinimum(.1);
1088c->SetAlias("cut","TOFcut2&&TPCcut&&NClcut")
1089c->SetLineColor(kBlue);
1090c->Draw("M>>hM2(50,1.99,3.99)","cut","esame");
1091c1->Modified();
1092c1->Update();
1093c1->Print("pics/M_nClCut_tpc.png");
1094
1095//========
1096//Inv Mass: different cuts
1097//=======
1098
1099c1->SetLogx(0)
1100c1->SetLogy(0)
1101c1->SetLogz(0)
1102c->SetAlias("cut","Ptcut&&TPCcut&&NClcut")
1103c->SetLineColor(kMagenta);
1104c->SetMarkerColor(kMagenta);
1105c->SetMarkerStyle(22);
1106c->Draw("M>>hM(50,1.99,3.99)","cut","e");
1107
1108c->SetAlias("cut","Pcut&&TPCcut&&NClcut")
1109c->SetLineColor(kCyan+1);
1110c->SetMarkerColor(kCyan+1);
1111c->SetMarkerStyle(21);
1112c->Draw("M>>hM2(50,1.99,3.99)","cut","esame");
1113
1114c->SetAlias("cut","TOFcut2&&TPCcut&&NClcut")
1115c->SetMarkerStyle(20);
1116c->SetLineColor(kBlue);
1117c->SetMarkerColor(kBlue);
1118c->Draw("M>>hM3(50,1.99,3.99)","cut","esame");
1119
1120c1->Modified();
1121c1->Update();
1122c1->Print("pics/M_nClCut_tpc_tof.png");
1123
1124
1125
8df8e382 1126
2a14a7b1 1127c->SetAlias("NClcut","Leg1_NclsTPC>90&&Leg2_NclsTPC>90");
1128c->SetAlias("Ptcut","Leg1_Pt>1&&Leg2_Pt>1")
1129c->SetAlias("PairT","PairType==1");
1130c->SetAlias("cut","NClcut&&Ptcut&&PairT")
1131
1132
1133c->SetAlias("cutE","abs(Leg1_TPC_nSigma_Electrons)<3&&abs(Leg2_TPC_nSigma_Electrons)<3");
1134c->SetAlias("cutPi","abs(Leg1_TPC_nSigma_Pions)>3&&abs(Leg2_TPC_nSigma_Pions)>3");
1135c->SetAlias("cutP","(Leg1_TPC_nSigma_Protons)>3&&(Leg2_TPC_nSigma_Protons)>3");
1136c->SetAlias("cut","NClcut&&Ptcut&&PairT&&cutE&&cutPi&&cutP")
1137
1138c->SetAlias("eta","abs(Leg1_Eta)<0.88&&abs(Leg1_Eta)<0.88");
1139c->SetAlias("rap","abs(Y)<0.88");
1140c->SetAlias("spdFirst","(Leg1_ITS_clusterMap&1)==1&&(Leg2_ITS_clusterMap&1)==1")
1141c->SetAlias("cut","NClcut&&Ptcut&&PairT&&cutE&&cutPi&&cutP&&eta&&rap&&spdFirst")
1142 c->Draw("M>>hM(50,1.98,1.98+50*.04)","cut","e")
1143
1144
1145c->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))
1146
1147c->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))
1148
1149c->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))
1150
1151c->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))
1152
1153c->SetAlias("cutPro","cutProL1&&cutProL2");
1154c->SetAlias("cutPio","cutPioL1&&cutPioL2");
8df8e382 1155*/
1156