36 #include <TGraphAsymmErrors.h>
63 fCentralityClass(
"0100"),
64 fRapidityRange(
"0101"),
66 fIsLowPtAnalysis(false),
67 fIsPass4Analysis(false),
69 fIsRapidityScan(false)
101 AliFatal(
"Only settings for 2010 and the low energy runs are implemented so far");
120 else AliFatal(
"Not yet implemented");
134 else AliFatal(
"Not yet implemented");
165 else AliFatal(
"Not yet implemented");
185 else AliFatal(
"Not yet implemented");
197 else AliFatal(
"Not yet implemented");
228 else AliFatal(
"Not yet implemented");
246 else AliFatal(
"Not yet implemented");
258 else AliFatal(
"Not yet implemented");
287 else AliFatal(
"Not yet implemented");
296 else AliFatal(
"Not yet implemented");
301 else AliFatal(
"Not yet implemented");
305 else AliFatal(
"Not yet implemented");
309 else AliFatal(
"Not yet implemented");
313 printf(
"Invalid decay type: %d\n",decay);
327 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
328 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.035);
331 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
332 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.012);
335 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
336 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.08);
339 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
343 for(Int_t i=4;i<=7;i++)
fRawYield->SetBinContent(i,0.04);
344 for(Int_t i=8;i<=12;i++)
fRawYield->SetBinContent(i,0.07);
345 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.10);
346 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,1);
349 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
350 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.10);
353 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
354 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.03);
355 fPIDEff->SetBinContent(2,0.05);
358 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
359 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0);
378 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
379 for(Int_t i=1;i<=20;i++)
fNorm->SetBinContent(i,0.05);
382 fBR =
new TH1F(
"fBR",
"fBR",20,0,20);
383 for(Int_t i=1;i<=20;i++)
fBR->SetBinContent(i,0.012);
386 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",20,0,20);
387 for(Int_t i=1;i<=20;i++)
fTrackingEff->SetBinContent(i,0.10);
390 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
394 for(Int_t i=4;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
395 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.10);
398 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
403 for(Int_t i=5;i<=16;i++)
fCutsEff->SetBinContent(i,0.10);
404 for(Int_t i=17;i<=20;i++)
fCutsEff->SetBinContent(i,0.);
407 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
408 for(Int_t i=3;i<=16;i++)
fPIDEff->SetBinContent(i,0.05);
411 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",20,0,20);
412 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0.01);
415 for(Int_t i=13;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
436 fNorm =
new TH1F(
"fNorm",
"fNorm",36,0,36);
437 for(Int_t i=1;i<36;i++)
fNorm->SetBinContent(i,0.048);
440 fBR =
new TH1F(
"fBR",
"fBR",36,0,36);
441 for(Int_t i=1;i<=36;i++)
fBR->SetBinContent(i,0.012);
444 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",36,0,36);
445 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.10);
446 for(Int_t i=25;i<=36;i++)
fTrackingEff->SetBinContent(i,0.);
449 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",36,0,36);
453 for(Int_t i=4;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
454 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.10);
455 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.30);
456 for(Int_t i=25;i<=36;i++)
fRawYield->SetBinContent(i,0.);
459 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",36,0,36);
465 for(Int_t i=6;i<=24;i++)
fCutsEff->SetBinContent(i,0.06);
466 for(Int_t i=25;i<=36;i++)
fCutsEff->SetBinContent(i,0.0);
469 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",36,0,36);
470 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.05);
471 for(Int_t i=25;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
474 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",36,0,36);
475 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
479 for(Int_t i=13;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
480 for(Int_t i=17;i<=24;i++)
fMCPtShape->SetBinContent(i,0.05);
481 for(Int_t i=25;i<=36;i++)
fMCPtShape->SetBinContent(i,0.0);
507 for(Int_t i=5;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
508 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
509 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.15);
510 for(Int_t i=17;i<=36;i++)
fRawYield->SetBinContent(i,0.);
520 for(Int_t i=7;i<=8;i++)
fCutsEff->SetBinContent(i,0.15);
521 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.15);
522 for(Int_t i=25;i<=36;i++)
fCutsEff->SetBinContent(i,0.0);
537 for(Int_t i=4;i<=6;i++)
fRawYield->SetBinContent(i,0.07);
538 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
539 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.10);
540 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.15);
541 for(Int_t i=17;i<=36;i++)
fRawYield->SetBinContent(i,0.);
551 for(Int_t i=7;i<=8;i++)
fCutsEff->SetBinContent(i,0.15);
552 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.15);
553 for(Int_t i=17;i<=36;i++)
fCutsEff->SetBinContent(i,0.0);
564 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
565 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.07);
568 fBR =
new TH1F(
"fBR",
"fBR",20,0,20);
569 for(Int_t i=1;i<=20;i++)
fBR->SetBinContent(i,0.012);
572 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",20,0,20);
573 for(Int_t i=1;i<=20;i++)
fTrackingEff->SetBinContent(i,0.10);
577 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
580 for(Int_t i=3;i<=16;i++)
fRawYield->SetBinContent(i,0.05);
584 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
589 for(Int_t i=5;i<=16;i++)
fCutsEff->SetBinContent(i,0.10);
590 for(Int_t i=17;i<=20;i++)
fCutsEff->SetBinContent(i,0.);
593 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
596 for(Int_t i=3;i<=16;i++)
fPIDEff->SetBinContent(i,0.05);
599 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",20,0,20);
600 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0.01);
603 for(Int_t i=13;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
619 AliInfo(
" Settings for D0 --> Kpi Pb-Pb collisions at 2.76 TeV - 0-10 centrality - newww!!");
622 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
623 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.048);
626 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
627 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.012);
630 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
631 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.10);
635 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
642 for (Int_t i=7; i<=8; i++)
fRawYield->SetBinContent(i,0.06);
643 for (Int_t i=9; i<=12; i++)
fRawYield->SetBinContent(i,0.07);
644 for (Int_t i=13; i<=16; i++)
fRawYield->SetBinContent(i,0.11);
645 for (Int_t i=17; i<=24; i++)
fRawYield->SetBinContent(i,0.30);
647 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
654 for (Int_t i=7; i<=8; i++)
fCutsEff->SetBinContent(i,0.05);
655 for (Int_t i=9; i<=12; i++)
fCutsEff->SetBinContent(i,0.05);
656 for (Int_t i=13; i<=16; i++)
fCutsEff->SetBinContent(i,0.05);
657 for (Int_t i=17; i<=24; i++)
fCutsEff->SetBinContent(i,0.05);
661 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
664 for(Int_t i=2;i<=16;i++)
fPIDEff->SetBinContent(i,0.05);
667 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
672 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.);
689 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
690 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.07);
693 fBR =
new TH1F(
"fBR",
"fBR",20,0,20);
694 for(Int_t i=1;i<=20;i++)
fBR->SetBinContent(i,0.012);
697 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",20,0,20);
698 for(Int_t i=1;i<=20;i++)
fTrackingEff->SetBinContent(i,0.10);
702 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
716 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.112);
717 for(Int_t i=17;i<=20;i++)
fRawYield->SetBinContent(i,0);
719 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
726 for(Int_t i=7;i<=8;i++)
fCutsEff->SetBinContent(i,0.15);
727 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.15);
731 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
734 for(Int_t i=2;i<=16;i++)
fPIDEff->SetBinContent(i,0.05);
737 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",0,0,20);
738 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0.01);
742 for(Int_t i=13;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
743 for(Int_t i=17;i<=20;i++)
fMCPtShape->SetBinContent(i,0.);
758 AliInfo(
" Settings for D0 --> K pi, p-p collisions at 2.76 TeV");
761 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
762 for(Int_t i=1;i<=20;i++)
fNorm->SetBinContent(i,0.019);
765 fBR =
new TH1F(
"fBR",
"fBR",20,0,20);
766 for(Int_t i=1;i<=20;i++)
fBR->SetBinContent(i,0.012);
769 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",20,0,20);
770 for(Int_t i=1;i<=20;i++)
fTrackingEff->SetBinContent(i,0.10);
773 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
775 for(Int_t i=1;i<=20;i++)
fRawYield->SetBinContent(i,0.15);
778 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
779 for(Int_t i=1;i<=20;i++)
fCutsEff->SetBinContent(i,0.10);
781 for(Int_t i=7;i<=20;i++)
fCutsEff->SetBinContent(i,0.15);
785 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
786 for(Int_t i=1;i<=20;i++)
fPIDEff->SetBinContent(i,0.15);
788 for(Int_t i=7;i<=20;i++)
fPIDEff->SetBinContent(i,0.05);
791 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",20,0,20);
792 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0.01);
812 AliInfo(
" Settings for D0 --> K pi, pp collisions at 7 TeV, analysis without topological cuts");
815 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
816 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.035);
819 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
820 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.0129);
823 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
824 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.08);
827 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
828 for(Int_t i=1;i<=24;i++)
fRawYield->SetBinContent(i,0.2);
837 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.2);
840 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
841 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.);
844 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
845 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.03);
846 fPIDEff->SetBinContent(1,0.05);
847 fPIDEff->SetBinContent(2,0.05);
850 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
851 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0);
864 AliInfo(
" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV");
866 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
867 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.02);
870 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
871 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.013);
874 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
875 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.06);
878 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
881 for(Int_t i=3;i<=12;i++)
fRawYield->SetBinContent(i,0.04);
882 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.07);
883 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.11);
886 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
889 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
892 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
893 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
896 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
899 for(Int_t i=3;i<=24;i++)
fMCPtShape->SetBinContent(i,0.0);
917 AliInfo(
" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV, analysis without topological cuts");
920 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
921 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.035);
924 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
925 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.0129);
928 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
929 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.06);
932 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
933 for(Int_t i=1;i<=24;i++)
fRawYield->SetBinContent(i,0.2);
942 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.11);
945 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
946 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.);
949 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
950 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.);
953 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
954 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0);
974 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
975 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.035);
978 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
979 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.021);
982 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
983 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.12);
987 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
994 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
995 for(Int_t i=12;i<=24;i++)
fRawYield->SetBinContent(i,0.10);
998 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
999 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.10);
1002 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
1003 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.05);
1004 fPIDEff->SetBinContent(1,0.15);
1005 fPIDEff->SetBinContent(2,0.15);
1006 fPIDEff->SetBinContent(3,0.15);
1007 fPIDEff->SetBinContent(4,0.15);
1008 for(Int_t i=12;i<=16;i++)
fPIDEff->SetBinContent(i,0.10);
1011 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
1012 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0);
1037 fNorm =
new TH1F(
"fNorm",
"fNorm",36,0,36);
1038 for(Int_t i=1;i<=36;i++)
fNorm->SetBinContent(i,0.035);
1041 fBR =
new TH1F(
"fBR",
"fBR",36,0,36);
1042 for(Int_t i=1;i<=36;i++)
fBR->SetBinContent(i,0.021);
1045 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",36,0,36);
1046 for(Int_t i=1;i<=36;i++)
fTrackingEff->SetBinContent(i,0.12);
1050 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",36,0,36);
1054 for(Int_t i=4;i<=12;i++)
fRawYield->SetBinContent(i,0.03);
1055 for(Int_t i=13;i<=24;i++)
fRawYield->SetBinContent(i,0.07);
1056 for(Int_t i=25;i<=36;i++)
fRawYield->SetBinContent(i,0.20);
1059 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",36,0,36);
1060 for(Int_t i=1;i<=36;i++)
fCutsEff->SetBinContent(i,0.05);
1063 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",36,0,36);
1064 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.05);
1067 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",36,0,36);
1068 for(Int_t i=1;i<=36;i++)
fMCPtShape->SetBinContent(i,0);
1092 fNorm =
new TH1F(
"fNorm",
"fNorm",12,0,12);
1093 for(Int_t i=1;i<=12;i++)
fNorm->SetBinContent(i,0.035);
1096 fBR =
new TH1F(
"fBR",
"fBR",12,0,12);
1097 for(Int_t i=1;i<=12;i++)
fBR->SetBinContent(i,0.06);
1100 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",12,0,12);
1101 for(Int_t i=1;i<=12;i++)
fTrackingEff->SetBinContent(i,0.12);
1105 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",12,0,12);
1120 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",12,0,12);
1121 for(Int_t i=1;i<=12;i++)
fCutsEff->SetBinContent(i,0.15);
1124 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",12,0,12);
1125 for(Int_t i=1;i<=12;i++)
fPIDEff->SetBinContent(i,0.07);
1128 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",12,0,12);
1129 for(Int_t i=1; i<=2; i++)
fMCPtShape->SetBinContent(i,1.);
1130 for(Int_t i=3; i<=4; i++)
fMCPtShape->SetBinContent(i,0.03);
1131 for(Int_t i=5; i<=6; i++)
fMCPtShape->SetBinContent(i,0.03);
1132 for(Int_t i=7; i<=8; i++)
fMCPtShape->SetBinContent(i,0.02);
1133 for(Int_t i=9; i<=12; i++)
fMCPtShape->SetBinContent(i,0.02);
1156 fNorm =
new TH1F(
"fNorm",
"fNorm",36,0,36);
1157 for(Int_t i=1;i<=36;i++)
fNorm->SetBinContent(i,0.048);
1160 fBR =
new TH1F(
"fBR",
"fBR",36,0,36);
1161 for(Int_t i=1;i<=36;i++)
fBR->SetBinContent(i,0.021);
1164 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",36,0,36);
1165 for(Int_t i=1;i<=36;i++)
fTrackingEff->SetBinContent(i,0.15);
1168 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",36,0,36);
1169 for(Int_t i=1;i<=36;i++)
fRawYield->SetBinContent(i,.10);
1173 for(Int_t i=7; i<=8; i++)
fRawYield->SetBinContent(i,0.10);
1174 for(Int_t i=9; i<=12; i++)
fRawYield->SetBinContent(i,0.08);
1175 for(Int_t i=13; i<=16; i++)
fRawYield->SetBinContent(i,0.05);
1176 for(Int_t i=17; i<=24; i++)
fRawYield->SetBinContent(i,0.08);
1177 for(Int_t i=25; i<=36; i++)
fRawYield->SetBinContent(i,0.20);
1181 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",36,0,36);
1182 for(Int_t i=1; i<=12; i++)
fCutsEff->SetBinContent(i,0.10);
1183 for(Int_t i=13; i<=36; i++)
fCutsEff->SetBinContent(i,0.08);
1187 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",36,0,36);
1188 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.05);
1191 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",36,0,36);
1192 for(Int_t iBin=1; iBin<=5; iBin++)
fMCPtShape->SetBinContent(iBin,0.05);
1193 for(Int_t iBin=6; iBin<=36; iBin++)
fMCPtShape->SetBinContent(iBin,0.03);
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<=3;i++)
fRawYield->SetBinContent(i,0.);
1233 for(Int_t i=7; i<=8; i++)
fRawYield->SetBinContent(i,0.08);
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.08);
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.08);
1241 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",36,0,36);
1242 for(Int_t i=1; i<=36; i++)
fCutsEff->SetBinContent(i,0.10);
1246 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",36,0,36);
1247 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.05);
1250 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",36,0,36);
1251 for(Int_t i=1;i<=3;i++)
fMCPtShape->SetBinContent(i,0.);
1254 for (Int_t i=6; i<=36; i++)
fMCPtShape->SetBinContent(i,0.01);
1276 AliInfo(
" Settings for D+ --> K pi pi, Pb-Pb collisions for 30-50 CC");
1278 fNorm =
new TH1F(
"fNorm",
"fNorm",16,0,16);
1279 for(Int_t i=1;i<=16;i++)
fNorm->SetBinContent(i,0.048);
1282 fBR =
new TH1F(
"fBR",
"fBR",16,0,16);
1283 for(Int_t i=1;i<=16;i++)
fBR->SetBinContent(i,0.021);
1286 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",16,0,16);
1287 for(Int_t i=1;i<=36;i++)
fTrackingEff->SetBinContent(i,0.15);
1290 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",16,0,16);
1291 for(Int_t i=1;i<=12;i++)
fRawYield->SetBinContent(i,0.09);
1292 for(Int_t i=13; i<=16; i++)
fRawYield->SetBinContent(i,0.12);
1296 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",16,0,16);
1297 for(Int_t i=1; i<=12; i++)
fCutsEff->SetBinContent(i,0.10);
1298 for(Int_t i=12; i<=16; i++)
fCutsEff->SetBinContent(i,0.14);
1301 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",16,0,16);
1302 for(Int_t i=1;i<=16;i++)
fPIDEff->SetBinContent(i,0.05);
1305 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",16,0,16);
1308 for(Int_t iBin=5; iBin<=16; iBin++)
fMCPtShape->SetBinContent(iBin,0.02);
1330 AliInfo(
" Settings for Ds --> K K pi, Pb-Pb collisions for 0-7.5 CC");
1333 fNorm =
new TH1F(
"fNorm",
"fNorm",12,0,12);
1334 for(Int_t i=1;i<=12;i++)
fNorm->SetBinContent(i,0.048);
1337 fBR =
new TH1F(
"fBR",
"fBR",12,0,12);
1338 for(Int_t i=1;i<=12;i++)
fBR->SetBinContent(i,0.053);
1341 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",12,0,12);
1342 for(Int_t i=1;i<=12;i++)
fTrackingEff->SetBinContent(i,0.15);
1345 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",12,0,12);
1346 for(Int_t i=1;i<=6;i++)
fRawYield->SetBinContent(i,.30);
1347 for(Int_t i=7; i<=12; i++)
fRawYield->SetBinContent(i,0.20);
1352 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",12,0,12);
1353 for(Int_t i=1;i<=12;i++)
fCutsEff->SetBinContent(i,0.20);
1356 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",12,0,12);
1357 for(Int_t i=1;i<=12;i++)
fPIDEff->SetBinContent(i,0.1);
1360 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",12,0,12);
1361 for(Int_t i=1; i<=2; i++)
fMCPtShape->SetBinContent(i,1.);
1362 for(Int_t i=3; i<=4; i++)
fMCPtShape->SetBinContent(i,0.03);
1363 for(Int_t i=5; i<=6; i++)
fMCPtShape->SetBinContent(i,0.03);
1364 for(Int_t i=7; i<=8; i++)
fMCPtShape->SetBinContent(i,0.02);
1365 for(Int_t i=9; i<=12; i++)
fMCPtShape->SetBinContent(i,0.02);
1386 AliInfo(
" Settings for Ds --> K K pi, Pb-Pb collisions for 0-10 CC");
1388 fNorm =
new TH1F(
"fNorm",
"fNorm",12,0,12);
1389 for(Int_t i=1;i<=12;i++)
fNorm->SetBinContent(i,0.048);
1392 fBR =
new TH1F(
"fBR",
"fBR",12,0,12);
1393 for(Int_t i=1;i<=12;i++)
fBR->SetBinContent(i,0.045);
1396 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",12,0,12);
1397 for(Int_t i=1;i<=12;i++)
fTrackingEff->SetBinContent(i,0.15);
1400 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",12,0,12);
1401 for(Int_t i=1;i<=12;i++)
fRawYield->SetBinContent(i,0.08);
1405 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",12,0,12);
1406 for(Int_t i=1;i<=12;i++)
fCutsEff->SetBinContent(i,0.20);
1409 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",12,0,12);
1410 for(Int_t i=1;i<=12;i++)
fPIDEff->SetBinContent(i,0.07);
1413 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",12,0,12);
1414 for(Int_t i=1; i<=6; i++)
fMCPtShape->SetBinContent(i,0.02);
1415 for(Int_t i=7; i<=8; i++)
fMCPtShape->SetBinContent(i,0.01);
1416 for(Int_t i=9; i<=12; i++)
fMCPtShape->SetBinContent(i,0.01);
1428 AliInfo(
" Settings for Ds --> K K pi, Pb-Pb collisions for 20-50 CC");
1430 fNorm =
new TH1F(
"fNorm",
"fNorm",12,0,12);
1431 for(Int_t i=1;i<=12;i++)
fNorm->SetBinContent(i,0.051);
1434 fBR =
new TH1F(
"fBR",
"fBR",12,0,12);
1435 for(Int_t i=1;i<=12;i++)
fBR->SetBinContent(i,0.045);
1438 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",12,0,12);
1439 for(Int_t i=1;i<=12;i++)
fTrackingEff->SetBinContent(i,0.15);
1442 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",12,0,12);
1443 for(Int_t i=1;i<=12;i++)
fRawYield->SetBinContent(i,0.08);
1447 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",12,0,12);
1448 for(Int_t i=1;i<=12;i++)
fCutsEff->SetBinContent(i,0.20);
1451 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",12,0,12);
1452 for(Int_t i=1;i<=12;i++)
fPIDEff->SetBinContent(i,0.07);
1455 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",12,0,12);
1456 for(Int_t i=1; i<=12; i++)
fMCPtShape->SetBinContent(i,0.01);
1465 AliInfo(
" Settings for Ds --> K K pi, p-Pb collisions at 5.023 TeV");
1468 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
1469 for(Int_t i=3;i<=12;i++)
fNorm->SetBinContent(i,0.00);
1472 fBR =
new TH1F(
"fBR",
"fBR",12,0,12);
1473 for(Int_t i=3;i<=12;i++)
fBR->SetBinContent(i,0.053);
1476 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
1477 for(Int_t i=3;i<=12;i++)
fTrackingEff->SetBinContent(i,0.09);
1480 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
1481 for(Int_t i=3; i<=4; i++)
fRawYield->SetBinContent(i,0.10);
1482 for(Int_t i=5; i<=6; i++)
fRawYield->SetBinContent(i,0.08);
1483 for(Int_t i=7; i<=8; i++)
fRawYield->SetBinContent(i,0.05);
1484 for(Int_t i=9; i<=12; i++)
fRawYield->SetBinContent(i,0.08);
1487 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
1488 for(Int_t i=3;i<=12;i++)
fCutsEff->SetBinContent(i,0.15);
1491 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
1492 for(Int_t i=3;i<=4;i++)
fPIDEff->SetBinContent(i,0.10);
1493 for(Int_t i=5;i<=12;i++)
fPIDEff->SetBinContent(i,0.05);
1496 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
1497 for(Int_t i=3;i<=12;i++)
fMCPtShape->SetBinContent(i,0.04);
1519 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
1520 for(Int_t i=1;i<=20;i++)
fNorm->SetBinContent(i,0.05);
1523 fBR =
new TH1F(
"fBR",
"fBR",20,0,20);
1524 for(Int_t i=1;i<=20;i++)
fBR->SetBinContent(i,0.021);
1527 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",20,0,20);
1528 for(Int_t i=1;i<=20;i++)
fTrackingEff->SetBinContent(i,0.15);
1531 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
1532 for(Int_t i=1;i<=20;i++)
fRawYield->SetBinContent(i,.10);
1543 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
1544 for(Int_t i=1;i<=20;i++)
fCutsEff->SetBinContent(i,0.15);
1547 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
1548 for(Int_t i=1;i<=20;i++)
fPIDEff->SetBinContent(i,0.05);
1551 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",20,0,20);
1552 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0.);
1553 for(Int_t iBin=7; iBin<=8; iBin++)
fMCPtShape->SetBinContent(iBin,0.01);
1554 for(Int_t iBin=9; iBin<=12; iBin++)
fMCPtShape->SetBinContent(iBin,0.05);
1555 for(Int_t iBin=13; iBin<=16; iBin++)
fMCPtShape->SetBinContent(iBin,0.05);
1579 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
1580 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.07);
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);
1592 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
1601 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.15);
1602 for(Int_t i=13;i<=20;i++)
fRawYield->SetBinContent(i,1);
1605 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
1606 for(Int_t i=1;i<=20;i++)
fCutsEff->SetBinContent(i,0.10);
1609 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
1610 for(Int_t i=1;i<=20;i++)
fPIDEff->SetBinContent(i,0.05);
1611 fPIDEff->SetBinContent(3,0.13);
1615 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",20,0,20);
1616 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0);
1617 for(Int_t iBin=4; iBin<=8; iBin++)
fMCPtShape->SetBinContent(iBin,0.01);
1618 for(Int_t iBin=9; iBin<=12; iBin++)
fMCPtShape->SetBinContent(iBin,0.03);
1619 for(Int_t iBin=13; iBin<=16; iBin++)
fMCPtShape->SetBinContent(iBin,0.03);
1640 AliInfo(
" Settings for D+ --> K pi pi p-p collisions at 2.76 TeV");
1643 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
1644 for(Int_t i=1;i<=20;i++)
fNorm->SetBinContent(i,0.019);
1647 fBR =
new TH1F(
"fBR",
"fBR",20,0,20);
1648 for(Int_t i=1;i<=20;i++)
fBR->SetBinContent(i,0.021);
1651 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",20,0,20);
1652 for(Int_t i=1;i<=20;i++)
fTrackingEff->SetBinContent(i,0.15);
1655 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
1658 for(Int_t i=3;i<=6;i++)
fRawYield->SetBinContent(i,0.10);
1661 for(Int_t i=9;i<=20;i++)
fRawYield->SetBinContent(i,0.055);
1664 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
1665 for(Int_t i=1;i<=20;i++)
fCutsEff->SetBinContent(i,0.15);
1668 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
1669 for(Int_t i=1;i<=20;i++)
fPIDEff->SetBinContent(i,0.05);
1670 fPIDEff->SetBinContent(3,0.10);
1671 fPIDEff->SetBinContent(4,0.10);
1674 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",20,0,20);
1675 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0.01);
1686 AliInfo(
" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV");
1688 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
1689 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.02);
1691 fBR =
new TH1F(
"fBR",
"fBR",20,0,24);
1692 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.021);
1695 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
1696 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.09);
1699 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
1708 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
1709 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
1710 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
1713 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
1722 for(Int_t i=9;i<=12;i++)
fCutsEff->SetBinContent(i,0.05);
1723 for(Int_t i=13;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
1724 for(Int_t i=17;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
1727 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
1729 fPIDEff->SetBinContent(2,0.1);
1730 for(Int_t i=3;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
1733 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
1734 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0.0);
1757 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
1758 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.035);
1761 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
1762 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.015);
1765 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
1770 for(Int_t i=4;i<=24;i++)
fTrackingEff->SetBinContent(i,0.12);
1774 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
1783 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.04);
1784 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
1785 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.2);
1788 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
1790 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.10);
1793 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
1794 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.04);
1798 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
1799 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0);
1814 AliInfo(
" Settings for D*+ --> D0 pi p-p collisions at 2.76 TeV");
1817 fNorm =
new TH1F(
"fNorm",
"fNorm",20,0,20);
1818 for(Int_t i=1;i<=20;i++)
fNorm->SetBinContent(i,0.019);
1821 fBR =
new TH1F(
"fBR",
"fBR",20,0,20);
1822 for(Int_t i=1;i<=20;i++)
fBR->SetBinContent(i,0.015);
1825 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",20,0,20);
1826 for(Int_t i=1;i<=20;i++)
fTrackingEff->SetBinContent(i,0.15);
1829 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",20,0,20);
1842 for(Int_t i=9;i<=20;i++)
fRawYield->SetBinContent(i,0.065);
1845 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",20,0,20);
1846 for(Int_t i=1;i<=20;i++)
fCutsEff->SetBinContent(i,0.10);
1859 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",20,0,20);
1860 for(Int_t i=1;i<=20;i++)
fPIDEff->SetBinContent(i,0.05);
1863 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",20,0,20);
1864 for(Int_t i=1;i<=20;i++)
fMCPtShape->SetBinContent(i,0.01);
1878 AliInfo(
" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-20 centrality - DUMMY");
1881 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
1882 for(Int_t i=1;i<=20;i++)
fNorm->SetBinContent(i,0.05);
1885 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
1886 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.015);
1889 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);;
1890 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.15);
1894 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
1895 for(Int_t i=1;i<=24;i++)
fRawYield->SetBinContent(i,0.1);
1902 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
1903 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.10);
1909 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
1910 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.05);
1914 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
1915 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0.045);
1936 AliInfo(
" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-7.5 centrality - DUMMY");
1940 fNorm =
new TH1F(
"fNorm",
"fNorm",36,0,36);
1941 for(Int_t i=1;i<=36;i++)
fNorm->SetBinContent(i,0.048);
1944 fBR =
new TH1F(
"fBR",
"fBR",36,0,36);
1945 for(Int_t i=1;i<=36;i++)
fBR->SetBinContent(i,0.015);
1948 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",36,0,36);;
1949 for(Int_t i=1;i<=36;i++)
fTrackingEff->SetBinContent(i,0.15);
1953 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",36,0,36);
1954 for(Int_t i=1;i<=36;i++)
fRawYield->SetBinContent(i,0.05);
1959 for(Int_t i=25;i<=36;i++)
fRawYield->SetBinContent(i,0.15);
1962 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",36,0,36);
1963 for(Int_t i=1;i<=36;i++)
fCutsEff->SetBinContent(i,0.10);
1966 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",36,0,36);
1967 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.05);
1968 fPIDEff->SetBinContent(4,0.09);
1971 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",36,0,36);
1972 for(Int_t i=1;i<=36;i++)
fMCPtShape->SetBinContent(i,0.035);
1995 AliInfo(
" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 20-40 centrality - DUMMY");
1998 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
1999 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.10);
2002 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
2003 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.015);
2006 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);;
2007 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.15);
2011 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
2012 for(Int_t i=1;i<=24;i++)
fRawYield->SetBinContent(i,0.15);
2015 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
2016 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.10);
2019 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
2020 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.04);
2024 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
2025 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0.);
2040 AliInfo(
" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 40-80 centrality - DUMMY");
2043 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
2044 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.07);
2047 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
2048 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.015);
2051 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);;
2052 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.15);
2056 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
2057 for(Int_t i=1;i<=24;i++)
fRawYield->SetBinContent(i,0.2);
2077 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
2078 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.10);
2081 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
2082 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.05);
2086 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
2087 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
2107 AliInfo(
" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-10% centrality ");
2111 fNorm =
new TH1F(
"fNorm",
"fNorm",36,0,36);
2112 for(Int_t i=1;i<=36;i++)
fNorm->SetBinContent(i,0.048);
2115 fBR =
new TH1F(
"fBR",
"fBR",36,0,36);
2116 for(Int_t i=1;i<=36;i++)
fBR->SetBinContent(i,0.015);
2119 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",36,0,36);;
2120 for(Int_t i=1;i<=36;i++)
fTrackingEff->SetBinContent(i,0.15);
2124 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",36,0,36);
2125 for(Int_t i=1;i<=3;i++)
fRawYield->SetBinContent(i, 0.);
2129 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2130 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
2131 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.05);
2132 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.05);
2133 for(Int_t i=25;i<=36;i++)
fRawYield->SetBinContent(i,0.10);
2136 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",36,0,36);
2137 for(Int_t i=1;i<=3;i++)
fCutsEff->SetBinContent(i, 0.);
2139 for(Int_t i=5;i<=36;i++)
fCutsEff->SetBinContent(i,0.10);
2142 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",36,0,36);
2143 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.05);
2146 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",36,0,36);
2147 for(Int_t i=1;i<=3;i++)
fMCPtShape->SetBinContent(i, 0.);
2150 for(Int_t i=5;i<=36;i++)
fMCPtShape->SetBinContent(i,0.01);
2163 AliInfo(
" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 30-50 centrality - DUMMY");
2166 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
2167 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.07);
2170 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
2171 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.015);
2174 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);;
2175 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.15);
2179 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
2180 for(Int_t i=1;i<=24;i++)
fRawYield->SetBinContent(i,0.1);
2200 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
2201 for(Int_t i=1;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
2209 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
2210 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.05);
2214 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
2215 for(Int_t i=1;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
2230 AliInfo(
" Settings for D* --> D0 pi, p-Pb collisions at 5.023 TeV");
2233 fNorm =
new TH1F(
"fNorm",
"fNorm",24,0,24);
2234 for(Int_t i=1;i<=24;i++)
fNorm->SetBinContent(i,0.07);
2237 fBR =
new TH1F(
"fBR",
"fBR",24,0,24);
2238 for(Int_t i=1;i<=24;i++)
fBR->SetBinContent(i,0.015);
2241 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",24,0,24);
2242 for(Int_t i=1;i<=24;i++)
fTrackingEff->SetBinContent(i,0.09);
2245 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
2254 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
2255 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.06);
2256 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.06);
2259 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",24,0,24);
2262 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
2265 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",24,0,24);
2266 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
2269 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",24,0,24);
2270 for(Int_t i=1;i<=3;i++)
fMCPtShape->SetBinContent(i,0.03);
2271 for(Int_t i=4;i<=6;i++)
fMCPtShape->SetBinContent(i,0.);
2272 for(Int_t i=7;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
2291 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
2292 for(Int_t i=3;i<=5;i++)
fMCPtShape->SetBinContent(i,0.17);
2293 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2299 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
2300 for(Int_t i=3;i<=5;i++)
fMCPtShape->SetBinContent(i,0.17);
2301 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2307 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
2308 for(Int_t i=3;i<=5;i++)
fMCPtShape->SetBinContent(i,0.14);
2309 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2315 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
2316 for(Int_t i=3;i<=5;i++)
fMCPtShape->SetBinContent(i,0.11);
2317 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2323 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.08);
2324 for(Int_t i=3;i<=5;i++)
fMCPtShape->SetBinContent(i,0.08);
2325 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2331 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.18);
2332 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.09);
2339 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.23);
2340 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2346 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.08);
2347 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.095);
2353 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.08);
2354 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2360 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.15);
2361 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.07);
2368 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.16);
2369 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.15);
2375 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.05);
2376 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.15);
2382 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.10);
2383 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.08);
2389 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.10);
2390 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.045);
2396 for(Int_t i=7;i<=12;i++)
fRawYield->SetBinContent(i,0.10);
2397 for(Int_t i=7;i<=12;i++)
fMCPtShape->SetBinContent(i,0.045);
2442 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.15);
2443 for(Int_t i=5;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2444 for(Int_t i=9;i<=17;i++)
fRawYield->SetBinContent(i,0.15);
2446 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2447 for(Int_t i=5;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2448 for(Int_t i=9;i<=17;i++)
fMCPtShape->SetBinContent(i,0.05);
2450 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.15);
2451 for(Int_t i=5;i<=8;i++)
fCutsEff->SetBinContent(i,0.20);
2452 for(Int_t i=9;i<=17;i++)
fCutsEff->SetBinContent(i,0.10);
2460 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.10);
2461 for(Int_t i=5;i<=8;i++)
fRawYield->SetBinContent(i,0.10);
2462 for(Int_t i=9;i<=17;i++)
fRawYield->SetBinContent(i,0.15);
2464 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2465 for(Int_t i=5;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2466 for(Int_t i=9;i<=17;i++)
fMCPtShape->SetBinContent(i,0.05);
2468 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.15);
2469 for(Int_t i=5;i<=8;i++)
fCutsEff->SetBinContent(i,0.20);
2470 for(Int_t i=9;i<=17;i++)
fCutsEff->SetBinContent(i,0.15);
2478 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.10);
2479 for(Int_t i=5;i<=8;i++)
fRawYield->SetBinContent(i,0.10);
2480 for(Int_t i=9;i<=17;i++)
fRawYield->SetBinContent(i,0.08);
2482 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2483 for(Int_t i=5;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2484 for(Int_t i=9;i<=17;i++)
fMCPtShape->SetBinContent(i,0.05);
2486 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.20);
2487 for(Int_t i=5;i<=8;i++)
fCutsEff->SetBinContent(i,0.15);
2488 for(Int_t i=9;i<=17;i++)
fCutsEff->SetBinContent(i,0.20);
2496 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.20);
2497 for(Int_t i=5;i<=8;i++)
fRawYield->SetBinContent(i,0.10);
2498 for(Int_t i=9;i<=17;i++)
fRawYield->SetBinContent(i,0.20);
2500 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2501 for(Int_t i=5;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2502 for(Int_t i=9;i<=17;i++)
fMCPtShape->SetBinContent(i,0.05);
2504 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.20);
2505 for(Int_t i=5;i<=8;i++)
fCutsEff->SetBinContent(i,0.20);
2506 for(Int_t i=9;i<=17;i++)
fCutsEff->SetBinContent(i,0.20);
2514 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.10);
2515 for(Int_t i=5;i<=8;i++)
fRawYield->SetBinContent(i,0.10);
2516 for(Int_t i=9;i<=17;i++)
fRawYield->SetBinContent(i,0.05);
2518 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2519 for(Int_t i=5;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2520 for(Int_t i=9;i<=17;i++)
fMCPtShape->SetBinContent(i,0.05);
2522 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.20);
2523 for(Int_t i=5;i<=8;i++)
fCutsEff->SetBinContent(i,0.15);
2524 for(Int_t i=9;i<=17;i++)
fCutsEff->SetBinContent(i,0.20);
2532 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.05);
2533 for(Int_t i=5;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2534 for(Int_t i=9;i<=17;i++)
fRawYield->SetBinContent(i,0.10);
2536 for(Int_t i=1;i<=17;i++)
fMCPtShape->SetBinContent(i,0.05);
2537 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.20);
2538 for(Int_t i=5;i<=8;i++)
fCutsEff->SetBinContent(i,0.10);
2539 for(Int_t i=9;i<=17;i++)
fCutsEff->SetBinContent(i,0.20);
2583 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
2585 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.04);
2586 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.06);
2587 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.04);
2589 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2591 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2592 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2593 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2595 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2596 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2597 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2605 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
2607 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.04);
2608 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.08);
2609 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.04);
2611 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2613 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2614 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2615 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2617 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2618 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2619 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2625 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
2627 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.04);
2628 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.12);
2629 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
2631 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2633 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2634 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2635 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2637 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2638 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2639 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2646 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
2648 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.04);
2649 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.03);
2650 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
2652 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2654 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2655 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2656 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2658 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2659 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2660 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2667 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
2669 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.04);
2670 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.03);
2671 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.04);
2673 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2675 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2676 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2677 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2679 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2680 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2681 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2689 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
2691 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.07);
2692 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2693 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.10);
2695 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2697 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.07);
2698 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.08);
2699 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.08);
2701 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2702 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.10);
2703 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.10);
2709 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
2711 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.05);
2712 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2713 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.10);
2715 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2717 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.07);
2718 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.08);
2719 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.08);
2721 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2722 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2723 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2729 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
2731 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.05);
2732 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2733 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.05);
2735 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2737 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.07);
2738 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.08);
2739 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.08);
2741 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2742 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2743 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2749 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
2751 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.05);
2752 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2753 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.05);
2755 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2757 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.07);
2758 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.08);
2759 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.08);
2761 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2762 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2763 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2770 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
2772 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.07);
2773 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2774 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.07);
2776 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2778 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.07);
2779 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.08);
2780 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.08);
2782 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2783 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2784 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2793 AliInfo(
"Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
2795 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.06);
2796 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.03);
2797 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.07);
2799 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2801 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2802 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2803 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2805 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2806 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2807 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2814 AliInfo(
"Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
2816 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.04);
2817 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.05);
2818 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.05);
2820 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2822 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2823 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2824 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2826 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2827 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2828 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2835 AliInfo(
"Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
2837 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.05);
2838 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.03);
2839 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.05);
2841 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2843 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2844 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2845 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2847 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2848 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2849 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2856 AliInfo(
"Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
2858 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.04);
2859 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.03);
2860 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.05);
2862 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2864 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2865 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2866 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2868 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2869 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2870 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2876 AliInfo(
"Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
2878 for(Int_t i=2;i<=5;i++)
fRawYield->SetBinContent(i,0.06);
2879 for(Int_t i=6;i<=8;i++)
fRawYield->SetBinContent(i,0.07);
2880 for(Int_t i=9;i<=16;i++)
fRawYield->SetBinContent(i,0.07);
2882 for(Int_t i=1;i<=36;i++)
fPIDEff->SetBinContent(i,0.0);
2884 for(Int_t i=2;i<=5;i++)
fCutsEff->SetBinContent(i,0.05);
2885 for(Int_t i=6;i<=8;i++)
fCutsEff->SetBinContent(i,0.05);
2886 for(Int_t i=9;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
2888 for(Int_t i=2;i<=5;i++)
fMCPtShape->SetBinContent(i,0.10);
2889 for(Int_t i=6;i<=8;i++)
fMCPtShape->SetBinContent(i,0.05);
2890 for(Int_t i=9;i<=16;i++)
fMCPtShape->SetBinContent(i,0.05);
2898 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
2904 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
2905 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
2906 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
2907 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
2908 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
2909 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
2914 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
2917 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
2923 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
2933 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
2940 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
2941 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
2942 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
2943 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
2944 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
2945 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
2950 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
2954 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
2960 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
2968 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
2974 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
2975 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
2976 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
2977 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
2978 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
2979 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
2984 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
2988 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
2994 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3003 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
3008 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3009 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3010 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3011 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3012 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3013 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3019 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3023 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
3031 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.02);
3041 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult ZNA");
3046 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3047 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3048 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3049 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3050 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3051 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3056 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3060 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3069 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3078 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult ZNA");
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);
3094 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);
3106 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3114 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult ZNA");
3119 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3120 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3121 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3122 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3123 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3124 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3129 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.00);
3141 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3151 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult ZNA");
3156 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3157 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3158 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3159 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3160 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3161 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3166 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.0);
3177 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3186 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult CL1");
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);
3202 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3205 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3211 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3220 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult CL1");
3225 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3226 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3227 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3228 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3229 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3230 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3236 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3240 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3246 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3253 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult CL1");
3259 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3260 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3261 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3262 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3263 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3264 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3269 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3273 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3279 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);}
3285 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult CL1");
3290 for(Int_t i=3;i<=4;i++)
fRawYield->SetBinContent(i,0.042);
3291 for(Int_t i=5;i<=6;i++)
fRawYield->SetBinContent(i,0.042);
3292 for(Int_t i=7;i<=8;i++)
fRawYield->SetBinContent(i,0.058);
3293 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.058);
3294 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.09);
3295 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.09);
3301 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3305 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
3313 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.02);
3321 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
3332 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3333 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3334 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
3339 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3342 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3349 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3359 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
3370 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3371 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3372 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
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);
3405 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
3416 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3417 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3418 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
3423 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3428 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3434 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3452 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
3463 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3464 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3465 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
3471 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3475 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
3484 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.02);
3501 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
3512 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3513 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3514 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
3519 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3523 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.);
3532 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3551 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
3562 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3563 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3564 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
3570 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3574 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
3583 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3601 AliInfo(
"Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
3612 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3613 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3614 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
3621 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3625 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
3634 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3655 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
3666 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3667 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3668 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
3673 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3677 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
3685 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
3703 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult CL1");
3714 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3715 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3716 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
3721 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3724 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3730 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3740 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult CL1");
3751 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3752 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3753 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
3758 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3762 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3769 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3786 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 40-60% mult CL1");
3797 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3798 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3799 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.0);
3804 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3808 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3816 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3834 AliInfo(
"Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult CL1");
3845 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3846 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.00);
3847 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.0);
3852 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3856 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
3864 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.02);
3881 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
3892 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3893 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3894 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
3899 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3906 for(Int_t i=4;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3909 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3918 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
3930 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3931 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3932 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
3938 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
3942 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
3949 for(Int_t i=4;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
3966 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
3979 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
3980 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
3981 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
3989 for(Int_t i=3;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
3993 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
4000 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
4018 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
4030 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4031 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4032 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
4039 for(Int_t i=3;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
4043 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
4052 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.02);
4070 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 0-20% mult ZNA");
4082 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4083 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4084 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
4090 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
4094 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.);
4102 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
4119 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 20-40% mult ZNA");
4131 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4132 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4133 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
4139 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
4144 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
4152 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
4171 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 40-60% mult ZNA");
4183 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4184 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4185 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
4192 for(Int_t i=3;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
4197 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
4205 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
4224 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 60-100% mult ZNA");
4236 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4237 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4238 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.00);
4244 for(Int_t i=3;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
4248 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
4256 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.03);
4274 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 0-20% mult CL1");
4285 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4286 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4287 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.0);
4292 for(Int_t i=3;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
4298 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
4301 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
4310 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 20-40% mult CL1");
4322 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4323 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4324 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.0);
4331 for(Int_t i=3;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
4335 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
4341 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
4358 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 40-60% mult CL1");
4370 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4371 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4372 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.0);
4378 for(Int_t i=3;i<=16;i++)
fCutsEff->SetBinContent(i,0.05);
4382 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.00);
4388 for(Int_t i=5;i<=24;i++)
fMCPtShape->SetBinContent(i,0.01);
4407 AliInfo(
"Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 60-100% mult CL1");
4419 for(Int_t i=9;i<=12;i++)
fRawYield->SetBinContent(i,0.06);
4420 for(Int_t i=13;i<=16;i++)
fRawYield->SetBinContent(i,0.08);
4421 for(Int_t i=17;i<=24;i++)
fRawYield->SetBinContent(i,0.08);
4428 for(Int_t i=3;i<=24;i++)
fCutsEff->SetBinContent(i,0.05);
4432 for(Int_t i=1;i<=24;i++)
fPIDEff->SetBinContent(i,0.0);
4441 for(Int_t i=7;i<=16;i++)
fMCPtShape->SetBinContent(i,0.02);
4463 fNorm =
new TH1F(
"fNorm",
"fNorm",4,2,6);
4464 for(Int_t i=1;i<=4;i++)
fNorm->SetBinContent(i,0.035);
4467 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",4,2,6);
4468 for(Int_t i=1;i<=4;i++)
fTrackingEff->SetBinContent(i,0.12);
4471 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",24,0,24);
4477 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",4,2,6);
4484 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",4,2,4);
4485 fPIDEff->SetBinContent(1,0.13);
4486 fPIDEff->SetBinContent(2,0.30);
4487 fPIDEff->SetBinContent(3,0.5);
4488 fPIDEff->SetBinContent(4,0.30);
4491 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",12,0,12);
4492 for(Int_t i=1; i<=2; i++)
fMCPtShape->SetBinContent(i,1.);
4493 for(Int_t i=3; i<=4; i++)
fMCPtShape->SetBinContent(i,0.03);
4494 for(Int_t i=5; i<=6; i++)
fMCPtShape->SetBinContent(i,0.03);
4495 for(Int_t i=7; i<=8; i++)
fMCPtShape->SetBinContent(i,0.02);
4496 for(Int_t i=9; i<=12; i++)
fMCPtShape->SetBinContent(i,0.02);
4508 fNorm =
new TH1F(
"fNorm",
"fNorm",12,0.,12.);
4509 for(Int_t i=1;i<=12;i++)
fNorm->SetBinContent(i,0.05);
4512 fTrackingEff =
new TH1F(
"fTrackingEff",
"fTrackingEff",12,0.,12.);
4513 for(Int_t i=1;i<=12;i++)
fTrackingEff->SetBinContent(i,0.1);
4516 fRawYield =
new TH1F(
"fRawYield",
"fRawYield",12,0.,12.);
4517 for(Int_t i=1;i<=2;i++)
fRawYield->SetBinContent(i,1.);
4518 for(Int_t i=3;i<=12;i++)
fRawYield->SetBinContent(i,0.1);
4520 fCutsEff =
new TH1F(
"fCutsEff",
"fCutsEff",12,0.,12.);
4521 for(Int_t i=1;i<=12;i++)
fCutsEff->SetBinContent(i,0.1);
4524 fPIDEff =
new TH1F(
"fPIDEff",
"fPIDEff",12,0.,12.);
4525 for(Int_t i=1;i<=12;i++)
fPIDEff->SetBinContent(i,0.1);
4528 fMCPtShape =
new TH1F(
"fMCPtShape",
"fMCPtShape",12,0.,12.);
4529 for(Int_t i=1; i<=2;i++)
fMCPtShape->SetBinContent(i,1.);
4530 for(Int_t i=3; i<=12;i++)
fMCPtShape->SetBinContent(i,0.1);
4542 return fCutsEff->GetBinContent(bin);
4562 return TMath::Sqrt(err);
4570 Int_t bin=
fPIDEff->FindBin(pt);
4572 return fPIDEff->GetBinContent(bin);
4620 err += feeddownErr*feeddownErr;
4622 return TMath::Sqrt(err);
4629 gStyle->SetOptStat(0);
4631 TCanvas *cSystErr =
new TCanvas(
"cSystErr",
"Systematic Errors",300,80,640,500);
4632 cSystErr->Range(0.20,-0.5,18.4,0.34);
4633 cSystErr->SetRightMargin(0.318);
4634 cSystErr->SetFillColor(0);
4636 TH2F *hFrame =
new TH2F(
"hFrame",
"Systematic errors; p_{t} (GeV/c); Relative Error",40,0,40,100,-1,+1);
4637 hFrame->SetAxisRange(1.,35.9,
"X");
4639 hFrame->SetAxisRange(-0.5,0.5,
"Y");
4642 TLegend *leg =
new TLegend(0.69,0.44,0.98,0.86,NULL,
"brNDC");
4643 leg->SetTextSize(0.03601695);
4644 leg->SetFillStyle(0);
4645 leg->SetBorderSize(0);
4647 TH1F *hTotErr=
new TH1F(
"hTotErr",
"",36,0,36);
4649 TGraphAsymmErrors *gTotErr =
new TGraphAsymmErrors(nbins);
4650 for(Int_t i=1;i<=36;i++) {
4651 Double_t pt = hTotErr->GetBinCenter(i);
4652 Double_t ptwidth = hTotErr->GetBinWidth(i);
4655 Double_t x=0., y=0., errxl=0., errxh=0., erryl=0., erryh=0.;
4656 Double_t toterryl=0., toterryh=0.;
4657 for(Int_t j=0; j<grErrFeeddown->GetN(); j++) {
4658 grErrFeeddown->GetPoint(j,x,y);
4659 errxh = grErrFeeddown->GetErrorXhigh(j);
4660 errxl = grErrFeeddown->GetErrorXlow(j);
4661 if ( ( (x-errxl) <= pt) && ( (x+errxl) >= pt) ) {
4662 erryh = grErrFeeddown->GetErrorYhigh(j);
4663 erryl = grErrFeeddown->GetErrorYlow(j);
4671 hTotErr->SetBinContent(i,toterryh);
4672 gTotErr->SetPoint(i,pt,0.);
4673 gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,toterryl,toterryh);
4677 gTotErr->SetPoint(i,pt,0.);
4682 gTotErr->SetLineColor(kBlack);
4683 gTotErr->SetFillColor(kRed);
4684 gTotErr->SetFillStyle(3002);
4686 leg->AddEntry(gTotErr,
"Total (excl. norm.)",
"f");
4693 fNorm->SetFillColor(1);
4694 fNorm->SetFillStyle(3002);
4698 Double_t norm =
fNorm->GetBinContent(1)*100;
4699 leg->AddEntry(
fNorm,Form(
"Normalization (%.1f%s)",norm,
"%"),
"");
4702 grErrFeeddown->SetFillColor(kTeal-8);
4703 grErrFeeddown->SetFillStyle(3001);
4704 grErrFeeddown->Draw(
"2");
4705 leg->AddEntry(grErrFeeddown,
"Feed-down from B",
"f");
4712 hTrackingEffRefl->Draw(
"same");
4716 fBR->SetFillColor(6);
4717 fBR->SetFillStyle(3005);
4721 hBRRefl->Draw(
"same");
4722 leg->AddEntry(
fBR,
"Branching ratio",
"f");
4726 ci = TColor::GetColor(
"#00cc00");
4732 hRawYieldRefl->Draw(
"same");
4733 leg->AddEntry(
fRawYield,
"Yield extraction",
"l");
4740 hCutsEffRefl->Draw(
"same");
4741 leg->AddEntry(
fCutsEff,
"Cuts efficiency",
"l");
4748 hPIDEffRefl->Draw(
"same");
4749 leg->AddEntry(
fPIDEff,
"PID efficiency",
"l");
4752 Int_t ci = TColor::GetColor(
"#9933ff");
4758 hMCPtShapeRefl->Draw(
"same");
4759 leg->AddEntry(
fMCPtShape,
"MC p_{t} shape",
"l");
4762 Int_t ci = TColor::GetColor(
"#ff6600");
4768 hPartAntipartRefl->Draw(
"same");
4775 cSystErr->SaveAs(
"RelativeSystematics.eps");
4784 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 InitDplustoKpipi2010ppPass4()
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()
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 InitDplustoKpipi2013pPb60100V0A()
void InitD0toKpi2011PbPb3050InPlane()
void InitDplustoKpipi2010PbPb5080CentScan()
Double_t GetPIDEffErr(Double_t pt) const