caching of the correction coeficients from the graph

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

diff --git a/TPC/AliTPCcalibCosmic.h b/TPC/AliTPCcalibCosmic.h
index 38c4af8d797b9814be8d2cd74a03e4a1d2e99034..ffee9e5e13eccc559f3582fa31b3ab67acd38da6 100644 (file)
--- a/TPC/AliTPCcalibCosmic.h
+++ b/TPC/AliTPCcalibCosmic.h
@@ -5,15 +5,12 @@
  * See cxx source for full Copyright notice                               */
 
 #include "AliTPCcalibBase.h"
-#include "AliTPCCalPad.h"
-#include "TH2F.h"
+class TH2F;
 class TH1F;
 class TList;
 class AliESDEvent;
 class AliESDtrack;
-
-#include "TTreeStream.h"
-
+class THnSparse;
 
 class AliTPCcalibCosmic:public AliTPCcalibBase {
 public:
@@ -22,33 +19,56 @@ public:
   virtual ~AliTPCcalibCosmic();
   
   virtual void      Process(AliESDEvent *event);
-  virtual Long64_t  Merge(TCollection *li);
+  virtual Long64_t  Merge(TCollection *const li);
   virtual void      Analyze();
+  void              Add(const AliTPCcalibCosmic* cosmic);
   //
+  //
+  void              Init();
   void              FindPairs(AliESDEvent *event);
   Bool_t            IsPair(AliExternalTrackParam *tr0, AliExternalTrackParam *tr1);
-  void              SetGainMap(AliTPCCalPad *GainMap){fGainMap = GainMap;};
   static void       CalculateBetheParams(TH2F *hist, Double_t * initialParam);
   static Double_t   CalculateMIPvalue(TH1F * hist);
-  AliExternalTrackParam *Invert(AliExternalTrackParam *input);
+  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(const AliExternalTrackParam *par0, const AliExternalTrackParam *par1, const AliExternalTrackParam *inner0, const AliExternalTrackParam *inner1, AliTPCseed *seed0,  AliTPCseed *seed1, const AliExternalTrackParam *param0Combined);
+  void MaterialBudgetDump(AliExternalTrackParam *const par0, AliExternalTrackParam *const par1, const AliExternalTrackParam *inner0, const AliExternalTrackParam *inner1, AliTPCseed *const seed0,  AliTPCseed *const seed1, AliExternalTrackParam *const param0Combined, AliExternalTrackParam *const param1Combined);
+
+
   //
-  TH1F   *          GetHistNTracks(){return fHistNTracks;};
-  TH1F   *          GetHistClusters(){return fClusters;};
-  TH2F   *          GetHistAcorde(){return fModules;};
-  TH1F   *          GetHistPt(){return fHistPt;};
-  TH2F   *          GetHistDeDx(){return fDeDx;};
-  TH1F   *          GetHistMIP(){return fDeDxMIP;};
+  TH1F   *          GetHistNTracks() const {return fHistNTracks;};
+  TH1F   *          GetHistClusters() const {return fClusters;};
+  TH2F   *          GetHistAcorde()const {return fModules;};
+  TH1F   *          GetHistPt() const {return fHistPt;};
+  TH2F   *          GetHistDeDx() const {return fDeDx;};
+  TH1F   *          GetHistMIP() const {return fDeDxMIP;};
   //
-  Double_t          GetMIPvalue(){return fMIPvalue;};
+  Double_t          GetMIPvalue()const {return fMIPvalue;};
   //
-  static void       BinLogX(TH1 * h);   // method for correct histogram binning
+  static void       BinLogX(TH1 *const h);   // method for correct histogram binning
+  static void       BinLogX(THnSparse *const h, Int_t axisDim);   // method for correct histogram binning
 
+  void     Process(AliESDtrack *const track, Int_t runNo=-1) {AliTPCcalibBase::Process(track,runNo);};
+  void     Process(AliTPCseed *const track)  {return AliTPCcalibBase::Process(track);}
 
-private:
+public:  
+  //
+  // 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:
+  
   void              FillAcordeHist(AliESDtrack *upperTrack);
 
-  AliTPCCalPad *fGainMap;         //  gain map from Krypton calibration
+  
+
   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
@@ -68,7 +88,7 @@ private:
   AliTPCcalibCosmic(const AliTPCcalibCosmic&); 
   AliTPCcalibCosmic& operator=(const AliTPCcalibCosmic&); 
 
-  ClassDef(AliTPCcalibCosmic, 1); 
+  ClassDef(AliTPCcalibCosmic, 2); 
 };
 
 #endif