Coverity and bug fix: prabhat n sjena
authorsjena <sjena@cern.ch>
Wed, 10 Dec 2014 14:47:01 +0000 (15:47 +0100)
committersjena <sjena@cern.ch>
Wed, 10 Dec 2014 14:47:01 +0000 (15:47 +0100)
PWGCF/Correlations/DPhi/AliAnalysisTaskpxpx.cxx
PWGCF/Correlations/DPhi/AliAnalysisTaskpxpy.cxx
PWGCF/Correlations/DPhi/AliAnalysisTaskpypy.cxx
PWGCF/Correlations/DPhi/AliAnalysisTaskpzpz.cxx
PWGCF/EBYE/PIDFluctuation/macros/AddAliEbyENetChargeFluctuationTask.C
PWGCF/EBYE/PIDFluctuation/task/AliEbyENetChargeFluctuationTask.cxx
PWGCF/EBYE/PIDFluctuation/task/AliEbyENetChargeFluctuationTask.h
PWGCF/EBYE/PIDFluctuation/task/AliEbyEParticleRatioFluctuationTask.cxx

index 9484ecd..1dc1087 100644 (file)
@@ -1063,8 +1063,8 @@ void  AliAnalysisTaskpxpx::UserExec(Option_t */*option*/)
   float trkCentr = -999.;
   float spdCentr = -999.;
   
-  float vertexX  = -999;
-  float vertexY  = -999;
+  //float vertexX  = -999;
+  //float vertexY  = -999;
   float vertexZ  = -999;
   //float vertexXY = -999;
   //float dcaZ     = -999;
@@ -1128,8 +1128,8 @@ void  AliAnalysisTaskpxpx::UserExec(Option_t */*option*/)
            {
              if(fCov[5] != 0)
                {
-                 vertexX = vertex->GetX();
-                 vertexY = vertex->GetY();
+                 //vertexX = vertex->GetX();
+                 //vertexY = vertex->GetY();
                  vertexZ = vertex->GetZ();
                  
                  if(TMath::Abs(vertexZ) > 10)
@@ -1190,12 +1190,15 @@ void  AliAnalysisTaskpxpx::UserExec(Option_t */*option*/)
          charge = int(q);
          phi    = t->Phi();
          pt     = t->Pt(); 
-         px     = t->Px();
-         py     = t->Py();
+         //px     = t->Px();
+         //py     = t->Py();
          pz     = t->Pz();
          eta    = t->Eta();
          dedx   = t->GetTPCsignal();
                  
+         px     = pt*cos(phi);
+         py     = pt*sin(phi);
+         
          //for Global tracks
          Double_t nsigmaelectron = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(newAodTrack,(AliPID::EParticleType)AliPID::kElectron));
          Double_t nsigmapion = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(newAodTrack,(AliPID::EParticleType)AliPID::kPion));
index 067cf86..6c72e0d 100644 (file)
@@ -1063,8 +1063,8 @@ void  AliAnalysisTaskpxpy::UserExec(Option_t */*option*/)
   float trkCentr = -999.;
   float spdCentr = -999.;
   
-  float vertexX  = -999;
-  float vertexY  = -999;
+  //float vertexX  = -999;
+  //float vertexY  = -999;
   float vertexZ  = -999;
   //float vertexXY = -999;
   //float dcaZ     = -999;
@@ -1128,8 +1128,8 @@ void  AliAnalysisTaskpxpy::UserExec(Option_t */*option*/)
            {
              if(fCov[5] != 0)
                {
-                 vertexX = vertex->GetX();
-                 vertexY = vertex->GetY();
+                 //vertexX = vertex->GetX();
+                 //vertexY = vertex->GetY();
                  vertexZ = vertex->GetZ();
                  
                  if(TMath::Abs(vertexZ) > 10)
@@ -1190,12 +1190,13 @@ void  AliAnalysisTaskpxpy::UserExec(Option_t */*option*/)
          charge = int(q);
          phi    = t->Phi();
          pt     = t->Pt(); 
-         px     = t->Px();
-         py     = t->Py();
+         //px     = t->Px();
+         //py     = t->Py();
          pz     = t->Pz();
          eta    = t->Eta();
          dedx   = t->GetTPCsignal();
-                 
+         px     = pt*cos(phi);
+          py     = pt*sin(phi);
          //for Global tracks
          Double_t nsigmaelectron = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(newAodTrack,(AliPID::EParticleType)AliPID::kElectron));
          Double_t nsigmapion = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(newAodTrack,(AliPID::EParticleType)AliPID::kPion));
