From 99b8e903fa1c52ec68b01acf27e6c805bebfbd80 Mon Sep 17 00:00:00 2001 From: gconesab Date: Fri, 23 Sep 2011 18:02:02 +0000 Subject: [PATCH] Add invariant mass histograms to check from MC the mass from conversions, check needed for clusterization studies --- PWG4/PartCorrDep/AliAnaPi0.cxx | 42 ++++++++++++++++++++++++++++++++-- PWG4/PartCorrDep/AliAnaPi0.h | 10 ++++++-- 2 files changed, 48 insertions(+), 4 deletions(-) diff --git a/PWG4/PartCorrDep/AliAnaPi0.cxx b/PWG4/PartCorrDep/AliAnaPi0.cxx index dd26370b94b..78d12b7b9be 100755 --- a/PWG4/PartCorrDep/AliAnaPi0.cxx +++ b/PWG4/PartCorrDep/AliAnaPi0.cxx @@ -87,7 +87,8 @@ fhPrimEtaPhi(0x0), fhPrimEtaAccPhi(0x0), fhPrimPi0PtOrigin(0x0 fhMCOrgMass(), fhMCOrgAsym(), fhMCOrgDeltaEta(), fhMCOrgDeltaPhi(), fhMCPi0MassPtRec(), fhMCPi0MassPtTrue(), fhMCPi0PtTruePtRec(), fhMCEtaMassPtRec(), fhMCEtaMassPtTrue(), fhMCEtaPtTruePtRec(), -fhMCPi0PtOrigin(0x0), fhMCEtaPtOrigin(0x0) +fhMCPi0PtOrigin(0x0), fhMCEtaPtOrigin(0x0), +fhReMCFromConversion(0), fhReMCFromNotConversion(0), fhReMCFromMixConversion(0) { //Default Ctor InitParameters(); @@ -625,6 +626,25 @@ TList * AliAnaPi0::GetCreateOutputObjects() //Histograms filled only if MC data is requested if(IsDataMC()){ + + fhReMCFromConversion = new TH2F("hReMCFromConversion","Invariant mass of 2 clusters originated in conversions", + nptbins,ptmin,ptmax,nmassbins,massmin,massmax); + fhReMCFromConversion->SetXTitle("p_{T} (GeV/c)"); + fhReMCFromConversion->SetYTitle("m_{#gamma,#gamma} (GeV/c^{2})"); + outputContainer->Add(fhReMCFromConversion) ; + + fhReMCFromNotConversion = new TH2F("hReMCNotFromConversion","Invariant mass of 2 clusters not originated in conversions", + nptbins,ptmin,ptmax,nmassbins,massmin,massmax); + fhReMCFromNotConversion->SetXTitle("p_{T} (GeV/c)"); + fhReMCFromNotConversion->SetYTitle("m_{#gamma,#gamma} (GeV/c^{2})"); + outputContainer->Add(fhReMCFromNotConversion) ; + + fhReMCFromMixConversion = new TH2F("hReMCFromMixConversion","Invariant mass of 2 clusters one from conversion and the other not", + nptbins,ptmin,ptmax,nmassbins,massmin,massmax); + fhReMCFromMixConversion->SetXTitle("p_{T} (GeV/c)"); + fhReMCFromMixConversion->SetYTitle("m_{#gamma,#gamma} (GeV/c^{2})"); + outputContainer->Add(fhReMCFromMixConversion) ; + //Pi0 fhPrimPi0Pt = new TH1F("hPrimPi0Pt","Primary pi0 pt, Y<1",nptbins,ptmin,ptmax) ; fhPrimPi0AccPt = new TH1F("hPrimPi0AccPt","Primary pi0 pt with both photons in acceptance",nptbins,ptmin,ptmax) ; @@ -1920,7 +1940,25 @@ void AliAnaPi0::MakeAnalysisFillHistograms() // MC data //--------- //Do some MC checks on the origin of the pair, is there any common ancestor and if there is one, who? - if(IsDataMC()) FillMCVersusRecDataHistograms(p1->GetLabel(), p2->GetLabel(),p1->Pt(), p2->Pt(),ncell1, ncell2, m, pt, a,deta, dphi); + if(IsDataMC()) { + + if(GetMCAnalysisUtils()->CheckTagBit(p1->GetTag(),AliMCAnalysisUtils::kMCConversion) && + GetMCAnalysisUtils()->CheckTagBit(p2->GetTag(),AliMCAnalysisUtils::kMCConversion)) + { + fhReMCFromConversion->Fill(pt,m); + } + else if(!GetMCAnalysisUtils()->CheckTagBit(p1->GetTag(),AliMCAnalysisUtils::kMCConversion) && + !GetMCAnalysisUtils()->CheckTagBit(p2->GetTag(),AliMCAnalysisUtils::kMCConversion)) + { + fhReMCFromNotConversion->Fill(pt,m); + } + else + { + fhReMCFromMixConversion->Fill(pt,m); + } + + FillMCVersusRecDataHistograms(p1->GetLabel(), p2->GetLabel(),p1->Pt(), p2->Pt(),ncell1, ncell2, m, pt, a,deta, dphi); + } //----------------------- //Multi cuts analysis diff --git a/PWG4/PartCorrDep/AliAnaPi0.h b/PWG4/PartCorrDep/AliAnaPi0.h index 1814d4a16c1..33ac1f476d4 100755 --- a/PWG4/PartCorrDep/AliAnaPi0.h +++ b/PWG4/PartCorrDep/AliAnaPi0.h @@ -228,7 +228,7 @@ class AliAnaPi0 : public AliAnaPartCorrBaseClass { //Multiple cuts: Assymmetry, pt, n cells, PID TH2F ** fhRePtNCellAsymCuts ; //![fNPtCuts*fNAsymCuts*fNCellNCuts*] REAL two-photon invariant mass distribution for different pt cut, n cell cuts and assymetry TH2F ** fhMiPtNCellAsymCuts ; //![fNPtCuts*fNAsymCuts*fNCellNCuts] Mixed two-photon invariant mass distribution for different pt cut, n cell cuts and assymetry - TH2F ** fhRePtNCellAsymCutsSM[12] ; //![fNPtCuts*fNAsymCuts*fNCellNCutsfNModules] REAL two-photon invariant mass distribution for different pt cut, n cell cuts and assymetry for each module + TH2F ** fhRePtNCellAsymCutsSM[12] ; //![fNPtCuts*fNAsymCuts*fNCellNCutsfNModules] REAL two-photon invariant mass distribution for different pt cut, n cell cuts and assymetry for each module TH2F ** fhRePIDBits ; //![fNPIDBits] REAL two-photon invariant mass distribution for different PID bits TH3F ** fhRePtMult ; //![fNAsymCuts] REAL two-photon invariant mass distribution for different track multiplicity and assymetry cuts @@ -294,7 +294,13 @@ class AliAnaPi0 : public AliAnaPartCorrBaseClass { TH2F * fhMCPi0PtOrigin ; //! Mass of reoconstructed pi0 pairs in calorimeter vs mother TH2F * fhMCEtaPtOrigin ; //! Mass of reoconstructed pi0 pairs in calorimeter vs mother - ClassDef(AliAnaPi0,20) + TH2F * fhReMCFromConversion ; //! Invariant mass of 2 clusters originated in conversions + TH2F * fhReMCFromNotConversion ; //! Invariant mass of 2 clusters not originated in conversions + TH2F * fhReMCFromMixConversion ; //! Invariant mass of 2 clusters one from conversion and the other not + + + + ClassDef(AliAnaPi0,21) } ; -- 2.43.5