]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWGLF/totEt/AliAnalysisEtReconstructed.h
add default options for kepler gpu family
[u/mrichter/AliRoot.git] / PWGLF / totEt / AliAnalysisEtReconstructed.h
index 718487d36476fed2972c9cde4578a98d986e1e11..96f6cbad9625b88467401f974dfb6738c63b0015 100644 (file)
@@ -1,10 +1,6 @@
 #ifndef ALIANALYSISETRECONSTRUCTED_H
 #define ALIANALYSISETRECONSTRUCTED_H
-
-class TH2F;
-class TH2D;
-class AliPHOSGeometry;
-class AliPHOSGeoUtils;//_________________________________________________________________________
+//_________________________________________________________________________
 //  Utility Class for transverse energy studies
 //  Base class for ESD analysis
 //  - reconstruction output
@@ -13,7 +9,9 @@ class AliPHOSGeoUtils;//________________________________________________________
 //_________________________________________________________________________
 
 #include "AliAnalysisEt.h"
-
+class TH2F;
+class TH3F;
+class TH2D;
 class TH2I;
 class AliVParticle;
 class AliESDEvent;
@@ -46,9 +44,8 @@ protected:
 
     bool CheckGoodVertex(AliVParticle *track);
     virtual bool TrackHitsCalorimeter(AliVParticle *track, Double_t magField);
-    virtual Bool_t TooCloseToBadChannel(const AliESDCaloCluster &cluster) const;
 
-    //Bool_t TooCloseToBadChannel(const AliESDCaloCluster &cluster) const;
+    virtual Double_t GetCorrectionModification(const AliESDCaloCluster& cluster,Int_t nonLinCorr, Int_t effCorr, Int_t cent);//nonLinCorr 0 = nominal 1 = high -1 = low, effCorr  0 = nominal 1 = high -1 = low
 
     AliAnalysisHadEtCorrections *fCorrections;//corrections needed for hadronic et
 
@@ -66,21 +63,59 @@ protected:
     Double_t fEMinCorrection; // Emin correction
     
     Double_t fRecEffCorrection; // Eff correction
-    Double_t CalcTrackClusterDistance(const Float_t pos[3],Int_t *trkMatchId, const AliESDEvent *event);
-    
-//    AliPHOSGeoUtils *fGeoUtils;
-    AliPHOSGeometry *fGeoUtils;
-    TH2I *fBadMapM2; // Bad map
-    TH2I *fBadMapM3; // Bad map
-    TH2I *fBadMapM4; // Bad map
-
-    
-    TH2D *fClusterPosition; // Position of clusters
     
-    TH2D *fHistChargedEnergyRemoved;
-    TH2D *fHistNeutralEnergyRemoved;
-    TH2D *fHistGammaEnergyAdded;
+    TH2D *fClusterPositionAccepted; // Position of clusters
+    TH2D *fClusterPositionAll; // Position of clusters
+    TH2D *fClusterPositionAcceptedEnergy; // Position of clusters
+    TH2D *fClusterPositionAllEnergy; // Position of clusters
+    TH1F *fClusterEnergy; // Distribution of cluster energies
+    TH2F *fClusterEnergyCent; // Distribution of cluster energies vs centrality bin
+    TH2F *fClusterEnergyCentMatched; // Distribution of cluster energies vs centrality bin
+    TH2F *fClusterEnergyCentNotMatched; // Distribution of cluster energies vs centrality bin
+    TH1F *fClusterEt; // Distribution of cluster energies
     
+    TH2D *fHistChargedEnergyRemoved; // Charged energy removed
+    TH2D *fHistNeutralEnergyRemoved; // Neutral energy removed
+    TH2D *fHistGammaEnergyAdded; // gamma energy added
+
+    TH3F *fHistMatchedTracksEvspTvsCent;   //For measuring hadron deposits
+    TH3F *fHistMatchedTracksEvspTvsCentEffCorr;   //For measuring hadron deposits
+    TH3F *fHistMatchedTracksEvspTvsCentEffTMCorr;   //For measuring hadron deposits
+    TH3F *fHistMatchedTracksEvspTvsCentEffTMCorr500MeV;   //For measuring hadron deposits
+    TH2F *fHistFoundHadronsvsCent;   //For measuring hadron deposits
+    TH2F *fHistNotFoundHadronsvsCent;   //For measuring hadron deposits
+    TH2F *fHistFoundHadronsEtvsCent;   //For measuring hadron deposits
+    TH2F *fHistNotFoundHadronsEtvsCent;   //For measuring hadron deposits
+    TH2F *fHistFoundHadronsvsCent500MeV;   //For measuring hadron deposits
+    TH2F *fHistNotFoundHadronsvsCent500MeV;   //For measuring hadron deposits
+    TH2F *fHistFoundHadronsEtvsCent500MeV;   //For measuring hadron deposits
+    TH2F *fHistNotFoundHadronsEtvsCent500MeV;   //For measuring hadron deposits
+    TH2D *fHistNominalRawEt;//Total ET from clusters with nominal reconstruction efficiency and nonlinearity correction vs centrality
+    TH2D *fHistNominalNonLinHighEt;//Total ET from clusters with nominal reconstruction efficiency and high bound of nonlinearity correction vs centrality
+    TH2D *fHistNominalNonLinLowEt;//Total ET from clusters with nominal reconstruction efficiency and low bound of nonlinearity correction vs centrality
+    TH2D *fHistNominalEffHighEt;//Total ET from clusters with high bound on reconstruction efficiency and nominal nonlinearity correction vs centrality
+    TH2D *fHistNominalEffLowEt;//Total ET from clusters with low bound on reconstruction efficiency and nominal nonlinearity correction vs centrality
+
+    TH2F *fHistTotRawEtEffCorr; // gamma efficiency applied
+    TH2F *fHistTotRawEt; //no  gamma efficiency applied
+    TH2F *fHistTotRawEtEffCorr500MeV;//Total ET from clusters with nominal reconstruction efficiency and nonlinearity correction vs centrality
+    TH2F *fHistTotAllRawEt; // all clusters no reco eff
+    TH2F *fHistTotAllRawEtEffCorr; // all clusters reco eff applied
+    Double_t ApplyModifiedCorrections(const AliESDCaloCluster& cluster,Int_t nonLinCorr, Int_t effCorr, Int_t cent);//nonLinCorr 0 = nominal 1 = high -1 = low, effCorr  0 = nominal 1 = high -1 = low
+
+    TH3F *fHistNClustersPhosVsEmcal; // all clusters no reco eff
+    TH2F *fHistClusterSizeVsCent; // all clusters no reco eff
+    TH2F *fHistMatchedClusterSizeVsCent; // all clusters no reco eff
+    TH2F *fHistTotAllRawEtVsTotalPt; // all clusters no reco eff
+    //fHistTotAllRawEtVsTotalPtCent
+    TH3F *fHistTotAllRawEtVsTotalPtVsCent; // all clusters no reco eff
+    TH3F *fHistTotMatchedRawEtVsTotalPtVsCent; // all clusters no reco eff
+    TH2F *fHistPIDProtonsTrackMatchedDepositedVsNch;
+    TH2F *fHistPIDAntiProtonsTrackMatchedDepositedVsNch;
+    TH2F *fHistPIDProtonsTrackMatchedDepositedVsNcl;
+    TH2F *fHistPIDAntiProtonsTrackMatchedDepositedVsNcl;
+    TH2F *fHistPiKPTrackMatchedDepositedVsNch;
+    TH3F *fHistCentVsNchVsNclReco;
 
 private: