]> git.uio.no Git - u/mrichter/AliRoot.git/blob - TPC/AliTPCcalibTime.h
adding merging procedure for distortion maps
[u/mrichter/AliRoot.git] / TPC / AliTPCcalibTime.h
1 #ifndef ALITPCCALIBTIME_H
2 #define ALITPCCALIBTIME_H
3
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5  * See cxx source for full Copyright notice                               */
6
7 #include "AliTPCcalibBase.h"
8 class TObjArray;
9
10 class TH1F;
11 class TH3F;
12 class TH2F;
13 class THnSparse;
14 class TList;
15 class AliESDEvent;
16 class AliESDtrack;
17 class AliTPCcalibLaser;
18 class TGraphErrors;
19 class AliSplineFit;
20 class AliESDfriendTrack;
21
22 class AliTPCcalibTime:public AliTPCcalibBase {
23 public:
24   AliTPCcalibTime(); 
25   AliTPCcalibTime(const Text_t *name, const Text_t *title, UInt_t StartTime, UInt_t EndTime, Int_t deltaIntegrationTimeVdrift);
26   virtual ~AliTPCcalibTime();
27   
28   virtual void           Process(AliESDEvent *event);
29   virtual Long64_t       Merge(TCollection *const li);
30   virtual void           Analyze();
31   static Bool_t          IsLaser      (const AliESDEvent *const event);
32   static Bool_t          IsCosmics    (const AliESDEvent *const event);
33   static Bool_t          IsBeam       (const AliESDEvent *const event);
34   void                   ProcessLaser (AliESDEvent *event);
35   void                   ProcessCosmic(const AliESDEvent *const event);
36   void                   ProcessBeam  (const AliESDEvent *const event);
37   Bool_t                 IsPair(AliExternalTrackParam *tr0, AliExternalTrackParam *tr1);
38   Bool_t                 IsCross(AliESDtrack *const tr0, AliESDtrack *const tr1);
39   Bool_t                 IsSame (AliESDtrack *const tr0, AliESDtrack *const tr1);
40   void                   ProcessSame(AliESDtrack *const track, AliESDfriendTrack *const friendTrack, const AliESDEvent *const event);
41   void                   ProcessAlignITS(AliESDtrack *const track, AliESDfriendTrack *const friendTrack, const AliESDEvent *const event, AliESDfriend *const ESDfriend);
42   void                   ProcessAlignTRD(AliESDtrack* const track, AliESDfriendTrack *const friendTrack);
43   void                   ProcessAlignTOF(AliESDtrack* const track, AliESDfriendTrack *const friendTrack);
44
45   THnSparse*    GetHistVdriftLaserA(Int_t index=1) const {return fHistVdriftLaserA[index];};
46   THnSparse*    GetHistVdriftLaserC(Int_t index=1) const {return fHistVdriftLaserC[index];};
47   THnSparse*    GetHistoDrift(const char* name) const;
48   TObjArray*    GetHistoDrift() const;
49   TGraphErrors* GetGraphDrift(const char* name);
50   TObjArray*    GetGraphDrift();
51   AliSplineFit* GetFitDrift(const char* name);
52 //  TObjArray*    GetFitDrift();
53   TH1F*         GetCosmiMatchingHisto(Int_t index=0) const {return fCosmiMatchingHisto[index];};
54   
55   void     Process(AliESDtrack *track, Int_t runNo=-1){AliTPCcalibBase::Process(track,runNo);};
56   void     Process(AliTPCseed *track){return AliTPCcalibBase::Process(track);}
57   TObjArray* GetAlignITSTPC() const {return fAlignITSTPC;}              // alignemnt array ITS TPC match
58   TObjArray* GetAlignTRDTPC() const {return fAlignTRDTPC;}              // alignemnt array TRD TPC match
59   TObjArray* GetAlignTOFTPC() const {return fAlignTOFTPC;}              // alignemnt array TOF TPC match
60
61   THnSparse*  GetResHistoTPCITS(Int_t index) const { return (index<5) ? fResHistoTPCITS[index]:0;}        //TPC-ITS    matching map
62   THnSparse*  GetResHistoTPCvertex(Int_t index)      const { return (index<5) ? fResHistoTPCvertex[index]   :0;}        //TPC vertex matching map
63   THnSparse*  GetResHistoTPCTRD(Int_t index)   const { return (index<5) ? fResHistoTPCTRD[index]:0;}        //TPC-TRD    matching map
64
65   void        BookDistortionMaps();      // book histograms
66   void        FillResHistoTPCITS(const AliExternalTrackParam * pTPCIn, const AliExternalTrackParam * pITSOut );       // fill residual histo
67   void        FillResHistoTPC(const AliESDtrack * pTrack);
68   void        FillResHistoTPCTRD(const AliExternalTrackParam * pTPCOut, const AliExternalTrackParam * pTRDIn );
69
70 private:
71   void ResetCurrent();                  // reset current values
72
73   AliTPCcalibLaser * fLaser;            //! laser calibration
74   //
75   // current information
76   //
77   Float_t fDz;          //! current delta z
78   
79   // cuts
80   //
81   Float_t fCutMaxD;     // maximal distance in rfi ditection
82   Float_t fCutMaxDz;    // maximal distance in z ditection
83   Float_t fCutTheta;    // maximal distance in theta ditection
84   Float_t fCutMinDir;   // direction vector products
85   Int_t   fCutTracks;   // maximal number of tracks
86  
87   AliTPCcalibTime(const AliTPCcalibTime&); 
88   AliTPCcalibTime& operator=(const AliTPCcalibTime&); 
89
90   TH1F* fCosmiMatchingHisto[10];
91   //
92   // distortion maps
93   //
94   THnSparse*  fResHistoTPCITS[5];        //TPC-ITS    matching map
95   THnSparse*  fResHistoTPCvertex[5];           //TPC-ITS    vertex matching map
96   THnSparse*  fResHistoTPCTRD[5];        //TPC-TRD    matching map
97   // laser histo
98   THnSparse * fHistVdriftLaserA[3];     //Histograms for V drift from laser
99   THnSparse * fHistVdriftLaserC[3];     //Histograms for V drift from laser
100   // DELTA Z histo
101   TObjArray* fArrayDz;                  // array of DZ histograms for different triggers
102   TObjArray* fAlignITSTPC;              // alignemnt array ITS TPC match
103   TObjArray* fAlignTRDTPC;              // alignemnt array TRD TPC match
104   TObjArray* fAlignTOFTPC;              // alignemnt array TOF TPC match
105   Int_t    fTimeBins;                   //Bins time
106   Double_t fTimeStart;                  //Start time
107   Double_t fTimeEnd;                    //End time
108   Int_t    fPtBins;                     //Bins pt
109   Double_t fPtStart;                    //Start pt
110   Double_t fPtEnd;                      //End pt
111   Int_t    fVdriftBins;                 //Bins vdrift
112   Double_t fVdriftStart;                //Start vdrift
113   Double_t fVdriftEnd;                  //End vdrift
114   Int_t    fRunBins;                    //Bins run
115   Double_t fRunStart;                   //Start run
116   Double_t fRunEnd;                     //End run
117   Int_t    fBinsVdrift[4];              //Bins for vdrift
118   Double_t fXminVdrift[4];              //Xmax for vdrift
119   Double_t fXmaxVdrift[4];              //Xmin for vdrift
120   ClassDef(AliTPCcalibTime, 3); 
121 };
122
123 #endif
124
125