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
  15 #include "TTreeStream.h"
  16
  17
  18 class AliTPCcalibCosmic:public AliTPCcalibBase {
  19 public:
  20   AliTPCcalibCosmic();
  21   AliTPCcalibCosmic(const Text_t *name, const Text_t *title);
  22   virtual ~AliTPCcalibCosmic();
  23
  24   virtual void      Process(AliESDEvent *event);
  25   virtual Long64_t  Merge(TCollection *li);
  26   virtual void      Analyze();
  27   //
  28   void              FindPairs(AliESDEvent *event);
  29   Bool_t            IsPair(AliExternalTrackParam *tr0, AliExternalTrackParam *tr1);
  30   void              SetGainMap(AliTPCCalPad *GainMap){fGainMap = GainMap;};
  31   static void       CalculateBetheParams(TH2F *hist, Double_t * initialParam);
  32   static Double_t   CalculateMIPvalue(TH1F * hist);
  33   AliExternalTrackParam *Invert(AliExternalTrackParam *input);
  34   //
  35   TH1F   *          GetHistNTracks(){return fHistNTracks;};
  36   TH1F   *          GetHistClusters(){return fClusters;};
  37   TH2F   *          GetHistAcorde(){return fModules;};
  38   TH1F   *          GetHistPt(){return fHistPt;};
  39   TH2F   *          GetHistDeDx(){return fDeDx;};
  40   TH1F   *          GetHistMIP(){return fDeDxMIP;};
  41   //
  42   Double_t          GetMIPvalue(){return fMIPvalue;};
  43   //
  44   static void       BinLogX(TH1 * h);   // method for correct histogram binning
  45
  46
  47 private:
  48
  49   void              FillAcordeHist(AliESDtrack *upperTrack);
  50
  51   AliTPCCalPad *fGainMap;         //  gain map from Krypton calibration
  52   TH1F  *fHistNTracks;            //  histogram showing number of ESD tracks per event
  53   TH1F  *fClusters;               //  histogram showing the number of clusters per track
  54   TH2F  *fModules;                //  2d histogram of tracks which are propagated to the ACORDE scintillator array
  55   TH1F  *fHistPt;                 //  Pt histogram of reconstructed tracks
  56   TH2F  *fDeDx;                   //  dEdx spectrum showing the different particle types
  57   TH1F  *fDeDxMIP;                //  TPC signal close to the MIP region of muons 0.4 < p < 0.45 GeV
  58
  59   Double_t fMIPvalue;             //  MIP value calculated via a fit to fDeDxMIP
  60
  61   // cuts
  62   //
  63   Float_t fCutMaxD;     // maximal distance in rfi ditection
  64   Float_t fCutTheta;    // maximal distance in theta ditection
  65   Float_t fCutMinDir;   // direction vector products
  66
  67   AliTPCcalibCosmic(const AliTPCcalibCosmic&);
  68   AliTPCcalibCosmic& operator=(const AliTPCcalibCosmic&);
  69
  70   ClassDef(AliTPCcalibCosmic, 1);
  71 };
  72
  73 #endif
  74