3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
9 //_________________________________________________________________________
10 // Implementation version v1 of the PHOS particle identifier
11 // Identification is based on information from CPV and EMC
13 //*-- Author: Yves Schutz (SUBATECH), Gustavo Conesa.
15 // --- ROOT system ---
20 // --- Standard library ---
22 // --- AliRoot header files ---
23 class AliPHOSEmcRecPoint ;
24 class AliPHOSRecPoint ;
26 #include "AliPHOSPID.h"
28 class AliPHOSPIDv1 : public AliPHOSPID {
32 AliPHOSPIDv1() ; // ctor
33 AliPHOSPIDv1(const TString alirunFileNameFile, const TString eventFolderName = AliConfig::fgkDefaultEventFolderName) ;
34 AliPHOSPIDv1(const AliPHOSPIDv1 & pid) ; // cpy ctor
36 virtual ~AliPHOSPIDv1() ; // dtor
38 virtual void Exec(Option_t *option); // Does the job
40 //Get file name that contain the PCA
41 const TString GetFileNamePrincipal(TString particle) const;
43 //Get file name that contain PID parameters
44 const TString GetFileNameParameters() const {return fFileNameParameters ;}
46 // Get number of rec.particles in this run
47 virtual Int_t GetRecParticlesInRun() const {return fRecParticlesInRun ;}
49 // Get PID parameters as they are defined in fParameters
50 Float_t GetParameterCalibration (Int_t i) const;
51 Float_t GetParameterCpv2Emc (Int_t i, TString axis) const;
52 Float_t GetParameterTimeGate (Int_t i) const;
53 Float_t GetParameterToCalculateEllipse(TString particle, TString param, Int_t i) const ;
54 Float_t GetParameterPhotonBoundary (Int_t i) const;
55 Float_t GetParameterPi0Boundary (Int_t i) const;
57 // Get energy-dependent PID parameters
58 Float_t GetCalibratedEnergy (Float_t e) const;
59 Float_t GetCpv2EmcDistanceCut (TString axis, Float_t e) const ;
60 Float_t GetEllipseParameter (TString particle, TString param, Float_t e) const;
62 // Set PID parameters to change appropriate element of fParameters
63 void SetParameterCalibration (Int_t i, Float_t param);
64 void SetParameterCpv2Emc (Int_t i, TString axis, Float_t cut) ;
65 void SetParameterTimeGate (Int_t i, Float_t gate) ;
66 void SetParameterToCalculateEllipse(TString particle, TString param, Int_t i, Float_t value) ;
67 void SetParameterPhotonBoundary(Int_t i, Float_t param);
68 void SetParameterPi0Boundary (Int_t i, Float_t param);
72 virtual const char * Version() const { return "pid-v1" ; }
74 AliPHOSPIDv1 & operator = (const AliPHOSPIDv1 & /*pid*/) { return *this ;}
78 const TString BranchName() const ;
80 virtual void InitParameters() ;
81 void MakeRecParticles(void ) ;
82 // Relative Distance CPV-EMC
83 Float_t GetDistance (AliPHOSEmcRecPoint * emc, AliPHOSRecPoint * cpv, Option_t * axis)const ;
84 Int_t GetCPVBit (AliPHOSEmcRecPoint * emc, AliPHOSRecPoint * cpv, Int_t EffPur, Float_t e) const;
85 Int_t GetPrincipalBit (TString particle, const Double_t* P, Int_t EffPur, Float_t e)const ; //Principal cut
86 Int_t GetHardPhotonBit(AliPHOSEmcRecPoint * emc) const;
87 Int_t GetHardPi0Bit (AliPHOSEmcRecPoint * emc) const;
88 TVector3 GetMomentumDirection(AliPHOSEmcRecPoint * emc, AliPHOSRecPoint * cpv)const ;
89 void PrintRecParticles(Option_t * option) ;
90 virtual void WriteRecParticles() ;
91 void SetParameters() ; //Fills the matrix of parameters
96 Bool_t fDefaultInit; //! kTRUE if the task was created by defaut ctor (only parameters are initialized)
97 Int_t fNEvent ; //! current event number
98 TString fFileNamePrincipalPhoton ; // File name of the photon principals
99 TString fFileNamePrincipalPi0 ; // File name of the pi0 principals
100 TString fFileNameParameters ; // File name with PID parameters
101 TPrincipal *fPrincipalPhoton ; //! TPrincipal from photon pca file
102 TPrincipal *fPrincipalPi0 ; //! TPrincipal from pi0 pca file
103 Double_t *fX ; //! Shower shape for the principal data
104 Double_t *fPPhoton ; //! Principal photon eigenvalues
105 Double_t *fPPi0 ; //! Principal pi0 eigenvalues
106 Int_t fRecParticlesInRun ; //! Total number of recparticles in one run
107 TMatrix *fParameters; //! Matrix of identification Parameters
110 ClassDef( AliPHOSPIDv1,9) // Particle identifier implementation version 1
114 #endif // AliPHOSPIDV1_H