#include "TH2D.h"
#include "AliAnaConvCorrPion.h"
//#include "AliAODTrack.h"
-#include "TClonesArray.h"
#include "AliAODConversionParticle.h"
-#include <AliLog.h>
-//#include "AliAODConversionMother.h"
-//#include "AliAODConversionPhoton.h"
-//#include "THnSparse.h"
-//#include "TH2F.h"
-#include <iostream>
using namespace std;
AliAnaConvCorrPion::AliAnaConvCorrPion() :
AliAnaConvCorrBase("pion_hadron_corr", "Pion dPhi"),
//hTriggerPtvsMass(NULL),
- fAxisM()
+ hTriggerPtvsMass(NULL),
+ fAxisM()
{
//consctructor
InitMassAxis();
AliAnaConvCorrPion::AliAnaConvCorrPion(TString name, TString title = "Pion Corr") :
AliAnaConvCorrBase(name, title),
//hTriggerPtvsMass(NULL),
+ hTriggerPtvsMass(NULL),
fAxisM()
{
//consctructor
Double_t mbins[7] = {0.1, 0.11, 0.12, 0.15, 0.16, 0.18, 0.2};
fAxisM.Set(6, mbins);
fAxisM.SetNameTitle("InvMass", "invariant mass");
- GetAxisList().AddAt(&fAxisM, 4);
- GetTrackAxisList().AddAt(&fAxisM, 3);
- GetTrigAxisList().AddAt(&fAxisM, 2);
+ GetAxisList().AddAt(&fAxisM, 6);
+ GetTrackAxisList().AddAt(&fAxisM, 5);
+ GetTrigAxisList().AddAt(&fAxisM, 4);
}
///________________________________________________________________________________
//Create histograms
CreateBaseHistograms();
- hTriggerPtvsMass[0] = new TH2D(Form("hTriggerPtvsMass_all_%s", GetName()), "Pt vs Mass all pizero", 400, 0, .400, GetAxistPt().GetNbins(), GetAxistPt().GetXbins()->GetArray());
- hTriggerPtvsMass[1] = new TH2D("hTriggerPtvsMass_leadingcone", "Pt vs Mass leading cone", 1, 0, .400, 1, 0, 100);
- hTriggerPtvsMass[2] = new TH2D("hTriggerPtvsMass_leadingevent", "Pt vs Mass leading event", 1, 0, .400, 1, 0, 100);
- GetHistograms()->Add(hTriggerPtvsMass[0]);
- //GetHistograms()->Add(hTriggerPtvsMass[1]);
- //GetHistograms()->Add(hTriggerPtvsMass[2]);
+ hTriggerPtvsMass = new TH2D(Form("hTriggerPtvsMass_all_%s", GetName()), "Pt vs Mass all pizero",
+ 400, 0, .400, GetAxistPt().GetNbins(), GetAxistPt().GetXbins()->GetArray());
+ GetHistograms()->Add(hTriggerPtvsMass);
}
///________________________________________________________________________________
-void AliAnaConvCorrPion::FillTriggerCounters(const AliAODConversionParticle * particle, Int_t leading) {
- //Fill histograms counting triggers
- //fHNTriggers[leading]->Fill(particle->Pt());
- AliDebug(AliLog::kDebug + 5, Form("Fill trigger countder %f %f", particle->M(), particle->Pt()));
- hTriggerPtvsMass[leading]->Fill(particle->M(), particle->Pt());
+void AliAnaConvCorrPion::FillTriggerCounters(const AliAODConversionParticle * particle) {
+ hTriggerPtvsMass->Fill(particle->M(), particle->Pt());
}
-//________________________________________________________________________________
-// void AliAnaConvCorrPion::Process(TClonesArray * pions, TClonesArray * photons, TClonesArray * tracks) {
-
-// for(Int_t ip = 0; ip < pions->GetEntriesFast(); ip++) {
-
-// AliAODConversionParticle * pion = static_cast<AliAODConversionParticle*>(pions->UncheckedAt(ip));
-
-// Int_t tIDs[4] = {-1, -1, -1, -1};
-// AliAODConversionParticle * photon1 = static_cast<AliAODConversionParticle*>(photons->UncheckedAt(pion->GetLabel(0)));
-// tIDs[0] = photon1->GetLabel(0);
-// tIDs[1] = photon1->GetLabel(1);
-// AliAODConversionParticle * photon2 = static_cast<AliAODConversionParticle*>(photons->UncheckedAt(pion->GetLabel(1)));
-// tIDs[2] = photon2->GetLabel(0);
-// tIDs[3] = photon2->GetLabel(1);
-
-// CorrelateWithTracks(static_cast<AliAODConversionParticle*>(pion), tracks, tIDs, kFALSE);
-// }
-// }
-