]> git.uio.no Git - u/mrichter/AliRoot.git/blob - TPC/AliTPCcalibCosmic.h
Updated flags for low flux case (A. Dainese)
[u/mrichter/AliRoot.git] / TPC / AliTPCcalibCosmic.h
1 #ifndef ALITPCCALIBCOSMIC_H
2 #define ALITPCCALIBCOSMIC_H
3
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5  * See cxx source for full Copyright notice                               */
6
7 #include "AliTPCcalibBase.h"
8 #include "AliTPCCalPad.h"
9 #include "TH2F.h"
10 class TH1F;
11 class TList;
12 class AliESDEvent;
13 class AliESDtrack;
14 class AliExternalComparison;
15
16 #include "TTreeStream.h"
17
18
19 class AliTPCcalibCosmic:public AliTPCcalibBase {
20 public:
21   AliTPCcalibCosmic(); 
22   AliTPCcalibCosmic(const Text_t *name, const Text_t *title);
23   virtual ~AliTPCcalibCosmic();
24   
25   virtual void      Process(AliESDEvent *event);
26   virtual Long64_t  Merge(TCollection *li);
27   virtual void      Analyze();
28   //
29   void ProcessTree(TTree * tree, AliExternalComparison *comp=0);
30   //
31   void              FindPairs(AliESDEvent *event);
32   Bool_t            IsPair(AliExternalTrackParam *tr0, AliExternalTrackParam *tr1);
33   void              SetGainMap(AliTPCCalPad *GainMap){fGainMap = GainMap;};
34   static void       CalculateBetheParams(TH2F *hist, Double_t * initialParam);
35   static Double_t   CalculateMIPvalue(TH1F * hist);
36   AliExternalTrackParam *Invert(AliExternalTrackParam *input);
37   AliExternalTrackParam *MakeTrack(const AliExternalTrackParam *track0, const AliExternalTrackParam *track1);
38   void UpdateTrack(AliExternalTrackParam &track0, const AliExternalTrackParam &track1);
39   //
40   TH1F   *          GetHistNTracks(){return fHistNTracks;};
41   TH1F   *          GetHistClusters(){return fClusters;};
42   TH2F   *          GetHistAcorde(){return fModules;};
43   TH1F   *          GetHistPt(){return fHistPt;};
44   TH2F   *          GetHistDeDx(){return fDeDx;};
45   TH1F   *          GetHistMIP(){return fDeDxMIP;};
46   //
47   Double_t          GetMIPvalue(){return fMIPvalue;};
48   //
49   static void       BinLogX(TH1 * h);   // method for correct histogram binning
50
51   void     Process(AliESDtrack *track, Int_t runNo=-1){AliTPCcalibBase::Process(track,runNo);};
52   void     Process(AliTPCseed *track){return AliTPCcalibBase::Process(track);}
53
54 private:
55
56   void              FillAcordeHist(AliESDtrack *upperTrack);
57
58   AliTPCCalPad *fGainMap;         //  gain map from Krypton calibration
59   TH1F  *fHistNTracks;            //  histogram showing number of ESD tracks per event
60   TH1F  *fClusters;               //  histogram showing the number of clusters per track
61   TH2F  *fModules;                //  2d histogram of tracks which are propagated to the ACORDE scintillator array
62   TH1F  *fHistPt;                 //  Pt histogram of reconstructed tracks
63   TH2F  *fDeDx;                   //  dEdx spectrum showing the different particle types
64   TH1F  *fDeDxMIP;                //  TPC signal close to the MIP region of muons 0.4 < p < 0.45 GeV
65
66   Double_t fMIPvalue;             //  MIP value calculated via a fit to fDeDxMIP
67   //
68   // cuts
69   //
70   Float_t fCutMaxD;     // maximal distance in rfi ditection
71   Float_t fCutMaxDz;     // maximal distance in z ditection
72   Float_t fCutTheta;    // maximal distance in theta ditection
73   Float_t fCutMinDir;   // direction vector products
74
75   AliTPCcalibCosmic(const AliTPCcalibCosmic&); 
76   AliTPCcalibCosmic& operator=(const AliTPCcalibCosmic&); 
77
78   ClassDef(AliTPCcalibCosmic, 1); 
79 };
80
81 #endif
82