StRoot: StFlowAnalysisMaker/StFlowAnalysisMaker.h Source File

00001 
00002 //
00003 // $Id: StFlowAnalysisMaker.h,v 1.51 2011/03/10 18:56:25 posk Exp $
00004 //
00005 // Authors: Art Poskanzer and Raimond Snellings, LBNL, Aug 1999
00006 //          FTPC added by Markus Oldenburg, MPI, Dec 2000
00007 //
00009 //
00010 // Description:  Maker to analyze Flow using StFlowEvent
00011 //
00013 
00014 #ifndef StFlowAnalysisMaker_H
00015 #define StFlowAnalysisMaker_H
00016 #include <Stiostream.h>
00017 #include "StMaker.h"
00018 #include "StFlowMaker/StFlowConstants.h"
00019 #include "TVector2.h"
00020 #include "TString.h"
00021 class StFlowMaker;
00022 class StFlowEvent;
00023 class StFlowSelection;
00024 class TH1F;
00025 class TH1D;
00026 class TH2F;
00027 class TH2D;
00028 class TH3F;
00029 class TProfile;
00030 class TProfile2D;
00031 
00032 class StFlowAnalysisMaker : public StMaker {
00033 
00042 public:
00043 
00045            StFlowAnalysisMaker(const Char_t* name="FlowAnalysis");
00047           StFlowAnalysisMaker(const Char_t* name,
00048                                const StFlowSelection& pFlowSelect);
00049            StFlowAnalysisMaker(const StFlowAnalysisMaker &from) {};
00050   virtual  ~StFlowAnalysisMaker();
00051 
00052   Int_t    Init();
00053   Int_t    Make();
00054   Int_t    Finish();
00055   void     SetHistoRanges(Bool_t ftpc_included = kFALSE);
00056   void     SetPtRange_for_vEta(Float_t lo, Float_t hi);
00057   void     SetEtaRange_for_vPt(Float_t lo, Float_t hi);
00058   void     SetV1Ep1Ep2(Bool_t v1Ep1Ep2 = kTRUE);
00059   virtual  const char *GetCVS() const {static const char cvs[]=
00060     "Tag $Name:  $ $Id: StFlowAnalysisMaker.h,v 1.51 2011/03/10 18:56:25 posk Exp $ built "__DATE__" "__TIME__ ;
00061     return cvs;}
00062 
00063 private:
00064 
00065   Bool_t   FillFromFlowEvent();
00066   void     FillEventHistograms();
00067   void     FillParticleHistograms();
00068   Bool_t   mCalcReCentPars;
00069 #ifndef __CINT__
00070   TVector2 mQ[Flow::nSels][Flow::nHars];                     
00071   Float_t  mPsi[Flow::nSels][Flow::nHars];                   
00072   UInt_t   mMult[Flow::nSels][Flow::nHars];                  
00073   Float_t  m_q[Flow::nSels][Flow::nHars];                    
00074   TVector2 mQSub[Flow::nSels*Flow::nSubs][Flow::nHars];      
00075   Float_t  mPsiSub[Flow::nSels*Flow::nSubs][Flow::nHars];    
00076   UInt_t   mMultSub[Flow::nSels*Flow::nSubs][Flow::nHars];   
00077   Float_t  mRes[Flow::nSels][Flow::nHars];      
00078   Float_t  mResErr[Flow::nSels][Flow::nHars];   
00079   Float_t  mZDCSMD_e_PsiWgt,mZDCSMD_w_PsiWgt,mZDCSMD_f_PsiWgt;   
00080   Float_t  mFlowWeight;                         
00081 #endif /*__CINT__*/
00082   TString          xLabel;      
00083   StFlowMaker*     pFlowMaker;  
00084   StFlowEvent*     pFlowEvent;  
00085   StFlowSelection* pFlowSelect; 
00086 
00087   // for single histograms
00088   TH1F*     mHistTrigger;              
00089   TH1F*     mHistChargeFtpc;           
00090   TH1F*     mHistDcaGlobalTpc;         
00091   TH1F*     mHistDcaGlobalFtpc;        
00092   TH1F*     mHistDcaTpc;               
00093   TH1F*     mHistDcaFtpc;              
00094   TH1F*     mHistChi2Tpc;              
00095   TH1F*     mHistChi2Ftpc;             
00096   TH1F*     mHistFitPtsTpc;            
00097   TH1F*     mHistMaxPtsTpc;            
00098   TH1F*     mHistFitOverMaxTpc;        
00099   TH1F*     mHistFitPtsFtpc;           
00100   TH1F*     mHistMaxPtsFtpc;           
00101   TH1F*     mHistFitOverMaxFtpc;       
00102   TH1F*     mHistOrigMult;             
00103   TH1F*     mHistMultEta;              
00104   TH1F*     mHistMult;                 
00105   TH1F*     mHistMultOverOrig;         
00106   TH1F*     mHistMultPart;             
00107   TH1F*     mHistVertexZ;              
00108   TH2F*     mHistVertexXY2D;           
00109   TH1F*     mHistEtaSymTpc;            
00110   TH1F*     mHistEtaSymFtpc;           
00111   TH3F*     mHistEtaPtPhi3D;           
00112   TH2D*     mHistYieldAll2D;           
00113   TH2D*     mHistYieldPart2D;          
00114   TProfile* mHistBinEta;               
00115   TProfile* mHistBinPt;                
00116   TH1F*     mHistPidPiPlus;            
00117   TH1F*     mHistPidPiMinus;           
00118   TH1F*     mHistPidProton;            
00119   TH1F*     mHistPidAntiProton;        
00120   TH1F*     mHistPidKplus;             
00121   TH1F*     mHistPidKminus;            
00122   TH1F*     mHistPidDeuteron;          
00123   TH1F*     mHistPidAntiDeuteron;      
00124   TH1F*     mHistPidElectron;          
00125   TH1F*     mHistPidPositron;          
00126   TH1F*     mHistPidPiPlusPart;        
00127   TH1F*     mHistPidPiMinusPart;       
00128   TH1F*     mHistPidProtonPart;        
00129   TH1F*     mHistPidAntiProtonPart;    
00130   TH1F*     mHistPidKplusPart;         
00131   TH1F*     mHistPidKminusPart;        
00132   TH1F*     mHistPidDeuteronPart;      
00133   TH1F*     mHistPidAntiDeuteronPart;  
00134   TH1F*     mHistPidElectronPart;      
00135   TH1F*     mHistPidPositronPart;      
00136   TProfile* mHistPidMult;              
00137   TH1F*     mHistCent;                 
00138   TH2F*     mHistEtaSymVerZ2DTpc;      
00139   TH2F*     mHistEtaSymVerZ2DFtpc;     
00140   TH2F*     mHistCTBvsZDC2D;           
00141   TH2F*     mHistMeanDedxPos2D;        
00142   TH2F*     mHistMeanDedxNeg2D;        
00143   TH2F*     mHistMeanDedxPiPlus2D;     
00144   TH2F*     mHistMeanDedxPiMinus2D;    
00145   TH2F*     mHistMeanDedxProton2D;     
00146   TH2F*     mHistMeanDedxPbar2D;       
00147   TH2F*     mHistMeanDedxKplus2D;      
00148   TH2F*     mHistMeanDedxKminus2D;     
00149   TH2F*     mHistMeanDedxDeuteron2D;   
00150   TH2F*     mHistMeanDedxAntiDeuteron2D;
00151   TH2F*     mHistMeanDedxPositron2D;   
00152   TH2F*     mHistMeanDedxElectron2D;   
00153   TH1F*     mZDC_SMD_west_vert;        
00154   TH1F*     mZDC_SMD_west_hori;        
00155   TH1F*     mZDC_SMD_east_vert;        
00156   TH1F*     mZDC_SMD_east_hori;        
00157   TH1D*     mHistZDCSMDPsiWgtEast;     
00158   TH1D*     mHistZDCSMDPsiWgtWest;     
00159   TH1D*     mHistZDCSMDPsiWgtTest;     
00160   TH1D*     mHistZDCSMDPsiWgtFull;     
00161   TH1D*     mHistZDCSMDPsiCorTest;     
00162   TH1D*     mHistZDCSMDPsiCorFull;     
00163   
00164   // for each harmonic, each selection, and each sub-event
00165   struct histSubHars {
00166     TH1F*     mHistPsiSubs;
00167   };
00168   struct histSubs;      
00169   friend struct histSubs;
00170   struct histSubs {
00171     struct histSubHars histSubHar[Flow::nHars];
00172   };
00173   struct histSubs histSub[Flow::nSels*Flow::nSubs]; 
00174 
00175   // for each harmonic and each selection
00176   struct histFullHars {
00177     TH1F*       mHistPhiCorr;
00178     TH1F*       mHistPhiLab;
00179     TH1F*       mHistPsiSubCorr;
00180     TH1F*       mHistPsiSubCorrDiff;
00181     TH1F*       mHistPsi;
00182     TProfile*   mHistReCentX;
00183     TProfile*   mHistReCentY;
00184     TProfile*   mHistQreCent;
00185     TH2D*       mHistQXY2D;
00186     TH2D*       mHistQFTPCSubXY2D;
00187     TH2D*       mHistQTPCSubXY2D;
00188     TH1F*       mHistPsi_Diff;
00189     TH1F*       mHistMult;
00190     TH1F*       mHist_q;
00191     TH2D*       mHistYield2D;
00192     TProfile2D* mHist_vObs2D;
00193     TProfile*   mHist_vObsEta;
00194     TProfile*   mHist_vObsPt;
00195     TH2D*       mHist_v2D;
00196     TH1D*       mHist_vEta;
00197     TH1D*       mHist_vPt;
00198   };
00199 
00200   // for two harmonics and each selection
00201   struct histTwoHars {
00202     TH1D*       mHistPhiFarEast;
00203     TH1D*       mHistPhiEast;
00204     TH1D*       mHistPhiWest;
00205     TH1D*       mHistPhiFarWest;
00206     TH1D*       mHistPhiFtpcFarEast;
00207     TH1D*       mHistPhiFtpcEast;
00208     TH1D*       mHistPhiFtpcWest;
00209     TH1D*       mHistPhiFtpcFarWest;
00210     TH1D*       mHistPhiWgtFarEast;
00211     TH1D*       mHistPhiWgtEast;
00212     TH1D*       mHistPhiWgtWest;
00213     TH1D*       mHistPhiWgtFarWest;
00214     TH1D*       mHistPhiWgtFtpcFarEast;
00215     TH1D*       mHistPhiWgtFtpcEast;
00216     TH1D*       mHistPhiWgtFtpcWest;
00217     TH1D*       mHistPhiWgtFtpcFarWest;
00218     TH1D*       mHistPhiFlatFarEast;
00219     TH1D*       mHistPhiFlatEast;
00220     TH1D*       mHistPhiFlatWest;
00221     TH1D*       mHistPhiFlatFarWest;
00222     TH1D*       mHistPhiFlatFtpcFarEast;
00223     TH1D*       mHistPhiFlatFtpcEast;
00224     TH1D*       mHistPhiFlatFtpcWest;
00225     TH1D*       mHistPhiFlatFtpcFarWest;
00226   };
00227 
00228   // for each selection
00229   struct histFulls;     
00230   friend struct histFulls;
00231   struct histFulls {
00232     TProfile* mHistCos;
00233     TH1F*     mHistRes;
00234     TProfile* mHist_vObs;
00235     TH1D*     mHist_v;
00236     struct histFullHars histFullHar[Flow::nHars];
00237     struct histTwoHars histTwoHar[2];
00238   };
00239   struct histFulls histFull[Flow::nSels]; 
00240 
00241   TString      MakerName;
00242 
00243   Float_t mEtaMin;
00244   Float_t mEtaMax;
00245     Int_t mNEtaBins;
00246 
00247   Float_t mPtRange_for_vEta[2];
00248   Float_t mEtaRange_for_vPt[2];
00249 
00250   Bool_t  mV1Ep1Ep2;
00251 
00252   ClassDef(StFlowAnalysisMaker,0)              // macro for rootcint
00253 };
00254 #endif
00255 
00257 //
00258 // $Log: StFlowAnalysisMaker.h,v $
00259 // Revision 1.51  2011/03/10 18:56:25  posk
00260 // Added histogram for laboratory azimuthal distribution of particles.
00261 //
00262 // Revision 1.50  2010/09/30 19:28:12  posk
00263 // Instead of reversing the weight for negative pseudrapidity for odd harmonics,
00264 // it is now done only for the first harmonic.
00265 // Recentering is now done for all harmonics.
00266 //
00267 // Revision 1.49  2009/11/24 19:29:12  posk
00268 // Added reCenter to remove acceptance correlations as an option instead of phiWgt.
00269 //
00270 // Revision 1.48  2006/02/22 19:36:27  posk
00271 // Minor updates.
00272 //
00273 // Revision 1.47  2004/12/17 22:33:33  aihong
00274 // add in full Psi weight for ZDC SMD and fix a few bugs, done by Gang
00275 //
00276 // Revision 1.46  2004/12/09 23:47:07  posk
00277 // Minor changes in code formatting.
00278 // Added hist for TPC primary dca to AnalysisMaker.
00279 //
00280 // Revision 1.45  2004/12/07 23:10:21  posk
00281 // Only odd and even phiWgt hists. If the old phiWgt file contains more than
00282 // two harmonics, only the first two are read. Now writes only the first two.
00283 //
00284 // Revision 1.44  2004/05/05 21:13:47  aihong
00285 // Gang's code for ZDC-SMD added
00286 //
00287 // Revision 1.43  2003/11/14 20:00:43  oldi
00288 // Implementation of v1{EP1,EP2}. This method is set to be the default for v1 now!
00289 // Minor code clean-ups.
00290 //
00291 // Revision 1.42  2003/09/10 19:47:14  perev
00292 // ansi corrs
00293 //
00294 // Revision 1.41  2003/09/02 17:58:10  perev
00295 // gcc 3.2 updates + WarnOff
00296 //
00297 // Revision 1.40  2003/08/06 20:54:09  oldi
00298 // Introduction of possibility to exclude pt ranges for v(eta) and eta regions
00299 // for v(pt) histograms. Default behavior stays the same (all available tracks
00300 // are included in v(pt) and v(eta)).
00301 //
00302 // Revision 1.39  2003/05/06 21:33:06  posk
00303 // Removed some histograms.
00304 //
00305 // Revision 1.38  2003/05/06 18:38:06  posk
00306 // Removed StFlowTagMaker.
00307 //
00308 // Revision 1.37  2003/01/10 16:40:33  oldi
00309 // Several changes to comply with FTPC tracks:
00310 // - Switch to include/exclude FTPC tracks introduced.
00311 //   The same switch changes the range of the eta histograms.
00312 // - Eta symmetry plots for FTPC tracks added and separated from TPC plots.
00313 // - PhiWgts and related histograms for FTPC tracks split in FarEast, East,
00314 //   West, FarWest (depending on vertex.z()).
00315 // - Psi_Diff plots for 2 different selections and the first 2 harmonics added.
00316 // - Cut to exclude mu-events with no primary vertex introduced.
00317 //   (This is possible for UPC events and FTPC tracks.)
00318 // - Global DCA cut for FTPC tracks added.
00319 // - Global DCA cuts for event plane selection separated for TPC and FTPC tracks.
00320 // - Charge cut for FTPC tracks added.
00321 //
00322 // Revision 1.36  2002/11/26 22:11:53  posk
00323 // First use of doxygen.
00324 //
00325 // Revision 1.35  2002/10/28 19:45:55  posk
00326 // Eliminate events with Psi=0.
00327 //
00328 // Revision 1.34  2001/12/11 22:03:56  posk
00329 // Four sets of phiWgt histograms.
00330 // StFlowMaker StFlowEvent::PhiWeight() changes.
00331 // Cumulant histogram names changed.
00332 //
00333 // Revision 1.33  2001/11/09 21:14:46  posk
00334 // Switched from CERNLIB to TMath. Using global dca instead of dca.
00335 //
00336 // Revision 1.32  2001/08/02 17:41:50  snelling
00337 // Added trigger histogram
00338 //
00339 // Revision 1.31  2001/05/22 20:11:15  posk
00340 // Changed dEdx graphs.
00341 //
00342 // Revision 1.30  2001/04/25 17:45:52  perev
00343 // HPcorrs
00344 //
00345 // Revision 1.29  2000/12/12 15:01:11  posk
00346 // Put log comments at end of file.
00347 //
00348 // Revision 1.28  2000/12/10 02:02:02  oldi
00349 // A new member (StTrackTopologyMap mTopology) was added to StFlowPicoTrack.
00350 // The evaluation of either a track originates from the FTPC or not is
00351 // unambiguous now. The evaluation itself is easily extendible for other
00352 // detectors (e.g. SVT+TPC). Old flowpicoevent.root files are treated as if
00353 // they contain TPC tracks only (backward compatibility).
00354 //
00355 // Revision 1.27  2000/12/08 17:04:09  oldi
00356 // Phi weights for both FTPCs included.
00357 //
00358 // Revision 1.24  2000/09/22 22:01:40  posk
00359 // Doubly integrated v now contains resolution error.
00360 //
00361 // Revision 1.23  2000/09/16 22:23:06  snelling
00362 // Auto magically switch to rapidity when identified particles are used
00363 //
00364 // Revision 1.22  2000/09/15 22:52:55  posk
00365 // Added Pt weighting for event plane calculation.
00366 //
00367 // Revision 1.21  2000/09/12 01:31:00  snelling
00368 // Added pid histograms for e- e+ and dbar
00369 //
00370 // Revision 1.20  2000/08/31 18:50:30  posk
00371 // Added plotCen.C to plot from a series of files with different centralities.
00372 //
00373 // Revision 1.19  2000/08/09 21:38:59  snelling
00374 // Added monitor histograms
00375 //
00376 // Revision 1.18  2000/08/01 21:51:20  posk
00377 // Added doubly integrated v.
00378 //
00379 // Revision 1.17  2000/07/03 02:07:49  perev
00380 // StEvent: vector<TObject*>
00381 //
00382 // Revision 1.16  2000/06/30 14:51:19  posk
00383 // Using MessageMgr. Added graph for Eta Symmetry vs. Vertex Z.
00384 //
00385 // Revision 1.15  2000/05/26 21:25:22  posk
00386 // Use TProfile2D class and profile projection methods.
00387 // Correction needed for >2 subevents.
00388 //
00389 // Revision 1.14  2000/05/16 17:30:21  snelling
00390 // removed the dependencies cint did not like
00391 //
00392 // Revision 1.13  2000/05/03 16:38:35  posk
00393 // Compatable with ROOT 2.24/02.
00394 //
00395 // Revision 1.12  2000/04/13 22:34:15  posk
00396 // Resolution correction is now made.
00397 //
00398 // Revision 1.11  2000/03/21 00:24:45  posk
00399 // Added GetCVS and changed some plot names.
00400 //
00401 // Revision 1.10  2000/03/15 23:32:04  posk
00402 // Added StFlowSelection.
00403 //
00404 // Revision 1.9  2000/03/02 22:55:34  posk
00405 // Changed header file extensions from .hh to .h .
00406 //
00407 // Revision 1.11  2000/02/29 21:55:14  posk
00408 // Removed static const int& statements.
00409 //
00410 // Revision 1.10  2000/02/18 23:44:53  posk
00411 // Added PID and centrality.
00412 //
00413 // Revision 1.9  2000/02/04 16:26:42  posk
00414 // Added correct calculation of event plane resolution for large flow.
00415 //
00416 // Revision 1.7  2000/01/14 01:35:52  snelling
00417 // changed include path ../FlowMaker/ to FlowMaker/
00418 //
00419 // Revision 1.3  1999/12/04 00:15:40  posk
00420 // Works with StFlowEvent which works with the new StEvent
00421 //
00422 // Revision 1.2  1999/11/24 18:14:07  posk
00423 // Now reads event quantities with StFlowEvent methods
00424 //
00425 // Revision 1.7  1999/11/05 00:02:03  posk
00426 // Changed the flow vector, Q, to a TVector2.
00427 //
00428 // Revision 1.6  1999/10/05 16:54:11  posk
00429 // Added getPhiWeight method for making the event plane isotropic.
00430 //