Changes in drawing macros:
authormiweber <m.weber@cern.ch>
Fri, 17 Oct 2014 18:25:14 +0000 (20:25 +0200)
committermiweber <m.weber@cern.ch>
Fri, 17 Oct 2014 18:27:16 +0000 (20:27 +0200)
- remove range restriction in Delta Eta when projecting
- automatic calculation of range for medium and high momentum = 3*sigma of Gauss fit (only if not above safety margin, to be checked)

PWGCF/EBYE/macros/drawBalanceFunction2DPsi.C
PWGCF/EBYE/macros/drawCorrelationFunctionPsi.C

index db7024a..997dbe4 100644 (file)
@@ -829,7 +829,7 @@ void drawBFPsi2D(const char* lhcPeriod = "LHC11h",
   gHistBalanceFunction->SetTitle("");
   gHistBalanceFunction->GetYaxis()->SetTitleOffset(1.4);
   gHistBalanceFunction->GetYaxis()->SetNdivisions(10);
-  gHistBalanceFunction->GetXaxis()->SetRangeUser(-1.4,1.4); 
+  //gHistBalanceFunction->GetXaxis()->SetRangeUser(-1.4,1.4); 
   gHistBalanceFunction->GetXaxis()->SetNdivisions(10);
   gHistBalanceFunction->GetYaxis()->SetTitle("#Delta #varphi (rad)");
   gHistBalanceFunction->DrawCopy("lego2");
@@ -1099,14 +1099,26 @@ void drawProjections(TH2D *gHistBalanceFunction2D = 0x0,
     // need to restrict to near side in case of dphi
     if(!kProjectInEta) gHistBalanceFunctionSubtracted->GetXaxis()->SetRangeUser(-TMath::Pi()/2,TMath::Pi()/2);
 
-    if(bReduceRangeForMoments){
+   if(bReduceRangeForMoments){
 
+      // safety check:
       // default (for 2<pT,assoc<3<pT,trig<4)
-      Double_t rangeReduced = 1.0;
-
+      Double_t rangeReduced = 1.2;
       //for 3<pT,assoc<8<pT,trig<15
       if(ptTriggerMax>10){
-       rangeReduced = 0.3;
+       rangeReduced = 0.35;
+      }
+
+      // new define range by fit!
+      gHistBalanceFunctionSubtracted->Fit("gaus","","");
+      Double_t sigmaGaus = gHistBalanceFunctionSubtracted->GetFunction("gaus")->GetParameter(2);
+
+      // if safety check OK, set rangeReduced to 3 sigma of gauss fit
+      if(3*sigmaGaus > rangeReduced){
+       cout<<"RANGE REDUCE ERROR: "<<3*sigmaGaus<<" > "<<rangeReduced<<endl;
+      }
+      else{
+       rangeReduced = 3*sigmaGaus;
       }
 
       cout<<"Use reduced range = "<<rangeReduced<<endl;
@@ -1432,6 +1444,16 @@ void drawBFPsi2DFromCorrelationFunctions(const char* lhcPeriod = "LHC10h",
   TH2D *gHistNN = dynamic_cast<TH2D *>(f->Get("gHistNNCorrelationFunctions"));
   if(!gHistNN) return;
 
+  // in order to get unzoomed (in older versions used smaller user ranger)
+  Int_t binMinX = gHistPN->GetXaxis()->FindBin(gHistPN->GetXaxis()->GetXmin()+0.00001);
+  Int_t binMaxX = gHistPN->GetXaxis()->FindBin(gHistPN->GetXaxis()->GetXmax()-0.00001);
+  Int_t binMinY = gHistPN->GetYaxis()->FindBin(gHistPN->GetYaxis()->GetXmin()+0.00001);
+  Int_t binMaxY = gHistPN->GetYaxis()->FindBin(gHistPN->GetYaxis()->GetXmax()-0.00001);
+  gHistPN->GetXaxis()->SetRange(binMinX,binMaxX); gHistPN->GetYaxis()->SetRange(binMinY,binMaxY);
+  gHistNP->GetXaxis()->SetRange(binMinX,binMaxX); gHistNP->GetYaxis()->SetRange(binMinY,binMaxY);
+  gHistPP->GetXaxis()->SetRange(binMinX,binMaxX); gHistPP->GetYaxis()->SetRange(binMinY,binMaxY);
+  gHistNN->GetXaxis()->SetRange(binMinX,binMaxX); gHistNN->GetYaxis()->SetRange(binMinY,binMaxY);
+
   gHistPN->Sumw2();
   gHistPP->Sumw2();
   gHistPN->Add(gHistPP,-1);
@@ -1455,7 +1477,7 @@ void drawBFPsi2DFromCorrelationFunctions(const char* lhcPeriod = "LHC10h",
   gHistBalanceFunction2D->GetZaxis()->SetNdivisions(10);
   gHistBalanceFunction2D->GetYaxis()->SetTitleOffset(1.4);
   gHistBalanceFunction2D->GetYaxis()->SetNdivisions(10);
-  gHistBalanceFunction2D->GetXaxis()->SetRangeUser(-1.4,1.4); 
+  //gHistBalanceFunction2D->GetXaxis()->SetRangeUser(-1.4,1.4); 
   gHistBalanceFunction2D->GetXaxis()->SetNdivisions(10);
   gHistBalanceFunction2D->DrawCopy("lego2");
   gPad->SetTheta(30); // default is 30
@@ -1570,8 +1592,6 @@ void mergeBFPsi2D(TString momDirectory = "./",
   TGaxis::SetMaxDigits(3);
   gStyle->SetPalette(55,0);
 
-  cout<<"REDUCE"<<bReduceRangeForMoments<<"  "<<endl;
-
   const Int_t nMaxDirectories = 4; // maximum number of directories to merge (set to 4 for now)
   TString sDirectory[nMaxDirectories];
   Int_t nDirectories = nMaxDirectories;
index ebcea28..182f7bc 100644 (file)
@@ -460,7 +460,7 @@ void draw(TList *list, TList *listBFShuffled, TList *listBFMixed,
       gHistPP[2]->Scale(1./gHistPP[2]->GetBinContent(gHistPP[2]->FindBin(0,0)));
     }
   }
-  
+
   //(+-)
   if(eventClass == "Centrality"){
     histoTitle = "(+-) | Centrality: ";
@@ -586,7 +586,8 @@ void draw(TList *list, TList *listBFShuffled, TList *listBFMixed,
 
     //Correlation function (+-)
     gHistPN[3] = b->GetCorrelationFunction("PN",psiMin,psiMax,vertexZMin,vertexZMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax,bMixed,normToTrig,normalizationRangePhi);
-    gHistPN[3]->GetXaxis()->SetRangeUser(-1.5,1.5);
+    //gHistPN[3]->GetXaxis()->SetRangeUser(-1.5,1.5);
+
     if(normToTrig)
       gHistPN[3]->GetZaxis()->SetTitle("#frac{1}{N_{trig}}#frac{d^{2}N_{assoc}}{d#Delta#eta#Delta#varphi} (rad^{-1})");
     else
@@ -653,7 +654,7 @@ void draw(TList *list, TList *listBFShuffled, TList *listBFMixed,
     //Correlation function (+-)
     gHistPN[3] = dynamic_cast<TH2D *>(gHistPN[0]->Clone());
     gHistPN[3]->Divide(gHistPN[2]);
-    gHistPN[3]->GetXaxis()->SetRangeUser(-1.5,1.5);
+    //gHistPN[3]->GetXaxis()->SetRangeUser(-1.5,1.5);
     if(normToTrig)
       gHistPN[3]->GetZaxis()->SetTitle("#frac{1}{N_{trig}}#frac{d^{2}N_{assoc}}{d#Delta#eta#Delta#varphi} (rad^{-1})");
     else
@@ -800,7 +801,7 @@ void draw(TList *list, TList *listBFShuffled, TList *listBFMixed,
 
     //Correlation function (-+)
     gHistNP[3] = b->GetCorrelationFunction("NP",psiMin,psiMax,vertexZMin,vertexZMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax,bMixed,normToTrig,normalizationRangePhi);
-    gHistNP[3]->GetXaxis()->SetRangeUser(-1.5,1.5);
+    //gHistNP[3]->GetXaxis()->SetRangeUser(-1.5,1.5);
     if(normToTrig)
       gHistNP[3]->GetZaxis()->SetTitle("#frac{1}{N_{trig}}#frac{d^{2}N_{assoc}}{d#Delta#eta#Delta#varphi} (rad^{-1})");
     else
@@ -867,7 +868,7 @@ void draw(TList *list, TList *listBFShuffled, TList *listBFMixed,
     //Correlation function (-+)
     gHistNP[3] = dynamic_cast<TH2D *>(gHistNP[0]->Clone());
     gHistNP[3]->Divide(gHistNP[2]);
-    gHistNP[3]->GetXaxis()->SetRangeUser(-1.5,1.5);
+    //gHistNP[3]->GetXaxis()->SetRangeUser(-1.5,1.5);
     if(normToTrig)
       gHistNP[3]->GetZaxis()->SetTitle("#frac{1}{N_{trig}}#frac{d^{2}N_{assoc}}{d#Delta#eta#Delta#varphi} (rad^{-1})");
     else
@@ -1015,7 +1016,7 @@ void draw(TList *list, TList *listBFShuffled, TList *listBFMixed,
 
     //Correlation function (++)
     gHistPP[3] = b->GetCorrelationFunction("PP",psiMin,psiMax,vertexZMin,vertexZMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax,bMixed,normToTrig,normalizationRangePhi);
-    gHistPP[3]->GetXaxis()->SetRangeUser(-1.5,1.5);
+    //gHistPP[3]->GetXaxis()->SetRangeUser(-1.5,1.5);
     if(normToTrig)
       gHistPP[3]->GetZaxis()->SetTitle("#frac{1}{N_{trig}}#frac{d^{2}N_{assoc}}{d#Delta#eta#Delta#varphi} (rad^{-1})");
     else
@@ -1081,7 +1082,7 @@ void draw(TList *list, TList *listBFShuffled, TList *listBFMixed,
     //Correlation function (++)
     gHistPP[3] = dynamic_cast<TH2D *>(gHistPP[0]->Clone());
     gHistPP[3]->Divide(gHistPP[2]);
-    gHistPP[3]->GetXaxis()->SetRangeUser(-1.5,1.5);
+    //gHistPP[3]->GetXaxis()->SetRangeUser(-1.5,1.5);
     if(normToTrig)
       gHistPP[3]->GetZaxis()->SetTitle("#frac{1}{N_{trig}}#frac{d^{2}N_{assoc}}{d#Delta#eta#Delta#varphi} (rad^{-1})");
     else
@@ -1228,7 +1229,7 @@ void draw(TList *list, TList *listBFShuffled, TList *listBFMixed,
 
     //Correlation function (--)
     gHistNN[3] = b->GetCorrelationFunction("NN",psiMin,psiMax,vertexZMin,vertexZMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax,bMixed,normToTrig,normalizationRangePhi);
-    gHistNN[3]->GetXaxis()->SetRangeUser(-1.5,1.5);
+    //gHistNN[3]->GetXaxis()->SetRangeUser(-1.5,1.5);
     if(normToTrig)
       gHistNN[3]->GetZaxis()->SetTitle("#frac{1}{N_{trig}}#frac{d^{2}N_{assoc}}{d#Delta#eta#Delta#varphi} (rad^{-1})");
     else
@@ -1294,7 +1295,7 @@ void draw(TList *list, TList *listBFShuffled, TList *listBFMixed,
     //Correlation function (--)
     gHistNN[3] = dynamic_cast<TH2D *>(gHistNN[0]->Clone());
     gHistNN[3]->Divide(gHistNN[2]);
-    gHistNN[3]->GetXaxis()->SetRangeUser(-1.5,1.5);
+    //gHistNN[3]->GetXaxis()->SetRangeUser(-1.5,1.5);
     if(normToTrig)
       gHistNN[3]->GetZaxis()->SetTitle("#frac{1}{N_{trig}}#frac{d^{2}N_{assoc}}{d#Delta#eta#Delta#varphi} (rad^{-1})");
     else
@@ -1501,7 +1502,7 @@ void drawCorrelationFunctions(const char* lhcPeriod = "LHC10h",
   TH2D *gHistPN = dynamic_cast<TH2D *>(f->Get("gHistPNCorrelationFunctions"));
   gHistPN->SetStats(kFALSE);
   gHistPN->SetTitle("C_{+-}(#Delta#eta,#Delta#varphi)");
-  gHistPN->GetXaxis()->SetRangeUser(-1.4,1.4);
+  //gHistPN->GetXaxis()->SetRangeUser(-1.4,1.4);
   gHistPN->GetXaxis()->CenterTitle();
   gHistPN->GetXaxis()->SetTitleOffset(1.2);
   gHistPN->GetYaxis()->CenterTitle();
@@ -1545,7 +1546,7 @@ void drawCorrelationFunctions(const char* lhcPeriod = "LHC10h",
   TH2D *gHistNP = dynamic_cast<TH2D *>(f->Get("gHistNPCorrelationFunctions"));
   gHistNP->SetStats(kFALSE);
   gHistNP->SetTitle("C_{-+}(#Delta#eta,#Delta#varphi)");
-  gHistNP->GetXaxis()->SetRangeUser(-1.4,1);
+  //gHistNP->GetXaxis()->SetRangeUser(-1.4,1.4);
   gHistNP->GetXaxis()->CenterTitle();
   gHistNP->GetXaxis()->SetTitleOffset(1.2);
   gHistNP->GetYaxis()->CenterTitle();
@@ -1590,7 +1591,7 @@ void drawCorrelationFunctions(const char* lhcPeriod = "LHC10h",
   TH2D *gHistPP = dynamic_cast<TH2D *>(f->Get("gHistPPCorrelationFunctions"));
   gHistPP->SetStats(kFALSE);
   gHistPP->SetTitle("C_{++}(#Delta#eta,#Delta#varphi)");
-  gHistPP->GetXaxis()->SetRangeUser(-1.4,1.4);
+  //gHistPP->GetXaxis()->SetRangeUser(-1.4,1.4);
   gHistPP->GetXaxis()->CenterTitle();
   gHistPP->GetXaxis()->SetTitleOffset(1.2);
   gHistPP->GetYaxis()->CenterTitle();
@@ -1635,7 +1636,7 @@ void drawCorrelationFunctions(const char* lhcPeriod = "LHC10h",
   TH2D *gHistNN = dynamic_cast<TH2D *>(f->Get("gHistNNCorrelationFunctions"));
   gHistNN->SetStats(kFALSE);
   gHistNN->SetTitle("C_{--}(#Delta#eta,#Delta#varphi)");
-  gHistNN->GetXaxis()->SetRangeUser(-1.4,1.4);
+  //gHistNN->GetXaxis()->SetRangeUser(-1.4,1.4);
   gHistNN->GetXaxis()->CenterTitle();
   gHistNN->GetXaxis()->SetTitleOffset(1.2);
   gHistNN->GetYaxis()->CenterTitle();
@@ -1782,7 +1783,7 @@ void drawProjections(Bool_t kProjectInEta = kFALSE,
   TH2D *gHistPN = dynamic_cast<TH2D *>(f->Get("gHistPNCorrelationFunctions"));
   gHistPN->SetStats(kFALSE);
   gHistPN->SetTitle("C_{+-}(#Delta#eta,#Delta#varphi)");
-  gHistPN->GetXaxis()->SetRangeUser(-1.4,1.4);
+  //gHistPN->GetXaxis()->SetRangeUser(-1.4,1.4);
   gHistPN->GetXaxis()->CenterTitle();
   gHistPN->GetXaxis()->SetTitleOffset(1.2);
   gHistPN->GetYaxis()->CenterTitle();
@@ -1816,7 +1817,7 @@ void drawProjections(Bool_t kProjectInEta = kFALSE,
       gHistPNprojection->SetBinContent(iBinX,sum);
       gHistPNprojection->SetBinError(iBinX,gError);
     }
-    gHistPNprojection->GetXaxis()->SetRangeUser(-1.4,1.4);
+    //gHistPNprojection->GetXaxis()->SetRangeUser(-1.4,1.4);
     if(kUseZYAM) 
       gHistPNprojection->GetYaxis()->SetTitle("#frac{1}{N_{trig}}#frac{dN_{assoc}}{#Delta#eta} - b_{ZYAM}");
     else 
@@ -1882,7 +1883,7 @@ void drawProjections(Bool_t kProjectInEta = kFALSE,
   TH2D *gHistNP = dynamic_cast<TH2D *>(f->Get("gHistNPCorrelationFunctions"));
   gHistNP->SetStats(kFALSE);
   gHistNP->SetTitle("C_{-+}(#Delta#eta,#Delta#varphi)");
-  gHistNP->GetXaxis()->SetRangeUser(-1.4,1.4);
+  //gHistNP->GetXaxis()->SetRangeUser(-1.4,1.4);
   gHistNP->GetXaxis()->CenterTitle();
   gHistNP->GetXaxis()->SetTitleOffset(1.2);
   gHistNP->GetYaxis()->CenterTitle();
@@ -1915,7 +1916,7 @@ void drawProjections(Bool_t kProjectInEta = kFALSE,
       gHistNPprojection->SetBinContent(iBinX,sum);
       gHistNPprojection->SetBinError(iBinX,gError);
     }
-    gHistNPprojection->GetXaxis()->SetRangeUser(-1.4,1.4);
+    //gHistNPprojection->GetXaxis()->SetRangeUser(-1.4,1.4);
     if(kUseZYAM) 
       gHistNPprojection->GetYaxis()->SetTitle("#frac{1}{N_{trig}}#frac{dN_{assoc}}{#Delta#eta} - b_{ZYAM}");
     else
@@ -1979,7 +1980,7 @@ void drawProjections(Bool_t kProjectInEta = kFALSE,
   TH2D *gHistPP = dynamic_cast<TH2D *>(f->Get("gHistPPCorrelationFunctions"));
   gHistPP->SetStats(kFALSE);
   gHistPP->SetTitle("C_{++}(#Delta#eta,#Delta#varphi)");
-  gHistPP->GetXaxis()->SetRangeUser(-1.4,1.4);
+  //gHistPP->GetXaxis()->SetRangeUser(-1.4,1.4);
   gHistPP->GetXaxis()->CenterTitle();
   gHistPP->GetXaxis()->SetTitleOffset(1.2);
   gHistPP->GetYaxis()->CenterTitle();
@@ -2012,7 +2013,7 @@ void drawProjections(Bool_t kProjectInEta = kFALSE,
       gHistPPprojection->SetBinContent(iBinX,sum);
       gHistPPprojection->SetBinError(iBinX,gError);
     }
-    gHistPPprojection->GetXaxis()->SetRangeUser(-1.4,1.4);
+    //gHistPPprojection->GetXaxis()->SetRangeUser(-1.4,1.4);
     if(kUseZYAM) 
       gHistPPprojection->GetYaxis()->SetTitle("#frac{1}{N_{trig}}#frac{dN_{assoc}}{#Delta#eta} - b_{ZYAM}");
     else
@@ -2076,7 +2077,7 @@ void drawProjections(Bool_t kProjectInEta = kFALSE,
   TH2D *gHistNN = dynamic_cast<TH2D *>(f->Get("gHistNNCorrelationFunctions"));
   gHistNN->SetStats(kFALSE);
   gHistNN->SetTitle("C_{--}(#Delta#eta,#Delta#varphi)");
-  gHistNN->GetXaxis()->SetRangeUser(-1.4,1.4);
+  //gHistNN->GetXaxis()->SetRangeUser(-1.4,1.4);
   gHistNN->GetXaxis()->CenterTitle();
   gHistNN->GetXaxis()->SetTitleOffset(1.2);
   gHistNN->GetYaxis()->CenterTitle();
@@ -2109,7 +2110,7 @@ void drawProjections(Bool_t kProjectInEta = kFALSE,
       gHistNNprojection->SetBinContent(iBinX,sum);
       gHistNNprojection->SetBinError(iBinX,gError);
     }
-    gHistNNprojection->GetXaxis()->SetRangeUser(-1.4,1.4);
+    //gHistNNprojection->GetXaxis()->SetRangeUser(-1.4,1.4);
     if(kUseZYAM) 
       gHistNNprojection->GetYaxis()->SetTitle("#frac{1}{N_{trig}}#frac{dN_{assoc}}{#Delta#eta} - b_{ZYAM}");
     else