StRoot: StHbtMaker/Infrastructure/StHbtReactionPlaneAnalysis.cxx Source File

00001 /***************************************************************************
00002  *
00003  * $Id: StHbtReactionPlaneAnalysis.cxx,v 1.6 2005/06/28 23:12:53 chajecki Exp $
00004  *
00005  * Author: Randall Wells, Ohio State, rcwells@mps.ohio-state.edu
00006  ***************************************************************************
00007  *
00008  * Description: part of STAR HBT Framework: StHbtMaker package
00009  *      This is the Class for Analysis objects.  Each of the simultaneous
00010  *      Analyses running should have one of these instantiated.  They link
00011  *      into the Manager in an Analysis Collection.
00012  *
00013  ***************************************************************************
00014  *
00015  * $Log: StHbtReactionPlaneAnalysis.cxx,v $
00016  * Revision 1.6  2005/06/28 23:12:53  chajecki
00017  * UncorrectedNumberOfNegativePrimaries() -> UncorrectedNumberOfPrimaries()
00018  *
00019  * For data taken in Y2 and later the centrality definition bases
00020  * on UncorrectedNumberOfPrimaries() while for Y1(AuAu@130)
00021  * it based on UncorrectedNumberOfNegativePrimaries().
00022  * But in many places of HBT code the number of negative primaries
00023  * was used as a multiplicity for all productions.
00024  * This has been fixed.
00025  *
00026  * Revision 1.5  2004/04/12 14:05:46  magestro
00027  * Added Vz dimension to event-mixing
00028  *
00029  * Revision 1.4  2002/05/28 14:04:07  rcwells
00030  * Added multiplicity binning to StHbtReactionPlaneAnalysis
00031  *
00032  * Revision 1.2  2001/07/20 20:03:53  rcwells
00033  * Added pT weighting and moved event angle cal. to event cut
00034  *
00035  * Revision 1.1  2001/07/13 20:03:14  rcwells
00036  * Adding reaction plane analysis
00037  *
00038  *
00039  **************************************************************************/
00040 
00041 #include "StHbtMaker/Infrastructure/StHbtReactionPlaneAnalysis.h"
00042 #include "StHbtMaker/Infrastructure/StHbtParticleCollection.hh"
00043 #include "StHbtMaker/Base/StHbtTrackCut.h"
00044 #include "StHbtMaker/Base/StHbtV0Cut.h"
00045 #include "StHbtMaker/Infrastructure/StHbtPicoEventCollectionVector.hh"
00046 #include "StHbtMaker/Infrastructure/StHbtPicoEventCollectionVectorHideAway.hh"
00047 #include "PhysicalConstants.h"
00048 
00049 #ifdef __ROOT__ 
00050 ClassImp(StHbtReactionPlaneAnalysis)
00051 #endif
00052 
00053 extern void FillHbtParticleCollection(StHbtParticleCut*         partCut,
00054                                      StHbtEvent*               hbtEvent,
00055                                      StHbtParticleCollection*  partCollection);
00056 
00057 
00058 //____________________________
00059 StHbtReactionPlaneAnalysis::StHbtReactionPlaneAnalysis(int PtWgt, unsigned int binsRP, double minRP, double maxRP,
00060                                                        unsigned int binsMult, double minMult, double maxMult, 
00061                    unsigned int binsVert, double minVert, double maxVert){
00062   //  mControlSwitch     = 0;
00063   mEventCut          = 0;
00064   mFirstParticleCut  = 0;
00065   mSecondParticleCut = 0;
00066   mPairCut           = 0;
00067   mCorrFctnCollection= 0;
00068   mCorrFctnCollection = new StHbtCorrFctnCollection;
00069   mReactionPlaneBins = binsRP;
00070   mReactionPlane[0] = minRP;
00071   mReactionPlane[1] = maxRP;
00072   mUnderFlow = 0; 
00073   mOverFlow = 0; 
00074   if (mMixingBuffer) delete mMixingBuffer;
00075   mPicoEventCollectionVectorHideAway = new StHbtPicoEventCollectionVectorHideAway(mReactionPlaneBins,mReactionPlane[0],mReactionPlane[1],binsMult,minMult,maxMult
00076                   ,binsVert,minVert,maxVert);
00077   mPtWgt = PtWgt;
00078 };
00079 //____________________________
00080 
00081 StHbtReactionPlaneAnalysis::StHbtReactionPlaneAnalysis(const StHbtReactionPlaneAnalysis& a) : StHbtAnalysis() {
00082   //StHbtReactionPlaneAnalysis();
00083   mPtWgt = a.mPtWgt;
00084   mEventCut          = 0;
00085   mFirstParticleCut  = 0;
00086   mSecondParticleCut = 0;
00087   mPairCut           = 0;
00088   mCorrFctnCollection= 0;
00089   mCorrFctnCollection = new StHbtCorrFctnCollection;
00090   mReactionPlaneBins = a.mReactionPlaneBins; 
00091   mReactionPlane[0] = a.mReactionPlane[0]; 
00092   mReactionPlane[1] = a.mReactionPlane[1];
00093   mUnderFlow = 0; 
00094   mOverFlow = 0; 
00095   if (mMixingBuffer) delete mMixingBuffer;
00096   mPicoEventCollectionVectorHideAway = new StHbtPicoEventCollectionVectorHideAway(mReactionPlaneBins,mReactionPlane[0],mReactionPlane[1]);
00097 
00098   // find the right event cut
00099   mEventCut = a.mEventCut->Clone();
00100   // find the right first particle cut
00101   mFirstParticleCut = a.mFirstParticleCut->Clone();
00102   // find the right second particle cut
00103   if (a.mFirstParticleCut==a.mSecondParticleCut) 
00104     SetSecondParticleCut(mFirstParticleCut); // identical particle hbt
00105   else
00106   mSecondParticleCut = a.mSecondParticleCut->Clone();
00107 
00108   mPairCut = a.mPairCut->Clone();
00109   
00110   if ( mEventCut ) {
00111       SetEventCut(mEventCut); // this will set the myAnalysis pointer inside the cut
00112       cout << " StHbtReactionPlaneAnalysis::StHbtReactionPlaneAnalysis(const StHbtReactionPlaneAnalysis& a) - event cut set " << endl;
00113   }
00114   if ( mFirstParticleCut ) {
00115       SetFirstParticleCut(mFirstParticleCut); // this will set the myAnalysis pointer inside the cut
00116       cout << " StHbtReactionPlaneAnalysis::StHbtReactionPlaneAnalysis(const StHbtReactionPlaneAnalysis& a) - first particle cut set " << endl;
00117   }
00118   if ( mSecondParticleCut ) {
00119       SetSecondParticleCut(mSecondParticleCut); // this will set the myAnalysis pointer inside the cut
00120       cout << " StHbtReactionPlaneAnalysis::StHbtReactionPlaneAnalysis(const StHbtReactionPlaneAnalysis& a) - second particle cut set " << endl;
00121   }  if ( mPairCut ) {
00122       SetPairCut(mPairCut); // this will set the myAnalysis pointer inside the cut
00123       cout << " StHbtReactionPlaneAnalysis::StHbtReactionPlaneAnalysis(const StHbtReactionPlaneAnalysis& a) - pair cut set " << endl;
00124   }
00125 
00126   StHbtCorrFctnIterator iter;
00127   for (iter=a.mCorrFctnCollection->begin(); iter!=a.mCorrFctnCollection->end();iter++){
00128     cout << " StHbtReactionPlaneAnalysis::StHbtReactionPlaneAnalysis(const StHbtReactionPlaneAnalysis& a) - looking for correlation functions " << endl;
00129     StHbtCorrFctn* fctn = (*iter)->Clone();
00130     if (fctn) AddCorrFctn(fctn);
00131     else cout << " StHbtReactionPlaneAnalysis::StHbtReactionPlaneAnalysis(const StHbtReactionPlaneAnalysis& a) - correlation function not found " << endl;
00132   }
00133 
00134   mNumEventsToMix = a.mNumEventsToMix;
00135 
00136   cout << " StHbtReactionPlaneAnalysis::StHbtReactionPlaneAnalysis(const StHbtReactionPlaneAnalysis& a) - analysis copied " << endl;
00137 
00138 }
00139 //____________________________
00140 StHbtReactionPlaneAnalysis::~StHbtReactionPlaneAnalysis(){
00141   // now delete every PicoEvent in the EventMixingBuffer and then the Buffer itself
00142   delete mPicoEventCollectionVectorHideAway;
00143 }
00144 
00145 //____________________________
00146 StHbtString StHbtReactionPlaneAnalysis::Report()
00147 {
00148   cout << "StHbtReactionPlaneAnalysis - constructing Report..."<<endl;
00149   char Ctemp[200];
00150   StHbtString temp = "-----------\nHbt StHbtReactionPlaneAnalysis Report:\n";
00151   sprintf(Ctemp,"Events are mixed in %d bins in the range %E cm to %E cm.\n",mReactionPlaneBins,mReactionPlane[0],mReactionPlane[1]);
00152   temp += Ctemp;
00153   sprintf(Ctemp,"Events underflowing: %d\n",mUnderFlow);
00154   temp += Ctemp;
00155   sprintf(Ctemp,"Events overflowing: %d\n",mOverFlow);
00156   temp += Ctemp;
00157   sprintf(Ctemp,"Now adding StHbtAnalysis(base) Report\n");
00158   temp += Ctemp;
00159   temp += StHbtAnalysis::Report();
00160   StHbtString returnThis=temp;
00161   return returnThis;
00162 }
00163 //_________________________
00164 void StHbtReactionPlaneAnalysis::ProcessEvent(const StHbtEvent* hbtEvent) {
00165   // get right mixing buffer
00166   mVertexZ = hbtEvent->PrimVertPos().z();
00167   if (mPtWgt) mReactionPlaneAngle = hbtEvent->ReactionPlane(1);
00168   else mReactionPlaneAngle = hbtEvent->ReactionPlane(0);
00169   // cout << "Reaction Plane " << mReactionPlaneAngle << endl; 
00170   if (mReactionPlaneAngle<-990) {
00171     cout << "No event plane!" << endl;
00172     return;
00173   }
00174   if ( mReactionPlaneAngle<0.0 ) mReactionPlaneAngle+=2*pi;
00175   int multiplicity = hbtEvent->UncorrectedNumberOfPrimaries();
00176   mMixingBuffer = mPicoEventCollectionVectorHideAway->PicoEventCollection(mReactionPlaneAngle,(double)multiplicity, mVertexZ); 
00177   if (!mMixingBuffer) {
00178     if ( mReactionPlaneAngle < mReactionPlane[0] ) mUnderFlow++;
00179     if ( mReactionPlaneAngle > mReactionPlane[1] ) mOverFlow++;
00180     return;
00181   }
00182   // call ProcessEvent() from StHbtAnalysis-bas
00183   StHbtAnalysis::ProcessEvent(hbtEvent);
00184 }