ClassImp(AliFlowEventSimpleMaker)
//-----------------------------------------------------------------------
AliFlowEventSimpleMaker::AliFlowEventSimpleMaker() :
- fMCReactionPlaneAngle(0.)
+ fMCReactionPlaneAngle(0.),
+ fCount(0)
{
//constructor
}
Int_t iPIDDiff = diffCuts->GetPID();
Int_t iNumberOfInputTracks = anInput->GetEntries() ;
- //cerr<<"iNumberOfInputTracks = "<<iNumberOfInputTracks<<endl;
+
TParticle* pParticle = new TParticle();
anInput->SetBranchAddress("Particles",&pParticle);
// AliFlowEventSimple* pEvent = new AliFlowEventSimple(iNumberOfInputTracks);
AliFlowEventSimple* pEvent = new AliFlowEventSimple(10);
- //cerr<<pEvent<<" pEvent "<<endl;
Int_t iN = iNumberOfInputTracks; // additional variable to artificially fix the number of tracks
// Int_t iN = 576; //multiplicity for chi=1.5
pEvent->SetEventNSelTracksIntFlow(iSelParticlesInt);
pEvent->SetNumberOfTracks(iGoodTracks);//tracks used either for int. or for diff. flow
-
pEvent->SetMCReactionPlaneAngle(fMCReactionPlaneAngle);
+
if (!fMCReactionPlaneAngle == 0) cout<<" MC Reaction Plane Angle = "<< fMCReactionPlaneAngle << endl;
else cout<<" MC Reaction Plane Angle = unknown "<< endl;
cout<<" iGoodTracks = "<< iGoodTracks << endl;
cout<<" # of selected tracks for int. flow = "<<iSelParticlesInt<<endl;
cout<<" # of selected tracks for diff. flow = "<<iSelParticlesDiff<<endl;
+ cout << "# " << ++fCount << " events processed" << endl;
delete pParticle;
return pEvent;
//Fills the event from the MC kinematic information
Int_t iNumberOfInputTracks = anInput->GetNumberOfTracks() ;
- cerr<<"anInput->GetNumberOfTracks() = "<<iNumberOfInputTracks<<endl;
AliFlowEventSimple* pEvent = new AliFlowEventSimple(10);
cout<<" iGoodTracks = "<<iGoodTracks<<endl;
cout<<" # of selected tracks for int. flow = "<<iSelParticlesInt<<endl;
cout<<" # of selected tracks for diff. flow = "<<iSelParticlesDiff<<endl;
+ cout << "# " << ++fCount << " events processed" << endl;
- //cout << " iSelectedTracksInt = " << iSelParticlesInt << endl;
return pEvent;
}
//Fills the event from the MC kinematic information
Int_t iNumberOfInputTracks = anInput->GetNumberOfTracks() ;
- cerr<<"anInput->GetNumberOfTracks() = "<<iNumberOfInputTracks<<endl;
+
if (iNumberOfInputTracks==-1) {
cout<<"Skipping Event -- No MC information available for this event"<<endl;
return 0;
cout<<" iGoodTracks = "<<iGoodTracks<<endl;
cout<<" # of selected tracks for int. flow = "<<iSelParticlesInt<<endl;
cout<<" # of selected tracks for diff. flow = "<<iSelParticlesDiff<<endl;
+ cout << "# " << ++fCount << " events processed" << endl;
- //cout << " iSelectedTracksInt = " << iSelParticlesInt << endl;
return pEvent;
}
//Fills the event from the ESD
Int_t iNumberOfInputTracks = anInput->GetNumberOfTracks() ;
- cerr<<"anInput->GetNumberOfTracks() = "<<iNumberOfInputTracks<<endl;
AliFlowEventSimple* pEvent = new AliFlowEventSimple(10);
cout<<" iGoodTracks = "<<iGoodTracks<<endl;
cout<<" # of selected tracks for int. flow = "<<iSelParticlesInt<<endl;
cout<<" # of selected tracks for diff. flow = "<<iSelParticlesDiff<<endl;
+ cout << "# " << ++fCount << " events processed" << endl;
- //cout << " iSelectedTracksInt = " << iSelParticlesInt << endl;
return pEvent;
-
-
}
//Fills the event from the ESD
Int_t iNumberOfInputTracks = anInput->GetNumberOfTracks() ;
- cerr<<"anInput->GetNumberOfTracks() = "<<iNumberOfInputTracks<<endl;
AliFlowEventSimple* pEvent = new AliFlowEventSimple(10);
cout<<" iGoodTracks = "<<iGoodTracks<<endl;
cout<<" # of selected tracks for int. flow = "<<iSelParticlesInt<<endl;
cout<<" # of selected tracks for diff. flow = "<<iSelParticlesDiff<<endl;
+ cout << "# " << ++fCount << " events processed" << endl;
- //cout << " iSelectedTracksInt = " << iSelParticlesInt << endl;
return pEvent;
-
-
}
//-----------------------------------------------------------------------
//Fills the event from the AOD
Int_t iNumberOfInputTracks = anInput->GetNumberOfTracks() ;
- cerr<<"anInput->GetNumberOfTracks() = "<<iNumberOfInputTracks<<endl;
AliFlowEventSimple* pEvent = new AliFlowEventSimple(10);
cout<<" iGoodTracks = "<<iGoodTracks<<endl;
cout<<" # of selected tracks for int. flow = "<<iSelParticlesInt<<endl;
cout<<" # of selected tracks for diff. flow = "<<iSelParticlesDiff<<endl;
+ cout << "# " << ++fCount << " events processed" << endl;
- //cout << " iSelectedTracksInt = " << iSelParticlesInt << endl;
return pEvent;
-
}
//Fills the event from the AOD
Int_t iNumberOfInputTracks = anInput->GetNumberOfTracks() ;
- cerr<<"anInput->GetNumberOfTracks() = "<<iNumberOfInputTracks<<endl;
AliFlowEventSimple* pEvent = new AliFlowEventSimple(10);
cout<<" iGoodTracks = "<<iGoodTracks<<endl;
cout<<" # of selected tracks for int. flow = "<<iSelParticlesInt<<endl;
cout<<" # of selected tracks for diff. flow = "<<iSelParticlesDiff<<endl;
+ cout << "# " << ++fCount << " events processed" << endl;
- //cout << " iSelectedTracksInt = " << iSelParticlesInt << endl;
return pEvent;
-
}
//-----------------------------------------------------------------------
}
Int_t iNumberOfInputTracks = anInput->GetNumberOfTracks() ;
- cerr<<"anInput->GetNumberOfTracks() = "<<iNumberOfInputTracks<<endl;
AliFlowEventSimple* pEvent = new AliFlowEventSimple(10);
cout<<" iGoodTracks = "<<iGoodTracks<<endl;
cout<<" # of selected tracks for int. flow = "<<iSelParticlesInt<<endl;
cout<<" # of selected tracks for diff. flow = "<<iSelParticlesDiff<<endl;
+ cout << "# " << ++fCount << " events processed" << endl;
- //cout << " iSelectedTracksInt = " << iSelParticlesInt << endl;
return pEvent;
-
-
}
//-----------------------------------------------------------------------
}
Int_t iNumberOfInputTracks = anInput->GetNumberOfTracks() ;
- cerr<<"anInput->GetNumberOfTracks() = "<<iNumberOfInputTracks<<endl;
+
Int_t iNumberOfInputTracksMC = anInputMc->GetNumberOfTracks() ;
- cerr<<"anInputMc->GetNumberOfTracks() = "<<iNumberOfInputTracksMC<<endl;
if (iNumberOfInputTracksMC==-1) {
cout<<"Skipping Event -- No MC information available for this event"<<endl;
return 0;
pEvent-> SetEventNSelTracksIntFlow(iSelParticlesInt);
pEvent->SetNumberOfTracks(iGoodTracks);
+ cout << " Number of MC input tracks = " << iNumberOfInputTracksMC << endl;
+
pEvent->SetMCReactionPlaneAngle(fMCReactionPlaneAngle);
if (!fMCReactionPlaneAngle == 0) cout<<" MC Reaction Plane Angle = "<< fMCReactionPlaneAngle << endl;
else cout<<" MC Reaction Plane Angle = unknown "<< endl;
cout<<" iGoodTracks = "<<iGoodTracks<<endl;
cout<<" # of selected tracks for int. flow = "<<iSelParticlesInt<<endl;
cout<<" # of selected tracks for diff. flow = "<<iSelParticlesDiff<<endl;
+ cout << "# " << ++fCount << " events processed" << endl;
- //cout << " iSelectedTracksInt = " << iSelParticlesInt << endl;
return pEvent;
-
}
-
-
-
-
-