StRoot: StStrangeMuDstMaker/StStrangeEvMuDst.hh Source File

00001 
00010 #ifndef StStrangeEvMuDst_hh
00011 #define StStrangeEvMuDst_hh
00012 #include "StStrangeMuDst.hh"
00013 #include "TMath.h"
00014 #include "stdlib.h"
00015 
00016 
00017 class StEvent;
00018 class StMcEvent;
00019 
00020 class StStrangeEvMuDst : public StStrangeMuDst {
00021   friend class StMuMomentumShiftMaker;
00022 
00023 public:
00024   StStrangeEvMuDst();
00025   ~StStrangeEvMuDst();
00026   StStrangeEvMuDst(StEvent&);
00027   StStrangeEvMuDst(StMcEvent&);
00028   void Clear();
00029   void Clear(Option_t* option) {TObject::Clear(option);}
00030 
00031   static void SetCorrectionFile(char*); 
00032   static void SetFractionFile(char*);
00033 
00035 
00036 
00037   Int_t   run()                    const;
00039   Int_t   event()                  const;
00041   Bool_t  unbiasedTrigger()        const;
00043   UInt_t  l0TriggerWord()          const;
00045   Bool_t  ftpcGoodRun();
00046   Float_t primaryVertexX()         const;
00047   Float_t primaryVertexY()         const; 
00048   Float_t primaryVertexZ()         const;
00050   Float_t magneticField()          const;
00052 
00054 
00055   Int_t   globalTracks()           const;
00056   Int_t   tpcTracks()              const;
00057   Int_t   ftpcEastTracks()         const;
00058   Int_t   ftpcWestTracks()         const;
00059   Int_t   primaryTracks()          const;
00060   Int_t   tpcPrimTracks()          const;
00061   Int_t   ftpcEastPrimTracks()     const;
00062   Int_t   ftpcWestPrimTracks()     const;
00063   Int_t   primaryNegTracks()       const;
00064   Float_t primaryCorrectedTracks() const;
00066   UInt_t  zdcWestADC()             const;
00068   Float_t fractionSigma()          const;
00070 
00072 
00073   void Fill(StEvent&);
00074   void Fill(StMcEvent&);
00075   void setMagneticField(Float_t);
00076   void setL0TriggerWord(UInt_t);
00077   void setGlobalTracks(Int_t);
00078   void setPrimaryTracks(Int_t);
00080   void tagRunNumber(Int_t);
00082 
00083 
00084 protected:
00085                                   // These are written out
00087   Int_t   mRun;
00089   Int_t   mEvent;
00090   Float_t mPrimaryVertexX;
00091   Float_t mPrimaryVertexY;
00092   Float_t mPrimaryVertexZ;
00093   Int_t   mGlobalTracks;
00094   Int_t   mPrimaryTracks;
00095   Int_t   mPrimaryNegTracks;
00096   Float_t mMagneticField;
00097   UInt_t  mL0TriggerWord;
00098 
00099 private:
00100 
00101   Int_t FTPC(const Int_t) const;
00102 
00103   ClassDef(StStrangeEvMuDst,8)
00104 };
00105 
00106 inline         StStrangeEvMuDst::StStrangeEvMuDst(StEvent& event)
00107                { Fill(event); }
00108 inline         StStrangeEvMuDst::StStrangeEvMuDst(StMcEvent& event)
00109                { Fill(event); }
00110 
00111 inline void    StStrangeEvMuDst::tagRunNumber(Int_t r)
00112                { mRun = - TMath::Abs(r); }
00113 inline Bool_t  StStrangeEvMuDst::ftpcGoodRun()
00114                { return (mRun > 0); }
00115 
00116 inline Int_t   StStrangeEvMuDst::run() const
00117                { return TMath::Abs(mRun); }
00118 inline Int_t   StStrangeEvMuDst::event() const
00119                { return TMath::Abs(mEvent); }
00120 inline Bool_t  StStrangeEvMuDst::unbiasedTrigger() const
00121                { return (mEvent > 0); }
00122 inline UInt_t  StStrangeEvMuDst::l0TriggerWord() const
00123                { return ((run() < 3999999) ? mL0TriggerWord :
00124                                              (UInt_t) mL0TriggerWord%100000); }
00125 inline UInt_t  StStrangeEvMuDst::zdcWestADC() const
00126                { return ((run() < 3999999) ? 0U :
00127                                              (UInt_t) mL0TriggerWord/100000); }
00128 inline Float_t StStrangeEvMuDst::magneticField() const 
00129                { return mMagneticField; }
00130 
00131 inline Float_t StStrangeEvMuDst::primaryVertexX() const 
00132                { return mPrimaryVertexX; }
00133 inline Float_t StStrangeEvMuDst::primaryVertexY() const 
00134                { return mPrimaryVertexY; }
00135 inline Float_t StStrangeEvMuDst::primaryVertexZ() const 
00136                { return mPrimaryVertexZ; }
00137 
00138 // Store FTPCwest*724 + FTPCeast in most significant 19 bits (max = 723 each)
00139 // Store main TPC in least significant 13 bits (max = 8191)
00140 inline Int_t   StStrangeEvMuDst::FTPC(const Int_t x) const
00141                { return ((Int_t) (((UInt_t) x)>>13)); }
00142 
00143 inline Int_t   StStrangeEvMuDst::globalTracks() const
00144                { div_t ftpc = div(FTPC(mGlobalTracks),724);
00145                  return (tpcTracks() + ftpc.quot + ftpc.rem); }
00146 inline Int_t   StStrangeEvMuDst::tpcTracks() const
00147                { return mGlobalTracks & 0x1fff; }
00148 inline Int_t   StStrangeEvMuDst::ftpcEastTracks() const
00149                { return FTPC(mGlobalTracks)%724; }
00150 inline Int_t   StStrangeEvMuDst::ftpcWestTracks() const
00151                { return FTPC(mGlobalTracks)/724; }
00152 
00153 inline Int_t   StStrangeEvMuDst::primaryTracks() const
00154                { div_t ftpc = div(FTPC(mPrimaryTracks),724);
00155                  return (tpcPrimTracks() + ftpc.quot + ftpc.rem); }
00156 inline Int_t   StStrangeEvMuDst::tpcPrimTracks() const
00157                { return mPrimaryTracks & 0x1fff; }
00158 inline Int_t   StStrangeEvMuDst::ftpcEastPrimTracks() const
00159                { return FTPC(mPrimaryTracks)%724; }
00160 inline Int_t   StStrangeEvMuDst::ftpcWestPrimTracks() const
00161                { return FTPC(mPrimaryTracks)/724; }
00162 inline Int_t   StStrangeEvMuDst::primaryNegTracks() const
00163                { return mPrimaryNegTracks; }
00164 
00165 inline void    StStrangeEvMuDst::setMagneticField(Float_t b)
00166                { mMagneticField = b; }
00167 inline void    StStrangeEvMuDst::setL0TriggerWord(UInt_t trgWrd)
00168                { mL0TriggerWord = trgWrd; }
00169 inline void    StStrangeEvMuDst::setGlobalTracks(Int_t newGlob)
00170                { mGlobalTracks = newGlob;}
00171 inline void    StStrangeEvMuDst::setPrimaryTracks(Int_t newPrim)
00172                { mPrimaryTracks = newPrim;}
00173 
00174 #endif
00175 
00176 
00177 /***********************************************************************
00178  * $Id: StStrangeEvMuDst.hh,v 3.9 2005/03/17 05:06:51 genevb Exp $
00179  * $Log: StStrangeEvMuDst.hh,v $
00180  * Revision 3.9  2005/03/17 05:06:51  genevb
00181  * Add StMuMomentumShiftMaker friend
00182  *
00183  * Revision 3.8  2003/05/30 21:20:19  genevb
00184  * doxygen savvy, encoding of FTPC mults, change virtual funcs
00185  *
00186  * Revision 3.7  2003/02/10 16:00:29  genevb
00187  * Implement cleared events
00188  *
00189  * Revision 3.6  2002/05/17 14:05:28  genevb
00190  * Added L3 unbiased trigger info
00191  *
00192  * Revision 3.5  2002/04/30 16:02:48  genevb
00193  * Common muDst, improved MC code, better kinks, StrangeCuts now a branch
00194  *
00195  * Revision 3.4  2001/11/05 23:41:06  genevb
00196  * Add more dEdx, B field info, careful of changes to TTree unrolling
00197  *
00198  * Revision 3.3  2001/09/14 21:22:26  genevb
00199  * Avoid hiding TObject::Clear()
00200  *
00201  * Revision 3.2  2001/08/23 13:20:54  genevb
00202  * Many bug workarounds...
00203  *
00204  * Revision 3.1  2001/05/04 20:15:14  genevb
00205  * Common interfaces and reorganization of components, add MC event info
00206  *
00207  * Revision 3.0  2000/07/14 12:56:49  genevb
00208  * Revision 3 has event multiplicities and dedx information for vertex tracks
00209  *
00210  * Revision 2.1  2000/06/09 22:17:10  genevb
00211  * Allow MC data to be copied between DSTs, other small improvements
00212  *
00213  * Revision 2.0  2000/06/02 22:11:54  genevb
00214  * New version of Strangeness micro DST package
00215  *
00216  * Revision 1.3  2000/03/31 03:20:24  jones
00217  * Added topology map to V0/Xi; access funcs for each data member
00218  *
00219  * Revision 1.2  2000/03/29 20:52:13  genevb
00220  * Added StKinkMuDst, replaced arrays
00221  *
00222  * Revision 1.1  2000/03/29 03:10:07  genevb
00223  * Int_troduction of Strangeness Micro DST package
00224  *
00225  ***********************************************************************/