/* 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();
+ void SetGainMap(AliTPCCalPad *GainMap){fGainMap = GainMap;};
//
- // Tracks manipulation
+ // Tracks and cluster manipulation
//
- void Process(AliTPCseed* seed);
+ Float_t GetMaxNorm(AliTPCclusterMI * cl);
+ Float_t GetQNorm(AliTPCclusterMI * cl);
+ Float_t GetGain(AliTPCclusterMI* cl);
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 AddTracklet(UInt_t sector, UInt_t padType,TVectorD &dedxQ, TVectorD &dedxM,TVectorD& parY, TVectorD& parZ, TVectorD& meanPos);
void AddCluster(AliTPCclusterMI* cluster);
+ void Add(AliTPCcalibTracksGain* cal);
//
- // Merging of the component
+ // Debug stream analyze part
//
- Long64_t Merge(TCollection *list);
- void Add(AliTPCcalibTracksGain* cal);
+ static TVectorD * MakeQPosNorm(TTree * chain, Int_t ipad, Bool_t isMax, Int_t maxPoints=1000000, Int_t verbose=0);
+
+ static void MakeQPosNormAll(TTree * chain, AliTPCClusterParam * param, Int_t maxPoints=1000000, Int_t verbose=0);
+
//
// Histogram part
//
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 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
+ AliTPCCalPad *fGainMap; // gain map to be applied
+ //
//
// Simple Profiles and histograms - per chambers + 1 total
//
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
+ //
+ // 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 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
//