StRoot: StHbtMaker/Infrastructure/StHbtParticle.hh Source File

00001 /***************************************************************************
00002  *
00003  * $Id: StHbtParticle.hh,v 1.19 2003/01/14 09:41:16 renault Exp $
00004  *
00005  * Author: Mike Lisa, Ohio State, lisa@mps.ohio-state.edu
00006  ***************************************************************************
00007  *
00008  * Description: part of STAR HBT Framework: StHbtMaker package
00009  *   Particle objects are part of the PicoEvent, which is what is
00010  *   stored in the EventMixingBuffers
00011  *   A Track object gets converted to a Particle object if it
00012  *   passes the ParticleCut of an Analysis
00013  *
00014  ***************************************************************************
00015  *
00016  * $Log: StHbtParticle.hh,v $
00017  * Revision 1.19  2003/01/14 09:41:16  renault
00018  * changes on average separation calculation, hit shared finder and memory optimisation
00019  * for Z,U and Sectors variables.
00020  *
00021  * Revision 1.18  2002/12/12 17:01:50  kisiel
00022  * Hidden Information handling and purity calculation
00023  *
00024  * Revision 1.17  2002/11/19 23:35:52  renault
00025  * Enable calculation of exit/entrance separation for V0 daughters
00026  *
00027  * Revision 1.16  2001/12/14 23:11:30  fretiere
00028  * Add class HitMergingCut. Add class fabricesPairCut = HitMerginCut + pair purity cuts. Add TpcLocalTransform function which convert to local tpc coord (not pretty). Modify StHbtTrack, StHbtParticle, StHbtHiddenInfo, StHbtPair to handle the hit information and cope with my code
00029  *
00030  * Revision 1.15  2001/05/25 23:23:59  lisa
00031  * Added in StHbtKink stuff
00032  *
00033  * Revision 1.14  2001/05/23 00:19:05  lisa
00034  * Add in Smearing classes and methods needed for momentum resolution studies and correction
00035  *
00036  * Revision 1.13  2001/04/03 21:04:36  kisiel
00037  *
00038  *
00039  *   Changes needed to make the Theoretical code
00040  *   work. The main code is the ThCorrFctn directory.
00041  *   The most visible change is the addition of the
00042  *   HiddenInfo to StHbtPair.
00043  *
00044  * Revision 1.12  2000/10/05 23:09:05  lisa
00045  * Added kT-dependent radii to mixed-event simulator AND implemented AverageSeparation Cut and CorrFctn
00046  *
00047  * Revision 1.11  2000/07/17 20:03:17  lisa
00048  * Implemented tools for addressing and assessing trackmerging
00049  *
00050  * Revision 1.10  2000/07/16 21:38:23  laue
00051  * StHbtCoulomb.cxx StHbtSectoredAnalysis.cxx : updated for standalone version
00052  * StHbtV0.cc StHbtV0.hh : some cast to prevent compiling warnings
00053  * StHbtParticle.cc StHbtParticle.hh : pointers mTrack,mV0 initialized to 0
00054  * StHbtIOBinary.cc : some printouts in #ifdef STHBTDEBUG
00055  * StHbtEvent.cc : B-Field set to 0.25Tesla, we have to think about a better
00056  *                 solution
00057  *
00058  * Revision 1.9  2000/05/03 17:44:43  laue
00059  * StHbtEvent, StHbtTrack & StHbtV0 declared friend to StHbtIOBinary
00060  * StHbtParticle updated for V0 pos,neg track Id
00061  *
00062  * Revision 1.8  2000/04/03 16:21:51  laue
00063  * some include files changed
00064  * Multi track cut added
00065  *
00066  * Revision 1.6  2000/02/26 19:04:52  laue
00067  * Some unnecessary includes removed.
00068  * StThreeVectorD replace by StHbtThreeVector.
00069  * StHbtCoulomb modified to compile without Root (ClassDef embraced into
00070  *   #ifdef __ROOT__  ..... #endif)
00071  * StHbtParticle now returns references (FourMomentum(),Helix(),
00072  *   DecayVertexPosiion())
00073  *
00074  * Revision 1.5  1999/12/11 15:58:29  lisa
00075  * Add vertex decay position datum and accessor to StHbtParticle to allow pairwise cuts on seperation of V0s
00076  *
00077  * Revision 1.4  1999/09/17 22:38:02  lisa
00078  * first full integration of V0s into StHbt framework
00079  *
00080  * Revision 1.3  1999/09/01 19:04:54  lisa
00081  * update Particle class AND add parity cf and Randys Coulomb correction
00082  *
00083  * Revision 1.2  1999/07/06 22:33:23  lisa
00084  * Adjusted all to work in pro and new - dev itself is broken
00085  *
00086  * Revision 1.1.1.1  1999/06/29 16:02:57  lisa
00087  * Installation of StHbtMaker
00088  *
00089  **************************************************************************/
00090 
00091 #ifndef StHbtParticle_hh
00092 #define StHbtParticle_hh
00093 
00094 #include "StHbtMaker/Infrastructure/StHbtTypes.hh"
00095 #include "StHbtMaker/Infrastructure/StHbtTrack.hh"
00096 #include "StHbtMaker/Infrastructure/StHbtV0.hh"
00097 #include "StHbtMaker/Infrastructure/StHbtKink.hh"
00098 #include "StHbtMaker/Infrastructure/StHbtXi.hh"
00099 #include "StPhysicalHelixD.hh"
00100 // ***
00101 class StHbtHiddenInfo;
00102 // ***
00103 class StHbtParticle{
00104 public:
00105   StHbtParticle();
00106   StHbtParticle(const StHbtTrack* const hbtTrack, const double& mass);
00107   StHbtParticle(const StHbtV0* const hbtV0, const double& mass);
00108   StHbtParticle(const StHbtKink* const hbtKink, const double& mass);
00109   StHbtParticle(const StHbtXi* const hbtXi, const double& mass);
00110   ~StHbtParticle();
00111 
00112   const StHbtLorentzVector& FourMomentum() const;
00113 
00114   StPhysicalHelixD& Helix();
00115 
00116   const StHbtThreeVector DecayVertexPosition() const;
00117 
00118   unsigned long TopologyMap(const int word) const;
00119   int NumberOfHits() const;
00120 
00121   unsigned long TrackId() const;         // only for particles from tracks 
00122   unsigned short   NegTrackId() const;   // only for particles from v0 
00123   unsigned short   PosTrackId() const;   // only for particles from v0 
00124 
00125   StHbtTrack* Track() const;
00126   StHbtV0* V0() const;
00127   StHbtKink* Kink() const;
00128 
00129   const StHbtThreeVector& NominalTpcExitPoint() const;     // position track exits TPC assuming start at (0,0,0)
00130   const StHbtThreeVector& NominalTpcEntrancePoint() const; // position track crosses IFC assuming start at (0,0,0)
00131   const StHbtThreeVector& TpcV0PosExitPoint() const;  
00132   const StHbtThreeVector& TpcV0PosEntrancePoint() const;
00133   const StHbtThreeVector& TpcV0NegExitPoint() const;  
00134   const StHbtThreeVector& TpcV0NegEntrancePoint() const;
00135 
00136   // the following method is for explicit internal calculation to fill datamembers.
00137   // It is invoked automatically if StHbtParticle constructed from StHbtTrack
00138   //void CalculateNominalTpcExitAndEntrancePoints(); // NOTE - this requires the mHelix, so be sure this is filled
00139 
00140 
00141   StHbtThreeVector mNominalPosSample[11];  // I make this public for convenience and speed of StHbtPair()
00142   float mZ[45];
00143   float mU[45];
00144   int mSect[45];
00145 
00146   void ResetFourMomentum(const StHbtLorentzVector& fourMomentum);
00147 
00148   const StHbtHiddenInfo*  HiddenInfo() const;
00149   // Fab private
00150   StHbtHiddenInfo*  getHiddenInfo() const;
00151   void SetHiddenInfo(StHbtHiddenInfo* aHiddenInfo);
00152   void CalculatePurity();
00153   double GetPionPurity();
00154   double GetKaonPurity();
00155   double GetProtonPurity();
00156   void CalculateTpcExitAndEntrancePoints( StPhysicalHelixD* tHelix,
00157                                           StHbtThreeVector* PrimVert,
00158                                           StHbtThreeVector* SecVert,
00159                                           StHbtThreeVector* tmpTpcEntrancePoint,
00160                                           StHbtThreeVector* tmpTpcExitPoint,
00161                                           StHbtThreeVector* tmpPosSample,
00162                                           float* tmpZ,float* tmpU,int* tmpSect);
00163 
00164   // For V0 Neg Daugthers TpcEntrance/ExitPoints
00165   StHbtThreeVector* mTpcV0NegPosSample;
00166   float* mV0NegZ;
00167   float* mV0NegU;
00168   int* mV0NegSect;
00169  
00170 private:
00171   StHbtTrack* mTrack;  // copy of the track the particle was formed of, else Null
00172   StHbtV0* mV0;        // copy of the v0 the particle was formed of, else Null
00173   StHbtKink* mKink;        // copy of the v0 the particle was formed of, else Null
00174   StHbtXi* mXi;
00175 
00176   StHbtLorentzVector mFourMomentum;
00177   StPhysicalHelixD mHelix;
00178   unsigned long  mMap[2]; 
00179   int mNhits;
00180   StHbtThreeVector mNominalTpcExitPoint;
00181   StHbtThreeVector mNominalTpcEntrancePoint;
00182   StHbtHiddenInfo* mHiddenInfo;  // Fab private
00183 
00184   double mPurity[6];
00185 
00186   static double mPrimPimPar0;
00187   static double mPrimPimPar1;
00188   static double mPrimPimPar2;
00189   static double mPrimPipPar0;
00190   static double mPrimPipPar1;
00191   static double mPrimPipPar2;
00192   static double mPrimPmPar0;
00193   static double mPrimPmPar1;
00194   static double mPrimPmPar2;
00195   static double mPrimPpPar0;
00196   static double mPrimPpPar1;
00197   static double mPrimPpPar2;
00198 
00199    // For V0 Daugthers TpcEntrance/ExitPoints
00200   StHbtThreeVector mPrimaryVertex;
00201   StHbtThreeVector mSecondaryVertex;
00202 
00203   StPhysicalHelixD mHelixV0Pos;
00204   StHbtThreeVector mTpcV0PosEntrancePoint;
00205   StHbtThreeVector mTpcV0PosExitPoint;
00206 
00207   StPhysicalHelixD mHelixV0Neg;
00208   StHbtThreeVector mTpcV0NegEntrancePoint;
00209   StHbtThreeVector mTpcV0NegExitPoint;
00210 };
00211 
00212 inline StHbtTrack* StHbtParticle::Track() const { return mTrack; }
00213 inline unsigned long  StHbtParticle::TrackId() const { return mTrack->TrackId(); }; 
00214 inline const StHbtLorentzVector& StHbtParticle::FourMomentum() const {return mFourMomentum;}
00215 inline StPhysicalHelixD& StHbtParticle::Helix() {return mHelix;}
00216 inline unsigned long StHbtParticle::TopologyMap(const int word) const {return mMap[word];}
00217 inline int StHbtParticle::NumberOfHits() const {return mNhits;}
00218 inline StHbtV0* StHbtParticle::V0() const { return mV0; }
00219 inline unsigned short StHbtParticle::NegTrackId() const { return mV0->idNeg(); }
00220 inline unsigned short StHbtParticle::PosTrackId() const { return mV0->idPos(); }
00221 inline const StHbtThreeVector StHbtParticle::DecayVertexPosition() const {return mV0->decayVertexV0(); }
00222 // ***
00223 inline StHbtHiddenInfo* StHbtParticle::getHiddenInfo() const
00224 {return mHiddenInfo;}
00225 inline const StHbtHiddenInfo* StHbtParticle::HiddenInfo() const
00226 {return mHiddenInfo;}
00227 inline void StHbtParticle::SetHiddenInfo(StHbtHiddenInfo* aHiddenInfo)
00228 { mHiddenInfo = aHiddenInfo->clone();}
00229 // ***
00230 
00231 inline void StHbtParticle::ResetFourMomentum(const StHbtLorentzVector& vec){mFourMomentum = vec;}
00232 
00233 inline StHbtKink* StHbtParticle::Kink() const { return mKink; }
00234 
00235 #endif