StRoot: StMuDSTMaker/COMMON/StMuEvent.h Source File

00001 /***************************************************************************
00002  *
00003  * $Id: StMuEvent.h,v 1.30 2010/05/28 19:47:51 tone421 Exp $
00004  * Author: Frank Laue, BNL, laue@bnl.gov
00005  *
00006  ***************************************************************************/
00007 #ifndef StMuEvent_h
00008 #define StMuEvent_h
00009 
00010 #include "TObject.h"
00011 #include "TArrayI.h"
00012 #include "StMuL3EventSummary.h"
00013 #include "StMuEmcCollection.h"
00014 #include "StMuFmsCollection.h"
00015 #include "StMuTriggerIdCollection.h"
00016 
00017 #include "StEvent/StEventInfo.h"
00018 #include "StEvent/StRunInfo.h"
00019 #include "StEvent/StEventSummary.h"
00020 #include "StEvent/StVpdTriggerDetector.h"
00021 #include "StEvent/StMtdTriggerDetector.h"
00022 #include "StEvent/StCtbTriggerDetector.h"
00023 #include "StEvent/StZdcTriggerDetector.h"
00024 #include "StEvent/StBbcTriggerDetector.h"
00025 #include "StEvent/StEmcTriggerDetector.h"
00026 #include "StEvent/StFpdTriggerDetector.h"
00027 #include "StEvent/StFmsTriggerDetector.h"
00028 #include "StEvent/StFpdCollection.h"
00029 #include "StEvent/StL0Trigger.h"
00030 #include "StEvent/StTriggerIdCollection.h"
00031 #include "StEvent/StDetectorState.h"
00032 
00033 #include "StStrangeMuDstMaker/StStrangeMuDst.hh"
00034 //#include "StarClassLibrary/StThreeVectorD.hh"
00035 
00036 class StEvent;
00037 class StMuCut;
00038 class StTofCollection; // calibrated vpd for TOF
00039 class StTriggerData;
00040 
00047 class StMuEvent : public TObject {
00048  public:
00049   StMuEvent();
00050   StMuEvent(const StEvent*); 
00051   virtual ~StMuEvent();
00052 
00053   int eventId();
00054   int eventNumber();
00055   int runId();
00056   int runNumber();
00057   // classes taken straight from StEvent
00058   StRunInfo& runInfo();
00059   StEventInfo& eventInfo();
00060   StEventSummary& eventSummary();
00061   StVpdTriggerDetector& vpdTriggerDetector();
00062   StMtdTriggerDetector& mtdTriggerDetector();
00063   StCtbTriggerDetector& ctbTriggerDetector();
00064   StZdcTriggerDetector& zdcTriggerDetector();
00065   StBbcTriggerDetector& bbcTriggerDetector();
00066   StEmcTriggerDetector& emcTriggerDetector();
00067   StFpdTriggerDetector& fpdTriggerDetector();
00068   StFmsTriggerDetector& fmsTriggerDetector();
00069   StFpdCollection& fpdCollection(); 
00070   StL0Trigger& l0Trigger(); 
00071   // Special classes for the muDst
00072   StMuL3EventSummary& l3EventSummary();
00073   StMuTriggerIdCollection& triggerIdCollection();
00074   const StTriggerData* triggerData() const;
00075 
00077   unsigned short refMultPos(int vtx_id = -1);
00079   unsigned short refMultNeg(int vtx_id = -1);
00081   unsigned short refMult(int vtx_id = -1);
00083   unsigned short refMultFtpcEast(int vtx_id = -1);
00085   unsigned short refMultFtpcWest(int vtx_id = -1);
00087   unsigned short refMultFtpc(int vtx_id = -1);
00088   unsigned short grefmult(int vtx_id=-1);
00089   unsigned short btofTrayMultiplicity();
00090   float nearestVertexZ(int vtx_id=-1);
00091 
00093   double reactionPlane(unsigned short);
00094   void   setReactionPlane(unsigned short, double v);
00096   double reactionPlanePtWgt(unsigned short);
00097   void   setReactionPlanePtWgt(unsigned short, double v);
00098   double magneticField();
00099   double zdcAdcAttentuatedSumWest();
00100   double zdcAdcAttentuatedSumEast();
00101   double ctbMultiplicity();
00103   StThreeVectorF primaryVertexPosition(int vtx_id = -1);
00104   StThreeVectorF primaryVertexErrors(int vtx_id = -1);
00105   TArrayI& L2Result(); // Raw L2Result[] array
00106 
00107   // Calibrated VPD info from StTofCollection in StEvent
00108   unsigned int numberOfVpdEastHits();  
00109   unsigned int numberOfVpdWestHits();
00110   float vpdTstart();
00111   float vpdTdiff(); 
00112   float vpdVz();
00113 
00114 
00115  protected:
00116   void clear();
00117   void fill(const StEvent*);
00118 
00119   // classes that we just takes from StEvent
00120   StRunInfo mRunInfo;
00121   StEventInfo mEventInfo;
00122   StEventSummary mEventSummary;
00123   StVpdTriggerDetector mVpdTriggerDetector;
00124   StMtdTriggerDetector mMtdTriggerDetector;
00125   StCtbTriggerDetector mCtbTriggerDetector;
00126   StZdcTriggerDetector mZdcTriggerDetector;
00127   StBbcTriggerDetector mBbcTriggerDetector;
00128   StEmcTriggerDetector mEmcTriggerDetector;
00129   StFpdTriggerDetector mFpdTriggerDetector;
00130   StFmsTriggerDetector mFmsTriggerDetector;
00131   StFpdCollection mFpdCollection; 
00132   StL0Trigger mL0Trigger; 
00133   // special classes from MuDst
00134   StMuL3EventSummary mL3EventSummary;
00135   StMuTriggerIdCollection mTriggerIdCollection;
00136   StTriggerData* mTriggerData;
00137 
00138   UShort_t mRefMultPos;
00139   UShort_t mRefMultNeg;
00140   UShort_t mRefMultFtpcEast;
00141   UShort_t mRefMultFtpcWest;
00142   Float_t mReactionPlane[2];              
00143   Float_t mReactionPlanePtWgt[2];
00144 
00145   StThreeVectorF mPrimaryVertexError;
00146 
00147   TArrayI mL2Result; // Raw L2 info
00148 
00149   UInt_t  mVpdEast;    // VPD East Hit pattern 
00150   UInt_t  mVpdWest;    // VPD West Hit pattern
00151   Float_t mVpdTstart;  // VPD start time (calibrated)
00152   Float_t mVpdTdiff;   // VPD time difference (calibrated)
00153   Float_t mVpdVz;      // VPD vertex z
00154 
00155   friend class StMuDst;
00156   friend class StMuDstMaker;
00157   friend class StMuMomentumShiftMaker;
00158   friend class StMuL3EventSummary;
00159   ClassDef(StMuEvent,13)
00160 };
00161 
00162 inline int StMuEvent::eventId() { return mEventInfo.id();}
00163 inline int StMuEvent::eventNumber() { return mEventInfo.id();}
00164 inline int StMuEvent::runId() { return mEventInfo.runId();}
00165 inline int StMuEvent::runNumber() { return mEventInfo.runId();}
00166 inline StRunInfo& StMuEvent::runInfo() {return mRunInfo;}
00167 inline StEventInfo& StMuEvent::eventInfo() {return mEventInfo;}
00168 inline StEventSummary& StMuEvent::eventSummary() {return mEventSummary;}
00169 inline StVpdTriggerDetector& StMuEvent::vpdTriggerDetector() {return mVpdTriggerDetector;}
00170 inline StMtdTriggerDetector& StMuEvent::mtdTriggerDetector() {return mMtdTriggerDetector;}
00171 inline StCtbTriggerDetector& StMuEvent::ctbTriggerDetector() {return mCtbTriggerDetector;}
00172 inline StZdcTriggerDetector& StMuEvent::zdcTriggerDetector() {return mZdcTriggerDetector;}
00173 inline StBbcTriggerDetector& StMuEvent::bbcTriggerDetector() {return mBbcTriggerDetector;}
00174 inline StEmcTriggerDetector& StMuEvent::emcTriggerDetector() {return mEmcTriggerDetector;}
00175 inline StFpdTriggerDetector& StMuEvent::fpdTriggerDetector() {return mFpdTriggerDetector;}
00176 inline StFmsTriggerDetector& StMuEvent::fmsTriggerDetector() {return mFmsTriggerDetector;}
00177 inline StFpdCollection& StMuEvent::fpdCollection() {return mFpdCollection;} 
00178 inline StL0Trigger& StMuEvent::l0Trigger() {return mL0Trigger;} 
00179 // special classes for muDst
00180 inline StMuL3EventSummary& StMuEvent::l3EventSummary() {return mL3EventSummary;}
00181 inline StMuTriggerIdCollection& StMuEvent::triggerIdCollection(){return mTriggerIdCollection;}
00182 //inline const StTriggerData* StMuEvent::triggerData() const { if(mTriggerData!=0) {mTriggerData->setDebug(0); return mTriggerData; } else return 0; }
00183 inline const StTriggerData* StMuEvent::triggerData() const { return mTriggerData; }
00184 inline double StMuEvent::reactionPlane(unsigned short s) {return (s==0) ? mReactionPlane[0] : mReactionPlane[1];}
00185 inline void StMuEvent::setReactionPlane(unsigned short s, double v) {(s==0) ? mReactionPlane[0]=v : mReactionPlane[1]=v;}
00186 inline double StMuEvent::reactionPlanePtWgt(unsigned short s) {return (s==0) ? mReactionPlanePtWgt[0] : mReactionPlanePtWgt[1];}
00187 inline void StMuEvent::setReactionPlanePtWgt(unsigned short s, double v) {(s==0) ? mReactionPlanePtWgt[0]=v : mReactionPlanePtWgt[1]=v;}
00188 inline double StMuEvent::magneticField() { return mEventSummary.magneticField();}
00189 inline double StMuEvent::zdcAdcAttentuatedSumWest() { return mZdcTriggerDetector.adc(10);}
00190 inline double StMuEvent::zdcAdcAttentuatedSumEast() { return mZdcTriggerDetector.adc(13);}
00191 inline double StMuEvent::ctbMultiplicity() { 
00192   double ctb=0;
00193   for (unsigned int slat = 0; slat < mCtbTriggerDetector.numberOfSlats(); slat++) {
00194     for (unsigned int tray = 0; tray < mCtbTriggerDetector.numberOfTrays(); tray++) {
00195       ctb += mCtbTriggerDetector.mips(tray,slat,0);
00196     }
00197   }
00198   return ctb;
00199 }
00200 inline TArrayI &StMuEvent::L2Result() { return mL2Result; }
00201 inline unsigned int StMuEvent::numberOfVpdEastHits() {  
00202   unsigned int num = 0;
00203   for(int i=0;i<32;i++) {
00204     num += mVpdEast>>i & 1;
00205   }
00206   return num;
00207 }
00208 inline unsigned int StMuEvent::numberOfVpdWestHits() {
00209   unsigned int num = 0;
00210   for(int i=0;i<32;i++) {
00211     num += mVpdWest>>i & 1;
00212   }
00213   return num;
00214 }
00215 inline float StMuEvent::vpdTstart() { return mVpdTstart; }
00216 inline float StMuEvent::vpdTdiff() { return mVpdTdiff; }
00217 inline float StMuEvent::vpdVz() { return mVpdVz; }
00218 #endif
00219 /***************************************************************************
00220  *
00221  * $Log: StMuEvent.h,v $
00222  * Revision 1.30  2010/05/28 19:47:51  tone421
00223  * Removed a cout needed for test purposes in StMuDstMaker. Made sure StTriggerData objects copied into the MuDst have a debug value of 0..
00224  *
00225  * Revision 1.29  2010/05/26 17:34:59  tone421
00226  * Added const protection to StTriggerData* StMuEvent::triggerData()
00227  *
00228  * Revision 1.27  2010/02/03 17:16:22  tone421
00229  * Added function StMuEvent::nearestVertexZ(int vtx_id) which returns the z distance of the nearest vertex in relation to vertex vtx_id
00230  *
00231  * Revision 1.26  2010/02/03 04:54:45  tone421
00232  * Added StMuEvent::btofTrayMultiplicity() to return only TOF hits from trays. Should be looked at instead of ctbSum for run 9 and beyond.
00233  *
00234  * Revision 1.25  2010/01/25 03:57:39  tone421
00235  * Added FMS and Roman pot arrays
00236  *
00237  * Revision 1.24  2009/08/25 15:49:49  fine
00238  * fix the compilation issues under SL5_64_bits  gcc 4.3.2
00239  *
00240  * Revision 1.23  2009/01/09 19:43:47  tone421
00241  * OAdded gremult in StMuEvent (globals tracks with DCA < 3cm, >= 10 TPC fit hits and |eta| < 0.5)
00242  *
00243  * Revision 1.22  2008/11/18 15:34:33  tone421
00244  * 2 changes. The first ensures StMuEvent::primaryVertexPosition() returns the position of current vertex (set by StMuDst::setVertexIndex(Int_t vtx_id)): previously it returned the position of best ranked vertex. The second insures events with no vertex have a PVx=PYy=PYz=-999 rather than 0.
00245  *
00246  * Revision 1.21  2008/06/26 15:41:29  tone421
00247  *
00248  * Add getter and setter for vpd z vertex position
00249  *
00250  * Revision 1.20  2008/02/20 09:00:48  mvl
00251  * Included FMS data (StFMSTriggerDetector) (code by Akio)
00252  *
00253  * Revision 1.19  2007/09/21 02:27:12  mvl
00254  * Added calibrated VPD info from StTofCollection (run-8 prep)
00255  *
00256  * Revision 1.18  2007/09/05 23:21:21  mvl
00257  * Added StMtdTriggerDetector
00258  *
00259  * Revision 1.17  2007/08/02 20:46:46  mvl
00260  * Switch off Q-vector branhces in StMuDstMaker and increase version number in StMuEvent.
00261  * This is to avoid wranings when reading P07ib data which has Q-vector information stored with more recent libraries.
00262  *
00263  * Revision 1.16  2007/04/20 06:25:21  mvl
00264  * Removed Q-vectors (will implement utility class).
00265  * Added Vpd info.
00266  *
00267  * Revision 1.14  2006/09/20 01:50:35  mvl
00268  * Added data member and code for L2Result array (TArrayI).
00269  *
00270  * Revision 1.13  2005/08/19 19:46:05  mvl
00271  * Further updates for multiple vertices. The main changes are:
00272  * 1) StMudst::primaryTracks() now returns a list (TObjArray*) of tracks
00273  *    belonging to the 'current' primary vertex. The index number of the
00274  *    'current' vertex can be set using StMuDst::setCurrentVertex().
00275  *    This also affects StMuDst::primaryTracks(int i) and
00276  *    StMuDst::numberOfprimaryTracks().
00277  * 2) refMult is now stored for all vertices, in StMuPrimaryVertex. The
00278  *    obvious way to access these numbers is from the StMuprimaryVertex structures,
00279  *    but for ebakcward compatibility a function is provided in StMuEvent as well
00280  *    (this is the only function taht works for existing MuDst)
00281  *
00282  * As an aside, I've also changes the internals of StMuDst::createStEvent and
00283  * StMuDst::fixTrackIndices() to be able to deal with a larger range of index numbers for tracks as generated by Ittf.
00284  *
00285  * BIG FAT WARNING: StMudst2StEventMaker and StMuDstFilterMaker
00286  * do not fully support the multiple vertex functionality yet.
00287  *
00288  * Revision 1.12  2005/03/17 21:55:00  mvl
00289  * Added StMuMomentumShiftMaker for applying a magnetic field scaling to the reconstructed MuDst. This class accesses StMuTrack, StMuEvent and StMuHelix and some Strangeness MuDst data members as 'friend'
00290  *
00291  * Revision 1.11  2004/12/02 00:19:52  mvl
00292  * Added error on primary vertex
00293  *
00294  * Revision 1.10  2004/08/04 17:57:13  mvl
00295  * Added EMC trigger information and fpd trigger (tower) information
00296  *
00297  * Revision 1.9  2004/05/02 04:10:14  perev
00298  * private => protected
00299  *
00300  * Revision 1.8  2003/07/22 19:14:40  laue
00301  * multiplicities for FTPC added
00302  *
00303  * Revision 1.7  2003/02/20 15:29:42  laue
00304  * StMuTriggerIdCollection added
00305  *
00306  * Revision 1.6  2003/02/19 13:52:11  laue
00307  * added the StTriggerIdCollection
00308  *
00309  * Revision 1.5  2003/01/09 18:59:46  laue
00310  * initial check in of new EMC classes and the changes required
00311  *
00312  * Revision 1.4  2002/08/23 17:30:18  laue
00313  * additional member functions added (Helen Caines' request)
00314  *
00315  * Revision 1.3  2002/08/20 19:55:49  laue
00316  * Doxygen comments added
00317  *
00318  * Revision 1.2  2002/03/20 16:04:11  laue
00319  * minor changes, mostly added access functions
00320  *
00321  * Revision 1.1  2002/03/08 17:04:17  laue
00322  * initial revision
00323  *
00324  *
00325  **************************************************************************/