2 #ifndef ALIPERFANALYZEINVPT_H
3 #define ALIPERFANALYZEINVPT_H
5 //----------------------------------------------------------------------------------------------------
6 // Class to analyse output of AliPerformancePtCalib.cxx and AliPerformancePtCalibMC.cxx:
7 // Projection of charge/pt vs theta and vs phi resp. histoprams will be fitted with either
8 // polynomial or gaussian fit function to extract minimum position of charge/pt.
9 // Fit options and theta, phi bins can be set by user.
10 // Attention: use the Set* functions of AliPerformancePtCalib.h and AliPerformancePtCalibMC.h
11 // when running AliPerformancePtCalib*::Analyse().
13 // Author: S. Schuchmann 11/13/2009
14 //----------------------------------------------------------------------------------------------------
29 class AliPerfAnalyzeInvPt : public TNamed {
32 AliPerfAnalyzeInvPt(const Char_t* name="AliPerfAnalyzeInvPt", const Char_t* title="AliPerfAnalyzeInvPt");
33 virtual ~AliPerfAnalyzeInvPt(){;}
34 void InitGraphs(Double_t *binsXTheta,Double_t *fitParamTheta,Double_t *errFitParamTheta,Double_t *binsXPhi,Double_t *fitParamPhi,Double_t *errFitParamPhi);
36 void StartAnalysis(const TH2F *histThetaInvPt, const TH2F *histPhiInvPt, TObjArray* aFolderObj);
37 void SetProjBinsPhi(const Double_t *pBins,const Int_t sizep);
38 void SetProjBinsTheta(const Double_t *tBins, const Int_t sizet);
39 void SetMakeFitOption(const Bool_t setGausFit, const Double_t exclusionR,const Double_t fitR );
40 void SetDoRebin(const Int_t rebin){if(rebin) {fDoRebin = kTRUE; fRebin = rebin;}}
43 Double_t fThetaBins[100];// array of theta bins for projection and fit
44 Double_t fPhiBins[100];// array of phi bins for projection and fit
46 Int_t fNThetaBins;// number of theta bins for projection and fit
47 Int_t fNPhiBins; // number of phi bins for projection and fit
48 Double_t fRange;// fit range of 1/pt spectrum
49 Double_t fExclRange ;// range of exlusion of points around 0 when fitting 1/pt
50 Bool_t fFitGaus;// set this flag for usage of gaussian fit function instead of polynomial (default)
54 TH1D *fHistFitTheta[100];// projection histos for analysis in theta bins
55 TH1D *fHistFitPhi[100];// projection histos for analysis in phi bins
57 //histos for copy of input histos
58 TH2F *fHistH2InvPtTheta; // for copy of input histos
59 TH2F *fHistH2InvPtPhi; // for copy of input histos
61 //graphs for displaying the fit parameters (minimum position of 1/pt and error)
62 TGraphErrors *fGrMinPosTheta;
63 TGraphErrors *fGrMinPosPhi;
69 TF1 *fFitMinPos;//fit function for first fit with polynomial
70 TF1 *fFitMinPosRejP;//fit function for second fit with polynomial to improve result
71 TF1 *fFitInvGauss;//fit functionfor first fit with gaussian
72 TF1 *fFitInvGaussRejP;//fit function for second fit with gaussian to improve result
74 // infput for fit functions
75 static Double_t Polynomial(Double_t *x, const Double_t *par);
76 static Double_t PolynomialRejP(Double_t *x, const Double_t *par);
77 static Double_t InvGauss(Double_t *x, const Double_t *par);
78 static Double_t InvGaussRejP(Double_t *x, const Double_t *par);
80 //functions for fitting procedure
81 void MakeFit(TH1D *dmproy, TF1 * fitpb, Double_t &mean, Double_t &ErrMean, Double_t &excl,Double_t &range);
82 void MakeFitBetter(TH1D *dmproy, TF1 * fitpb2, Double_t &mean, Double_t &ErrMean, Double_t &f, Double_t &excl, Double_t &range);
83 void MakeFitInvGauss(TH1D *dmproy, TF1 * fitpb2, Double_t &mean, Double_t &ErrMean,Double_t &excl , Double_t &range);
84 void MakeFitInvGaussBetter(TH1D *dmproy, TF1 * fitpb2, Double_t &mean, Double_t &ErrMean, Double_t &f,Double_t &excl, Double_t &range);
89 AliPerfAnalyzeInvPt(const AliPerfAnalyzeInvPt&); // not implemented
90 AliPerfAnalyzeInvPt& operator=(const AliPerfAnalyzeInvPt&); // not implemented
92 ClassDef(AliPerfAnalyzeInvPt, 1);