fCommonHists(0),
fCommonHistsRes(0),
fEventNumber(0),
- fQX(0),
- fQY(0),
+// fQX(0),
+// fQY(0),
fQ2sum(0),
fQtheta(0),
- fEvent(0x0),
- fTrack(0x0),
- fLYZEP(0x0)
+// fEvent(0x0),
+ fTrack(NULL),
+ fLYZEP(NULL)
{
// Constructor.
//-----------------------------------------------------------------------
-void AliFlowAnalysisWithLYZEventPlane::Make(AliFlowEventSimple* fEvent, AliFlowLYZEventPlane* fLYZEP) {
+void AliFlowAnalysisWithLYZEventPlane::Make(AliFlowEventSimple* anEvent, AliFlowLYZEventPlane* fLYZEP) {
+ //void AliFlowAnalysisWithLYZEventPlane::Make(AliFlowEventSimple* anEvent) {
//Get the event plane and weight for each event
- if (fEvent) {
+ if (anEvent) {
//fill control histograms
- fCommonHists->FillControlHistograms(fEvent);
+ fCommonHists->FillControlHistograms(anEvent);
//get the Q vector from the FlowEvent
- fQ = fEvent->GetQ();
+ fQ = anEvent->GetQ();
//Weight with the multiplicity
Double_t fQX = 0.;
Double_t fQY = 0.;
//calculate flow
//loop over the tracks of the event
- Int_t fNumberOfTracks = fEvent->NumberOfTracks();
+ Int_t fNumberOfTracks = anEvent->NumberOfTracks();
for (Int_t i=0;i<fNumberOfTracks;i++)
{
- fTrack = fEvent->GetTrack(i) ;
+ fTrack = anEvent->GetTrack(i) ;
if (fTrack){
if (fTrack->UseForDifferentialFlow()) {
Double_t fPhi = fTrack->Phi();
virtual void Init();
virtual void Make(AliFlowEventSimple* fEvent, AliFlowLYZEventPlane* fLYZEP);
+ // virtual void Make(AliFlowEventSimple* anEvent);
virtual void Finish();
// input files
private:
- //AliFlowAnalysisWithLYZEventPlane(const AliFlowAnalysisWithLYZEventPlane& lyz) {} // AliFlowAnalysisWithLYZEventPlane object cannot be copied
- //void operator=(const AliFlowAnalysisWithLYZEventPlane &lyz) {}
+ // AliFlowAnalysisWithLYZEventPlane(const AliFlowAnalysisWithLYZEventPlane& lyz); // AliFlowAnalysisWithLYZEventPlane object cannot be copied
+ // void operator=(const AliFlowAnalysisWithLYZEventPlane &lyz);
TFile* fOutFile; //!
TFile* fFirstRunFile ; //! pointer to file from first run
AliFlowCommonHistResults* fCommonHistsRes; //!
Int_t fEventNumber; // event counter
- Double_t fQX; // local flow vector
- Double_t fQY; // local flow vector
+ // Double_t fQX; // local flow vector
+ // Double_t fQY; // local flow vector
AliFlowVector fQ; // flow vector
TVector2 fQsum; // flow vector sum
Double_t fQ2sum; // flow vector sum squared
Double_t fQtheta; // flow vector projected on ref. angle theta
- AliFlowEventSimple* fEvent ; //!
+ // AliFlowEventSimple* fEvent ; //!
AliFlowTrackSimple* fTrack ; //!
AliFlowLYZEventPlane* fLYZEP ; //!
fEventNumber(0),
fMult(0),
fNbins(0),
- fEvent(0x0),
- fTrack(0x0),
+ fTrack(NULL),
fDebug(kFALSE),
fHistFileName(0),
fHistFile(0),
- fCommonHists(0),
- fCommonHistsRes(0),
- fHistProFlow(0),
- fHistRP(0)
+ fCommonHists(NULL),
+ fCommonHistsRes(NULL),
+ fHistProFlow(NULL),
+ fHistRP(NULL)
{
//-----------------------------------------------------------------------
-void AliFlowAnalysisWithMCEventPlane::Make(AliFlowEventSimple* fEvent, Double_t fRP) {
+void AliFlowAnalysisWithMCEventPlane::Make(AliFlowEventSimple* anEvent, Double_t fRP) {
//Calculate v2 from the MC reaction plane
- if (fEvent) {
+ if (anEvent) {
//fill control histograms
- fCommonHists->FillControlHistograms(fEvent);
+ fCommonHists->FillControlHistograms(anEvent);
//get the Q vector from the FlowEvent
- fQ = fEvent->GetQ();
+ fQ = anEvent->GetQ();
//cout<<"fQ.Mod() = " << fQ.Mod() << endl;
//for chi calculation:
fQsum += fQ;
//calculate flow
//loop over the tracks of the event
- Int_t fNumberOfTracks = fEvent->NumberOfTracks();
+ Int_t fNumberOfTracks = anEvent->NumberOfTracks();
for (Int_t i=0;i<fNumberOfTracks;i++)
{
- fTrack = fEvent->GetTrack(i) ;
+ fTrack = anEvent->GetTrack(i) ;
if (fTrack){
if (fTrack->UseForDifferentialFlow()) {
Double_t fPhi = fTrack->Phi();
virtual ~AliFlowAnalysisWithMCEventPlane(); //destructor
void Init(); //defines variables and histograms
- void Make(AliFlowEventSimple* fEvent, Double_t fRP); //calculates variables and fills histograms
+ void Make(AliFlowEventSimple* anEvent, Double_t fRP); //calculates variables and fills histograms
void Finish(); //saves histograms
void SetDebug(Bool_t kt) { this->fDebug = kt ; }
Int_t fMult; // multiplicity
Int_t fNbins; // number of bins
- AliFlowEventSimple* fEvent ; //!
AliFlowTrackSimple* fTrack ; //!
Bool_t fDebug ; //! flag for lyz analysis: more print statements
AliFlowAnalysisWithScalarProduct::AliFlowAnalysisWithScalarProduct():
fEventNumber(0),
- fEvent(0x0),
- fTrack(0x0),
+ fTrack(NULL),
fDebug(kFALSE),
fHistFileName(0),
- fHistFile(0),
- fHistProUQ(0),
- fCommonHists(0)
+ fHistFile(NULL),
+ fHistProUQ(NULL),
+ fCommonHists(NULL)
{
// Constructor.
//-----------------------------------------------------------------------
-void AliFlowAnalysisWithScalarProduct::Make(AliFlowEventSimple* fEvent) {
+void AliFlowAnalysisWithScalarProduct::Make(AliFlowEventSimple* anEvent) {
//Fill histogram
- if (fEvent) {
+ if (anEvent) {
//fill control histograms
- fCommonHists->FillControlHistograms(fEvent);
+ fCommonHists->FillControlHistograms(anEvent);
//get the Q vector from the FlowEvent
- fQ = fEvent->GetQ();
+ fQ = anEvent->GetQ();
//Double_t fMult = fQ.GetMult();
//loop over the tracks of the event
- Int_t fNumberOfTracks = fEvent->NumberOfTracks();
+ Int_t fNumberOfTracks = anEvent->NumberOfTracks();
for (Int_t i=0;i<fNumberOfTracks;i++)
{
- fTrack = fEvent->GetTrack(i) ;
+ fTrack = anEvent->GetTrack(i) ;
if (fTrack){
if (fTrack->UseForDifferentialFlow()) {
Double_t fPhi = fTrack->Phi();
cout<<".....finished"<<endl;
}
+
+
virtual ~AliFlowAnalysisWithScalarProduct(); //destructor
void Init(); //defines variables and histograms
- void Make(AliFlowEventSimple* fEvent); //calculates variables and fills histograms
+ void Make(AliFlowEventSimple* anEvent); //calculates variables and fills histograms
void Finish(); //saves histograms
void SetDebug(Bool_t kt) { this->fDebug = kt ; }
Int_t fEventNumber; // event counter
- AliFlowEventSimple* fEvent ; //!
AliFlowTrackSimple* fTrack ; //!
Bool_t fDebug ; //! flag for lyz analysis: more print statements
//-----------------------------------------------------------------------
AliFlowCommonHist::AliFlowCommonHist(TString input):
- fEvent(0),
+// fEvent(0),
fTrack(0),
fHistMultOrig(0),
fHistMultInt(0),
//-----------------------------------------------------------------------
-Bool_t AliFlowCommonHist::FillControlHistograms(AliFlowEventSimple* fEvent)
+Bool_t AliFlowCommonHist::FillControlHistograms(AliFlowEventSimple* Event)
{
//Fills the control histograms
- if (!fEvent){
+ if (!Event){
cout<<"##### FillControlHistograms: FlowEvent pointer null"<<endl;
return kFALSE;
}
//fill the histograms
- Int_t fNumberOfTracks = fEvent->NumberOfTracks();
+ Int_t fNumberOfTracks = Event->NumberOfTracks();
fHistMultOrig->Fill(fNumberOfTracks);
- AliFlowVector fQ = fEvent->GetQ();
+ AliFlowVector fQ = Event->GetQ();
//weight by the Multiplicity
Double_t fQX = fQ.X()/fQ.GetMult();
Double_t fQY = fQ.Y()/fQ.GetMult();
Int_t fMultInt = 0;
Int_t fMultDiff = 0;
for (Int_t i=0;i<fNumberOfTracks;i++) {
- fTrack = fEvent->GetTrack(i);
+ fTrack = Event->GetTrack(i);
if (fTrack ) {
if (fTrack->UseForIntegratedFlow()){
fPt = fTrack->Pt();
virtual ~AliFlowCommonHist();
//make fill methods here
- Bool_t FillControlHistograms(AliFlowEventSimple* fEvent);
+ Bool_t FillControlHistograms(AliFlowEventSimple* Event);
//make get methods here
Double_t GetEntriesInPtBin(Int_t fBin); //gets entries from fHistPtDiff
private:
- AliFlowEventSimple* fEvent; //!
+ // AliFlowEventSimple* fEvent; //!
AliFlowTrackSimple* fTrack; //!
//define histograms here
//normal loop
while (fGoodTracks < fN && ftrkN < fNumberOfInputTracks) {
- AliMCParticle* fParticle = fInput->GetTrack(ftrkN); //get input particle
+ AliMCParticle* myParticle = fInput->GetTrack(ftrkN); //get input particle
//cut on tracks
- if (TMath::Abs(fParticle->Eta()) < 0.2)
+ if (TMath::Abs(myParticle->Eta()) < 0.2)
{
if(
- TMath::Abs(fParticle->Particle()->GetPdgCode()) == 211
- // TMath::Abs(fParticle->Particle()->GetPdgCode()) == 211 ||
- // TMath::Abs(fParticle->Particle()->GetPdgCode()) == 321 ||
- // TMath::Abs(fParticle->Particle()->GetPdgCode()) == 2212
+ TMath::Abs(myParticle->Particle()->GetPdgCode()) == 211
+ // TMath::Abs(myParticle->Particle()->GetPdgCode()) == 211 ||
+ // TMath::Abs(myParticle->Particle()->GetPdgCode()) == 321 ||
+ // TMath::Abs(myParticle->Particle()->GetPdgCode()) == 2212
)
{
fTrack = new AliFlowTrackSimple() ;
- fTrack->SetPt(fParticle->Pt() );
- fTrack->SetEta(fParticle->Eta() );
- fTrack->SetPhi(fParticle->Phi() );
+ fTrack->SetPt(myParticle->Pt() );
+ fTrack->SetEta(myParticle->Eta() );
+ fTrack->SetPhi(myParticle->Phi() );
fTrack->SetForIntegratedFlow(kTRUE);
fTrack->SetForDifferentialFlow(kTRUE);
fEvent->TrackCollection()->Add(fTrack) ;
}
/* else if(
- TMath::Abs(fParticle->Particle()->GetPdgCode()) == 211
+ TMath::Abs(myParticle->Particle()->GetPdgCode()) == 211
)
{
fTrack = new AliFlowTrackSimple();
- fTrack->SetPt(fParticle->Pt() );
- fTrack->SetEta(fParticle->Eta() );
- fTrack->SetPhi(fParticle->Phi() );
+ fTrack->SetPt(myParticle->Pt() );
+ fTrack->SetEta(myParticle->Eta() );
+ fTrack->SetPhi(myParticle->Phi() );
fTrack->SetForIntegratedFlow(kFALSE);
fTrack->SetForDifferentialFlow(kTRUE);
//normal loop
while (fGoodTracks < fN && ftrkN < fNumberOfInputTracks) {
- AliESDtrack* fParticle = fInput->GetTrack(ftrkN); //get input particle
+ AliESDtrack* myParticle = fInput->GetTrack(ftrkN); //get input particle
//cut on tracks
- if (TMath::Abs(fParticle->Eta()) < 0.2)
+ if (TMath::Abs(myParticle->Eta()) < 0.2)
{
fTrack = new AliFlowTrackSimple() ;
- fTrack->SetPt(fParticle->Pt() );
- fTrack->SetEta(fParticle->Eta() );
- fTrack->SetPhi(fParticle->Phi() );
+ fTrack->SetPt(myParticle->Pt() );
+ fTrack->SetEta(myParticle->Eta() );
+ fTrack->SetPhi(myParticle->Phi() );
fTrack->SetForIntegratedFlow(kTRUE);
fTrack->SetForDifferentialFlow(kTRUE);
//normal loop
while (fGoodTracks < fN && ftrkN < fNumberOfInputTracks) {
- AliAODTrack* fParticle = fInput->GetTrack(ftrkN); //get input particle
+ AliAODTrack* myParticle = fInput->GetTrack(ftrkN); //get input particle
//cut on tracks
- if (TMath::Abs(fParticle->Eta()) < 0.2)
+ if (TMath::Abs(myParticle->Eta()) < 0.2)
{
fTrack = new AliFlowTrackSimple() ;
- fTrack->SetPt(fParticle->Pt() );
- fTrack->SetEta(fParticle->Eta() );
- fTrack->SetPhi(fParticle->Phi() );
+ fTrack->SetPt(myParticle->Pt() );
+ fTrack->SetEta(myParticle->Eta() );
+ fTrack->SetPhi(myParticle->Phi() );
fTrack->SetForIntegratedFlow(kTRUE);
fTrack->SetForDifferentialFlow(kTRUE);
//normal loop
while (fGoodTracks < fN && ftrkN < fNumberOfInputTracks) {
- AliESDtrack* fParticle = fInput->GetTrack(ftrkN); //get input particle
+ AliESDtrack* myParticle = fInput->GetTrack(ftrkN); //get input particle
//get Label
- Int_t fLabel = fParticle->GetLabel();
+ Int_t fLabel = myParticle->GetLabel();
//match to mc particle
AliMCParticle* fMcParticle = fInputMc->GetTrack(TMath::Abs(fLabel));
//check
- if (TMath::Abs(fParticle->GetLabel())!=fMcParticle->Label()) cout<<"fParticle->GetLabel()!=fMcParticle->Label() "<<fParticle->GetLabel()<<" "<<fMcParticle->Label()<<endl;
+ if (TMath::Abs(myParticle->GetLabel())!=fMcParticle->Label()) cout<<"myParticle->GetLabel()!=fMcParticle->Label() "<<myParticle->GetLabel()<<" "<<fMcParticle->Label()<<endl;
//cut on tracks
- if (TMath::Abs(fParticle->Eta()) < 0.2)
+ if (TMath::Abs(myParticle->Eta()) < 0.2)
{
if(
TMath::Abs(fMcParticle->Particle()->GetPdgCode()) == 211 //pions
{
if(fOption == 0) { //take the PID from the MC & the kinematics from the ESD
fTrack = new AliFlowTrackSimple() ;
- fTrack->SetPt(fParticle->Pt() );
- fTrack->SetEta(fParticle->Eta() );
- fTrack->SetPhi(fParticle->Phi() );
+ fTrack->SetPt(myParticle->Pt() );
+ fTrack->SetEta(myParticle->Eta() );
+ fTrack->SetPhi(myParticle->Phi() );
fTrack->SetForIntegratedFlow(kTRUE);
fTrack->SetForDifferentialFlow(kTRUE);
}
-/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+/* copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
* See cxx source for full Copyright notice */
/* $Id$ */