// --- ROOT system ---
#include "TParticle.h"
#include "TH2F.h"
+#include "TDatabasePDG.h"
//---- AliRoot system ----
#include "AliAnaChargedParticles.h"
#include "AliFiducialCut.h"
#include "AliVTrack.h"
#include "AliAODMCParticle.h"
+#include "AliAODTrack.h"
+#include "AliAODEvent.h"
+#include "AliESDEvent.h"
ClassImp(AliAnaChargedParticles)
//__________________________________________________
- AliAnaChargedParticles::AliAnaChargedParticles() :
+ AliAnaChargedParticles::AliAnaChargedParticles() :
AliAnaCaloTrackCorrBaseClass(),
- fPdg(0),
- fhNtracks(0), fhPt(0),
- fhPhiNeg(0), fhEtaNeg(0),
- fhPhiPos(0), fhEtaPos(0),
- fhEtaPhiPos(0), fhEtaPhiNeg(0),
- //MC
- fhPtPion(0), fhPhiPion(0), fhEtaPion(0),
- fhPtProton(0), fhPhiProton(0), fhEtaProton(0),
- fhPtElectron(0),fhPhiElectron(0), fhEtaElectron(0),
- fhPtKaon(0), fhPhiKaon(0), fhEtaKaon(0),
- fhPtUnknown(0), fhPhiUnknown(0), fhEtaUnknown(0)
+ fFillPileUpHistograms(0), fFillTrackBCHistograms(0),
+ fFillVertexBC0Histograms(0),
+ //Histograms
+ fhNtracks(0), fhPt(0), fhPtNoCut(0),
+ fhPtCutDCA(0), fhPtCutDCABCOK(0),
+ fhPhiNeg(0), fhEtaNeg(0),
+ fhPhiPos(0), fhEtaPos(0),
+ fhEtaPhiPos(0), fhEtaPhiNeg(0),
+ fhPtVtxOutBC0(0), fhEtaPhiVtxOutBC0(0),
+ fhPtVtxInBC0(0), fhEtaPhiVtxInBC0(0),
+ fhPtSPDRefit(0), fhPtNoSPDRefit(0), fhPtNoSPDNoRefit(0),
+ fhEtaPhiSPDRefitPt02(0), fhEtaPhiNoSPDRefitPt02(0), fhEtaPhiNoSPDNoRefitPt02(0),
+ fhEtaPhiSPDRefitPt3(0), fhEtaPhiNoSPDRefitPt3(0), fhEtaPhiNoSPDNoRefitPt3(0),
+
+ //TOF
+ fhTOFSignal(0), fhTOFSignalPtCut(0), fhTOFSignalBCOK(0),
+ fhPtTOFSignal(0), fhPtTOFSignalDCACut(0),
+ fhPtTOFSignalVtxOutBC0(0), fhPtTOFSignalVtxInBC0(0),
+ fhPtTOFStatus0(0), fhEtaPhiTOFStatus0(0),
+ fhEtaPhiTOFBC0(0), fhEtaPhiTOFBCPlus(0), fhEtaPhiTOFBCMinus(0),
+ fhEtaPhiTOFBC0PileUpSPD(0),
+ fhEtaPhiTOFBCPlusPileUpSPD(0),
+ fhEtaPhiTOFBCMinusPileUpSPD(0),
+ fhProductionVertexBC(0),
+ fhPtNPileUpSPDVtx(0), fhPtNPileUpTrkVtx(0),
+ fhPtNPileUpSPDVtxBC0(0), fhPtNPileUpTrkVtxBC0(0)
+
{
//Default Ctor
+ for(Int_t i = 0; i < 7; i++)
+ {
+ fhPtPileUp [i] = 0;
+ fhPtTOFSignalPileUp[i] = 0;
+ fhPtTOFSignalVtxOutBC0PileUp[i] = 0;
+ fhPtTOFSignalVtxInBC0PileUp [i] = 0;
+ fhProductionVertexBCPileUp [i] = 0;
+ }
+
+ for(Int_t i = 0; i < 3; i++)
+ {
+ fhPtDCA [i] = 0 ;
+
+ fhPtDCASPDRefit [i] = 0 ;
+ fhPtDCANoSPDRefit [i] = 0 ;
+ fhPtDCANoSPDNoRefit [i] = 0 ;
+
+ fhPtDCAPileUp [i] = 0 ;
+ fhPtDCATOFBC0 [i] = 0 ;
+ fhPtDCATOFBCOut [i] = 0 ;
+ fhPtDCAPileUpTOFBC0 [i] = 0 ;
+ fhPtDCANoTOFHit [i] = 0 ;
+ fhPtDCAPileUpNoTOFHit [i] = 0 ;
+
+ fhPtDCAVtxOutBC0 [i] = 0 ;
+ fhPtDCAVtxInBC0 [i] = 0 ;
+ fhPtDCAVtxOutBC0PileUp[i] = 0 ;
+ fhPtDCAVtxInBC0PileUp [i] = 0 ;
+ fhPtDCAVtxOutBC0NoTOFHit[i] = 0 ;
+ fhPtDCAVtxInBC0NoTOFHit [i] = 0 ;
+ fhPtDCAVtxOutBC0PileUpNoTOFHit[i] = 0 ;
+ fhPtDCAVtxInBC0PileUpNoTOFHit [i] = 0 ;
+ }
+
+ // MC
+ for(Int_t imcPart = 0; imcPart < 6; imcPart++)
+ {
+ fhPtMCPart [imcPart] = 0;
+ fhPtMCPrimPart [imcPart] = 0;
+ fhPhiMCPart [imcPart] = 0;
+ fhPhiMCPrimPart[imcPart] = 0;
+ fhEtaMCPart [imcPart] = 0;
+ fhEtaMCPrimPart[imcPart] = 0;
+ }
+
//Initialize parameters
InitParameters();
}
+//__________________________________________________
+void AliAnaChargedParticles::FillPrimaryHistograms()
+{
+ // Fill primary generated particles histograms if MC data is available
+
+ Int_t pdg = 0 ;
+ Int_t nprim = 0 ;
+
+ TParticle * primStack = 0;
+ AliAODMCParticle * primAOD = 0;
+ TLorentzVector lv;
+
+ // Get the ESD MC particles container
+ AliStack * stack = 0;
+ if( GetReader()->ReadStack() )
+ {
+ stack = GetMCStack();
+ if(!stack ) return;
+ nprim = stack->GetNtrack();
+ }
+
+ // Get the AOD MC particles container
+ TClonesArray * mcparticles = 0;
+ if( GetReader()->ReadAODMCParticles() )
+ {
+ mcparticles = GetReader()->GetAODMCParticles();
+ if( !mcparticles ) return;
+ nprim = mcparticles->GetEntriesFast();
+ }
+
+ for(Int_t i=0 ; i < nprim; i++)
+ {
+ if(GetReader()->AcceptOnlyHIJINGLabels() && !GetReader()->IsHIJINGLabel(i)) continue ;
+
+ if(GetReader()->ReadStack())
+ {
+ primStack = stack->Particle(i) ;
+ if(!primStack)
+ {
+ printf("AliAnaChargedParticles::FillPrimaryMCHistograms() - ESD primaries pointer not available!!\n");
+ continue;
+ }
+
+ if( primStack->GetStatusCode() != 1 ) continue;
+
+ Int_t charge = (Int_t )TDatabasePDG::Instance()->GetParticle(primStack->GetPdgCode())->Charge();
+ if( TMath::Abs(charge) == 0 ) continue;
+
+ pdg = TMath::Abs(primStack->GetPdgCode());
+
+ if(primStack->Energy() == TMath::Abs(primStack->Pz())) continue ; //Protection against floating point exception
+
+ //printf("i %d, %s %d %s %d \n",i, stack->Particle(i)->GetName(), stack->Particle(i)->GetPdgCode(),
+ // prim->GetName(), prim->GetPdgCode());
+
+ //Charged kinematics
+ primStack->Momentum(lv);
+ }
+ else
+ {
+ primAOD = (AliAODMCParticle *) mcparticles->At(i);
+ if(!primAOD)
+ {
+ printf("AliAnaChargedParticles::FillPrimaryHistograms() - AOD primaries pointer not available!!\n");
+ continue;
+ }
+
+ if( primAOD->GetStatus() != 1 ) continue;
+
+ if(TMath::Abs(primAOD->Charge()) == 0 ) continue;
+
+ pdg = TMath::Abs(primAOD->GetPdgCode());
+
+ if(primAOD->E() == TMath::Abs(primAOD->Pz())) continue ; //Protection against floating point exception
+
+ //Charged kinematics
+ lv.SetPxPyPzE(primAOD->Px(),primAOD->Py(),primAOD->Pz(),primAOD->E());
+ }
+
+ Int_t mcType = kmcUnknown;
+ if (pdg==211 ) mcType = kmcPion;
+ else if(pdg==2212) mcType = kmcProton;
+ else if(pdg==321 ) mcType = kmcKaon;
+ else if(pdg==11 ) mcType = kmcElectron;
+ else if(pdg==13 ) mcType = kmcMuon;
+
+ //printf("pdg %d, charge %f, mctype %d\n",pdg, charge, mcType);
+
+ Bool_t in = GetFiducialCut()->IsInFiducialCut(lv,"CTS") ;
+
+ Float_t ptPrim = lv.Pt();
+ Float_t etaPrim = lv.Eta();
+ Float_t phiPrim = lv.Phi();
+ if(phiPrim < 0) phiPrim+=TMath::TwoPi();
+
+ if(in)
+ fhPtMCPrimPart[mcType]->Fill(ptPrim);
+ fhEtaMCPrimPart[mcType]->Fill(ptPrim,etaPrim);
+ fhPhiMCPrimPart[mcType]->Fill(ptPrim,phiPrim);
+ }
+}
+
//_______________________________________________________
TList * AliAnaChargedParticles::GetCreateOutputObjects()
{
TList * outputContainer = new TList() ;
- outputContainer->SetName("ExampleHistos") ;
+ outputContainer->SetName("ChargedParticleHistos") ;
Int_t nptbins = GetHistogramRanges()->GetHistoPtBins(); Int_t nphibins = GetHistogramRanges()->GetHistoPhiBins(); Int_t netabins = GetHistogramRanges()->GetHistoEtaBins();
Float_t ptmax = GetHistogramRanges()->GetHistoPtMax(); Float_t phimax = GetHistogramRanges()->GetHistoPhiMax(); Float_t etamax = GetHistogramRanges()->GetHistoEtaMax();
fhNtracks->SetXTitle("# of tracks");
outputContainer->Add(fhNtracks);
- fhPt = new TH1F ("hPt","p_T distribution", nptbins,ptmin,ptmax);
- fhPt->SetXTitle("p_{T} (GeV/c)");
+ fhPt = new TH1F ("hPt","#it{p}_{T} distribution", nptbins,ptmin,ptmax);
+ fhPt->SetXTitle("#it{p}_{T} (GeV/#it{c})");
outputContainer->Add(fhPt);
+ fhPtNoCut = new TH1F ("hPtNoCut","#it{p}_{T} distribution, raw tracks", nptbins,ptmin,ptmax);
+ fhPtNoCut->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhPtNoCut);
+
+ if(!GetReader()->IsDCACutOn())
+ {
+ fhPtCutDCA = new TH1F ("hPtCutDCA","#it{p}_{T} distribution, cut DCA", nptbins,ptmin,ptmax);
+ fhPtCutDCA->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhPtCutDCA);
+ }
+
+ if(fFillVertexBC0Histograms && !GetReader()->IsDCACutOn())
+ {
+ fhPtCutDCABCOK = new TH1F ("hPtCutDCABCOK","#it{p}_{T} distribution, DCA cut, track BC=0 or -100", nptbins,ptmin,ptmax);
+ fhPtCutDCABCOK->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhPtCutDCABCOK);
+ }
+
fhPhiNeg = new TH2F ("hPhiNegative","#phi of negative charges distribution",
nptbins,ptmin,ptmax, nphibins,phimin,phimax);
fhPhiNeg->SetYTitle("#phi (rad)");
- fhPhiNeg->SetXTitle("p_{T} (GeV/c)");
+ fhPhiNeg->SetXTitle("#it{p}_{T} (GeV/#it{c})");
outputContainer->Add(fhPhiNeg);
fhEtaNeg = new TH2F ("hEtaNegative","#eta of negative charges distribution",
nptbins,ptmin,ptmax, netabins,etamin,etamax);
fhEtaNeg->SetYTitle("#eta ");
- fhEtaNeg->SetXTitle("p_{T} (GeV/c)");
+ fhEtaNeg->SetXTitle("#it{p}_{T} (GeV/#it{c})");
outputContainer->Add(fhEtaNeg);
fhPhiPos = new TH2F ("hPhiPositive","#phi of positive charges distribution",
nptbins,ptmin,ptmax, nphibins,phimin,phimax);
fhPhiPos->SetYTitle("#phi (rad)");
- fhPhiPos->SetXTitle("p_{T} (GeV/c)");
+ fhPhiPos->SetXTitle("#it{p}_{T} (GeV/#it{c})");
outputContainer->Add(fhPhiPos);
fhEtaPos = new TH2F ("hEtaPositive","#eta of positive charges distribution",
nptbins,ptmin,ptmax, netabins,etamin,etamax);
fhEtaPos->SetYTitle("#eta ");
- fhEtaPos->SetXTitle("p_{T} (GeV/c)");
+ fhEtaPos->SetXTitle("#it{p}_{T} (GeV/#it{c})");
outputContainer->Add(fhEtaPos);
- fhEtaPhiPos = new TH2F ("hEtaPhiPositive","pt/eta/phi of positive charge",netabins,etamin,etamax, nphibins,phimin,phimax);
+ fhEtaPhiPos = new TH2F ("hEtaPhiPositive","pt/eta/phi of positive charge",netabins,etamin,etamax, nphibins,phimin,phimax);
fhEtaPhiPos->SetXTitle("#eta ");
fhEtaPhiPos->SetYTitle("#phi (rad)");
outputContainer->Add(fhEtaPhiPos);
- fhEtaPhiNeg = new TH2F ("hEtaPhiNegative","eta vs phi of negative charge",netabins,etamin,etamax, nphibins,phimin,phimax);
+ fhEtaPhiNeg = new TH2F ("hEtaPhiNegative","#eta vs #phi of negative charge",netabins,etamin,etamax, nphibins,phimin,phimax);
fhEtaPhiNeg->SetXTitle("#eta ");
fhEtaPhiNeg->SetYTitle("#phi (rad)");
outputContainer->Add(fhEtaPhiNeg);
- if(IsDataMC()){
+ if(fFillVertexBC0Histograms)
+ {
+ fhPtVtxOutBC0 = new TH1F ("hPtVtxOutBC0","#it{p}_{T} distribution, vertex in BC=0", nptbins,ptmin,ptmax);
+ fhPtVtxOutBC0->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhPtVtxOutBC0);
- fhPtPion = new TH1F ("hPtMCPion","p_T distribution from #pi", nptbins,ptmin,ptmax);
- fhPtPion->SetXTitle("p_{T} (GeV/c)");
- outputContainer->Add(fhPtPion);
+ fhEtaPhiVtxOutBC0 = new TH2F ("hEtaPhiVtxOutBC0","#eta vs #phi of all charges with vertex in BC=0",netabins,etamin,etamax, nphibins,phimin,phimax);
+ fhEtaPhiVtxOutBC0->SetXTitle("#eta ");
+ fhEtaPhiVtxOutBC0->SetYTitle("#phi (rad)");
+ outputContainer->Add(fhEtaPhiVtxOutBC0);
- fhPhiPion = new TH2F ("hPhiMCPion","#phi distribution from #pi",nptbins,ptmin,ptmax, nphibins,phimin,phimax);
- fhPhiPion->SetXTitle("#phi (rad)");
- outputContainer->Add(fhPhiPion);
+ fhPtVtxInBC0 = new TH1F ("hPtVtxInBC0","#it{p}_{T} distribution, vertex in BC=0", nptbins,ptmin,ptmax);
+ fhPtVtxInBC0->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhPtVtxInBC0);
- fhEtaPion = new TH2F ("hEtaMCPion","#eta distribution from #pi",nptbins,ptmin,ptmax, netabins,etamin,etamax);
- fhEtaPion->SetXTitle("#eta ");
- outputContainer->Add(fhEtaPion);
-
- fhPtProton = new TH1F ("hPtMCProton","p_T distribution from proton", nptbins,ptmin,ptmax);
- fhPtProton->SetXTitle("p_{T} (GeV/c)");
- outputContainer->Add(fhPtProton);
-
- fhPhiProton = new TH2F ("hPhiMCProton","#phi distribution from proton",nptbins,ptmin,ptmax, nphibins,phimin,phimax);
- fhPhiProton->SetXTitle("#phi (rad)");
- outputContainer->Add(fhPhiProton);
+ fhEtaPhiVtxInBC0 = new TH2F ("hEtaPhiVtxInBC0","#eta vs #phi of all charges with vertex in BC=0",netabins,etamin,etamax, nphibins,phimin,phimax);
+ fhEtaPhiVtxInBC0->SetXTitle("#eta ");
+ fhEtaPhiVtxInBC0->SetYTitle("#phi (rad)");
+ outputContainer->Add(fhEtaPhiVtxInBC0);
+ }
+
+ fhPtSPDRefit = new TH1F ("hPtSPDRefit","#it{p}_{T} distribution of tracks with SPD and ITS refit", nptbins,ptmin,ptmax);
+ fhPtSPDRefit->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhPtSPDRefit);
+
+ fhEtaPhiSPDRefitPt02 = new TH2F ("hEtaPhiSPDRefitPt02","#eta vs #phi of tracks with SPD and ITS refit, #it{p}_{T}< 2 GeV/#it{c}",
+ netabins,etamin,etamax, nphibins,phimin,phimax);
+ fhEtaPhiSPDRefitPt02->SetXTitle("#eta ");
+ fhEtaPhiSPDRefitPt02->SetYTitle("#phi (rad)");
+ outputContainer->Add(fhEtaPhiSPDRefitPt02);
+
+ fhEtaPhiSPDRefitPt3 = new TH2F ("hEtaPhiSPDRefitPt3","#eta vs #phi of tracks with SPD and ITS refit, #it{p}_{T}> 3 GeV/#it{c}",
+ netabins,etamin,etamax, nphibins,phimin,phimax);
+ fhEtaPhiSPDRefitPt3->SetXTitle("#eta ");
+ fhEtaPhiSPDRefitPt3->SetYTitle("#phi (rad)");
+ outputContainer->Add(fhEtaPhiSPDRefitPt3);
+
+ fhPtNoSPDRefit = new TH1F ("hPtNoSPDRefit","#it{p}_{T} distribution of constrained tracks no SPD and with ITSRefit",
+ nptbins,ptmin,ptmax);
+ fhPtNoSPDRefit->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhPtNoSPDRefit);
+
+ fhEtaPhiNoSPDRefitPt02 = new TH2F ("hEtaPhiNoSPDRefitPt02","#eta vs #phi of constrained tracks no SPD and with ITSRefit, #it{p}_{T}< 2 GeV/#it{c}",
+ netabins,etamin,etamax, nphibins,phimin,phimax);
+ fhEtaPhiNoSPDRefitPt02->SetXTitle("#eta ");
+ fhEtaPhiNoSPDRefitPt02->SetYTitle("#phi (rad)");
+ outputContainer->Add(fhEtaPhiNoSPDRefitPt02);
+
+ fhEtaPhiNoSPDRefitPt3 = new TH2F ("hEtaPhiNoSPDRefitPt3","#eta vs #phi of of constrained tracks no SPD and with ITSRefit, #it{p}_{T}> 3 GeV/#it{c}",
+ netabins,etamin,etamax, nphibins,phimin,phimax);
+ fhEtaPhiNoSPDRefitPt3->SetXTitle("#eta ");
+ fhEtaPhiNoSPDRefitPt3->SetYTitle("#phi (rad)");
+ outputContainer->Add(fhEtaPhiNoSPDRefitPt3);
+
+ fhPtNoSPDNoRefit = new TH1F ("hPtNoSPDNoRefit","#it{p}_{T} distribution of constrained tracks with no SPD requierement and without ITSRefit",
+ nptbins,ptmin,ptmax);
+ fhPtNoSPDNoRefit->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhPtNoSPDNoRefit);
+
+ fhEtaPhiNoSPDNoRefitPt02 = new TH2F ("hEtaPhiNoSPDNoRefitPt02",
+ "#eta vs #phi of constrained tracks with no SPD requierement and without ITSRefit, #it{p}_{T}< 2 GeV/#it{c}",
+ netabins,etamin,etamax, nphibins,phimin,phimax);
+ fhEtaPhiNoSPDNoRefitPt02->SetXTitle("#eta ");
+ fhEtaPhiNoSPDNoRefitPt02->SetYTitle("#phi (rad)");
+ outputContainer->Add(fhEtaPhiNoSPDNoRefitPt02);
+
+ fhEtaPhiNoSPDNoRefitPt3 = new TH2F ("hEtaPhiNoSPDNoRefitPt3",
+ "#eta vs #phi of constrained tracks with no SPD requierement and without ITSRefit, #it{p}_{T}> 3 GeV/#it{c}",
+ netabins,etamin,etamax, nphibins,phimin,phimax);
+ fhEtaPhiNoSPDNoRefitPt3->SetXTitle("#eta ");
+ fhEtaPhiNoSPDNoRefitPt3->SetYTitle("#phi (rad)");
+ outputContainer->Add(fhEtaPhiNoSPDNoRefitPt3);
+
+ if(fFillVertexBC0Histograms)
+ {
+ fhProductionVertexBC = new TH1F("hProductionVertexBC", "tracks production vertex bunch crossing ", 41 , -20 , 20 ) ;
+ fhProductionVertexBC->SetYTitle("# tracks");
+ fhProductionVertexBC->SetXTitle("Bunch crossing");
+ outputContainer->Add(fhProductionVertexBC);
+ }
+
+ Int_t ntofbins = 1000;
+ Int_t mintof = -500;
+ Int_t maxtof = 500;
+
+ fhTOFSignal = new TH1F ("hTOFSignal","TOF signal", ntofbins,mintof,maxtof);
+ fhTOFSignal->SetXTitle("TOF signal (ns)");
+ outputContainer->Add(fhTOFSignal);
+
+ if(fFillTrackBCHistograms)
+ {
+ fhTOFSignalBCOK = new TH1F ("hTOFSignalBCOK","TOF signal", ntofbins,mintof,maxtof);
+ fhTOFSignalBCOK->SetXTitle("TOF signal (ns)");
+ outputContainer->Add(fhTOFSignalBCOK);
+ }
+
+ fhTOFSignalPtCut = new TH1F ("hTOFSignalPtCut","TOF signal", ntofbins,mintof,maxtof);
+ fhTOFSignalPtCut->SetXTitle("TOF signal (ns)");
+ outputContainer->Add(fhTOFSignalPtCut);
+
+ fhPtTOFSignal = new TH2F ("hPtTOFSignal","TOF signal", nptbins,ptmin,ptmax,ntofbins,mintof,maxtof);
+ fhPtTOFSignal->SetYTitle("TOF signal (ns)");
+ fhPtTOFSignal->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhPtTOFSignal);
+
+ if(!GetReader()->IsDCACutOn())
+ {
+ fhPtTOFSignalDCACut = new TH2F ("hPtTOFSignalDCACut","TOF signal after DCA cut", nptbins,ptmin,ptmax,ntofbins,mintof,maxtof);
+ fhPtTOFSignalDCACut->SetYTitle("TOF signal (ns)");
+ fhPtTOFSignalDCACut->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhPtTOFSignalDCACut);
+ }
+
+ if(fFillVertexBC0Histograms)
+ {
+ fhPtTOFSignalVtxOutBC0 = new TH2F ("hPtTOFSignalVtxOutBC0","TOF signal, vtx BC!=0", nptbins,ptmin,ptmax,ntofbins,mintof,maxtof);
+ fhPtTOFSignalVtxOutBC0->SetYTitle("TOF signal (ns)");
+ fhPtTOFSignalVtxOutBC0->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhPtTOFSignalVtxOutBC0);
- fhEtaProton = new TH2F ("hEtaMCProton","#eta distribution from proton",nptbins,ptmin,ptmax, netabins,etamin,etamax);
- fhEtaProton->SetXTitle("#eta ");
- outputContainer->Add(fhEtaProton);
+ fhPtTOFSignalVtxInBC0 = new TH2F ("hPtTOFSignalVtxInBC0","TOF signal, vtx BC=0", nptbins,ptmin,ptmax,ntofbins,mintof,maxtof);
+ fhPtTOFSignalVtxInBC0->SetYTitle("TOF signal (ns)");
+ fhPtTOFSignalVtxInBC0->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhPtTOFSignalVtxInBC0);
+ }
+
+ if(fFillPileUpHistograms)
+ {
+ TString pileUpName[] = {"SPD","EMCAL","SPDOrEMCAL","SPDAndEMCAL","SPDAndNotEMCAL","EMCALAndNotSPD","NotSPDAndNotEMCAL"} ;
- fhPtKaon = new TH1F ("hPtMCKaon","p_T distribution from kaon", nptbins,ptmin,ptmax);
- fhPtKaon->SetXTitle("p_{T} (GeV/c)");
- outputContainer->Add(fhPtKaon);
+ for(Int_t i = 0 ; i < 7 ; i++)
+ {
+ fhPtPileUp[i] = new TH1F(Form("hPtPileUp%s",pileUpName[i].Data()),
+ Form("Track #it{p}_{T}distribution, %s Pile-Up event",pileUpName[i].Data()),
+ nptbins,ptmin,ptmax);
+ fhPtPileUp[i]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhPtPileUp[i]);
+
+ fhPtTOFSignalPileUp[i] = new TH2F(Form("hPtTOFSignalPileUp%s",pileUpName[i].Data()),
+ Form("Track TOF vs #it{p}_{T}distribution, %s Pile-Up event",pileUpName[i].Data()),
+ nptbins,ptmin,ptmax,ntofbins,mintof,maxtof);
+ fhPtTOFSignalPileUp[i]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ fhPtTOFSignalPileUp[i]->SetXTitle("TOF signal (ns)");
+ outputContainer->Add(fhPtTOFSignalPileUp[i]);
+
+ if(fFillVertexBC0Histograms)
+ {
+ fhPtTOFSignalVtxOutBC0PileUp[i] = new TH2F(Form("hPtTOFSignalVtxOutBC0PileUp%s",pileUpName[i].Data()),
+ Form("Track TOF vs #it{p}_{T}distribution, %s Pile-Up event, vtx BC!=0",pileUpName[i].Data()),
+ nptbins,ptmin,ptmax,ntofbins,mintof,maxtof);
+ fhPtTOFSignalVtxOutBC0PileUp[i]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ fhPtTOFSignalVtxOutBC0PileUp[i]->SetXTitle("TOF signal (ns)");
+ outputContainer->Add(fhPtTOFSignalVtxOutBC0PileUp[i]);
+
+ fhPtTOFSignalVtxInBC0PileUp[i] = new TH2F(Form("hPtTOFSignalVtxInBC0PileUp%s",pileUpName[i].Data()),
+ Form("Track TOF vs #it{p}_{T}distribution, %s Pile-Up event, vtx BC=0",pileUpName[i].Data()),
+ nptbins,ptmin,ptmax,ntofbins,mintof,maxtof);
+ fhPtTOFSignalVtxInBC0PileUp[i]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ fhPtTOFSignalVtxInBC0PileUp[i]->SetXTitle("TOF signal (ns)");
+ outputContainer->Add(fhPtTOFSignalVtxInBC0PileUp[i]);
+ }
+
+ if(fFillVertexBC0Histograms)
+ {
+ fhProductionVertexBCPileUp[i] = new TH1F(Form("hProductionVertexBCPileUp%s",pileUpName[i].Data()),
+ Form("tracks production vertex bunch crossing, %s Pile-Up event",pileUpName[i].Data()),
+ 41 , -20 , 20 ) ;
+ fhProductionVertexBCPileUp[i]->SetYTitle("# tracks");
+ fhProductionVertexBCPileUp[i]->SetXTitle("Bunch crossing");
+ outputContainer->Add(fhProductionVertexBCPileUp[i]);
+ }
+ }
+
- fhPhiKaon = new TH2F ("hPhiMCKaon","#phi distribution from kaon",nptbins,ptmin,ptmax, nphibins,phimin,phimax);
- fhPhiKaon->SetXTitle("#phi (rad)");
- outputContainer->Add(fhPhiKaon);
+ if(fFillTrackBCHistograms)
+ {
+ fhEtaPhiTOFBC0 = new TH2F ("hEtaPhiTOFBC0","eta-phi for tracks with hit on TOF, and tof corresponding to BC=0",netabins,etamin,etamax, nphibins,phimin,phimax);
+ fhEtaPhiTOFBC0->SetXTitle("#eta ");
+ fhEtaPhiTOFBC0->SetYTitle("#phi (rad)");
+ outputContainer->Add(fhEtaPhiTOFBC0);
+
+ fhEtaPhiTOFBCPlus = new TH2F ("hEtaPhiTOFBCPlus","eta-phi for tracks with hit on TOF, and tof corresponding to BC>0",netabins,etamin,etamax, nphibins,phimin,phimax);
+ fhEtaPhiTOFBCPlus->SetXTitle("#eta ");
+ fhEtaPhiTOFBCPlus->SetYTitle("#phi (rad)");
+ outputContainer->Add(fhEtaPhiTOFBCPlus);
+
+ fhEtaPhiTOFBCMinus = new TH2F ("hEtaPhiTOFBCMinus","eta-phi for tracks with hit on TOF, and tof corresponding to BC<0",netabins,etamin,etamax, nphibins,phimin,phimax);
+ fhEtaPhiTOFBCMinus->SetXTitle("#eta ");
+ fhEtaPhiTOFBCMinus->SetYTitle("#phi (rad)");
+ outputContainer->Add(fhEtaPhiTOFBCMinus);
+
+ if(fFillPileUpHistograms)
+ {
+ fhEtaPhiTOFBC0PileUpSPD = new TH2F ("hEtaPhiTOFBC0PileUpSPD","eta-phi for tracks with hit on TOF, and tof corresponding to BC=0, SPD pile-up",netabins,etamin,etamax, nphibins,phimin,phimax);
+ fhEtaPhiTOFBC0PileUpSPD->SetXTitle("#eta ");
+ fhEtaPhiTOFBC0PileUpSPD->SetYTitle("#phi (rad)");
+ outputContainer->Add(fhEtaPhiTOFBC0PileUpSPD);
+
+ fhEtaPhiTOFBCPlusPileUpSPD = new TH2F ("hEtaPhiTOFBCPlusPileUpSPD","eta-phi for tracks with hit on TOF, and tof corresponding to BC>0, SPD pile-up",netabins,etamin,etamax, nphibins,phimin,phimax);
+ fhEtaPhiTOFBCPlusPileUpSPD->SetXTitle("#eta ");
+ fhEtaPhiTOFBCPlusPileUpSPD->SetYTitle("#phi (rad)");
+ outputContainer->Add(fhEtaPhiTOFBCPlusPileUpSPD);
+
+ fhEtaPhiTOFBCMinusPileUpSPD = new TH2F ("hEtaPhiTOFBCMinusPileUpSPD","eta-phi for tracks with hit on TOF, and tof corresponding to BC<0, SPD pile-up",netabins,etamin,etamax, nphibins,phimin,phimax);
+ fhEtaPhiTOFBCMinusPileUpSPD->SetXTitle("#eta ");
+ fhEtaPhiTOFBCMinusPileUpSPD->SetYTitle("#phi (rad)");
+ outputContainer->Add(fhEtaPhiTOFBCMinusPileUpSPD);
+ }
+ }
- fhEtaKaon = new TH2F ("hEtaMCKaon","#eta distribution from kaon",nptbins,ptmin,ptmax, netabins,etamin,etamax);
- fhEtaKaon->SetXTitle("#eta ");
- outputContainer->Add(fhEtaKaon);
+ }
+
+ fhPtTOFStatus0 = new TH1F ("hPtTOFStatus0","#it{p}_{T} distribution of tracks not hitting TOF", nptbins,ptmin,ptmax);
+ fhPtTOFStatus0->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhPtTOFStatus0);
+
+
+ fhEtaPhiTOFStatus0 = new TH2F ("hEtaPhiTOFStatus0","eta-phi for tracks without hit on TOF",netabins,etamin,etamax, nphibins,phimin,phimax);
+ fhEtaPhiTOFStatus0->SetXTitle("#eta ");
+ fhEtaPhiTOFStatus0->SetYTitle("#phi (rad)");
+ outputContainer->Add(fhEtaPhiTOFStatus0);
+
+ TString dcaName[] = {"xy","z","Cons"} ;
+ Int_t ndcabins = 800;
+ Int_t mindca = -4;
+ Int_t maxdca = 4;
+
+ for(Int_t i = 0 ; i < 3 ; i++)
+ {
- fhPtElectron = new TH1F ("hPtMCElectron","p_T distribution from electron", nptbins,ptmin,ptmax);
- fhPtElectron->SetXTitle("p_{T} (GeV/c)");
- outputContainer->Add(fhPtElectron);
+ fhPtDCA[i] = new TH2F(Form("hPtDCA%s",dcaName[i].Data()),
+ Form("Track DCA%s vs #it{p}_{T}distribution",dcaName[i].Data()),
+ nptbins,ptmin,ptmax,ndcabins,mindca,maxdca);
+ fhPtDCA[i]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ fhPtDCA[i]->SetYTitle(Form("DCA_{%s}",dcaName[i].Data()));
+ outputContainer->Add(fhPtDCA[i]);
- fhPhiElectron = new TH2F ("hPhiMCElectron","#phi distribution from electron",nptbins,ptmin,ptmax, nphibins,phimin,phimax);
- fhPhiElectron->SetXTitle("#phi (rad)");
- outputContainer->Add(fhPhiElectron);
+ fhPtDCASPDRefit[i] = new TH2F(Form("hPtDCA%sSPDRefit",dcaName[i].Data()),
+ Form("Track DCA%s vs #it{p}_{T}distribution of tracks with SPD and ITS refit",dcaName[i].Data()),
+ nptbins,ptmin,ptmax,ndcabins,mindca,maxdca);
+ fhPtDCASPDRefit[i]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ fhPtDCASPDRefit[i]->SetYTitle(Form("DCA_{%s}",dcaName[i].Data()));
+ outputContainer->Add(fhPtDCASPDRefit[i]);
+
+ fhPtDCANoSPDRefit[i] = new TH2F(Form("hPtDCA%sNoSPDRefit",dcaName[i].Data()),
+ Form("Track DCA%s vs #it{p}_{T}distributionof constrained tracks no SPD and with ITSRefit",dcaName[i].Data()),
+ nptbins,ptmin,ptmax,ndcabins,mindca,maxdca);
+ fhPtDCANoSPDRefit[i]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ fhPtDCANoSPDRefit[i]->SetYTitle(Form("DCA_{%s}",dcaName[i].Data()));
+ outputContainer->Add(fhPtDCANoSPDRefit[i]);
+
+ fhPtDCANoSPDNoRefit[i] = new TH2F(Form("hPtDCA%sNoSPDNoRefit",dcaName[i].Data()),
+ Form("Track DCA%s vs #it{p}_{T}distribution, constrained tracks with no SPD requierement and without ITSRefit",dcaName[i].Data()),
+ nptbins,ptmin,ptmax,ndcabins,mindca,maxdca);
+ fhPtDCANoSPDNoRefit[i]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ fhPtDCANoSPDNoRefit[i]->SetYTitle(Form("DCA_{%s}",dcaName[i].Data()));
+ outputContainer->Add(fhPtDCANoSPDNoRefit[i]);
- fhEtaElectron = new TH2F ("hEtaMCElectron","#eta distribution from electron",nptbins,ptmin,ptmax, netabins,etamin,etamax);
- fhEtaElectron->SetXTitle("#eta ");
- outputContainer->Add(fhEtaElectron);
+ if(fFillTrackBCHistograms)
+ {
+ fhPtDCATOFBC0[i] = new TH2F(Form("hPtDCA%sTOFBC0",dcaName[i].Data()),
+ Form("Track DCA%s vs #it{p}_{T}distribution, BC=0",dcaName[i].Data()),
+ nptbins,ptmin,ptmax,ndcabins,mindca,maxdca);
+ fhPtDCATOFBC0[i]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ fhPtDCATOFBC0[i]->SetYTitle(Form("DCA_{%s}",dcaName[i].Data()));
+ outputContainer->Add(fhPtDCATOFBC0[i]);
+
+ fhPtDCATOFBCOut[i] = new TH2F(Form("hPtDCA%sTOFBCOut",dcaName[i].Data()),
+ Form("Track DCA%s vs #it{p}_{T}distribution, BC!=0",dcaName[i].Data()),
+ nptbins,ptmin,ptmax,ndcabins,mindca,maxdca);
+ fhPtDCATOFBCOut[i]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ fhPtDCATOFBCOut[i]->SetYTitle(Form("DCA_{%s}",dcaName[i].Data()));
+ outputContainer->Add(fhPtDCATOFBCOut[i]);
+ }
- fhPtUnknown = new TH1F ("hPtMCUnknown","p_T distribution from unknown", nptbins,ptmin,ptmax);
- fhPtUnknown->SetXTitle("p_{T} (GeV/c)");
- outputContainer->Add(fhPtUnknown);
+ fhPtDCANoTOFHit[i] = new TH2F(Form("hPtDCA%sNoTOFHit",dcaName[i].Data()),
+ Form("Track (no TOF hit) DCA%s vs #it{p}_{T}distribution",dcaName[i].Data()),
+ nptbins,ptmin,ptmax,ndcabins,mindca,maxdca);
+ fhPtDCANoTOFHit[i]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ fhPtDCANoTOFHit[i]->SetYTitle(Form("DCA_{%s}",dcaName[i].Data()));
+ outputContainer->Add(fhPtDCANoTOFHit[i]);
+
+ if(fFillVertexBC0Histograms)
+ {
+ fhPtDCAVtxOutBC0[i] = new TH2F(Form("hPtDCA%sVtxOutBC0",dcaName[i].Data()),
+ Form("Track DCA%s vs #it{p}_{T}distribution, vertex with BC!=0",dcaName[i].Data()),
+ nptbins,ptmin,ptmax,ndcabins,mindca,maxdca);
+ fhPtDCAVtxOutBC0[i]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ fhPtDCAVtxOutBC0[i]->SetYTitle(Form("DCA_{%s}",dcaName[i].Data()));
+ outputContainer->Add(fhPtDCAVtxOutBC0[i]);
+
+ fhPtDCAVtxOutBC0NoTOFHit[i] = new TH2F(Form("hPtDCA%sVtxOutBC0NoTOFHit",dcaName[i].Data()),
+ Form("Track (no TOF hit) DCA%s vs #it{p}_{T}distribution, vertex with BC!=0",dcaName[i].Data()),
+ nptbins,ptmin,ptmax,ndcabins,mindca,maxdca);
+ fhPtDCAVtxOutBC0NoTOFHit[i]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ fhPtDCAVtxOutBC0NoTOFHit[i]->SetYTitle(Form("DCA_{%s}",dcaName[i].Data()));
+ outputContainer->Add(fhPtDCAVtxOutBC0NoTOFHit[i]);
+
+ fhPtDCAVtxInBC0[i] = new TH2F(Form("hPtDCA%sVtxInBC0",dcaName[i].Data()),
+ Form("Track DCA%s vs #it{p}_{T}distribution, vertex with BC==0",dcaName[i].Data()),
+ nptbins,ptmin,ptmax,ndcabins,mindca,maxdca);
+ fhPtDCAVtxInBC0[i]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ fhPtDCAVtxInBC0[i]->SetYTitle(Form("DCA_{%s}",dcaName[i].Data()));
+ outputContainer->Add(fhPtDCAVtxInBC0[i]);
+
+ fhPtDCAVtxInBC0NoTOFHit[i] = new TH2F(Form("hPtDCA%sVtxInBC0NoTOFHit",dcaName[i].Data()),
+ Form("Track (no TOF hit) DCA%s vs #it{p}_{T}distribution, vertex with BC==0",dcaName[i].Data()),
+ nptbins,ptmin,ptmax,ndcabins,mindca,maxdca);
+ fhPtDCAVtxInBC0NoTOFHit[i]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ fhPtDCAVtxInBC0NoTOFHit[i]->SetYTitle(Form("DCA_{%s}",dcaName[i].Data()));
+ outputContainer->Add(fhPtDCAVtxInBC0NoTOFHit[i]);
+ }
- fhPhiUnknown = new TH2F ("hPhiMCUnknown","#phi distribution from unknown",nptbins,ptmin,ptmax, nphibins,phimin,phimax);
- fhPhiUnknown->SetXTitle("#phi (rad)");
- outputContainer->Add(fhPhiUnknown);
+ if(fFillPileUpHistograms)
+ {
+ fhPtDCAPileUp[i] = new TH2F(Form("hPtDCA%sPileUp",dcaName[i].Data()),
+ Form("Track DCA%s vs #it{p}_{T}distribution, SPD Pile-Up",dcaName[i].Data()),
+ nptbins,ptmin,ptmax,ndcabins,mindca,maxdca);
+ fhPtDCAPileUp[i]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ fhPtDCAPileUp[i]->SetYTitle(Form("DCA_{%s}",dcaName[i].Data()));
+ outputContainer->Add(fhPtDCAPileUp[i]);
+
+ if(fFillTrackBCHistograms)
+ {
+ fhPtDCAPileUpTOFBC0[i] = new TH2F(Form("hPtDCA%sPileUpTOFBC0",dcaName[i].Data()),
+ Form("Track DCA%s vs #it{p}_{T}distribution",dcaName[i].Data()),
+ nptbins,ptmin,ptmax,ndcabins,mindca,maxdca);
+ fhPtDCAPileUpTOFBC0[i]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ fhPtDCAPileUpTOFBC0[i]->SetYTitle(Form("DCA_{%s}",dcaName[i].Data()));
+ outputContainer->Add(fhPtDCAPileUpTOFBC0[i]);
+ }
+
+ fhPtDCAPileUpNoTOFHit[i] = new TH2F(Form("hPtDCA%sPileUpNoTOFHit",dcaName[i].Data()),
+ Form("Track (no TOF hit) DCA%s vs #it{p}_{T}distribution, SPD Pile-Up, vertex with BC!=0",dcaName[i].Data()),
+ nptbins,ptmin,ptmax,ndcabins,mindca,maxdca);
+ fhPtDCAPileUpNoTOFHit[i]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ fhPtDCAPileUpNoTOFHit[i]->SetYTitle(Form("DCA_{%s}",dcaName[i].Data()));
+ outputContainer->Add(fhPtDCAPileUpNoTOFHit[i]);
+
+ if(fFillVertexBC0Histograms)
+ {
+ fhPtDCAVtxOutBC0PileUp[i] = new TH2F(Form("hPtDCA%sPileUpVtxOutBC0",dcaName[i].Data()),
+ Form("Track DCA%s vs #it{p}_{T}distribution, SPD Pile-Up, vertex with BC!=0",dcaName[i].Data()),
+ nptbins,ptmin,ptmax,ndcabins,mindca,maxdca);
+ fhPtDCAVtxOutBC0PileUp[i]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ fhPtDCAVtxOutBC0PileUp[i]->SetYTitle(Form("DCA_{%s}",dcaName[i].Data()));
+ outputContainer->Add(fhPtDCAVtxOutBC0PileUp[i]);
+
+ fhPtDCAVtxOutBC0PileUpNoTOFHit[i] = new TH2F(Form("hPtDCA%sVtxOutBC0PileUpNoTOFHit",dcaName[i].Data()),
+ Form("Track (no TOF hit) DCA%s vs #it{p}_{T}distribution, SPD Pile-Up, vertex with BC!=0",dcaName[i].Data()),
+ nptbins,ptmin,ptmax,ndcabins,mindca,maxdca);
+ fhPtDCAVtxOutBC0PileUpNoTOFHit[i]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ fhPtDCAVtxOutBC0PileUpNoTOFHit[i]->SetYTitle(Form("DCA_{%s}",dcaName[i].Data()));
+ outputContainer->Add(fhPtDCAVtxOutBC0PileUpNoTOFHit[i]);
+
+ fhPtDCAVtxInBC0PileUp[i] = new TH2F(Form("hPtDCA%sPileUpVtxInBC0",dcaName[i].Data()),
+ Form("Track DCA%s vs #it{p}_{T}distribution, SPD Pile-Up,vertex with BC==0",dcaName[i].Data()),
+ nptbins,ptmin,ptmax,ndcabins,mindca,maxdca);
+ fhPtDCAVtxInBC0PileUp[i]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ fhPtDCAVtxInBC0PileUp[i]->SetYTitle(Form("DCA_{%s}",dcaName[i].Data()));
+ outputContainer->Add(fhPtDCAVtxInBC0PileUp[i]);
+
+ fhPtDCAVtxInBC0PileUpNoTOFHit[i] = new TH2F(Form("hPtDCA%sVtxInBC0PileUpNoTOFHit",dcaName[i].Data()),
+ Form("Track (no TOF hit) DCA%s vs #it{p}_{T}distribution, SPD Pile-Up, vertex with BC==0",dcaName[i].Data()),
+ nptbins,ptmin,ptmax,ndcabins,mindca,maxdca);
+ fhPtDCAVtxInBC0PileUpNoTOFHit[i]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ fhPtDCAVtxInBC0PileUpNoTOFHit[i]->SetYTitle(Form("DCA_{%s}",dcaName[i].Data()));
+ outputContainer->Add(fhPtDCAVtxInBC0PileUpNoTOFHit[i]);
+
+ }
+ }
+ }
+
+ if(IsDataMC())
+ {
+ //enum mvType{kmcPion = 0, kmcProton = 1, kmcKaon = 2, kmcMuon = 3, kmcElectron = 4, kmcUnknown = 4 };
+
+ TString histoName[] = {"Pion","Proton","Kaon","Muon","Electron","Unknown"};
+ TString titleName[] = {"#pi^{#pm}","p^{#pm}","K^{#pm}","#mu^{#pm}","e^{#pm}","x^{#pm}"};
- fhEtaUnknown = new TH2F ("hEtaMCUnknown","#eta distribution from unknown",nptbins,ptmin,ptmax, netabins,etamin,etamax);
- fhEtaUnknown->SetXTitle("#eta ");
- outputContainer->Add(fhEtaUnknown);
+ for(Int_t imcPart = 0; imcPart < 6; imcPart++)
+ {
+ fhPtMCPart[imcPart] = new TH1F (Form("hPtMC%s",histoName[imcPart].Data()),
+ Form("reconstructed #it{p}_{T} distribution from %s",titleName[imcPart].Data()),
+ nptbins,ptmin,ptmax);
+ fhPtMCPart[imcPart]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhPtMCPart[imcPart]);
+
+ fhPhiMCPart[imcPart] = new TH2F (Form("hPhiMC%s",histoName[imcPart].Data()),
+ Form("reconstructed #phi vs #it{p}_{T} distribution from %s",titleName[imcPart].Data()),
+ nptbins,ptmin,ptmax, nphibins,phimin,phimax);
+ fhPhiMCPart[imcPart]->SetYTitle("#phi (rad)");
+ fhPhiMCPart[imcPart]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhPhiMCPart[imcPart]);
+
+ fhEtaMCPart[imcPart] = new TH2F (Form("hEtaMC%s",histoName[imcPart].Data()),
+ Form("reconstructed #eta vs #it{p}_{T} distribution from %s",titleName[imcPart].Data()),
+ nptbins,ptmin,ptmax, netabins,etamin,etamax);
+ fhEtaMCPart[imcPart]->SetYTitle("#eta ");
+ fhEtaMCPart[imcPart]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhEtaMCPart[imcPart]);
+
+ fhPtMCPrimPart[imcPart] = new TH1F (Form("hPtMCPrimary%s",histoName[imcPart].Data()),
+ Form("generated #it{p}_{T} distribution from %s",titleName[imcPart].Data()),
+ nptbins,ptmin,ptmax);
+ fhPtMCPrimPart[imcPart]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhPtMCPrimPart[imcPart]);
+
+ fhPhiMCPrimPart[imcPart] = new TH2F (Form("hPhiMCPrimary%s",histoName[imcPart].Data()),
+ Form("generated #phi vs #it{p}_{T} distribution from %s",titleName[imcPart].Data()),
+ nptbins,ptmin,ptmax, nphibins,phimin,phimax);
+ fhPhiMCPrimPart[imcPart]->SetYTitle("#phi (rad)");
+ fhPhiMCPrimPart[imcPart]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhPhiMCPrimPart[imcPart]);
+
+ fhEtaMCPrimPart[imcPart] = new TH2F (Form("hEtaMCPrimary%s",histoName[imcPart].Data()),
+ Form("generated #eta vs #it{p}_{T} distribution from %s",titleName[imcPart].Data()),
+ nptbins,ptmin,ptmax, netabins,etamin,etamax);
+ fhEtaMCPrimPart[imcPart]->SetYTitle("#eta ");
+ fhEtaMCPrimPart[imcPart]->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhEtaMCPrimPart[imcPart]);
+ }
}
+ fhPtNPileUpSPDVtx = new TH2F ("hPt_NPileUpVertSPD","pT of cluster vs N pile-up SPD vertex",
+ nptbins,ptmin,ptmax,20,0,20);
+ fhPtNPileUpSPDVtx->SetYTitle("# vertex ");
+ fhPtNPileUpSPDVtx->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhPtNPileUpSPDVtx);
+
+ fhPtNPileUpTrkVtx = new TH2F ("hPt_NPileUpVertTracks","pT of cluster vs N pile-up Tracks vertex",
+ nptbins,ptmin,ptmax, 20,0,20 );
+ fhPtNPileUpTrkVtx->SetYTitle("# vertex ");
+ fhPtNPileUpTrkVtx->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhPtNPileUpTrkVtx);
+
+ if(fFillVertexBC0Histograms)
+ {
+ fhPtNPileUpSPDVtxBC0 = new TH2F ("hPt_NPileUpVertSPD_BC0","pT of cluster vs N pile-up SPD vertex",
+ nptbins,ptmin,ptmax,20,0,20);
+ fhPtNPileUpSPDVtxBC0->SetYTitle("# vertex ");
+ fhPtNPileUpSPDVtxBC0->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhPtNPileUpSPDVtxBC0);
+
+ fhPtNPileUpTrkVtxBC0 = new TH2F ("hPt_NPileUpVertTracks_BC0","pT of cluster vs N pile-up Tracks vertex",
+ nptbins,ptmin,ptmax, 20,0,20 );
+ fhPtNPileUpTrkVtxBC0->SetYTitle("# vertex ");
+ fhPtNPileUpTrkVtxBC0->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+ outputContainer->Add(fhPtNPileUpTrkVtxBC0);
+ }
+
return outputContainer;
}
+
//___________________________________________
void AliAnaChargedParticles::InitParameters()
{
SetOutputAODName("PWG4Particle");
AddToHistogramsName("AnaCharged_");
-
- fPdg = -1; //Select all tracks
}
printf("Min Pt = %3.2f\n", GetMinPt());
printf("Max Pt = %3.2f\n", GetMaxPt());
- printf("Select clusters with pdg %d \n",fPdg);
}
{
//Init
//Do some checks
- if(!GetReader()->IsCTSSwitchedOn()){
- printf("AliAnaChargedParticles::Init() - STOP!: You want to use CTS tracks in analysis but not read!! \n!!Check the configuration file!!\n");
- abort();
- }
+
+ if(!GetReader()->IsCTSSwitchedOn())
+ AliFatal("STOP!: You want to use CTS tracks in analysis but not read!! \n!!Check the configuration file!!\n");
+
}
if(GetDebug() > 0)
printf("AliAnaChargedParticles::MakeAnalysisFillAOD() - In CTS aod entries %d\n", ntracks);
+ AliVEvent * event = GetReader()->GetInputEvent();
+ AliESDEvent* esdEv = dynamic_cast<AliESDEvent*> (event);
+ AliAODEvent* aodEv = dynamic_cast<AliAODEvent*> (event);
+
+ Int_t vtxBC = GetReader()->GetVertexBC();
+ if(!GetReader()->IsDCACutOn()) vtxBC = GetReader()->GetVertexBC(event->GetPrimaryVertex());
+
+ if(fFillVertexBC0Histograms)
+ {
+ fhProductionVertexBC->Fill(vtxBC);
+ if(fFillPileUpHistograms)
+ {
+ if(GetReader()->IsPileUpFromSPD()) fhProductionVertexBCPileUp[0]->Fill(vtxBC);
+ if(GetReader()->IsPileUpFromEMCal()) fhProductionVertexBCPileUp[1]->Fill(vtxBC);
+ if(GetReader()->IsPileUpFromSPDOrEMCal()) fhProductionVertexBCPileUp[2]->Fill(vtxBC);
+ if(GetReader()->IsPileUpFromSPDAndEMCal()) fhProductionVertexBCPileUp[3]->Fill(vtxBC);
+ if(GetReader()->IsPileUpFromSPDAndNotEMCal()) fhProductionVertexBCPileUp[4]->Fill(vtxBC);
+ if(GetReader()->IsPileUpFromEMCalAndNotSPD()) fhProductionVertexBCPileUp[5]->Fill(vtxBC);
+ if(GetReader()->IsPileUpFromNotSPDAndNotEMCal()) fhProductionVertexBCPileUp[6]->Fill(vtxBC);
+ }
+ }
+
+ // N pile up vertices
+ Int_t nVtxSPD = -1;
+ Int_t nVtxTrk = -1;
+
+ if (esdEv)
+ {
+ nVtxSPD = esdEv->GetNumberOfPileupVerticesSPD();
+ nVtxTrk = esdEv->GetNumberOfPileupVerticesTracks();
+
+ }//ESD
+ else if (aodEv)
+ {
+ nVtxSPD = aodEv->GetNumberOfPileupVerticesSPD();
+ nVtxTrk = aodEv->GetNumberOfPileupVerticesTracks();
+ }//AOD
+
+
+ //printf("AliAnaChargedParticles::MakeAnalysisFillAOD() - primary vertex BC %d\n",vtxBC);
+
+ Double_t bz = event->GetMagneticField();
+
//Fill AODParticle with CTS aods
- TVector3 p3;
+ Float_t pt = 0;
+ Float_t phi = 0;
+ Float_t eta = 0;
Int_t evtIndex = 0;
- for(Int_t i = 0; i < ntracks; i++){
-
+ for(Int_t i = 0; i < ntracks; i++)
+ {
AliVTrack * track = (AliVTrack*) (GetCTSTracks()->At(i));
+
+ pt = track->Pt();
+ eta = track->Eta();
+ phi = track->Phi();
+
+ fhPtNoCut->Fill(pt);
+
+ fhPtNPileUpSPDVtx->Fill(pt,nVtxSPD);
+ fhPtNPileUpTrkVtx->Fill(pt,nVtxTrk);
+
+ AliAODTrack * aodTrack = dynamic_cast<AliAODTrack*>(track);
+ AliESDtrack * esdTrack = dynamic_cast<AliESDtrack*>(track);
+
+ //TOF
+ ULong_t status = track->GetStatus();
+ Bool_t okTOF = (status & AliVTrack::kTOFout) == AliVTrack::kTOFout ;
+ Double32_t tof = track->GetTOFsignal()*1e-3;
+
+ //DCA
+ Double_t dcaCons = -999;
+ if(aodTrack)
+ {
+ dcaCons = aodTrack->DCA();
+ //vtxBC = aodTrack->GetProdVertex()->GetBC();
+ }
+
+ Double_t dca[2] = {1e6,1e6};
+ Double_t covar[3] = {1e6,1e6,1e6};
+ track->PropagateToDCA(GetReader()->GetInputEvent()->GetPrimaryVertex(),bz,100.,dca,covar);
+
+ Float_t trackDCA = dca[0];
+
+ if(dcaCons == -999)
+ {
+ fhPtDCA[0]->Fill(pt, dca[0]);
+ fhPtDCA[1]->Fill(pt, dca[1]);
+ }
+ else
+ {
+ trackDCA = dcaCons;
+ fhPtDCA[2]->Fill(pt, dcaCons);
+ }
+
+ if(GetReader()->AcceptDCA(pt,trackDCA) && !GetReader()->IsDCACutOn() ) fhPtCutDCA->Fill(pt);
+
+ if(fFillVertexBC0Histograms)
+ {
+ if(TMath::Abs(vtxBC) > 0 && vtxBC!=AliVTrack::kTOFBCNA)
+ {
+ fhPtVtxOutBC0->Fill(pt);
+ fhEtaPhiVtxOutBC0->Fill(eta,phi);
+
+ if(dcaCons == -999)
+ {
+ fhPtDCAVtxOutBC0[0]->Fill(pt, dca[0]);
+ fhPtDCAVtxOutBC0[1]->Fill(pt, dca[1]);
+ }
+ else
+ fhPtDCAVtxOutBC0[2]->Fill(pt, dcaCons);
+ }
+ else
+ {
+ fhPtVtxInBC0->Fill(pt);
+ fhEtaPhiVtxInBC0->Fill(eta,phi);
+
+ fhPtNPileUpSPDVtxBC0->Fill(pt,nVtxSPD);
+ fhPtNPileUpTrkVtxBC0->Fill(pt,nVtxTrk);
+
+ if(GetReader()->AcceptDCA(pt,trackDCA) && !GetReader()->IsDCACutOn() ) fhPtCutDCABCOK->Fill(pt);
+
+ if(dcaCons == -999)
+ {
+ fhPtDCAVtxInBC0[0]->Fill(pt, dca[0]);
+ fhPtDCAVtxInBC0[1]->Fill(pt, dca[1]);
+ }
+ else
+ fhPtDCAVtxInBC0[2]->Fill(pt, dcaCons);
+
+ }
+ }
+
+ if(fFillPileUpHistograms && GetReader()->IsPileUpFromSPD())
+ {
+ if(dcaCons == -999)
+ {
+ fhPtDCAPileUp[0]->Fill(pt, dca[0]);
+ fhPtDCAPileUp[1]->Fill(pt, dca[1]);
+ }
+ else
+ fhPtDCAPileUp[2]->Fill(pt, dcaCons);
+
+ if(fFillVertexBC0Histograms)
+ {
+ if(TMath::Abs(vtxBC) > 0 && vtxBC!=AliVTrack::kTOFBCNA)
+ {
+ if(dcaCons == -999)
+ {
+ fhPtDCAVtxOutBC0PileUp[0]->Fill(pt, dca[0]);
+ fhPtDCAVtxOutBC0PileUp[1]->Fill(pt, dca[1]);
+ }
+ else fhPtDCAVtxOutBC0PileUp[2]->Fill(pt, dcaCons);
+ }
+ else
+ {
+ if(dcaCons == -999)
+ {
+ fhPtDCAVtxInBC0PileUp[0]->Fill(pt, dca[0]);
+ fhPtDCAVtxInBC0PileUp[1]->Fill(pt, dca[1]);
+ }
+ else fhPtDCAVtxInBC0PileUp[2]->Fill(pt, dcaCons);
+ }
+ }
+ }
+
+ if(!okTOF)
+ {
+ if(dcaCons == -999)
+ {
+ fhPtDCANoTOFHit[0]->Fill(pt, dca[0]);
+ fhPtDCANoTOFHit[1]->Fill(pt, dca[1]);
+ }
+ else
+ fhPtDCANoTOFHit[2]->Fill(pt, dcaCons);
+
+ if(fFillVertexBC0Histograms)
+ {
+ if(TMath::Abs(vtxBC) > 0 && vtxBC!=AliVTrack::kTOFBCNA)
+ {
+ if(dcaCons == -999)
+ {
+ fhPtDCAVtxOutBC0NoTOFHit[0]->Fill(pt, dca[0]);
+ fhPtDCAVtxOutBC0NoTOFHit[1]->Fill(pt, dca[1]);
+ }
+ else
+ fhPtDCAVtxOutBC0NoTOFHit[2]->Fill(pt, dcaCons);
+ }
+ else
+ {
+ if(dcaCons == -999)
+ {
+ fhPtDCAVtxInBC0NoTOFHit[0]->Fill(pt, dca[0]);
+ fhPtDCAVtxInBC0NoTOFHit[1]->Fill(pt, dca[1]);
+ }
+ else
+ fhPtDCAVtxInBC0NoTOFHit[2]->Fill(pt, dcaCons);
+
+ }
+ }
+
+ if(fFillPileUpHistograms && GetReader()->IsPileUpFromSPD())
+ {
+ if(dcaCons == -999)
+ {
+ fhPtDCAPileUpNoTOFHit[0]->Fill(pt, dca[0]);
+ fhPtDCAPileUpNoTOFHit[1]->Fill(pt, dca[1]);
+ }
+ else
+ fhPtDCAPileUpNoTOFHit[2]->Fill(pt, dcaCons);
+
+ if(fFillVertexBC0Histograms)
+ {
+ if(TMath::Abs(vtxBC) > 0 && vtxBC!=AliVTrack::kTOFBCNA)
+ {
+ if(dcaCons == -999)
+ {
+ fhPtDCAVtxOutBC0PileUpNoTOFHit[0]->Fill(pt, dca[0]);
+ fhPtDCAVtxOutBC0PileUpNoTOFHit[1]->Fill(pt, dca[1]);
+ }
+ else
+ fhPtDCAVtxOutBC0PileUpNoTOFHit[2]->Fill(pt, dcaCons);
+ }
+ else
+ {
+ if(dcaCons == -999)
+ {
+ fhPtDCAVtxInBC0PileUpNoTOFHit[0]->Fill(pt, dca[0]);
+ fhPtDCAVtxInBC0PileUpNoTOFHit[1]->Fill(pt, dca[1]);
+ }
+ else
+ fhPtDCAVtxInBC0PileUpNoTOFHit[2]->Fill(pt, dcaCons);
+
+ }
+ }
+ }
+ }
+
+ //printf("track pT %2.2f, DCA Cons %f, DCA1 %f, DCA2 %f, TOFBC %d, oktof %d, tof %f\n",
+ // pt,dcaCons,dca[0],dca[1],track->GetTOFBunchCrossing(bz),okTOF, tof);
+
+
+// if( vtxBC == 0 && trackBC !=0 && trackBC!=AliVTrack::kTOFBCNA)
+// printf("TOF Signal %e, BC %d, pt %f, dca_xy %f, dca_z %f, dca_tpc %f \n", tof,trackBC, pt,dca[0],dca[1],dcaCons);
+
- //Fill AODParticle after some selection
+ if(okTOF)
+ {
+ fhTOFSignal ->Fill(tof);
+ fhPtTOFSignal->Fill(pt, tof);
+ if(GetReader()->AcceptDCA(pt,trackDCA) && !GetReader()->IsDCACutOn() ) fhPtTOFSignalDCACut->Fill(pt, tof);
+
+ if(fFillVertexBC0Histograms)
+ {
+ if(TMath::Abs(vtxBC) > 0 && vtxBC!=AliVTrack::kTOFBCNA)
+ fhPtTOFSignalVtxOutBC0->Fill(pt, tof);
+ else
+ fhPtTOFSignalVtxInBC0->Fill(pt, tof);
+ }
+
+ Int_t trackBC = 1000;
+
+ if(fFillTrackBCHistograms)
+ {
+ trackBC = track->GetTOFBunchCrossing(bz);
+
+ if(trackBC==0)
+ {
+ fhTOFSignalBCOK->Fill(tof);
+
+ if(dcaCons == -999)
+ {
+ fhPtDCATOFBC0[0]->Fill(pt, dca[0]);
+ fhPtDCATOFBC0[1]->Fill(pt, dca[1]);
+ }
+ else
+ fhPtDCATOFBC0[2]->Fill(pt, dcaCons);
+
+ if(fFillPileUpHistograms && GetReader()->IsPileUpFromSPD())
+ {
+ if(dcaCons == -999)
+ {
+ fhPtDCAPileUpTOFBC0[0]->Fill(pt, dca[0]);
+ fhPtDCAPileUpTOFBC0[1]->Fill(pt, dca[1]);
+ }
+ else
+ fhPtDCAPileUpTOFBC0[2]->Fill(pt, dcaCons);
+ }
+ }
+ else if(trackBC!=AliVTrack::kTOFBCNA)
+ {
+ if(dcaCons == -999)
+ {
+ fhPtDCATOFBCOut[0]->Fill(pt, dca[0]);
+ fhPtDCATOFBCOut[1]->Fill(pt, dca[1]);
+ }
+ else
+ fhPtDCATOFBCOut[2]->Fill(pt, dcaCons);
+
+ }
+ }
+
+ if(fFillPileUpHistograms)
+ {
+ if(GetReader()->IsPileUpFromSPD()) fhPtTOFSignalPileUp[0]->Fill(pt, tof);
+ if(GetReader()->IsPileUpFromEMCal()) fhPtTOFSignalPileUp[1]->Fill(pt, tof);
+ if(GetReader()->IsPileUpFromSPDOrEMCal()) fhPtTOFSignalPileUp[2]->Fill(pt, tof);
+ if(GetReader()->IsPileUpFromSPDAndEMCal()) fhPtTOFSignalPileUp[3]->Fill(pt, tof);
+ if(GetReader()->IsPileUpFromSPDAndNotEMCal()) fhPtTOFSignalPileUp[4]->Fill(pt, tof);
+ if(GetReader()->IsPileUpFromEMCalAndNotSPD()) fhPtTOFSignalPileUp[5]->Fill(pt, tof);
+ if(GetReader()->IsPileUpFromNotSPDAndNotEMCal()) fhPtTOFSignalPileUp[6]->Fill(pt, tof);
+
+ if(fFillTrackBCHistograms)
+ {
+ if (trackBC ==0) { fhEtaPhiTOFBC0 ->Fill(eta,phi); if(fFillPileUpHistograms && GetReader()->IsPileUpFromSPD()) fhEtaPhiTOFBC0PileUpSPD ->Fill(eta,phi); }
+ else if (trackBC < 0) { fhEtaPhiTOFBCPlus ->Fill(eta,phi); if(fFillPileUpHistograms && GetReader()->IsPileUpFromSPD()) fhEtaPhiTOFBCPlusPileUpSPD ->Fill(eta,phi); }
+ else if (trackBC > 0) { fhEtaPhiTOFBCMinus->Fill(eta,phi); if(fFillPileUpHistograms && GetReader()->IsPileUpFromSPD()) fhEtaPhiTOFBCMinusPileUpSPD->Fill(eta,phi); }
+ }
+
+ if(fFillVertexBC0Histograms)
+ {
+ if(TMath::Abs(vtxBC) > 0 && vtxBC!=AliVTrack::kTOFBCNA)
+ {
+ if(GetReader()->IsPileUpFromSPD()) fhPtTOFSignalVtxOutBC0PileUp[0]->Fill(pt, tof);
+ if(GetReader()->IsPileUpFromEMCal()) fhPtTOFSignalVtxOutBC0PileUp[1]->Fill(pt, tof);
+ if(GetReader()->IsPileUpFromSPDOrEMCal()) fhPtTOFSignalVtxOutBC0PileUp[2]->Fill(pt, tof);
+ if(GetReader()->IsPileUpFromSPDAndEMCal()) fhPtTOFSignalVtxOutBC0PileUp[3]->Fill(pt, tof);
+ if(GetReader()->IsPileUpFromSPDAndNotEMCal()) fhPtTOFSignalVtxOutBC0PileUp[4]->Fill(pt, tof);
+ if(GetReader()->IsPileUpFromEMCalAndNotSPD()) fhPtTOFSignalVtxOutBC0PileUp[5]->Fill(pt, tof);
+ if(GetReader()->IsPileUpFromNotSPDAndNotEMCal()) fhPtTOFSignalVtxOutBC0PileUp[6]->Fill(pt, tof);
+ }
+ else
+ {
+ if(GetReader()->IsPileUpFromSPD()) fhPtTOFSignalVtxInBC0PileUp[0]->Fill(pt, tof);
+ if(GetReader()->IsPileUpFromEMCal()) fhPtTOFSignalVtxInBC0PileUp[1]->Fill(pt, tof);
+ if(GetReader()->IsPileUpFromSPDOrEMCal()) fhPtTOFSignalVtxInBC0PileUp[2]->Fill(pt, tof);
+ if(GetReader()->IsPileUpFromSPDAndEMCal()) fhPtTOFSignalVtxInBC0PileUp[3]->Fill(pt, tof);
+ if(GetReader()->IsPileUpFromSPDAndNotEMCal()) fhPtTOFSignalVtxInBC0PileUp[4]->Fill(pt, tof);
+ if(GetReader()->IsPileUpFromEMCalAndNotSPD()) fhPtTOFSignalVtxInBC0PileUp[5]->Fill(pt, tof);
+ if(GetReader()->IsPileUpFromNotSPDAndNotEMCal()) fhPtTOFSignalVtxInBC0PileUp[6]->Fill(pt, tof);
+ }
+ }
+ }
+ }
+
+ //Fill AODParticle after some selection
Double_t mom[3] = {track->Px(),track->Py(),track->Pz()};
- p3.SetXYZ(mom[0],mom[1],mom[2]);
Bool_t in = GetFiducialCut()->IsInFiducialCut(mom,"CTS") ;
if(GetDebug() > 1)
- printf("AliAnaChargedParticles::MakeAnalysisFillAOD() - Track pt %2.2f, phi %2.2f, eta %2.2f in fiducial cut %d\n",p3.Pt(), p3.Phi(), p3.Eta(),in);
+ printf("AliAnaChargedParticles::MakeAnalysisFillAOD() - Track pt %2.2f, eta %2.2f, phi %2.2f in fiducial cut %d\n",pt,eta,phi,in);
//Acceptance selection
if(IsFiducialCutOn() && ! in ) continue ;
// Momentum selection
- if(p3.Pt() < GetMinPt() || p3.Pt() > GetMaxPt()) continue;
+ if(pt < GetMinPt() || pt > GetMaxPt()) continue;
+
+ if(okTOF) fhTOFSignalPtCut->Fill(tof);
+ else
+ {
+ fhPtTOFStatus0 ->Fill(pt);
+ fhEtaPhiTOFStatus0->Fill(eta,phi);
+ }
- //Keep only particles identified with fPdg
- //Selection not done for the moment
- //Should be done here.
+ Bool_t bITSRefit = (status & AliVTrack::kITSrefit) == AliVTrack::kITSrefit;
+ Bool_t bConstrained = kFALSE;
+ if (aodTrack) bConstrained = aodTrack->IsGlobalConstrained();
+ else if(esdTrack) bConstrained = (!esdTrack->HasPointOnITSLayer(0) && !esdTrack->HasPointOnITSLayer(1));
+ //printf("Track %d, pt %2.2f, eta %2.2f, phi %2.2f, SPDRefit %d, refit %d, dcaCons %2.2f\n",
+ // i, pt, eta, phi, bConstrained, bITSRefit, dcaCons);
+ if(bConstrained)
+ {
+ if(bITSRefit)
+ {
+ fhPtNoSPDRefit->Fill(pt);
+ if(pt < 2)fhEtaPhiNoSPDRefitPt02->Fill(eta,phi);
+ if(pt > 3)fhEtaPhiNoSPDRefitPt3 ->Fill(eta,phi);
+
+ if(dcaCons == -999)
+ {
+ fhPtDCANoSPDRefit[0]->Fill(pt, dca[0]);
+ fhPtDCANoSPDRefit[1]->Fill(pt, dca[1]);
+ }
+ else
+ fhPtDCANoSPDRefit[2]->Fill(pt, dcaCons);
+
+ }
+ else
+ {
+ fhPtNoSPDNoRefit->Fill(pt);
+ if(pt < 2)fhEtaPhiNoSPDNoRefitPt02->Fill(eta,phi);
+ if(pt > 3)fhEtaPhiNoSPDNoRefitPt3 ->Fill(eta,phi);
+ if(dcaCons == -999)
+ {
+ fhPtDCANoSPDNoRefit[0]->Fill(pt, dca[0]);
+ fhPtDCANoSPDNoRefit[1]->Fill(pt, dca[1]);
+ }
+ else
+ fhPtDCANoSPDNoRefit[2]->Fill(pt, dcaCons);
+
+ }
+ }
+ else
+ {
+ fhPtSPDRefit->Fill(pt);
+ if(pt < 2)fhEtaPhiSPDRefitPt02->Fill(eta,phi);
+ if(pt > 3)fhEtaPhiSPDRefitPt3 ->Fill(eta,phi);
+ if(dcaCons == -999)
+ {
+ fhPtDCASPDRefit[0]->Fill(pt, dca[0]);
+ fhPtDCASPDRefit[1]->Fill(pt, dca[1]);
+ }
+ else
+ fhPtDCASPDRefit[2]->Fill(pt, dcaCons);
+ }
+
// Mixed event
- if (GetMixedEvent()){
+ if (GetMixedEvent())
+ {
evtIndex = GetMixedEvent()->EventIndex(track->GetID()) ;
- }
+ }
+
GetVertex(vert,evtIndex);
if(TMath::Abs(vert[2])> GetZvertexCut()) return;
}
//__________________________________________________________________
-void AliAnaChargedParticles::MakeAnalysisFillHistograms()
+void AliAnaChargedParticles::MakeAnalysisFillHistograms()
{
//Do analysis and fill histograms
+ if(IsDataMC()) FillPrimaryHistograms();
+
//Loop on stored AODParticles
Int_t naod = GetOutputAODBranch()->GetEntriesFast();
fhNtracks->Fill(GetReader()->GetTrackMultiplicity()) ;
- if(GetDebug() > 0)
+ if(GetDebug() > 0)
printf("AliAnaChargedParticles::MakeAnalysisFillHistograms() - aod branch entries %d\n", naod);
- for(Int_t iaod = 0; iaod < naod ; iaod++){
- AliAODPWG4Particle* tr = (AliAODPWG4Particle*) (GetOutputAODBranch()->At(iaod));
-
- fhPt->Fill(tr->Pt());
+ Float_t pt = 0;
+ Float_t phi = 0;
+ Float_t eta = 0;
+
+ for(Int_t iaod = 0; iaod < naod ; iaod++)
+ {
+ AliAODPWG4Particle* track = (AliAODPWG4Particle*) (GetOutputAODBranch()->At(iaod));
+
+ pt = track->Pt();
+ eta = track->Eta();
+ phi = track->Phi();
- if(tr->GetChargedBit()){
- fhPhiPos ->Fill(tr->Pt(), tr->Phi());
- fhEtaPos ->Fill(tr->Pt(), tr->Eta());
- fhEtaPhiPos->Fill(tr->Eta(),tr->Phi());
+ fhPt->Fill(pt);
+
+ if(track->GetChargedBit())
+ {
+ fhPhiPos ->Fill(pt, phi);
+ fhEtaPos ->Fill(pt, eta);
+ fhEtaPhiPos->Fill(eta,phi);
}
- else{
- fhPhiNeg ->Fill(tr->Pt(), tr->Phi());
- fhEtaNeg ->Fill(tr->Pt(), tr->Eta());
- fhEtaPhiNeg->Fill(tr->Eta(),tr->Phi());
+ else
+ {
+ fhPhiNeg ->Fill(pt, phi);
+ fhEtaNeg ->Fill(pt, eta);
+ fhEtaPhiNeg->Fill(eta,phi);
}
- if(IsDataMC()){
- //Play with the MC stack if available
+ if(fFillPileUpHistograms)
+ {
+ if(GetReader()->IsPileUpFromSPD()) {fhPtPileUp[0]->Fill(pt);}
+ if(GetReader()->IsPileUpFromEMCal()) {fhPtPileUp[1]->Fill(pt);}
+ if(GetReader()->IsPileUpFromSPDOrEMCal()) {fhPtPileUp[2]->Fill(pt);}
+ if(GetReader()->IsPileUpFromSPDAndEMCal()) {fhPtPileUp[3]->Fill(pt);}
+ if(GetReader()->IsPileUpFromSPDAndNotEMCal()) {fhPtPileUp[4]->Fill(pt);}
+ if(GetReader()->IsPileUpFromEMCalAndNotSPD()) {fhPtPileUp[5]->Fill(pt);}
+ if(GetReader()->IsPileUpFromNotSPDAndNotEMCal()) {fhPtPileUp[6]->Fill(pt);}
+ }
+
+
+ if(IsDataMC())
+ {
+ //Play with the MC stack if available
Int_t mompdg = -1;
- Int_t label = TMath::Abs(tr->GetLabel());
- if(GetReader()->ReadStack()){
- TParticle * mom = GetMCStack()->Particle(label);
- mompdg =TMath::Abs(mom->GetPdgCode());
- }
- else if(GetReader()->ReadAODMCParticles()){
- AliAODMCParticle * aodmom = 0;
- //Get the list of MC particles
- aodmom = (AliAODMCParticle*) (GetReader()->GetAODMCParticles(tr->GetInputFileIndex()))->At(label);
- mompdg =TMath::Abs(aodmom->GetPdgCode());
- }
+ Int_t label = track->GetLabel();
- if(mompdg==211){
- fhPtPion ->Fill(tr->Pt());
- fhPhiPion->Fill(tr->Pt(), tr->Phi());
- fhEtaPion->Fill(tr->Pt(), tr->Eta());
- }
- else if(mompdg==2212){
- fhPtProton ->Fill(tr->Pt());
- fhPhiProton->Fill(tr->Pt(), tr->Phi());
- fhEtaProton->Fill(tr->Pt(), tr->Eta());
- }
- else if(mompdg==321){
- fhPtKaon ->Fill(tr->Pt());
- fhPhiKaon->Fill(tr->Pt(), tr->Phi());
- fhEtaKaon->Fill(tr->Pt(), tr->Eta());
- }
- else if(mompdg==11){
- fhPtElectron ->Fill(tr->Pt());
- fhPhiElectron->Fill(tr->Pt(), tr->Phi());
- fhEtaElectron->Fill(tr->Pt(), tr->Eta());
- }
- else {
- //printf("unknown pdg %d\n",mompdg);
- fhPtUnknown ->Fill(tr->Pt());
- fhPhiUnknown->Fill(tr->Pt(), tr->Phi());
- fhEtaUnknown->Fill(tr->Pt(), tr->Eta());
+ if(label >= 0)
+ {
+ if( GetReader()->ReadStack() && label < GetMCStack()->GetNtrack())
+ {
+ TParticle * mom = GetMCStack()->Particle(label);
+ mompdg =TMath::Abs(mom->GetPdgCode());
+ }
+ else if(GetReader()->ReadAODMCParticles())
+ {
+ AliAODMCParticle * aodmom = 0;
+ //Get the list of MC particles
+ aodmom = (AliAODMCParticle*) (GetReader()->GetAODMCParticles())->At(label);
+ mompdg =TMath::Abs(aodmom->GetPdgCode());
+ }
}
+ Int_t mcType = kmcUnknown;
+ if (mompdg==211 ) mcType = kmcPion;
+ else if(mompdg==2212) mcType = kmcProton;
+ else if(mompdg==321 ) mcType = kmcKaon;
+ else if(mompdg==11 ) mcType = kmcElectron;
+ else if(mompdg==13 ) mcType = kmcMuon;
+
+ fhPtMCPart [mcType]->Fill(pt);
+ fhEtaMCPart[mcType]->Fill(pt,eta);
+ fhPhiMCPart[mcType]->Fill(pt,phi);
+
}//Work with stack also
}// aod branch loop
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff