Disable Evchar if no TR read.

[u/mrichter/AliRoot.git] / TPC / AliTPCcalibCosmic.h

diff --git a/TPC/AliTPCcalibCosmic.h b/TPC/AliTPCcalibCosmic.h
index 5b7729f2662c11e68fc4e5c6e7bc9384f46c581f..6977bbd4c0d3acd067382257de1fae8dac75554a 100644 (file)
--- a/TPC/AliTPCcalibCosmic.h
+++ b/TPC/AliTPCcalibCosmic.h
@@ -5,10 +5,13 @@
  * See cxx source for full Copyright notice                               */
 
 #include "AliTPCcalibBase.h"
+#include "AliTPCCalPad.h"
+#include "TH2F.h"
 class TH1F;
-class TH2F;
 class TList;
 class AliESDEvent;
+class AliESDtrack;
+class THnSparse;
 
 #include "TTreeStream.h"
 
@@ -19,24 +22,75 @@ public:
   AliTPCcalibCosmic(const Text_t *name, const Text_t *title);
   virtual ~AliTPCcalibCosmic();
   
-  virtual void     Process(AliESDEvent *event);
-  virtual Long64_t Merge(TCollection *li);
+  virtual void      Process(AliESDEvent *event);
+  virtual Long64_t  Merge(TCollection *li);
+  virtual void      Analyze();
+  void              Add(const AliTPCcalibCosmic* cosmic);
+  //
+  //
+  void              Init();
+  void              FindPairs(AliESDEvent *event);
+  Bool_t            IsPair(AliExternalTrackParam *tr0, AliExternalTrackParam *tr1);
+  static void       CalculateBetheParams(TH2F *hist, Double_t * initialParam);
+  static Double_t   CalculateMIPvalue(TH1F * hist);
+  AliExternalTrackParam *MakeTrack(const AliExternalTrackParam *track0, const AliExternalTrackParam *track1);
+  AliExternalTrackParam *MakeCombinedTrack(const AliExternalTrackParam *track0, const AliExternalTrackParam *track1);
+
+  void UpdateTrack(AliExternalTrackParam &track0, const AliExternalTrackParam &track1);
+  //
+  void FillHistoPerformance(AliExternalTrackParam *par0, AliExternalTrackParam *par1, AliExternalTrackParam *inner0, AliExternalTrackParam *inner1, AliTPCseed *seed0,  AliTPCseed *seed1);
+
+
+  //
+  TH1F   *          GetHistNTracks(){return fHistNTracks;};
+  TH1F   *          GetHistClusters(){return fClusters;};
+  TH2F   *          GetHistAcorde(){return fModules;};
+  TH1F   *          GetHistPt(){return fHistPt;};
+  TH2F   *          GetHistDeDx(){return fDeDx;};
+  TH1F   *          GetHistMIP(){return fDeDxMIP;};
+  //
+  Double_t          GetMIPvalue(){return fMIPvalue;};
+  //
+  static void       BinLogX(TH1 * h);   // method for correct histogram binning
+  static void       BinLogX(THnSparse * h, Int_t axisDim);   // method for correct histogram binning
+
+  void     Process(AliESDtrack *track, Int_t runNo=-1){AliTPCcalibBase::Process(track,runNo);};
+  void     Process(AliTPCseed *track){return AliTPCcalibBase::Process(track);}
   
+  //
+  // Performance histograms
+  //
+  THnSparse   *fHistoDelta[6];  // histograms of tracking performance delta
+  THnSparse   *fHistoPull[6];   // histograms of tracking performance pull
+  THnSparse   *fHistodEdxMax[4];   // histograms of dEdx perfomance - max charge
+  THnSparse   *fHistodEdxTot[4];   // histograms of dEdx perfomance - tot charge
+
+private:
   
-public:
-  static void BinLogX(TH1 * h);   // method for correct histogram binning
+  void              FillAcordeHist(AliESDtrack *upperTrack);
+
+  
+
+  TH1F  *fHistNTracks;            //  histogram showing number of ESD tracks per event
+  TH1F  *fClusters;               //  histogram showing the number of clusters per track
+  TH2F  *fModules;                //  2d histogram of tracks which are propagated to the ACORDE scintillator array
+  TH1F  *fHistPt;                 //  Pt histogram of reconstructed tracks
+  TH2F  *fDeDx;                   //  dEdx spectrum showing the different particle types
+  TH1F  *fDeDxMIP;                //  TPC signal close to the MIP region of muons 0.4 < p < 0.45 GeV
 
-  TH1F  *fHistNTracks;
-  TH1F  *fClusters;
-  TH2F  *fModules;
-  TH1F  *fHistPt;
-  TH1F  *fPtResolution;
-  TH2F  *fDeDx;
+  Double_t fMIPvalue;             //  MIP value calculated via a fit to fDeDxMIP
+  //
+  // cuts
+  //
+  Float_t fCutMaxD;     // maximal distance in rfi ditection
+  Float_t fCutMaxDz;     // maximal distance in z ditection
+  Float_t fCutTheta;    // maximal distance in theta ditection
+  Float_t fCutMinDir;   // direction vector products
 
   AliTPCcalibCosmic(const AliTPCcalibCosmic&); 
   AliTPCcalibCosmic& operator=(const AliTPCcalibCosmic&); 
 
-  ClassDef(AliTPCcalibCosmic, 1); 
+  ClassDef(AliTPCcalibCosmic, 2); 
 };
 
 #endif