PHOS/AliPHOSPIDv1.h

   1 #ifndef ALIPHOSPIDV1_H
   2 #define ALIPHOSPIDV1_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 v1 of the PHOS particle identifier
  11 // Identification is based on information from CPV and EMC
  12 // Oh yeah
  13 //*-- Author: Yves Schutz (SUBATECH), Gustavo Conesa.
  14
  15 // --- ROOT system ---
  16 //class TFormula ;
  17 class TVector3 ;
  18 class TMatrixD ;
  19 class TPrincipal ;
  20
  21 // --- Standard library ---
  22
  23 // --- AliRoot header files ---
  24 class AliPHOSEmcRecPoint ;
  25 class AliPHOSRecPoint ;
  26
  27 #include "AliPHOSPID.h"
  28
  29 class  AliPHOSPIDv1 : public AliPHOSPID {
  30
  31 public:
  32
  33   AliPHOSPIDv1() ;          // ctor
  34   AliPHOSPIDv1(const char* headerFile, const char * tsBranch = "Default", const char * from = 0) ;
  35
  36   virtual ~AliPHOSPIDv1() ; // dtor
  37
  38   virtual void Exec(Option_t * option) ;
  39   virtual char * GetRecParticlesBranch()const {return (char*) fRecParticlesTitle.Data() ;}
  40   virtual char * GetTrackSegmentsBranch()const{return (char*) fTrackSegmentsTitle.Data(); }
  41   virtual const Int_t GetRecParticlesInRun() const  {return fRecParticlesInRun ;}
  42
  43   virtual void Print(Option_t * option) const {}
  44   void Print() ;
  45   // Get CpvtoEmcDistanceCut and TimeGate parameters depending on the custer energy and
  46   // Purity-Efficiency point (possible options "HIGH EFFICIENCY" "MEDIUM EFFICIENCY" "LOW
  47   // EFFICIENCY" and 3 more options changing EFFICIENCY by PURITY)
  48   Double_t GetCpvtoEmcDistanceCut(const Float_t Cluster_En, const TString Eff_Pur)const  ;
  49   Double_t GetTimeGate(const Float_t Cluster_En, const TString Eff_Pur)const  ;
  50   const TString GetPrincipalFile( )const {return fFileName ;}
  51   const TString GetPrincipalFilePar( )const {return fFileNamePar ;}
  52   void  SetPrincipalFileOptions(TString OptFileName) ;
  53
  54   // Set all parameters necessary in the PID depending on the custer energy and
  55   // Purity-Efficiency point (possible options "HIGH EFFICIENCY" "MEDIUM EFFICIENCY" "LOW
  56   // EFFICIENCY" and 3 more options changing EFFICIENCY by PURITY)
  57   void SetCpvtoEmcDistanceCut(Float_t Cluster_En, TString Eff_Pur, Float_t cut)  ;
  58   void SetTimeGate(Float_t Cluster_En, TString Eff_Pur, Float_t gate)  ;
  59   void SetEllipseXCenter(Float_t Cluster_En, TString Eff_Pur, Float_t x)  ;
  60   void SetEllipseYCenter(Float_t Cluster_En, TString Eff_Pur, Float_t y)  ;
  61   void SetEllipseAParameter(Float_t Cluster_En, TString Eff_Pur, Float_t a)  ;
  62   void SetEllipseBParameter(Float_t Cluster_En, TString Eff_Pur, Float_t b)  ;
  63   void SetEllipseAngle(Float_t Cluster_En, TString Eff_Pur, Float_t angle)  ;
  64   void SetEllipseParameters(Float_t Cluster_En, TString Eff_Pur, Float_t x, Float_t y,Float_t a, Float_t b,Float_t angle) ;
  65
  66
  67   virtual void SetTrackSegmentsBranch(const char* title) { fTrackSegmentsTitle = title;}
  68   virtual void SetRecParticlesBranch (const char* title) { fRecParticlesTitle = title;}
  69   virtual const char * Version() const { return "pid-v1" ; }
  70
  71  private:
  72
  73   const TString AliPHOSPIDv1::BranchName() const ;
  74   virtual void Init() ;
  75   void     MakeRecParticles(void ) ;
  76   Float_t  GetDistance(AliPHOSEmcRecPoint * emc, AliPHOSRecPoint * cpv, Option_t * Axis)const ; // Relative Distance CPV-EMC
  77   Int_t    GetPrincipalSign(Double_t* P, Int_t ell, Int_t eff_pur)const ; //Principal cut
  78   TVector3 GetMomentumDirection(AliPHOSEmcRecPoint * emc, AliPHOSRecPoint * cpv)const ;
  79   void     PrintRecParticles(Option_t * option) ;
  80   void     SetParameters() ;
  81   virtual void WriteRecParticles(Int_t event) ;
  82
  83  private:
  84
  85   TString                fFileName ;          // Name of the file which contains the Principal file
  86   TString                fFileNamePar ;       // Name of the file which contains the parameters
  87   TString                fOptFileName ;       // choose de parameters and principal file
  88   TString                fFrom ;              // name of Recpoints and TrackSegments
  89   TString                fHeaderFileName ;    // file name with event header
  90   TString                fTrackSegmentsTitle; // branch name with track segments
  91   TString                fRecPointsTitle ;    // branch name with rec points
  92   TString                fRecParticlesTitle ; // branch name with rec particles
  93
  94   Int_t                      fNEvent ;            //! current event number
  95   AliPHOSClusterizer *       fClusterizer ;       //! clusterizer
  96   AliPHOSTrackSegmentMaker * fTSMaker ;           //! track segment maker
  97   TPrincipal *               fPrincipal ;         //! TPrincipal from fFileName
  98   Int_t                      fRecParticlesInRun ; //! Total number of recparticles in one run
  99   Double_t *                 fX ;                 //! Principal data
 100   Double_t *                 fP ;                 //! Principal eigenvalues
 101   TMatrixD *                 fParameters ;        //! Matrix of all identification Parameters
 102
 103   ClassDef( AliPHOSPIDv1,2)  // Particle identifier implementation version 1
 104
 105 };
 106
 107 #endif // AliPHOSPIDV1_H