StRoot: StMiniMcEvent/StMiniMcEvent.h Source File

00001 
00016 #ifndef StMiniMcEvent_H
00017 #define StMiniMcEvent_H
00018 
00019 #include "TObject.h"
00020 #include "TClonesArray.h"
00021 #include "StMiniMcPair.h"
00022 #include "StContamPair.h"
00023 
00024 enum Category { MC,MATCHED,MERGED,SPLIT,CONTAM,GHOST,MATGLOB};
00025 
00026 class StMiniMcEvent : public TObject {
00027  public:
00028   StMiniMcEvent();
00029   virtual ~StMiniMcEvent();
00030   void Clear(Option_t *option=""); // clear the tracks
00031 
00032   StTinyMcTrack* addMcTrack(StTinyMcTrack* mc=0);
00033   StMiniMcPair*  addTrackPair(StMiniMcPair* pair=0, Category K=MATCHED);
00034   StContamPair*  addContamPair(StContamPair* pair=0);
00035 
00036   //  TClonesArray* tracks(Category);
00037   TClonesArray* tracks(Category) const;
00038 
00039   Int_t         eventId()                       const { return mEventId; }
00040   Int_t         runId()                         const { return mRunId; }
00041   Int_t         originMult()                    const { return mOriginMult; }
00042   Int_t         centralMult()                   const { return mCentralMult; }
00043   Int_t         centrality()                    const { return mCentrality; }
00044   Int_t         nUncorrectedNegativePrimaries() const { return mNUncorrectedNegativePrimaries; }
00045   Int_t         nUncorrectedPrimaries()         const { return mNUncorrectedPrimaries; }
00046   Int_t         nUncorrectedGlobals()           const { return mNUncorrectedGlobals; }
00047   Int_t         nFtpcWUncorrectedPrimaries()    const { return mNFtpcWUncorrectedPrimaries; }
00048   Int_t         nFtpcEUncorrectedPrimaries()    const { return mNFtpcEUncorrectedPrimaries; }
00049   Int_t         mcMult()                        const { return mMcMult; }
00050   Int_t         nMcNch()                        const { return mNMcNch; }
00051   Int_t         nMcFtpcWNch()                   const { return mNMcFtpcWNch; }
00052   Int_t         nMcFtpcENch()                   const { return mNMcFtpcENch; }
00053   Int_t         nMcHminus()                     const { return mNMcHminus; }
00054   Int_t         nMcGlobal()                     const { return mNMcGlobal; }
00055   Int_t         nMcGoodGlobal20()               const { return mNMcGoodGlobal20; }
00056   Int_t         nRcGlobal()                     const { return mNRcGlobal; }
00057   Int_t         nRcGoodGlobal20()               const { return mNRcGoodGlobal20; }
00058   Float_t       vertexX()                       const { return mVertexX; }  
00059   Float_t       vertexY()                       const { return mVertexY; }
00060   Float_t       vertexZ()                       const { return mVertexZ; }
00061   const Float_t*vertexCovMatrix()               const { return &mVertexCovMatrix[0];}
00062   Float_t       mcVertexX()                     const { return mMcVertexX       ; }
00063   Float_t       mcVertexY()                     const { return mMcVertexY       ; }
00064   Float_t       mcVertexZ()                     const { return mMcVertexZ       ; }
00065   Float_t       centerOfMassEnergy()            const { return mCenterOfMassEnergy; }
00066   Float_t       magneticField()                 const { return mMagField        ; }
00067   Float_t       backgroundRate()                const { return mBackgroundRate  ; }
00068   Short_t       beamMassNumberEast()            const { return mBeamMassNumberEast; }
00069   Short_t       beamMassNumberWest()            const { return mBeamMassNumberWest; }
00070   Float_t       ctb()                           const { return mCtb             ; }
00071   Float_t       zdcE()                          const { return mZdcE            ; }
00072   Float_t       zdcW()                          const { return mZdcW            ; }
00073   Int_t         nMcTrack()                      const { return mNMcTrack        ; }     
00074   Int_t         nMatchedPair()                  const { return mNMatchedPair    ; }
00075   Int_t         nMergedPair()                   const { return mNMergedPair     ; }
00076   Int_t         nSplitPair()                    const { return mNSplitPair      ; }
00077   Int_t         nGhostPair()                    const { return mNGhostPair      ; }
00078   Int_t         nContamPair()                   const { return mNContamPair     ; }
00079 
00080   float impact()                                const { return mImpact          ; }
00081   float impactPhi()                             const { return mImpactPhi       ; }
00082   float timeOffset()                            const { return mTimeOffset      ; }
00083 
00084   void setEventId(Int_t val)                            { mEventId=val; }
00085   void setRunId(Int_t val)                              { mRunId=val; }
00086   void setOriginMult(Int_t val)                         { mOriginMult=val; }
00087   void setCentralMult(Int_t val)                        { mCentralMult=val; }
00088   void setCentrality(Int_t val)                         { mCentrality=val; }
00089   void setNUncorrectedNegativePrimaries(Int_t val)      { mNUncorrectedNegativePrimaries=val; }
00090   void setNUncorrectedPrimaries(Int_t val)              { mNUncorrectedPrimaries=val; }
00091   void setNUncorrectedGlobals(Int_t val)   { mNUncorrectedGlobals=val; }
00092   void setNFtpcWUncorrectedPrimaries(Int_t val)         { mNFtpcWUncorrectedPrimaries=val; }
00093   void setNFtpcEUncorrectedPrimaries(Int_t val)         { mNFtpcEUncorrectedPrimaries=val; }
00094   void setMcMult(Int_t val)                             { mMcMult=val; }
00095   void setNMcNch(Int_t val)                             { mNMcNch=val; }
00096   void setNMcFtpcWNch(Int_t val)                        { mNMcFtpcWNch=val; }
00097   void setNMcFtpcENch(Int_t val)                        { mNMcFtpcENch=val; }
00098   void setNMcHminus(Int_t val)                          { mNMcHminus=val; }
00099   void setNMcGlobal(Int_t val)                          { mNMcGlobal=val; }
00100   void setNMcGoodGlobal20(Int_t val)                    { mNMcGoodGlobal20=val; }
00101   void setNRcGlobal(Int_t val)                          { mNRcGlobal=val; }
00102   void setNRcGoodGlobal20(Int_t val)                    { mNRcGoodGlobal20=val; }
00103   void setVertexX(Float_t val)                          { mVertexX=val; }
00104   void setVertexY(Float_t val)                          { mVertexY=val; }
00105   void setVertexZ(Float_t val)                          { mVertexZ=val; }
00106   void setVertexCovMatrix(Float_t *cov) {for (Int_t i = 0; i < 6; i++) mVertexCovMatrix[i] = cov[i];}
00107   void setMcVertexX(Float_t val)                        { mMcVertexX=val; }
00108   void setMcVertexY(Float_t val)                        { mMcVertexY=val; }
00109   void setMcVertexZ(Float_t val)                        { mMcVertexZ=val; } 
00110   void setCenterOfMassEnergy(Float_t val)               { mCenterOfMassEnergy=val; }
00111   void setMagField(Float_t val)                         { mMagField=val; }
00112   void setBackgroundRate(Float_t val)                   { mBackgroundRate=val; }
00113   void setBeamMassNumberEast(Short_t val)               { mBeamMassNumberEast=val; }
00114   void setBeamMassNumberWest(Short_t val)               { mBeamMassNumberWest=val; }
00115   void setCtb(Float_t val)                              { mCtb=val; }
00116   void setZdcE(Float_t val)                             { mZdcE=val; }
00117   void setZdcW(Float_t val)                             { mZdcW=val; }
00118   void setNMcTrack(Int_t val)                           { mNMcTrack     =val; }
00119   void setNMatchedPair(Int_t val)                       { mNMatchedPair =val; }
00120   void setNMergedPair(Int_t val)                        { mNMergedPair  =val; }
00121   void setNSplitPair(Int_t val)                         { mNSplitPair   =val; }
00122   void setNGhostPair(Int_t val)                         { mNGhostPair   =val; }
00123   void setNContamPair(Int_t val)                        { mNContamPair  =val; }
00124 
00125   void setImpact(float imp)                             { mImpact       =imp  ;}
00126   void setImpactPhi(float imphi)                        { mImpactPhi    =imphi;}
00127   void setTimeOffset(float time)                        { mTimeOffset   =time ;}
00128 
00129 
00130   virtual void Print(Option_t *option="") const;
00131 private:
00132   //
00133   // data members
00134   //
00135   Int_t         mEventId;
00136   Int_t         mRunId;       // set to 0 for simulations
00137   Int_t         mOriginMult;  // StEvent::primaryVertex(0)->numberOfDaughters()
00138   Int_t         mCentralMult; // reco, primary trk, flag>0, |eta|<0.75 
00139   Int_t         mCentrality;  // centrality bin, Nch cuts, P02gd, 2k2
00140   Int_t         mNUncorrectedNegativePrimaries; // from StuRefMult
00141   Int_t         mNUncorrectedPrimaries; // from StuRefMult
00142   Int_t         mNUncorrectedGlobals; // reco, globals, glTrk->helix->dist(vtx)<3, gl.fitPts>=10, -0.5 < eta < 0.5  
00143   Int_t         mNFtpcWUncorrectedPrimaries; // reco, primaries, flag>0, glTrk->helix->dist(vtx)<3, prim.pt<3, gl.fitPts>=5, 2.8 < eta < 3.8
00144   Int_t         mNFtpcEUncorrectedPrimaries; // reco, primaries, flag>0, glTrk->helix->dist(vtx)<3, prim.pt<3, gl.fitPts>=5,-2.8 > eta >-3.8
00145   Int_t         mMcMult;      // embedding: n mc tracks; (StMcEvent::numberOfPrimaryTracks()
00146                               // simulation: same as mOriginMult 
00147   Int_t         mNMcNch;      // mc, primary, charge!=0, |eta|<0.5
00148   Int_t         mNMcFtpcWNch; // mc, primary, charge!=0, 2.8 < eta < 3.8
00149   Int_t         mNMcFtpcENch; // mc, primary, charge!=0,-2.8 > eta >-3.8
00150   Int_t         mNMcHminus;   // mc, primary, charg  <0, |eta|<0.5
00151 
00152   Int_t         mNMcGlobal;             // mc, primary, |eta|<4
00153   Int_t         mNMcGoodGlobal20;       // mc, primary, |eta|<4, MC TPC hits >=20
00154   Int_t         mNRcGlobal;             // reco, primaries flag > 0
00155   Int_t         mNRcGoodGlobal20;       // reco, primaries flag > 0, 20 fit hits
00156 
00157   Float_t       mImpact;        // Impact parameter
00158   Float_t       mImpactPhi;     // Phi Impact parameter, azimuth of reaction plane 
00159   Float_t       mTimeOffset;    // time offset in seconds wrt trigger event
00160 
00161   Float_t       mVertexX;  
00162   Float_t       mVertexY;     // 
00163   Float_t       mVertexZ;
00164   Float_t       mVertexCovMatrix[6];
00165   Float_t       mMcVertexX;
00166   Float_t       mMcVertexY;
00167   Float_t       mMcVertexZ;
00168     
00169   Float_t       mMagField;    // in kGauss
00170 
00171   Float_t   mCenterOfMassEnergy;
00172   Float_t   mBackgroundRate;
00173   Short_t   mBeamMassNumberEast;
00174   Short_t   mBeamMassNumberWest;
00175 
00176   Float_t       mCtb;
00177   Float_t       mZdcE;
00178   Float_t       mZdcW;
00179 
00180   //
00181   // for root reasons, even though ghosts are just the rc tracks,
00182   // save it as a StMiniMcPair ... 
00183 
00184   Int_t mNMcTrack;     
00185   Int_t mNMatchedPair;
00186   Int_t mNMergedPair;
00187   Int_t mNSplitPair;
00188   Int_t mNGhostPair;
00189   Int_t mNContamPair;
00190   Int_t mNMatGlobPair;
00191     
00192   TClonesArray* mMcTracks;
00193   TClonesArray* mMatchedPairs;
00194   TClonesArray* mMergedPairs;
00195   TClonesArray* mSplitPairs;
00196   TClonesArray* mGhostPairs; 
00197   TClonesArray* mContamPairs;
00198   TClonesArray* mMatGlobPairs;
00199   static Int_t mSFirst; 
00200 
00201   ClassDef(StMiniMcEvent,5)
00202 };
00203 
00204 #endif
00205 
00206   
00207 //
00208 // $Log: StMiniMcEvent.h,v $
00209 // Revision 1.9  2012/03/15 23:37:20  perev
00210 // Uncorrected globals added(Chris)
00211 //
00212 // Revision 1.8  2011/03/22 00:31:48  perev
00213 // Added impact,phi impact & trigger time
00214 //
00215 // Revision 1.7  2007/05/21 16:17:16  fisyak
00216 // Increament ClassDef, thanks to Adam Kocoloski
00217 //
00218 // Revision 1.6  2007/02/23 17:07:00  fisyak
00219 // Add Ssd and DCA
00220 //
00221 // Revision 1.5  2004/01/26 13:58:18  calderon
00222 // Introduction of global matched branch.
00223 //
00224 // Revision 1.4  2003/07/09 01:05:51  calderon
00225 // Added the FTPC reference multiplicity, East and West
00226 // Added various other multiplicity definitions, good globals reco and MC,
00227 // and MC Nch and Nh-.
00228 // Keep data members private and provide setters and getters
00229 //
00230 // Revision 1.3  2003/05/08 02:09:20  calderon
00231 // Added data members for svt and ftpc fit points for StTinyRcTrack.
00232 // Added data members for svt and ftpc hits for StTinyMcTrack.
00233 // Added methods to calculate px, py, and p from the available pt,  phi and pz, for
00234 // global and primary momenta and also for monte carlo momentum.
00235 // Cleaned up includes in StMiniMcEvent.
00236 //
00237 // Revision 1.2  2002/06/06 15:14:13  calderon
00238 // the comment about the magnetic field said it is in tesla, but it is actually
00239 // in kGauss, so I modified the comment to reflect this
00240 //
00241 // Revision 1.1  2002/05/30 01:20:57  calderon
00242 // Classes for use in a general framework for extracting efficiencies
00243 // from both embedding and full simulations
00244 // (after GSTAR+TRS+StEvent+StMcEvent+StAssociationMaker)
00245 // so that the information of the track matches gets stored persistently.
00246 //
00247 //