PHOS/AliPHOSAnalyze.h

   1 #ifndef ALIPHOSANALYZE_H
   2 #define ALIPHOSANALYZE_H
   3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   4  * See cxx source for full Copyright notice                               */
   5
   6 /* $Id$ */
   7
   8 //_________________________________________________________________________
   9 //  Algorithm class to analyze PHOSv1 events
  10 //*-- Author : Yves Schutz (SUBATECH)
  11
  12 // --- ROOT system ---
  13
  14 #include "TFile.h"
  15 #include "TH1.h"
  16 #include "TH2.h"
  17
  18 // --- Standard library ---
  19
  20 // --- AliRoot header files ---
  21
  22 #include "AliPHOSv1.h"
  23 #include "AliPHOSGeometry.h"
  24 #include "AliPHOSPID.h"
  25 #include "AliPHOSIndexToObject.h"
  26
  27 class AliPHOSAnalyze : public TObject {
  28
  29 public:
  30
  31   AliPHOSAnalyze() ;              // ctor
  32   AliPHOSAnalyze(Text_t * name) ; // ctor
  33   AliPHOSAnalyze(const AliPHOSAnalyze & ana) ; // cpy ctor
  34   virtual ~AliPHOSAnalyze() ;     // dtor
  35
  36   void ActivePPSD(Int_t Nevents) ;
  37   void AnalyzeManyEvents(Int_t Nevtents = 100, Int_t Module=0) ;  // analyzes many events   ;
  38   void ReconstructCPV(Int_t Nevents = 1) ;  // reconstruct points in EMC and CPV
  39   void AnalyzeResolutions(Int_t Nevtents) ; // analyzes Energy and Position resolutions   ;
  40   void ReadAndPrintCPV(Int_t Nevents);      // Read & print generated and reconstructed hits in CPV
  41   void AnalyzeCPV(Int_t Nevents);           // analyzes various CPV characteristics
  42   void InvariantMass(Int_t Nevents = 100) ;
  43   void Reconstruct(Int_t Nevtents = 100,Int_t FirstEvent = 0) ;
  44   void AnalyzeResolutions(Int_t Nevtents) ;  // analyzes Energy and Position resolutions   ;
  45   void BookingHistograms() ;                // booking histograms for the ManyEvent analysis ;
  46   void BookResolutionHistograms() ;         // booking histograms for the Resoluion analysis ;
  47   void Copy(TObject & obj) ;                // copies an analysis into an other one
  48   Float_t CorrectEnergy(Float_t ERecPart) ;   //Corrects reconstracted energy
  49   Bool_t Init(Int_t evt) ;                  // does various initialisations
  50   void DisplayKineEvent(Int_t evt = -999) ; // displays the Kine events in ALICE coordinate
  51   void DisplayRecParticles() ;              // displays RecParticles in ALICE coordinate
  52   void DisplayRecPoints() ;                 // displays RecPoints in module coordinate
  53   void DisplayTrackSegments() ;             // displays TrackSegments in module coordinate
  54   Bool_t OpenRootFile(Text_t * name) ;      // opens the root file
  55   void SaveHistograms() ;                   // Save histograms in a root file
  56   void ResetHistograms() ;                  //
  57   AliPHOSAnalyze & operator = (const AliPHOSAnalyze & rvalue)  {
  58     // assignement operator requested by coding convention
  59     // but not needed
  60     assert(0==1) ;
  61     return *this ;
  62   }
  63   void SetDebugLevel(Int_t flag) { fDebugLevel = flag; }
  64
  65  private:
  66
  67   AliPHOSClusterizer * fClu ;         // a clusterizer
  68   Int_t fEvt ;                        // the evt number being processed
  69   AliPHOSGeometry * fGeom ;           // the PHOS Geometry object
  70   AliPHOSIndexToObject * fObjGetter ; // provides methods to retrieve objects from their index in a list
  71   AliPHOSv1 * fPHOS ;                 // the PHOS object from the root file
  72   AliPHOSPID * fPID ;                 // a particle identifier
  73   AliPHOSReconstructioner * fRec ;    // a reconstructioner
  74   TFile * fRootFile ;                 // the root file that contains the data
  75   AliPHOSTrackSegmentMaker * fTrs ;   // a tracksegmentmaker ;
  76
  77   Int_t fDebugLevel;                  // debug level for analysis
  78
  79   TH2F * fhEnergyCorrelations ;     //Energy correlations between Eloss in Convertor and PPSD(2)
  80
  81
  82   TH1F * fhEmcDigit ;               // Histo of digit energies in the Emc
  83   TH1F * fhVetoDigit ;              // Histo of digit energies in the Veto
  84   TH1F * fhConvertorDigit ;         // Histo of digit energies in the Convertor
  85   TH1F * fhEmcCluster ;             // Histo of Cluster energies in Emc
  86   TH1F * fhVetoCluster ;            // Histo of Cluster energies in Veto
  87   TH1F * fhConvertorCluster ;       // Histo of Cluster energies in Convertor
  88   TH2F * fhConvertorEmc ;           // 2d Convertor versus Emc energies
  89
  90   TH2F * fhAllEnergy ;           // Spectrum of detected photons with photon primary
  91   TH2F * fhPhotEnergy ;        // Total spectrum of detected photons
  92   TH2F * fhEMEnergy ;        // Spectrum of detected neutral EM with EM primary
  93   TH2F * fhPPSDEnergy ;
  94
  95   TH2F * fhAllPosition ;        // Position Resolution of  photons with photon primary
  96   TH2F * fhPhotPosition ;     // Position Resolution of  photons
  97   TH2F * fhEMPosition ;     // Position Resolution of neutral EM with EM primary
  98   TH2F * fhPPSDPosition ;  // Position Resolution of neutral EM
  99
 100   TH1F * fhPhotonPositionY ;        // Y distribution of detected photons
 101   TH1F * fhElectronPositionY ;      // Y distribution of detected electrons
 102   TH1F * fhNeutralHadronPositionY ; // Y distribution of detected neutral hadron
 103   TH1F * fhNeutralEMPositionY ;     // Y distribution of detected neutral EM
 104   TH1F * fhChargedHadronPositionY ; // Y distribution of detected charged
 105   TH1F * fhPhotonHadronPositionY ;  // Y distribution of detected Photon-Hadron
 106
 107   TH1F * fhPhotReg ;
 108   TH1F * fhAllReg ;
 109   TH1F * fhNReg ;
 110   TH1F * fhNBarReg ;
 111   TH1F * fhChargedReg ;
 112   TH1F * fhPhotEM ;
 113   TH1F * fhAllEM ;
 114   TH1F * fhNEM ;
 115   TH1F * fhNBarEM ;
 116   TH1F * fhChargedEM ;
 117   TH1F * fhPhotPPSD ;
 118   TH1F * fhAllPPSD ;
 119   TH1F * fhNPPSD ;
 120   TH1F * fhNBarPPSD ;
 121   TH1F * fhChargedPPSD ;
 122
 123   TH1F * fhPrimary ;
 124   TH1F * fhAllRP ;
 125   TH1F * fhPPSD ;
 126   TH1F * fhShape ;
 127   TH1F * fhVeto ;
 128
 129   TH1F * fhPhotPhot ;
 130   TH1F * fhPhotElec ;
 131   TH1F * fhPhotNeuH ;
 132   TH1F * fhPhotNuEM ;
 133   TH1F * fhPhotChHa ;
 134   TH1F * fhPhotGaHa ;
 135
 136
 137
 138 ClassDef(AliPHOSAnalyze,1)  // PHOSv1 event analyzis algorithm
 139
 140 };
 141
 142 #endif // AliPHOSANALYZE_H