From 03f5696bd10312e86654be612458bd9492c4096c Mon Sep 17 00:00:00 2001 From: miweber Date: Fri, 30 Aug 2013 15:47:46 +0000 Subject: [PATCH] adding mass squared to conversion cut control histograms (MW) --- PWGCF/EBYE/BalanceFunctions/AliBalancePsi.cxx | 9 +++++---- PWGCF/EBYE/BalanceFunctions/AliBalancePsi.h | 8 ++++---- 2 files changed, 9 insertions(+), 8 deletions(-) diff --git a/PWGCF/EBYE/BalanceFunctions/AliBalancePsi.cxx b/PWGCF/EBYE/BalanceFunctions/AliBalancePsi.cxx index 1695d852127..1c9b7199677 100644 --- a/PWGCF/EBYE/BalanceFunctions/AliBalancePsi.cxx +++ b/PWGCF/EBYE/BalanceFunctions/AliBalancePsi.cxx @@ -396,8 +396,8 @@ void AliBalancePsi::InitHistograms() { // QA histograms fHistHBTbefore = new TH2D("fHistHBTbefore","before HBT cut",200,0,2,200,0,2.*TMath::Pi()); fHistHBTafter = new TH2D("fHistHBTafter","after HBT cut",200,0,2,200,0,2.*TMath::Pi()); - fHistConversionbefore = new TH2D("fHistConversionbefore","before Conversion cut",200,0,2,200,0,2.*TMath::Pi()); - fHistConversionafter = new TH2D("fHistConversionafter","after Conversion cut",200,0,2,200,0,2.*TMath::Pi()); + fHistConversionbefore = new TH3D("fHistConversionbefore","before Conversion cut;#Delta#eta;#Delta#phi;M_{inv}^{2}",50,-2.0,2.0,50,-TMath::Pi()/2.,3.*TMath::Pi()/2.,300,0,1.5); + fHistConversionafter = new TH3D("fHistConversionafter","after Conversion cut;#Delta#eta;#Delta#phi;M_{inv}^{2}",50,-2.0,2.0,50,-TMath::Pi()/2.,3.*TMath::Pi()/2.,300,0,1.5); fHistPsiMinusPhi = new TH2D("fHistPsiMinusPhi","",4,-0.5,3.5,100,0,2.*TMath::Pi()); fHistResonancesBefore = new TH3D("fHistResonancesBefore","before resonance cut;#Delta#eta;#Delta#phi;M_{inv}",50,-2.0,2.0,50,-TMath::Pi()/2.,3.*TMath::Pi()/2.,300,0,1.5); fHistResonancesRho = new TH3D("fHistResonancesRho","after #rho resonance cut;#Delta#eta;#Delta#phi;M_{inv}",50,-2.0,2.0,50,-TMath::Pi()/2.,3.*TMath::Pi()/2.,300,0,1.5); @@ -640,7 +640,6 @@ void AliBalancePsi::CalculateBalance(Double_t gReactionPlane, if (charge1 * charge2 < 0) { Double_t deta = firstEta - secondEta[j]; Double_t dphi = firstPhi - secondPhi[j]; - fHistConversionbefore->Fill(deta,dphi); Float_t m0 = 0.510e-3; Float_t tantheta1 = 1e10; @@ -662,12 +661,14 @@ void AliBalancePsi::CalculateBalance(Double_t gReactionPlane, Float_t e2squ = m0 * m0 + secondPt[j] * secondPt[j] * (1.0 + 1.0 / tantheta2 / tantheta2); Float_t masssqu = 2 * m0 * m0 + 2 * ( TMath::Sqrt(e1squ * e2squ) - ( firstPt * secondPt[j] * ( TMath::Cos(phi1rad - phi2rad) + 1.0 / tantheta1 / tantheta2 ) ) ); + + fHistConversionbefore->Fill(deta,dphi,masssqu); if (masssqu < fInvMassCutConversion*fInvMassCutConversion){ //AliInfo(Form("Conversion: Removed track pair %d %d with [[%f %f] %f %f] %d %d <- %f %f %f %f %f %f ", i, j, deta, dphi, masssqu, charge1, charge2,eta1,eta2,phi1,phi2,pt1,pt2)); continue; } - fHistConversionafter->Fill(deta,dphi); + fHistConversionafter->Fill(deta,dphi,masssqu); } }//conversion cut diff --git a/PWGCF/EBYE/BalanceFunctions/AliBalancePsi.h b/PWGCF/EBYE/BalanceFunctions/AliBalancePsi.h index 79883ccda91..27cb676a0c2 100644 --- a/PWGCF/EBYE/BalanceFunctions/AliBalancePsi.h +++ b/PWGCF/EBYE/BalanceFunctions/AliBalancePsi.h @@ -193,8 +193,8 @@ class AliBalancePsi : public TObject { TH2D *GetQAHistHBTbefore() {return fHistHBTbefore;} TH2D *GetQAHistHBTafter() {return fHistHBTafter;} - TH2D *GetQAHistConversionbefore() {return fHistConversionbefore;} - TH2D *GetQAHistConversionafter() {return fHistConversionafter;} + TH3D *GetQAHistConversionbefore() {return fHistConversionbefore;} + TH3D *GetQAHistConversionafter() {return fHistConversionafter;} TH2D *GetQAHistPsiMinusPhi() {return fHistPsiMinusPhi;} TH3D *GetQAHistResonancesBefore() {return fHistResonancesBefore;} TH3D *GetQAHistResonancesRho() {return fHistResonancesRho;} @@ -233,8 +233,8 @@ class AliBalancePsi : public TObject { //QA histograms TH2D *fHistHBTbefore; // Delta Eta vs. Delta Phi before HBT inspired cuts TH2D *fHistHBTafter; // Delta Eta vs. Delta Phi after HBT inspired cuts - TH2D *fHistConversionbefore; // Delta Eta vs. Delta Phi before Conversion cuts - TH2D *fHistConversionafter; // Delta Eta vs. Delta Phi before Conversion cuts + TH3D *fHistConversionbefore; // 3D histogram (Deta,Dphi,Invmass) before Conversion cuts + TH3D *fHistConversionafter; // 3D histogram (Deta,Dphi,Invmass) before Conversion cuts TH2D *fHistPsiMinusPhi;// psi - phi QA histogram TH3D *fHistResonancesBefore; // 3D histogram (Deta,Dphi,Invmass) before resonance cuts TH3D *fHistResonancesRho; // 3D histogram (Deta,Dphi,Invmass) after removing rho -- 2.39.3