36 #include <TGraphAsymmErrors.h>
63 fCentralityClass(
"0100"),
64 fRapidityRange(
"0101"),
66 fIsLowPtAnalysis(false),
67 fIsPass4Analysis(false),
69 fIsRapidityScan(false)
112 }
else AliFatal(
"Not yet implemented");
123 else AliFatal(
"Not yet implemented");
137 else AliFatal(
"Not yet implemented");
168 else AliFatal(
"Not yet implemented");
180 }
else AliFatal(
"Not yet implemented");
192 else AliFatal(
"Not yet implemented");
204 else AliFatal(
"Not yet implemented");
235 else AliFatal(
"Not yet implemented");
246 }
else AliFatal(
"Not yet implemented");
258 else AliFatal(
"Not yet implemented");
270 else AliFatal(
"Not yet implemented");
299 else AliFatal(
"Not yet implemented");
308 else AliFatal(
"Not yet implemented");
313 else AliFatal(
"Not yet implemented");
317 else AliFatal(
"Not yet implemented");
322 else AliFatal(
"Not yet implemented");
326 printf(
"Invalid decay type: %d\n",decay);
340 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
341 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.035);
344 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
345 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.012);
348 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
349 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.08);
352 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
356 for(Int_t i=4;i<=7;i++)
fRawYield->SetBinContent(i,0.04);
357 for(Int_t i=8;i<=12;i++)
fRawYield->SetBinContent(i,0.07);
358 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.10);
359 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,1);
362 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
363 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.10);
366 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
367 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.03);
368 fPIDEff->SetBinContent(2,0.05);
371 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
372 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0);
391 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
392 for(Int_t i=1;i<=20;i++)
fNorm->SetBinContent(i,0.05);
395 fBR =
new TH1F(
"fBR",
"fBR",20,0,20);
396 for(Int_t i=1;i<=20;i++)
fBR->SetBinContent(i,0.012);
399 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",20,0,20);
400 for(Int_t i=1;i<=20;i++)
fTrackingEff->SetBinContent(i,0.10);
403 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
407 for(Int_t i=4;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
408 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.10);
411 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
416 for(Int_t i=5;i<=16;i++)
fCutsEff->SetBinContent(i,0.10);
417 for(Int_t i=17;i<=20;i++)
fCutsEff->SetBinContent(i,0.);
420 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
421 for(Int_t i=3;i<=16;i++)
fPIDEff->SetBinContent(i,0.05);
424 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",20,0,20);
425 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0.01);
428 for(Int_t i=13;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
449 fNorm =
new TH1F(
"fNorm",
"fNorm",36,0,36);
450 for(Int_t i=1;i<36;i++)
fNorm->SetBinContent(i,0.048);
453 fBR =
new TH1F(
"fBR",
"fBR",36,0,36);
454 for(Int_t i=1;i<=36;i++)
fBR->SetBinContent(i,0.012);
457 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",36,0,36);
458 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.10);
459 for(Int_t i=25;i<=36;i++)
fTrackingEff->SetBinContent(i,0.);
462 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",36,0,36);
466 for(Int_t i=4;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
467 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.10);
468 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.30);
469 for(Int_t i=25;i<=36;i++)
fRawYield->SetBinContent(i,0.);
472 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",36,0,36);
478 for(Int_t i=6;i<=24;i++)
fCutsEff->SetBinContent(i,0.06);
479 for(Int_t i=25;i<=36;i++)
fCutsEff->SetBinContent(i,0.0);
482 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",36,0,36);
483 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.05);
484 for(Int_t i=25;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
487 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",36,0,36);
488 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
492 for(Int_t i=13;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
493 for(Int_t i=17;i<=24;i++)
fMCPtShape->SetBinContent(i,0.05);
494 for(Int_t i=25;i<=36;i++)
fMCPtShape->SetBinContent(i,0.0);
520 for(Int_t i=5;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
521 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
522 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.15);
523 for(Int_t i=17;i<=36;i++)
fRawYield->SetBinContent(i,0.);
533 for(Int_t i=7;i<=8;i++)
fCutsEff->SetBinContent(i,0.15);
534 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.15);
535 for(Int_t i=25;i<=36;i++)
fCutsEff->SetBinContent(i,0.0);
550 for(Int_t i=4;i<=6;i++)
fRawYield->SetBinContent(i,0.07);
551 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
552 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.10);
553 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.15);
554 for(Int_t i=17;i<=36;i++)
fRawYield->SetBinContent(i,0.);
564 for(Int_t i=7;i<=8;i++)
fCutsEff->SetBinContent(i,0.15);
565 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.15);
566 for(Int_t i=17;i<=36;i++)
fCutsEff->SetBinContent(i,0.0);
577 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
578 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.07);
581 fBR =
new TH1F(
"fBR",
"fBR",20,0,20);
582 for(Int_t i=1;i<=20;i++)
fBR->SetBinContent(i,0.012);
585 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",20,0,20);
586 for(Int_t i=1;i<=20;i++)
fTrackingEff->SetBinContent(i,0.10);
590 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
593 for(Int_t i=3;i<=16;i++)
fRawYield->SetBinContent(i,0.05);
597 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
602 for(Int_t i=5;i<=16;i++)
fCutsEff->SetBinContent(i,0.10);
603 for(Int_t i=17;i<=20;i++)
fCutsEff->SetBinContent(i,0.);
606 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
609 for(Int_t i=3;i<=16;i++)
fPIDEff->SetBinContent(i,0.05);
612 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",20,0,20);
613 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0.01);
616 for(Int_t i=13;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
632 AliInfo(
" Settings for D0 --> Kpi Pb-Pb collisions at 2.76 TeV - 0-10 centrality - newww!!");
635 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
636 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.048);
639 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
640 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.012);
643 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
644 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.10);
648 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
655 for (Int_t i=7; i<=8; i++)
fRawYield->SetBinContent(i,0.06);
656 for (Int_t i=9; i<=12; i++)
fRawYield->SetBinContent(i,0.07);
657 for (Int_t i=13; i<=16; i++)
fRawYield->SetBinContent(i,0.11);
658 for (Int_t i=17; i<=24; i++)
fRawYield->SetBinContent(i,0.30);
660 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
667 for (Int_t i=7; i<=8; i++)
fCutsEff->SetBinContent(i,0.05);
668 for (Int_t i=9; i<=12; i++)
fCutsEff->SetBinContent(i,0.05);
669 for (Int_t i=13; i<=16; i++)
fCutsEff->SetBinContent(i,0.05);
670 for (Int_t i=17; i<=24; i++)
fCutsEff->SetBinContent(i,0.05);
674 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
677 for(Int_t i=2;i<=16;i++)
fPIDEff->SetBinContent(i,0.05);
680 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
685 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.);
702 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
703 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.07);
706 fBR =
new TH1F(
"fBR",
"fBR",20,0,20);
707 for(Int_t i=1;i<=20;i++)
fBR->SetBinContent(i,0.012);
710 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",20,0,20);
711 for(Int_t i=1;i<=20;i++)
fTrackingEff->SetBinContent(i,0.10);
715 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
729 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.112);
730 for(Int_t i=17;i<=20;i++)
fRawYield->SetBinContent(i,0);
732 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
739 for(Int_t i=7;i<=8;i++)
fCutsEff->SetBinContent(i,0.15);
740 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.15);
744 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
747 for(Int_t i=2;i<=16;i++)
fPIDEff->SetBinContent(i,0.05);
750 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",0,0,20);
751 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0.01);
755 for(Int_t i=13;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
756 for(Int_t i=17;i<=20;i++)
fMCPtShape->SetBinContent(i,0.);
771 AliInfo(
" Settings for D0 --> K pi, p-p collisions at 2.76 TeV");
774 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
775 for(Int_t i=1;i<=20;i++)
fNorm->SetBinContent(i,0.019);
778 fBR =
new TH1F(
"fBR",
"fBR",20,0,20);
779 for(Int_t i=1;i<=20;i++)
fBR->SetBinContent(i,0.012);
782 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",20,0,20);
783 for(Int_t i=1;i<=20;i++)
fTrackingEff->SetBinContent(i,0.10);
786 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
788 for(Int_t i=1;i<=20;i++)
fRawYield->SetBinContent(i,0.15);
791 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
792 for(Int_t i=1;i<=20;i++)
fCutsEff->SetBinContent(i,0.10);
794 for(Int_t i=7;i<=20;i++)
fCutsEff->SetBinContent(i,0.15);
798 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
799 for(Int_t i=1;i<=20;i++)
fPIDEff->SetBinContent(i,0.15);
801 for(Int_t i=7;i<=20;i++)
fPIDEff->SetBinContent(i,0.05);
804 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",20,0,20);
805 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0.01);
825 AliInfo(
" Settings for D0 --> K pi, pp collisions at 7 TeV, analysis without topological cuts");
828 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
829 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.035);
832 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
833 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.0129);
836 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
837 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.08);
840 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
841 for(Int_t i=1;i<=24;i++)
fRawYield->SetBinContent(i,0.2);
850 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.2);
853 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
854 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.);
857 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
858 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.03);
859 fPIDEff->SetBinContent(1,0.05);
860 fPIDEff->SetBinContent(2,0.05);
863 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
864 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0);
877 AliInfo(
" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV");
879 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
880 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.02);
883 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
884 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.013);
887 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
888 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.06);
891 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
894 for(Int_t i=3;i<=12;i++)
fRawYield->SetBinContent(i,0.04);
895 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.07);
896 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.11);
899 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
902 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
905 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
906 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
909 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
912 for(Int_t i=3;i<=24;i++)
fMCPtShape->SetBinContent(i,0.0);
930 AliInfo(
" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV, analysis without topological cuts");
933 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
934 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.035);
937 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
938 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.0129);
941 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
942 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.06);
945 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
946 for(Int_t i=1;i<=24;i++)
fRawYield->SetBinContent(i,0.2);
955 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.11);
958 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
959 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.);
962 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
963 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.);
966 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
967 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0);
987 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
988 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.035);
991 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
992 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.021);
995 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
996 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.12);
1000 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
1007 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
1008 for(Int_t i=12;i<=24;i++)
fRawYield->SetBinContent(i,0.10);
1011 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
1012 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.10);
1015 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
1016 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.05);
1017 fPIDEff->SetBinContent(1,0.15);
1018 fPIDEff->SetBinContent(2,0.15);
1019 fPIDEff->SetBinContent(3,0.15);
1020 fPIDEff->SetBinContent(4,0.15);
1021 for(Int_t i=12;i<=16;i++)
fPIDEff->SetBinContent(i,0.10);
1024 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
1025 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0);
1050 fNorm =
new TH1F(
"fNorm",
"fNorm",36,0,36);
1051 for(Int_t i=1;i<=36;i++)
fNorm->SetBinContent(i,0.035);
1054 fBR =
new TH1F(
"fBR",
"fBR",36,0,36);
1055 for(Int_t i=1;i<=36;i++)
fBR->SetBinContent(i,0.021);
1058 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",36,0,36);
1059 for(Int_t i=1;i<=36;i++)
fTrackingEff->SetBinContent(i,0.12);
1063 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",36,0,36);
1067 for(Int_t i=4;i<=12;i++)
fRawYield->SetBinContent(i,0.03);
1068 for(Int_t i=13;i<=24;i++)
fRawYield->SetBinContent(i,0.07);
1069 for(Int_t i=25;i<=36;i++)
fRawYield->SetBinContent(i,0.20);
1072 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",36,0,36);
1073 for(Int_t i=1;i<=36;i++)
fCutsEff->SetBinContent(i,0.05);
1076 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",36,0,36);
1077 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.05);
1080 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",36,0,36);
1081 for(Int_t i=1;i<=36;i++)
fMCPtShape->SetBinContent(i,0);
1105 fNorm =
new TH1F(
"fNorm",
"fNorm",36,0,36);
1106 for(Int_t i=1;i<=36;i++)
fNorm->SetBinContent(i,0.05);
1109 fBR =
new TH1F(
"fBR",
"fBR",36,0,36);
1110 for(Int_t i=1;i<=36;i++)
fBR->SetBinContent(i,0.025);
1113 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",36,0,36);
1114 for(Int_t i=1;i<=36;i++)
fTrackingEff->SetBinContent(i,0.09);
1117 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",36,0,36);
1121 for(Int_t i=4;i<=36;i++)
fRawYield->SetBinContent(i,0.04);
1124 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",36,0,36);
1127 for(Int_t i=3;i<=36;i++)
fCutsEff->SetBinContent(i,0.05);
1130 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",36,0,36);
1131 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.00);
1134 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",36,0,36);
1135 for(Int_t i=1;i<=36;i++)
fMCPtShape->SetBinContent(i,0);
1136 for(Int_t i=2;i<=3;i++)
fMCPtShape->SetBinContent(i,0.02);
1137 for(Int_t i=4;i<=5;i++)
fMCPtShape->SetBinContent(i,0.01);
1152 fNorm =
new TH1F(
"fNorm",
"fNorm",12,0,12);
1153 for(Int_t i=1;i<=12;i++)
fNorm->SetBinContent(i,0.035);
1156 fBR =
new TH1F(
"fBR",
"fBR",12,0,12);
1157 for(Int_t i=1;i<=12;i++)
fBR->SetBinContent(i,0.06);
1160 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",12,0,12);
1161 for(Int_t i=1;i<=12;i++)
fTrackingEff->SetBinContent(i,0.12);
1165 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",12,0,12);
1180 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",12,0,12);
1181 for(Int_t i=1;i<=12;i++)
fCutsEff->SetBinContent(i,0.15);
1184 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",12,0,12);
1185 for(Int_t i=1;i<=12;i++)
fPIDEff->SetBinContent(i,0.07);
1188 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",12,0,12);
1189 for(Int_t i=1; i<=2; i++)
fMCPtShape->SetBinContent(i,1.);
1190 for(Int_t i=3; i<=4; i++)
fMCPtShape->SetBinContent(i,0.03);
1191 for(Int_t i=5; i<=6; i++)
fMCPtShape->SetBinContent(i,0.03);
1192 for(Int_t i=7; i<=8; i++)
fMCPtShape->SetBinContent(i,0.02);
1193 for(Int_t i=9; i<=12; i++)
fMCPtShape->SetBinContent(i,0.02);
1216 fNorm =
new TH1F(
"fNorm",
"fNorm",36,0,36);
1217 for(Int_t i=1;i<=36;i++)
fNorm->SetBinContent(i,0.048);
1220 fBR =
new TH1F(
"fBR",
"fBR",36,0,36);
1221 for(Int_t i=1;i<=36;i++)
fBR->SetBinContent(i,0.021);
1224 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",36,0,36);
1225 for(Int_t i=1;i<=36;i++)
fTrackingEff->SetBinContent(i,0.15);
1228 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",36,0,36);
1229 for(Int_t i=1;i<=36;i++)
fRawYield->SetBinContent(i,.10);
1233 for(Int_t i=7; i<=8; i++)
fRawYield->SetBinContent(i,0.10);
1234 for(Int_t i=9; i<=12; i++)
fRawYield->SetBinContent(i,0.08);
1235 for(Int_t i=13; i<=16; i++)
fRawYield->SetBinContent(i,0.05);
1236 for(Int_t i=17; i<=24; i++)
fRawYield->SetBinContent(i,0.08);
1237 for(Int_t i=25; i<=36; i++)
fRawYield->SetBinContent(i,0.20);
1241 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",36,0,36);
1242 for(Int_t i=1; i<=12; i++)
fCutsEff->SetBinContent(i,0.10);
1243 for(Int_t i=13; i<=36; i++)
fCutsEff->SetBinContent(i,0.08);
1247 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",36,0,36);
1248 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.05);
1251 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",36,0,36);
1252 for(Int_t iBin=1; iBin<=5; iBin++)
fMCPtShape->SetBinContent(iBin,0.05);
1253 for(Int_t iBin=6; iBin<=36; iBin++)
fMCPtShape->SetBinContent(iBin,0.03);
1276 fNorm =
new TH1F(
"fNorm",
"fNorm",36,0,36);
1277 for(Int_t i=1;i<=36;i++)
fNorm->SetBinContent(i,0.048);
1280 fBR =
new TH1F(
"fBR",
"fBR",36,0,36);
1281 for(Int_t i=1;i<=36;i++)
fBR->SetBinContent(i,0.021);
1284 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",36,0,36);
1285 for(Int_t i=1;i<=36;i++)
fTrackingEff->SetBinContent(i,0.15);
1288 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",36,0,36);
1289 for(Int_t i=1;i<=3;i++)
fRawYield->SetBinContent(i,0.);
1293 for(Int_t i=7; i<=8; i++)
fRawYield->SetBinContent(i,0.08);
1294 for(Int_t i=9; i<=12; i++)
fRawYield->SetBinContent(i,0.08);
1295 for(Int_t i=13; i<=16; i++)
fRawYield->SetBinContent(i,0.08);
1296 for(Int_t i=17; i<=24; i++)
fRawYield->SetBinContent(i,0.08);
1297 for(Int_t i=25; i<=36; i++)
fRawYield->SetBinContent(i,0.08);
1301 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",36,0,36);
1302 for(Int_t i=1; i<=36; i++)
fCutsEff->SetBinContent(i,0.10);
1306 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",36,0,36);
1307 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.05);
1310 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",36,0,36);
1311 for(Int_t i=1;i<=3;i++)
fMCPtShape->SetBinContent(i,0.);
1314 for (Int_t i=6; i<=36; i++)
fMCPtShape->SetBinContent(i,0.01);
1336 AliInfo(
" Settings for D+ --> K pi pi, Pb-Pb collisions for 30-50 CC");
1338 fNorm =
new TH1F(
"fNorm",
"fNorm",16,0,16);
1339 for(Int_t i=1;i<=16;i++)
fNorm->SetBinContent(i,0.048);
1342 fBR =
new TH1F(
"fBR",
"fBR",16,0,16);
1343 for(Int_t i=1;i<=16;i++)
fBR->SetBinContent(i,0.021);
1346 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",16,0,16);
1347 for(Int_t i=1;i<=36;i++)
fTrackingEff->SetBinContent(i,0.15);
1350 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",16,0,16);
1351 for(Int_t i=1;i<=12;i++)
fRawYield->SetBinContent(i,0.09);
1352 for(Int_t i=13; i<=16; i++)
fRawYield->SetBinContent(i,0.12);
1356 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",16,0,16);
1357 for(Int_t i=1; i<=12; i++)
fCutsEff->SetBinContent(i,0.10);
1358 for(Int_t i=12; i<=16; i++)
fCutsEff->SetBinContent(i,0.14);
1361 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",16,0,16);
1362 for(Int_t i=1;i<=16;i++)
fPIDEff->SetBinContent(i,0.05);
1365 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",16,0,16);
1368 for(Int_t iBin=5; iBin<=16; iBin++)
fMCPtShape->SetBinContent(iBin,0.02);
1390 AliInfo(
" Settings for Ds --> K K pi, Pb-Pb collisions for 0-7.5 CC");
1393 fNorm =
new TH1F(
"fNorm",
"fNorm",12,0,12);
1394 for(Int_t i=1;i<=12;i++)
fNorm->SetBinContent(i,0.048);
1397 fBR =
new TH1F(
"fBR",
"fBR",12,0,12);
1398 for(Int_t i=1;i<=12;i++)
fBR->SetBinContent(i,0.053);
1401 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",12,0,12);
1402 for(Int_t i=1;i<=12;i++)
fTrackingEff->SetBinContent(i,0.15);
1405 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",12,0,12);
1406 for(Int_t i=1;i<=6;i++)
fRawYield->SetBinContent(i,.30);
1407 for(Int_t i=7; i<=12; i++)
fRawYield->SetBinContent(i,0.20);
1412 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",12,0,12);
1413 for(Int_t i=1;i<=12;i++)
fCutsEff->SetBinContent(i,0.20);
1416 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",12,0,12);
1417 for(Int_t i=1;i<=12;i++)
fPIDEff->SetBinContent(i,0.1);
1420 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",12,0,12);
1421 for(Int_t i=1; i<=2; i++)
fMCPtShape->SetBinContent(i,1.);
1422 for(Int_t i=3; i<=4; i++)
fMCPtShape->SetBinContent(i,0.03);
1423 for(Int_t i=5; i<=6; i++)
fMCPtShape->SetBinContent(i,0.03);
1424 for(Int_t i=7; i<=8; i++)
fMCPtShape->SetBinContent(i,0.02);
1425 for(Int_t i=9; i<=12; i++)
fMCPtShape->SetBinContent(i,0.02);
1446 AliInfo(
" Settings for Ds --> K K pi, Pb-Pb collisions for 0-10 CC");
1448 fNorm =
new TH1F(
"fNorm",
"fNorm",12,0,12);
1449 for(Int_t i=1;i<=12;i++)
fNorm->SetBinContent(i,0.048);
1452 fBR =
new TH1F(
"fBR",
"fBR",12,0,12);
1453 for(Int_t i=1;i<=12;i++)
fBR->SetBinContent(i,0.045);
1456 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",12,0,12);
1457 for(Int_t i=1;i<=12;i++)
fTrackingEff->SetBinContent(i,0.15);
1460 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",12,0,12);
1461 for(Int_t i=1;i<=12;i++)
fRawYield->SetBinContent(i,0.08);
1465 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",12,0,12);
1466 for(Int_t i=1;i<=12;i++)
fCutsEff->SetBinContent(i,0.20);
1469 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",12,0,12);
1470 for(Int_t i=1;i<=12;i++)
fPIDEff->SetBinContent(i,0.07);
1473 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",12,0,12);
1474 for(Int_t i=1; i<=6; i++)
fMCPtShape->SetBinContent(i,0.02);
1475 for(Int_t i=7; i<=8; i++)
fMCPtShape->SetBinContent(i,0.01);
1476 for(Int_t i=9; i<=12; i++)
fMCPtShape->SetBinContent(i,0.01);
1488 AliInfo(
" Settings for Ds --> K K pi, Pb-Pb collisions for 20-50 CC");
1490 fNorm =
new TH1F(
"fNorm",
"fNorm",12,0,12);
1491 for(Int_t i=1;i<=12;i++)
fNorm->SetBinContent(i,0.051);
1494 fBR =
new TH1F(
"fBR",
"fBR",12,0,12);
1495 for(Int_t i=1;i<=12;i++)
fBR->SetBinContent(i,0.045);
1498 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",12,0,12);
1499 for(Int_t i=1;i<=12;i++)
fTrackingEff->SetBinContent(i,0.15);
1502 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",12,0,12);
1503 for(Int_t i=1;i<=12;i++)
fRawYield->SetBinContent(i,0.08);
1507 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",12,0,12);
1508 for(Int_t i=1;i<=12;i++)
fCutsEff->SetBinContent(i,0.20);
1511 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",12,0,12);
1512 for(Int_t i=1;i<=12;i++)
fPIDEff->SetBinContent(i,0.07);
1515 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",12,0,12);
1516 for(Int_t i=1; i<=12; i++)
fMCPtShape->SetBinContent(i,0.01);
1525 AliInfo(
" Settings for Ds --> K K pi, p-Pb collisions at 5.023 TeV");
1528 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
1529 for(Int_t i=3;i<=12;i++)
fNorm->SetBinContent(i,0.00);
1532 fBR =
new TH1F(
"fBR",
"fBR",12,0,12);
1533 for(Int_t i=3;i<=12;i++)
fBR->SetBinContent(i,0.053);
1536 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
1537 for(Int_t i=3;i<=12;i++)
fTrackingEff->SetBinContent(i,0.09);
1540 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
1541 for(Int_t i=3; i<=4; i++)
fRawYield->SetBinContent(i,0.10);
1542 for(Int_t i=5; i<=6; i++)
fRawYield->SetBinContent(i,0.08);
1543 for(Int_t i=7; i<=8; i++)
fRawYield->SetBinContent(i,0.05);
1544 for(Int_t i=9; i<=12; i++)
fRawYield->SetBinContent(i,0.08);
1547 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
1548 for(Int_t i=3;i<=12;i++)
fCutsEff->SetBinContent(i,0.15);
1551 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
1552 for(Int_t i=3;i<=4;i++)
fPIDEff->SetBinContent(i,0.10);
1553 for(Int_t i=5;i<=12;i++)
fPIDEff->SetBinContent(i,0.05);
1556 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
1557 for(Int_t i=3;i<=12;i++)
fMCPtShape->SetBinContent(i,0.04);
1579 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
1580 for(Int_t i=1;i<=20;i++)
fNorm->SetBinContent(i,0.05);
1583 fBR =
new TH1F(
"fBR",
"fBR",20,0,20);
1584 for(Int_t i=1;i<=20;i++)
fBR->SetBinContent(i,0.021);
1587 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",20,0,20);
1588 for(Int_t i=1;i<=20;i++)
fTrackingEff->SetBinContent(i,0.15);
1591 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
1592 for(Int_t i=1;i<=20;i++)
fRawYield->SetBinContent(i,.10);
1603 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
1604 for(Int_t i=1;i<=20;i++)
fCutsEff->SetBinContent(i,0.15);
1607 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
1608 for(Int_t i=1;i<=20;i++)
fPIDEff->SetBinContent(i,0.05);
1611 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",20,0,20);
1612 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0.);
1613 for(Int_t iBin=7; iBin<=8; iBin++)
fMCPtShape->SetBinContent(iBin,0.01);
1614 for(Int_t iBin=9; iBin<=12; iBin++)
fMCPtShape->SetBinContent(iBin,0.05);
1615 for(Int_t iBin=13; iBin<=16; iBin++)
fMCPtShape->SetBinContent(iBin,0.05);
1639 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
1640 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.07);
1643 fBR =
new TH1F(
"fBR",
"fBR",20,0,20);
1644 for(Int_t i=1;i<=20;i++)
fBR->SetBinContent(i,0.021);
1647 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",20,0,20);
1648 for(Int_t i=1;i<=20;i++)
fTrackingEff->SetBinContent(i,0.15);
1652 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
1661 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.15);
1662 for(Int_t i=13;i<=20;i++)
fRawYield->SetBinContent(i,1);
1665 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
1666 for(Int_t i=1;i<=20;i++)
fCutsEff->SetBinContent(i,0.10);
1669 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
1670 for(Int_t i=1;i<=20;i++)
fPIDEff->SetBinContent(i,0.05);
1671 fPIDEff->SetBinContent(3,0.13);
1675 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",20,0,20);
1676 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0);
1677 for(Int_t iBin=4; iBin<=8; iBin++)
fMCPtShape->SetBinContent(iBin,0.01);
1678 for(Int_t iBin=9; iBin<=12; iBin++)
fMCPtShape->SetBinContent(iBin,0.03);
1679 for(Int_t iBin=13; iBin<=16; iBin++)
fMCPtShape->SetBinContent(iBin,0.03);
1700 AliInfo(
" Settings for D+ --> K pi pi p-p collisions at 2.76 TeV");
1703 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
1704 for(Int_t i=1;i<=20;i++)
fNorm->SetBinContent(i,0.019);
1707 fBR =
new TH1F(
"fBR",
"fBR",20,0,20);
1708 for(Int_t i=1;i<=20;i++)
fBR->SetBinContent(i,0.021);
1711 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",20,0,20);
1712 for(Int_t i=1;i<=20;i++)
fTrackingEff->SetBinContent(i,0.15);
1715 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
1718 for(Int_t i=3;i<=6;i++)
fRawYield->SetBinContent(i,0.10);
1721 for(Int_t i=9;i<=20;i++)
fRawYield->SetBinContent(i,0.055);
1724 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
1725 for(Int_t i=1;i<=20;i++)
fCutsEff->SetBinContent(i,0.15);
1728 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
1729 for(Int_t i=1;i<=20;i++)
fPIDEff->SetBinContent(i,0.05);
1730 fPIDEff->SetBinContent(3,0.10);
1731 fPIDEff->SetBinContent(4,0.10);
1734 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",20,0,20);
1735 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0.01);
1746 AliInfo(
" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV");
1748 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
1749 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.02);
1751 fBR =
new TH1F(
"fBR",
"fBR",20,0,24);
1752 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.021);
1755 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
1756 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.09);
1759 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
1768 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
1769 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
1770 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
1773 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
1782 for(Int_t i=9;i<=12;i++)
fCutsEff->SetBinContent(i,0.05);
1783 for(Int_t i=13;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
1784 for(Int_t i=17;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
1787 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
1789 fPIDEff->SetBinContent(2,0.1);
1790 for(Int_t i=3;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
1793 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
1794 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0.0);
1817 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
1818 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.035);
1821 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
1822 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.015);
1825 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
1830 for(Int_t i=4;i<=24;i++)
fTrackingEff->SetBinContent(i,0.12);
1834 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
1843 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.04);
1844 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
1845 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.2);
1848 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
1850 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.10);
1853 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
1854 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.04);
1858 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
1859 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0);
1874 AliInfo(
" Settings for D*+ --> D0 pi p-p collisions at 2.76 TeV");
1877 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
1878 for(Int_t i=1;i<=20;i++)
fNorm->SetBinContent(i,0.019);
1881 fBR =
new TH1F(
"fBR",
"fBR",20,0,20);
1882 for(Int_t i=1;i<=20;i++)
fBR->SetBinContent(i,0.015);
1885 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",20,0,20);
1886 for(Int_t i=1;i<=20;i++)
fTrackingEff->SetBinContent(i,0.15);
1889 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
1902 for(Int_t i=9;i<=20;i++)
fRawYield->SetBinContent(i,0.065);
1905 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
1906 for(Int_t i=1;i<=20;i++)
fCutsEff->SetBinContent(i,0.10);
1919 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
1920 for(Int_t i=1;i<=20;i++)
fPIDEff->SetBinContent(i,0.05);
1923 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",20,0,20);
1924 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0.01);
1940 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
1941 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.05);
1944 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
1945 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.013);
1948 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
1949 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.09);
1952 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
1956 for(Int_t i=4;i<=8;i++)
fRawYield->SetBinContent(i,0.03);
1957 for(Int_t i=9;i<=24;i++)
fRawYield->SetBinContent(i,0.02);
1960 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
1962 for(Int_t i=2;i<=3;i++)
fCutsEff->SetBinContent(i,0.10);
1963 for(Int_t i=4;i<=7;i++)
fCutsEff->SetBinContent(i,0.05);
1964 for(Int_t i=8;i<=24;i++)
fCutsEff->SetBinContent(i,0.01);
1967 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
1968 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
1971 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
1972 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0);
1973 for(Int_t i=2;i<=3;i++)
fMCPtShape->SetBinContent(i,0.02);
1974 for(Int_t i=4;i<=5;i++)
fMCPtShape->SetBinContent(i,0.01);
1987 AliInfo(
" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-20 centrality - DUMMY");
1990 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
1991 for(Int_t i=1;i<=20;i++)
fNorm->SetBinContent(i,0.05);
1994 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
1995 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.015);
1998 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);;
1999 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.15);
2003 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
2004 for(Int_t i=1;i<=24;i++)
fRawYield->SetBinContent(i,0.1);
2011 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
2012 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.10);
2018 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
2019 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.05);
2023 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
2024 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0.045);
2045 AliInfo(
" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-7.5 centrality - DUMMY");
2049 fNorm =
new TH1F(
"fNorm",
"fNorm",36,0,36);
2050 for(Int_t i=1;i<=36;i++)
fNorm->SetBinContent(i,0.048);
2053 fBR =
new TH1F(
"fBR",
"fBR",36,0,36);
2054 for(Int_t i=1;i<=36;i++)
fBR->SetBinContent(i,0.015);
2057 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",36,0,36);;
2058 for(Int_t i=1;i<=36;i++)
fTrackingEff->SetBinContent(i,0.15);
2062 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",36,0,36);
2063 for(Int_t i=1;i<=36;i++)
fRawYield->SetBinContent(i,0.05);
2068 for(Int_t i=25;i<=36;i++)
fRawYield->SetBinContent(i,0.15);
2071 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",36,0,36);
2072 for(Int_t i=1;i<=36;i++)
fCutsEff->SetBinContent(i,0.10);
2075 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",36,0,36);
2076 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.05);
2077 fPIDEff->SetBinContent(4,0.09);
2080 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",36,0,36);
2081 for(Int_t i=1;i<=36;i++)
fMCPtShape->SetBinContent(i,0.035);
2104 AliInfo(
" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 20-40 centrality - DUMMY");
2107 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
2108 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.10);
2111 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
2112 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.015);
2115 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);;
2116 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.15);
2120 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
2121 for(Int_t i=1;i<=24;i++)
fRawYield->SetBinContent(i,0.15);
2124 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
2125 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.10);
2128 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
2129 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.04);
2133 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
2134 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0.);
2149 AliInfo(
" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 40-80 centrality - DUMMY");
2152 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
2153 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.07);
2156 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
2157 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.015);
2160 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);;
2161 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.15);
2165 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
2166 for(Int_t i=1;i<=24;i++)
fRawYield->SetBinContent(i,0.2);
2186 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
2187 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.10);
2190 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
2191 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.05);
2195 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
2196 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
2216 AliInfo(
" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-10% centrality ");
2220 fNorm =
new TH1F(
"fNorm",
"fNorm",36,0,36);
2221 for(Int_t i=1;i<=36;i++)
fNorm->SetBinContent(i,0.048);
2224 fBR =
new TH1F(
"fBR",
"fBR",36,0,36);
2225 for(Int_t i=1;i<=36;i++)
fBR->SetBinContent(i,0.015);
2228 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",36,0,36);;
2229 for(Int_t i=1;i<=36;i++)
fTrackingEff->SetBinContent(i,0.15);
2233 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",36,0,36);
2234 for(Int_t i=1;i<=3;i++)
fRawYield->SetBinContent(i, 0.);
2238 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2239 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
2240 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.05);
2241 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.05);
2242 for(Int_t i=25;i<=36;i++)
fRawYield->SetBinContent(i,0.10);
2245 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",36,0,36);
2246 for(Int_t i=1;i<=3;i++)
fCutsEff->SetBinContent(i, 0.);
2248 for(Int_t i=5;i<=36;i++)
fCutsEff->SetBinContent(i,0.10);
2251 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",36,0,36);
2252 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.05);
2255 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",36,0,36);
2256 for(Int_t i=1;i<=3;i++)
fMCPtShape->SetBinContent(i, 0.);
2259 for(Int_t i=5;i<=36;i++)
fMCPtShape->SetBinContent(i,0.01);
2272 AliInfo(
" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 30-50 centrality - DUMMY");
2275 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
2276 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.07);
2279 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
2280 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.015);
2283 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);;
2284 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.15);
2288 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
2289 for(Int_t i=1;i<=24;i++)
fRawYield->SetBinContent(i,0.1);
2309 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
2310 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
2318 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
2319 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.05);
2323 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
2324 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
2339 AliInfo(
" Settings for D* --> D0 pi, p-Pb collisions at 5.023 TeV");
2342 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
2343 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.07);
2346 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
2347 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.015);
2350 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
2351 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.09);
2354 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
2363 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
2364 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.06);
2365 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.06);
2368 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
2371 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
2374 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
2375 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
2378 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
2379 for(Int_t i=1;i<=3;i++)
fMCPtShape->SetBinContent(i,0.03);
2380 for(Int_t i=4;i<=6;i++)
fMCPtShape->SetBinContent(i,0.);
2381 for(Int_t i=7;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
2400 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
2401 for(Int_t i=3;i<=5;i++)
fMCPtShape->SetBinContent(i,0.17);
2402 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2408 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
2409 for(Int_t i=3;i<=5;i++)
fMCPtShape->SetBinContent(i,0.17);
2410 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2416 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
2417 for(Int_t i=3;i<=5;i++)
fMCPtShape->SetBinContent(i,0.14);
2418 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2424 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
2425 for(Int_t i=3;i<=5;i++)
fMCPtShape->SetBinContent(i,0.11);
2426 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2432 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.08);
2433 for(Int_t i=3;i<=5;i++)
fMCPtShape->SetBinContent(i,0.08);
2434 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2440 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.18);
2441 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.09);
2448 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.23);
2449 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2455 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.08);
2456 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.095);
2462 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.08);
2463 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2469 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.15);
2470 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.07);
2477 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.16);
2478 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.15);
2484 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
2485 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.15);
2491 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.10);
2492 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2498 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.10);
2499 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.045);
2505 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.10);
2506 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.045);
2551 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.15);
2552 for(Int_t i=5;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2553 for(Int_t i=9;i<=17;i++)
fRawYield->SetBinContent(i,0.15);
2555 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2556 for(Int_t i=5;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2557 for(Int_t i=9;i<=17;i++)
fMCPtShape->SetBinContent(i,0.05);
2559 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.15);
2560 for(Int_t i=5;i<=8;i++)
fCutsEff->SetBinContent(i,0.20);
2561 for(Int_t i=9;i<=17;i++)
fCutsEff->SetBinContent(i,0.10);
2569 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.10);
2570 for(Int_t i=5;i<=8;i++)
fRawYield->SetBinContent(i,0.10);
2571 for(Int_t i=9;i<=17;i++)
fRawYield->SetBinContent(i,0.15);
2573 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2574 for(Int_t i=5;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2575 for(Int_t i=9;i<=17;i++)
fMCPtShape->SetBinContent(i,0.05);
2577 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.15);
2578 for(Int_t i=5;i<=8;i++)
fCutsEff->SetBinContent(i,0.20);
2579 for(Int_t i=9;i<=17;i++)
fCutsEff->SetBinContent(i,0.15);
2587 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.10);
2588 for(Int_t i=5;i<=8;i++)
fRawYield->SetBinContent(i,0.10);
2589 for(Int_t i=9;i<=17;i++)
fRawYield->SetBinContent(i,0.08);
2591 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2592 for(Int_t i=5;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2593 for(Int_t i=9;i<=17;i++)
fMCPtShape->SetBinContent(i,0.05);
2595 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.20);
2596 for(Int_t i=5;i<=8;i++)
fCutsEff->SetBinContent(i,0.15);
2597 for(Int_t i=9;i<=17;i++)
fCutsEff->SetBinContent(i,0.20);
2605 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.20);
2606 for(Int_t i=5;i<=8;i++)
fRawYield->SetBinContent(i,0.10);
2607 for(Int_t i=9;i<=17;i++)
fRawYield->SetBinContent(i,0.20);
2609 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2610 for(Int_t i=5;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2611 for(Int_t i=9;i<=17;i++)
fMCPtShape->SetBinContent(i,0.05);
2613 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.20);
2614 for(Int_t i=5;i<=8;i++)
fCutsEff->SetBinContent(i,0.20);
2615 for(Int_t i=9;i<=17;i++)
fCutsEff->SetBinContent(i,0.20);
2623 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.10);
2624 for(Int_t i=5;i<=8;i++)
fRawYield->SetBinContent(i,0.10);
2625 for(Int_t i=9;i<=17;i++)
fRawYield->SetBinContent(i,0.05);
2627 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2628 for(Int_t i=5;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2629 for(Int_t i=9;i<=17;i++)
fMCPtShape->SetBinContent(i,0.05);
2631 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.20);
2632 for(Int_t i=5;i<=8;i++)
fCutsEff->SetBinContent(i,0.15);
2633 for(Int_t i=9;i<=17;i++)
fCutsEff->SetBinContent(i,0.20);
2641 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.05);
2642 for(Int_t i=5;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2643 for(Int_t i=9;i<=17;i++)
fRawYield->SetBinContent(i,0.10);
2645 for(Int_t i=1;i<=17;i++)
fMCPtShape->SetBinContent(i,0.05);
2646 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.20);
2647 for(Int_t i=5;i<=8;i++)
fCutsEff->SetBinContent(i,0.10);
2648 for(Int_t i=9;i<=17;i++)
fCutsEff->SetBinContent(i,0.20);
2692 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
2694 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.05);
2695 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
2696 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.058);
2698 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2700 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2701 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2702 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2704 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2705 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.028);
2706 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.064);
2714 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
2716 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.042);
2717 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
2718 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.058);
2720 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2722 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2723 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2724 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2726 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2727 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.028);
2728 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.064);
2734 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
2736 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.104);
2737 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.104);
2738 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.104);
2740 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2742 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2743 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2744 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2746 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2747 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.028);
2748 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.064);
2755 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
2757 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.042);
2758 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
2759 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.058);
2761 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2763 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2764 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2765 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2767 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2768 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.028);
2769 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.064);
2776 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
2778 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.05);
2779 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
2780 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.058);
2782 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2784 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2785 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2786 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2788 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2789 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.028);
2790 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.064);
2798 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
2800 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.07);
2801 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2802 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.10);
2804 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2806 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.07);
2807 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.08);
2808 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.08);
2810 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.117);
2811 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.054);
2812 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.064);
2819 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
2821 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.05);
2822 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2823 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.10);
2825 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2827 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.07);
2828 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.08);
2829 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.08);
2831 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.117);
2832 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.054);
2833 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.064);
2839 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
2841 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.05);
2842 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2843 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.05);
2845 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2847 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.07);
2848 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.08);
2849 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.08);
2851 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.117);
2852 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.054);
2853 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.064);
2859 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
2861 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.05);
2862 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2863 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.05);
2865 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2867 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.07);
2868 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.08);
2869 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.08);
2871 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.117);
2872 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.054);
2873 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.064);
2879 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
2881 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.07);
2882 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2883 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.07);
2885 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2887 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.07);
2888 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.08);
2889 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.08);
2891 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.117);
2892 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.054);
2893 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.064);
2902 AliInfo(
"Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
2904 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.06);
2905 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.03);
2906 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.07);
2908 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2910 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2911 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2912 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2914 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.117);
2915 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.054);
2916 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.064);
2922 AliInfo(
"Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
2924 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.04);
2925 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2926 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.05);
2928 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2930 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2931 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2932 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2934 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.117);
2935 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.054);
2936 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.064);
2943 AliInfo(
"Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
2945 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.05);
2946 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.03);
2947 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.05);
2949 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2951 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2952 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2953 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2955 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.117);
2956 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.054);
2957 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.064);
2964 AliInfo(
"Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
2966 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.04);
2967 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.03);
2968 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.05);
2970 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2972 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2973 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2974 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2976 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.117);
2977 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.054);
2978 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.064);
2986 AliInfo(
"Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
2988 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.06);
2989 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.07);
2990 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.07);
2992 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2994 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2995 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2996 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2998 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.117);
2999 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.054);
3000 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.064);
3007 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
3013 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3014 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3015 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3016 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3017 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3018 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3023 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3026 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3032 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3042 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
3049 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3050 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3051 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3052 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3053 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3054 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3059 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3063 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3069 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3077 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
3083 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3084 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3085 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3086 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3087 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3088 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3093 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3097 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3103 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3112 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
3117 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3118 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3119 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3120 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3121 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3122 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3128 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3132 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
3140 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.02);
3150 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult ZNA");
3155 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3156 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3157 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3158 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3159 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3160 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3165 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3169 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3178 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3187 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult ZNA");
3192 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3193 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3194 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3195 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3196 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3197 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3203 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3206 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3215 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3223 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult ZNA");
3228 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3229 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3230 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3231 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3232 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3233 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3238 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3241 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3250 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3260 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult ZNA");
3265 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3266 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3267 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3268 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3269 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3270 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3275 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3278 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
3286 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3295 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult CL1");
3301 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3302 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3303 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3304 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3305 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3306 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3311 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3314 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3320 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3329 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult CL1");
3334 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3335 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3336 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3337 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3338 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3339 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3345 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3349 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3355 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3362 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult CL1");
3368 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3369 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3370 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3371 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3372 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3373 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3378 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3382 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3388 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);}
3394 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult CL1");
3399 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3400 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3401 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3402 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3403 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3404 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3410 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3414 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
3422 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.02);
3430 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
3441 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3442 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3443 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
3448 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3451 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3458 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3468 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
3479 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3480 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3481 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
3487 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3491 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3497 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3514 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
3525 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3526 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3527 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
3532 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3537 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3543 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3561 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
3572 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3573 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3574 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
3580 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3584 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
3593 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.02);
3610 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
3621 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3622 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3623 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
3628 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3632 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.);
3641 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3660 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
3671 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3672 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3673 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
3679 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3683 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
3692 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3710 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
3721 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3722 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3723 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
3730 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3734 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
3743 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3764 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
3775 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3776 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3777 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
3782 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3786 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
3794 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3812 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult CL1");
3823 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3824 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3825 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
3830 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3833 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3839 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3849 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult CL1");
3860 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3861 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3862 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
3867 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3871 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3878 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3895 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 40-60% mult CL1");
3906 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3907 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3908 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.0);
3913 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3917 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3925 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3943 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult CL1");
3954 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3955 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.00);
3956 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.0);
3961 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3965 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
3973 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.02);
3990 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
4001 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4002 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4003 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
4008 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
4015 for(Int_t i=4;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
4018 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
4027 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
4039 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4040 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4041 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
4047 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
4051 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
4058 for(Int_t i=4;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
4075 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
4088 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4089 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4090 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
4098 for(Int_t i=3;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
4102 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
4109 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
4127 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
4139 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4140 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4141 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
4148 for(Int_t i=3;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
4152 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
4161 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.02);
4179 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 0-20% mult ZNA");
4191 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4192 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4193 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
4199 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
4203 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.);
4211 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
4228 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 20-40% mult ZNA");
4240 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4241 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4242 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
4248 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
4253 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
4261 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
4280 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 40-60% mult ZNA");
4292 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4293 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4294 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
4301 for(Int_t i=3;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
4306 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
4314 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
4333 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 60-100% mult ZNA");
4345 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4346 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4347 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
4353 for(Int_t i=3;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
4357 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
4365 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
4383 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 0-20% mult CL1");
4394 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4395 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4396 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.0);
4401 for(Int_t i=3;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
4407 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
4410 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
4419 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 20-40% mult CL1");
4431 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4432 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4433 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.0);
4440 for(Int_t i=3;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
4444 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
4450 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
4467 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 40-60% mult CL1");
4479 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4480 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4481 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.0);
4487 for(Int_t i=3;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
4491 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
4497 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
4516 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 60-100% mult CL1");
4528 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4529 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4530 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
4537 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
4541 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
4550 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.02);
4572 fNorm =
new TH1F(
"fNorm",
"fNorm",4,2,6);
4573 for(Int_t i=1;i<=4;i++)
fNorm->SetBinContent(i,0.035);
4576 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",4,2,6);
4577 for(Int_t i=1;i<=4;i++)
fTrackingEff->SetBinContent(i,0.12);
4580 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
4586 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",4,2,6);
4593 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",4,2,4);
4594 fPIDEff->SetBinContent(1,0.13);
4595 fPIDEff->SetBinContent(2,0.30);
4596 fPIDEff->SetBinContent(3,0.5);
4597 fPIDEff->SetBinContent(4,0.30);
4600 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",12,0,12);
4601 for(Int_t i=1; i<=2; i++)
fMCPtShape->SetBinContent(i,1.);
4602 for(Int_t i=3; i<=4; i++)
fMCPtShape->SetBinContent(i,0.03);
4603 for(Int_t i=5; i<=6; i++)
fMCPtShape->SetBinContent(i,0.03);
4604 for(Int_t i=7; i<=8; i++)
fMCPtShape->SetBinContent(i,0.02);
4605 for(Int_t i=9; i<=12; i++)
fMCPtShape->SetBinContent(i,0.02);
4616 fNorm =
new TH1F(
"fNorm",
"fNorm",10,2.,12.);
4617 for(Int_t i=1;i<=10;i++)
fNorm->SetBinContent(i,0.05);
4620 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",10,2.,12.);
4621 for(Int_t i=1;i<=10;i++)
fTrackingEff->SetBinContent(i,0.12);
4624 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",10,2.,12.);
4625 for(Int_t i=1;i<=10;i++)
fRawYield->SetBinContent(i,0.11);
4627 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",10,2.,12.);
4628 for(Int_t i=1;i<=2;i++)
fCutsEff->SetBinContent(i,0.09);
4629 for(Int_t i=3;i<=10;i++)
fCutsEff->SetBinContent(i,0.08);
4632 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",10,2.,12.);
4633 for(Int_t i=1;i<=2;i++)
fPIDEff->SetBinContent(i,0.04);
4634 for(Int_t i=3;i<=12;i++)
fPIDEff->SetBinContent(i,0.05);
4637 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",10,2.,12.);
4638 for(Int_t i=1; i<=10;i++)
fMCPtShape->SetBinContent(i,0.005);
4641 fBR =
new TH1F(
"fBR",
"fBR",10,2,12);
4642 for(Int_t i=1;i<=10;i++)
fBR->SetBinContent(i,0.093);
4654 fNorm =
new TH1F(
"fNorm",
"fNorm",10,2.,12.);
4655 for(Int_t i=1;i<=10;i++)
fNorm->SetBinContent(i,0.05);
4658 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",10,2.,12.);
4659 for(Int_t i=1;i<=10;i++)
fTrackingEff->SetBinContent(i,0.09);
4662 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",10,2.,12.);
4663 for(Int_t i=1;i<=2;i++)
fRawYield->SetBinContent(i,0.11);
4664 for(Int_t i=3;i<=10;i++)
fRawYield->SetBinContent(i,0.11);
4666 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",10,2.,12.);
4667 for(Int_t i=1;i<=2;i++)
fCutsEff->SetBinContent(i,0.08);
4668 for(Int_t i=3;i<=10;i++)
fCutsEff->SetBinContent(i,0.09);
4671 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",10,2.,12.);
4672 for(Int_t i=1;i<=2;i++)
fPIDEff->SetBinContent(i,0.04);
4673 for(Int_t i=3;i<=10;i++)
fPIDEff->SetBinContent(i,0.05);
4676 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",10,2.,12.);
4677 for(Int_t i=1; i<=2;i++)
fMCPtShape->SetBinContent(i,0.005);
4678 for(Int_t i=3; i<=10;i++)
fMCPtShape->SetBinContent(i,0.005);
4681 fBR =
new TH1F(
"fBR",
"fBR",10,2,12);
4682 for(Int_t i=1;i<=10;i++)
fBR->SetBinContent(i,0.093);
4694 return fCutsEff->GetBinContent(bin);
4714 return TMath::Sqrt(err);
4722 Int_t bin=
fPIDEff->FindBin(pt);
4724 return fPIDEff->GetBinContent(bin);
4772 err += feeddownErr*feeddownErr;
4774 return TMath::Sqrt(err);
4781 gStyle->SetOptStat(0);
4783 TCanvas *cSystErr =
new TCanvas(
"cSystErr",
"Systematic Errors",300,80,640,500);
4784 cSystErr->Range(0.20,-0.5,18.4,0.34);
4785 cSystErr->SetRightMargin(0.318);
4786 cSystErr->SetFillColor(0);
4788 TH2F *hFrame =
new TH2F(
"hFrame",
"Systematic errors; p_{t} (GeV/c); Relative Error",40,0,40,100,-1,+1);
4789 hFrame->SetAxisRange(1.,35.9,
"X");
4791 hFrame->SetAxisRange(-0.5,0.5,
"Y");
4794 TLegend *leg =
new TLegend(0.69,0.44,0.98,0.86,NULL,
"brNDC");
4795 leg->SetTextSize(0.03601695);
4796 leg->SetFillStyle(0);
4797 leg->SetBorderSize(0);
4799 TH1F *hTotErr=
new TH1F(
"hTotErr",
"",36,0,36);
4801 TGraphAsymmErrors *gTotErr =
new TGraphAsymmErrors(nbins);
4802 for(Int_t i=1;i<=36;i++) {
4803 Double_t pt = hTotErr->GetBinCenter(i);
4804 Double_t ptwidth = hTotErr->GetBinWidth(i);
4807 Double_t x=0., y=0., errxl=0., errxh=0., erryl=0., erryh=0.;
4808 Double_t toterryl=0., toterryh=0.;
4809 for(Int_t j=0; j<grErrFeeddown->GetN(); j++) {
4810 grErrFeeddown->GetPoint(j,x,y);
4811 errxh = grErrFeeddown->GetErrorXhigh(j);
4812 errxl = grErrFeeddown->GetErrorXlow(j);
4813 if ( ( (x-errxl) <= pt) && ( (x+errxl) >= pt) ) {
4814 erryh = grErrFeeddown->GetErrorYhigh(j);
4815 erryl = grErrFeeddown->GetErrorYlow(j);
4823 hTotErr->SetBinContent(i,toterryh);
4824 gTotErr->SetPoint(i,pt,0.);
4825 gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,toterryl,toterryh);
4829 gTotErr->SetPoint(i,pt,0.);
4834 gTotErr->SetLineColor(kBlack);
4835 gTotErr->SetFillColor(kRed);
4836 gTotErr->SetFillStyle(3002);
4838 leg->AddEntry(gTotErr,
"Total (excl. norm.)",
"f");
4845 fNorm->SetFillColor(1);
4846 fNorm->SetFillStyle(3002);
4850 Double_t norm =
fNorm->GetBinContent(1)*100;
4851 leg->AddEntry(
fNorm,Form(
"Normalization (%.1f%s)",norm,
"%"),
"");
4854 grErrFeeddown->SetFillColor(kTeal-8);
4855 grErrFeeddown->SetFillStyle(3001);
4856 grErrFeeddown->Draw(
"2");
4857 leg->AddEntry(grErrFeeddown,
"Feed-down from B",
"f");
4864 hTrackingEffRefl->Draw(
"same");
4868 fBR->SetFillColor(6);
4869 fBR->SetFillStyle(3005);
4873 hBRRefl->Draw(
"same");
4874 leg->AddEntry(
fBR,
"Branching ratio",
"f");
4878 ci = TColor::GetColor(
"#00cc00");
4884 hRawYieldRefl->Draw(
"same");
4885 leg->AddEntry(
fRawYield,
"Yield extraction",
"l");
4892 hCutsEffRefl->Draw(
"same");
4893 leg->AddEntry(
fCutsEff,
"Cuts efficiency",
"l");
4900 hPIDEffRefl->Draw(
"same");
4901 leg->AddEntry(
fPIDEff,
"PID efficiency",
"l");
4904 Int_t ci = TColor::GetColor(
"#9933ff");
4910 hMCPtShapeRefl->Draw(
"same");
4911 leg->AddEntry(
fMCPtShape,
"MC p_{t} shape",
"l");
4914 Int_t ci = TColor::GetColor(
"#ff6600");
4920 hPartAntipartRefl->Draw(
"same");
4927 cSystErr->SaveAs(
"RelativeSystematics.eps");
4936 TH1F *hout=(TH1F*)hin->Clone(
"hout");
void InitDstartoD0pi2010PbPb2040CentScan()
void InitDstartoD0pi2011PbPb2030CentScan()
void InitDplustoKpipi2011PbPb3040CentScan()
void InitDstartoD0pi2013pPb2040ZNA()
void InitD0toKpi2010PbPb6080CentScan()
void InitDplustoKpipi2010ppLowEn()
void InitDplustoKpipi2013pPb020ZNA()
void InitD0toKpi2013pPb0100()
void InitD0toKpi2013pPb60100CL1()
void InitDplustoKpipi2013pPb60100CL1()
void InitDstartoD0pi2013pPb2040V0A()
TString fRapidityRange
MB:0100, 0-10:010, 0-20:020 ...40-80:4080...
void InitDstartoD0pi2010PbPb6080CentScan()
void InitDstartoD0pi2013pPb0100RapScan0104()
void InitDplustoKpipi2011PbPb3050()
void InitDplustoKpipi2013pPb0100RapScan0408()
void InitD0toKpi2013pPb2040CL1()
void InitD0toKpi2013pPb4060CL1()
TH1F * fPIDEff
cuts efficiency
void InitDplustoKpipi2011PbPb010()
void InitDstartoD0pi2013pPb60100CL1()
void InitDplustoKpipi2011PbPb010CentScan()
void InitD0toKpi2011PbPb1020CentScan()
Bool_t fIsLowEnergy
Rapidity range fot y measurements.
void InitDplustoKpipi2013pPb0100RapScan0401()
void InitDstartoD0pi2011PbPb3040CentScan()
void InitD0toKpi2013pPb0100RapScan0408()
TH1F * fBR
tracking efficiency
void InitD0toKpi2010PbPb020()
void InitDstartoD0pi2013pPb0100()
void InitDstartoD0pi2012pp()
TH1F * fMCPtShape
PID efficiency.
void InitDplustoKpipi2013pPb4060V0A()
void InitLctopK0S2010pp()
void InitDstartoD0pi2010PbPb2040()
void InitDstartoD0pi2010PbPb4060CentScan()
void InitDplustoKpipi2010PbPb1020CentScan()
void InitDstoKKpi2011PbPb2050()
Double_t GetCutsEffErr(Double_t pt) const
void InitDstartoD0pi2011PbPb010()
void InitDstartoD0pi2010PbPb020()
void InitDstartoD0pi2013pPb60100ZNA()
void InitDplustoKpipi2013pPb0100RapScan0101()
void InitDplustoKpipi2010PbPb020()
void InitDstartoD0pi2013pPb0100RapScan0408()
TH1F * ReflectHisto(TH1F *hin) const
void InitD0toKpi2011PbPb010CentScan()
void InitDplustoKpipi2013pPb60100ZNA()
void InitDplustoKpipi2013pPb0100RapScan0804()
void InitD0toKpi2011PbPb3050()
void InitDstartoD0pi2011PbPb1020CentScan()
Double_t GetPartAntipartErr(Double_t pt) const
void InitDstartoD0pi2011PbPb4050CentScan()
void InitDplustoKpipi2010PbPb4060CentScan()
void InitD0toKpi2013pPb2040V0A()
void InitDplustoKpipi2010PbPb4080()
void InitD0toKpi2013pPb4060ZNA()
void InitD0toKpi2010PbPb2040CentScan()
Double_t GetSeleEffErr(Double_t pt) const
void InitD0toKpi2013pPb0100RapScan0401()
void InitD0toKpi2013pPb020V0A()
void InitD0toKpi2010PbPb010CentScan()
void InitD0toKpi2013pPb0100RapScan0804()
void InitD0toKpi2013pPb2040ZNA()
void InitDstartoD0pi2010ppLowEn()
void InitLctopKpi2010pp()
void InitDstartoD0pi2013pPb020V0A()
void InitLctopK0S2013pPb()
void InitD0toKpi2010PbPb5080CentScan()
void InitD0toKpi2013pPb0100RapScan0104()
void InitDstartoD0pi2010PbPb1020CentScan()
void InitD0toKpi2013pPb020CL1()
void InitDstartoD0pi2010PbPb010CentScan()
Double_t GetTrackingEffErr(Double_t pt) const
void InitD0toKpi2013pPb60100V0A()
void InitDstartoD0pi2011PbPb07half()
void InitD0toKpi2011PbPb07half()
void InitDstartoD0pi2013pPb4060ZNA()
void InitDstartoD0pi2013pPb0100RapScan0401()
TH1F * fRawYield
normalization
void InitDstartoD0pi2010PbPb5080CentScan()
void InitDplustoKpipi2013pPb2040CL1()
void InitDstartoD0pi2010PbPb4080()
void InitDplustoKpipi2013pPb020V0A()
void InitDstartoD0pi2013pPb4060V0A()
AliHFSystErr(const Char_t *name="HFSystErr", const Char_t *title="")
void InitDstartoD0pi2013pPb0100RapScan0101()
void InitDplustoKpipi2010PbPb6080CentScan()
TH1F * fPartAntipart
MC dNdpt.
void InitDstartoD0pi2013pPb020ZNA()
Bool_t fIsRapidityScan
flag fot the PbPb centrality scan
void InitD0toKpi2013pPb4060V0A()
void InitDstoKKpi2013pPb0100()
void InitDplustoKpipi2013pPb4060ZNA()
TH1F * fTrackingEff
raw yield
void InitDplustoKpipi2013pPb0100()
TString fCentralityClass
Collision type: pp=0, PbPb=1.
void InitDplustoKpipi2010pp()
Double_t GetTotalSystErr(Double_t pt, Double_t feeddownErr=0) const
void InitD0toKpi2011PbPb3040CentScan()
void InitDplustoKpipi2010PbPb010CentScan()
void InitDplustoKpipi2011PbPb4050CentScan()
void InitDstartoD0pi2013pPb020CL1()
void InitD0toKpi2013pPb60100ZNA()
void InitD0toKpi2013pPb020ZNA()
void InitD0toKpi2010ppLowEn()
Bool_t fIsLowPtAnalysis
flag for the low energy (2.76TeV) run
void InitDstoKKpi2011PbPb07half()
void InitD0toKpi2010PbPb4060CentScan()
void InitDplustoKpipi2010ppPass4()
void Init(Int_t decay)
Function to initialize the variables/histograms.
void InitD0toKpi2010PbPb1020CentScan()
void InitDstartoD0pi2011PbPb010CentScan()
Double_t GetRawYieldErr(Double_t pt) const
ClassImp(AliAnalysisTaskCRC) AliAnalysisTaskCRC
void InitDstoKKpi2010pp()
void InitDplustoKpipi2011PbPb1020CentScan()
void InitDplustoKpipi2013pPb020CL1()
void InitD0toKpi2010ppLowPtAn()
void InitD0toKpi2010PbPb4080()
void InitDstartoD0pi2013pPb2040CL1()
void InitD0toKpi2011PbPb2030CentScan()
void InitDplustoKpipi2011PbPb07half()
void InitDstartoD0pi2010pp()
void InitDplustoKpipi2013pPb2040ZNA()
void InitD0toKpi2011PbPb010()
void InitDplustoKpipi2013pPb2040V0A()
void InitD0toKpi2011PbPb4050CentScan()
void InitDstartoD0pi2013pPb4060CL1()
void InitDplustoKpipi2013pPb4060CL1()
void InitDplustoKpipi2010PbPb2040CentScan()
void InitD0toKpi2013pPb0100RapScan0101()
Bool_t fIsPass4Analysis
flag for the low pt analysis (no topological cuts)
void DrawErrors(TGraphAsymmErrors *grErrFeeddown=0) const
TH1F * fCutsEff
branching ratio
void InitDstoKKpi2011PbPb010()
Double_t GetMCPtShapeErr(Double_t pt) const
Int_t fCollisionType
Run Number (year)
void InitDplustoKpipi2013pPb0100RapScan0104()
void InitDstartoD0pi2013pPb60100V0A()
void InitD0toKpi2011PbPb3050OutOfPlane()
void InitDstartoD0pi2013pPb0100RapScan0804()
void InitD0toKpi2013pPb0100LowPtAn()
void InitDplustoKpipi2011PbPb2030CentScan()
Int_t fRunNumber
particle=antiparticle
void InitDstartoD0pi2011PbPb3050()
Bool_t fIsCentScan
flag for the pass4 analysis
void InitDplustoKpipi2012pp()
void InitDplustoKpipi2013pPb60100V0A()
void InitD0toKpi2011PbPb3050InPlane()
void InitDplustoKpipi2010PbPb5080CentScan()