* 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;
class AliTPCcalibLaser;
+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 deltaIntegrationTimeDeDx, 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 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 ProcessCosmic(const AliESDEvent *const event);
+ void ProcessBeam (const AliESDEvent *const event);
Bool_t IsPair(AliExternalTrackParam *tr0, AliExternalTrackParam *tr1);
- //
- THnSparse * GetHistVdrift(){return (THnSparse*) fHistVdrift;};
- THnSparse * GetHistDeDxVsTgl(){return (THnSparse*) fHistDeDxTgl;};
- THnSparse * GetHistDeDx(){return (THnSparse*) fHistDeDx;};
- THnSparse * GetHistVdriftLaserA(Int_t index=1){return (THnSparse*) fHistVdriftLaserA[index];};
- THnSparse * GetHistVdriftLaserC(Int_t index=1){return (THnSparse*) fHistVdriftLaserC[index];};
- TMap * GetMapDz(){return (TMap *) fMapDz;};
+ 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) 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) const {return fCosmiMatchingHisto[index];};
void Process(AliESDtrack *track, Int_t runNo=-1){AliTPCcalibBase::Process(track,runNo);};
void Process(AliTPCseed *track){return AliTPCcalibBase::Process(track);}
-private:
- void ResetCurrent(); // reset current values
+ 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];}
+
+ 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
- THnSparse * fHistDeDxTgl; // dEdx vs. dip angle vs time histogram
- THnSparse * fHistDeDx; // dEdx vs. time histogram (cosmics: all particles on Fermi plateau)
- THnSparse * fHistVdrift; // drift velocity vs time histogram
+ 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 );
- Float_t fIntegrationTimeDeDx; // required statistics for each dEdx time bin
- Float_t fIntegrationTimeVdrift; // required statistics for each Vdrift time bin
+ 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
//
Float_t fDz; //! current delta z
- Float_t fdEdx; //! current dEdx
- Float_t fdEdxRatio; //! current dEdx ratio
- Float_t fTl; //! current tan(lambda)
// cuts
//
Float_t fCutMaxD; // maximal distance in rfi ditection
- Float_t fCutMaxDz; // maximal distance in z ditection
+ Float_t fCutMaxDz; // maximal distance in z ditection
Float_t fCutTheta; // maximal distance in theta ditection
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]; //NEW! Histograms for V drift from laser
- THnSparse * fHistVdriftLaserC[3]; //NEW! Histograms for V drift from laser
+ 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; //NEW! 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
+ Int_t fPtBins; //Bins pt
+ Double_t fPtStart; //Start pt
+ Double_t fPtEnd; //End pt
+ Int_t fVdriftBins; //Bins vdrift
+ Double_t fVdriftStart; //Start vdrift
+ Double_t fVdriftEnd; //End vdrift
+ Int_t fRunBins; //Bins run
+ Double_t fRunStart; //Start run
+ Double_t fRunEnd; //End run
+ Int_t fBinsVdrift[4]; //Bins for vdrift
+ Double_t fXminVdrift[4]; //Xmax for vdrift
+ Double_t fXmaxVdrift[4]; //Xmin for vdrift
+
+private:
+ AliTPCcalibTime(const AliTPCcalibTime&);
+ AliTPCcalibTime& operator=(const AliTPCcalibTime&);
- Int_t fTimeBins; //NEW! Bins time
- Double_t fTimeStart; //NEW! Start time
- Double_t fTimeEnd; //NEW! End time
- Int_t fPtBins; //NEW! Bins pt
- Double_t fPtStart; //NEW! Start pt
- Double_t fPtEnd; //NEW! End pt
- Int_t fVdriftBins; //NEW! Bins vdrift
- Double_t fVdriftStart; //NEW! Start vdrift
- Double_t fVdriftEnd; //NEW! End vdrift
- Int_t fRunBins; //NEW! Bins run
- Double_t fRunStart; //NEW! Start run
- Double_t fRunEnd; //NEW! End run
- Int_t fBinsVdrift[4]; //NEW! Bins for vdrift
- Double_t fXminVdrift[4]; //NEW! Xmax for vdrift
- Double_t fXmaxVdrift[4]; //NEW! Xmin for vdrift
- ClassDef(AliTPCcalibTime, 1);
+ ClassDef(AliTPCcalibTime, 10);
};
#endif