Implemented projections for per trigger yield (Alis)

diff --git a/PWGCF/EBYE/macros/drawBalanceFunction2DPsi.C b/PWGCF/EBYE/macros/drawBalanceFunction2DPsi.C
index 5b529d23f25051f7d1c3d990e210cb6afb690fc7..03b95d4b23f478974c0c1010a3bd4c3debb183d2 100644 (file)
--- a/PWGCF/EBYE/macros/drawBalanceFunction2DPsi.C
+++ b/PWGCF/EBYE/macros/drawBalanceFunction2DPsi.C
@@ -445,7 +445,7 @@ void draw(TList *listBF, TList *listBFShuffled, TList *listBFMixed,
     c4a->SetLeftMargin(0.15);
     gHistBalanceFunctionSubtracted->DrawCopy("colz");
 
-    fitbalanceFunction(gCentrality, psiMin = -0.5, psiMax,
+    fitbalanceFunction(gCentrality, psiMin , psiMax,
                       ptTriggerMin, ptTriggerMax,
                       ptAssociatedMin, ptAssociatedMax,
                       gHistBalanceFunctionSubtracted,k2pMethod, eventClass);
@@ -606,6 +606,9 @@ void fitbalanceFunction(Int_t gCentrality = 1,
 
 //____________________________________________________________//
 void drawBFPsi2D(const char* lhcPeriod = "LHC11h",
+                const char* gCentralityEstimator = "V0M",
+                Int_t gBit = 128,
+                const char* gEventPlaneEstimator = "VZERO",
                 Int_t gCentrality = 1,
                 Bool_t kShowShuffled = kFALSE, 
                 Bool_t kShowMixed = kFALSE, 
@@ -613,14 +616,19 @@ void drawBFPsi2D(const char* lhcPeriod = "LHC11h",
                 Double_t ptTriggerMin = -1.,
                 Double_t ptTriggerMax = -1.,
                 Double_t ptAssociatedMin = -1.,
-                Double_t ptAssociatedMax = -1.) {
+                Double_t ptAssociatedMax = -1.,
+                Bool_t k2pMethod = kTRUE) {
   //Macro that draws the BF distributions for each centrality bin
   //for reaction plane dependent analysis
   //Author: Panos.Christakoglou@nikhef.nl
   TGaxis::SetMaxDigits(3);
 
   //Get the input file
-  TString filename = lhcPeriod; filename +="/PttFrom";
+  TString filename = lhcPeriod; 
+  filename += "/Centrality"; filename += gCentralityEstimator;
+  filename += "_Bit"; filename += gBit;
+  filename += "_"; filename += gEventPlaneEstimator;
+  filename +="/PttFrom";
   filename += Form("%.1f",ptTriggerMin); filename += "To"; 
   filename += Form("%.1f",ptTriggerMax); filename += "PtaFrom";
   filename += Form("%.1f",ptAssociatedMin); filename += "To"; 
@@ -635,7 +643,9 @@ void drawBFPsi2D(const char* lhcPeriod = "LHC11h",
   filename += Form("%.1f",ptTriggerMin); filename += "To"; 
   filename += Form("%.1f",ptTriggerMax); filename += "PtaFrom";
   filename += Form("%.1f",ptAssociatedMin); filename += "To"; 
-  filename += Form("%.1f",ptAssociatedMax);  filename += ".root";  
+  filename += Form("%.1f",ptAssociatedMax);  
+  if(k2pMethod) filename += "_2pMethod";
+  filename += ".root";  
 
   //Open the file
   TFile *f = TFile::Open(filename.Data());
@@ -736,7 +746,7 @@ void drawBFPsi2D(const char* lhcPeriod = "LHC11h",
   TCanvas *c1 = new TCanvas("c1","Raw balance function 2D",0,0,600,500);
   c1->SetFillColor(10); c1->SetHighLightColor(10);
   c1->SetLeftMargin(0.17); c1->SetTopMargin(0.05);
-  gHistBalanceFunction->SetTitle("Raw balance function");
+  gHistBalanceFunction->SetTitle("");
   gHistBalanceFunction->GetYaxis()->SetTitleOffset(1.4);
   gHistBalanceFunction->GetYaxis()->SetNdivisions(10);
   gHistBalanceFunction->GetXaxis()->SetNdivisions(10);
@@ -750,6 +760,23 @@ void drawBFPsi2D(const char* lhcPeriod = "LHC11h",
   latexInfo1->DrawLatex(0.64,0.76,pttLatex.Data());
   latexInfo1->DrawLatex(0.64,0.70,ptaLatex.Data());
 
+  TString pngName = "BalanceFunction2D."; 
+  pngName += "Centrality";
+  pngName += gCentrality; 
+  if((psiMin == -0.5)&&(psiMax == 0.5)) pngName += "InPlane.Ptt";
+  else if((psiMin == 0.5)&&(psiMax == 1.5)) pngName += "Intermediate.Ptt";
+  else if((psiMin == 1.5)&&(psiMax == 2.5)) pngName += "OutOfPlane.Ptt";
+  else if((psiMin == 2.5)&&(psiMax == 3.5)) pngName += "Rest.PttFrom";
+  else pngName += "All.PttFrom";  
+  pngName += Form("%.1f",ptTriggerMin); pngName += "To"; 
+  pngName += Form("%.1f",ptTriggerMax); pngName += "PtaFrom";
+  pngName += Form("%.1f",ptAssociatedMin); pngName += "To"; 
+  pngName += Form("%.1f",ptAssociatedMax); 
+  if(k2pMethod) pngName += "_2pMethod";
+  pngName += ".png";
+
+  c1->SaveAs(pngName.Data());
+
   if(kShowShuffled) {
     TCanvas *c2 = new TCanvas("c2","Shuffled balance function 2D",100,100,600,500);
     c2->SetFillColor(10); c2->SetHighLightColor(10);

diff --git a/PWGCF/EBYE/macros/drawBalanceFunctionPsi.C b/PWGCF/EBYE/macros/drawBalanceFunctionPsi.C
index 8d033a38328f3326be245468b19b65c69e3d908d..dcff5bc2fbce2651002a2961c8e0e08be2c835de 100644 (file)
--- a/PWGCF/EBYE/macros/drawBalanceFunctionPsi.C
+++ b/PWGCF/EBYE/macros/drawBalanceFunctionPsi.C
@@ -34,7 +34,7 @@ void drawBalanceFunctionPsi(const char* filename = "AnalysisResultsPsi.root",
 
   //Prepare the objects and return them
   TList *listBF = GetListOfObjects(filename,gCentrality,gBit,gCentralityEstimator,0);
-  TList *listBFShuffled = GetListOfObjects(filename,gCentrality,gBit,gCentralityEstimator,1);
+  TList *listBFShuffled = 0x0;//GetListOfObjects(filename,gCentrality,gBit,gCentralityEstimator,1);
   TList *listBFMixed = GetListOfObjects(filename,gCentrality,gBit,gCentralityEstimator,2);
   if(!listBF) {
     Printf("The TList object was not created");
@@ -240,22 +240,24 @@ void draw(TList *listBF, TList *listBFShuffled, TList *listBFMixed,
   AliTHn *hNPShuffled = NULL;
   AliTHn *hPPShuffled = NULL;
   AliTHn *hNNShuffled = NULL;
-  //listBFShuffled->ls();
-  hPShuffled = (AliTHn*) listBFShuffled->FindObject("fHistP_shuffleV0M");
-  hNShuffled = (AliTHn*) listBFShuffled->FindObject("fHistN_shuffleV0M");
-  hPNShuffled = (AliTHn*) listBFShuffled->FindObject("fHistPN_shuffleV0M");
-  hNPShuffled = (AliTHn*) listBFShuffled->FindObject("fHistNP_shuffleV0M");
-  hPPShuffled = (AliTHn*) listBFShuffled->FindObject("fHistPP_shuffleV0M");
-  hNNShuffled = (AliTHn*) listBFShuffled->FindObject("fHistNN_shuffleV0M");
-
-  AliBalancePsi *bShuffled = new AliBalancePsi();
-  bShuffled->SetEventClass(eventClass);
-  bShuffled->SetHistNp(hPShuffled);
-  bShuffled->SetHistNn(hNShuffled);
-  bShuffled->SetHistNpn(hPNShuffled);
-  bShuffled->SetHistNnp(hNPShuffled);
-  bShuffled->SetHistNpp(hPPShuffled);
-  bShuffled->SetHistNnn(hNNShuffled);
+  if(listBFShuffled) {
+    //listBFShuffled->ls();
+    hPShuffled = (AliTHn*) listBFShuffled->FindObject("fHistP_shuffleV0M");
+    hNShuffled = (AliTHn*) listBFShuffled->FindObject("fHistN_shuffleV0M");
+    hPNShuffled = (AliTHn*) listBFShuffled->FindObject("fHistPN_shuffleV0M");
+    hNPShuffled = (AliTHn*) listBFShuffled->FindObject("fHistNP_shuffleV0M");
+    hPPShuffled = (AliTHn*) listBFShuffled->FindObject("fHistPP_shuffleV0M");
+    hNNShuffled = (AliTHn*) listBFShuffled->FindObject("fHistNN_shuffleV0M");
+    
+    AliBalancePsi *bShuffled = new AliBalancePsi();
+    bShuffled->SetEventClass(eventClass);
+    bShuffled->SetHistNp(hPShuffled);
+    bShuffled->SetHistNn(hNShuffled);
+    bShuffled->SetHistNpn(hPNShuffled);
+    bShuffled->SetHistNnp(hNPShuffled);
+    bShuffled->SetHistNpp(hPPShuffled);
+    bShuffled->SetHistNnn(hNNShuffled);
+  }
 
   //balance function mixing
   AliTHn *hPMixed = NULL;
@@ -348,31 +350,31 @@ void draw(TList *listBF, TList *listBFShuffled, TList *listBFMixed,
   gHistBalanceFunction->SetName("gHistBalanceFunction");
   
   //Shuffled balance function
-  if(k2pMethod){ 
-    if(bMixed)
-      gHistBalanceFunctionShuffled = bShuffled->GetBalanceFunctionHistogram2pMethod(0,gDeltaEtaDeltaPhi,psiMin,psiMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax,bMixed);
-    else{
-      cerr<<"SHUFFLE: NO MIXED BF BUT REQUESTED CORRECTING WITH IT! --> FAIL"<<endl;
-      return;
-    }
-  }
-  else if(k2pMethod2D){ 
-    if(bMixed){
-      if(gDeltaEtaDeltaPhi==1) //Delta eta
-       gHistBalanceFunctionShuffled = bShuffled->GetBalanceFunction1DFrom2D2pMethod(0,psiMin,psiMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax,bMixed);
-      else //Delta phi
-       gHistBalanceFunctionShuffled = bShuffled->GetBalanceFunction1DFrom2D2pMethod(1,psiMin,psiMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax,bMixed);
-    }
-    else{
-      cerr<<"SHUFFLE: NO MIXED BF BUT REQUESTED CORRECTING WITH IT! --> FAIL"<<endl;
-      return;
-    }
-  }
-  else
-    gHistBalanceFunctionShuffled = bShuffled->GetBalanceFunctionHistogram(0,gDeltaEtaDeltaPhi,psiMin,psiMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax);
-  gHistBalanceFunctionShuffled->SetMarkerStyle(24);
-  gHistBalanceFunctionShuffled->SetName("gHistBalanceFunctionShuffled");
-
+  //if(k2pMethod){ 
+  //if(bMixed)
+      //gHistBalanceFunctionShuffled = bShuffled->GetBalanceFunctionHistogram2pMethod(0,gDeltaEtaDeltaPhi,psiMin,psiMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax,bMixed);
+  //else{
+  //cerr<<"SHUFFLE: NO MIXED BF BUT REQUESTED CORRECTING WITH IT! --> FAIL"<<endl;
+  //return;
+  //}
+  //}
+  //else if(k2pMethod2D){ 
+  //if(bMixed){
+  //  if(gDeltaEtaDeltaPhi==1) //Delta eta
+  //gHistBalanceFunctionShuffled = bShuffled->GetBalanceFunction1DFrom2D2pMethod(0,psiMin,psiMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax,bMixed);
+  //  else //Delta phi
+  //gHistBalanceFunctionShuffled = bShuffled->GetBalanceFunction1DFrom2D2pMethod(1,psiMin,psiMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax,bMixed);
+  //}
+  //else{
+  //  cerr<<"SHUFFLE: NO MIXED BF BUT REQUESTED CORRECTING WITH IT! --> FAIL"<<endl;
+  //  return;
+  //}
+  //}
+  //else
+  //gHistBalanceFunctionShuffled = bShuffled->GetBalanceFunctionHistogram(0,gDeltaEtaDeltaPhi,psiMin,psiMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax);
+  //gHistBalanceFunctionShuffled->SetMarkerStyle(24);
+  //gHistBalanceFunctionShuffled->SetName("gHistBalanceFunctionShuffled");
+  
   //Mixed balance function
   if(k2pMethod){ 
     if(bMixed)
@@ -414,7 +416,7 @@ void draw(TList *listBF, TList *listBFShuffled, TList *listBFMixed,
   c1->SetHighLightColor(10);
   c1->SetLeftMargin(0.15);
   gHistBalanceFunction->DrawCopy("E");
-  gHistBalanceFunctionShuffled->DrawCopy("ESAME");
+  //gHistBalanceFunctionShuffled->DrawCopy("ESAME");
   gHistBalanceFunctionMixed->DrawCopy("ESAME");
   
   legend = new TLegend(0.18,0.62,0.45,0.82,"","brNDC");
@@ -426,7 +428,7 @@ void draw(TList *listBF, TList *listBFShuffled, TList *listBFMixed,
   legend->SetMargin(0.25); 
   legend->SetShadowColor(10);
   legend->AddEntry(gHistBalanceFunction,"Data","lp");
-  legend->AddEntry(gHistBalanceFunctionShuffled,"Shuffled data","lp");
+  //legend->AddEntry(gHistBalanceFunctionShuffled,"Shuffled data","lp");
   legend->AddEntry(gHistBalanceFunctionMixed,"Mixed data","lp");
   legend->Draw();
   
@@ -463,7 +465,7 @@ void draw(TList *listBF, TList *listBFShuffled, TList *listBFMixed,
   c1->SaveAs(pngName.Data());
   
   GetWeightedMean(gHistBalanceFunction);
-  GetWeightedMean(gHistBalanceFunctionShuffled);
+  //GetWeightedMean(gHistBalanceFunctionShuffled);
   
   TString meanLatex, rmsLatex, skewnessLatex, kurtosisLatex;
   meanLatex = "#mu = "; 
@@ -538,7 +540,7 @@ void draw(TList *listBF, TList *listBFShuffled, TList *listBFMixed,
   TFile *fOutput = new TFile(newFileName.Data(),"recreate");
   fOutput->cd();
   gHistBalanceFunction->Write();
-  gHistBalanceFunctionShuffled->Write();
+  //gHistBalanceFunctionShuffled->Write();
   gHistBalanceFunctionMixed->Write();
   gHistBalanceFunctionSubtracted->Write();
   fOutput->Close();