PWG3/dielectron/macros/MakeDataReport.C

   1 void SetupStyle();
   2
   3 void 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
  54 void 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
  86 void 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
 203   Double_t res=5.2e-2;
 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
 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;
 217
 218
 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);
 222
 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);
 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
 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};
 332
 333   Double_t sigmas[5]={3,3,3,3,3};
 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;
 346   Double_t mip=50;
 347
 348   Color_t color=kRed;
 349   Int_t lineWidth=2;
 350
 351 TF1 *fBethe[5];
 352 TF1 *fBetheP[5];
 353 TF1 *fBetheM[5];
 354
 355 for (Int_t i=0; i<5; ++i){
 356   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);
 357   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);
 358   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);
 359
 360   fBethe[i]->SetParameters(alephParameters);
 361   fBetheP[i]->SetParameters(alephParameters);
 362   fBetheM[i]->SetParameters(alephParameters);
 363
 364   fBethe[i]->SetLineColor(baseColors[i]);
 365   fBetheP[i]->SetLineColor(baseColors[i]-sigmaColorOffset);
 366   fBetheM[i]->SetLineColor(baseColors[i]-sigmaColorOffset);
 367
 368   fBethe[i]->SetLineWidth(lineWidth);
 369   fBetheP[i]->SetLineWidth(lineWidth);
 370   fBetheM[i]->SetLineWidth(lineWidth);
 371
 372   fBetheP[i]->SetLineStyle(2);
 373   fBetheM[i]->SetLineStyle(2);
 374
 375   fBethe[i]->SetNpx(200);
 376   fBetheP[i]->SetNpx(200);
 377   fBetheM[i]->SetNpx(200);
 378 }
 379
 380 for (Int_t i=0; i<5; ++i){
 381   fBethe[i]->Draw("same");
 382 //   fBetheP[i]->Draw("same");
 383 //   fBetheM[i]->Draw("same");
 384 }
 385 }
 386
 387 */
 388
 389
 390 /*
 391
 392     Double_t resolution=0.052;
 393     Double_t nSigma=3.;
 394     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.);
 395     ffPio->SetParameters(0.0283086/0.97,2.63394e+01,5.04114e-11,2.12543e+00,4.88663e+00);
 396 fPPio->Draw("same");
 397
 398   TF1 f("fP","-8*exp(-0.6*x)",0,40);
 399 f.Draw("same")
 400
 401
 402 */
 403
 404
 405 /*
 406
 407
 408 c->SetAlias("cut","PairType==1")
 409
 410
 411 c->SetAlias("nCls","Leg1_NclsTPC>90&&Leg2_NclsTPC>90");
 412
 413 //--------PID
 414
 415 //-- nsigma
 416 c->SetAlias("cutE","abs(Leg1_TPC_nSigma_Electrons)<3&&abs(Leg2_TPC_nSigma_Electrons)<3");
 417 c->SetAlias("cutE","Leg1_TPC_nSigma_Electrons>-1 && Leg2_TPC_nSigma_Electrons>-1");
 418 c->SetAlias("cutPi","abs(Leg1_TPC_nSigma_Pions)>3&&abs(Leg2_TPC_nSigma_Pions)>3");
 419 c->SetAlias("cutP","(Leg1_TPC_nSigma_Protons)>3.&&(Leg2_TPC_nSigma_Protons)>3.");
 420 c->SetAlias("pidSig","cutE&&cutPi&&cutP");
 421
 422 //-- Pi param
 423 // 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)");
 424 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)");
 425 c->SetAlias("pidParam","eleParam&&cutP");
 426
 427
 428
 429 c->SetAlias("LegEta","abs(Leg1_Eta)<0.9&&abs(Leg2_Eta<0.9)");
 430 c->SetAlias("Rap","abs(Y)<0.9");
 431 c->SetAlias("QA","LegNcl&&LegEta&&Rap");
 432 c->SetAlias("spdFirst","(Leg1_ITS_clusterMap&1)==1 && (Leg2_ITS_clusterMap&1)==1");
 433
 434
 435 c->SetAlias("cut","PairType==1&&nCls&&pidSig&&LegEta&&Rap")
 436
 437 c->SetMarkerStyle(20);
 438 c->SetMarkerSize(.8);
 439 c->SetMarkerColor(kBlack);
 440 c->SetLineColor(kBlack);
 441
 442 // c->SetAlias("nCls","Leg1_NclsTPC>120&&Leg2_NclsTPC>120");
 443 c->Draw("M>>hM(50,2,4)","cut","e");
 444
 445 c->SetMarkerColor(kBlue);
 446 c->SetLineColor(kBlue);
 447 c->SetAlias("cut","PairType==1&&nCls&&pidParam&&LegEta&&Rap")
 448 // c->SetAlias("nCls","Leg1_NclsTPC>140&&Leg2_NclsTPC>140");
 449 c->Draw("M>>hM2(50,2,4)","cut","esame");
 450
 451 c->SetMarkerColor(kGreen);
 452 c->SetLineColor(kGreen);
 453 c->SetAlias("nCls","Leg1_NclsTPC>150&&Leg2_NclsTPC>150");
 454 c->Draw("M>>hM3(50,2,4)","cut","esame");
 455
 456
 457
 458
 459 c->SetAlias("cut","PairType==1&&nCls&&LegEta&&Rap")
 460 // histos
 461 AliDielectronHistos h("h","h");
 462 h.AddClass("TPCsignal");
 463 h.UserHistogram("TPCsignal","sigTPC","TPC signal;P [GeV];TPC signal [arb. Units]",400,.3,40,400,0.,200.,0,0,kTRUE,kFALSE)
 464 h.GetHistogram("TPCsignal","sigTPC")->SetDirectory(gDirectory)
 465
 466 h.UserHistogram("TPCsignal","nSigE","TPC n #sigma Electrons;P [GeV];TPC n #sigma Electrons",200,.3,40.,100,-10.,10.,0,0,kTRUE,kFALSE)
 467 h.GetHistogram("TPCsignal","nSigE")->SetDirectory(gDirectory)
 468 h.UserHistogram("TPCsignal","nSigMu","TPC n #sigma Muons;P [GeV];TPC n #sigma Muons",400,.3,40.,500,-10.,10.,0,0,kTRUE,kFALSE)
 469 h.GetHistogram("TPCsignal","nSigMu")->SetDirectory(gDirectory)
 470 h.UserHistogram("TPCsignal","nSigPi","TPC n #sigma Pions;P [GeV];TPC n #sigma Pions",400,.3,40.,500,-10.,10.,0,0,kTRUE,kFALSE)
 471 h.GetHistogram("TPCsignal","nSigPi")->SetDirectory(gDirectory)
 472 h.UserHistogram("TPCsignal","nSigK","TPC n #sigma Kaons;P [GeV];TPC n #sigma Kaons",400,.3,40,500,-10,10,0,0,kTRUE,kFALSE)
 473 h.GetHistogram("TPCsignal","nSigK")->SetDirectory(gDirectory)
 474 h.UserHistogram("TPCsignal","nSigP","TPC n #sigma Protons;P [GeV];TPC n #sigma Protons",400,.3,40.,500,-10,10.,0,0,kTRUE,kFALSE)
 475 h.GetHistogram("TPCsignal","nSigP")->SetDirectory(gDirectory)
 476
 477
 478 c->Draw("Leg1_TPC_signal:Leg1_P_InnerParam>>sigTPC","cut","colz")
 479 c->Draw("Leg2_TPC_signal:Leg2_P_InnerParam>>+sigTPC","cut","colz")
 480
 481
 482 c->Draw("Leg1_TPC_nSigma_Electrons:Leg1_P_InnerParam>>nSigE","cut","colz")
 483 c->Draw("Leg2_TPC_nSigma_Electrons:Leg2_P_InnerParam>>+nSigE","cut","colz")
 484
 485 c->Draw("Leg1_TPC_nSigma_Muos:Leg1_P_InnerParam>>nSigMu","cut","goff")
 486 c->Draw("Leg2_TPC_nSigma_Muons:Leg1_P_InnerParam>>nSigMu","cut","goff")
 487 c->Draw("Leg2_TPC_nSigma_Muons:Leg2_P_InnerParam>>+nSigMu","cut","colz")
 488
 489 c->Draw("Leg1_TPC_nSigma_Pions:Leg1_P_InnerParam>>nSigPi","cut","goff")
 490 c->Draw("Leg2_TPC_nSigma_Pions:Leg2_P_InnerParam>>+nSigPi","cut","colz")
 491
 492 c->Draw("Leg1_TPC_nSigma_Kaons:Leg1_P_InnerParam>>nSigK","cut","goff")
 493 c->Draw("Leg2_TPC_nSigma_Kaons:Leg2_P_InnerParam>>+nSigK","cut","colz")
 494
 495 c->Draw("Leg1_TPC_nSigma_Protons+.5:Leg1_P_InnerParam>>nSigP","cut","goff")
 496 c->Draw("Leg2_TPC_nSigma_Protons+.5:Leg2_P_InnerParam>>+nSigP","cut","colz")
 497
 498
 499
 500 c->Draw("Leg1_TOF_nSigma_Electrons:Leg1_P_InnerParam>>nSigE","cut","colz")
 501 c->Draw("Leg2_TOF_nSigma_Electrons:Leg2_P_InnerParam>>+nSigE","cut","colz")
 502
 503
 504 c->Draw("Leg1_TOF_nSigma_Protons:Leg1_P_InnerParam>>nSigP","cut","goff")
 505 c->Draw("Leg2_TOF_nSigma_Protons:Leg2_P_InnerParam>>+nSigP","cut","colz")
 506
 507
 508 Pair->SetScanField(0)
 509 Pair->Scan("EventInFile:File.GetString()","","colsize=1 col=3.d:100.s")
 510
 511
 512 AliDielectronSignalFunc sig;
 513 sig.SetDefaults(1);
 514
 515
 516
 517 //WooJins cuts:
 518 c->SetAlias("cutE","abs(Leg1_TPC_nSigma_Electrons)<3&&abs(Leg2_TPC_nSigma_Electrons)<3");
 519 c->SetAlias("cutPi","abs(Leg1_TPC_nSigma_Pions)>5&&abs(Leg2_TPC_nSigma_Pions)>5");
 520 // c->SetAlias("cutPi","Leg1_TPC_signal>65&&Leg2_TPC_signal>65");
 521 c->SetAlias("cutP","(Leg1_TPC_nSigma_Protons)>5&&(Leg2_TPC_nSigma_Protons)>5");
 522 c->SetAlias("cutK","abs(Leg1_TPC_nSigma_Kaons)>2&&abs(Leg2_TPC_nSigma_Kaons)>2");
 523 c->SetAlias("pid","cutE&&cutPi&&cutP");
 524 c->SetAlias("etaLeg","abs(Leg1_Eta)<.9&&abs(Leg2_Eta)<.9");
 525 c->SetAlias("rap","abs(Y)<.9");
 526
 527
 528 c->SetAlias("cutAdd","PairType==1&&abs(Leg1_ImpactParXY)<.02&&abs(Leg2_ImpactParXY)<.02&&Leg2_Pt>1.")
 529 c->Draw("M>>hM(50,2,4)","PairType==1&&pid","e");
 530 h.Rebin();
 531 h.Rebin();
 532 h.Rebin();
 533 sig.Process(hM);
 534 sig.Draw("samestat");
 535
 536
 537
 538
 539 //test
 540 c->SetAlias("cutAdd","PairType==1&&abs(Leg1_ImpactParXY)<.03&&abs(Leg2_ImpactParXY)<.03&&Leg2_Pt>1")
 541 c->SetAlias("cut","cutAdd&&pid&etaLeg&&rap")
 542
 543 c->Draw("M>>hM(50,1.99,3.99)","cut","e");
 544 sig.Process(hM);
 545 sig.Draw("samestat");
 546
 547
 548 c->SetAlias("cut","PairType==1")
 549 c->Draw("Leg1_TPC_signal:Leg1_P_InnerParam>>sigTPC","cut","colz")
 550 c->Draw("Leg2_TPC_signal:Leg2_P_InnerParam>>+sigTPC","cut","colz")
 551
 552 ////
 553 c->SetAlias("cutE","abs(Leg1_TPC_nSigma_Electrons)<2&&abs(Leg2_TPC_nSigma_Electrons)<2");
 554 c->SetAlias("cutPi","abs(Leg1_TPC_nSigma_Pions)>2&&abs(Leg2_TPC_nSigma_Pions)>2");
 555 c->SetAlias("cutPi2","Leg1_TPC_signal>65&&Leg2_TPC_signal>65");
 556 c->SetAlias("cutPi3","abs(Leg1_TPC_nSigma_Pions)>2.5&&abs(Leg2_TPC_nSigma_Pions)>2.5");
 557 c->SetAlias("cutP","abs(Leg1_TPC_nSigma_Protons)>2&&abs(Leg2_TPC_nSigma_Protons)>2");
 558 c->SetAlias("cutP2","abs(Leg1_TPC_nSigma_Protons)>1.5&&abs(Leg2_TPC_nSigma_Protons)>1.5");
 559 c->SetAlias("cutK","abs(Leg1_TPC_nSigma_Kaons)>2&&abs(Leg2_TPC_nSigma_Kaons)>2");
 560 c->SetAlias("cutdXY","abs(Leg1_ImpactParXY)<.02&&abs(Leg2_ImpactParXY)<.02")
 561 c->SetAlias("cutPt","Leg2_Pt>1.")
 562 c->SetAlias("cutPt2","Leg2_Pt>1.2")
 563
 564 //----
 565 c->SetMarkerSize(0.7);
 566
 567 c->SetMarkerStyle(20);
 568 c->SetLineColor(kRed);
 569 c->SetMarkerColor(kRed);
 570 c->Draw("M>>hM(100,2,4)","cutPi&&cutPt","e");
 571
 572 c->SetMarkerStyle(21);
 573 c->SetLineColor(kRed-1);
 574 c->SetMarkerColor(kRed-2);
 575 c->Draw("M>>hM2(100,2,4)","cutPi2&&cutPt","esame");
 576
 577 c->SetMarkerStyle(22);
 578 c->SetLineColor(kRed-2);
 579 c->SetMarkerColor(kRed-2);
 580 c->Draw("M>>hM3(100,2,4)","cutPi3&&cutPt","esame");
 581
 582 //----
 583 c->SetMarkerStyle(20);
 584 c->SetLineColor(kBlue);
 585 c->SetMarkerColor(kBlue);
 586 c->Draw("M>>hM4(100,2,4)","cutPi&&cutPt&&cutP","esame");
 587
 588 c->SetMarkerStyle(21);
 589 c->SetLineColor(kBlue-1);
 590 c->SetMarkerColor(kBlue-1);
 591 c->Draw("M>>hM5(100,2,4)","cutPi2&&cutPt&&cutP","esame");
 592
 593 c->SetMarkerStyle(22);
 594 c->SetLineColor(kBlue-2);
 595 c->SetMarkerColor(kBlue-2);
 596 c->Draw("M>>hM6(100,2,4)","cutPi3&&cutPt&&cutP","esame");
 597
 598 //----
 599
 600 c->SetMarkerStyle(20);
 601 c->SetLineColor(kGreen);
 602 c->SetMarkerColor(kGreen);
 603 c->Draw("M>>hM7(100,2,4)","cutPi&&cutPt&&cutP2","esame");
 604
 605 c->SetMarkerStyle(21);
 606 c->SetLineColor(kGreen-1);
 607 c->SetMarkerColor(kGreen-1);
 608 c->Draw("M>>hM8(100,2,4)","cutPi2&&cutPt&&cutP2","esame");
 609
 610 c->SetMarkerStyle(22);
 611 c->SetLineColor(kGreen-2);
 612 c->SetMarkerColor(kGreen-2);
 613 c->Draw("M>>hM9(100,2,4)","cutPi3&&cutPt&&cutP2","esame");
 614
 615 //----
 616
 617
 618 c->SetMarkerStyle(20);
 619 c->SetLineColor(kMagentha);
 620 c->SetMarkerColor(kMagentha);
 621 c->Draw("M>>hM7(100,2,4)","cutPi&&cutPt&&cutP2","esame");
 622
 623 c->SetMarkerStyle(21);
 624 c->SetLineColor(kMagentha-1);
 625 c->SetMarkerColor(kMagentha-1);
 626 c->Draw("M>>hM8(100,2,4)","cutPi2&&cutPt&&cutP2","esame");
 627
 628 c->SetMarkerStyle(22);
 629 c->SetLineColor(kMagentha-2);
 630 c->SetMarkerColor(kMagentha-2);
 631 c->Draw("M>>hM9(100,2,4)","cutPi3&&cutPt&&cutP2","esame");
 632
 633
 634 c->SetLineColor(kBlack);
 635 c->SetMarkerColor(kBlue);
 636 c->Draw("M>>hM4(100,2,4)","cutE&&cutPi&&cutK&&cutP&&cutdXY&&cutPt","esame");
 637
 638
 639
 640
 641 //
 642 c->SetAlias("cutE","Leg1_TPC_nSigma_Electrons>-1.5&&Leg2_TPC_nSigma_Electrons>-1.5");
 643 // c->SetAlias("cutE","Leg1_TPC_signal>60&&Leg2_TPC_signal>60");
 644 c->SetAlias("cutP","abs(Leg1_TPC_nSigma_Protons)>3&&abs(Leg2_TPC_nSigma_Protons)>3")
 645
 646 c->SetAlias("cutAdd","PairType==1&&abs(Leg1_ImpactParXY)<.03&&abs(Leg2_ImpactParXY)<.03&&Leg2_Pt>0")
 647 c->SetAlias("cut","Leg2_Pt>1&&cutE&&cutP")
 648
 649 c->Draw("M>>hM(50,1.99,3.99)","cut","e");
 650
 651 c->SetAlias("cutAdd","PairType==1&&abs(Leg1_ImpactParXY)<.03&&abs(Leg2_ImpactParXY)<.03&&Leg2_Pt>.8")
 652 c->Draw("M>>hM2(50,1.99,3.99)","cut","esame");
 653
 654 c->SetAlias("cutAdd","PairType==1&&abs(Leg1_ImpactParXY)<.03&&abs(Leg2_ImpactParXY)<.03&&Leg2_Pt>1")
 655 c->Draw("M>>hM3(50,1.99,3.99)","cut","esame");
 656
 657 c->SetAlias("cutAdd","PairType==1&&abs(Leg1_ImpactParXY)<.03&&abs(Leg2_ImpactParXY)<.03&&Leg2_Pt>1.2")
 658 c->Draw("M>>hM4(50,1.99,3.99)","cut","esame");
 659
 660
 661 c->Draw("Leg1_TPC_signal:Leg1_P_InnerParam>>sigTPC","cut","goff")
 662 c->Draw("Leg2_TPC_signal:Leg2_P_InnerParam>>+sigTPC","cut","goff")
 663 c1->Modified();c1->Update()
 664
 665
 666 c->SetAlias("cutE","Leg1_TPC_nSigma_Electrons>-1.5&&Leg2_TPC_nSigma_Electrons>-1.5");
 667 c->SetAlias("cut","Leg2_P_InnerParam>1.5&&cutE")
 668 c->SetMarkerStyle(21);
 669 c->Draw("M>>hM(50,1.99,3.99)","cut","e");
 670
 671 c->SetAlias("cutE","Leg1_TPC_nSigma_Electrons>-1.5+.8&&Leg2_TPC_nSigma_Electrons>-1.5+.8");
 672 c->SetAlias("cut","Leg2_P_InnerParam>1.5&&cutE")
 673 c->SetMarkerStyle(20);
 674 c->SetMarkerColor(kRed);
 675 c->Draw("M>>hM3(50,1.99,3.99)","cut","esame");
 676
 677
 678
 679
 680 //=================== pass 1 pass 2 comparison =======================
 681
 682
 683 //------ cuts ------------
 684
 685
 686 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");
 687
 688 c->SetAlias("Pcut","Leg1_P_InnerParam>1.3&&Leg2_P_InnerParam>1.3")
 689 c->SetAlias("Ptcut","Leg1_Pt>1&&Leg2_Pt>1")
 690 c->SetAlias("TOFcut","abs(Leg1_TOF_nSigma_Electrons)<3&&abs(Leg2_TOF_nSigma_Electrons)<3");
 691 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)");
 692 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))))")
 693 c->SetAlias("NClcut","Leg1_NclsTPC>90&&Leg2_NclsTPC>90");
 694
 695 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)")
 696 c->SetAlias("cut","PairType==1&&eleParam&&Run<127719")
 697 c->SetAlias("cut","1==1")
 698 c->SetAlias("cut","NClcut")
 699
 700 c->SetAlias("cut","TOFcut&&TPCcut&&NClcut")
 701
 702 c->SetAlias("cut","TOFcut2&&TPCcut&&NClcut")
 703
 704 c->SetAlias("cut","cutPro&&TPCcut&&NClcut")
 705
 706 c->SetAlias("cut","Pcut&&TPCcut&&NClcut")
 707
 708 c->SetAlias("cut","Ptcut&&TPCcut&&NClcut")
 709
 710
 711
 712 //------------ plots --------------
 713
 714 //no cut
 715 c->SetAlias("cut","1==1")
 716 c1->SetLogx(0)
 717 c1->SetLogz(0)
 718 c->SetLineColor(kBlack);
 719 c->Draw("M>>hM(50,1.99,3.99)","cut","e");
 720 hM->SetTitle(";Inv. Mass [GeV]; Pair (e+e-) / 40GeV")
 721 c1->Modified()
 722 c1->Update();
 723 c1->Print("pics/M_noCut.png");
 724
 725 //=============
 726 //ncl: No cut
 727 //=============
 728 c->SetAlias("cut","1==1")
 729 c1->SetLogx(1)
 730 c1->SetLogz(0)
 731 gStyle->SetOptStat(0);
 732 c->Draw("Leg1_NclsTPC:Leg1_TPC_signal:Leg1_P_InnerParam>>test(1000,0,40,200,60,100)","cut","profcolz")
 733 c->Draw("Leg2_NclsTPC:Leg2_TPC_signal:Leg2_P_InnerParam>>+test","cut","profcolz")
 734 test->SetMinimum(80)
 735 test->SetTitle("mean TPC number of clusters;P_{TPC} [GeV]; TPC signal [arb. units]")
 736 c1->Modified()
 737 c1->Update();
 738 c1->Print("pics/TPCnCl_P.png");
 739
 740 //=============
 741 //TPCsignal: ncl cut
 742 //=============
 743 c->SetAlias("cut","NClcut")
 744 c->Draw("Leg1_NclsTPC:Leg1_TPC_signal:Leg1_P_InnerParam>>test(1000,0,40,200,60,100)","cut","profcolz")
 745 c->Draw("Leg2_NclsTPC:Leg2_TPC_signal:Leg2_P_InnerParam>>+test","cut","profcolz")
 746 test->SetMinimum(80)
 747 test->SetTitle("mean TPC number of clusters;P_{TPC} [GeV]; TPC signal [arb. units]")
 748 c1->Modified()
 749 c1->Update();
 750 c1->Print("pics/TPCnCl_P_cutNcl.png");
 751
 752
 753 //=============
 754 //tpc signal + signal cut
 755 //=============
 756 c1->SetLogx(1)
 757 c1->SetLogy(0)
 758 c1->SetLogz(1)
 759 h.GetHistogram("TPCsignal","sigTPC")->GetYaxis()->SetRangeUser(60,100);
 760 c->SetAlias("cut","NClcut")
 761 c->Draw("Leg1_TPC_signal:Leg1_P_InnerParam>>sigTPC","cut","colz")
 762 c->Draw("Leg2_TPC_signal:Leg2_P_InnerParam>>+sigTPC","cut","colz")
 763 TF1 f("f1","[0]*(1-exp(-[1]*(x-[2])))",0.3,40);
 764 f.SetParameters(70,1,-2);
 765 f.Draw("same");
 766 c1->Modified();
 767 c1->Update();
 768 c1->Print("pics/TPCsignal_P_cutNcl.png");
 769
 770 //------- Mass
 771
 772 c1->SetLogx(0)
 773 c1->SetLogy(1)
 774 c1->SetLogz(0)
 775 c->SetAlias("cut","1==1")
 776 c->SetLineColor(kBlack);
 777 c->SetMarkerColor(kBlack);
 778 c->Draw("M>>hM(50,1.99,3.99)","cut","e");
 779 hM->SetTitle(";Inv. Mass [GeV]; Pair (e+e-) / 40GeV")
 780 hM->SetMinimum(5e2);
 781 c->SetAlias("cut","NClcut")
 782 c->SetLineColor(kRed);
 783 c->SetMarkerColor(kRed);
 784 c->Draw("M>>hM2(50,1.99,3.99)","cut","esame");
 785 c1->Modified();
 786 c1->Update();
 787 c1->Print("pics/M_nClCut.png");
 788
 789
 790 //==========
 791 //tpc signal: ncl + tpc cut
 792 //==========
 793 c1->SetLogx(1)
 794 c1->SetLogy(0)
 795 c1->SetLogz(1)
 796 c->SetAlias("cut","TPCcut&&NClcut")
 797 c->Draw("Leg1_TPC_signal:Leg1_P_InnerParam>>sigTPC","cut","colz")
 798 c->Draw("Leg2_TPC_signal:Leg2_P_InnerParam>>+sigTPC","cut","colz")
 799 c1->Modified();
 800 c1->Update();
 801 c1->Print("pics/TPCsignal_P_cutNcl_tpc.png");
 802
 803 /--- Mass
 804
 805 c1->SetLogx(0)
 806 c1->SetLogy(1)
 807 c1->SetLogz(0)
 808 c->SetAlias("cut","1==1")
 809 c->SetLineColor(kBlack);
 810 c->SetMarkerColor(kBlack);
 811 c->Draw("M>>hM(50,1.99,3.99)","cut","e");
 812 hM->SetTitle(";Inv. Mass [GeV]; Pair (e+e-) / 40GeV")
 813 hM->SetMinimum(5);
 814 c->SetAlias("cut","NClcut")
 815 c->SetLineColor(kRed);
 816 c->SetMarkerColor(kRed);
 817 c->Draw("M>>hM2(50,1.99,3.99)","cut","esame");
 818 c->SetAlias("cut","TPCcut&&NClcut")
 819 c->SetLineColor(kGreen);
 820 c->SetMarkerColor(kGreen);
 821 c->Draw("M>>hM3(50,1.99,3.99)","cut","esame");
 822 c1->Modified();
 823 c1->Update();
 824 c1->Print("pics/M_nClCut_tpc.png");
 825
 826
 827 //========
 828 //TPC signal: ncl + tpc +  tof cut
 829 //=======
 830 c1->SetLogx(1)
 831 c1->SetLogy(0)
 832 c1->SetLogz(1)
 833 c->SetAlias("cut","TOFcut2&&TPCcut&&NClcut")
 834 c->Draw("Leg1_TPC_signal:Leg1_P_InnerParam>>sigTPC","cut","colz")
 835 c->Draw("Leg2_TPC_signal:Leg2_P_InnerParam>>+sigTPC","cut","colz")
 836 c1->Modified();
 837 c1->Update();
 838 c1->Print("pics/TPCsignal_P_cutNcl_tpc.png");
 839
 840 //--- Mass
 841
 842 c1->SetLogx(0)
 843 c1->SetLogy(0)
 844 c1->SetLogz(0)
 845 c->SetAlias("cut","1==1")
 846 c->SetAlias("cut","TPCcut&&NClcut")
 847 c->SetLineColor(kGreen);
 848 c->Draw("M>>hM(50,1.99,3.99)","cut","e");
 849 hM->SetTitle(";Inv. Mass [GeV]; Pair (e+e-) / 40GeV")
 850 hM->SetMinimum(.1);
 851 c->SetAlias("cut","TOFcut2&&TPCcut&&NClcut")
 852 c->SetLineColor(kBlue);
 853 c->Draw("M>>hM2(50,1.99,3.99)","cut","esame");
 854 c1->Modified();
 855 c1->Update();
 856 c1->Print("pics/M_nClCut_tpc.png");
 857
 858 //========
 859 //Inv Mass: different cuts
 860 //=======
 861
 862 c1->SetLogx(0)
 863 c1->SetLogy(0)
 864 c1->SetLogz(0)
 865 c->SetAlias("cut","Ptcut&&TPCcut&&NClcut")
 866 c->SetLineColor(kMagenta);
 867 c->SetMarkerColor(kMagenta);
 868 c->SetMarkerStyle(22);
 869 c->Draw("M>>hM(50,1.99,3.99)","cut","e");
 870
 871 c->SetAlias("cut","Pcut&&TPCcut&&NClcut")
 872 c->SetLineColor(kCyan+1);
 873 c->SetMarkerColor(kCyan+1);
 874 c->SetMarkerStyle(21);
 875 c->Draw("M>>hM2(50,1.99,3.99)","cut","esame");
 876
 877 c->SetAlias("cut","TOFcut2&&TPCcut&&NClcut")
 878 c->SetMarkerStyle(20);
 879 c->SetLineColor(kBlue);
 880 c->SetMarkerColor(kBlue);
 881 c->Draw("M>>hM3(50,1.99,3.99)","cut","esame");
 882
 883 c1->Modified();
 884 c1->Update();
 885 c1->Print("pics/M_nClCut_tpc_tof.png");
 886
 887
 888
 889
 890 c->SetAlias("NClcut","Leg1_NclsTPC>90&&Leg2_NclsTPC>90");
 891 c->SetAlias("Ptcut","Leg1_Pt>1&&Leg2_Pt>1")
 892 c->SetAlias("PairT","PairType==1");
 893 c->SetAlias("cut","NClcut&&Ptcut&&PairT")
 894
 895
 896 c->SetAlias("cutE","abs(Leg1_TPC_nSigma_Electrons)<3&&abs(Leg2_TPC_nSigma_Electrons)<3");
 897 c->SetAlias("cutPi","abs(Leg1_TPC_nSigma_Pions)>3&&abs(Leg2_TPC_nSigma_Pions)>3");
 898 c->SetAlias("cutP","(Leg1_TPC_nSigma_Protons)>3&&(Leg2_TPC_nSigma_Protons)>3");
 899 c->SetAlias("cut","NClcut&&Ptcut&&PairT&&cutE&&cutPi&&cutP")
 900
 901 c->SetAlias("eta","abs(Leg1_Eta)<0.88&&abs(Leg1_Eta)<0.88");
 902 c->SetAlias("rap","abs(Y)<0.88");
 903 c->SetAlias("spdFirst","(Leg1_ITS_clusterMap&1)==1&&(Leg2_ITS_clusterMap&1)==1")
 904 c->SetAlias("cut","NClcut&&Ptcut&&PairT&&cutE&&cutPi&&cutP&&eta&&rap&&spdFirst")
 905  c->Draw("M>>hM(50,1.98,1.98+50*.04)","cut","e")
 906
 907
 908 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))
 909
 910 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))
 911
 912 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))
 913
 914 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))
 915
 916 c->SetAlias("cutPro","cutProL1&&cutProL2");
 917 c->SetAlias("cutPio","cutPioL1&&cutPioL2");
 918 */
 919