totEt cuts into their own proper class

[u/mrichter/AliRoot.git] / PWG4 / totEt / AliAnalysisHadEt.cxx

diff --git a/PWG4/totEt/AliAnalysisHadEt.cxx b/PWG4/totEt/AliAnalysisHadEt.cxx
index 961a762c59b26a1126bac7958d9e49231ff07441..3d086c259a58268b01cfec791b32b7acd13fbb4d 100644 (file)
--- a/PWG4/totEt/AliAnalysisHadEt.cxx
+++ b/PWG4/totEt/AliAnalysisHadEt.cxx
@@ -1,8 +1,13 @@
-//Create by Christine Nattrass, Rebecca Scott, Irakli Martashvili
+//Created by Christine Nattrass, Rebecca Scott, Irakli Martashvili
 //University of Tennessee at Knoxville
-//This class is designed for the analysis of the hadronic component of transverse energy.  It is used by AliAnalysisTaskHadEt.
-//This gets information about the hadronic component of the transverse energy from tracks reconstructed in an event
-//it has daughters, AliAnalysisHadEtMonteCarlo and AliAnalysisHadEtReconstructed, which loop over either Monte Carlo data or real data to get Et
+//
+// This class is designed for the analysis of the hadronic component of 
+// transverse energy.  It is used by AliAnalysisTaskHadEt.
+// This gets information about the hadronic component of the transverse energy 
+// from tracks reconstructed in an event
+// it has daughters, AliAnalysisHadEtMonteCarlo and 
+// AliAnalysisHadEtReconstructed which loop over either Monte Carlo data or 
+// real data to get Et
 #include "AliAnalysisHadEt.h"
 #include "TMath.h"
 #include "TList.h"
@@ -38,6 +43,7 @@ Float_t AliAnalysisHadEt::fgPtAxis[117]=
 
 AliAnalysisHadEt::AliAnalysisHadEt() :
         fHistogramNameSuffix("")
+       ,fCuts(0)
         ,fPdgDB(0)
         ,fPiPlusCode(0)
         ,fPiMinusCode(0)
@@ -73,21 +79,9 @@ AliAnalysisHadEt::AliAnalysisHadEt() :
        ,fMultiplicity(0)
        ,fChargedMultiplicity(0)
        ,fNeutralMultiplicity(0)
-        ,fEtaCut(EtCommonCuts::kEtaCut)
-        ,fEtaCutAcc(0)
-                                 //,fPhiCutAccMin(0)
-                                 //,fPhiCutAccMax(360.)
-        ,fVertexXCut(0)
-        ,fVertexYCut(0)
-        ,fVertexZCut(0)
-                                 ,fIPxyCut(0)
-                                 ,fIPzCut(0)
-                                 //,fSingleCellEnergyCut(0)
-                                 //,fClusterEnergyCut(EtCommonCuts::kClusterEnergyCut)
-                                 //,fTrackPtCut(EtCommonCuts::kTrackPtCut)
-        ,ffesdtrackCutsITSTPC(0)
-        ,fesdtrackCutsTPC(0)
-       ,fesdtrackCutsITS(0)
+        ,fEsdtrackCutsITSTPC(0)
+        ,fEsdtrackCutsTPC(0)
+       ,fEsdtrackCutsITS(0)
         ,fhistoList(0)
 {//default constructor
 
@@ -129,31 +123,10 @@ void AliAnalysisHadEt::FillOutputList()
 
 void AliAnalysisHadEt::Init()
 {//Initiate member vaiables to reasonable values
+  if (!fCuts) fCuts = new AliAnalysisEtCuts();
+
   if(!fPdgDB) fPdgDB = new TDatabasePDG();
-  //the codes are defined in $ROOTSYS/etc/pdg_table.txt
-  fPionMass = fPdgDB->GetParticle("pi+")->Mass();
-  fPiPlusCode = fPdgDB->GetParticle("pi+")->PdgCode();
-    fPiMinusCode = fPdgDB->GetParticle("pi-")->PdgCode();
-    fKPlusCode = fPdgDB->GetParticle("K+")->PdgCode();
-    fKMinusCode = fPdgDB->GetParticle("K-")->PdgCode();
-    fProtonCode = fPdgDB->GetParticle("proton")->PdgCode();
-    fAntiProtonCode = fPdgDB->GetParticle("antiproton")->PdgCode();
-    fLambdaCode = fPdgDB->GetParticle("Lambda0")->PdgCode();
-    fAntiLambdaCode = fPdgDB->GetParticle("Lambda0_bar")->PdgCode();
-    fK0SCode = fPdgDB->GetParticle("K_S0")->PdgCode();
-    fOmegaCode = fPdgDB->GetParticle("Omega-")->PdgCode();
-    fAntiOmegaCode = fPdgDB->GetParticle("Omega+")->PdgCode();
-    fXi0Code = fPdgDB->GetParticle("Xi0")->PdgCode();
-    fAntiXi0Code = fPdgDB->GetParticle("Xi0_bar")->PdgCode();
-    fXiCode = fPdgDB->GetParticle("Xi-")->PdgCode();
-    fAntiXiCode = fPdgDB->GetParticle("Xi-_bar")->PdgCode();
-    fSigmaCode = fPdgDB->GetParticle("Sigma-")->PdgCode();
-    fAntiSigmaCode = fPdgDB->GetParticle("Sigma+")->PdgCode();
-    fK0LCode = fPdgDB->GetParticle("K_L0")->PdgCode();
-    fNeutronCode = fPdgDB->GetParticle("neutron")->PdgCode();
-    fAntiNeutronCode = fPdgDB->GetParticle("antineutron")->PdgCode();
-    fEPlusCode = fPdgDB->GetParticle("e+")->PdgCode();
-    fEMinusCode = fPdgDB->GetParticle("e-")->PdgCode();
+  SetParticleCodes();  
 }
 
 void AliAnalysisHadEt::CreateHistograms()
@@ -176,6 +149,34 @@ void AliAnalysisHadEt::ResetEventValues()
     fChargedMultiplicity = 0;
     fNeutralMultiplicity = 0;
 }
