| d02977ee |
1 | #ifndef AliStaObjects_H |
| 2 | #define AliStaObjects_H |
| 3 | |
| 4 | // $Id$ |
| 5 | |
| 6 | #include <TLorentzVector.h> |
| 7 | |
| 8 | class AliStaHeader |
| 9 | { |
| 10 | public: |
| 11 | AliStaHeader() : fRun(0), fOrbit(0), fPeriod(0), fBx(0), fL0(0), fL1(0), fL2(0), |
| 12 | fTrClassMask(0), fTrCluster(0), fOffTriggers(0), fFiredTriggers(), |
| 13 | fTcls(0), fV0And(0), fIsHT(0), fIsPileup(0), fIsPileup2(0), fIsPileup4(0), fIsPileup8(0), |
| 14 | fNSpdVertices(0), fNTpcVertices(0), fV0Cent(0), fV0(0), fCl1Cent(0), fCl1(0), fTrCent(0), |
| 15 | fTr(0), fCqual(-1), fPsi(0), fPsiRes(0), fNSelTr(0), fNSelPrimTr(0), fNSelPrimTr1(0), |
| 16 | fNSelPrimTr2(0), fNCells(0), fNCells0(0), fNCells01(0), fNCells03(0), |
| 17 | fNCells1(0), fNCells2(0), fNCells5(0), fNClus(0), fNClus1(0), fNClus2(0), fNClus5(0), |
| 18 | fMaxCellE(0), fMaxClusE(0), fMaxTrE(0), fNcSM0(0), fNcSM1(0), fNcSM2(0), fNcSM3(0), |
| 19 | fNcSM4(0), fNcSM5(0), fNcSM6(0),fNcSM7(0),fNcSM8(0),fNcSM9(0) {;} |
| 20 | virtual ~AliStaHeader() {;} |
| 21 | |
| 22 | ULong64_t GetEventId() const { |
| 23 | return (((ULong64_t)fPeriod << 36) | |
| 24 | ((ULong64_t)fOrbit << 12) | |
| 25 | (ULong64_t)fBx); |
| 26 | } |
| 27 | |
| 28 | public: |
| 29 | Int_t fRun; // run number |
| 30 | UInt_t fOrbit; // orbit number |
| 31 | UInt_t fPeriod; // period number |
| 32 | UShort_t fBx; // bunch crossing id |
| 33 | UInt_t fL0; // l0 trigger bits |
| 34 | UInt_t fL1; // l1 trigger bits |
| 35 | UShort_t fL2; // l2 trigger bits |
| 36 | ULong64_t fTrClassMask; // trigger class mask |
| 37 | UChar_t fTrCluster; // trigger cluster mask |
| 38 | UInt_t fOffTriggers; // fired offline triggers for this event |
| 39 | TString fFiredTriggers; // string with fired triggers |
| 40 | UInt_t fTcls; // custom trigger definition |
| 41 | Bool_t fV0And; // if V0AND (from AliTriggerAnalysis) |
| 42 | Bool_t fIsHT; // if EMCAL L0 (from AliTriggerAnalysis) |
| 43 | Bool_t fIsPileup; // indicate pileup from IsPileupFromSPD with 0.8 minzdist |
| 44 | Bool_t fIsPileup2; // indicate pileup from IsPileupFromSPD with 0.4 minzdist |
| 45 | Bool_t fIsPileup4; // indicate pileup from IsPileupFromSPD with 0.2 minzdist |
| 46 | Bool_t fIsPileup8; // indicate pileup from IsPileupFromSPD with 0.1 minzdist |
| 47 | UShort_t fNSpdVertices; // number of pileup vertices (spd) |
| 48 | UShort_t fNTpcVertices; // number of pileup vertices (tpc) |
| 49 | Double32_t fV0Cent; //[0,0,16] v0 cent |
| 50 | Double32_t fV0; //[0,0,16] v0 result used for cent |
| 51 | Double32_t fCl1Cent; //[0,0,16] cl1 cent |
| 52 | Double32_t fCl1; //[0,0,16] cl1 result used for cent |
| 53 | Double32_t fTrCent; //[0,0,16] tr cent |
| 54 | Double32_t fTr; //[0,0,16] tr result used for cent |
| 55 | Int_t fCqual; // centrality quality |
| 56 | Double32_t fPsi; //[0,0,16] event-plane angle |
| 57 | Double32_t fPsiRes; //[0,0,16] event-plane ange resolution |
| 58 | UShort_t fNSelTr; // # selected tracks |
| 59 | UShort_t fNSelPrimTr; // # selected tracks (primary) |
| 60 | UShort_t fNSelPrimTr1; // # selected tracks (primary) pt > 1 GeV/c |
| 61 | UShort_t fNSelPrimTr2; // # selected tracks (primary) pt > 2 GeV/c |
| 62 | UShort_t fNCells; // # cells |
| 63 | UShort_t fNCells0; // # cells > 0.45 GeV |
| 64 | UShort_t fNCells01; // # cells > 0.1 GeV |
| 65 | UShort_t fNCells03; // # cells > 0.3 GeV |
| 66 | UShort_t fNCells1; // # cells > 1 GeV |
| 67 | UShort_t fNCells2; // # cells > 2 GeV |
| 68 | UShort_t fNCells5; // # cells > 5 GeV |
| 69 | UShort_t fNClus; // # clus |
| 70 | UShort_t fNClus1; // # clus > 1 GeV |
| 71 | UShort_t fNClus2; // # clus > 2 GeV |
| 72 | UShort_t fNClus5; // # clus > 5 GeV |
| 73 | Double32_t fMaxCellE; //[0,0,16] maximum cell energy |
| 74 | Double32_t fMaxClusE; //[0,0,16] maximum clus energy |
| 75 | Double32_t fMaxTrE; //[0,0,16] maximum trigger energy |
| 76 | UShort_t fNcSM0; // # cells > 0.1 GeV in SM 0 |
| 77 | UShort_t fNcSM1; // # cells > 0.1 GeV in SM 1 |
| 78 | UShort_t fNcSM2; // # cells > 0.1 GeV in SM 2 |
| 79 | UShort_t fNcSM3; // # cells > 0.1 GeV in SM 3 |
| 80 | UShort_t fNcSM4; // # cells > 0.1 GeV in SM 4 |
| 81 | UShort_t fNcSM5; // # cells > 0.1 GeV in SM 5 |
| 82 | UShort_t fNcSM6; // # cells > 0.1 GeV in SM 6 |
| 83 | UShort_t fNcSM7; // # cells > 0.1 GeV in SM 7 |
| 84 | UShort_t fNcSM8; // # cells > 0.1 GeV in SM 8 |
| 85 | UShort_t fNcSM9; // # cells > 0.1 GeV in SM 9 |
| 86 | |
| 87 | ClassDef(AliStaHeader,6) // Header class |
| 88 | }; |
| 89 | |
| 90 | class AliStaVertex |
| 91 | { |
| 92 | public: |
| 93 | AliStaVertex(Double_t x=0, Double_t y=0, Double_t z=0) : fVx(x), fVy(y), fVz(z), fVc(-1), fDisp(0), fZres(0), |
| 94 | fChi2(0), fSt(0), fIs3D(0), fIsZ(0) {;} |
| 95 | virtual ~AliStaVertex() {;} |
| 96 | |
| 97 | public: |
| 98 | Double_t fVx; //[0,0,16] vertex x |
| 99 | Double_t fVy; //[0,0,16] vertex y |
| 100 | Double_t fVz; //[0,0,16] vertex z |
| 101 | Double_t fVc; //[0,0,16] number of contributors to vertex |
| 102 | Double_t fDisp; //[0,0,16] dispersion |
| 103 | Double_t fZres; //[0,0,16] z-resolution |
| 104 | Double_t fChi2; //[0,0,16] chi2 of fit |
| 105 | Bool_t fSt; // status bit |
| 106 | Bool_t fIs3D; // is vertex from 3D |
| 107 | Bool_t fIsZ; // is vertex from Z only |
| 108 | |
| 109 | ClassDef(AliStaVertex,1) // Vertex class |
| 110 | }; |
| 111 | |
| 112 | class AliStaCluster : public TObject |
| 113 | { |
| 114 | public: |
| 115 | AliStaCluster() : TObject(), |
| 116 | fE(0), fR(0), fEta(0), fPhi(0), fN(0), fN1(0), fN3(0), fIdMax(-1), fSM(-1), fEmax(0), fE2max(0), fEcross(0), |
| 117 | fTmax(0), fDbc(-1), fDisp(-1), fM20(-1), fM02(-1), fEcc(-1), fSig(-1), fSigEtaEta(-1), fSigPhiPhi(-1), |
| 118 | fIsTrackM(0), fTrDz(0), fTrDr(-1), fTrEp(0), fTrDedx(0), fTrIso(0), fTrIso1(0), fTrIso2(0), |
| 119 | fTrIsoD1(0), fTrIso1D1(0), fTrIso2D1(0), fTrIsoD3(0), fTrIso1D3(0), fTrIso2D3(0), |
| 120 | fTrIsoD4(0), fTrIso1D4(0), fTrIso2D4(0), fTrIsoStrip(0), fCeIso(0), fCeIso1(0), |
| 121 | fCeIso3(0), fCeIso4(0), fCeIso3x3(0), fCeIso4x4(0), fCeIso5x5(0), fCeCore(0), fCeIso3x22(0), |
| 122 | fIsShared(0), fTrigId(-1), fTrigE(0), fMcLabel(-1), fEmbE(0) {;} |
| 123 | |
| 124 | void GetMom(TLorentzVector& p, Double_t *vertex=0); |
| 125 | void GetMom(TLorentzVector& p, AliStaVertex *vertex); |
| 126 | |
| 127 | public: |
| 128 | Double32_t fE; //[0,0,16] energy |
| 129 | Double32_t fR; //[0,0,16] radius (cylinder) |
| 130 | Double32_t fEta; //[0,0,16] eta |
| 131 | Double32_t fPhi; //[0,0,16] phi |
| 132 | UChar_t fN; // number of cells |
| 133 | UChar_t fN1; // number of cells > 100 MeV |
| 134 | UChar_t fN3; // number of cells > 300 MeV |
| 135 | Short_t fIdMax; // id maximum cell |
| 136 | Char_t fSM; // super module number (from maximum cell) |
| 137 | Double32_t fEmax; //[0,0,16] energy of maximum cell |
| 138 | Double32_t fE2max; //[0,0,16] energy of second maximum cell |
| 139 | Double32_t fEcross; //[0,0,16] energy of the 4 adjacent cells around the seed |
| 140 | Double32_t fTmax; //[0,0,16] time of maximum cell |
| 141 | Double32_t fDbc; //[0,0,16] distance to nearest bad channel |
| 142 | Double32_t fDisp; //[0,0,16] cluster dispersion, for shape analysis |
| 143 | Double32_t fM20; //[0,0,16] 2-nd moment along the main eigen axis |
| 144 | Double32_t fM02; //[0,0,16] 2-nd moment along the second eigen axis |
| 145 | Double32_t fEcc; //[0,0,16] eccentricity |
| 146 | Double32_t fSig; //[0,0,16] sigma |
| 147 | Double32_t fSigEtaEta; //[0,0,16] sigma eta-eta |
| 148 | Double32_t fSigPhiPhi; //[0,0,16] sigma phi-phi |
| 149 | Bool_t fIsTrackM; // if true then track values are set |
| 150 | Double32_t fTrDz; //[0,0,16] dZ to nearest track |
| 151 | Double32_t fTrDr; //[0,0,16] dR to nearest track (in x,y) |
| 152 | Double32_t fTrEp; //[0,0,16] E/P to nearest track |
| 153 | Double32_t fTrDedx; //[0,0,16] dE/dx (TPC signal) to nearest track |
| 154 | Double32_t fTrIso; //[0,0,16] track isolation |
| 155 | Double32_t fTrIso1; //[0,0,16] track isolation (pt>1GeV/c) |
| 156 | Double32_t fTrIso2; //[0,0,16] track isolation (pt>2GeV/c) |
| 157 | Double32_t fTrIsoD1; //[0,0,16] track isolation, iso dist 0.25 |
| 158 | Double32_t fTrIso1D1; //[0,0,16] track isolation (pt>1GeV/c), iso dist 0.1 |
| 159 | Double32_t fTrIso2D1; //[0,0,16] track isolation (pt>2GeV/c), iso dist 0.1 |
| 160 | Double32_t fTrIsoD3; //[0,0,16] track isolation, iso dist 0.3 |
| 161 | Double32_t fTrIso1D3; //[0,0,16] track isolation (pt>1GeV/c), iso dist 0.3 |
| 162 | Double32_t fTrIso2D3; //[0,0,16] track isolation (pt>2GeV/c), iso dist 0.3 |
| 163 | Double32_t fTrIsoD4; //[0,0,16] track isolation, iso dist 0.4 |
| 164 | Double32_t fTrIso1D4; //[0,0,16] track isolation (pt>1GeV/c), iso dist 0.4 |
| 165 | Double32_t fTrIso2D4; //[0,0,16] track isolation (pt>2GeV/c), iso dist 0.4 |
| 166 | Double32_t fTrIsoStrip; //[0,0,16] track isolation strip, dEtaXdPhi=0.015x0.3 |
| 167 | Double32_t fCeIso; //[0,0,16] cell isolation in R=0.20 |
| 168 | Double32_t fCeIso1; //[0,0,16] cell isolation in R=0.10 |
| 169 | Double32_t fCeIso3; //[0,0,16] cell isolation in R=0.30 |
| 170 | Double32_t fCeIso4; //[0,0,16] cell isolation in R=0.40 |
| 171 | Double32_t fCeIso3x3; //[0,0,16] cell isolation in 3x3 cells |
| 172 | Double32_t fCeIso4x4; //[0,0,16] cell isolation in 4x4 cells |
| 173 | Double32_t fCeIso5x5; //[0,0,16] cell isolation in 5x5 cells |
| 174 | Double32_t fCeCore; //[0,0,16] cell content in R=0.05 |
| 175 | Double32_t fCeIso3x22; //[0,0,16] cell isolation in rectangular strip of dEtaXdPhi=0.042x0.308 |
| 176 | Bool_t fIsShared; // =true then extends across more than one super module |
| 177 | Short_t fTrigId; // index of matched trigger tower |
| 178 | Double32_t fTrigE; //[0,0,16] energy (FEE) of matched trigger tower |
| 179 | Short_t fMcLabel; // index of closest MC particle |
| 180 | Double32_t fEmbE; //[0,0,16] sum of energy of embedded (MC) cells in cluster |
| 181 | |
| 182 | ClassDef(AliStaCluster,10) // Cluster class |
| 183 | }; |
| 184 | |
| 185 | class AliStaTrigger : public TObject |
| 186 | { |
| 187 | public: |
| 188 | AliStaTrigger() : TObject(), fE(0), fEta(0), fPhi(0), fIdMax(-1) {} |
| 189 | |
| 190 | public: |
| 191 | Double32_t fE; //[0,0,16] energy |
| 192 | Double32_t fEta; //[0,0,16] eta |
| 193 | Double32_t fPhi; //[0,0,16] phi |
| 194 | Short_t fIdMax; // id maximum cell |
| 195 | |
| 196 | ClassDef(AliStaTrigger,2) // Trigger class |
| 197 | }; |
| 198 | |
| 199 | class AliStaPart : public TObject |
| 200 | { |
| 201 | public: |
| 202 | AliStaPart() : TObject(), fPt(0), fEta(0), fPhi(0), fVR(0), fVEta(0), fVPhi(0), fPid(0), fMo(-1), fDet(-2), |
| 203 | fLab(-1), fNs(0) { memset(fDs,-1,sizeof(Short_t)*99); } |
| 204 | |
| 205 | Int_t OnEmcal() const { return (fDet==8); } |
| 206 | Int_t IsSim() const { return (fDet!=-2); } |
| 207 | |
| 208 | public: |
| 209 | Double32_t fPt; //[0,0,16] pt |
| 210 | Double32_t fEta; //[0,0,16] eta |
| 211 | Double32_t fPhi; //[0,0,16] phi |
| 212 | Double32_t fVR; //[0,0,16] prod r (cylinder) |
| 213 | Double32_t fVEta; //[0,0,16] prod eta |
| 214 | Double32_t fVPhi; //[0,0,16] prod phi |
| 215 | Short_t fPid; // pid |
| 216 | Short_t fMo; // index of mother |
| 217 | Short_t fDet; // detector in which particle left trace (8 for EMCAL, see AliTrackReference.h) |
| 218 | // the following must be filled before first usage |
| 219 | Short_t fLab; //! label (index in array) |
| 220 | Short_t fNs; //! number of daughters |
| 221 | Short_t fDs[99]; //! daughters |
| 222 | |
| 223 | ClassDef(AliStaPart,1) // Particle class |
| 224 | }; |
| 225 | #endif |
git://git.uio.no / u / mrichter / AliRoot.git / blame