Forgotten commit

[u/mrichter/AliRoot.git] / PWG0 / multPbPb / AliAnalysisMultPbTrackHistoManager.h

diff --git a/PWG0/multPbPb/AliAnalysisMultPbTrackHistoManager.h b/PWG0/multPbPb/AliAnalysisMultPbTrackHistoManager.h
index 0a52c2e49540bae8045c29f6affa50bdc65c0cd1..bde7afffba977cf977bbada18a2660531ef0ae4b 100644 (file)
--- a/PWG0/multPbPb/AliAnalysisMultPbTrackHistoManager.h
+++ b/PWG0/multPbPb/AliAnalysisMultPbTrackHistoManager.h
@@ -9,6 +9,10 @@ class TH1D;
 class TH1I;
 class AliMCParticle;
 
+//#define WEIGHTED_DCA
+#define TRANSVERSE_DCA
+
+
 //-------------------------------------------------------------------------
 //                      AliAnalysisMultPbTrackHistoManager
 // 
@@ -23,10 +27,15 @@ class AliAnalysisMultPbTrackHistoManager : public AliHistoListWrapper {
 
 public:
 
-  typedef enum {kHistoGen, kHistoRec, kHistoRecPrim, kHistoRecSecWeak, kHistoRecSecMat, kHistoRecFake, kNHistos} Histo_t;
-  typedef enum {kStatAll, kStatPhysSel, kStatCentr, kStatVtx, kNStatBins} Stat_t;
+  typedef enum {kHistoGen, kHistoRec, kHistoRecPrim, kHistoRecSecWeak, kHistoRecSecMat, kHistoRecFake, kHistoRecHighestMeanPt, kHistoRecLowestMeanPt, kNHistos} Histo_t;
+  typedef enum {kStatAll,  kStatCentr, kStatPhysSel, kStatVtx, kStatZDCCut, kStatVtxRangeCut, kNStatBins} Stat_t;
   typedef enum {kPartPiPlus, kPartKPlus, kPartP, kPartLPlus, kPartPiMinus, kPartKMinus, kPartPBar, kPartLMinus, kPartOther, kNPart} Part_t;
 
+  // these bits define the behaviour at merging
+  // if kKeepMax(Mib)Mean is set, the histogram with the highest (lowest) mean is kept
+  enum {kKeepMaxMean = BIT(22), kKeepMinMean = BIT(23)};
+
+
 
   AliAnalysisMultPbTrackHistoManager();
   AliAnalysisMultPbTrackHistoManager(const char * name,const char * title);
@@ -41,15 +50,19 @@ public:
   TH1D * GetHistoPt (Histo_t id, Float_t minEta = -22222, Float_t maxEta = -22222, Float_t minVz  = -22222, Float_t maxVz  = -22222, Bool_t scaleWidth = kTRUE);
   TH1D * GetHistoEta(Histo_t id, Float_t minPt  = -22222, Float_t maxPt  = -22222, Float_t minVz  = -22222, Float_t maxVz  = -22222, Bool_t scaleWidth = kTRUE);
   TH1D * GetHistoVz (Histo_t id, Float_t minPt  = -22222, Float_t maxPt  = -22222, Float_t minEta = -22222, Float_t maxEta = -22222, Bool_t scaleWidth = kTRUE);
+  TH2D * GetHistoPtVz (Histo_t id, Float_t minEta = -22222, Float_t maxEta = -22222, Bool_t scaleWidth = kFALSE);
 
   TH1I * GetHistoStats();
-  TH1D * GetHistoDCA(Histo_t id);
+  TH2D * GetHistoDCA(Histo_t id);
   TH1D * GetHistoMult(Histo_t id);
+  TH2D * GetHistoV0vsNtracks(Histo_t id);
+
 
   TH1D * GetHistoSpecies(Histo_t id);
   TH1D * GetHistoProcess(Histo_t id);
   TH1D * GetHistoVzEvent(Histo_t id);
-
+  TH1D * GetHistoPtEvent(Histo_t id);
+  TH1D * GetHistoMeanPt (Histo_t id);
 
 
   // Misch utils
@@ -59,13 +72,18 @@ public:
 
   // Histo bookers
   TH3D * BookHistoPtEtaVz(const char * name, const char * title);
-  TH1D * BookHistoDCA(const char * name, const char * title);
+  TH2D * BookHistoDCA(const char * name, const char * title);
   TH1I * BookHistoStats();
   TH1D * BookHistoMult(const char * name, const char * title);
 
   // 
   TH1 * GetHisto(const char * name);
 
+  virtual void Print(Option_t* option = "") const {fList->Print();}
+
+
+  Long64_t Merge(TCollection* list); // need to override this because of the mean pts
+
 private:
 
   static const char * kStatStepNames[];       // names of the step hist