+
+void AliAnalysisHadEt::SetParticleCodes()
+{  //the codes are defined in $ROOTSYS/etc/pdg_table.txt
+  fPionMass = fPdgDB->GetParticle("pi+")->Mass();
+  fPiPlusCode = fPdgDB->GetParticle("pi+")->PdgCode();
+  fPiMinusCode = fPdgDB->GetParticle("pi-")->PdgCode();
+  fKPlusCode = fPdgDB->GetParticle("K+")->PdgCode();
+  fKMinusCode = fPdgDB->GetParticle("K-")->PdgCode();
+  fProtonCode = fPdgDB->GetParticle("proton")->PdgCode();
+  fAntiProtonCode = fPdgDB->GetParticle("antiproton")->PdgCode();
+  fLambdaCode = fPdgDB->GetParticle("Lambda0")->PdgCode();
+  fAntiLambdaCode = fPdgDB->GetParticle("Lambda0_bar")->PdgCode();
+  fK0SCode = fPdgDB->GetParticle("K_S0")->PdgCode();
+  fOmegaCode = fPdgDB->GetParticle("Omega-")->PdgCode();
+  fAntiOmegaCode = fPdgDB->GetParticle("Omega+")->PdgCode();
+  fXi0Code = fPdgDB->GetParticle("Xi0")->PdgCode();
+  fAntiXi0Code = fPdgDB->GetParticle("Xi0_bar")->PdgCode();
+  fXiCode = fPdgDB->GetParticle("Xi-")->PdgCode();
+  fAntiXiCode = fPdgDB->GetParticle("Xi-_bar")->PdgCode();
+  fSigmaCode = fPdgDB->GetParticle("Sigma-")->PdgCode();
+  fAntiSigmaCode = fPdgDB->GetParticle("Sigma+")->PdgCode();
+  fK0LCode = fPdgDB->GetParticle("K_L0")->PdgCode();
+  fNeutronCode = fPdgDB->GetParticle("neutron")->PdgCode();
+  fAntiNeutronCode = fPdgDB->GetParticle("antineutron")->PdgCode();
+  fEPlusCode = fPdgDB->GetParticle("e+")->PdgCode();
+  fEMinusCode = fPdgDB->GetParticle("e-")->PdgCode();
+}
+
 void AliAnalysisHadEt::CreateEtaPtHisto2D(TString name, TString title)
 {     //creates a 2-d histogram in eta and phi and adds it to the list of histograms to be saved
   TString *histoname   = new TString();