#ifndef ALITRDRESOLUTION_H
#define ALITRDRESOLUTION_H
#ifndef ALITRDRECOTASK_H
#include "AliTRDrecoTask.h"
#endif
#ifndef Root_TNamed
#include "TNamed.h"
#endif
class TH1;
class TH2;
class TF1;
class TGraphErrors;
class TObjArray;
class TDatabasePDG;
class AliTRDrecoParam;
class AliTRDseedV1;
class AliTRDtrackInfo;
class AliTrackPoint;
class AliTRDresolution : public AliTRDrecoTask
{
public:
enum ETRDresolutionSteer {
kVerbose = BIT(18)
,kVisual = BIT(19)
,kTrackRefit = BIT(20)
,kTrackSelect= BIT(21)
,kXchange = BIT(22)
};
enum ETRDresolutionSlots {
kClToTrk = 2
,kClToMC
,kTrkltToTrk
,kTrkltToMC
,kNSlots = 4
};
enum ETRDresolutionClass {
kDetector=0
,kCluster
,kTracklet
,kTrackIn
,kMCcluster
,kMCtracklet
,kMCtrackIn
,kMCtrack
,kNclasses
,kV0TrackIn
};
enum ETRDresolutionClassProjs {
kDetNproj=460
,kClNproj=2800
,kTrkltNproj=8000
,kTrkInNproj=1440
,kTrkNproj=2500
,kMCTrkInNproj=1260
};
enum ETRDresolutionProjs {
kBC = 0
,kPhi
,kEta
,kYrez
,kPrez
,kZrez
,kSpeciesChgRC
,kPt
,kNdim
,kNdimDet = 4
,kNdimCl = 4
,kNdimTrklt = 4
,kNdimTrkIn = 7
,kNbunchCross = 3
,kNspc = 3
,kNcharge = 2
,kNpads = 4
};
AliTRDresolution();
AliTRDresolution(char* name, Bool_t xchange=kTRUE);
AliTRDresolution(const AliTRDresolution&);
AliTRDresolution& operator=(const AliTRDresolution&);
virtual ~AliTRDresolution();
static Bool_t FitTrack(const Int_t np, AliTrackPoint *points, Float_t params[10]);
static Bool_t FitTracklet(const Int_t ly, const Int_t np, const AliTrackPoint *points, const Float_t trackPars[10], Float_t trackletPars[3]);
void UserCreateOutputObjects();
static void GetRangeZ(TH2 *h2, Float_t &m, Float_t &M);
Float_t GetPtThreshold() const {return fPtThreshold;}
Bool_t GetRefFigure(Int_t ifig);
virtual TObjArray* Histos();
void MakeSummary();
TObjArray* Results(ETRDresolutionClass c) const { if(!fProj) return NULL; return (TObjArray*)fProj->At(c);}
void UserExec(Option_t * opt);
void InitExchangeContainers();
Bool_t IsVerbose() const { return TestBit(kVerbose);}
Bool_t IsVisual() const { return TestBit(kVisual);}
Bool_t UseBCselectTOF() const { return fBCbinTOF>0;}
Bool_t UseBCselectFill() const { return fBCbinFill>0;}
Bool_t UseLYselectTrklt() const { return fLYselect>=0;}
Bool_t UseExchangeContainers() const { return TestBit(kXchange);}
Bool_t PostProcess();
TH1* DetCluster(const TObjArray *cl=NULL);
TH1* PlotCluster(const AliTRDtrackV1 *t=NULL);
TH1* PlotTracklet(const AliTRDtrackV1 *t=NULL);
TH1* PlotTrackIn(const AliTRDtrackV1 *t=NULL);
TH1* PlotMC(const AliTRDtrackV1 *t=NULL);
static Bool_t Process(TH2* const , TGraphErrors **, Int_t stat=100){ return Bool_t(stat);}
void SetPtThreshold(Float_t pt) { fPtThreshold = pt;}
void SetBCselectTOF(Int_t b=0) { fBCbinTOF = b;}
void SetBCselectFill(Int_t b=0) { fBCbinFill = b<0||b>3499?1:b+1;}
void SetBsign(Int_t b=0) { fBsign = Bool_t(b);}
void SetLYselect(Int_t ly=0) { fLYselect = ly;}
void SetProcesses(Bool_t det, Bool_t cl, Bool_t trklt, Bool_t trkin);
void SetDump3D(Bool_t det, Bool_t cl, Bool_t trklt, Bool_t trkin);
void SetVerbose(Bool_t v = kTRUE) { SetBit(kVerbose, v);}
void SetVisual(Bool_t v = kTRUE) { SetBit(kVisual, v);}
void SetUseExchangeContainers(Bool_t v = kTRUE) { SetBit(kXchange, v);}
void Terminate(Option_t * opt);
static Bool_t UseTrack(const Int_t np, const AliTrackPoint *points, Float_t params[10]);
void AdjustF1(TH1 *h, TF1 *f);
TH2* DrawSigma(TH2 *h2, const Char_t *t, Float_t m=0., Float_t M=-1., Float_t scale=1);
void GetLandauMpvFwhm(TF1 * const f, Float_t &mpv, Float_t &xm, Float_t &xM);
void GetRange(TH2 *h2, Char_t mod, Float_t *range);
protected:
Bool_t HasDump3DFor(ETRDresolutionClass cls) const { return TESTBIT(fSteer, 4+cls);}
Bool_t HasProcess(ETRDresolutionClass cls) const { return TESTBIT(fSteer, cls);}
Bool_t MakeProjectionDetector();
Bool_t MakeProjectionCluster(Bool_t mc=kFALSE);
Bool_t MakeProjectionTracklet(Bool_t mc=kFALSE);
Bool_t MakeProjectionTrackIn(Bool_t mc=kFALSE, Bool_t v0=kFALSE);
Bool_t MakeProjectionTrack();
Bool_t Process(TH2* const , TF1 *, Float_t , TGraphErrors **) { return kTRUE;}
Bool_t Pulls(Double_t dyz[2], Double_t cc[3], Double_t tilt) const;
UShort_t fSteer;
Float_t fPtThreshold;
Int_t fBCbinTOF;
Int_t fBCbinFill;
Int_t fLYselect;
Bool_t fBsign;
static Char_t const *fgPerformanceName[kNclasses];
static Int_t const fgkNbins[kNdim];
static Double_t const fgkMin[kNdim];
static Double_t const fgkMax[kNdim];
static Char_t const *fgkTitle[kNdim];
TObjArray *fProj;
TDatabasePDG *fDBPDG;
TObjArray *fCl;
TObjArray *fMCcl;
ClassDef(AliTRDresolution, 11)
};
#endif