/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
* See cxx source for full Copyright notice */
-
+
#include <TChain.h>
#include <TNamed.h>
#include <iostream>
#include <TH1F.h>
+#include <AliTPCcalibBase.h>
+
using namespace std;
class TTreeSRedirector;
class TGraph;
class AliTPCCClusterParam;
-class AliTPCcalibTracksGain : public TNamed {
+class AliTPCcalibTracksGain : public AliTPCcalibBase {
public:
enum {
kShortPads = 0,
AliTPCcalibTracksGain();
AliTPCcalibTracksGain(const AliTPCcalibTracksGain& obj);
- AliTPCcalibTracksGain(const char* name, const char* title, AliTPCcalibTracksCuts* cuts, TNamed* debugStreamPrefix = 0, AliTPCcalibTracksGain* prevIter = 0);
+ AliTPCcalibTracksGain(const char* name, const char* title, AliTPCcalibTracksCuts* cuts);
virtual ~AliTPCcalibTracksGain();
AliTPCcalibTracksGain& operator=(const AliTPCcalibTracksGain& rhs);
- static void AddInfo(TChain* chain, char* debugStreamPrefix = 0, char* prevIterFileName = 0);
- void Terminate();
+ virtual Long64_t Merge(TCollection *list);
+ virtual void Process(AliTPCseed* seed);
+ virtual void Terminate();
+ virtual void Analyze();
//
- // Tracks manipulation
+ // Tracks and cluster manipulation
//
- void Process(AliTPCseed* seed);
+ Float_t GetGain(AliTPCclusterMI *cl);
+ Float_t GetMaxNorm(AliTPCclusterMI * cl, Float_t ky, Float_t kz);
+ Float_t GetQNorm(AliTPCclusterMI * cl, Float_t ky, Float_t kz);
void AddTrack(AliTPCseed* seed);
- Bool_t AcceptTrack(AliTPCseed* track);
void DumpTrack(AliTPCseed* track);
Bool_t GetDedx(AliTPCseed* track, Int_t padType, Int_t*rows,
Int_t §or, Int_t& npoints,
void AddCluster(AliTPCclusterMI* cluster, Float_t momenta, Float_t mdedx, Int_t padType, Float_t xcenter, TVectorD &dedxQ, TVectorD & dedxM, Float_t fraction, Float_t fraction2, Float_t dedge, TVectorD& parY, TVectorD& parZ, TVectorD& meanPos);
void AddTracklet(UInt_t sector, UInt_t padType,TVectorD &dedxQ, TVectorD &dedxM,TVectorD& parY, TVectorD& parZ, TVectorD& meanPos);
- void AddCluster(AliTPCclusterMI* cluster);
-
- //
- // Merging of the component
- //
- Long64_t Merge(TCollection *list);
void Add(AliTPCcalibTracksGain* cal);
- //
- // Histogram part
- //
- TH1F * GetQM(Int_t sector=-1){return (TH1F*)(sector<0 ? fArrayQM->At(72): fArrayQM->At(sector));}
- TH1F * GetQT(Int_t sector=-1){return (TH1F*)(sector<0 ? fArrayQT->At(72): fArrayQT->At(sector));}
- TProfile* GetProfileQM(Int_t sector){return (TProfile*)(sector<0 ? fProfileArrayQM->At(36): fProfileArrayQM->At(sector));}
- TProfile* GetProfileQT(Int_t sector){return (TProfile*)(sector<0 ? fProfileArrayQT->At(36): fProfileArrayQT->At(sector));}
- TProfile2D* GetProfileQM2D(Int_t sector){return (TProfile2D*)(sector<0 ? fProfileArrayQM2D->At(36): fProfileArrayQM2D->At(sector));}
- TProfile2D* GetProfileQT2D(Int_t sector){return (TProfile2D*)(sector<0 ? fProfileArrayQT2D->At(36): fProfileArrayQT2D->At(sector));}
+
//
// Get Derived results - gain maps
//
AliTPCCalROC* CreateCombinedCalROC(const AliTPCCalROC* roc1, const AliTPCCalROC* roc2);
//
void Evaluate(Bool_t robust = kFALSE, Double_t frac = -1.);
- void GetParameters(UInt_t segment, UInt_t padType, UInt_t fitType, TVectorD &fitParam);
+ Bool_t GetParameters(UInt_t segment, UInt_t padType, UInt_t fitType, TVectorD &fitParam);
void GetErrors(UInt_t segment, UInt_t padType, UInt_t fitType, TVectorD &fitError);
Double_t GetRedChi2(UInt_t segment, UInt_t padType, UInt_t fitType);
void GetCovarianceMatrix(UInt_t segment, UInt_t padType, UInt_t fitType, TMatrixD& covMatrix);
TGraph * CreateAmpGraph(Int_t ipad, Bool_t qmax);
- TLinearFitter* GetFitter(UInt_t segment, UInt_t padType, UInt_t fitType);
+ TLinearFitter* GetFitter(UInt_t segment, UInt_t padType, UInt_t fitType);
+ void Process(AliESDEvent *event) {AliTPCcalibBase::Process(event);};
+ void Process(AliESDtrack *track, Int_t runNo=-1){AliTPCcalibBase::Process(track,runNo);};
+
public:
//
// Helper function
//
static Double_t GetPadLength(Double_t lx);
static Int_t GetPadType(Double_t lx);
- static Bool_t GetRowPad(Double_t lx, Double_t ly, Int_t& row, Int_t& pad); // just for debugging
//
//
- TTreeSRedirector* fDebugStream; //! debug stream for debugging
AliTPCcalibTracksCuts* fCuts; // cuts that are used for sieving the tracks used for calibration
//
- // Simple Profiles and histograms - per chambers + 1 total
+ // kalman fit - alignment of dEdx
//
- TObjArray* fArrayQM; // Qmax normalized
- TObjArray* fArrayQT; // Qtot normalized
- TObjArray* fProfileArrayQM; // Qmax normalized versus local X
- TObjArray* fProfileArrayQT; // Qtot normalized versus local X
- TObjArray* fProfileArrayQM2D; // Qmax normalized versus local X and phi
- TObjArray* fProfileArrayQT2D; // Qtot normalized versus local X and phi
//
// Fitters
//
AliTPCFitPad* fSimpleFitter; // simple fitter for short pads
AliTPCFitPad* fSqrtFitter; // sqrt fitter for medium pads
AliTPCFitPad* fLogFitter; // log fitter for long pads
+ //
TLinearFitter* fFitter0M; // fitting of the atenuation, angular correction, and mean chamber gain
TLinearFitter* fFitter1M; // fitting of the atenuation, angular correction, and mean chamber gain
TLinearFitter* fFitter2M; // fitting of the atenuation, angular correction, and mean chamber gain
TLinearFitter* fFitter0T; // fitting of the atenuation, angular correction, and mean chamber gain
TLinearFitter* fFitter1T; // fitting of the atenuation, angular correction, and mean chamber gain
TLinearFitter* fFitter2T; // fitting of the atenuation, angular correction, and mean chamber gain
- AliTPCFitPad* fSingleSectorFitter; // just for debugging
+ //
+ // angular adn diffusion effect fitter
+ //
+ TLinearFitter* fDFitter0M; // fitting of the atenuation, angular correction
+ TLinearFitter* fDFitter1M; // fitting of the atenuation, angular correction
+ TLinearFitter* fDFitter2M; // fitting of the atenuation, angular correction
+ TLinearFitter* fDFitter0T; // fitting of the atenuation, angular correction
+ TLinearFitter* fDFitter1T; // fitting of the atenuation, angular correction
+ TLinearFitter* fDFitter2T; // fitting of the atenuation, angular correction
+ //
+ AliTPCFitPad* fSingleSectorFitter; // just for debugging
//
// Conters
//
UInt_t fTotalTracks; // just for debugging
UInt_t fAcceptedTracks; // just for debugging
- AliTPCCalPad* fDebugCalPadRaw; // just for debugging
- AliTPCCalPad* fDebugCalPadCorr; // just for debugging
UInt_t fNShortClusters[36]; // number of clusters registered on short pads
UInt_t fNMediumClusters[36]; // number of clusters registered on medium pads
UInt_t fNLongClusters[36]; // number of clusters registered on long pads
//
//
- AliTPCcalibTracksGain* fPrevIter; // the calibration object in its previous iteration (will not be owned by the new object, don't forget to delete it!)
//
// Setup
//