index dddb8eb..fec8efd 100644 (file)
@@ -1054,8 +1054,8 @@ void  AliAnalysisTaskpypy::UserExec(Option_t */*option*/)
   float trkCentr = -999.;
   float spdCentr = -999.;
   
-  float vertexX  = -999;
-  float vertexY  = -999;
+  //float vertexX  = -999;
+  //float vertexY  = -999;
   float vertexZ  = -999;
   //float vertexXY = -999;
   //float dcaZ     = -999;
@@ -1119,8 +1119,8 @@ void  AliAnalysisTaskpypy::UserExec(Option_t */*option*/)
            {
              if(fCov[5] != 0)
                {
-                 vertexX = vertex->GetX();
-                 vertexY = vertex->GetY();
+                 //vertexX = vertex->GetX();
+                 //vertexY = vertex->GetY();
                  vertexZ = vertex->GetZ();
                  
                  if(TMath::Abs(vertexZ) > 10)
@@ -1181,12 +1181,15 @@ void  AliAnalysisTaskpypy::UserExec(Option_t */*option*/)
          charge = int(q);
          phi    = t->Phi();
          pt     = t->Pt(); 
-         px     = t->Px();
-         py     = t->Py();
+         //px     = t->Px();
+         //py     = t->Py();
          pz     = t->Pz();
          eta    = t->Eta();
          dedx   = t->GetTPCsignal();
                  
+         px     = pt*cos(phi);
+          py     = pt*sin(phi);
+
          //for Global tracks
          Double_t nsigmaelectron = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(newAodTrack,(AliPID::EParticleType)AliPID::kElectron));
          Double_t nsigmapion = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(newAodTrack,(AliPID::EParticleType)AliPID::kPion));
index df2c44c..fb219f3 100644 (file)
@@ -1056,8 +1056,8 @@ void  AliAnalysisTaskpzpz::UserExec(Option_t */*option*/)
   float trkCentr = -999.;
   float spdCentr = -999.;
   
-  float vertexX  = -999;
-  float vertexY  = -999;
+  //float vertexX  = -999;
+  //float vertexY  = -999;
   float vertexZ  = -999;
   //float vertexXY = -999;
   //float dcaZ     = -999;
@@ -1121,8 +1121,8 @@ void  AliAnalysisTaskpzpz::UserExec(Option_t */*option*/)
            {
              if(fCov[5] != 0)
                {
-                 vertexX = vertex->GetX();
-                 vertexY = vertex->GetY();
+                 //vertexX = vertex->GetX();
+                 //vertexY = vertex->GetY();
                  vertexZ = vertex->GetZ();
                  
                  if(TMath::Abs(vertexZ) > 10)
@@ -1183,12 +1183,14 @@ void  AliAnalysisTaskpzpz::UserExec(Option_t */*option*/)
          charge = int(q);
          phi    = t->Phi();
          pt     = t->Pt(); 
-         px     = t->Px();
-         py     = t->Py();
+         //px     = t->Px();
+         //py     = t->Py();
          pz     = t->Pz();
          eta    = t->Eta();
          dedx   = t->GetTPCsignal();
-                 
+         px     = pt*cos(phi);
+          py     = pt*sin(phi);                  
+
          //for Global tracks
          Double_t nsigmaelectron = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(newAodTrack,(AliPID::EParticleType)AliPID::kElectron));
          Double_t nsigmapion = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(newAodTrack,(AliPID::EParticleType)AliPID::kPion));
index 5c12368..798f9e2 100644 (file)
@@ -52,6 +52,9 @@ void AddAliEbyENetChargeFluctuationTask(const Char_t *taskname="TOFTPC",
     return NULL;
   }
 
+
+  Printf("============================== I am here very great ================================");
+
   if (isMC) task->SetIsMC();
 
   Int_t sysii = sysType;
