PHOS/AliPHOSPIDv0.h

   1 #ifndef ALIPHOSPIDV0_H
   2 #define ALIPHOSPIDV0_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 //_________________________________________________________________________
  10 // Implementation version v0 of the PHOS particle identifier
  11 // Identification is based on information from PPSD and EMC
  12 // Oh yeah
  13 //*-- Author: Yves Schutz (SUBATECH)
  14
  15 // --- ROOT system ---
  16 class TFormula ;
  17 class TVector3 ;
  18
  19 // --- Standard library ---
  20
  21 // --- AliRoot header files ---
  22 class AliPHOSEmcRecPoint ;
  23 class AliPHOSRecPoint ;
  24
  25 #include "AliPHOSPID.h"
  26
  27 class  AliPHOSPIDv0 : public AliPHOSPID {
  28
  29 public:
  30
  31   AliPHOSPIDv0() ;          // ctor
  32   AliPHOSPIDv0(const char* headerFile, const char * tsBranch = "Default", const Bool_t toSplit=kFALSE) ;
  33   AliPHOSPIDv0(AliPHOSPIDv0 & pid) ;          // cpy ctor
  34   virtual ~AliPHOSPIDv0() ; // dtor
  35
  36   virtual void Exec(Option_t * option);
  37   //  virtual char * GetRecParticlesBranch()const {return (char*) fRecParticlesTitle.Data() ;}
  38   //  virtual char * GetTrackSegmentsBranch()const{return (char*) fTrackSegmentsTitle.Data(); }
  39   virtual const Int_t GetRecParticlesInRun() const  {return fRecParticlesInRun ;}
  40
  41   virtual void PlotDispersionCuts()const ;
  42   virtual void Print(Option_t * option)const ;
  43   virtual void SetIdentificationMethod(char * option = "CPV DISP" ){fIDOptions = option ;}
  44   virtual void SetShowerProfileCut(char * formula = "0.35*0.35 - (x-1.386)*(x-1.386) - 1.707*1.707*(y-1.008)*(y-1.008)") ;
  45   virtual void SetDispersionCut(Float_t cut){fDispersion = cut ; }
  46   virtual void SetCpvtoEmcDistanceCut(Float_t cut )      {fCpvEmcDistance = cut ;}
  47   virtual void SetTimeGate(Float_t gate)                 {fTimeGate = gate ;}
  48   //  virtual void SetTrackSegmentsBranch(const char* title) { fTrackSegmentsTitle = title;}
  49   //  virtual void SetRecParticlesBranch (const char* title) { fRecParticlesTitle = title;}
  50   virtual const char * Version() const { return "pid-v0" ; }
  51   AliPHOSPIDv0 & operator = (const  AliPHOSPIDv0 & pid) { return *this ; }
  52  private:
  53
  54   virtual void Init() ;
  55   void     MakeRecParticles(void ) ;
  56   Float_t  GetDistance(AliPHOSEmcRecPoint * emc, AliPHOSRecPoint * cpv, Option_t * Axis)const ; // Relative Distance CPV-EMC
  57   TVector3 GetMomentumDirection(AliPHOSEmcRecPoint * emc, AliPHOSRecPoint * cpv)const ;
  58   void     PrintRecParticles(Option_t * option) ;
  59   virtual void WriteRecParticles(Int_t event) ;
  60
  61  private:
  62
  63   TString                fHeaderFileName ;    // file name with event header
  64   TString                fTrackSegmentsTitle; // branch name with track segments
  65   TString                fRecPointsTitle ;    // branch name with rec points
  66   TString                fRecParticlesTitle ; // branch name with rec particles
  67   TString                fIDOptions ;         // PID option
  68   Int_t                  fNEvent ;            // current event number
  69
  70   AliPHOSClusterizer   * fClusterizer ;       // !
  71   AliPHOSTrackSegmentMaker * fTSMaker ;       // !
  72
  73   TFormula             * fFormula ;           // formula to define cut on the shower elips axis
  74   Float_t                fDispersion ;        // dispersion cut
  75   Float_t                fCpvEmcDistance ;    // Max EMC-CPV distance
  76   Float_t                fTimeGate ;          // Time of the latest EmcRecPoint accepted as EM
  77   Int_t                  fRecParticlesInRun ; //! Total number of recparticles in one run
  78
  79   ClassDef( AliPHOSPIDv0,1)  // Particle identifier implementation version 1
  80
  81 };
  82
  83 #endif // AliPHOSPIDV0_H