ALIROOT-5420 Missing include

[u/mrichter/AliRoot.git] / PWGLF / totEt / AliAnalysisEtCuts.h

diff --git a/PWGLF/totEt/AliAnalysisEtCuts.h b/PWGLF/totEt/AliAnalysisEtCuts.h
index bb00ebf605b0dfdc7c5c64342fd804d8eb7f920f..6394116b736d4507c28547cb44351b1d3e66d1fc 100644 (file)
--- a/PWGLF/totEt/AliAnalysisEtCuts.h
+++ b/PWGLF/totEt/AliAnalysisEtCuts.h
@@ -9,6 +9,7 @@
 //_________________________________________________________________________
 
 #include "TNamed.h"
+#include <iostream>
 
 class AliAnalysisEtCuts : public TNamed
 {
@@ -46,14 +47,18 @@ class AliAnalysisEtCuts : public TNamed
   Int_t GetReconstructedNItsClustersCut() const { return fReconstructedNItsClustersCut; }
   Double_t GetReconstructedPidCut() const { return fReconstructedPidCut; }
   // ReconstructedPhos
-  Char_t GetReconstructedPhosClusterType() const { return fReconstructedPhosClusterType; }
+  Char_t GetPhosClusterType() const { return fReconstructedPhosClusterType; }
   Double_t GetReconstructedPhosClusterEnergyCut() const { return fReconstructedPhosClusterEnergyCut; }
   Double_t GetReconstructedPhosSingleCellEnergyCut() const { return fReconstructedPhosSingleCellEnergyCut; }
   Double_t GetPhosTrackDistanceCut() const { return fPhosTrackDistanceCut; }
   Double_t GetPhosTrackDxCut() const { return fPhosTrackDxCut; }
   Double_t GetPhosTrackDzCut() const { return fPhosTrackDzCut; }
+  Double_t GetPhosTrackRCut() const { return fPhosTrackRCut; }
+  
+  Double_t GetPhosBadDistanceCut() const { return fPhosBadDistanceCut; }
+  
   // ReconstructedEmcal
-  Char_t GetReconstructedEmcalClusterType() const { return fReconstructedEmcalClusterType; }
+  Char_t GetEmcalClusterType() const { return fReconstructedEmcalClusterType; }
   Double_t GetReconstructedEmcalClusterEnergyCut() const { return fReconstructedEmcalClusterEnergyCut; }
   Double_t GetReconstructedEmcalSingleCellEnergyCut() const { return fReconstructedEmcalSingleCellEnergyCut; }
   Double_t GetEmcalTrackDistanceCut() const { return fEmcalTrackDistanceCut; }
@@ -83,71 +88,84 @@ class AliAnalysisEtCuts : public TNamed
   Double_t GetHistMinParticlePt() const { return fHistMinParticlePt; }
   Double_t GetHistMaxParticlePt() const { return fHistMaxParticlePt; }
   
+  
+  
   Short_t GetDetectorPhos() const { return fgkDetectorPhos; }
   Short_t GetDetectorEmcal() const { return fgkDetectorEmcal; }
+  
+  Double_t GetPrimaryVertexCutXY() const { return fPrimaryVertexCutXY; }
+  Double_t GetPrimaryVertexCutZ() const { return fPrimaryVertexCutZ; }
+  
 
   // Setters
   // Common
-  void SetCommonEtaCut(const Double_t val) { fCommonEtaCut = val; }
-  void SetCommonClusterEnergyCut(const Double_t val) { fCommonClusterEnergyCut = val; }
-  void SetCommonTrackPtCut(const Double_t val) { fCommonTrackPtCut = val; }
-  void SetCommonSingleCell(const Int_t val) { fCommonSingleCell = val;}
+  void SetCommonEtaCut(Double_t val) { fCommonEtaCut = val; }
+  void SetCommonClusterEnergyCut(Double_t val) { fCommonClusterEnergyCut = val; }
+  void SetCommonTrackPtCut(Double_t val) { fCommonTrackPtCut = val; }
+  void SetCommonSingleCell(Int_t val) { fCommonSingleCell = val;}
   // GeometryPhos
-  void SetGeometryPhosEtaAccCut(const Double_t val) { fGeometryPhosEtaAccCut = val; }
-  void SetGeometryPhosPhiAccMinCut(const Double_t val) { fGeometryPhosPhiAccMinCut = val; }
-  void SetGeometryPhosPhiAccMaxCut(const Double_t val) { fGeometryPhosPhiAccMaxCut = val; }
-  void SetGeometryPhosDetectorRadius(const Double_t val) { fGeometryPhosDetectorRadius = val; }
+  void SetGeometryPhosEtaAccCut(Double_t val) { fGeometryPhosEtaAccCut = val; }
+  void SetGeometryPhosPhiAccMinCut(Double_t val) { fGeometryPhosPhiAccMinCut = val; }
+  void SetGeometryPhosPhiAccMaxCut(Double_t val) { fGeometryPhosPhiAccMaxCut = val; }
+  void SetGeometryPhosDetectorRadius(Double_t val) { fGeometryPhosDetectorRadius = val; }
   // GeometryEmcal
-  void SetGeometryEmcalEtaAccCut(const Double_t val) { fGeometryEmcalEtaAccCut = val; }
-  void SetGeometryEmcalPhiAccMinCut(const Double_t val) { fGeometryEmcalPhiAccMinCut = val; }
-  void SetGeometryEmcalPhiAccMaxCut(const Double_t val) { fGeometryEmcalPhiAccMaxCut = val; }
-  void SetGeometryEmcalDetectorRadius(const Double_t val) { fGeometryEmcalDetectorRadius = val; }
+  void SetGeometryEmcalEtaAccCut(Double_t val) { fGeometryEmcalEtaAccCut = val; }
+  void SetGeometryEmcalPhiAccMinCut(Double_t val) { fGeometryEmcalPhiAccMinCut = val; }
+  void SetGeometryEmcalPhiAccMaxCut(Double_t val) { fGeometryEmcalPhiAccMaxCut = val; }
+  void SetGeometryEmcalDetectorRadius(Double_t val) { fGeometryEmcalDetectorRadius = val; }
   // Reconstructed
-  void SetReconstructedVertexXCut(const Double_t val) { fReconstructedVertexXCut = val; }
-  void SetReconstructedVertexYCut(const Double_t val) { fReconstructedVertexYCut = val; }
-  void SetReconstructedVertexZCut(const Double_t val) { fReconstructedVertexZCut = val; }
-  void SetReconstructedIPxyCut(const Double_t val) { fReconstructedIPxyCut = val; }
-  void SetReconstructedIPzCut(const Double_t val) { fReconstructedIPzCut = val; }
-  void SetReconstructedNTpcClustersCut(const Int_t val) { fReconstructedNTpcClustersCut = val; }
-  void SetReconstructedNItsClustersCut(const Int_t val) { fReconstructedNItsClustersCut = val; }
-  void SetReconstrucedPidCut(const Double_t val) { fReconstructedPidCut = val; }
+  void SetReconstructedVertexXCut(Double_t val) { fReconstructedVertexXCut = val; }
+  void SetReconstructedVertexYCut(Double_t val) { fReconstructedVertexYCut = val; }
+  void SetReconstructedVertexZCut(Double_t val) { fReconstructedVertexZCut = val; }
+  void SetReconstructedIPxyCut(Double_t val) { fReconstructedIPxyCut = val; }
+  void SetReconstructedIPzCut(Double_t val) { fReconstructedIPzCut = val; }
+  void SetReconstructedNTpcClustersCut(Int_t val) { fReconstructedNTpcClustersCut = val; }
+  void SetReconstructedNItsClustersCut(Int_t val) { fReconstructedNItsClustersCut = val; }
+  void SetReconstrucedPidCut(Double_t val) { fReconstructedPidCut = val; }
   // ReconstructedPhos
-  void SetReconstructedPhosClusterType(const Char_t val) { fReconstructedPhosClusterType = val; }
-  void SetReconstructedPhosClusterEnergyCut(const Double_t val) { fReconstructedPhosClusterEnergyCut = val; }
-  void SetReconstructedPhosSingleCellEnergyCut(const Double_t val) { fReconstructedPhosSingleCellEnergyCut = val; }
-  void SetPhosTrackDistanceCut(const Double_t val) { fPhosTrackDistanceCut = val; }
-  void SetPhosTrackDxCut(const Double_t val) { fPhosTrackDxCut = val; }
-  void SetPhosTrackDzCut(const Double_t val) { fPhosTrackDzCut = val; }
+  void SetReconstructedPhosClusterType(Char_t val) { fReconstructedPhosClusterType = val; }
+  void SetReconstructedPhosClusterEnergyCut(Double_t val) { fReconstructedPhosClusterEnergyCut = val; }
+  void SetReconstructedPhosSingleCellEnergyCut(Double_t val) { fReconstructedPhosSingleCellEnergyCut = val; }
+  void SetPhosTrackDistanceCut(Double_t val) { fPhosTrackDistanceCut = val; }
+  void SetPhosTrackDxCut(Double_t val) { fPhosTrackDxCut = val; }
+  void SetPhosTrackDzCut(Double_t val) { fPhosTrackDzCut = val; }
+  void SetPhosTrackRCut(Double_t val) { std::cout << "Setting: " << val << std::endl; fPhosTrackRCut = val; }
+  
+  void SetPhosBadDistanceCut(Double_t val) { fPhosBadDistanceCut = val; }
   
   // ReconstructedEmcal
-  void SetReconstructedEmcalClusterType(const Char_t val) { fReconstructedEmcalClusterType = val; }
-  void SetReconstructedEmcalClusterEnergyCut(const Double_t val) { fReconstructedEmcalClusterEnergyCut = val; }
-  void SetReconstructedEmcalSingleCellEnergyCut(const Double_t val) { fReconstructedEmcalSingleCellEnergyCut = val; }
-  void SetEmcalTrackDistanceCut(const Double_t val) { fEmcalTrackDistanceCut = val; }
+  void SetReconstructedEmcalClusterType(Char_t val) { fReconstructedEmcalClusterType = val; }
+  void SetReconstructedEmcalClusterEnergyCut(Double_t val) { fReconstructedEmcalClusterEnergyCut = val; }
+  void SetReconstructedEmcalSingleCellEnergyCut(Double_t val) { fReconstructedEmcalSingleCellEnergyCut = val; }
+  void SetEmcalTrackDistanceCut(Double_t val) { fEmcalTrackDistanceCut = val; }
   // MonteCarlo
-  void SetMonteCarloSingleChargedParticle(const Double_t val) { fMonteCarloSingleChargedParticle = val; }
-  void SetMonteCarloNeutralParticle(const Double_t val) { fMonteCarloNeutralParticle = val; }
+  void SetMonteCarloSingleChargedParticle(Double_t val) { fMonteCarloSingleChargedParticle = val; }
+  void SetMonteCarloNeutralParticle(Double_t val) { fMonteCarloNeutralParticle = val; }
   // Hist: TTree and histogram info
-  void SetHistMakeTree(const Bool_t val) { fHistMakeTree = val; }
-  void SetHistMakeTreeDeposit(const Bool_t val) { fHistMakeTreeDeposit = val; }
+  void SetHistMakeTree(Bool_t val) { fHistMakeTree = val; }
+  void SetHistMakeTreeDeposit(Bool_t val) { fHistMakeTreeDeposit = val; }
   //
-  void SetHistNbinsMult(const Int_t val) { fHistNbinsMult = val; }
-  void SetHistMinMult(const Double_t val) { fHistMinMult = val; }
-  void SetHistMaxMult(const Double_t val) { fHistMaxMult = val; }
+  void SetHistNbinsMult(Int_t val) { fHistNbinsMult = val; }
+  void SetHistMinMult(Double_t val) { fHistMinMult = val; }
+  void SetHistMaxMult(Double_t val) { fHistMaxMult = val; }
   //
-  void SetHistNbinsTotEt(const Int_t val) { fHistNbinsTotEt = val; }
-  void SetHistMinTotEt(const Double_t val) { fHistMinTotEt = val; }
-  void SetHistMaxTotEt(const Double_t val) { fHistMaxTotEt = val; }
+  void SetHistNbinsTotEt(Int_t val) { fHistNbinsTotEt = val; }
+  void SetHistMinTotEt(Double_t val) { fHistMinTotEt = val; }
+  void SetHistMaxTotEt(Double_t val) { fHistMaxTotEt = val; }
   //
-  void SetHistNbinsParticleEt(const Int_t val) { fHistNbinsParticleEt = val; }
-  void SetHistMinParticleEt(const Double_t val) { fHistMinParticleEt = val; }
-  void SetHistMaxParticleEt(const Double_t val) { fHistMaxParticleEt = val; }
+  void SetHistNbinsParticleEt(Int_t val) { fHistNbinsParticleEt = val; }
+  void SetHistMinParticleEt(Double_t val) { fHistMinParticleEt = val; }
+  void SetHistMaxParticleEt(Double_t val) { fHistMaxParticleEt = val; }
   //
-  void SetHistNbinsParticlePt(const Int_t val) { fHistNbinsParticlePt = val; }
-  void SetHistMinParticlePt(const Double_t val) { fHistMinParticlePt = val; }
-  void SetHistMaxParticlePt(const Double_t val) { fHistMaxParticlePt = val; }
-
+  void SetHistNbinsParticlePt(Int_t val) { fHistNbinsParticlePt = val; }
+  void SetHistMinParticlePt(Double_t val) { fHistMinParticlePt = val; }
+  void SetHistMaxParticlePt(Double_t val) { fHistMaxParticlePt = val; }
 
+  void SetPrimaryVertexCutXY(Double_t val) { fPrimaryVertexCutXY = val; }
+  void SetPrimaryVertexCutZ(Double_t val) { fPrimaryVertexCutZ = val; }
+  
+  
+  
  protected:
 
   // Common   
@@ -162,7 +180,10 @@ class AliAnalysisEtCuts : public TNamed
   Double_t fPhosTrackDistanceCut; // PHOS track distance  
   Double_t fPhosTrackDxCut; // PHOS track distance in x
   Double_t fPhosTrackDzCut; // PHOS track distance  in z
-  
+  Double_t fPhosTrackRCut; // PHOS track distance  in r (using the parametrized track distance)
+ 
+ Double_t fPhosBadDistanceCut; // PHOS distance to bad channel 
+ 
   // GeometryPhos
   Double_t fGeometryPhosEtaAccCut; // PHOS Eta Acc cut
   Double_t fGeometryPhosPhiAccMinCut; // PHOS Phi Acc Min cut
@@ -228,6 +249,10 @@ class AliAnalysisEtCuts : public TNamed
 // Detector definition
   static const Short_t fgkDetectorPhos = -1; // PHOS 
   static const Short_t fgkDetectorEmcal = 1; // EMCAL 
+  
+  Double_t fPrimaryVertexCutXY; // Cut to decide if particle is from primary vertex
+  Double_t fPrimaryVertexCutZ; // Cut to decide if particle is from primary vertex
+  
 
 private:
   //Declare private to avoid compilation warning