index 318d38c..b60c2af 100644 (file)
@@ -106,9 +106,9 @@ AliEbyENetChargeFluctuationTask::AliEbyENetChargeFluctuationTask(const char *nam
   fRapMax(0.5),   
   fCentralityBin(-1.),    
   fCentralityPercentile(-1.),
-  fNp(NULL),
-  fMCNp(NULL),
-  fRedFactp(NULL),           
+// fNp(NULL),
+// fMCNp(NULL),
+// fRedFactp(NULL),           
    
   fMinTrackLengthMC(80), 
   fAODtrackCutBit(768),   
@@ -136,7 +136,26 @@ AliEbyENetChargeFluctuationTask::AliEbyENetChargeFluctuationTask(const char *nam
   fRan(0),             
   fRanIdx(0),                        
   fHelperPID(0x0) { 
+
+  Printf("Task is Initialized");
+  
+
+  //-- -  - - -  - - ---- -- - --- -- --- -- - --
+  for (Int_t ii = 0; ii <= fOrder; ++ii)
+    for (Int_t jj = 0; jj < 2; ++jj)
+      fRedFactp[ii][jj] = 1.;
+  
+  //-- -  - - -  - - ---- -- - --- -- --- -- - -- 
+  for (Int_t kk = 0; kk < 4; ++kk)
+    for (Int_t jj = 0; jj < 2; ++jj)
+      fNp[kk][jj] = 0;
   
+  //-- -  - - -  - - ---- -- - --- -- --- -- - --
+  for (Int_t kk = 0; kk < 4; ++kk)
+    for (Int_t jj = 0; jj < 2; ++jj)
+      fMCNp[kk][jj] = 0;  
+   
+
   DefineOutput(1, TList::Class()); 
   DefineOutput(2, TList::Class()); 
   DefineOutput(3, TList::Class()); 
@@ -149,25 +168,7 @@ AliEbyENetChargeFluctuationTask::~AliEbyENetChargeFluctuationTask() {
   if (fPhyList)   delete fPhyList;
   if (fEffList)   delete fEffList;
   if (fDcaList)   delete fDcaList;
-  
-
-
-  for (Int_t kk = 0; kk < 2; ++kk)
-    if (fNp[kk]) delete[] fNp[kk];
-  if (fNp) delete[] fNp;
-
-  for (Int_t kk = 0; kk < 2; ++kk)
-    if (fMCNp[kk]) delete[] fMCNp[kk];
-  if (fMCNp) delete[] fMCNp;
-
-
-  for (Int_t ii = 0; ii <= fOrder; ++ii)
-    if (fRedFactp[ii]) delete[] fRedFactp[ii];
-  if (fRedFactp) delete[] fRedFactp;  
-  
   return;
-
-
 }
 
 
@@ -208,8 +209,8 @@ const Char_t* fgkNetHistLatex[4]      = {"+ + +","+","-","+ - -"};
 
 //---------------------------------------------------------------------------------
 void AliEbyENetChargeFluctuationTask::UserCreateOutputObjects() {
-  Bool_t oldStatus = TH1::AddDirectoryStatus();
-  TH1::AddDirectory(kFALSE);
+  //Bool_t oldStatus = TH1::AddDirectoryStatus();
+  //TH1::AddDirectory(kFALSE);
   
   fQaList = new TList();
   fQaList->SetOwner(kTRUE);
@@ -223,16 +224,18 @@ void AliEbyENetChargeFluctuationTask::UserCreateOutputObjects() {
   fDcaList = new TList();
   fDcaList->SetOwner(kTRUE);
   
-  fRedFactp = new Double_t*[fOrder+1];
-  for (Int_t ii = 0 ; ii <= fOrder; ++ii)
-    fRedFactp[ii] = new Double_t[2];
-  
+    
   fRan = new TRandom3();
   fRan->SetSeed();
   
   fRanIdx = new TRandom3();
   fRanIdx->SetSeed();
   
+  Printf(" >>>%d %d %d %d %d %d %d %d %d %d", 
+        fIsAOD, fIsMC, fIsPhy, fIsEff, 
+        fIsDca, fIsQa, 
+        fIsRatio, fIsSub, fIsBS, fIsPer);
+  
   CreateBasicQA();
   if (fIsQa)  CreateQA();
   if (fIsPhy) InitPhy();
@@ -242,7 +245,7 @@ void AliEbyENetChargeFluctuationTask::UserCreateOutputObjects() {
 
   //  if (fIsQa) if (fESDtrackCuts) fQaList->Add(fESDtrackCuts);
 
-  TH1::AddDirectory(oldStatus);
+  // TH1::AddDirectory(oldStatus);
 
   PostData(1, fPhyList); 
   PostData(2, fQaList);
@@ -262,27 +265,14 @@ void AliEbyENetChargeFluctuationTask::UserExec( Option_t * ){
     PostData(4, fEffList);
     return;
   }
+   
   
-  // Printf("Number of Track %d",fNTracks);
   
-  fRedFactp = new Double_t*[fOrder+1];
-  for (Int_t ii = 0 ; ii <= fOrder; ++ii)
-    fRedFactp[ii] = new Double_t[2];
   //-- -  - - -  - - ---- -- - --- -- --- -- - --
   for (Int_t ii = 0; ii <= fOrder; ++ii)
     for (Int_t jj = 0; jj < 2; ++jj)
       fRedFactp[ii][jj] = 1.;
   
-
-  fNp = new Double_t*[4];
-  for (Int_t ii = 0 ; ii < 4; ++ii)
-    fNp[ii] = new Double_t[2];
-
-  fMCNp = new Double_t*[4];
-  for (Int_t ii = 0 ; ii < 4; ++ii)
-    fMCNp[ii] = new Double_t[2];
-
   //-- -  - - -  - - ---- -- - --- -- --- -- - -- 
   for (Int_t kk = 0; kk < 4; ++kk)
     for (Int_t jj = 0; jj < 2; ++jj)
@@ -294,15 +284,26 @@ void AliEbyENetChargeFluctuationTask::UserExec( Option_t * ){
       fMCNp[kk][jj] = 0;  
   
   
+
+  // Printf("Number of Track %d",fNTracks);
+
   fSubSampleIdx = fRanIdx->Integer(fNSubSamples);
+
+  
+
   fNTracks  = (fESD) ? fESD->GetNumberOfTracks() : fAOD->GetNumberOfTracks();  
 
+  
+
   if (fIsMC && fIsAOD) {
     fArrayMC = dynamic_cast<TClonesArray*>(fAOD->FindListObject(AliAODMCParticle::StdBranchName()));
     if (!fArrayMC)
       AliFatal("No array of MC particles found !!!"); 
   }
 
+  
+
   //-- -  - - -  - - ---- -- - --- -- --- -- - --
   if(fSystemType == 0)      { 
     Execpp(); 
@@ -321,6 +322,9 @@ void AliEbyENetChargeFluctuationTask::UserExec( Option_t * ){
   else 
     return;
   
+
+  return;
+
   //-- -  - - -  - - ---- -- - --- -- --- -- - --
   PostData(1, fPhyList); 
   PostData(2, fQaList);
@@ -1137,8 +1141,21 @@ return;
 
 //________________________________________________________________________
 void AliEbyENetChargeFluctuationTask::FillBasicHistos(const Char_t *name, Bool_t isMC)  {
-  Double_t **np = (isMC) ? fMCNp : fNp;
+  //Double_t **np = (isMC) ? fMCNp : fNp;
   
+  Double_t np[4][2];
+  if (isMC) {
+    for (Int_t i = 0; i < 4; i++) {
+      np[i][0] = fMCNp[i][0];
+      np[i][1] = fMCNp[i][1];
+    }
+  } else {
+    for (Int_t i = 0; i < 4; i++) {
+      np[i][0] = fNp[i][0];
+      np[i][1] = fNp[i][1];
+    }
+  }
+
   Float_t centralityBin = fCentralityBin;
   Float_t centralityPer = fCentralityPercentile;
   
@@ -1249,7 +1266,22 @@ void AliEbyENetChargeFluctuationTask::FillBasicHistos(const Char_t *name, Bool_t
 //________________________________________________________________________
 void AliEbyENetChargeFluctuationTask::FillRatioHistos(const Char_t *name, Bool_t isMC,Bool_t isPer)  {
    
-  Double_t **np = (isMC) ? fMCNp : fNp;
+  //  Double_t **np = (isMC) ? fMCNp : fNp;
+
+  Double_t np[4][2];
+  if (isMC) {
+    for (Int_t i = 0; i < 4; i++) {
+      np[i][0] = fMCNp[i][0];
+      np[i][1] = fMCNp[i][1];
+    }
+  } else {
+    for (Int_t i = 0; i < 4; i++) {
+      np[i][0] = fNp[i][0];
+      np[i][1] = fNp[i][1];
+    }
+  }
+
+
 
   Float_t centralityBin = (isPer) ? fCentralityPercentile : fCentralityBin;
 
@@ -1394,8 +1426,21 @@ void AliEbyENetChargeFluctuationTask::FillGroupHistos(const Char_t *name,Int_t i
   TList *list    = static_cast<TList*>(fPhyList->FindObject(Form("f%s",name)));
   TList *listSub = static_cast<TList*>(list->FindObject(Form("%s%02d",name, iSub)));
 
-  Double_t **np = (isMC) ? fMCNp : fNp;
-      
+  // Double_t **np = (isMC) ? fMCNp : fNp;
+     
+  Double_t np[4][2];
+  if (isMC) {
+    for (Int_t i = 0; i < 4; i++) {
+      np[i][0] = fMCNp[i][0];
+      np[i][1] = fMCNp[i][1];
+    }
+  } else {
+    for (Int_t i = 0; i < 4; i++) {
+      np[i][0] = fNp[i][0];
+      np[i][1] = fNp[i][1];
+    }
+  }
   for (Int_t iPid = 0; iPid < 4; ++iPid) {
     Int_t deltaNp = np[iPid][1]-np[iPid][0];  
     Double_t delta = 1.;
index 776c97c..0246161 100644 (file)
@@ -16,6 +16,8 @@
 #include "TProfile2D.h"
 #include "TRandom3.h"
 
+class TList;
+
 class AliESDtrack;
 class AliMCEvent;
 class AliStack;
@@ -29,7 +31,6 @@ class AliAODTrack;
 class AliAODMCParticle;
 class AliMCParticle;
 class TClonesArray;
-class TList;
 class AliHelperPID;
 
 #include "AliAnalysisTaskSE.h"
@@ -38,7 +39,7 @@ class AliEbyENetChargeFluctuationTask: public AliAnalysisTaskSE {
 
  public:
 
-  AliEbyENetChargeFluctuationTask( const char *name = "HigherMomentAnalysis");
+  AliEbyENetChargeFluctuationTask( const char *name = "NetChargeFluctuation");
   virtual ~AliEbyENetChargeFluctuationTask();
   
   virtual void   UserCreateOutputObjects();
@@ -60,7 +61,7 @@ class AliEbyENetChargeFluctuationTask: public AliAnalysisTaskSE {
   void SetNSubSamples(Int_t i) {fNSubSamples = i;}
   void Debug() {fDebug = kTRUE;}
   void SetIsAOD() {fIsAOD = kTRUE;}
-  void SetAnal(Int_t i);
+  void SetAnal(Int_t i);// {fAnalType = i; };
   void SetHelperPID(AliHelperPID* pid){ fHelperPID = pid; }
  
  
@@ -135,74 +136,74 @@ class AliEbyENetChargeFluctuationTask: public AliAnalysisTaskSE {
   
   /* Global Members Private */
 
-  AliInputEventHandler *fInputEventHandler; //! 
-  AliESDEvent          *fESD;               //! 
-  AliAODEvent          *fAOD;               //! 
-  AliMCEvent           *fMCEvent;           //! 
-  AliStack             *fStack;             //! 
-  AliAODInputHandler   *fAODHandler;        //! 
-  AliESDInputHandler   *fESDHandler;        //! 
-  AliStack             *fMCStack;           //! 
-  TClonesArray         *fArrayMC;           //! 
-  AliESDtrackCuts      *fESDtrackCuts;      //!
-
-  TList     *fQaList;                       //! Tree of QA
-  TList     *fPhyList;                      //! For Analysis
-  TList     *fDcaList;                      //! dca for Both Data and MC 
-  TList     *fEffList;                      //! Correction List
-
-  Int_t      fSystemType;                   //! "pp", "pA", "AA"
-  TString    fCentralityEstimator;          //! "V0M","TRK","TKL","ZDC","FMD"  
-
-  Double_t   fVxMax;                        //! X vertex  Range
-  Double_t   fVyMax;                        //! Y vertex Range
-  Double_t   fVzMax;                        //! Z vertex Range
-  Double_t   fPhiMin;                       //! Phi Minimum
-  Double_t   fPhiMax;                       //! Phi Maximum
-  Double_t   fPtMin;                        //! pT Minimum
-  Double_t   fPtMax;                        //! Pt Maximum
-  Double_t   fEtaMin;                       //! Eta Minimum
-  Double_t   fEtaMax;                       //! Eta Maximum
-  Double_t   fRapMin;                       //! Rapidity Minimum
-  Double_t   fRapMax;                       //! Rapidity Maximum
-  Double_t   fCentralityBin;                //  Centrality bin of current event within max centrality bin
-  Double_t   fCentralityPercentile;         //  Centrality percentile of current event
-  Double_t **fNp;                           //  Array of particle/anti-particle counts
-  Double_t **fMCNp;                         //  Array of MC particle/anti-particle counts
-  Double_t **fRedFactp;                     //  Array of particle/anti-particle reduced factorial
-  Double_t   fCurGen[8];                    // 
-  Double_t   fCurGenD[8];                   // 
-  Double_t   fCurRecD[7];                   // 
-  Double_t   fCurRec[5];                    // 
-  Double_t   fCurCont[6];                   //
+  AliInputEventHandler *fInputEventHandler; // Input Handler
+  AliESDEvent          *fESD;               // Current ESD events
+  AliAODEvent          *fAOD;               // Current AOD events
+  AliMCEvent           *fMCEvent;           // Current MC Event
+  AliStack             *fStack;             // Stak tree
+  AliAODInputHandler   *fAODHandler;        // AOD handler
+  AliESDInputHandler   *fESDHandler;        // ESD handler
+  AliStack             *fMCStack;           // MC Stack
+  TClonesArray         *fArrayMC;           // AOD MC stack
+  AliESDtrackCuts      *fESDtrackCuts;      // ESD Track Cuts
+
+  TList     *fQaList;                       // Tree of QA
+  TList     *fPhyList;                      // For Analysis
+  TList     *fDcaList;                      // dca for Both Data and MC 
+  TList     *fEffList;                      // Correction List
+
+  Int_t      fSystemType;                   // Translated to Int_t 0:"pp", 1:"pA", 2:"AA"
+  TString    fCentralityEstimator;          // "V0M","TRK","TKL","ZDC","FMD"  
  
-  Float_t    fMinTrackLengthMC;             //!  Min track length for MC tracks
-  Int_t      fAODtrackCutBit;               //!
-  Int_t      fNSubSamples;                  //!  N subsamples
-  Int_t      fSubSampleIdx;                 //  Subsample idx for current event
-  Int_t      fOrder;                        //!  Max order of higher order distributions
-  Int_t      fNTriggers;                    //!  N triggers used
-  Int_t      fHEventStatMax;                //!  Max N cuts to be included in HEventSta
-  Int_t      fNCentralityBins;              //!  N centrality bins used
-  Int_t      fCentralityBinMax;             //!
-  Int_t      fNTracks;                      //
-
-  Bool_t     fIsMC;                         //!  Is MC event 
-  Bool_t     fIsRatio;                      //!  Is Ratio
-  Bool_t     fIsAOD;                        //!  analysis mode: 0 = ESDs  | 1 = AODs
-  Bool_t     fIsSub;                        //!  analysis mode SS  
-  Bool_t     fIsBS;                         //!  analysis mode BS
-  Bool_t     fIsPer;                        //!  analysis mode PER
-  Bool_t     fIsEff;                        //!  analysis mode Eff
-  Bool_t     fDebug;                        //!  Check for QA
-  Bool_t     fIsQa;                         //!  Check for QA
-  Bool_t     fIsPhy;                        //!  Check for QA
-  Bool_t     fIsDca;                        //!  Check for QA
-
-  TRandom3  *fRan;                          //!  Radom Number BS
-  TRandom3  *fRanIdx;                       //!  Random Number SS
-
-  AliHelperPID *fHelperPID;       //! Customised HelperPID class
+  Double_t   fVxMax;                        // X vertex  Range
+  Double_t   fVyMax;                        // Y vertex Range
+  Double_t   fVzMax;                        // Z vertex Range
+  Double_t   fPhiMin;                       // Phi Minimum
+  Double_t   fPhiMax;                       // Phi Maximum
+  Double_t   fPtMin;                        // pT Minimum
+  Double_t   fPtMax;                        // Pt Maximum
+  Double_t   fEtaMin;                       // Eta Minimum
+  Double_t   fEtaMax;                       // Eta Maximum
+  Double_t   fRapMin;                       // Rapidity Minimum
+  Double_t   fRapMax;                       // Rapidity Maximum
+  Double_t   fCentralityBin;                // Centrality bin of current event within max centrality bin
+  Double_t   fCentralityPercentile;         // Centrality percentile of current event
+  Double_t   fNp[4][2];                     // Array of particle/anti-particle counts
+  Double_t   fMCNp[4][2];                   // Array of MC particle/anti-particle counts
+  Double_t   fRedFactp[9][2];               // Array of particle/anti-particle reduced factorial
+  Double_t   fCurGen[8];                    // Current Gen Track
+  Double_t   fCurGenD[8];                   // Current info on DGen
+  Double_t   fCurRecD[7];                   // Current DRec
+  Double_t   fCurRec[5];                    // Current Rec info
+  Double_t   fCurCont[6];                   // Current Contamination
+
+  Float_t    fMinTrackLengthMC;             // Min track length for MC tracks
+  Int_t      fAODtrackCutBit;               // AOD BITs
+  Int_t      fNSubSamples;                  // N subsamples
+  Int_t      fSubSampleIdx;                 // Subsample idx for current event
+  Int_t      fOrder;                        // Max order of higher order distributions
+  Int_t      fNTriggers;                    // N triggers used
+  Int_t      fHEventStatMax;                // Max N cuts to be included in HEventSta
+  Int_t      fNCentralityBins;              // N centrality bins used
+  Int_t      fCentralityBinMax;             // Max Cent
+  Int_t      fNTracks;                      // Number of Tracks of Current Events
+
+  Bool_t     fIsMC;                         // Is MC event - Auto set by Add Task
+  Bool_t     fIsRatio;                      // Is Ratio
+  Bool_t     fIsAOD;                        // analysis mode: 0 = ESDs  | 1 = AODs
+  Bool_t     fIsSub;                        // analysis mode SS  
+  Bool_t     fIsBS;                         // analysis mode BS
+  Bool_t     fIsPer;                        // analysis mode PER
+  Bool_t     fIsEff;                        // analysis mode Eff
+  Bool_t     fDebug;                        // Check for QA
+  Bool_t     fIsQa;                         // Check for QA
+  Bool_t     fIsPhy;                        // Check for QA
+  Bool_t     fIsDca;                        // Check for QA
+
+  TRandom3  *fRan;                          // Radom Number BS
+  TRandom3  *fRanIdx;                       // Random Number SS
+
+  AliHelperPID *fHelperPID;                 // Customised HelperPID class
 
   AliEbyENetChargeFluctuationTask(const AliEbyENetChargeFluctuationTask&);
   AliEbyENetChargeFluctuationTask& operator = (const AliEbyENetChargeFluctuationTask&);
index 6abe315..1c622bf 100644 (file)
@@ -145,7 +145,7 @@ void AliEbyEParticleRatioFluctuationTask::UserExec( Option_t * ){
   
   AliAODHeader *aodHeader = dynamic_cast<AliAODHeader*>(event->GetHeader());
   if(!aodHeader) AliFatal("Not a standard AOD");
-  gCent = (Int_t)aodHeader->GetCentralityP()->GetCentralityPercentile(fCentralityEstimator.Data());
+  gCent = (Int_t)aodHeader->GetCentralityP()->GetCentralityPercentile(fCentralityEstimator.Data()); // V0M`
   gRefMul = aodHeader->GetRefMultiplicity();
   if (gCent < 0 || gCent > 100) return;
   if (isQA) fEventCounter->Fill(2);