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