StRoot: StFlowMaker/StFlowEvent.h Source File

00001 
00002 //
00003 // $Id: StFlowEvent.h,v 1.56 2009/11/24 19:23:03 posk Exp $
00004 //
00005 // Author: Raimond Snellings and Art Poskanzer
00006 //          FTPC added by Markus Oldenburg, MPI, Dec 2000
00007 //          Cumulants added by Aihong Tang, KSU, Nov 2001
00008 //
00010 //
00011 // Description: A subset of StEvent with flow functions
00012 //
00014 
00015 #ifndef StFlowEvent_h
00016 #define StFlowEvent_h
00017 #include "StObject.h"
00018 #include "StFlowTrackCollection.h"
00019 #include "StTrackTopologyMap.h"
00020 #include "StThreeVectorF.hh"
00021 #include "StEnumerations.h"
00022 #include "Rtypes.h"
00023 #include "TVector2.h"
00024 #include "TComplex.h"
00025 class StFlowSelection;
00026 
00027 class StFlowEvent  : public StObject {
00028 
00029 public:
00030 
00031   StFlowEvent();
00032   virtual        ~StFlowEvent();
00033 
00034   Double_t       PhiWeight(Int_t selN, Int_t harN, StFlowTrack* pFlowTrack) const;
00035   Double_t       PhiWeightRaw(Int_t selN, Int_t harN, StFlowTrack* pFlowTrack) const;
00036   Double_t       Weight(Int_t selN, Int_t harN, StFlowTrack* pFlowTrack) const;
00037   Double_t       ZDCSMD_PsiWgtEast();
00038   Double_t       ZDCSMD_PsiWgtWest();
00039   Double_t       ZDCSMD_PsiWgtFull();
00040   Int_t          EventID() const;
00041   Int_t          RunID() const;
00042   Double_t       CenterOfMassEnergy() const;
00043   Double_t       MagneticField() const;
00044   Short_t        BeamMassNumberEast() const;
00045   Short_t        BeamMassNumberWest() const;
00046   UInt_t         OrigMult() const;
00047   UInt_t         L0TriggerWord() const;
00048   UInt_t         UncorrNegMult() const;
00049   UInt_t         UncorrPosMult() const;
00050   UInt_t         MultEta() const;
00051   UInt_t         FlowEventMult() const;
00052   UInt_t         Centrality() const;
00053   StThreeVectorF VertexPos() const;
00054   UInt_t         Mult(StFlowSelection*);
00055   UInt_t         MultPart(StFlowSelection*);
00056   TVector2       Q(StFlowSelection*);
00057   TVector2       NormQ(StFlowSelection* pFlowSelect);
00058   TVector2       QPart(StFlowSelection*);
00059   TVector2       ReCentEPPar(StFlowSelection*, char*); // for ana
00060   TVector2       ReCentPar(StFlowSelection*, char*); // for LYZ
00061   TVector2       ReCent(Int_t selN, Int_t harN, StFlowTrack* pFlowTrack) const;
00062   TVector2       ReCentEP(Int_t selN, Int_t harN, StFlowTrack* pFlowTrack) const;
00063   Float_t        q(StFlowSelection*);
00064   Float_t        MeanPt(StFlowSelection*);
00065   Float_t        Qtheta(StFlowSelection*, Float_t theta);
00066   TComplex       Grtheta(StFlowSelection*, Float_t r, Float_t theta);
00067   TComplex       GV1r0theta(StFlowSelection*, Float_t r, Float_t theta1, Float_t theta);
00068   TComplex       Gder_r0theta(StFlowSelection*, Float_t r, Float_t theta);
00069   Float_t        Psi(StFlowSelection*);
00070   Float_t        ZDCSMD_PsiCorr();
00071   Float_t        ZDCSMD_PsiEst();
00072   Float_t        ZDCSMD_PsiWst();
00073   Float_t        ZDCSMD_GetPosition(int eastwest,int verthori,int strip);
00074   Double_t       G_New(StFlowSelection* pFlowSelect, Double_t Zx, Double_t Zy);
00075   Double_t       G_Mix(StFlowSelection* pFlowSelect, Double_t Z1x, Double_t Z1y, Double_t Z2x, Double_t Z2y);
00076   Double_t       SumWeightSquare(StFlowSelection* pFlowSelect);
00077   Double_t       PtAbsWgtValue(Double_t pt) const;
00078   Double_t       EtaAbsWgtValue(Double_t eta) const;
00079   Float_t        CTB() const;
00080   Float_t        ZDCe() const;
00081   Float_t        ZDCw() const;
00082   Float_t        ZDCSMD(int eastwest,int verthori,int strip) const;
00083   Float_t        PtWgtSaturation() const;
00084   Bool_t         PtWgt() const;
00085   Bool_t         EtaWgt() const;
00086   Bool_t         FirstLastPhiWgt() const;
00087   Bool_t         FirstLastPoints() const;
00088   Bool_t         UseZDCSMD() const;
00089   Char_t*        Pid();
00090   Bool_t         EtaSubs() const;
00091   Bool_t         RanSubs() const;
00092 
00093   Float_t V1TPCDetctWgtG_Mix(Int_t selN) const;
00094   Float_t V1FtpcEastDetctWgtG_Mix(Int_t selN) const;
00095   Float_t V1FtpcWestDetctWgtG_Mix(Int_t selN) const;
00096 
00097   Float_t V2TPCDetctWgtG_Mix(Int_t selN) const;
00098   Float_t V2FtpcEastDetctWgtG_Mix(Int_t selN) const;
00099   Float_t V2FtpcWestDetctWgtG_Mix(Int_t selN) const;
00100 
00101 
00102   StFlowTrackCollection* TrackCollection() const;
00103 
00104   void SetSelections();
00105   void SetPid(const Char_t*);
00106   void SetPidsDeviant();
00107   void SetPidsProb();
00108   void PrintSelectionList();
00109   void MakeSubEvents();
00110   void MakeEtaSubEvents();
00111   void SetEventID(const Int_t&);
00112   void SetRunID(const Int_t&);
00113   void SetCenterOfMassEnergy(const Double_t&);
00114   void SetMagneticField(const Double_t&);
00115   void SetBeamMassNumberEast(const Short_t&);
00116   void SetBeamMassNumberWest(const Short_t&);
00117   void SetOrigMult(const UInt_t&);
00118   void SetL0TriggerWord(const UInt_t&);
00119   void SetUncorrPosMult(const UInt_t&);
00120   void SetUncorrNegMult(const UInt_t&);
00121   void SetMultEta(const UInt_t&);
00122   void SetCentrality();
00123   void SetVertexPos(const StThreeVectorF&);
00124   void SetCTB(const Float_t ctb);
00125   void SetZDCe(const Float_t zdce);
00126   void SetZDCw(const Float_t zdcw);
00127   void SetZDCSMD(int eastwest,int verthori,int strip,const Float_t zdcsmd);
00128 #ifndef __CINT__                
00129   void SetPhiWeightFarEast(const Flow::PhiWgt_t &pPhiWgt);
00130   void SetPhiWeightEast(const Flow::PhiWgt_t &pPhiWgt);
00131   void SetPhiWeightWest(const Flow::PhiWgt_t &pPhiWgt);
00132   void SetPhiWeightFarWest(const Flow::PhiWgt_t &pPhiWgt);
00133   void SetPhiWeightFtpcFarEast(const Flow::PhiWgtFtpc_t &pPhiWgt);
00134   void SetPhiWeightFtpcEast(const Flow::PhiWgtFtpc_t &pPhiWgt);
00135   void SetPhiWeightFtpcWest(const Flow::PhiWgtFtpc_t &pPhiWgt);
00136   void SetPhiWeightFtpcFarWest(const Flow::PhiWgtFtpc_t &pPhiWgt);
00137   void SetZDCSMD_PsiWeightWest(const Flow::ZDCSMD_PsiWgt_t&  ZDCSMD_PsiWgtWest);
00138   void SetZDCSMD_PsiWeightEast(const Flow::ZDCSMD_PsiWgt_t&  ZDCSMD_PsiWgtEast); 
00139   void SetZDCSMD_PsiWeightFull(const Flow::ZDCSMD_PsiWgt_t&  ZDCSMD_PsiWgtFull);
00140   void SetZDCSMD_BeamCenter(Double_t ex,Double_t ey,Double_t wx,Double_t wy);
00141   void SetReCentX(const Flow::ReCent_t &pReCentX);
00142   void SetReCentY(const Flow::ReCent_t &pReCentY);
00143 #endif
00144   static void SetEtaTpcCut(Float_t lo, Float_t hi, Int_t harN, Int_t selN);
00145   static void SetPtTpcCut(Float_t lo, Float_t hi, Int_t harN, Int_t selN);
00146   static void SetEtaFtpcCut(Float_t lo_neg, Float_t hi_neg, 
00147                             Float_t lo_pos, Float_t hi_pos, 
00148                             Int_t harN, Int_t selN);
00149   static void SetPtFtpcCut(Float_t lo, Float_t hi, Int_t harN, Int_t selN);
00150   static void SetPiPlusCut(Float_t lo, Float_t hi);
00151   static void SetPiMinusCut(Float_t lo, Float_t hi);
00152   static void SetProtonCut(Float_t lo, Float_t hi);
00153   static void SetAntiProtonCut(Float_t lo, Float_t hi);
00154   static void SetKPlusCut(Float_t lo, Float_t hi);
00155   static void SetKMinusCut(Float_t lo, Float_t hi);
00156   static void SetElectronCut(Float_t lo, Float_t hi);
00157   static void SetPositronCut(Float_t lo, Float_t hi);
00158   static void SetDeuteronCut(Float_t lo, Float_t hi);
00159   static void SetAntiDeuteronCut(Float_t lo, Float_t hi);
00160   static void SetDcaGlobalTpcCut(Float_t lo, Float_t hi);
00161   static void SetDcaGlobalFtpcCut(Float_t lo, Float_t hi);
00162   static void SetProbPid();
00163   static void SetEtaSubs();
00164   static void SetRanSubs();
00165   static void SetPtWgtSaturation(Float_t);
00166   static void SetPtWgt(Bool_t);
00167   static void SetEtaWgt(Bool_t);
00168   static void SetFirstLastPhiWgt();
00169   static void SetFirstLastPoints();
00170   static void SetUseZDCSMD(Bool_t);
00171 
00172   static void SetV1TPCDetctWgtG_Mix(Float_t val,  Int_t selN);
00173   static void SetV1FtpcEastDetctWgtG_Mix(Float_t val, Int_t selN);
00174   static void SetV1FtpcWestDetctWgtG_Mix(Float_t val,  Int_t selN);
00175 
00176 
00177   static void SetV2TPCDetctWgtG_Mix(Float_t val,  Int_t selN);
00178   static void SetV2FtpcEastDetctWgtG_Mix(Float_t val, Int_t selN);
00179   static void SetV2FtpcWestDetctWgtG_Mix(Float_t val,  Int_t selN);
00180   static Bool_t ProbPid();
00181 
00182 
00183 private:
00184 
00185   Int_t               mEventID;                                  // ID of the event
00186   Int_t               mRunID;                                    // ID of the run
00187   Double_t            mMagneticField;                            // 
00188   Double_t            mCenterOfMassEnergy;                       // 
00189   Short_t             mBeamMassNumberEast;                       //
00190   Short_t             mBeamMassNumberWest;                       //
00191   UInt_t              mOrigMult;                                 // number of tracks
00192   UInt_t              mL0TriggerWord;                            // L0 trigger word
00193   UInt_t              mUncorrNegMult;                            // number of h-
00194   UInt_t              mUncorrPosMult;                            // number of h+
00195   UInt_t              mMultEta;                                  // number of tracks
00196                                                 // with pos. flag in 1.5 unit of eta
00197   UInt_t              mCentrality;                               // centrality bin
00198   StThreeVectorF      mVertexPos;                                // primary vertex position
00199   Float_t             mCTB;                                      // CTB value sum
00200   Float_t             mZDCe;                                     // ZDC east
00201   Float_t             mZDCw;                                     // ZDC west
00202   Float_t             mZDCSMD[2][2][8];                          // ZDCSMD
00203   static Double_t     mZDCSMDCenterex,mZDCSMDCenterey;           // ZDCSMD Beam Center
00204   static Double_t     mZDCSMDCenterwx,mZDCSMDCenterwy;           // ZDCSMD Beam Center
00205   static Float_t      mEtaTpcCuts[2][2][Flow::nSels];            // range absolute values
00206   static Float_t      mEtaFtpcCuts[4][2][Flow::nSels];           // range values
00207   static Float_t      mPtTpcCuts[2][2][Flow::nSels];             // range
00208   static Float_t      mPtFtpcCuts[2][2][Flow::nSels];            // range
00209 
00210   static Float_t      mV1TPCDetctWgtG_Mix[Flow::nSels]; // detector wgt for G_Mix in v1{3} calc.
00211   static Float_t      mV1FtpcEastDetctWgtG_Mix[Flow::nSels]; 
00212   static Float_t      mV1FtpcWestDetctWgtG_Mix[Flow::nSels]; 
00213   static Float_t      mV2TPCDetctWgtG_Mix[Flow::nSels]; 
00214   static Float_t      mV2FtpcEastDetctWgtG_Mix[Flow::nSels]; 
00215   static Float_t      mV2FtpcWestDetctWgtG_Mix[Flow::nSels]; 
00216 
00217   Flow::PhiWgt_t      mPhiWgtFarEast;                            
00218   Flow::PhiWgt_t      mPhiWgtEast;                               
00219   Flow::PhiWgt_t      mPhiWgtWest;                               
00220   Flow::PhiWgt_t      mPhiWgtFarWest;                            
00221   Flow::PhiWgtFtpc_t  mPhiWgtFtpcFarEast;                        
00222   Flow::PhiWgtFtpc_t  mPhiWgtFtpcEast;                           
00223   Flow::PhiWgtFtpc_t  mPhiWgtFtpcWest;                           
00224   Flow::PhiWgtFtpc_t  mPhiWgtFtpcFarWest;                        
00225   Flow::ZDCSMD_PsiWgt_t  mZDCSMD_PsiWgtWest;                     
00226   Flow::ZDCSMD_PsiWgt_t  mZDCSMD_PsiWgtEast;                     
00227   Flow::ZDCSMD_PsiWgt_t  mZDCSMD_PsiWgtFull;                     
00228   Flow::ReCent_t      mReCentX;                                  
00229   Flow::ReCent_t      mReCentY;                                  
00230 
00231   static Float_t      mPiPlusCuts[2];                            // PID cuts
00232   static Float_t      mPtWgtSaturation;                          // saturation value for pt weighting
00233   static Bool_t       mPtWgt;                                    // flag for pt weighting
00234   static Bool_t       mEtaWgt;                                   // flag for y weighting for odd harmonics
00235   static Char_t       mPid[10];                                  // h+, h-, pi-, pi+, pi, k+, k-, k, pr+, pr-, pr, e+, e-, e
00236   static Bool_t       mProbPid;                                  // flag for probability pid
00237   static Bool_t       mEtaSubs;                                  // flag for eta subevents
00238   static Bool_t       mRanSubs;                                  // flag for random subevents
00239   static Bool_t       mFirstLastPhiWgt;                          // flag for using z of first and last points for reading phi weights
00240   static Bool_t       mFirstLastPoints;                          // flag for using z of first and last points for generating phi weights
00241   static Bool_t       mUseZDCSMD;                                // flag for using ZDC SMD for RP
00242   static Float_t      mPiMinusCuts[2];
00243   static Float_t      mProtonCuts[2];
00244   static Float_t      mKMinusCuts[2];
00245   static Float_t      mKPlusCuts[2];
00246   static Float_t      mAntiProtonCuts[2];
00247   static Float_t      mDeuteronCuts[2];
00248   static Float_t      mAntiDeuteronCuts[2];
00249   static Float_t      mElectronCuts[2];
00250   static Float_t      mPositronCuts[2];
00251   static Float_t      mDcaGlobalTpcCuts[2];
00252   static Float_t      mDcaGlobalFtpcCuts[2];
00253 
00254   StFlowTrackCollection* pTrackCollection;   
00255 
00256   ClassDef(StFlowEvent,1)                    // macro for rootcint
00257 
00258 };
00259 
00260 inline StFlowTrackCollection* StFlowEvent::TrackCollection() const {
00261   return pTrackCollection; }
00262 
00263 inline Int_t StFlowEvent::EventID() const { return mEventID; }
00264 
00265 inline Int_t StFlowEvent::RunID() const { return mRunID; }
00266 
00267 inline Double_t StFlowEvent::CenterOfMassEnergy() const { return mCenterOfMassEnergy; }
00268 
00269 inline Double_t StFlowEvent::MagneticField() const { return mMagneticField; }
00270 
00271 inline Short_t StFlowEvent::BeamMassNumberEast() const { return mBeamMassNumberEast; }
00272 
00273 inline Short_t StFlowEvent::BeamMassNumberWest() const { return mBeamMassNumberWest; }
00274 
00275 inline UInt_t StFlowEvent::OrigMult() const { return mOrigMult; }
00276 
00277 inline UInt_t StFlowEvent::L0TriggerWord() const { return mL0TriggerWord; }
00278 
00279 inline UInt_t StFlowEvent::UncorrNegMult() const { return mUncorrNegMult; }
00280 
00281 inline UInt_t StFlowEvent::UncorrPosMult() const { return mUncorrPosMult; }
00282 
00283 inline UInt_t StFlowEvent::MultEta() const { return mMultEta; }
00284 
00285 inline UInt_t StFlowEvent::FlowEventMult() const { return pTrackCollection->size(); }
00286 
00287 inline UInt_t StFlowEvent::Centrality() const { return mCentrality; }
00288 
00289 inline StThreeVectorF StFlowEvent::VertexPos() const { return mVertexPos; }
00290 
00291 inline Float_t  StFlowEvent::CTB() const { return mCTB; }
00292 
00293 inline Float_t  StFlowEvent::ZDCe() const { return mZDCe; }
00294 
00295 inline Float_t  StFlowEvent::ZDCw() const { return mZDCw; }  
00296 
00297 inline Float_t  StFlowEvent::ZDCSMD(int eastwest,int verthori,int strip) const {return mZDCSMD[eastwest][verthori][strip-1];}
00298 
00299 inline Float_t  StFlowEvent::PtWgtSaturation() const { return mPtWgtSaturation; }
00300 
00301 inline Bool_t   StFlowEvent::PtWgt() const { return mPtWgt; }
00302 
00303 inline Bool_t   StFlowEvent::EtaWgt() const { return mEtaWgt; }
00304 
00305 inline Bool_t   StFlowEvent::FirstLastPhiWgt() const { return mFirstLastPhiWgt; }
00306 
00307 inline Bool_t   StFlowEvent::FirstLastPoints() const { return mFirstLastPoints; }
00308 
00309 inline Char_t*  StFlowEvent::Pid() { return mPid; }
00310 
00311 inline Bool_t   StFlowEvent::ProbPid() { return mProbPid; }
00312 
00313 inline Bool_t   StFlowEvent::EtaSubs() const { return mEtaSubs; }
00314 
00315 inline Bool_t   StFlowEvent::RanSubs() const { return mRanSubs; }
00316 
00317 inline Bool_t   StFlowEvent::UseZDCSMD() const { return mUseZDCSMD;}
00318 
00319 inline   Float_t StFlowEvent::V1TPCDetctWgtG_Mix(Int_t selN) const {
00320   return mV1TPCDetctWgtG_Mix[selN]; }
00321 inline   Float_t StFlowEvent::V1FtpcEastDetctWgtG_Mix(Int_t selN) const {
00322   return mV1FtpcEastDetctWgtG_Mix[selN]; }
00323 inline   Float_t StFlowEvent::V1FtpcWestDetctWgtG_Mix(Int_t selN) const {
00324   return mV1FtpcWestDetctWgtG_Mix[selN]; }
00325 
00326 inline   Float_t StFlowEvent::V2TPCDetctWgtG_Mix(Int_t selN) const {
00327   return mV2TPCDetctWgtG_Mix[selN]; }
00328 inline   Float_t StFlowEvent::V2FtpcEastDetctWgtG_Mix(Int_t selN) const {
00329   return mV2FtpcEastDetctWgtG_Mix[selN]; }
00330 inline   Float_t StFlowEvent::V2FtpcWestDetctWgtG_Mix(Int_t selN) const {
00331   return mV2FtpcWestDetctWgtG_Mix[selN]; }
00332 
00333 #ifndef __CINT__
00334 inline void StFlowEvent::SetPhiWeightFarEast(const Flow::PhiWgt_t& pPhiWgtFarEast) {
00335   memcpy (mPhiWgtFarEast, pPhiWgtFarEast, sizeof(Flow::PhiWgt_t)); }
00336 
00337 inline void StFlowEvent::SetPhiWeightEast(const Flow::PhiWgt_t& pPhiWgtEast) {
00338   memcpy (mPhiWgtEast, pPhiWgtEast, sizeof(Flow::PhiWgt_t)); }
00339 
00340 inline void StFlowEvent::SetPhiWeightWest(const Flow::PhiWgt_t& pPhiWgtWest) {
00341   memcpy (mPhiWgtWest, pPhiWgtWest, sizeof(Flow::PhiWgt_t)); }
00342 
00343 inline void StFlowEvent::SetPhiWeightFarWest(const Flow::PhiWgt_t& pPhiWgtFarWest) {
00344   memcpy (mPhiWgtFarWest, pPhiWgtFarWest, sizeof(Flow::PhiWgt_t)); }
00345 
00346 inline void StFlowEvent::SetPhiWeightFtpcFarEast(const Flow::PhiWgtFtpc_t& pPhiWgtFtpcFarEast) {
00347   memcpy (mPhiWgtFtpcFarEast, pPhiWgtFtpcFarEast, sizeof(Flow::PhiWgtFtpc_t)); }
00348 
00349 inline void StFlowEvent::SetPhiWeightFtpcEast(const Flow::PhiWgtFtpc_t& pPhiWgtFtpcEast) {
00350   memcpy (mPhiWgtFtpcEast, pPhiWgtFtpcEast, sizeof(Flow::PhiWgtFtpc_t)); }
00351 
00352 inline void StFlowEvent::SetPhiWeightFtpcWest(const Flow::PhiWgtFtpc_t& pPhiWgtFtpcWest) {
00353   memcpy (mPhiWgtFtpcWest, pPhiWgtFtpcWest, sizeof(Flow::PhiWgtFtpc_t)); }
00354 
00355 inline void StFlowEvent::SetPhiWeightFtpcFarWest(const Flow::PhiWgtFtpc_t& pPhiWgtFtpcFarWest) {
00356   memcpy (mPhiWgtFtpcFarWest, pPhiWgtFtpcFarWest, sizeof(Flow::PhiWgtFtpc_t)); }
00357 
00358 inline void StFlowEvent::SetZDCSMD_PsiWeightEast(const Flow::ZDCSMD_PsiWgt_t&  ZDCSMD_PsiWgtEast) {
00359   memcpy (mZDCSMD_PsiWgtEast, ZDCSMD_PsiWgtEast, sizeof(Flow::ZDCSMD_PsiWgt_t)); }
00360 
00361 inline void StFlowEvent::SetZDCSMD_PsiWeightWest(const Flow::ZDCSMD_PsiWgt_t&  ZDCSMD_PsiWgtWest) {
00362   memcpy (mZDCSMD_PsiWgtWest, ZDCSMD_PsiWgtWest, sizeof(Flow::ZDCSMD_PsiWgt_t)); }
00363 
00364 inline void StFlowEvent::SetZDCSMD_PsiWeightFull(const Flow::ZDCSMD_PsiWgt_t&  ZDCSMD_PsiWgtFull) {
00365   memcpy (mZDCSMD_PsiWgtFull, ZDCSMD_PsiWgtFull, sizeof(Flow::ZDCSMD_PsiWgt_t)); }
00366  
00367 inline void StFlowEvent::SetZDCSMD_BeamCenter(Double_t ex,Double_t ey,Double_t wx,Double_t wy) {
00368   mZDCSMDCenterex = ex; mZDCSMDCenterey = ey; mZDCSMDCenterwx = wx; mZDCSMDCenterwy = wy;}
00369 
00370 inline void StFlowEvent::SetReCentX(const Flow::ReCent_t& pReCentX) {
00371   memcpy (mReCentX, pReCentX, sizeof(Flow::ReCent_t)); }
00372 
00373 inline void StFlowEvent::SetReCentY(const Flow::ReCent_t& pReCentY) {
00374   memcpy (mReCentY, pReCentY, sizeof(Flow::ReCent_t)); }
00375 #endif
00376 
00377 inline void StFlowEvent::SetEtaTpcCut(Float_t lo, Float_t hi, Int_t harN, Int_t selN)
00378 { mEtaTpcCuts[0][harN][selN] = lo; mEtaTpcCuts[1][harN][selN] = hi; }
00379 
00380 inline void StFlowEvent::SetPtTpcCut(Float_t lo, Float_t hi, Int_t harN, Int_t selN)
00381 { mPtTpcCuts[0][harN][selN] = lo; mPtTpcCuts[1][harN][selN] = hi; }
00382 
00383 inline void StFlowEvent::SetEtaFtpcCut(Float_t lo_neg, Float_t hi_neg, 
00384                                        Float_t lo_pos, Float_t hi_pos, 
00385                                        Int_t harN, Int_t selN)
00386 { mEtaFtpcCuts[0][harN][selN] = lo_neg; mEtaFtpcCuts[1][harN][selN] = hi_neg; 
00387  mEtaFtpcCuts[2][harN][selN] = lo_pos; mEtaFtpcCuts[3][harN][selN] = hi_pos; }
00388 
00389 inline void StFlowEvent::SetPtFtpcCut(Float_t lo, Float_t hi, Int_t harN, Int_t selN)
00390 { mPtFtpcCuts[0][harN][selN] = lo; mPtFtpcCuts[1][harN][selN] = hi; }
00391 
00392 inline void StFlowEvent::SetEventID(const Int_t& id) { mEventID = id; }
00393 
00394 inline void StFlowEvent::SetRunID(const Int_t& id) { mRunID = id; }
00395 
00396 inline void StFlowEvent::SetMagneticField(const Double_t& mf) { mMagneticField = mf; }
00397 
00398 inline void StFlowEvent::SetCenterOfMassEnergy(const Double_t& cms) { mCenterOfMassEnergy = cms; }
00399 
00400 inline void StFlowEvent::SetBeamMassNumberEast(const Short_t& bme) { mBeamMassNumberEast = bme; }
00401 
00402 inline void StFlowEvent::SetBeamMassNumberWest(const Short_t& bmw) { mBeamMassNumberWest = bmw; }
00403 
00404 inline void StFlowEvent::SetOrigMult(const UInt_t& tracks) {
00405   mOrigMult = tracks; }
00406 
00407 inline void StFlowEvent::SetL0TriggerWord(const UInt_t& trigger) {
00408   mL0TriggerWord = trigger; }
00409 
00410 inline void StFlowEvent::SetUncorrNegMult(const UInt_t& negtracks) {
00411   mUncorrNegMult = negtracks; }
00412 
00413 inline void StFlowEvent::SetUncorrPosMult(const UInt_t& postracks) {
00414   mUncorrPosMult = postracks; }
00415 
00416 inline void StFlowEvent::SetMultEta(const UInt_t& goodtracks) {
00417   mMultEta = goodtracks; }
00418 
00419 inline void StFlowEvent::SetVertexPos(const StThreeVectorF& vertexPos) {
00420   mVertexPos = vertexPos; }
00421 
00422 inline void StFlowEvent::SetPiPlusCut(Float_t lo, Float_t hi) { 
00423   mPiPlusCuts[0] = lo; mPiPlusCuts[1] = hi; }
00424 
00425 inline void StFlowEvent::SetPiMinusCut(Float_t lo, Float_t hi) { 
00426   mPiMinusCuts[0] = lo; mPiMinusCuts[1] = hi; }
00427 
00428 inline void StFlowEvent::SetProtonCut(Float_t lo, Float_t hi) { 
00429   mProtonCuts[0] = lo; mProtonCuts[1] = hi; }
00430 
00431 inline void StFlowEvent::SetKMinusCut(Float_t lo, Float_t hi) { 
00432   mKMinusCuts[0] = lo; mKMinusCuts[1] = hi; }
00433 
00434 inline void StFlowEvent::SetKPlusCut(Float_t lo, Float_t hi) { 
00435   mKPlusCuts[0] = lo; mKPlusCuts[1] = hi; }
00436 
00437 inline void StFlowEvent::SetAntiProtonCut(Float_t lo, Float_t hi) { 
00438   mAntiProtonCuts[0] = lo; mAntiProtonCuts[1] = hi; }
00439 
00440 inline void StFlowEvent::SetDeuteronCut(Float_t lo, Float_t hi) { 
00441   mDeuteronCuts[0] = lo; mDeuteronCuts[1] = hi; }
00442 
00443 inline void StFlowEvent::SetAntiDeuteronCut(Float_t lo, Float_t hi) { 
00444   mAntiDeuteronCuts[0] = lo; mAntiDeuteronCuts[1] = hi; }
00445 
00446 inline void StFlowEvent::SetElectronCut(Float_t lo, Float_t hi) { 
00447   mElectronCuts[0] = lo; mElectronCuts[1] = hi; }
00448 
00449 inline void StFlowEvent::SetPositronCut(Float_t lo, Float_t hi) { 
00450   mPositronCuts[0] = lo; mPositronCuts[1] = hi; }
00451 
00452 inline void StFlowEvent::SetDcaGlobalTpcCut(Float_t lo, Float_t hi) { 
00453   mDcaGlobalTpcCuts[0] = lo; mDcaGlobalTpcCuts[1] = hi; }
00454 
00455 inline void StFlowEvent::SetDcaGlobalFtpcCut(Float_t lo, Float_t hi) { 
00456   mDcaGlobalFtpcCuts[0] = lo; mDcaGlobalFtpcCuts[1] = hi; }
00457 
00458 inline void StFlowEvent::SetCTB(const Float_t ctb) { mCTB = ctb; }
00459 
00460 inline void StFlowEvent::SetZDCe(const Float_t zdce) { mZDCe = zdce; }
00461 
00462 inline void StFlowEvent::SetZDCw(const Float_t zdcw) { mZDCw = zdcw; }
00463 
00464 inline void StFlowEvent::SetZDCSMD(int eastwest,int verthori,int strip,const Float_t zdcsmd) {
00465   mZDCSMD[eastwest][verthori][strip-1] = (zdcsmd >0.)? zdcsmd:0.;}
00466 
00467 inline void StFlowEvent::SetPid(const Char_t* pid)  { 
00468   strncpy(mPid, pid, 9); mPid[9] = '\0'; }
00469 
00470 inline void StFlowEvent::SetProbPid() { mProbPid = kTRUE; }
00471 
00472 inline void StFlowEvent::SetEtaSubs() { mEtaSubs = kTRUE; }
00473 
00474 inline void StFlowEvent::SetRanSubs() { mRanSubs = kTRUE; }
00475 
00476 inline void StFlowEvent::SetFirstLastPhiWgt() { mFirstLastPhiWgt = kTRUE; }
00477 
00478 inline void StFlowEvent::SetFirstLastPoints() { mFirstLastPoints = kTRUE; }
00479 
00480 inline void StFlowEvent::SetUseZDCSMD(Bool_t UseZDCSMD) { mUseZDCSMD = UseZDCSMD; }
00481 
00482 inline void StFlowEvent::SetPtWgtSaturation(Float_t PtWgtSaturation) { mPtWgtSaturation = PtWgtSaturation; }
00483 
00484 inline void StFlowEvent::SetPtWgt(Bool_t PtWgt) { mPtWgt = PtWgt; }
00485 
00486 inline void StFlowEvent::SetEtaWgt(Bool_t EtaWgt) { mEtaWgt = EtaWgt; }
00487 
00488 inline void StFlowEvent::SetV1TPCDetctWgtG_Mix(Float_t val,  Int_t selN){
00489   mV1TPCDetctWgtG_Mix[selN]=val;}
00490 inline void StFlowEvent::SetV1FtpcEastDetctWgtG_Mix(Float_t val,  Int_t selN){
00491   mV1FtpcEastDetctWgtG_Mix[selN]=val;}
00492 inline void StFlowEvent::SetV1FtpcWestDetctWgtG_Mix(Float_t val,  Int_t selN){
00493   mV1FtpcWestDetctWgtG_Mix[selN]=val;}
00494 
00495 inline void StFlowEvent::SetV2TPCDetctWgtG_Mix(Float_t val,  Int_t selN){
00496   mV2TPCDetctWgtG_Mix[selN]=val;}
00497 inline void StFlowEvent::SetV2FtpcEastDetctWgtG_Mix(Float_t val,  Int_t selN){
00498   mV2FtpcEastDetctWgtG_Mix[selN]=val;}
00499 inline void StFlowEvent::SetV2FtpcWestDetctWgtG_Mix(Float_t val,  Int_t selN){
00500   mV2FtpcWestDetctWgtG_Mix[selN]=val;}
00501 
00502 #endif
00503 
00505 //
00506 // $Log: StFlowEvent.h,v $
00507 // Revision 1.56  2009/11/24 19:23:03  posk
00508 // Added reCenter option to remove acceptance correlations instead of phiWgt.
00509 //
00510 // Revision 1.55  2007/02/06 18:57:54  posk
00511 // In Lee Yang Zeros method, introduced recentering of Q vector.
00512 // Reactivated eta symmetry cut.
00513 //
00514 // Revision 1.54  2006/07/06 16:56:01  posk
00515 // Calculation of v1 for selection=2 is done with mixed harmonics.
00516 //
00517 // Revision 1.53  2006/02/22 19:29:16  posk
00518 // Additions needed for the StFlowLeeYangZerosMaker
00519 //
00520 // Revision 1.52  2005/02/10 21:04:57  aihong
00521 // test mProbPid of StFlowEvent before launch calculation pid on fly
00522 //
00523 // Revision 1.51  2004/12/17 22:33:08  aihong
00524 // add in full Psi weight for ZDC SMD and fix a few bugs, done by Gang
00525 //
00526 // Revision 1.50  2004/12/17 15:50:08  aihong
00527 // check in v1{3} code
00528 //
00529 // Revision 1.49  2004/12/07 17:04:46  posk
00530 // Eliminated the very old mOnePhiWgt, which used one phiWgt histogram for flttening
00531 // instead of four.
00532 //
00533 // Revision 1.48  2004/11/16 21:22:22  aihong
00534 // removed old cumulant method
00535 //
00536 // Revision 1.47  2004/05/05 21:13:45  aihong
00537 // Gang's code for ZDC-SMD added
00538 //
00539 // Revision 1.46  2004/03/11 17:58:42  posk
00540 // Added Random Subs analysis method.
00541 //
00542 // Revision 1.45  2003/07/30 22:00:40  oldi
00543 // Eta cuts for event plane selection separated for FTPC east and west.
00544 // PtWgtSaturation parameter introduced (default set to 2. -> no change of default behavior).
00545 //
00546 // Revision 1.44  2003/06/18 17:00:59  posk
00547 // Event plane cuts now only odd and even, instead of different for each harmonic.
00548 //
00549 // Revision 1.43  2003/04/01 00:27:07  posk
00550 // Little q is now unweighted by pt or eta. Big Q is unaffected.
00551 //
00552 // Revision 1.42  2003/01/10 16:42:15  oldi
00553 // Several changes to comply with FTPC tracks:
00554 // - Switch to include/exclude FTPC tracks introduced.
00555 //   The same switch changes the range of the eta histograms.
00556 // - Eta symmetry plots for FTPC tracks added and separated from TPC plots.
00557 // - PhiWgts and related histograms for FTPC tracks split in FarEast, East,
00558 //   West, FarWest (depending on vertex.z()).
00559 // - Psi_Diff plots for 2 different selections and the first 2 harmonics added.
00560 // - Cut to exclude mu-events with no primary vertex introduced.
00561 //   (This is possible for UPC events and FTPC tracks.)
00562 // - Global DCA cut for FTPC tracks added.
00563 // - Global DCA cuts for event plane selection separated for TPC and FTPC tracks.
00564 // - Charge cut for FTPC tracks added.
00565 //
00566 // Revision 1.41  2003/01/08 19:26:48  posk
00567 // PhiWgt hists sorted on sign of z of first and last points.
00568 // Version 6 of pico file.
00569 //
00570 // Revision 1.40  2002/05/23 18:54:11  posk
00571 // Moved centrality cuts into StFlowConstants
00572 //
00573 // Revision 1.39  2002/03/15 16:43:22  snelling
00574 // Added a method to recalculate the centrality in StFlowPicoEvent
00575 //
00576 // Revision 1.38  2002/03/14 18:51:50  snelling
00577 // Added new centralities
00578 //
00579 // Revision 1.37  2002/02/13 22:29:28  posk
00580 // Pt Weight now also weights Phi Weights. Added Eta Weight, default=FALSE.
00581 //
00582 // Revision 1.36  2002/02/01 23:06:26  snelling
00583 // Added entries for header information in flowPico (not everthing is available yet)
00584 //
00585 // Revision 1.35  2002/01/31 01:04:47  posk
00586 // *** empty log message ***
00587 //
00588 // Revision 1.34  2001/12/18 19:22:09  posk
00589 // "proton" and "antiproton" changed to "pr+" and "pr-".
00590 // Compiles on Solaris.
00591 //
00592 // Revision 1.33  2001/12/11 21:33:50  posk
00593 // Went from one to four sets of histograms for making the event plane isotropic.
00594 // StFlowEvent::PhiWeight() has changed arguments and return value.
00595 // The ptWgt saturates above 2 GeV/c.
00596 //
00597 // Revision 1.32  2001/11/09 21:10:42  posk
00598 // Switched from CERNLIB to TMath. Little q is now normalized.
00599 //
00600 // Revision 1.31  2001/11/02 04:49:52  aihong
00601 // add func. for cumulant maker
00602 //
00603 // Revision 1.30  2001/08/01 19:39:35  snelling
00604 // Added the trigger word
00605 //
00606 // Revision 1.29  2001/07/27 01:26:14  snelling
00607 // Added and changed variables for picoEvent. Changed trackCut class to StTrack
00608 //
00609 // Revision 1.28  2001/07/24 22:29:12  snelling
00610 // First attempt to get a standard root pico file again, added variables
00611 //
00612 // Revision 1.27  2001/06/07 20:06:20  posk
00613 // Global Dca cut for event plane particles.
00614 // Removed SetPtWgt().
00615 //
00616 // Revision 1.26  2001/06/06 13:02:58  rcwells
00617 // Added SetPtWgt(Bool_t) function to StFlowEvent
00618 //
00619 // Revision 1.25  2001/06/04 18:57:05  rcwells
00620 // Adding filling from HbtEvents
00621 //
00622 // Revision 1.24  2001/05/23 18:11:09  posk
00623 // Removed SetPids().
00624 //
00625 // Revision 1.23  2001/05/22 20:17:30  posk
00626 // Now can do pseudorapidity subevents.
00627 //
00628 // Revision 1.22  2001/04/03 17:47:23  oldi
00629 // Bug fix that excluded FTPC tracks from the determination of the reaction plane.
00630 //
00631 // Revision 1.21  2000/12/12 20:22:05  posk
00632 // Put log comments at end of files.
00633 // Deleted persistent StFlowEvent (old micro DST).
00634 //
00635 // Revision 1.20  2000/12/10 02:01:13  oldi
00636 // A new member (StTrackTopologyMap mTopology) was added to StFlowPicoTrack.
00637 // The evaluation of either a track originates from the FTPC or not is
00638 // unambiguous now. The evaluation itself is easily extendible for other
00639 // detectors (e.g. SVT+TPC). Old flowpicoevent.root files are treated as if
00640 // they contain TPC tracks only (backward compatibility).
00641 //
00642 // Revision 1.19  2000/12/08 17:03:38  oldi
00643 // Phi weights for both FTPCs included.
00644 //
00645 // Revision 1.18  2000/10/12 22:46:36  snelling
00646 // Added support for the new pDST's and the probability pid method
00647 //
00648 // Revision 1.16  2000/09/15 22:51:30  posk
00649 // Added pt weighting for event plane calcualtion.
00650 //
00651 // Revision 1.15  2000/09/05 16:11:32  snelling
00652 // Added global DCA, electron and positron
00653 //
00654 // Revision 1.14  2000/08/31 18:58:22  posk
00655 // For picoDST, added version number, runID, and multEta for centrality.
00656 // Added centrality cut when reading picoDST.
00657 // Added pt and eta selections for particles corr. wrt event plane.
00658 //
00659 // Revision 1.13  2000/08/09 21:38:23  snelling
00660 // PID added
00661 //
00662 // Revision 1.12  2000/08/05 22:07:18  fisyak
00663 // less restrictive selection for ROOTCINT
00664 //
00665 // Revision 1.11  2000/08/05 21:21:33  fisyak
00666 // hide from CINT inline functions
00667 //
00668 // Revision 1.10  2000/08/04 21:03:45  perev
00669 // Leaks + Clear() cleanup
00670 //
00671 // Revision 1.9  2000/06/30 14:48:32  posk
00672 // Using MessageMgr, changed Eta Symmetry cut.
00673 //
00674 // Revision 1.8  2000/06/20 16:34:25  snelling
00675 // fixed cout/streamer problem for mPhiWgt under Solaris
00676 //
00677 // Revision 1.7  2000/05/26 21:29:27  posk
00678 // Protected Track data members from overflow.
00679 //
00680 // Revision 1.5  2000/05/16 20:59:30  posk
00681 // Voloshin's flownanoevent.root added.
00682 //
00683 // Revision 1.4  2000/05/12 22:42:04  snelling
00684 // Additions for persistency and minor fix
00685 //
00686 // Revision 1.2  2000/03/15 23:28:51  posk
00687 // Added StFlowSelection.
00688 //
00689 // Revision 1.1  2000/03/02 23:02:50  posk
00690 // Changed extensions from .hh and .cc to .h and .cxx .
00691 //
00692 // Revision 1.10  2000/02/29 22:00:54  posk
00693 // Made SetPhiWeight inline, changed ImpactPar to Dca, etc.
00694 //
00695 // Revision 1.9  2000/02/18 22:49:55  posk
00696 // Added PID and centrality.
00697 //
00698 // Revision 1.5  1999/12/15 22:01:26  posk
00699 // Added StFlowConstants.hh
00700 //
00701 // Revision 1.4  1999/12/04 00:10:33  posk
00702 // Works with the new StEvent
00703 //
00704 // Revision 1.3  1999/11/30 18:52:52  snelling
00705 // First modification for the new StEvent
00706 //
00707 // Revision 1.2  1999/11/24 18:17:14  posk
00708 // Put the methods which act on the data in with the data in StFlowEvent.
00709 //
00710 // Revision 1.1  1999/11/04 19:02:06  snelling
00711 // First check in of StFlowMaker. It contains the common code from
00712 // StFlowTagMaker and StFlowAnalysisMaker.
00713 //