36 #include <TGraphAsymmErrors.h>
63 fCentralityClass(
"0100"),
64 fRapidityRange(
"0101"),
66 fIsLowPtAnalysis(false),
68 fIsRapidityScan(false)
100 AliFatal(
"Only settings for 2010 and the low energy runs are implemented so far");
119 else AliFatal(
"Not yet implemented");
133 else AliFatal(
"Not yet implemented");
164 else AliFatal(
"Not yet implemented");
183 else AliFatal(
"Not yet implemented");
195 else AliFatal(
"Not yet implemented");
226 else AliFatal(
"Not yet implemented");
244 else AliFatal(
"Not yet implemented");
256 else AliFatal(
"Not yet implemented");
285 else AliFatal(
"Not yet implemented");
294 else AliFatal(
"Not yet implemented");
299 else AliFatal(
"Not yet implemented");
303 else AliFatal(
"Not yet implemented");
307 else AliFatal(
"Not yet implemented");
311 printf(
"Invalid decay type: %d\n",decay);
325 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
326 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.035);
329 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
330 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.012);
333 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
334 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.08);
337 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
341 for(Int_t i=4;i<=7;i++)
fRawYield->SetBinContent(i,0.04);
342 for(Int_t i=8;i<=12;i++)
fRawYield->SetBinContent(i,0.07);
343 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.10);
344 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,1);
347 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
348 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.10);
351 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
352 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.03);
353 fPIDEff->SetBinContent(2,0.05);
356 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
357 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0);
376 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
377 for(Int_t i=1;i<=20;i++)
fNorm->SetBinContent(i,0.05);
380 fBR =
new TH1F(
"fBR",
"fBR",20,0,20);
381 for(Int_t i=1;i<=20;i++)
fBR->SetBinContent(i,0.012);
384 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",20,0,20);
385 for(Int_t i=1;i<=20;i++)
fTrackingEff->SetBinContent(i,0.10);
388 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
392 for(Int_t i=4;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
393 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.10);
396 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
401 for(Int_t i=5;i<=16;i++)
fCutsEff->SetBinContent(i,0.10);
402 for(Int_t i=17;i<=20;i++)
fCutsEff->SetBinContent(i,0.);
405 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
406 for(Int_t i=3;i<=16;i++)
fPIDEff->SetBinContent(i,0.05);
409 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",20,0,20);
410 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0.01);
413 for(Int_t i=13;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
434 fNorm =
new TH1F(
"fNorm",
"fNorm",36,0,36);
435 for(Int_t i=1;i<36;i++)
fNorm->SetBinContent(i,0.048);
438 fBR =
new TH1F(
"fBR",
"fBR",36,0,36);
439 for(Int_t i=1;i<=36;i++)
fBR->SetBinContent(i,0.012);
442 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",36,0,36);
443 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.10);
444 for(Int_t i=25;i<=36;i++)
fTrackingEff->SetBinContent(i,0.);
447 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",36,0,36);
451 for(Int_t i=4;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
452 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.10);
453 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.30);
454 for(Int_t i=25;i<=36;i++)
fRawYield->SetBinContent(i,0.);
457 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",36,0,36);
463 for(Int_t i=6;i<=24;i++)
fCutsEff->SetBinContent(i,0.06);
464 for(Int_t i=25;i<=36;i++)
fCutsEff->SetBinContent(i,0.0);
467 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",36,0,36);
468 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.05);
469 for(Int_t i=25;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
472 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",36,0,36);
473 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
477 for(Int_t i=13;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
478 for(Int_t i=17;i<=24;i++)
fMCPtShape->SetBinContent(i,0.05);
479 for(Int_t i=25;i<=36;i++)
fMCPtShape->SetBinContent(i,0.0);
505 for(Int_t i=5;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
506 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
507 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.15);
508 for(Int_t i=17;i<=36;i++)
fRawYield->SetBinContent(i,0.);
518 for(Int_t i=7;i<=8;i++)
fCutsEff->SetBinContent(i,0.15);
519 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.15);
520 for(Int_t i=25;i<=36;i++)
fCutsEff->SetBinContent(i,0.0);
535 for(Int_t i=4;i<=6;i++)
fRawYield->SetBinContent(i,0.07);
536 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
537 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.10);
538 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.15);
539 for(Int_t i=17;i<=36;i++)
fRawYield->SetBinContent(i,0.);
549 for(Int_t i=7;i<=8;i++)
fCutsEff->SetBinContent(i,0.15);
550 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.15);
551 for(Int_t i=17;i<=36;i++)
fCutsEff->SetBinContent(i,0.0);
562 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
563 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.07);
566 fBR =
new TH1F(
"fBR",
"fBR",20,0,20);
567 for(Int_t i=1;i<=20;i++)
fBR->SetBinContent(i,0.012);
570 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",20,0,20);
571 for(Int_t i=1;i<=20;i++)
fTrackingEff->SetBinContent(i,0.10);
575 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
578 for(Int_t i=3;i<=16;i++)
fRawYield->SetBinContent(i,0.05);
582 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
587 for(Int_t i=5;i<=16;i++)
fCutsEff->SetBinContent(i,0.10);
588 for(Int_t i=17;i<=20;i++)
fCutsEff->SetBinContent(i,0.);
591 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
594 for(Int_t i=3;i<=16;i++)
fPIDEff->SetBinContent(i,0.05);
597 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",20,0,20);
598 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0.01);
601 for(Int_t i=13;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
617 AliInfo(
" Settings for D0 --> Kpi Pb-Pb collisions at 2.76 TeV - 0-10 centrality - newww!!");
620 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
621 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.048);
624 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
625 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.012);
628 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
629 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.10);
633 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
640 for (Int_t i=7; i<=8; i++)
fRawYield->SetBinContent(i,0.06);
641 for (Int_t i=9; i<=12; i++)
fRawYield->SetBinContent(i,0.07);
642 for (Int_t i=13; i<=16; i++)
fRawYield->SetBinContent(i,0.11);
643 for (Int_t i=17; i<=24; i++)
fRawYield->SetBinContent(i,0.30);
645 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
652 for (Int_t i=7; i<=8; i++)
fCutsEff->SetBinContent(i,0.05);
653 for (Int_t i=9; i<=12; i++)
fCutsEff->SetBinContent(i,0.05);
654 for (Int_t i=13; i<=16; i++)
fCutsEff->SetBinContent(i,0.05);
655 for (Int_t i=17; i<=24; i++)
fCutsEff->SetBinContent(i,0.05);
659 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
662 for(Int_t i=2;i<=16;i++)
fPIDEff->SetBinContent(i,0.05);
665 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
670 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.);
687 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
688 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.07);
691 fBR =
new TH1F(
"fBR",
"fBR",20,0,20);
692 for(Int_t i=1;i<=20;i++)
fBR->SetBinContent(i,0.012);
695 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",20,0,20);
696 for(Int_t i=1;i<=20;i++)
fTrackingEff->SetBinContent(i,0.10);
700 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
714 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.112);
715 for(Int_t i=17;i<=20;i++)
fRawYield->SetBinContent(i,0);
717 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
724 for(Int_t i=7;i<=8;i++)
fCutsEff->SetBinContent(i,0.15);
725 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.15);
729 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
732 for(Int_t i=2;i<=16;i++)
fPIDEff->SetBinContent(i,0.05);
735 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",0,0,20);
736 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0.01);
740 for(Int_t i=13;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
741 for(Int_t i=17;i<=20;i++)
fMCPtShape->SetBinContent(i,0.);
756 AliInfo(
" Settings for D0 --> K pi, p-p collisions at 2.76 TeV");
759 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
760 for(Int_t i=1;i<=20;i++)
fNorm->SetBinContent(i,0.019);
763 fBR =
new TH1F(
"fBR",
"fBR",20,0,20);
764 for(Int_t i=1;i<=20;i++)
fBR->SetBinContent(i,0.012);
767 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",20,0,20);
768 for(Int_t i=1;i<=20;i++)
fTrackingEff->SetBinContent(i,0.10);
771 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
773 for(Int_t i=1;i<=20;i++)
fRawYield->SetBinContent(i,0.15);
776 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
777 for(Int_t i=1;i<=20;i++)
fCutsEff->SetBinContent(i,0.10);
779 for(Int_t i=7;i<=20;i++)
fCutsEff->SetBinContent(i,0.15);
783 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
784 for(Int_t i=1;i<=20;i++)
fPIDEff->SetBinContent(i,0.15);
786 for(Int_t i=7;i<=20;i++)
fPIDEff->SetBinContent(i,0.05);
789 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",20,0,20);
790 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0.01);
810 AliInfo(
" Settings for D0 --> K pi, pp collisions at 7 TeV, analysis without topological cuts");
813 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
814 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.035);
817 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
818 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.0129);
821 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
822 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.08);
825 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
826 for(Int_t i=1;i<=24;i++)
fRawYield->SetBinContent(i,0.2);
835 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.2);
838 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
839 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.);
842 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
843 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.03);
844 fPIDEff->SetBinContent(1,0.05);
845 fPIDEff->SetBinContent(2,0.05);
848 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
849 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0);
862 AliInfo(
" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV");
864 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
865 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.02);
868 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
869 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.013);
872 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
873 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.06);
876 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
879 for(Int_t i=3;i<=12;i++)
fRawYield->SetBinContent(i,0.04);
880 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.07);
881 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.11);
884 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
887 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
890 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
891 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
894 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
897 for(Int_t i=3;i<=24;i++)
fMCPtShape->SetBinContent(i,0.0);
915 AliInfo(
" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV, analysis without topological cuts");
918 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
919 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.035);
922 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
923 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.0129);
926 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
927 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.06);
930 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
931 for(Int_t i=1;i<=24;i++)
fRawYield->SetBinContent(i,0.2);
940 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.11);
943 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
944 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.);
947 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
948 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.);
951 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
952 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0);
972 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
973 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.035);
976 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
977 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.021);
980 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
981 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.12);
985 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
992 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
993 for(Int_t i=12;i<=24;i++)
fRawYield->SetBinContent(i,0.10);
996 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
997 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.10);
1000 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
1001 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.05);
1002 fPIDEff->SetBinContent(1,0.15);
1003 fPIDEff->SetBinContent(2,0.15);
1004 fPIDEff->SetBinContent(3,0.15);
1005 fPIDEff->SetBinContent(4,0.15);
1006 for(Int_t i=12;i<=16;i++)
fPIDEff->SetBinContent(i,0.10);
1009 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
1010 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0);
1035 fNorm =
new TH1F(
"fNorm",
"fNorm",12,0,12);
1036 for(Int_t i=1;i<=12;i++)
fNorm->SetBinContent(i,0.035);
1039 fBR =
new TH1F(
"fBR",
"fBR",12,0,12);
1040 for(Int_t i=1;i<=12;i++)
fBR->SetBinContent(i,0.06);
1043 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",12,0,12);
1044 for(Int_t i=1;i<=12;i++)
fTrackingEff->SetBinContent(i,0.12);
1048 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",12,0,12);
1063 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",12,0,12);
1064 for(Int_t i=1;i<=12;i++)
fCutsEff->SetBinContent(i,0.15);
1067 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",12,0,12);
1068 for(Int_t i=1;i<=12;i++)
fPIDEff->SetBinContent(i,0.07);
1071 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",12,0,12);
1072 for(Int_t i=1; i<=2; i++)
fMCPtShape->SetBinContent(i,1.);
1073 for(Int_t i=3; i<=4; i++)
fMCPtShape->SetBinContent(i,0.03);
1074 for(Int_t i=5; i<=6; i++)
fMCPtShape->SetBinContent(i,0.03);
1075 for(Int_t i=7; i<=8; i++)
fMCPtShape->SetBinContent(i,0.02);
1076 for(Int_t i=9; i<=12; i++)
fMCPtShape->SetBinContent(i,0.02);
1099 fNorm =
new TH1F(
"fNorm",
"fNorm",36,0,36);
1100 for(Int_t i=1;i<=36;i++)
fNorm->SetBinContent(i,0.048);
1103 fBR =
new TH1F(
"fBR",
"fBR",36,0,36);
1104 for(Int_t i=1;i<=36;i++)
fBR->SetBinContent(i,0.021);
1107 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",36,0,36);
1108 for(Int_t i=1;i<=36;i++)
fTrackingEff->SetBinContent(i,0.15);
1111 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",36,0,36);
1112 for(Int_t i=1;i<=36;i++)
fRawYield->SetBinContent(i,.10);
1116 for(Int_t i=7; i<=8; i++)
fRawYield->SetBinContent(i,0.10);
1117 for(Int_t i=9; i<=12; i++)
fRawYield->SetBinContent(i,0.08);
1118 for(Int_t i=13; i<=16; i++)
fRawYield->SetBinContent(i,0.05);
1119 for(Int_t i=17; i<=24; i++)
fRawYield->SetBinContent(i,0.08);
1120 for(Int_t i=25; i<=36; i++)
fRawYield->SetBinContent(i,0.20);
1124 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",36,0,36);
1125 for(Int_t i=1; i<=12; i++)
fCutsEff->SetBinContent(i,0.10);
1126 for(Int_t i=13; i<=36; i++)
fCutsEff->SetBinContent(i,0.08);
1130 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",36,0,36);
1131 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.05);
1134 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",36,0,36);
1135 for(Int_t iBin=1; iBin<=5; iBin++)
fMCPtShape->SetBinContent(iBin,0.05);
1136 for(Int_t iBin=6; iBin<=36; iBin++)
fMCPtShape->SetBinContent(iBin,0.03);
1159 fNorm =
new TH1F(
"fNorm",
"fNorm",36,0,36);
1160 for(Int_t i=1;i<=36;i++)
fNorm->SetBinContent(i,0.048);
1163 fBR =
new TH1F(
"fBR",
"fBR",36,0,36);
1164 for(Int_t i=1;i<=36;i++)
fBR->SetBinContent(i,0.021);
1167 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",36,0,36);
1168 for(Int_t i=1;i<=36;i++)
fTrackingEff->SetBinContent(i,0.15);
1171 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",36,0,36);
1172 for(Int_t i=1;i<=3;i++)
fRawYield->SetBinContent(i,0.);
1176 for(Int_t i=7; i<=8; i++)
fRawYield->SetBinContent(i,0.08);
1177 for(Int_t i=9; i<=12; i++)
fRawYield->SetBinContent(i,0.08);
1178 for(Int_t i=13; i<=16; i++)
fRawYield->SetBinContent(i,0.08);
1179 for(Int_t i=17; i<=24; i++)
fRawYield->SetBinContent(i,0.08);
1180 for(Int_t i=25; i<=36; i++)
fRawYield->SetBinContent(i,0.08);
1184 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",36,0,36);
1185 for(Int_t i=1; i<=36; i++)
fCutsEff->SetBinContent(i,0.10);
1189 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",36,0,36);
1190 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.05);
1193 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",36,0,36);
1194 for(Int_t i=1;i<=3;i++)
fMCPtShape->SetBinContent(i,0.);
1197 for (Int_t i=6; i<=36; i++)
fMCPtShape->SetBinContent(i,0.01);
1219 AliInfo(
" Settings for D+ --> K pi pi, Pb-Pb collisions for 30-50 CC");
1221 fNorm =
new TH1F(
"fNorm",
"fNorm",16,0,16);
1222 for(Int_t i=1;i<=16;i++)
fNorm->SetBinContent(i,0.048);
1225 fBR =
new TH1F(
"fBR",
"fBR",16,0,16);
1226 for(Int_t i=1;i<=16;i++)
fBR->SetBinContent(i,0.021);
1229 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",16,0,16);
1230 for(Int_t i=1;i<=36;i++)
fTrackingEff->SetBinContent(i,0.15);
1233 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",16,0,16);
1234 for(Int_t i=1;i<=12;i++)
fRawYield->SetBinContent(i,0.09);
1235 for(Int_t i=13; i<=16; i++)
fRawYield->SetBinContent(i,0.12);
1239 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",16,0,16);
1240 for(Int_t i=1; i<=12; i++)
fCutsEff->SetBinContent(i,0.10);
1241 for(Int_t i=12; i<=16; i++)
fCutsEff->SetBinContent(i,0.14);
1244 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",16,0,16);
1245 for(Int_t i=1;i<=16;i++)
fPIDEff->SetBinContent(i,0.05);
1248 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",16,0,16);
1251 for(Int_t iBin=5; iBin<=16; iBin++)
fMCPtShape->SetBinContent(iBin,0.02);
1273 AliInfo(
" Settings for Ds --> K K pi, Pb-Pb collisions for 0-7.5 CC");
1276 fNorm =
new TH1F(
"fNorm",
"fNorm",12,0,12);
1277 for(Int_t i=1;i<=12;i++)
fNorm->SetBinContent(i,0.048);
1280 fBR =
new TH1F(
"fBR",
"fBR",12,0,12);
1281 for(Int_t i=1;i<=12;i++)
fBR->SetBinContent(i,0.053);
1284 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",12,0,12);
1285 for(Int_t i=1;i<=12;i++)
fTrackingEff->SetBinContent(i,0.15);
1288 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",12,0,12);
1289 for(Int_t i=1;i<=6;i++)
fRawYield->SetBinContent(i,.30);
1290 for(Int_t i=7; i<=12; i++)
fRawYield->SetBinContent(i,0.20);
1295 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",12,0,12);
1296 for(Int_t i=1;i<=12;i++)
fCutsEff->SetBinContent(i,0.20);
1299 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",12,0,12);
1300 for(Int_t i=1;i<=12;i++)
fPIDEff->SetBinContent(i,0.1);
1303 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",12,0,12);
1304 for(Int_t i=1; i<=2; i++)
fMCPtShape->SetBinContent(i,1.);
1305 for(Int_t i=3; i<=4; i++)
fMCPtShape->SetBinContent(i,0.03);
1306 for(Int_t i=5; i<=6; i++)
fMCPtShape->SetBinContent(i,0.03);
1307 for(Int_t i=7; i<=8; i++)
fMCPtShape->SetBinContent(i,0.02);
1308 for(Int_t i=9; i<=12; i++)
fMCPtShape->SetBinContent(i,0.02);
1329 AliInfo(
" Settings for Ds --> K K pi, Pb-Pb collisions for 0-10 CC");
1331 fNorm =
new TH1F(
"fNorm",
"fNorm",12,0,12);
1332 for(Int_t i=1;i<=12;i++)
fNorm->SetBinContent(i,0.048);
1335 fBR =
new TH1F(
"fBR",
"fBR",12,0,12);
1336 for(Int_t i=1;i<=12;i++)
fBR->SetBinContent(i,0.045);
1339 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",12,0,12);
1340 for(Int_t i=1;i<=12;i++)
fTrackingEff->SetBinContent(i,0.15);
1343 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",12,0,12);
1344 for(Int_t i=1;i<=12;i++)
fRawYield->SetBinContent(i,0.08);
1348 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",12,0,12);
1349 for(Int_t i=1;i<=12;i++)
fCutsEff->SetBinContent(i,0.20);
1352 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",12,0,12);
1353 for(Int_t i=1;i<=12;i++)
fPIDEff->SetBinContent(i,0.07);
1356 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",12,0,12);
1357 for(Int_t i=1; i<=6; i++)
fMCPtShape->SetBinContent(i,0.02);
1358 for(Int_t i=7; i<=8; i++)
fMCPtShape->SetBinContent(i,0.01);
1359 for(Int_t i=9; i<=12; i++)
fMCPtShape->SetBinContent(i,0.01);
1371 AliInfo(
" Settings for Ds --> K K pi, Pb-Pb collisions for 20-50 CC");
1373 fNorm =
new TH1F(
"fNorm",
"fNorm",12,0,12);
1374 for(Int_t i=1;i<=12;i++)
fNorm->SetBinContent(i,0.051);
1377 fBR =
new TH1F(
"fBR",
"fBR",12,0,12);
1378 for(Int_t i=1;i<=12;i++)
fBR->SetBinContent(i,0.045);
1381 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",12,0,12);
1382 for(Int_t i=1;i<=12;i++)
fTrackingEff->SetBinContent(i,0.15);
1385 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",12,0,12);
1386 for(Int_t i=1;i<=12;i++)
fRawYield->SetBinContent(i,0.08);
1390 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",12,0,12);
1391 for(Int_t i=1;i<=12;i++)
fCutsEff->SetBinContent(i,0.20);
1394 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",12,0,12);
1395 for(Int_t i=1;i<=12;i++)
fPIDEff->SetBinContent(i,0.07);
1398 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",12,0,12);
1399 for(Int_t i=1; i<=12; i++)
fMCPtShape->SetBinContent(i,0.01);
1408 AliInfo(
" Settings for Ds --> K K pi, p-Pb collisions at 5.023 TeV");
1411 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
1412 for(Int_t i=3;i<=12;i++)
fNorm->SetBinContent(i,0.00);
1415 fBR =
new TH1F(
"fBR",
"fBR",12,0,12);
1416 for(Int_t i=3;i<=12;i++)
fBR->SetBinContent(i,0.053);
1419 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
1420 for(Int_t i=3;i<=12;i++)
fTrackingEff->SetBinContent(i,0.09);
1423 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
1424 for(Int_t i=3; i<=4; i++)
fRawYield->SetBinContent(i,0.10);
1425 for(Int_t i=5; i<=6; i++)
fRawYield->SetBinContent(i,0.08);
1426 for(Int_t i=7; i<=8; i++)
fRawYield->SetBinContent(i,0.05);
1427 for(Int_t i=9; i<=12; i++)
fRawYield->SetBinContent(i,0.08);
1430 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
1431 for(Int_t i=3;i<=12;i++)
fCutsEff->SetBinContent(i,0.15);
1434 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
1435 for(Int_t i=3;i<=4;i++)
fPIDEff->SetBinContent(i,0.10);
1436 for(Int_t i=5;i<=12;i++)
fPIDEff->SetBinContent(i,0.05);
1439 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
1440 for(Int_t i=3;i<=12;i++)
fMCPtShape->SetBinContent(i,0.04);
1462 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
1463 for(Int_t i=1;i<=20;i++)
fNorm->SetBinContent(i,0.05);
1466 fBR =
new TH1F(
"fBR",
"fBR",20,0,20);
1467 for(Int_t i=1;i<=20;i++)
fBR->SetBinContent(i,0.021);
1470 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",20,0,20);
1471 for(Int_t i=1;i<=20;i++)
fTrackingEff->SetBinContent(i,0.15);
1474 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
1475 for(Int_t i=1;i<=20;i++)
fRawYield->SetBinContent(i,.10);
1486 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
1487 for(Int_t i=1;i<=20;i++)
fCutsEff->SetBinContent(i,0.15);
1490 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
1491 for(Int_t i=1;i<=20;i++)
fPIDEff->SetBinContent(i,0.05);
1494 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",20,0,20);
1495 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0.);
1496 for(Int_t iBin=7; iBin<=8; iBin++)
fMCPtShape->SetBinContent(iBin,0.01);
1497 for(Int_t iBin=9; iBin<=12; iBin++)
fMCPtShape->SetBinContent(iBin,0.05);
1498 for(Int_t iBin=13; iBin<=16; iBin++)
fMCPtShape->SetBinContent(iBin,0.05);
1522 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
1523 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.07);
1526 fBR =
new TH1F(
"fBR",
"fBR",20,0,20);
1527 for(Int_t i=1;i<=20;i++)
fBR->SetBinContent(i,0.021);
1530 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",20,0,20);
1531 for(Int_t i=1;i<=20;i++)
fTrackingEff->SetBinContent(i,0.15);
1535 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
1544 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.15);
1545 for(Int_t i=13;i<=20;i++)
fRawYield->SetBinContent(i,1);
1548 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
1549 for(Int_t i=1;i<=20;i++)
fCutsEff->SetBinContent(i,0.10);
1552 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
1553 for(Int_t i=1;i<=20;i++)
fPIDEff->SetBinContent(i,0.05);
1554 fPIDEff->SetBinContent(3,0.13);
1558 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",20,0,20);
1559 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0);
1560 for(Int_t iBin=4; iBin<=8; iBin++)
fMCPtShape->SetBinContent(iBin,0.01);
1561 for(Int_t iBin=9; iBin<=12; iBin++)
fMCPtShape->SetBinContent(iBin,0.03);
1562 for(Int_t iBin=13; iBin<=16; iBin++)
fMCPtShape->SetBinContent(iBin,0.03);
1583 AliInfo(
" Settings for D+ --> K pi pi p-p collisions at 2.76 TeV");
1586 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
1587 for(Int_t i=1;i<=20;i++)
fNorm->SetBinContent(i,0.019);
1590 fBR =
new TH1F(
"fBR",
"fBR",20,0,20);
1591 for(Int_t i=1;i<=20;i++)
fBR->SetBinContent(i,0.021);
1594 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",20,0,20);
1595 for(Int_t i=1;i<=20;i++)
fTrackingEff->SetBinContent(i,0.15);
1598 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
1601 for(Int_t i=3;i<=6;i++)
fRawYield->SetBinContent(i,0.10);
1604 for(Int_t i=9;i<=20;i++)
fRawYield->SetBinContent(i,0.055);
1607 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
1608 for(Int_t i=1;i<=20;i++)
fCutsEff->SetBinContent(i,0.15);
1611 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
1612 for(Int_t i=1;i<=20;i++)
fPIDEff->SetBinContent(i,0.05);
1613 fPIDEff->SetBinContent(3,0.10);
1614 fPIDEff->SetBinContent(4,0.10);
1617 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",20,0,20);
1618 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0.01);
1629 AliInfo(
" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV");
1631 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
1632 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.02);
1634 fBR =
new TH1F(
"fBR",
"fBR",20,0,24);
1635 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.021);
1638 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
1639 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.09);
1642 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
1651 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
1652 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
1653 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
1656 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
1665 for(Int_t i=9;i<=12;i++)
fCutsEff->SetBinContent(i,0.05);
1666 for(Int_t i=13;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
1667 for(Int_t i=17;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
1670 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
1672 fPIDEff->SetBinContent(2,0.1);
1673 for(Int_t i=3;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
1676 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
1677 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0.0);
1700 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
1701 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.035);
1704 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
1705 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.015);
1708 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
1713 for(Int_t i=4;i<=24;i++)
fTrackingEff->SetBinContent(i,0.12);
1717 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
1726 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.04);
1727 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
1728 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.2);
1731 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
1733 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.10);
1736 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
1737 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.04);
1741 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
1742 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0);
1757 AliInfo(
" Settings for D*+ --> D0 pi p-p collisions at 2.76 TeV");
1760 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
1761 for(Int_t i=1;i<=20;i++)
fNorm->SetBinContent(i,0.019);
1764 fBR =
new TH1F(
"fBR",
"fBR",20,0,20);
1765 for(Int_t i=1;i<=20;i++)
fBR->SetBinContent(i,0.015);
1768 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",20,0,20);
1769 for(Int_t i=1;i<=20;i++)
fTrackingEff->SetBinContent(i,0.15);
1772 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
1785 for(Int_t i=9;i<=20;i++)
fRawYield->SetBinContent(i,0.065);
1788 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
1789 for(Int_t i=1;i<=20;i++)
fCutsEff->SetBinContent(i,0.10);
1802 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
1803 for(Int_t i=1;i<=20;i++)
fPIDEff->SetBinContent(i,0.05);
1806 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",20,0,20);
1807 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0.01);
1821 AliInfo(
" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-20 centrality - DUMMY");
1824 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
1825 for(Int_t i=1;i<=20;i++)
fNorm->SetBinContent(i,0.05);
1828 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
1829 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.015);
1832 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);;
1833 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.15);
1837 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
1838 for(Int_t i=1;i<=24;i++)
fRawYield->SetBinContent(i,0.1);
1845 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
1846 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.10);
1852 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
1853 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.05);
1857 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
1858 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0.045);
1879 AliInfo(
" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-7.5 centrality - DUMMY");
1883 fNorm =
new TH1F(
"fNorm",
"fNorm",36,0,36);
1884 for(Int_t i=1;i<=36;i++)
fNorm->SetBinContent(i,0.048);
1887 fBR =
new TH1F(
"fBR",
"fBR",36,0,36);
1888 for(Int_t i=1;i<=36;i++)
fBR->SetBinContent(i,0.015);
1891 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",36,0,36);;
1892 for(Int_t i=1;i<=36;i++)
fTrackingEff->SetBinContent(i,0.15);
1896 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",36,0,36);
1897 for(Int_t i=1;i<=36;i++)
fRawYield->SetBinContent(i,0.05);
1902 for(Int_t i=25;i<=36;i++)
fRawYield->SetBinContent(i,0.15);
1905 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",36,0,36);
1906 for(Int_t i=1;i<=36;i++)
fCutsEff->SetBinContent(i,0.10);
1909 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",36,0,36);
1910 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.05);
1911 fPIDEff->SetBinContent(4,0.09);
1914 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",36,0,36);
1915 for(Int_t i=1;i<=36;i++)
fMCPtShape->SetBinContent(i,0.035);
1938 AliInfo(
" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 20-40 centrality - DUMMY");
1941 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
1942 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.10);
1945 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
1946 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.015);
1949 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);;
1950 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.15);
1954 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
1955 for(Int_t i=1;i<=24;i++)
fRawYield->SetBinContent(i,0.15);
1958 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
1959 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.10);
1962 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
1963 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.04);
1967 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
1968 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0.);
1983 AliInfo(
" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 40-80 centrality - DUMMY");
1986 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
1987 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.07);
1990 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
1991 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.015);
1994 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);;
1995 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.15);
1999 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
2000 for(Int_t i=1;i<=24;i++)
fRawYield->SetBinContent(i,0.2);
2020 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
2021 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.10);
2024 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
2025 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.05);
2029 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
2030 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
2050 AliInfo(
" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-10% centrality ");
2054 fNorm =
new TH1F(
"fNorm",
"fNorm",36,0,36);
2055 for(Int_t i=1;i<=36;i++)
fNorm->SetBinContent(i,0.048);
2058 fBR =
new TH1F(
"fBR",
"fBR",36,0,36);
2059 for(Int_t i=1;i<=36;i++)
fBR->SetBinContent(i,0.015);
2062 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",36,0,36);;
2063 for(Int_t i=1;i<=36;i++)
fTrackingEff->SetBinContent(i,0.15);
2067 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",36,0,36);
2068 for(Int_t i=1;i<=3;i++)
fRawYield->SetBinContent(i, 0.);
2072 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2073 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
2074 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.05);
2075 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.05);
2076 for(Int_t i=25;i<=36;i++)
fRawYield->SetBinContent(i,0.10);
2079 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",36,0,36);
2080 for(Int_t i=1;i<=3;i++)
fCutsEff->SetBinContent(i, 0.);
2082 for(Int_t i=5;i<=36;i++)
fCutsEff->SetBinContent(i,0.10);
2085 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",36,0,36);
2086 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.05);
2089 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",36,0,36);
2090 for(Int_t i=1;i<=3;i++)
fMCPtShape->SetBinContent(i, 0.);
2093 for(Int_t i=5;i<=36;i++)
fMCPtShape->SetBinContent(i,0.01);
2106 AliInfo(
" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 30-50 centrality - DUMMY");
2109 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
2110 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.07);
2113 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
2114 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.015);
2117 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);;
2118 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.15);
2122 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
2123 for(Int_t i=1;i<=24;i++)
fRawYield->SetBinContent(i,0.1);
2143 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
2144 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
2152 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
2153 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.05);
2157 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
2158 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
2173 AliInfo(
" Settings for D* --> D0 pi, p-Pb collisions at 5.023 TeV");
2176 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
2177 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.07);
2180 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
2181 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.015);
2184 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
2185 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.09);
2188 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
2197 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
2198 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.06);
2199 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.06);
2202 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
2205 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
2208 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
2209 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
2212 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
2213 for(Int_t i=1;i<=3;i++)
fMCPtShape->SetBinContent(i,0.03);
2214 for(Int_t i=4;i<=6;i++)
fMCPtShape->SetBinContent(i,0.);
2215 for(Int_t i=7;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
2234 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
2235 for(Int_t i=3;i<=5;i++)
fMCPtShape->SetBinContent(i,0.17);
2236 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2242 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
2243 for(Int_t i=3;i<=5;i++)
fMCPtShape->SetBinContent(i,0.17);
2244 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2250 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
2251 for(Int_t i=3;i<=5;i++)
fMCPtShape->SetBinContent(i,0.14);
2252 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2258 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
2259 for(Int_t i=3;i<=5;i++)
fMCPtShape->SetBinContent(i,0.11);
2260 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2266 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.08);
2267 for(Int_t i=3;i<=5;i++)
fMCPtShape->SetBinContent(i,0.08);
2268 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2274 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.18);
2275 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.09);
2282 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.23);
2283 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2289 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.08);
2290 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.095);
2296 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.08);
2297 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2303 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.15);
2304 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.07);
2311 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.16);
2312 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.15);
2318 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
2319 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.15);
2325 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.10);
2326 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2332 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.10);
2333 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.045);
2339 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.10);
2340 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.045);
2385 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.15);
2386 for(Int_t i=5;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2387 for(Int_t i=9;i<=17;i++)
fRawYield->SetBinContent(i,0.15);
2389 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2390 for(Int_t i=5;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2391 for(Int_t i=9;i<=17;i++)
fMCPtShape->SetBinContent(i,0.05);
2393 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.15);
2394 for(Int_t i=5;i<=8;i++)
fCutsEff->SetBinContent(i,0.20);
2395 for(Int_t i=9;i<=17;i++)
fCutsEff->SetBinContent(i,0.10);
2403 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.10);
2404 for(Int_t i=5;i<=8;i++)
fRawYield->SetBinContent(i,0.10);
2405 for(Int_t i=9;i<=17;i++)
fRawYield->SetBinContent(i,0.15);
2407 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2408 for(Int_t i=5;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2409 for(Int_t i=9;i<=17;i++)
fMCPtShape->SetBinContent(i,0.05);
2411 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.15);
2412 for(Int_t i=5;i<=8;i++)
fCutsEff->SetBinContent(i,0.20);
2413 for(Int_t i=9;i<=17;i++)
fCutsEff->SetBinContent(i,0.15);
2421 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.10);
2422 for(Int_t i=5;i<=8;i++)
fRawYield->SetBinContent(i,0.10);
2423 for(Int_t i=9;i<=17;i++)
fRawYield->SetBinContent(i,0.08);
2425 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2426 for(Int_t i=5;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2427 for(Int_t i=9;i<=17;i++)
fMCPtShape->SetBinContent(i,0.05);
2429 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.20);
2430 for(Int_t i=5;i<=8;i++)
fCutsEff->SetBinContent(i,0.15);
2431 for(Int_t i=9;i<=17;i++)
fCutsEff->SetBinContent(i,0.20);
2439 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.20);
2440 for(Int_t i=5;i<=8;i++)
fRawYield->SetBinContent(i,0.10);
2441 for(Int_t i=9;i<=17;i++)
fRawYield->SetBinContent(i,0.20);
2443 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2444 for(Int_t i=5;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2445 for(Int_t i=9;i<=17;i++)
fMCPtShape->SetBinContent(i,0.05);
2447 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.20);
2448 for(Int_t i=5;i<=8;i++)
fCutsEff->SetBinContent(i,0.20);
2449 for(Int_t i=9;i<=17;i++)
fCutsEff->SetBinContent(i,0.20);
2457 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.10);
2458 for(Int_t i=5;i<=8;i++)
fRawYield->SetBinContent(i,0.10);
2459 for(Int_t i=9;i<=17;i++)
fRawYield->SetBinContent(i,0.05);
2461 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2462 for(Int_t i=5;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2463 for(Int_t i=9;i<=17;i++)
fMCPtShape->SetBinContent(i,0.05);
2465 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.20);
2466 for(Int_t i=5;i<=8;i++)
fCutsEff->SetBinContent(i,0.15);
2467 for(Int_t i=9;i<=17;i++)
fCutsEff->SetBinContent(i,0.20);
2475 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.05);
2476 for(Int_t i=5;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2477 for(Int_t i=9;i<=17;i++)
fRawYield->SetBinContent(i,0.10);
2479 for(Int_t i=1;i<=17;i++)
fMCPtShape->SetBinContent(i,0.05);
2480 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.20);
2481 for(Int_t i=5;i<=8;i++)
fCutsEff->SetBinContent(i,0.10);
2482 for(Int_t i=9;i<=17;i++)
fCutsEff->SetBinContent(i,0.20);
2526 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
2528 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.04);
2529 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.06);
2530 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.04);
2532 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2534 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2535 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2536 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2538 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2539 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2540 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2548 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
2550 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.04);
2551 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.08);
2552 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.04);
2554 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2556 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2557 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2558 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2560 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2561 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2562 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2568 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
2570 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.04);
2571 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.12);
2572 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
2574 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2576 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2577 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2578 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2580 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2581 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2582 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2589 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
2591 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.04);
2592 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.03);
2593 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
2595 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2597 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2598 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2599 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2601 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2602 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2603 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2610 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
2612 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.04);
2613 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.03);
2614 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.04);
2616 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2618 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2619 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2620 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2622 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2623 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2624 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2632 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
2634 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.07);
2635 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2636 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.10);
2638 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2640 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.07);
2641 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.08);
2642 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.08);
2644 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2645 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.10);
2646 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.10);
2652 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
2654 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.05);
2655 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2656 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.10);
2658 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2660 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.07);
2661 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.08);
2662 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.08);
2664 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2665 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2666 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2672 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
2674 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.05);
2675 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2676 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.05);
2678 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2680 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.07);
2681 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.08);
2682 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.08);
2684 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2685 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2686 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2692 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < 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.05);
2696 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.05);
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.07);
2701 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.08);
2702 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.08);
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.05);
2706 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2713 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
2715 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.07);
2716 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2717 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.07);
2719 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2721 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.07);
2722 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.08);
2723 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.08);
2725 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2726 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2727 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2736 AliInfo(
"Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
2738 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.06);
2739 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.03);
2740 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.07);
2742 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2744 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2745 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2746 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2748 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2749 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2750 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2757 AliInfo(
"Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
2759 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.04);
2760 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2761 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.05);
2763 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2765 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2766 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2767 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2769 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2770 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2771 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2778 AliInfo(
"Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
2780 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.05);
2781 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.03);
2782 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.05);
2784 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2786 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2787 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2788 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2790 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2791 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2792 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2799 AliInfo(
"Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
2801 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.04);
2802 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.03);
2803 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.05);
2805 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2807 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2808 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2809 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2811 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2812 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2813 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2819 AliInfo(
"Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
2821 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.06);
2822 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.07);
2823 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.07);
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.05);
2828 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2829 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2831 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2832 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2833 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2841 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
2847 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
2848 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
2849 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
2850 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
2851 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
2852 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
2857 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
2860 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
2866 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
2876 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
2883 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
2884 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
2885 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
2886 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
2887 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
2888 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
2893 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
2897 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
2903 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
2911 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
2917 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
2918 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
2919 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
2920 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
2921 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
2922 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
2927 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
2931 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
2937 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
2946 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
2951 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
2952 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
2953 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
2954 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
2955 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
2956 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
2962 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
2966 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
2974 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.02);
2984 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult ZNA");
2989 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
2990 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
2991 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
2992 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
2993 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
2994 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
2999 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3003 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3012 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3021 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult ZNA");
3026 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3027 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3028 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3029 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3030 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3031 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3037 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3040 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3049 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3057 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult ZNA");
3062 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3063 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3064 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3065 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3066 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3067 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3072 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3075 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3084 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3094 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult ZNA");
3099 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3100 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3101 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3102 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3103 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3104 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3109 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3112 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
3120 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3129 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult CL1");
3135 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3136 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3137 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3138 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3139 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3140 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3145 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3148 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3154 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3163 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult CL1");
3168 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3169 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3170 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3171 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3172 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3173 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3179 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3183 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3189 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3196 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult CL1");
3202 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3203 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3204 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3205 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3206 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3207 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3212 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3216 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3222 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);}
3228 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult CL1");
3233 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3234 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3235 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3236 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3237 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3238 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3244 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3248 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
3256 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.02);
3264 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
3275 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3276 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3277 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
3282 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3285 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3292 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3302 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
3313 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3314 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3315 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
3321 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3325 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3331 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3348 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
3359 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3360 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3361 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
3366 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3371 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3377 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3395 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
3406 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3407 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3408 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
3414 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3418 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
3427 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.02);
3444 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
3455 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3456 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3457 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
3462 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3466 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.);
3475 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3494 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
3505 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3506 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3507 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
3513 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3517 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
3526 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3544 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
3555 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3556 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3557 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
3564 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3568 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
3577 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3598 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
3609 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3610 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3611 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
3616 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3620 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
3628 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3646 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult CL1");
3657 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3658 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3659 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
3664 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3667 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3673 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3683 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult CL1");
3694 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3695 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3696 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
3701 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3705 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3712 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3729 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 40-60% mult CL1");
3740 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3741 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3742 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.0);
3747 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3751 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3759 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3777 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult CL1");
3788 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3789 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.00);
3790 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.0);
3795 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3799 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
3807 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.02);
3824 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
3835 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3836 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3837 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
3842 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3849 for(Int_t i=4;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3852 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3861 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
3873 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3874 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3875 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
3881 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3885 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3892 for(Int_t i=4;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3909 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
3922 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3923 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3924 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
3932 for(Int_t i=3;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
3936 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3943 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3961 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
3973 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3974 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3975 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
3982 for(Int_t i=3;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
3986 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
3995 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.02);
4013 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 0-20% mult ZNA");
4025 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4026 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4027 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
4033 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
4037 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.);
4045 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
4062 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 20-40% mult ZNA");
4074 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4075 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4076 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
4082 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
4087 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
4095 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
4114 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 40-60% mult ZNA");
4126 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4127 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4128 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
4135 for(Int_t i=3;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
4140 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
4148 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
4167 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 60-100% mult ZNA");
4179 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4180 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4181 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
4187 for(Int_t i=3;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
4191 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
4199 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
4217 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 0-20% mult CL1");
4228 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4229 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4230 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.0);
4235 for(Int_t i=3;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
4241 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
4244 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
4253 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 20-40% mult CL1");
4265 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4266 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4267 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.0);
4274 for(Int_t i=3;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
4278 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
4284 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
4301 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 40-60% mult CL1");
4313 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4314 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4315 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.0);
4321 for(Int_t i=3;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
4325 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
4331 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
4350 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 60-100% mult CL1");
4362 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4363 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4364 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
4371 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
4375 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
4384 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.02);
4406 fNorm =
new TH1F(
"fNorm",
"fNorm",4,2,6);
4407 for(Int_t i=1;i<=4;i++)
fNorm->SetBinContent(i,0.035);
4410 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",4,2,6);
4411 for(Int_t i=1;i<=4;i++)
fTrackingEff->SetBinContent(i,0.12);
4414 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
4420 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",4,2,6);
4427 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",4,2,4);
4428 fPIDEff->SetBinContent(1,0.13);
4429 fPIDEff->SetBinContent(2,0.30);
4430 fPIDEff->SetBinContent(3,0.5);
4431 fPIDEff->SetBinContent(4,0.30);
4434 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",12,0,12);
4435 for(Int_t i=1; i<=2; i++)
fMCPtShape->SetBinContent(i,1.);
4436 for(Int_t i=3; i<=4; i++)
fMCPtShape->SetBinContent(i,0.03);
4437 for(Int_t i=5; i<=6; i++)
fMCPtShape->SetBinContent(i,0.03);
4438 for(Int_t i=7; i<=8; i++)
fMCPtShape->SetBinContent(i,0.02);
4439 for(Int_t i=9; i<=12; i++)
fMCPtShape->SetBinContent(i,0.02);
4451 fNorm =
new TH1F(
"fNorm",
"fNorm",12,0.,12.);
4452 for(Int_t i=1;i<=12;i++)
fNorm->SetBinContent(i,0.05);
4455 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",12,0.,12.);
4456 for(Int_t i=1;i<=12;i++)
fTrackingEff->SetBinContent(i,0.1);
4459 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",12,0.,12.);
4460 for(Int_t i=1;i<=2;i++)
fRawYield->SetBinContent(i,1.);
4461 for(Int_t i=3;i<=12;i++)
fRawYield->SetBinContent(i,0.1);
4463 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",12,0.,12.);
4464 for(Int_t i=1;i<=12;i++)
fCutsEff->SetBinContent(i,0.1);
4467 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",12,0.,12.);
4468 for(Int_t i=1;i<=12;i++)
fPIDEff->SetBinContent(i,0.1);
4471 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",12,0.,12.);
4472 for(Int_t i=1; i<=2;i++)
fMCPtShape->SetBinContent(i,1.);
4473 for(Int_t i=3; i<=12;i++)
fMCPtShape->SetBinContent(i,0.1);
4485 return fCutsEff->GetBinContent(bin);
4505 return TMath::Sqrt(err);
4513 Int_t bin=
fPIDEff->FindBin(pt);
4515 return fPIDEff->GetBinContent(bin);
4563 err += feeddownErr*feeddownErr;
4565 return TMath::Sqrt(err);
4572 gStyle->SetOptStat(0);
4574 TCanvas *cSystErr =
new TCanvas(
"cSystErr",
"Systematic Errors",300,80,640,500);
4575 cSystErr->Range(0.20,-0.5,18.4,0.34);
4576 cSystErr->SetRightMargin(0.318);
4577 cSystErr->SetFillColor(0);
4579 TH2F *hFrame =
new TH2F(
"hFrame",
"Systematic errors; p_{t} (GeV/c); Relative Error",40,0,40,100,-1,+1);
4580 hFrame->SetAxisRange(1.,35.9,
"X");
4582 hFrame->SetAxisRange(-0.5,0.5,
"Y");
4585 TLegend *leg =
new TLegend(0.69,0.44,0.98,0.86,NULL,
"brNDC");
4586 leg->SetTextSize(0.03601695);
4587 leg->SetFillStyle(0);
4588 leg->SetBorderSize(0);
4590 TH1F *hTotErr=
new TH1F(
"hTotErr",
"",36,0,36);
4592 TGraphAsymmErrors *gTotErr =
new TGraphAsymmErrors(nbins);
4593 for(Int_t i=1;i<=36;i++) {
4594 Double_t pt = hTotErr->GetBinCenter(i);
4595 Double_t ptwidth = hTotErr->GetBinWidth(i);
4598 Double_t x=0., y=0., errxl=0., errxh=0., erryl=0., erryh=0.;
4599 Double_t toterryl=0., toterryh=0.;
4600 for(Int_t j=0; j<grErrFeeddown->GetN(); j++) {
4601 grErrFeeddown->GetPoint(j,x,y);
4602 errxh = grErrFeeddown->GetErrorXhigh(j);
4603 errxl = grErrFeeddown->GetErrorXlow(j);
4604 if ( ( (x-errxl) <= pt) && ( (x+errxl) >= pt) ) {
4605 erryh = grErrFeeddown->GetErrorYhigh(j);
4606 erryl = grErrFeeddown->GetErrorYlow(j);
4614 hTotErr->SetBinContent(i,toterryh);
4615 gTotErr->SetPoint(i,pt,0.);
4616 gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,toterryl,toterryh);
4620 gTotErr->SetPoint(i,pt,0.);
4625 gTotErr->SetLineColor(kBlack);
4626 gTotErr->SetFillColor(kRed);
4627 gTotErr->SetFillStyle(3002);
4629 leg->AddEntry(gTotErr,
"Total (excl. norm.)",
"f");
4636 fNorm->SetFillColor(1);
4637 fNorm->SetFillStyle(3002);
4641 Double_t norm =
fNorm->GetBinContent(1)*100;
4642 leg->AddEntry(
fNorm,Form(
"Normalization (%.1f%s)",norm,
"%"),
"");
4645 grErrFeeddown->SetFillColor(kTeal-8);
4646 grErrFeeddown->SetFillStyle(3001);
4647 grErrFeeddown->Draw(
"2");
4648 leg->AddEntry(grErrFeeddown,
"Feed-down from B",
"f");
4655 hTrackingEffRefl->Draw(
"same");
4659 fBR->SetFillColor(6);
4660 fBR->SetFillStyle(3005);
4664 hBRRefl->Draw(
"same");
4665 leg->AddEntry(
fBR,
"Branching ratio",
"f");
4669 ci = TColor::GetColor(
"#00cc00");
4675 hRawYieldRefl->Draw(
"same");
4676 leg->AddEntry(
fRawYield,
"Yield extraction",
"l");
4683 hCutsEffRefl->Draw(
"same");
4684 leg->AddEntry(
fCutsEff,
"Cuts efficiency",
"l");
4691 hPIDEffRefl->Draw(
"same");
4692 leg->AddEntry(
fPIDEff,
"PID efficiency",
"l");
4695 Int_t ci = TColor::GetColor(
"#9933ff");
4701 hMCPtShapeRefl->Draw(
"same");
4702 leg->AddEntry(
fMCPtShape,
"MC p_{t} shape",
"l");
4705 Int_t ci = TColor::GetColor(
"#ff6600");
4711 hPartAntipartRefl->Draw(
"same");
4718 cSystErr->SaveAs(
"RelativeSystematics.eps");
4727 TH1F *hout=(TH1F*)hin->Clone(
"hout");
void InitDstartoD0pi2010PbPb2040CentScan()
void InitDstartoD0pi2011PbPb2030CentScan()
void InitDplustoKpipi2011PbPb3040CentScan()
void InitDstartoD0pi2013pPb2040ZNA()
ClassImp(AliAnalysisTaskTriggerRates) AliAnalysisTaskTriggerRates
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()
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 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 Init(Int_t decay)
Function to initialize the variables/histograms.
void InitD0toKpi2010PbPb1020CentScan()
void InitDstartoD0pi2011PbPb010CentScan()
Double_t GetRawYieldErr(Double_t pt) const
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()
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 low pt analysis (no topological cuts)
void InitDplustoKpipi2013pPb60100V0A()
void InitD0toKpi2011PbPb3050InPlane()
void InitDplustoKpipi2010PbPb5080CentScan()
Double_t GetPIDEffErr(Double_t pt) const