PWGCF/Correlations/JCORRAN/AliJEbeHistos.cxx

   1 #include  "AliJEbeHistos.h"
   2 #include  "AliJCard.h"
   3
   4 //______________________________________________________________________________
   5 AliJEbeHistos::AliJEbeHistos(AliJCard* cardP):
   6         fcard(cardP),
   7         fmaxEtaRange(0),
   8         fmaxTriggEtaRange(0),
   9         ftriggFiducCut(0),
  10         fhistoList(NULL),
  11         fUseDirectory(true),
  12         fTopDirectory(NULL)
  13 {   // constructor
  14
  15         fcard=cardP;
  16         fmaxEtaRange = fcard->Get("EtaRange");
  17         ftriggFiducCut =  fcard->Get("TriggerFiducialEtaCut"); //FK// Fiduc cut
  18         fmaxTriggEtaRange =  fmaxEtaRange - ftriggFiducCut; //FK// Trigger range
  19
  20         fhistoList = new TList();
  21         fTopDirectory = gDirectory;
  22
  23         //fhtyp[1] = "Real";
  24         //fhtyp[2] = "Mixed";
  25         //fhtyp[3] = "Rap. Gap";
  26 }
  27
  28 //______________________________________________________________________________
  29 AliJEbeHistos::AliJEbeHistos(const AliJEbeHistos& obj):
  30         fcard(obj.fcard),
  31         fmaxEtaRange(obj.fmaxEtaRange),
  32         fmaxTriggEtaRange(obj.fmaxTriggEtaRange),
  33         ftriggFiducCut(obj.ftriggFiducCut),
  34         fhistoList(obj.fhistoList),
  35         fUseDirectory(obj.fUseDirectory),
  36         fTopDirectory(obj.fTopDirectory)
  37 {
  38         // copy constructor
  39 }
  40
  41 //______________________________________________________________________________
  42 AliJEbeHistos& AliJEbeHistos::operator=(const AliJEbeHistos& obj)
  43 {
  44
  45         JUNUSED(obj);
  46         // copy constructor
  47         return *this;
  48 }
  49
  50
  51 void AliJEbeHistos::CreateUnfoldingHistos(){
  52         // pi0mass histos
  53
  54         int    hic;
  55         float  b1 = 0, b2 = 0.50;
  56         int    bins=300;
  57         double lbin = 0.0;
  58         double hbin = 0.5;
  59         const int NPTBins = 51;
  60
  61         Double_t PTBINS[NPTBins+1] = {
  62                 0.0, 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45,
  63                 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95,
  64                 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9,
  65                 2.0, 2.2, 2.4, 2.6, 2.8, 3.0, 3.2, 3.4, 3.6, 3.8,
  66                 4.0, 4.5, 5.0, 5.5, 6.0, 6.5, 7.0, 8.0, 10.5, 14.0,
  67                 21.0, 32.0};
  68
  69         if(fcard->GetNoOfBins(kCentrType) > kMaxNoCentrBin ){
  70                 cout<<"ERROR: No of Centrality bins exceed max dim in AliJEbeHistos.cxx "<<endl;
  71                 exit(0);
  72         }
  73
  74         TDirectory * cwd = gDirectory;
  75         TDirectory * nwd = gDirectory->mkdir("UnfoldingHistos");
  76         if( nwd ) nwd->cd();
  77
  78         for( hic=0;hic<fcard->GetNoOfBins(kCentrType);hic++){
  79                 for(int ih = 1; ih < kNHarmonics; ih++){
  80                         b1 = fcard->GetBinBorder(kCentrType, hic);
  81                         b2 = fcard->GetBinBorder(kCentrType, hic+1);
  82                         //inclusive fpt spectra================================================
  83                         sprintf(fhname,"hVnObsVector%02d%02d", hic, ih);
  84                         sprintf(fhtit, "Observed v_{%d}: %3.0f -%3.0f ", ih, b1, b2);
  85                         fhVnObsVector[hic][ih] = new TH1D(fhname, fhtit, bins, lbin,hbin); fhVnObsVector[hic][ih]->Sumw2();
  86                         fhistoList->Add(fhVnObsVector[hic][ih]);
  87                         sprintf(fhname,"hResponseDist%02d%02d", hic,ih);
  88                         sprintf(fhtit, "Response Distribution v_{%d}: %3.0f -%3.0f ", ih, b1, b2);
  89                         fhResponseDist[hic][ih] = new TH2D(fhname, fhtit, bins, -0.5, 0.5, bins, -0.5,0.5); fhResponseDist[hic][ih]->Sumw2();
  90                         fhistoList->Add(fhResponseDist[hic][ih]);
  91                         sprintf(fhname,"hVnObsEP%02d%02d", hic,ih);
  92                         sprintf(fhtit, "Observed EP v_{%d}: %3.0f -%3.0f ", ih, b1, b2);
  93                         fhVnObsEP[hic][ih] = new TH1D(fhname, fhtit, bins, -0.2,0.5); fhVnObsEP[hic][ih]->Sumw2();
  94                         fhistoList->Add(fhVnObsEP[hic][ih]);
  95                         sprintf(fhname,"hCosndPhiPt%02d%02d", hic,ih);
  96                         sprintf(fhtit, "Cos%d#Delta#phi Pt: %3.0f -%3.0f ", ih, b1, b2);
  97                         fhCosndPhiPt[hic][ih] = new TH1D(fhname, fhtit, NPTBins, PTBINS); fhCosndPhiPt[hic][ih]->Sumw2();
  98                         fhistoList->Add(fhCosndPhiPt[hic][ih]);
  99                         sprintf(fhname,"heCosndPhiPt%02d%02d", hic,ih);
 100                         sprintf(fhtit, "eCos%d#Delta#phi Pt: %3.0f -%3.0f ", ih, b1, b2);
 101                         fheCosndPhiPt[hic][ih] = new TH1D(fhname, fhtit, NPTBins, PTBINS); fheCosndPhiPt[hic][ih]->Sumw2();
 102                         fhistoList->Add(fheCosndPhiPt[hic][ih]);
 103                         sprintf(fhname,"hCounter%02d%02d", hic,ih);
 104                         sprintf(fhtit, "Counter Pt %02d: %3.0f -%3.0f ", ih, b1, b2);
 105                         fhCounter[hic][ih] = new TH1D(fhname, fhtit, NPTBins, PTBINS); fhCounter[hic][ih]->Sumw2();
 106                         fhistoList->Add(fhCounter[hic][ih]);
 107
 108                         sprintf(fhname,"hEventPlane%02d%02d", hic,ih);
 109                         sprintf(fhtit, "Event Planes %02d: %3.0f -%3.0f ", ih, b1, b2);
 110                         fhEventPlane[hic][ih] = new TH1D(fhname, fhtit, 100, -6.5,6.5); fhEventPlane[hic][ih]->Sumw2();
 111                         fhistoList->Add(fhEventPlane[hic][ih]);
 112
 113                         sprintf(fhtit, "Q vector from V0 %02d: %3.0f -%3.0f ", ih, b1, b2);
 114                         fhQvectorV0[hic][ih] = new TH1D(Form("hQvectorV0_%02d%02d",hic,ih),fhtit ,400, 0, 100); fhQvectorV0[hic][ih]->Sumw2();
 115                         sprintf(fhtit, "Q vector from V0A %02d: %3.0f -%3.0f ", ih, b1, b2);
 116                         fhQvectorV0A[hic][ih] = new TH1D(Form("hQvectorV0A_%02d%02d",hic,ih),fhtit ,400, 0, 100); fhQvectorV0A[hic][ih]->Sumw2();
 117                         sprintf(fhtit, "Q vector from V0C %02d: %3.0f -%3.0f ", ih, b1, b2);
 118                         fhQvectorV0C[hic][ih] = new TH1D(Form("hQvectorV0C_%02d%02d",hic,ih), fhtit ,400, 0, 100); fhQvectorV0C[hic][ih]->Sumw2();
 119                         fhistoList->Add(fhQvectorV0[hic][ih]);
 120                         fhistoList->Add(fhQvectorV0A[hic][ih]);
 121                         fhistoList->Add(fhQvectorV0C[hic][ih]);
 122
 123                         sprintf(fhtit, "Q vector correlation V0A and V0C %02d: %3.0f -%3.0f ", ih, b1, b2);
 124                         fhQvectorCorrelation[hic][ih] = new TH2D(Form("hQvectorCorrelation_%02d%02d",hic,ih),fhtit , 400, 0, 100, 400, 0 ,100); fhQvectorCorrelation[hic][ih]->Sumw2();
 125                         fhistoList->Add(fhQvectorCorrelation[hic][ih]);
 126
 127                         sprintf(fhtit, "VnObs vs V0C Qvector correlation %02d: %3.0f -%3.0f ", ih, b1, b2);
 128                         fhVnObsVsQvectorCorrelation[hic][ih] = new TH2D(Form("hVnObsVsQvectorCorrelation_%02d%02d",hic,ih),fhtit , 400, 0, 0.5, 400, 0 ,100); fhVnObsVsQvectorCorrelation[hic][ih]->Sumw2();
 129                         fhistoList->Add(fhVnObsVsQvectorCorrelation[hic][ih]);
 130                 }
 131                 sprintf(fhname,"hMultiCount%02d", hic);
 132                 sprintf(fhtit, "Multiplicity count: %3.0f -%3.0f ", b1, b2);
 133                 fhMultiCount[hic] = new TH1D( fhname, fhtit, 1000, 0, 10000 );
 134                 fhistoList->Add(fhMultiCount[hic]);
 135         }
 136         // tmp histogram for e-b-e
 137         for(int ih = 1; ih < kNHarmonics; ih++){
 138                 fhEPCosndPhi[ih] = new TH1D(Form("hEPCosndPhi_%02d",ih), Form("cos%d#Delta#phi",ih), 400, -1.1, 1.1); fhEPCosndPhi[ih]->Sumw2();
 139                 fhEPCosndPhi2[ih] = new TH1D(Form("hEPCosndPhi2_%02d",ih), Form("cos%d#Delta#phi",ih), 400, -1.1, 1.1); fhEPCosndPhi2[ih]->Sumw2();
 140                 fhistoList->Add(fhEPCosndPhi[ih]);
 141                 fhistoList->Add(fhEPCosndPhi2[ih]);
 142
 143
 144         }
 145         cwd->cd();
 146 }
 147