The CTP (trigger) infor is written into the ESD. The corresponding entries are fTrigg...
[u/mrichter/AliRoot.git] / STEER / AliESD.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3  *                                                                        *
4  * Author: The ALICE Off-line Project.                                    *
5  * Contributors are mentioned in the code where appropriate.              *
6  *                                                                        *
7  * Permission to use, copy, modify and distribute this software and its   *
8  * documentation strictly for non-commercial purposes is hereby granted   *
9  * without fee, provided that the above copyright notice appears in all   *
10  * copies and that both the copyright notice and this permission notice   *
11  * appear in the supporting documentation. The authors make no claims     *
12  * about the suitability of this software for any purpose. It is          *
13  * provided "as is" without express or implied warranty.                  *
14  **************************************************************************/
15
16 /* $Id$ */
17
18 //-----------------------------------------------------------------
19 //           Implementation of the ESD class
20 //   This is the class to deal with during the phisical analysis of data
21 //   This class is generated directly by the reconstruction methods
22 //      Origin: Iouri Belikov, CERN, Jouri.Belikov@cern.ch
23 //-----------------------------------------------------------------
24
25 #include "AliESD.h"
26 #include "AliESDfriend.h"
27
28 ClassImp(AliESD)
29
30 //______________________________________________________________________________
31 AliESD::AliESD():
32   fEventNumber(0),
33   fRunNumber(0),
34   fTriggerMask(0),
35   fTriggerCluster(0),
36   fRecoVersion(0),
37   fMagneticField(0),
38   fZDCN1Energy(0),
39   fZDCP1Energy(0),
40   fZDCN2Energy(0),
41   fZDCP2Energy(0),
42   fZDCEMEnergy(0),
43   fZDCParticipants(0),
44   fT0zVertex(0),
45   fT0timeStart(0),
46   fPrimaryVertex(),
47   fTracks("AliESDtrack",15000),
48   fHLTConfMapTracks("AliESDHLTtrack",25000),
49   fHLTHoughTracks("AliESDHLTtrack",15000),
50   fMuonTracks("AliESDMuonTrack",30),
51   fPmdTracks("AliESDPmdTrack",3000),
52   fTrdTracks("AliESDTrdTrack",300),
53   fV0s("AliESDv0",200),  
54   fCascades("AliESDcascade",20),
55   fKinks("AliESDkink",4000),
56   fV0MIs("AliESDV0MI",4000),
57   fCaloClusters("AliESDCaloCluster",10000),
58   fEMCALClusters(0), 
59   fFirstEMCALCluster(-1),
60   fPHOSClusters(0), 
61   fFirstPHOSCluster(-1),
62   fESDFMD(0x0)
63 {
64   for (Int_t i=0; i<24; i++) {
65     fT0time[i] = 0;
66     fT0amplitude[i] = 0;
67   }
68 }
69
70 //______________________________________________________________________________
71 AliESD::~AliESD()
72 {
73   //
74   // Standard destructor
75   //
76   fTracks.Delete();
77   fHLTConfMapTracks.Delete();
78   fHLTHoughTracks.Delete();
79   fMuonTracks.Delete();
80   fPmdTracks.Delete();
81   fTrdTracks.Delete();
82   fV0s.Delete();
83   fCascades.Delete();
84   fKinks.Delete();
85   fV0MIs.Delete();
86   fCaloClusters.Delete();
87   delete fESDFMD;
88 }
89
90 void AliESD::UpdateV0PIDs()
91 {
92   //
93   //
94   //
95   Int_t nV0 = GetNumberOfV0MIs();
96   for (Int_t i=0;i<nV0;i++){
97     AliESDV0MI * v0 = GetV0MI(i);
98     AliESDtrack* tp = GetTrack(v0->GetIndex(0));
99     AliESDtrack* tm = GetTrack(v0->GetIndex(1));
100     if (!tm || !tp){
101       printf("BBBUUUUUUUGGGG\n");
102     }
103     Double_t pp[5],pm[5];
104     tp->GetESDpid(pp);
105     tm->GetESDpid(pm);
106     v0->UpdatePID(pp,pm);    
107   }
108 }
109
110 //______________________________________________________________________________
111 void AliESD::Reset()
112 {
113   fEventNumber=0;
114   fRunNumber=0;
115   fTriggerMask=0;
116   fTriggerCluster=0;
117   fRecoVersion=0;
118   fMagneticField=0;
119   fZDCN1Energy=0;
120   fZDCP1Energy=0;
121   fZDCN2Energy=0;
122   fZDCP2Energy=0;
123   fZDCEMEnergy=0;
124   fZDCParticipants=0;
125   fT0zVertex=0;
126   fT0timeStart = 0;
127   fPrimaryVertex.Reset();
128   fTracks.Clear();
129   fHLTConfMapTracks.Clear();
130   fHLTHoughTracks.Clear();
131   fMuonTracks.Clear();
132   fPmdTracks.Clear();
133   fTrdTracks.Clear();
134   fV0s.Clear();
135   fCascades.Clear();
136   fCaloClusters.Clear();
137   fEMCALClusters=0; 
138   fFirstEMCALCluster=-1; 
139   fPHOSClusters=0; 
140   fFirstPHOSCluster=-1; 
141   if (fESDFMD) fESDFMD->Clear();
142 }
143
144 //______________________________________________________________________________
145 void AliESD::Print(Option_t *) const 
146 {
147   //
148   // Print header information of the event
149   //
150   printf("ESD run information\n");
151   printf("Event # %d Run # %d Trigger %lld Magnetic field %f \n",
152          GetEventNumber(),
153          GetRunNumber(),
154          GetTriggerMask(),
155          GetMagneticField() );
156   printf("Vertex: (%.4f +- %.4f, %.4f +- %.4f, %.4f +- %.4f) cm\n",
157          fPrimaryVertex.GetXv(), fPrimaryVertex.GetXRes(),
158          fPrimaryVertex.GetYv(), fPrimaryVertex.GetYRes(),
159          fPrimaryVertex.GetZv(), fPrimaryVertex.GetZRes());
160   printf("Event from reconstruction version %d \n",fRecoVersion);
161   printf("Number of tracks: \n");
162   printf("                 charged   %d\n", GetNumberOfTracks());
163   printf("                 hlt CF    %d\n", GetNumberOfHLTConfMapTracks());
164   printf("                 hlt HT    %d\n", GetNumberOfHLTHoughTracks());
165   printf("                 phos      %d\n", GetNumberOfPHOSClusters());
166   printf("                 emcal     %d\n", GetNumberOfEMCALClusters());
167   printf("                 muon      %d\n", GetNumberOfMuonTracks());
168   printf("                 pmd       %d\n", GetNumberOfPmdTracks());
169   printf("                 trd       %d\n", GetNumberOfTrdTracks());
170   printf("                 v0        %d\n", GetNumberOfV0s());
171   printf("                 cascades  %d\n)", GetNumberOfCascades());
172   printf("                 kinks     %d\n)", GetNumberOfKinks());
173   printf("                 V0MIs     %d\n)", GetNumberOfV0MIs());
174   printf("                 CaloClusters %d\n)", GetNumberOfCaloClusters());
175   printf("                 FMD       %s\n)", (fESDFMD ? "yes" : "no"));
176 }
177
178 void AliESD::SetESDfriend(const AliESDfriend *ev) {
179   //
180   // Attaches the complementary info to the ESD
181   //
182   if (!ev) return;
183
184   Int_t ntrk=ev->GetNumberOfTracks();
185
186   for (Int_t i=0; i<ntrk; i++) {
187     const AliESDfriendTrack *f=ev->GetTrack(i);
188     GetTrack(i)->SetFriendTrack(f);
189   }
190 }
191
192 void AliESD::GetESDfriend(AliESDfriend *ev) const {
193   //
194   // Extracts the complementary info from the ESD
195   //
196   if (!ev) return;
197
198   Int_t ntrk=GetNumberOfTracks();
199
200   for (Int_t i=0; i<ntrk; i++) {
201     const AliESDtrack *t=GetTrack(i);
202     const AliESDfriendTrack *f=t->GetFriendTrack();
203     ev->AddTrack(f);
204   }
205 }