optimize range for LeeYangZeroes to find zero

diff --git a/PWG2/FLOW/AliFlowCommon/AliFlowAnalysisWithLeeYangZeros.cxx b/PWG2/FLOW/AliFlowCommon/AliFlowAnalysisWithLeeYangZeros.cxx
index 14ecb243339c8f417b61cc71f3f460a13d8c5a6c..e929d0fff559e397063480e304b07d15abb8133b 100644 (file)
--- a/PWG2/FLOW/AliFlowCommon/AliFlowAnalysisWithLeeYangZeros.cxx
+++ b/PWG2/FLOW/AliFlowCommon/AliFlowAnalysisWithLeeYangZeros.cxx
@@ -613,7 +613,9 @@ void AliFlowAnalysisWithLeeYangZeros::GetOutputHistograms(TList *outputListHisto
        else { cout<<"r0 is not found! Leaving LYZ analysis."<<endl; return kFALSE; }
 
        //for estimating systematic error resulting from d0
-       Double_t dBinsize = (AliFlowLYZConstants::fgMax)/(AliFlowLYZConstants::kNbins);
+       Double_t dBinsize =0.;
+       if (fUseSum){ dBinsize = (AliFlowLYZConstants::fgMaxSUM)/(AliFlowLYZConstants::kNbins);}
+       else { dBinsize = (AliFlowLYZConstants::fgMaxPROD)/(AliFlowLYZConstants::kNbins);}
        Double_t dVplus = -1.;
        Double_t dVmin  = -1.;
        if (dR0+dBinsize!=0.) {dVplus = dJ01/(dR0+dBinsize);}
@@ -645,14 +647,17 @@ void AliFlowAnalysisWithLeeYangZeros::GetOutputHistograms(TList *outputListHisto
     //get average value of fVtheta = fV
     dV /=iNtheta;
     if (!fUseSum) { if (dMultRP!=0.){dV /=dMultRP;}} //scale with multiplicity for PRODUCT
-
+    
     //sigma2 and chi 
     Double_t  dSigma2 = 0;
     Double_t  dChi= 0;
     if (fEventNumber!=0) {
       *fQsum /= fEventNumber;
+      //cout<<"fQsum is "<<fQsum->X()<<" "<<fQsum->Y()<<endl; 
       fQ2sum /= fEventNumber;
+      //cout<<"fQ2sum is "<<fQ2sum<<endl; 
       dSigma2 = fQ2sum - TMath::Power(fQsum->X(),2.) - TMath::Power(fQsum->Y(),2.) - TMath::Power(dV,2.);  //BP eq. 62
+      //cout<<"dSigma2 is "<<dSigma2<<endl; 
       if (dSigma2>0) dChi = dV/TMath::Sqrt(dSigma2);
       else dChi = -1.;
       fCommonHistsRes->FillChiRP(dChi);

diff --git a/PWG2/FLOW/AliFlowCommon/AliFlowLYZConstants.cxx b/PWG2/FLOW/AliFlowCommon/AliFlowLYZConstants.cxx
index 1dbd1c3dfcfafecac168e5e61cb835a894bbbbda..7b819753df25024b882ab6437cee1390d1493fee 100644 (file)
--- a/PWG2/FLOW/AliFlowCommon/AliFlowLYZConstants.cxx
+++ b/PWG2/FLOW/AliFlowCommon/AliFlowLYZConstants.cxx
@@ -25,8 +25,8 @@ $Log$
 // Author: Naomi van der Kolk (kolk@nikhef.nl)
 
   
-Double_t AliFlowLYZConstants::fgMin       =  0. ;
-Double_t AliFlowLYZConstants::fgMax       =  120.;
+Double_t AliFlowLYZConstants::fgMaxSUM       =  120.; 
+Double_t AliFlowLYZConstants::fgMaxPROD      =  1.; 
 
   
 

diff --git a/PWG2/FLOW/AliFlowCommon/AliFlowLYZConstants.h b/PWG2/FLOW/AliFlowCommon/AliFlowLYZConstants.h
index c75ca4b71ec6357c75cf3d5b4e56340595b5ae23..b2cb9686e9756b988011ddbb040613c2fd421f11 100644 (file)
--- a/PWG2/FLOW/AliFlowCommon/AliFlowLYZConstants.h
+++ b/PWG2/FLOW/AliFlowCommon/AliFlowLYZConstants.h
@@ -18,13 +18,13 @@ namespace AliFlowLYZConstants {
  // Enumerators
   enum {
     kTheta       = 5,     // number of reference angles theta
-    kNbins       = 8000   // number of bins in fHistGtheta (AliFlowLYZHist1)
+    kNbins       = 1200   // number of bins in fHistGtheta (AliFlowLYZHist1)
   };
  
 
  // Histograms limits
-  extern Double_t  fgMin ;   // lower limit for fHistGtheta (AliFlowLYZHist1)
-  extern Double_t  fgMax ;   // upper limit for fHistGtheta (AliFlowLYZHist1)
+  extern Double_t  fgMaxSUM ;   // upper limit for fHistGtheta (AliFlowLYZHist1)
+  extern Double_t  fgMaxPROD ;   // upper limit for fHistGtheta (AliFlowLYZHist1)
    
  
 }

diff --git a/PWG2/FLOW/AliFlowCommon/AliFlowLYZHist1.cxx b/PWG2/FLOW/AliFlowCommon/AliFlowLYZHist1.cxx
index b39ec0998d30d21bea62cafda8e862e2585cc960..9077e935a6606a53fdf0fc34ea9a9c0a755a97a6 100644 (file)
--- a/PWG2/FLOW/AliFlowCommon/AliFlowLYZHist1.cxx
+++ b/PWG2/FLOW/AliFlowCommon/AliFlowLYZHist1.cxx
@@ -56,9 +56,11 @@ ClassImp(AliFlowLYZHist1)
 {
 
   //constructor creating histograms 
-  Int_t iNbins = AliFlowLYZConstants::kNbins;
-  Double_t dMin = AliFlowLYZConstants::fgMin;
-  Double_t dMax = AliFlowLYZConstants::fgMax;
+  Int_t iNbins      = AliFlowLYZConstants::kNbins;
+  Double_t dMaxSUM  = AliFlowLYZConstants::fgMaxSUM;
+  Double_t dMaxPROD = AliFlowLYZConstants::fgMaxPROD;
+  Double_t dMin     = 0.;
+
   TString name, addlast;
  
   if (useSum) { addlast = "LYZSUM"; }
@@ -68,7 +70,8 @@ ClassImp(AliFlowLYZHist1)
   name = "First_Flow_Gtheta";
   name +=theta;
   name +=addlast;
-  fHistGtheta = new TH1D(name.Data(),name.Data(),iNbins,dMin,dMax);  
+  if (useSum) { fHistGtheta = new TH1D(name.Data(),name.Data(),iNbins,dMin,dMaxSUM); }
+  else { fHistGtheta = new TH1D(name.Data(),name.Data(),iNbins,dMin,dMaxPROD); } 
   fHistGtheta->SetXTitle("r");
   fHistGtheta->SetYTitle("|G^{#theta}(ir)|^{2}");
   
@@ -76,7 +79,8 @@ ClassImp(AliFlowLYZHist1)
   name = "First_FlowPro_ReGtheta";
   name +=theta;
   name +=addlast;
-  fHistProReGtheta = new TProfile(name.Data(),name.Data(),iNbins,dMin,dMax);
+  if (useSum) { fHistProReGtheta = new TProfile(name.Data(),name.Data(),iNbins,dMin,dMaxSUM); }
+  else { fHistProReGtheta = new TProfile(name.Data(),name.Data(),iNbins,dMin,dMaxPROD); }
   fHistProReGtheta->SetXTitle("r");
   fHistProReGtheta->SetYTitle("Re G^{#theta}(ir)");
   
@@ -84,7 +88,8 @@ ClassImp(AliFlowLYZHist1)
   name = "First_FlowPro_ImGtheta";
   name +=theta;
   name +=addlast;
-  fHistProImGtheta = new TProfile(name.Data(),name.Data(),iNbins,dMin,dMax);
+  if (useSum) { fHistProImGtheta = new TProfile(name.Data(),name.Data(),iNbins,dMin,dMaxSUM); }
+  else { fHistProImGtheta = new TProfile(name.Data(),name.Data(),iNbins,dMin,dMaxPROD); }
   fHistProImGtheta->SetXTitle("r");
   fHistProImGtheta->SetYTitle("Im G^{#theta}(ir)");