1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
19 // lu@physi.uni-heidelberg.de
23 grep " AliTRDdEdxBaseUtils::" AliTRDdEdxBaseUtils.cxx | grep "=" -v | grep -v "[6]" | grep -v printf |wc
24 grep "(" AliTRDdEdxBaseUtils.h | grep ";" | grep -v grep | grep -v ClassDef | grep -v "{" | grep -v typedef | wc
28 #ifndef ALITRDDEDXBASEUTILS_H
29 #define ALITRDDEDXBASEUTILS_H
40 #include "TTreeStream.h"
53 class AliTRDdEdxBaseUtils
56 //===================================================================================
58 //===================================================================================
59 static void BinLogX(TAxis *axis);
60 static void GetCDFCuts(const TH1D *hh, const Int_t ncut, Double_t cuts[], const Double_t cdfs[], const Double_t thres);
61 static Double_t GetMeanRMS(const Double_t nn, const Double_t sum, const Double_t w2s, Double_t * grms=0x0, Double_t * gerr=0x0);
62 static Double_t TruncatedMean(const Int_t nx, const Double_t xdata[], const Double_t lowfrac, const Double_t highfrac, Double_t * grms=0x0, Double_t * gerr=0x0, Double_t *wws=0x0);
63 static Double_t TruncatedMean(const TH1 *hh, const Double_t lowfrac, const Double_t highfrac, Double_t * grms=0x0, Double_t * gerr=0x0);
64 static void FitSlicesY(const TH2D *hh, TH1D *&hnor, TH1D *&hmpv, TH1D *&hwid, TH1D *&hres, const Double_t thres, const Double_t lowfrac, const Double_t highfrac);
66 //===================================================================================
67 // TRD Analysis Fast Tool
68 //===================================================================================
69 static Int_t GetNtracklet(const AliESDEvent *esd);
70 static AliTRDtrackV1 * GetTRDtrackV1(const AliESDtrack * esdtrack);
71 static Bool_t IsInSameStack(const AliTRDtrackV1 *trdtrack);
72 static AliTRDseedV1 * GetLastTracklet(const AliTRDtrackV1 *trdtrack);
73 static AliTRDseedV1 * GetFirstSectorStackMomentum(const AliTRDtrackV1 *trdtrack, Int_t & isec, Int_t & istk, Double_t & mom);
74 static Double_t GetDeltaPhi(const AliTRDseedV1 *tracklet);
76 //===================================================================================
77 // Detector, Data and Control Constant
78 //===================================================================================
80 static Int_t NTRDchamber(){return 18*5*6;} //540
81 static Int_t NTRDtimebin(){return NTRDchamber()*31;} //16740
82 static Int_t ToDetector(const Int_t gtb);
83 static Int_t ToTimeBin(const Int_t gtb);
84 static Int_t ToSector(const Int_t gtb);
85 static Int_t ToStack(const Int_t gtb);
86 static Int_t ToLayer(const Int_t gtb);
88 static TString GetRunType(const Int_t run);
90 static void SetQ0Frac(const Double_t q0){ fgQ0Frac = q0; }
91 static void SetQ1Frac(const Double_t q1){ fgQ1Frac = q1; }
92 static void SetTimeBinCountCut(const Double_t tbc){ fgTimeBinCountCut = tbc; }
93 static void SetCalibTPCnclsCut(const Int_t tpc){ fgCalibTPCnclsCut = tpc; }
94 static void SetExBOn(const Bool_t kon){ fgExBOn = kon; }
95 static void SetPadGainOn(const Bool_t kon){ fgPadGainOn = kon; }
96 static void SetQScale(const Double_t scale){ fgQScale = scale; }
98 static Double_t Q0Frac(){return fgQ0Frac;}
99 static Double_t Q1Frac(){return fgQ1Frac;}
100 static Double_t TimeBinCountCut(){return fgTimeBinCountCut;}
101 static Int_t CalibTPCnclsCut(){return fgCalibTPCnclsCut;}
102 static Bool_t IsExBOn(){return fgExBOn;}
103 static Bool_t IsPadGainOn(){return fgPadGainOn;}
104 static Double_t QScale(){return fgQScale;}
106 static void PrintControl();
108 //===================================================================================
109 // dEdx Parameterization
110 //===================================================================================
111 static void FastFitdEdxTR(TH1 * hh);
113 static Double_t MeandEdx(const Double_t * xx, const Double_t * par);
114 static Double_t MeanTR(const Double_t * xx, const Double_t * par);
115 static Double_t MeandEdxTR(const Double_t * xx, const Double_t * par);
117 static Double_t QMeanTPC(const Double_t bg);
118 static Double_t Q0MeanTRDpp(const Double_t bg);
119 static Double_t Q1MeanTRDpp(const Double_t bg);
120 static Double_t Q0MeanTRDPbPb(const Double_t bg);
121 static Double_t Q1MeanTRDPbPb(const Double_t bg);
123 typedef Double_t (*FFunc)(const Double_t *xx, const Double_t *par);
125 static Double_t MeandEdxLogx(const Double_t * xx, const Double_t * par){return ToLogx(MeandEdx, xx, par);}
126 static Double_t MeanTRLogx(const Double_t * xx, const Double_t * par){return ToLogx(MeanTR, xx, par);}
127 static Double_t MeandEdxTRLogx(const Double_t * xx, const Double_t * par){return ToLogx(MeandEdxTR, xx, par);}
130 //dEdx Parameterization
131 static Double_t ToLogx(FFunc func, const Double_t * xx, const Double_t * par);
134 static Double_t fgQ0Frac; //q0frac
135 static Double_t fgQ1Frac; //q1frac
136 static Double_t fgTimeBinCountCut; //tbcut
137 static Int_t fgCalibTPCnclsCut; //tpccut
138 static Bool_t fgExBOn; //exbon
139 static Bool_t fgPadGainOn; //pad gain
140 static Double_t fgQScale; //Qscale