3 /* Copyright(c) 1998-2007, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
8 //-------------------------------------------------------------------------
9 // AOD event header class
10 // Author: Markus Oldenburg, CERN
11 //-------------------------------------------------------------------------
13 #include "AliVHeader.h"
14 #include "AliAODVertex.h"
16 #include "AliCentrality.h"
22 class AliAODHeader : public AliVHeader {
27 AliAODHeader(Int_t nRun, UShort_t nBunchX, UInt_t nOrbit, UInt_t nPeriod, const Char_t *title="");
28 AliAODHeader(Int_t nRun,
36 Double_t muonMagFieldScale,
43 ULong64_t fTriggerMask,
44 UChar_t fTriggerCluster,
46 const Char_t *title="",
50 virtual ~AliAODHeader();
51 AliAODHeader(const AliAODHeader& evt);
52 AliAODHeader& operator=(const AliAODHeader& evt);
54 Int_t GetRunNumber() const { return fRunNumber;}
55 Int_t GetEventNumberESDFile() const { return fEventNumberESDFile;}
56 UShort_t GetBunchCrossNumber() const { return fBunchCrossNumber; }
57 UInt_t GetOrbitNumber() const { return fOrbitNumber; }
58 UInt_t GetPeriodNumber() const { return fPeriodNumber; }
59 ULong64_t GetTriggerMask() const { return fTriggerMask; }
60 UChar_t GetTriggerCluster() const { return fTriggerCluster; }
61 TString GetFiredTriggerClasses()const { return fFiredTriggers;}
62 UInt_t GetEventType() const { return fEventType; }
63 Double_t GetMagneticField() const { return fMagneticField; }
64 Double_t GetMuonMagFieldScale() const { return fMuonMagFieldScale; }
66 Double_t GetCentrality() const { return fCentrality; }
67 Double_t GetZDCN1Energy() const { return fZDCN1Energy; }
68 Double_t GetZDCP1Energy() const { return fZDCP1Energy; }
69 Double_t GetZDCN2Energy() const { return fZDCN2Energy; }
70 Double_t GetZDCP2Energy() const { return fZDCP2Energy; }
71 Double_t GetZDCEMEnergy(Int_t i) const { return fZDCEMEnergy[i]; }
72 Int_t GetRefMultiplicity() const { return fRefMult; }
73 Int_t GetRefMultiplicityPos() const { return fRefMultPos; }
74 Int_t GetRefMultiplicityNeg() const { return fRefMultNeg; }
75 Int_t GetNumberOfMuons() const { return fNMuons; }
76 Int_t GetNumberOfDimuons() const { return fNDimuons; }
78 Double_t GetQTheta(UInt_t i) const;
79 UInt_t GetNQTheta() const { return (UInt_t)fNQTheta; }
81 Double_t GetDiamondX() const {return fDiamondXY[0];}
82 Double_t GetDiamondY() const {return fDiamondXY[1];}
83 Double_t GetDiamondZ() const {return fDiamondZ;}
84 Double_t GetSigma2DiamondX() const {return fDiamondCovXY[0];}
85 Double_t GetSigma2DiamondY() const {return fDiamondCovXY[2];}
86 Double_t GetSigma2DiamondZ() const {return fDiamondSig2Z;}
87 void GetDiamondCovXY(Float_t cov[3]) const {
88 for(Int_t i=0;i<3;i++) cov[i]=fDiamondCovXY[i]; return;
90 UInt_t GetL0TriggerInputs() const {return fL0TriggerInputs;}
91 UInt_t GetL1TriggerInputs() const {return fL1TriggerInputs;}
92 UShort_t GetL2TriggerInputs() const {return fL2TriggerInputs;}
93 AliCentrality* GetCentralityP() const { return fCentralityP; }
96 void SetRunNumber(Int_t nRun) { fRunNumber = nRun; }
97 void SetEventNumberESDFile(Int_t n) { fEventNumberESDFile=n; }
98 void SetBunchCrossNumber(UShort_t nBx) { fBunchCrossNumber = nBx; }
99 void SetOrbitNumber(UInt_t nOr) { fOrbitNumber = nOr; }
100 void SetPeriodNumber(UInt_t nPer) { fPeriodNumber = nPer; }
101 void SetTriggerMask(ULong64_t trigMsk) { fTriggerMask = trigMsk; }
102 void SetFiredTriggerClasses(TString trig) { fFiredTriggers = trig;}
103 void SetTriggerCluster(UChar_t trigClus) { fTriggerCluster = trigClus; }
104 void SetEventType(UInt_t evttype) { fEventType = evttype; }
105 void SetMagneticField(Double_t magFld) { fMagneticField = magFld; }
106 void SetMuonMagFieldScale(Double_t magFldScl){ fMuonMagFieldScale = magFldScl; }
108 void SetCentrality(AliCentrality* cent) {
110 if(fCentralityP)*fCentralityP = *cent;
111 else fCentralityP = new AliCentrality(*cent);
112 fCentrality = cent->GetCentralityPercentile("V0M");
118 void SetZDCN1Energy(Double_t n1Energy) { fZDCN1Energy = n1Energy; }
119 void SetZDCP1Energy(Double_t p1Energy) { fZDCP1Energy = p1Energy; }
120 void SetZDCN2Energy(Double_t n2Energy) { fZDCN2Energy = n2Energy; }
121 void SetZDCP2Energy(Double_t p2Energy) { fZDCP2Energy = p2Energy; }
122 void SetZDCEMEnergy(Double_t emEnergy1, Double_t emEnergy2)
123 { fZDCEMEnergy[0] = emEnergy1; fZDCEMEnergy[1] = emEnergy2;}
124 void SetRefMultiplicity(Int_t refMult) { fRefMult = refMult; }
125 void SetRefMultiplicityPos(Int_t refMultPos) { fRefMultPos = refMultPos; }
126 void SetRefMultiplicityNeg(Int_t refMultNeg) { fRefMultNeg = refMultNeg; }
127 void SetNumberOfMuons(Int_t nMuons) { fNMuons = nMuons; }
128 void SetNumberOfDimuons(Int_t nDimuons) { fNDimuons = nDimuons; }
130 void SetQTheta(Double_t *QTheta, UInt_t size = 5);
133 void SetDiamond(Float_t xy[2],Float_t cov[3]) {
134 for(Int_t i=0;i<3;i++) {fDiamondCovXY[i] = cov[i];}
135 for(Int_t i=0;i<2;i++) {fDiamondXY[i] = xy[i] ;}
137 void SetDiamondZ(Float_t z, Float_t sig2z){
138 fDiamondZ=z; fDiamondSig2Z=sig2z;
140 void SetL0TriggerInputs(UInt_t n) {fL0TriggerInputs=n;}
141 void SetL1TriggerInputs(UInt_t n) {fL1TriggerInputs=n;}
142 void SetL2TriggerInputs(UShort_t n) {fL2TriggerInputs=n;}
143 void SetESDFileName(TString name) {fESDFileName = name;}
144 void Print(Option_t* option = "") const;
146 void SetPHOSMatrix(TGeoHMatrix*matrix, Int_t i) {
147 if ((i >= 0) && (i < kNPHOSMatrix)) fPHOSMatrix[i] = matrix;
149 const TGeoHMatrix* GetPHOSMatrix(Int_t i) const {
150 return ((i >= 0) && (i < kNPHOSMatrix)) ? fPHOSMatrix[i] : NULL;
153 void SetEMCALMatrix(TGeoHMatrix*matrix, Int_t i) {
154 if ((i >= 0) && (i < kNEMCALMatrix)) fEMCALMatrix[i] = matrix;
156 const TGeoHMatrix* GetEMCALMatrix(Int_t i) const {
157 return ((i >= 0) && (i < kNEMCALMatrix)) ? fEMCALMatrix[i] : NULL;
160 UInt_t GetOfflineTrigger() { return fOfflineTrigger; }
161 void SetOfflineTrigger(UInt_t trigger) { fOfflineTrigger = trigger; }
162 TString GetESDFileName() {return fESDFileName;}
163 void Clear(Option_t* = "");
164 enum {kNPHOSMatrix = 5};
165 enum {kNEMCALMatrix = 12};
169 Double32_t fMagneticField; // Solenoid Magnetic Field in kG
170 Double32_t fMuonMagFieldScale; // magnetic field scale of muon arm magnet
171 Double32_t fCentrality; // Centrality
172 Double32_t fZDCN1Energy; // reconstructed energy in the neutron1 ZDC
173 Double32_t fZDCP1Energy; // reconstructed energy in the proton1 ZDC
174 Double32_t fZDCN2Energy; // reconstructed energy in the neutron2 ZDC
175 Double32_t fZDCP2Energy; // reconstructed energy in the proton2 ZDC
176 Double32_t fZDCEMEnergy[2]; // reconstructed energy in the electromagnetic ZDCs
177 Int_t fNQTheta; // number of QTheta elements
178 Double32_t *fQTheta; // [fNQTheta] values to store Lee-Yang-Zeros
179 ULong64_t fTriggerMask; // Trigger Type (mask)
180 TString fFiredTriggers; // String with fired triggers
181 Int_t fRunNumber; // Run Number
182 Int_t fRefMult; // reference multiplicity
183 Int_t fRefMultPos; // reference multiplicity of positive particles
184 Int_t fRefMultNeg; // reference multiplicity of negative particles
185 Int_t fNMuons; // number of muons in the forward spectrometer
186 Int_t fNDimuons; // number of dimuons in the forward spectrometer
187 UInt_t fEventType; // Type of Event
188 UInt_t fOrbitNumber; // Orbit Number
189 UInt_t fPeriodNumber; // Period Number
190 UShort_t fBunchCrossNumber; // BunchCrossingNumber
191 UChar_t fTriggerCluster; // Trigger cluster (mask)
192 Double32_t fDiamondXY[2]; // Interaction diamond (x,y) in RUN
193 Double32_t fDiamondCovXY[3]; // Interaction diamond covariance (x,y) in RUN
194 Double32_t fDiamondZ; // Interaction diamond (z) in RUN
195 Double32_t fDiamondSig2Z; // Interaction diamond sigma^2 (z) in RUN
196 TGeoHMatrix* fPHOSMatrix[kNPHOSMatrix]; //PHOS module position and orientation matrices
197 TGeoHMatrix* fEMCALMatrix[kNEMCALMatrix]; //EMCAL supermodule position and orientation matrices
198 UInt_t fOfflineTrigger; // fired offline triggers for this event
199 TString fESDFileName; // ESD file name to which this event belongs
200 Int_t fEventNumberESDFile; // Event number in ESD file
201 UInt_t fL0TriggerInputs; // L0 Trigger Inputs (mask)
202 UInt_t fL1TriggerInputs; // L1 Trigger Inputs (mask)
203 UShort_t fL2TriggerInputs; // L2 Trigger Inputs (mask)
204 AliCentrality* fCentralityP; // Pointer to full centrality information
205 ClassDef(AliAODHeader, 14);