26d4b141 |
1 | #ifndef ALIPHOSPIDV1_H |
2 | #define ALIPHOSPIDV1_H |
6ad0bfa0 |
3 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * |
4 | * See cxx source for full Copyright notice */ |
5 | |
6 | /* $Id$ */ |
7 | |
b2a60966 |
8 | |
9 | //_________________________________________________________________________ |
10 | // Implementation version v1 of the PHOS particle identifier |
11 | // Identification is based on information from PPSD and EMC |
2f04ed65 |
12 | // Oh yeah |
b2a60966 |
13 | //*-- Author: Yves Schutz (SUBATECH) |
6ad0bfa0 |
14 | |
15 | // --- ROOT system --- |
7acf6008 |
16 | class TFormula ; |
17 | class TVector3 ; |
6ad0bfa0 |
18 | |
19 | // --- Standard library --- |
20 | |
21 | // --- AliRoot header files --- |
7acf6008 |
22 | class AliPHOSEmcRecPoint ; |
23 | class AliPHOSRecPoint ; |
6ad0bfa0 |
24 | |
26d4b141 |
25 | #include "AliPHOSPID.h" |
6ad0bfa0 |
26 | |
26d4b141 |
27 | class AliPHOSPIDv1 : public AliPHOSPID { |
6ad0bfa0 |
28 | |
29 | public: |
30 | |
7acf6008 |
31 | AliPHOSPIDv1() ; // ctor |
32 | AliPHOSPIDv1(const char* headerFile, const char * tsBranch = 0) ; |
33 | virtual ~AliPHOSPIDv1() ; // dtor |
34 | |
35 | virtual void Exec(Option_t * option); |
36 | virtual char * GetRecParticlesBranch()const {return (char*) fRecparticlesTitle.Data() ;} |
37 | virtual char * GetTrackSegmentsBranch()const{return (char*) fTSTitle.Data(); } |
38 | |
39 | virtual void Init() ; |
40 | |
41 | virtual void PlotDispersionCuts()const ; |
42 | |
43 | virtual void Print(Option_t * option)const ; |
44 | |
45 | virtual Bool_t ReadTrackSegments() ; |
46 | |
47 | virtual void SetIdentificationMethod(char * option = "CPV DISP" ){fIDOptions = option ;} |
48 | |
49 | virtual void SetShowerProfileCut(char * formula = |
50 | "0.35*0.35 - (x-1.386)*(x-1.386) - 1.707*1.707*(y-1.008)*(y-1.008)") ; |
51 | |
52 | virtual void SetDispersionCut(Float_t cut){fDispersion = cut ; } |
53 | virtual void SetCpvtoEmcDistanceCut(Float_t cut ) {fCpvEmcDistance = cut ;} |
54 | virtual void SetTrackSegmentsBranch(const char* title) { fTSTitle = title;} |
55 | virtual void SetRecParticlesBranch (const char* title) { fRecparticlesTitle = title;} |
56 | |
57 | virtual void WriteRecParticles() ; |
69183710 |
58 | |
6ad0bfa0 |
59 | |
7acf6008 |
60 | private: |
61 | void MakeRecParticles(void ) ; |
62 | Float_t GetDistance(AliPHOSEmcRecPoint * emc, AliPHOSRecPoint * cpv, Option_t * Axis)const ; |
63 | // Relative Distance PPSD-EMC |
64 | TVector3 GetMomentumDirection(AliPHOSEmcRecPoint * emc, AliPHOSRecPoint * cpv, AliPHOSRecPoint * ppsd)const ; |
69183710 |
65 | |
7acf6008 |
66 | void PrintRecParticles(Option_t * option) ; |
09fc14a0 |
67 | |
68 | private: |
69 | |
baef0810 |
70 | TString fHeaderFileName ; // file name with event header |
71 | TString fTSTitle; // branch name with track segments |
72 | TString fRecPointsTitle ; // branch name with rec points |
73 | TString fRecparticlesTitle ; // branch name with rec particles |
7acf6008 |
74 | |
baef0810 |
75 | TString fIDOptions ; // PID option |
7acf6008 |
76 | |
baef0810 |
77 | Int_t fNEvent ; // current event number |
78 | TObjArray * fEmcRecPoints ; // ! initial EMC RecPoints |
79 | TObjArray * fCpvRecPoints ; // ! initial CPV RecPoints |
80 | TClonesArray * fTrackSegments; // ! initial list of TrackSegments |
81 | TClonesArray * fRecParticles ; // ! output |
7acf6008 |
82 | |
baef0810 |
83 | AliPHOSClusterizer * fClusterizer ; // ! |
84 | AliPHOSTrackSegmentMaker * fTSMaker ; // ! |
7acf6008 |
85 | |
baef0810 |
86 | AliPHOSGeometry * fGeom ; // !pointer to PHOS geometry |
87 | TFormula * fFormula ; // formula to define cut on the shouer elips axis |
88 | Float_t fDispersion ; // dispersion cut |
89 | Float_t fCpvEmcDistance ; // Max EMC-CPV distance |
7acf6008 |
90 | |
baef0810 |
91 | Bool_t fIsInitialized ; // kTRUE is inifialized |
7acf6008 |
92 | |
6ad0bfa0 |
93 | |
b2a60966 |
94 | ClassDef( AliPHOSPIDv1,1) // Particle identifier implementation version 1 |
6ad0bfa0 |
95 | |
96 | }; |
97 | |
26d4b141 |
98 | #endif // AliPHOSPIDV1_H |