StRoot: StFlowMaker/StFlowCutTrack.h Source File

00001 
00002 //
00003 // $Id: StFlowCutTrack.h,v 1.20 2003/09/02 17:58:11 perev Exp $
00004 //
00005 // Author: Art Poskanzer and Raimond Snellings, LBNL, Nov 1999
00006 //          FTPC added by Markus Oldenburg, MPI, Dec 2000
00007 //          MuDst enabled by Kirill Filimonov, LBNL, Jun 2002
00008 //
00009 // Description:  Class for applying flow track cuts
00010 //               If lo >= hi no cut is applied
00011 //               All functions and data members are static
00012 //               Therefor, no need to instantiate
00013 //               Just use StFlowCutTrack::func();
00014 //
00016 
00017 #ifndef _StFlowCutTrack_INCLUDED_
00018 #define _StFlowCutTrack_INCLUDED_
00019 #include <Stiostream.h>
00020 #include <stdlib.h>
00021 #include "Rtypes.h"
00022 class StTrack;
00023 class StFlowPicoTrack;
00024 class StMuTrack;
00025 
00026 class StFlowCutTrack {
00027 
00028  public:
00029 
00030                  StFlowCutTrack();
00031   virtual        ~StFlowCutTrack();
00032 
00033   static Int_t   CheckTrack(StTrack* pTrack);
00034   static Int_t   CheckTrack(StFlowPicoTrack* pPicoTrack);
00035   static Int_t   CheckTrack(StMuTrack* pMuTrack);
00036   static void    PrintCutList();
00037   static UInt_t  EtaSymPosTpc();
00038   static UInt_t  EtaSymNegTpc();
00039   static UInt_t  EtaSymPosFtpc();
00040   static UInt_t  EtaSymNegFtpc();
00041   static void    EtaSymClear();
00042   static void    SetFitPtsTpc(Int_t lo, Int_t hi);
00043   static void    SetFitPtsFtpc(Int_t lo, Int_t hi);
00044   static void    SetFitOverMaxPts(Float_t lo, Float_t hi);
00045   static void    SetChiSqTpc(Float_t lo, Float_t hi);
00046   static void    SetChiSqFtpc(Float_t lo, Float_t hi);
00047   static void    SetDcaFtpc(Float_t lo, Float_t hi);
00048   static void    SetDcaGlobalFtpc(Float_t lo, Float_t hi);
00049   static void    SetPtTpc(Float_t lo, Float_t hi);
00050   static void    SetPtFtpc(Float_t lo, Float_t hi);
00051   static void    SetEtaTpc(Float_t lo, Float_t hi);
00052   static void    SetEtaFtpc(Float_t lo_neg, Float_t hi_neg, Float_t lo_pos, Float_t hi_pos);
00053   static void    SetChgTpc(Int_t lo, Int_t hi);
00054   static void    SetChgFtpc(Int_t lo, Int_t hi);
00055   static void    IncludeTpcTracks(Bool_t tpc_include);
00056   static void    IncludeFtpcTracks(Bool_t ftpc_include);
00057 
00058  private:
00059 
00060   static UInt_t  mTrackN;                    // number of tracks
00061   static UInt_t  mTpcTrackN;                 // number of Tpc tracks
00062   static UInt_t  mFtpcTrackN;                // number of Ftpc tracks
00063   static UInt_t  mFtpcEastTrackN;            // number of Ftpc east tracks
00064   static UInt_t  mFtpcWestTrackN;            // number of Ftpc west tracks
00065 
00066   static UInt_t  mGoodTrackN;                // number of accepted tracks   
00067   static UInt_t  mGoodTpcTrackN;             // number of accepted Tpc tracks   
00068   static UInt_t  mGoodFtpcTrackN;            // number of accepted Ftpc tracks
00069    
00070   static UInt_t  mEtaSymPosTpcN;             // number of positive Eta Tpc tracks
00071   static UInt_t  mEtaSymNegTpcN;             // number of negative Eta Tpc tracks
00072   static UInt_t  mEtaSymPosFtpcN;            // number of positive Eta Ftpc tracks
00073   static UInt_t  mEtaSymNegFtpcN;            // number of negative Eta Ftpc tracks
00074                                                 
00075   static UInt_t  mFitPtsTpcCutN;             // number not accepted
00076   static Int_t   mFitPtsTpcCuts[2];          // range
00077 
00078   static UInt_t  mFitPtsFtpcCutN;            // number not accepted
00079   static Int_t   mFitPtsFtpcCuts[2];         // range
00080 
00081   static UInt_t  mFitOverMaxCutN;            // number not accepted (all)
00082   static UInt_t  mFitOverMaxTpcCutN;         // number not accepted (Tpc)
00083   static UInt_t  mFitOverMaxFtpcCutN;        // number not accepted (Ftpc)
00084   static Float_t mFitOverMaxCuts[2];         // range
00085 
00086   static UInt_t  mChiSqTpcCutN;              // number not accepted
00087   static Float_t mChiSqTpcCuts[2];           // range
00088 
00089   static UInt_t  mChiSqFtpcCutN;             // number not accepted
00090   static Float_t mChiSqFtpcCuts[2];          // range
00091 
00092   static UInt_t  mDcaFtpcCutN;               // number not accepted
00093   static Float_t mDcaFtpcCuts[2];            // range
00094 
00095   static UInt_t  mDcaGlobalFtpcCutN;         // number not accepted
00096   static Float_t mDcaGlobalFtpcCuts[2];      // range
00097 
00098   static UInt_t  mPtTpcCutN;                 // number not accepted
00099   static Float_t mPtTpcCuts[2];              // range
00100 
00101   static UInt_t  mPtFtpcCutN;                // number not accepted
00102   static Float_t mPtFtpcCuts[2];             // range
00103 
00104   static UInt_t  mEtaTpcCutN;                // number not accepted
00105   static Float_t mEtaTpcCuts[2];             // range
00106 
00107   static UInt_t  mEtaFtpcCutN;               // number not accepted
00108   static Float_t mEtaFtpcCuts[4];            // range
00109 
00110   static UInt_t  mChgTpcCutN;                // number not accepted
00111   static Int_t   mChgTpcCuts[2];             // range
00112 
00113   static UInt_t  mChgFtpcCutN;               // number not accepted
00114   static Int_t   mChgFtpcCuts[2];            // range
00115 
00116   static UInt_t  mTpcTrackCutN;              // number not accepted
00117   static Bool_t  mTpcTrackCut;               // cut or don't cut (kTRUE, kFALSE)
00118 
00119   static UInt_t  mFtpcTrackCutN;             // number not accepted
00120   static Bool_t  mFtpcTrackCut;              // cut or don't cut (kTRUE, kFALSE)
00121 
00122   ClassDef(StFlowCutTrack,1)                 // macro for rootcint
00123 }; 
00124 
00125 inline UInt_t StFlowCutTrack::EtaSymPosTpc() { return mEtaSymPosTpcN; }
00126 
00127 inline UInt_t StFlowCutTrack::EtaSymNegTpc() { return mEtaSymNegTpcN; }
00128 
00129 inline UInt_t StFlowCutTrack::EtaSymPosFtpc() { return mEtaSymPosFtpcN; }
00130 
00131 inline UInt_t StFlowCutTrack::EtaSymNegFtpc() { return mEtaSymNegFtpcN; }
00132 
00133 inline void StFlowCutTrack::EtaSymClear() { mEtaSymPosTpcN = 0; mEtaSymNegTpcN = 0; 
00134                                             mEtaSymPosFtpcN = 0; mEtaSymNegFtpcN = 0; }
00135 
00136 inline void StFlowCutTrack::SetFitPtsTpc(Int_t lo, Int_t hi) {
00137   mFitPtsTpcCuts[0] = lo; mFitPtsTpcCuts[1] = hi; }
00138 
00139 inline void StFlowCutTrack::SetFitPtsFtpc(Int_t lo, Int_t hi) {
00140   mFitPtsFtpcCuts[0] = lo; mFitPtsFtpcCuts[1] = hi; }
00141 
00142 inline void StFlowCutTrack::SetFitOverMaxPts(Float_t lo, Float_t hi) {
00143   mFitOverMaxCuts[0] = lo; mFitOverMaxCuts[1] = hi; }
00144 
00145 inline void StFlowCutTrack::SetChiSqTpc(Float_t lo, Float_t hi) {
00146   mChiSqTpcCuts[0] = lo; mChiSqTpcCuts[1] = hi; }
00147 
00148 inline void StFlowCutTrack::SetChiSqFtpc(Float_t lo, Float_t hi) {
00149   mChiSqFtpcCuts[0] = lo; mChiSqFtpcCuts[1] = hi; }
00150 
00151 inline void StFlowCutTrack::SetDcaFtpc(Float_t lo, Float_t hi) {
00152   mDcaFtpcCuts[0] = lo; mDcaFtpcCuts[1] = hi; }
00153 
00154 inline void StFlowCutTrack::SetDcaGlobalFtpc(Float_t lo, Float_t hi) {
00155   mDcaGlobalFtpcCuts[0] = lo; mDcaGlobalFtpcCuts[1] = hi; }
00156 
00157 inline void StFlowCutTrack::SetPtTpc(Float_t lo, Float_t hi) {
00158   mPtTpcCuts[0] = lo; mPtTpcCuts[1] = hi; }
00159 
00160 inline void StFlowCutTrack::SetPtFtpc(Float_t lo, Float_t hi) {
00161   mPtFtpcCuts[0] = lo; mPtFtpcCuts[1] = hi; }
00162 
00163 inline void StFlowCutTrack::SetEtaTpc(Float_t lo, Float_t hi) {
00164   mEtaTpcCuts[0] = lo; mEtaTpcCuts[1] = hi; }
00165 
00166 inline void StFlowCutTrack::SetEtaFtpc(Float_t lo_neg, Float_t hi_neg,
00167                                        Float_t lo_pos, Float_t hi_pos) {
00168   mEtaFtpcCuts[0] = lo_neg; mEtaFtpcCuts[1] = hi_neg; 
00169   mEtaFtpcCuts[2] = lo_pos; mEtaFtpcCuts[3] = hi_pos; }
00170 
00171 inline void StFlowCutTrack::SetChgTpc(Int_t lo, Int_t hi) {
00172   mChgTpcCuts[0] = lo; mChgTpcCuts[1] = hi; }
00173 
00174 inline void StFlowCutTrack::SetChgFtpc(Int_t lo, Int_t hi) {
00175   mChgFtpcCuts[0] = lo; mChgFtpcCuts[1] = hi; }
00176 
00177 inline void StFlowCutTrack::IncludeTpcTracks(Bool_t tpc_include) {
00178     mTpcTrackCut = !tpc_include; }
00179 
00180 inline void StFlowCutTrack::IncludeFtpcTracks(Bool_t ftpc_include) {
00181     mFtpcTrackCut = !ftpc_include; }
00182 
00183 #endif
00184 
00186 //
00187 // $Log: StFlowCutTrack.h,v $
00188 // Revision 1.20  2003/09/02 17:58:11  perev
00189 // gcc 3.2 updates + WarnOff
00190 //
00191 // Revision 1.19  2003/01/14 14:14:16  oldi
00192 // Possibility to exclude TPC tracks completely (= FTPC only).
00193 //
00194 // Revision 1.18  2003/01/10 16:42:05  oldi
00195 // Several changes to comply with FTPC tracks:
00196 // - Switch to include/exclude FTPC tracks introduced.
00197 //   The same switch changes the range of the eta histograms.
00198 // - Eta symmetry plots for FTPC tracks added and separated from TPC plots.
00199 // - PhiWgts and related histograms for FTPC tracks split in FarEast, East,
00200 //   West, FarWest (depending on vertex.z()).
00201 // - Psi_Diff plots for 2 different selections and the first 2 harmonics added.
00202 // - Cut to exclude mu-events with no primary vertex introduced.
00203 //   (This is possible for UPC events and FTPC tracks.)
00204 // - Global DCA cut for FTPC tracks added.
00205 // - Global DCA cuts for event plane selection separated for TPC and FTPC tracks.
00206 // - Charge cut for FTPC tracks added.
00207 //
00208 // Revision 1.17  2002/06/10 22:50:59  posk
00209 // pt and eta weighting now default.
00210 // DcaGlobalPart default now 0 to 1 cm.
00211 // Event cut order changed.
00212 //
00213 // Revision 1.16  2002/06/07 22:18:40  kirill
00214 // Introduced MuDst reader
00215 //
00216 // Revision 1.15  2002/02/13 22:29:17  posk
00217 // Pt Weight now also weights Phi Weights. Added Eta Weight, default=FALSE.
00218 //
00219 // Revision 1.14  2002/01/31 21:43:14  aihong
00220 // add SetChgTpc()
00221 //
00222 // Revision 1.13  2001/11/09 21:10:34  posk
00223 // Switched from CERNLIB to TMath. Little q is now normalized.
00224 //
00225 // Revision 1.12  2001/07/27 01:26:10  snelling
00226 // Added and changed variables for picoEvent. Changed trackCut class to StTrack
00227 //
00228 // Revision 1.11  2001/05/22 20:17:23  posk
00229 // Now can do pseudorapidity subevents.
00230 //
00231 // Revision 1.10  2000/12/12 20:22:05  posk
00232 // Put log comments at end of files.
00233 // Deleted persistent StFlowEvent (old micro DST).
00234 //
00235 // Revision 1.9  2000/12/10 02:01:13  oldi
00236 // A new member (StTrackTopologyMap mTopology) was added to StFlowPicoTrack.
00237 // The evaluation of either a track originates from the FTPC or not is
00238 // unambiguous now. The evaluation itself is easily extendible for other
00239 // detectors (e.g. SVT+TPC). Old flowpicoevent.root files are treated as if
00240 // they contain TPC tracks only (backward compatibility).
00241 //
00242 // Revision 1.7  2000/12/08 17:03:38  oldi
00243 // Phi weights for both FTPCs included.
00244 //
00245 // Revision 1.6  2000/12/06 15:38:46  oldi
00246 // Including FTPC.
00247 //
00248 // Revision 1.5  2000/10/12 22:46:35  snelling
00249 // Added support for the new pDST's and the probability pid method
00250 //
00251 // Revision 1.4  2000/08/31 18:58:20  posk
00252 // For picoDST, added version number, runID, and multEta for centrality.
00253 // Added centrality cut when reading picoDST.
00254 // Added pt and eta selections for particles corr. wrt event plane.
00255 //
00256 // Revision 1.3  2000/08/10 23:00:21  posk
00257 // New centralities. pt and eta cuts.
00258 //
00259 // Revision 1.2  2000/07/12 17:54:36  posk
00260 // Added chi2 and dca cuts. Multiplied EtaSym by ::sqrt(mult).
00261 // Apply cuts when reading picoevent file.
00262 //
00263 // Revision 1.1  2000/03/02 23:02:43  posk
00264 // Changed extensions from .hh and .cc to .h and .cxx .
00265 //
00266 // Revision 1.4  1999/12/15 22:01:24  posk
00267 // Added StFlowConstants.hh
00268 //
00269 // Revision 1.3  1999/11/30 18:52:50  snelling
00270 // First modification for the new StEvent
00271 //
00272 // Revision 1.2  1999/11/24 18:17:11  posk
00273 // Put the methods which act on the data in with the data in StFlowEvent.
00274 //
00275 // Revision 1.1  1999/11/05 00:06:45  posk
00276 // First versions of Flow cut classes.
00277 //