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"
19 class AliAODHeader : public AliVHeader {
24 AliAODHeader(Int_t nRun, UShort_t nBunchX, UInt_t nOrbit, UInt_t nPeriod, const Char_t *title="");
25 AliAODHeader(Int_t nRun,
33 Double_t muonMagFieldScale,
40 ULong64_t fTriggerMask,
41 UChar_t fTriggerCluster,
43 const Char_t *title="",
47 virtual ~AliAODHeader();
48 AliAODHeader(const AliAODHeader& evt);
49 AliAODHeader& operator=(const AliAODHeader& evt);
51 Int_t GetRunNumber() const { return fRunNumber; }
52 UShort_t GetBunchCrossNumber() const { return fBunchCrossNumber; }
53 UInt_t GetOrbitNumber() const { return fOrbitNumber; }
54 UInt_t GetPeriodNumber() const { return fPeriodNumber; }
55 ULong64_t GetTriggerMask() const { return fTriggerMask; }
56 UChar_t GetTriggerCluster() const { return fTriggerCluster; }
57 TString GetFiredTriggerClasses()const { return fFiredTriggers;}
58 UInt_t GetEventType() const { return fEventType; }
59 Double_t GetMagneticField() const { return fMagneticField; }
60 Double_t GetMuonMagFieldScale() const { return fMuonMagFieldScale; }
62 Double_t GetCentrality() const { return fCentrality; }
63 Double_t GetZDCN1Energy() const { return fZDCN1Energy; }
64 Double_t GetZDCP1Energy() const { return fZDCP1Energy; }
65 Double_t GetZDCN2Energy() const { return fZDCN2Energy; }
66 Double_t GetZDCP2Energy() const { return fZDCP2Energy; }
67 Double_t GetZDCEMEnergy(Int_t i) const { return fZDCEMEnergy[i]; }
68 Int_t GetRefMultiplicity() const { return fRefMult; }
69 Int_t GetRefMultiplicityPos() const { return fRefMultPos; }
70 Int_t GetRefMultiplicityNeg() const { return fRefMultNeg; }
71 Int_t GetNumberOfMuons() const { return fNMuons; }
72 Int_t GetNumberOfDimuons() const { return fNDimuons; }
74 Double_t GetQTheta(UInt_t i) const;
75 UInt_t GetNQTheta() const { return (UInt_t)fNQTheta; }
77 Double_t GetDiamondX() const {return fDiamondXY[0];}
78 Double_t GetDiamondY() const {return fDiamondXY[1];}
79 Double_t GetSigma2DiamondX() const {return fDiamondCovXY[0];}
80 Double_t GetSigma2DiamondY() const {return fDiamondCovXY[2];}
81 void GetDiamondCovXY(Float_t cov[3]) const {
82 for(Int_t i=0;i<3;i++) cov[i]=fDiamondCovXY[i]; return;
85 void SetRunNumber(Int_t nRun) { fRunNumber = nRun; }
86 void SetBunchCrossNumber(UShort_t nBx) { fBunchCrossNumber = nBx; }
87 void SetOrbitNumber(UInt_t nOr) { fOrbitNumber = nOr; }
88 void SetPeriodNumber(UInt_t nPer) { fPeriodNumber = nPer; }
89 void SetTriggerMask(ULong64_t trigMsk) { fTriggerMask = trigMsk; }
90 void SetFiredTriggerClasses(TString trig) { fFiredTriggers = trig;}
91 void SetTriggerCluster(UChar_t trigClus) { fTriggerCluster = trigClus; }
92 void SetEventType(UInt_t evttype) { fEventType = evttype; }
93 void SetMagneticField(Double_t magFld) { fMagneticField = magFld; }
94 void SetMuonMagFieldScale(Double_t magFldScl){ fMuonMagFieldScale = magFldScl; }
96 void SetCentrality(Double_t cent) { fCentrality = cent; }
97 void SetZDCN1Energy(Double_t n1Energy) { fZDCN1Energy = n1Energy; }
98 void SetZDCP1Energy(Double_t p1Energy) { fZDCP1Energy = p1Energy; }
99 void SetZDCN2Energy(Double_t n2Energy) { fZDCN2Energy = n2Energy; }
100 void SetZDCP2Energy(Double_t p2Energy) { fZDCP2Energy = p2Energy; }
101 void SetZDCEMEnergy(Double_t emEnergy1, Double_t emEnergy2)
102 { fZDCEMEnergy[0] = emEnergy1; fZDCEMEnergy[1] = emEnergy2;}
103 void SetRefMultiplicity(Int_t refMult) { fRefMult = refMult; }
104 void SetRefMultiplicityPos(Int_t refMultPos) { fRefMultPos = refMultPos; }
105 void SetRefMultiplicityNeg(Int_t refMultNeg) { fRefMultNeg = refMultNeg; }
106 void SetNumberOfMuons(Int_t nMuons) { fNMuons = nMuons; }
107 void SetNumberOfDimuons(Int_t nDimuons) { fNDimuons = nDimuons; }
109 void SetQTheta(Double_t *QTheta, UInt_t size = 5);
112 void SetDiamond(Float_t xy[2],Float_t cov[3]) {
113 for(Int_t i=0;i<3;i++) {if(i<2) fDiamondXY[i]=xy[i]; fDiamondCovXY[i]=cov[i];}
116 void Print(Option_t* option = "") const;
118 void SetPHOSMatrix(TGeoHMatrix*matrix, Int_t i) {
119 if ((i >= 0) && (i < kNPHOSMatrix)) fPHOSMatrix[i] = matrix;
121 const TGeoHMatrix* GetPHOSMatrix(Int_t i) const {
122 return ((i >= 0) && (i < kNPHOSMatrix)) ? fPHOSMatrix[i] : NULL;
125 void SetEMCALMatrix(TGeoHMatrix*matrix, Int_t i) {
126 if ((i >= 0) && (i < kNEMCALMatrix)) fEMCALMatrix[i] = matrix;
128 const TGeoHMatrix* GetEMCALMatrix(Int_t i) const {
129 return ((i >= 0) && (i < kNEMCALMatrix)) ? fEMCALMatrix[i] : NULL;
132 enum {kNPHOSMatrix = 5};
133 enum {kNEMCALMatrix = 12};
137 Double32_t fMagneticField; // Solenoid Magnetic Field in kG
138 Double32_t fMuonMagFieldScale; // magnetic field scale of muon arm magnet
139 Double32_t fCentrality; // Centrality
140 Double32_t fZDCN1Energy; // reconstructed energy in the neutron1 ZDC
141 Double32_t fZDCP1Energy; // reconstructed energy in the proton1 ZDC
142 Double32_t fZDCN2Energy; // reconstructed energy in the neutron2 ZDC
143 Double32_t fZDCP2Energy; // reconstructed energy in the proton2 ZDC
144 Double32_t fZDCEMEnergy[2]; // reconstructed energy in the electromagnetic ZDCs
145 Int_t fNQTheta; // number of QTheta elements
146 Double32_t *fQTheta; // [fNQTheta] values to store Lee-Yang-Zeros
147 ULong64_t fTriggerMask; // Trigger Type (mask)
148 TString fFiredTriggers; // String with fired triggers
149 Int_t fRunNumber; // Run Number
150 Int_t fRefMult; // reference multiplicity
151 Int_t fRefMultPos; // reference multiplicity of positive particles
152 Int_t fRefMultNeg; // reference multiplicity of negative particles
153 Int_t fNMuons; // number of muons in the forward spectrometer
154 Int_t fNDimuons; // number of dimuons in the forward spectrometer
155 UInt_t fEventType; // Type of Event
156 UInt_t fOrbitNumber; // Orbit Number
157 UInt_t fPeriodNumber; // Period Number
158 UShort_t fBunchCrossNumber; // BunchCrossingNumber
159 UChar_t fTriggerCluster; // Trigger cluster (mask)
160 Double32_t fDiamondXY[2]; // Interaction diamond (x,y) in RUN
161 Double32_t fDiamondCovXY[3]; // Interaction diamond covariance (x,y) in RUN
162 TGeoHMatrix* fPHOSMatrix[kNPHOSMatrix]; //PHOS module position and orientation matrices
163 TGeoHMatrix* fEMCALMatrix[kNEMCALMatrix]; //EMCAL supermodule position and orientation matrices
165 ClassDef(AliAODHeader,10);