From 95bd6d604cf04556f070b80b8e699b0a1f4df490 Mon Sep 17 00:00:00 2001 From: fbock Date: Wed, 20 Aug 2014 15:46:13 +0200 Subject: [PATCH] - added histo for validated two pion mass for same mother --- ...alysisTaskNeutralMesonToPiPlPiMiPiZero.cxx | 21 +++++++++++++------ ...AnalysisTaskNeutralMesonToPiPlPiMiPiZero.h | 1 + 2 files changed, 16 insertions(+), 6 deletions(-) diff --git a/PWGGA/GammaConv/AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.cxx b/PWGGA/GammaConv/AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.cxx index b6c145e30d8..a4f0166dcca 100644 --- a/PWGGA/GammaConv/AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.cxx +++ b/PWGGA/GammaConv/AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.cxx @@ -123,6 +123,7 @@ AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::AliAnalysisTaskNeutralMesonToPiPlPi fHistoTrueNegPionPt(NULL), fHistoTrueNegPionFromNeutralMesonPt(NULL), fHistoTruePionPionInvMassPt(NULL), + fHistoTruePionPionFromSameMotherInvMassPt(NULL), fHistoTruePionPionFromNeutralMesonInvMassPt(NULL), fHistoNEvents(NULL), fHistoNGoodESDTracks(NULL), @@ -210,6 +211,7 @@ AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::AliAnalysisTaskNeutralMesonToPiPlPi fHistoTrueNegPionPt(NULL), fHistoTrueNegPionFromNeutralMesonPt(NULL), fHistoTruePionPionInvMassPt(NULL), + fHistoTruePionPionFromSameMotherInvMassPt(NULL), fHistoTruePionPionFromNeutralMesonInvMassPt(NULL), fHistoNEvents(NULL), fHistoNGoodESDTracks(NULL), @@ -507,6 +509,7 @@ void AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::UserCreateOutputObjects() // if (fDoMesonQA){ fHistoTruePionPionInvMassPt = new TH2F*[fnCuts]; + fHistoTruePionPionFromSameMotherInvMassPt = new TH2F*[fnCuts]; fHistoTruePionPionFromNeutralMesonInvMassPt = new TH2F*[fnCuts]; // } @@ -588,6 +591,8 @@ void AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::UserCreateOutputObjects() // if (fDoMesonQA){ fHistoTruePionPionInvMassPt[iCut] = new TH2F("ESD_TruePiPlusPiNeg_InvMassPt","ESD_TruePiPlusPiNeg_InvMassPt",2000,0.,2.,200,0.,20.); fTrueList[iCut]->Add(fHistoTruePionPionInvMassPt[iCut]); + fHistoTruePionPionFromSameMotherInvMassPt[iCut] = new TH2F("ESD_TruePiPlusPiNegFromSameMother_InvMassPt","ESD_TruePiPlusPiNegFromSameMother_InvMassPt",2000,0.,2.,200,0.,20.); + fTrueList[iCut]->Add(fHistoTruePionPionFromSameMotherInvMassPt[iCut]); fHistoTruePionPionFromNeutralMesonInvMassPt[iCut] = new TH2F("ESD_TruePiPlusPiNegFromNeutralMeson_InvMassPt","ESD_TruePiPlusPiNegFromNeutralMeson_InvMassPt",2000,0.,2.,200,0.,20.); fTrueList[iCut]->Add(fHistoTruePionPionFromNeutralMesonInvMassPt[iCut]); // } @@ -1194,13 +1199,15 @@ void AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::ProcessPionCandidates(){ Int_t labeln=0; Int_t labelp=0; - + Int_t motherlabelp = 0; + Int_t motherlabeln = 0; if( fMCEvent ) { labeln=TMath::Abs(negPionCandidate->GetLabel()); labelp=TMath::Abs(posPionCandidate->GetLabel()); TParticle *fNegativeMCParticle = fMCStack->Particle(labeln); TParticle *fPositiveMCParticle = fMCStack->Particle(labelp); - + motherlabeln = fNegativeMCParticle->GetMother(0); + motherlabelp = fPositiveMCParticle->GetMother(0); if( fPositiveMCParticle && fNegativeMCParticle) { virtualPhoton->SetMCLabelPositive(labelp); virtualPhoton->SetMCLabelNegative(labeln); @@ -1218,10 +1225,12 @@ void AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::ProcessPionCandidates(){ } if (fMCEvent){ - if( (IsEtaPiPlPiMiPiZeroDaughter(labeln) && IsEtaPiPlPiMiPiZeroDaughter(labelp) ) || (IsOmegaPiPlPiMiPiZeroDaughter(labeln) && IsOmegaPiPlPiMiPiZeroDaughter(labelp)) ) { - fHistoTruePionPionFromNeutralMesonInvMassPt[fiCut]->Fill(vParticle->GetMass(),vParticle->Pt()); - } - + if (motherlabeln == motherlabelp){ + fHistoTruePionPionFromSameMotherInvMassPt[fiCut]->Fill(vParticle->GetMass(),vParticle->Pt()); + if( IsEtaPiPlPiMiPiZeroDaughter(labeln) || IsOmegaPiPlPiMiPiZeroDaughter(labeln) ) { + fHistoTruePionPionFromNeutralMesonInvMassPt[fiCut]->Fill(vParticle->GetMass(),vParticle->Pt()); + } + } } diff --git a/PWGGA/GammaConv/AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.h b/PWGGA/GammaConv/AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.h index eb33dc04991..2d1b1eab899 100644 --- a/PWGGA/GammaConv/AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.h +++ b/PWGGA/GammaConv/AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.h @@ -162,6 +162,7 @@ class AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero: public AliAnalysisTaskSE TH1F **fHistoTrueNegPionPt; // histos with reconstructed validated negative pion, pT TH1F **fHistoTrueNegPionFromNeutralMesonPt; // histos with reconstructed validated negative pion from eta or omega, pT TH2F **fHistoTruePionPionInvMassPt; // histos with reconstructed validated two pion, invariant mass, pT + TH2F **fHistoTruePionPionFromSameMotherInvMassPt;// histos with reconstructed validated two pion from same mother, invariant mass, pT TH2F **fHistoTruePionPionFromNeutralMesonInvMassPt;// histos with reconstructed validated two pion from eta or omega, invariant mass, pT // Event properties TH1I **fHistoNEvents; // histo for event counting -- 2.39.3