]> git.uio.no Git - u/mrichter/AliRoot.git/blob - STEER/AliESD.h
Reconstruction of multiplicity. Ipdated AlITSvertexerZ (Massimo)
[u/mrichter/AliRoot.git] / STEER / AliESD.h
1 // -*- mode: C++ -*- 
2 #ifndef ALIESD_H
3 #define ALIESD_H
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5  * See cxx source for full Copyright notice                               */
6
7
8 /* $Id$ */
9
10 //-------------------------------------------------------------------------
11 //                          Class AliESD
12 //   This is the class to deal with during the physical analysis of data
13 //      
14 //         Origin: Iouri Belikov, CERN, Jouri.Belikov@cern.ch 
15 //-------------------------------------------------------------------------
16
17 #include <TClonesArray.h>
18 #include <TObject.h>
19
20 #include "AliESDMuonTrack.h"
21 #include "AliESDPmdTrack.h"
22 #include "AliESDTrdTrack.h"
23 #include "AliESDVertex.h"
24 #include "AliESDcascade.h"
25 #include "AliESDkink.h"
26 #include "AliESDtrack.h"
27 #include "AliESDHLTtrack.h"
28 #include "AliESDCaloCluster.h"
29 #include "AliESDv0.h"
30 #include "AliESDV0MI.h"
31 #include "AliESDFMD.h"
32 #include "AliMultiplicity.h"
33
34 class AliESDfriend;
35
36 class AliESD : public TObject {
37 public:
38   AliESD();
39   virtual ~AliESD(); 
40
41   void SetESDfriend(const AliESDfriend *f);
42   void GetESDfriend(AliESDfriend *f) const;
43
44   void SetEventNumber(Int_t n) {fEventNumber=n;}
45   void SetRunNumber(Int_t n) {fRunNumber=n;}
46   void SetTriggerMask(ULong64_t n) {fTriggerMask=n;}
47   void SetTriggerCluster(UChar_t n) {fTriggerCluster = n;}
48   void SetMagneticField(Float_t mf){fMagneticField = mf;}
49   Float_t GetMagneticField() const {return fMagneticField;}
50   
51   AliESDtrack *GetTrack(Int_t i) const {
52     return (AliESDtrack *)fTracks.UncheckedAt(i);
53   }
54   AliESDHLTtrack *GetHLTConfMapTrack(Int_t i) const {
55     return (AliESDHLTtrack *)fHLTConfMapTracks.UncheckedAt(i);
56   }
57   AliESDHLTtrack *GetHLTHoughTrack(Int_t i) const {
58     return (AliESDHLTtrack *)fHLTHoughTracks.UncheckedAt(i);
59   }
60   AliESDMuonTrack *GetMuonTrack(Int_t i) const {
61     return (AliESDMuonTrack *)fMuonTracks.UncheckedAt(i);
62   }
63   AliESDPmdTrack *GetPmdTrack(Int_t i) const {
64     return (AliESDPmdTrack *)fPmdTracks.UncheckedAt(i);
65   }
66   AliESDTrdTrack *GetTrdTrack(Int_t i) const {
67     return (AliESDTrdTrack *)fTrdTracks.UncheckedAt(i);
68   }
69
70   Int_t  AddTrack(const AliESDtrack *t) {
71     AliESDtrack * track = new(fTracks[fTracks.GetEntriesFast()]) AliESDtrack(*t);
72     track->SetID(fTracks.GetEntriesFast()-1);
73     return  track->GetID();
74     
75   }
76   void AddHLTConfMapTrack(const AliESDHLTtrack *t) {
77     new(fHLTConfMapTracks[fHLTConfMapTracks.GetEntriesFast()]) AliESDHLTtrack(*t);
78   }
79   void AddHLTHoughTrack(const AliESDHLTtrack *t) {
80     new(fHLTHoughTracks[fHLTHoughTracks.GetEntriesFast()]) AliESDHLTtrack(*t);
81   }
82   void AddMuonTrack(const AliESDMuonTrack *t) {
83     new(fMuonTracks[fMuonTracks.GetEntriesFast()]) AliESDMuonTrack(*t);
84   }
85   void AddPmdTrack(const AliESDPmdTrack *t) {
86     new(fPmdTracks[fPmdTracks.GetEntriesFast()]) AliESDPmdTrack(*t);
87   }
88   void AddTrdTrack(const AliESDTrdTrack *t) {
89     new(fTrdTracks[fTrdTracks.GetEntriesFast()]) AliESDTrdTrack(*t);
90   }
91
92   AliESDv0 *GetV0(Int_t i) const {
93     return (AliESDv0 *)fV0s.UncheckedAt(i);
94   }
95   void AddV0(const AliESDv0 *v) {
96     new(fV0s[fV0s.GetEntriesFast()]) AliESDv0(*v);
97   }
98   void UpdateV0PIDs();
99
100   AliESDcascade *GetCascade(Int_t i) const {
101     return (AliESDcascade *)fCascades.UncheckedAt(i);
102   }
103   void AddCascade(const AliESDcascade *c) {
104     new(fCascades[fCascades.GetEntriesFast()]) AliESDcascade(*c);
105   }
106
107   AliESDkink *GetKink(Int_t i) const {
108     return (AliESDkink *)fKinks.UncheckedAt(i);
109   }
110   Int_t AddKink(const AliESDkink *c) {
111     AliESDkink * kink = new(fKinks[fKinks.GetEntriesFast()]) AliESDkink(*c);
112     kink->SetID(fKinks.GetEntriesFast());
113     return fKinks.GetEntriesFast()-1;
114   }
115
116   AliESDV0MI *GetV0MI(Int_t i) const {
117     return (AliESDV0MI *)fV0MIs.UncheckedAt(i);
118   }
119   Int_t AddV0MI(const AliESDV0MI *c) {
120     AliESDV0MI * v0 = new(fV0MIs[fV0MIs.GetEntriesFast()]) AliESDV0MI(*c);
121     v0->SetID(fV0MIs.GetEntriesFast()-1);
122     return fV0MIs.GetEntriesFast()-1;
123   }
124
125   AliESDCaloCluster *GetCaloCluster(Int_t i) const {
126     return (AliESDCaloCluster *)fCaloClusters.UncheckedAt(i);
127   }
128   Int_t AddCaloCluster(const AliESDCaloCluster *c) {
129     AliESDCaloCluster *clus = new(fCaloClusters[fCaloClusters.GetEntriesFast()]) AliESDCaloCluster(*c);
130     clus->SetID(fCaloClusters.GetEntriesFast()-1);
131     return fCaloClusters.GetEntriesFast()-1;
132   }
133     
134   void SetVertex(const AliESDVertex *vertex) {
135      new (&fSPDVertex) AliESDVertex(*vertex);
136   }
137   const AliESDVertex *GetVertex() const {return &fSPDVertex;}
138
139   void SetMultiplicity(const AliMultiplicity *mul) {
140      new (&fSPDMult) AliMultiplicity(*mul);
141   }
142   const AliMultiplicity *GetMultiplicity() const {return &fSPDMult;}
143
144   void SetPrimaryVertex(const AliESDVertex *vertex) {
145      new (&fPrimaryVertex) AliESDVertex(*vertex);
146   }
147   const AliESDVertex *GetPrimaryVertex() const {return &fPrimaryVertex;}
148
149   Int_t  GetEventNumber() const {return fEventNumber;}
150   Int_t  GetRunNumber() const {return fRunNumber;}
151   ULong64_t GetTriggerMask() const {return fTriggerMask;}
152   UChar_t  GetTriggerCluster() const {return fTriggerCluster;}
153   
154   Int_t GetNumberOfTracks()     const {return fTracks.GetEntriesFast();}
155   Int_t GetNumberOfHLTConfMapTracks()     const {return fHLTConfMapTracks.GetEntriesFast();}
156   Int_t GetNumberOfHLTHoughTracks()     const {return fHLTHoughTracks.GetEntriesFast();}
157   Int_t GetNumberOfMuonTracks() const {return fMuonTracks.GetEntriesFast();}
158   Int_t GetNumberOfPmdTracks() const {return fPmdTracks.GetEntriesFast();}
159   Int_t GetNumberOfTrdTracks() const {return fTrdTracks.GetEntriesFast();}
160   Int_t GetNumberOfV0s()      const {return fV0s.GetEntriesFast();}
161   Int_t GetNumberOfCascades() const {return fCascades.GetEntriesFast();}
162   Int_t GetNumberOfKinks() const {return fKinks.GetEntriesFast();}
163   Int_t GetNumberOfV0MIs() const {return fV0MIs.GetEntriesFast();}
164   Int_t GetNumberOfCaloClusters() const {return fCaloClusters.GetEntriesFast();}
165
166   Int_t GetNumberOfEMCALClusters() const {return fEMCALClusters;}
167   void  SetNumberOfEMCALClusters(Int_t clus) {fEMCALClusters = clus;}
168   Int_t GetFirstEMCALCluster() const {return fFirstEMCALCluster;}
169   void  SetFirstEMCALCluster(Int_t index) {fFirstEMCALCluster = index;}
170
171   Int_t GetNumberOfPHOSClusters() const {return fPHOSClusters;}
172   void  SetNumberOfPHOSClusters(Int_t part) { fPHOSClusters = part ; }
173   void  SetFirstPHOSCluster(Int_t index) { fFirstPHOSCluster = index ; } 
174   Int_t GetFirstPHOSCluster() const  { return fFirstPHOSCluster ; }
175
176   Float_t GetT0zVertex() const {return fT0zVertex;}
177   void SetT0zVertex(Float_t z) {fT0zVertex=z;}
178   Float_t GetT0() const {return fT0timeStart;}
179   void SetT0(Float_t timeStart) {fT0timeStart = timeStart;}
180   const Float_t * GetT0time() const {return fT0time;}
181   void SetT0time(Float_t time[24]) {
182     for (Int_t i=0; i<24; i++) fT0time[i] = time[i];
183   }
184   const Float_t * GetT0amplitude() const {return fT0amplitude;}
185   void SetT0amplitude(Float_t amp[24]) {
186     for (Int_t i=0; i<24; i++) fT0amplitude[i] = amp[i];
187   }
188
189   Float_t GetZDCN1Energy() const {return fZDCN1Energy;}
190   Float_t GetZDCP1Energy() const {return fZDCP1Energy;}
191   Float_t GetZDCN2Energy() const {return fZDCN2Energy;}
192   Float_t GetZDCP2Energy() const {return fZDCP2Energy;}
193   Float_t GetZDCEMEnergy() const {return fZDCEMEnergy;}
194   Int_t   GetZDCParticipants() const {return fZDCParticipants;}
195   void    SetZDC(Float_t n1Energy, Float_t p1Energy, Float_t emEnergy,
196                  Float_t n2Energy, Float_t p2Energy, Int_t participants) 
197    {fZDCN1Energy=n1Energy; fZDCP1Energy=p1Energy; fZDCEMEnergy=emEnergy;
198     fZDCN2Energy=n2Energy; fZDCP2Energy=p2Energy; fZDCParticipants=participants;}
199
200   void ResetV0s() { fV0s.Clear(); }
201   void ResetCascades() { fCascades.Clear(); }
202   void Reset();
203
204   void  Print(Option_t *option="") const;
205
206   void SetFMDData(AliESDFMD * obj) {
207     fESDFMD = new AliESDFMD(*obj);
208   }
209
210   AliESDFMD * GetFMDData(){ return fESDFMD;}
211    
212 protected:
213
214   // Event Identification
215   Int_t        fEventNumber;     // Event Number
216   Int_t        fRunNumber;       // Run Number
217   ULong64_t    fTriggerMask;     // Trigger Type (mask)
218   UChar_t      fTriggerCluster;  // Trigger cluster (mask)
219   Int_t        fRecoVersion;     // Version of reconstruction 
220   Float_t      fMagneticField;   // Solenoid Magnetic Field in kG : for compatibility with AliMagF
221
222   Float_t      fZDCN1Energy;      // reconstructed energy in the neutron ZDC
223   Float_t      fZDCP1Energy;      // reconstructed energy in the proton ZDC
224   Float_t      fZDCN2Energy;      // reconstructed energy in the neutron ZDC
225   Float_t      fZDCP2Energy;      // reconstructed energy in the proton ZDC
226   Float_t      fZDCEMEnergy;     // reconstructed energy in the electromagnetic ZDC
227   Int_t        fZDCParticipants; // number of participants estimated by the ZDC
228
229   Float_t      fT0zVertex;       // vertex z position estimated by the START
230   AliESDVertex fSPDVertex;       // Primary vertex estimated by the SPD
231   AliESDVertex fPrimaryVertex;   // Primary vertex estimated using ESD tracks
232   AliMultiplicity fSPDMult;      // SPD tracklet multiplicity
233
234   Float_t      fT0timeStart;     // interaction time estimated by the START
235   Float_t      fT0time[24];      // best TOF on each START PMT
236   Float_t      fT0amplitude[24]; // number of particles(MIPs) on each START PMT
237
238   TClonesArray fTracks;          // ESD tracks
239   TClonesArray fHLTConfMapTracks;// HLT ESD tracks from Conformal Mapper method
240   TClonesArray fHLTHoughTracks;  // HLT ESD tracks from Hough Transform method
241   TClonesArray fMuonTracks;      // MUON ESD tracks
242   TClonesArray fPmdTracks;       // PMD ESD tracks
243   TClonesArray fTrdTracks;       // TRD ESD tracks (triggered)
244   TClonesArray fV0s;             // V0 vertices
245   TClonesArray fCascades;        // Cascade vertices
246   TClonesArray fKinks;           // Kinks
247   TClonesArray fV0MIs;           // V0MI
248   TClonesArray fCaloClusters;    // Calorimeter clusters for PHOS/EMCAL
249   Int_t        fEMCALClusters;   // Number of EMCAL clusters (subset of caloclusters)
250   Int_t        fFirstEMCALCluster; // First EMCAL cluster in the fCaloClusters list 
251
252   Int_t        fPHOSClusters;     // Number of PHOS clusters (subset of caloclusters)
253   Int_t        fFirstPHOSCluster; // First PHOS cluster in the fCaloClusters list 
254  
255   AliESDFMD *  fESDFMD; // FMD object containing rough multiplicity
256
257   ClassDef(AliESD,13)  //ESD class 
258 };
259 #endif 
260