Update TPCCEda to write output file in parts (to avoid too big files produced in...
[u/mrichter/AliRoot.git] / TPC / AliTPCcalibTime.h
index 69c5691..1811a95 100644 (file)
@@ -5,15 +5,14 @@
  * See cxx source for full Copyright notice                               */
 
 #include "AliTPCcalibBase.h"
-#include "TH2F.h"
-#include "TF1.h"
-#include "TArrayD.h"
-#include "TObjArray.h"
-
+#include "THnSparse.h"           // Temporary
+#include "TH1D.h"                // Temporary make code compiling for HLT in the 
+class TObjArray;
 class TH1F;
 class TH3F;
 class TH2F;
 class THnSparse;
+class TH1D;
 class TList;
 class AliESDEvent;
 class AliESDtrack;
@@ -22,52 +21,66 @@ class TGraphErrors;
 class AliSplineFit;
 class AliESDfriendTrack;
 
-#include "TTreeStream.h"
-#include "TMap.h"
 class AliTPCcalibTime:public AliTPCcalibBase {
 public:
   AliTPCcalibTime(); 
-  AliTPCcalibTime(const Text_t *name, const Text_t *title, UInt_t StartTime, UInt_t EndTime, Int_t deltaIntegrationTimeVdrift);
+  AliTPCcalibTime(const Text_t *name, const Text_t *title, UInt_t StartTime, UInt_t EndTime, Int_t deltaIntegrationTimeVdrift, Int_t memoryMode=2);
   virtual ~AliTPCcalibTime();
   
   virtual void           Process(AliESDEvent *event);
-  virtual Long64_t       Merge(TCollection *li);
+  virtual Long64_t       Merge(TCollection *const li);
   virtual void           Analyze();
-  static Bool_t          IsLaser      (AliESDEvent *event);
-  static Bool_t          IsCosmics    (AliESDEvent *event);
-  static Bool_t          IsBeam       (AliESDEvent *event);
+  static Bool_t          IsLaser      (const AliESDEvent *const event);
+  static Bool_t          IsCosmics    (const AliESDEvent *const event);
+  static Bool_t          IsBeam       (const AliESDEvent *const event);
   void                   ProcessLaser (AliESDEvent *event);
-  void                   ProcessCosmic(AliESDEvent *event);
-  void                   ProcessBeam  (AliESDEvent *event);
+  void                   ProcessCosmic(const AliESDEvent *const event);
+  void                   ProcessBeam  (const AliESDEvent *const event);
   Bool_t                 IsPair(AliExternalTrackParam *tr0, AliExternalTrackParam *tr1);
-  Bool_t                 IsCross(AliESDtrack *tr0, AliESDtrack *tr1);
-  Bool_t                 IsSame (AliESDtrack *tr0, AliESDtrack *tr1);
-  void                   ProcessSame(AliESDtrack* track, AliESDfriendTrack *friendTrack,AliESDEvent *event);
-  void                   ProcessAlignITS(AliESDtrack* track, AliESDfriendTrack *friendTrack);
-  void                   ProcessAlignTRD(AliESDtrack* track, AliESDfriendTrack *friendTrack);
-  void                   ProcessAlignTOF(AliESDtrack* track, AliESDfriendTrack *friendTrack);
+  Bool_t                 IsCross(AliESDtrack *const tr0, AliESDtrack *const tr1);
+  Bool_t                 IsSame (AliESDtrack *const tr0, AliESDtrack *const tr1);
+  void                   ProcessSame(AliESDtrack *const track, AliESDfriendTrack *const friendTrack, const AliESDEvent *const event);
+  void                   ProcessAlignITS(AliESDtrack *const track, AliESDfriendTrack *const friendTrack, const AliESDEvent *const event, AliESDfriend *const ESDfriend);
+  void                   ProcessAlignTRD(AliESDtrack* const track, AliESDfriendTrack *const friendTrack);
+  void                   ProcessAlignTOF(AliESDtrack* const track, AliESDfriendTrack *const friendTrack);
 
-  THnSparse*    GetHistVdriftLaserA(Int_t index=1){return fHistVdriftLaserA[index];};
-  THnSparse*    GetHistVdriftLaserC(Int_t index=1){return fHistVdriftLaserC[index];};
-  THnSparse*    GetHistoDrift(const char* name);
-  TObjArray*    GetHistoDrift();
+  THnSparse*    GetHistVdriftLaserA(Int_t index=1) const {return fHistVdriftLaserA[index];};
+  THnSparse*    GetHistVdriftLaserC(Int_t index=1) const {return fHistVdriftLaserC[index];};
+  THnSparse*    GetHistoDrift(const char* name) const;
+  TObjArray*    GetHistoDrift() const;
   TGraphErrors* GetGraphDrift(const char* name);
   TObjArray*    GetGraphDrift();
   AliSplineFit* GetFitDrift(const char* name);
 //  TObjArray*    GetFitDrift();
-  TH1F*         GetCosmiMatchingHisto(Int_t index=0){return fCosmiMatchingHisto[index];};
-  
+  TH1F*         GetCosmiMatchingHisto(Int_t index=0) const {return fCosmiMatchingHisto[index];};
   void     Process(AliESDtrack *track, Int_t runNo=-1){AliTPCcalibBase::Process(track,runNo);};
   void     Process(AliTPCseed *track){return AliTPCcalibBase::Process(track);}
-  TObjArray* GetAlignITSTPC() {return fAlignITSTPC;}              // alignemnt array ITS TPC match
-  TObjArray* GetAlignTRDTPC() {return fAlignTRDTPC;}              // alignemnt array TRD TPC match
-  TObjArray* GetAlignTOFTPC() {return fAlignTOFTPC;}              // alignemnt array TOF TPC match
+  TObjArray* GetAlignITSTPC() const {return fAlignITSTPC;}              // alignemnt array ITS TPC match
+  TObjArray* GetAlignTRDTPC() const {return fAlignTRDTPC;}              // alignemnt array TRD TPC match
+  TObjArray* GetAlignTOFTPC() const {return fAlignTOFTPC;}              // alignemnt array TOF TPC match
 
+  THnSparse * GetTPCVertexHisto(Int_t index) { return fTPCVertex[index%12];}
+  THnSparse * GetTPCVertexHistoCorrelation(Int_t index) { return fTPCVertexCorrelation[index%5];}
 
-private:
-  void ResetCurrent();                  // reset current values
+  THnSparse*  GetResHistoTPCCE(Int_t index) const { return (index<5) ? fResHistoTPCCE[index]:0;}        //TPC-CE    matching map
+  THnSparse*  GetResHistoTPCITS(Int_t index) const { return (index<5) ? fResHistoTPCITS[index]:0;}        //TPC-ITS    matching map
+  THnSparse*  GetResHistoTPCvertex(Int_t index)      const { return (index<5) ? fResHistoTPCvertex[index]   :0;}        //TPC vertex matching map
+  THnSparse*  GetResHistoTPCTRD(Int_t index)   const { return (index<5) ? fResHistoTPCTRD[index]:0;}        //TPC-TRD    matching map
+  THnSparse*  GetResHistoTPCTOF(Int_t index)   const { return (index<5) ? fResHistoTPCTOF[index]:0;}        //TPC-TOF    matching map
 
+  void        BookDistortionMaps();      // book histograms
+  void        FillResHistoTPCCE(const AliExternalTrackParam * pTPCIn, const AliExternalTrackParam * pTPCOut );       // fill residual histo
+  void        FillResHistoTPCITS(const AliExternalTrackParam * pTPCIn, const AliExternalTrackParam * pITSOut );       // fill residual histo
+  void        FillResHistoTPC(const AliESDtrack * pTrack);
+  void        FillResHistoTPCTRD(const AliExternalTrackParam * pTPCOut, const AliExternalTrackParam * pTRDIn );
+  void        FillResHistoTPCTOF(const AliExternalTrackParam * pTPCOut, const AliExternalTrackParam * pTOFIn );
+
+  TObjArray * GetLaserArrayA() const { return fArrayLaserA;}
+  TObjArray * GetLaserArrayC() const { return fArrayLaserC;}
+
+protected:
+  void ResetCurrent();                  // reset current values
+  Int_t              fMemoryMode;       // 0 -do not fill THnSparse with residuals  1- fill only important QA THn 2 - Fill all THnsparse for calibration
   AliTPCcalibLaser * fLaser;            //! laser calibration
   //
   // current information
@@ -82,20 +95,32 @@ private:
   Float_t fCutMinDir;   // direction vector products
   Int_t   fCutTracks;   // maximal number of tracks
  
-  AliTPCcalibTime(const AliTPCcalibTime&); 
-  AliTPCcalibTime& operator=(const AliTPCcalibTime&); 
 
   TH1F* fCosmiMatchingHisto[10];
-
+  //
+  // distortion maps
+  //
+  THnSparse*  fResHistoTPCCE[5];        //TPC-TPCE matching map
+  THnSparse*  fResHistoTPCITS[5];        //TPC-ITS    matching map
+  THnSparse*  fResHistoTPCvertex[5];           //TPC-ITS    vertex matching map
+  THnSparse*  fResHistoTPCTRD[5];        //TPC-TRD    matching map
+  THnSparse*  fResHistoTPCTOF[5];        //TPC-TRD    matching map
   // laser histo
   THnSparse * fHistVdriftLaserA[3];    //Histograms for V drift from laser
   THnSparse * fHistVdriftLaserC[3];    //Histograms for V drift from laser
+  TObjArray *fArrayLaserA;              //Object array of driftvelocity laserA
+  TObjArray *fArrayLaserC;              //Object array of driftvelocity laserC
+  //
+  // TPC vertex A side C side histo
+  //
+  THnSparse * fTPCVertex[12];           // TPC vertex histograms A side c side - A+C -ESD
+  THnSparse * fTPCVertexCorrelation[5];       // TPC vertex correlation A side C side with TPC vertex and ITS vertex     
   // DELTA Z histo
-//  TMap*      fMapDz;                 //Tmap of V drifts for different triggers
   TObjArray* fArrayDz;                  // array of DZ histograms for different triggers
   TObjArray* fAlignITSTPC;              // alignemnt array ITS TPC match
   TObjArray* fAlignTRDTPC;              // alignemnt array TRD TPC match
   TObjArray* fAlignTOFTPC;              // alignemnt array TOF TPC match
+  Int_t      fTimeKalmanBin;            // width of Kalman bin - time in seconds
   Int_t    fTimeBins;                  //Bins time
   Double_t fTimeStart;                 //Start time
   Double_t fTimeEnd;                   //End time
@@ -111,7 +136,12 @@ private:
   Int_t    fBinsVdrift[4];             //Bins for vdrift
   Double_t fXminVdrift[4];             //Xmax for vdrift
   Double_t fXmaxVdrift[4];             //Xmin for vdrift
-  ClassDef(AliTPCcalibTime, 2); 
+
+private:
+  AliTPCcalibTime(const AliTPCcalibTime&); 
+  AliTPCcalibTime& operator=(const AliTPCcalibTime&); 
+
+  ClassDef(AliTPCcalibTime, 10); 
 };
 
 #endif