From 23c982671a2163e880ba2ce1c75ca71cc0c3348c Mon Sep 17 00:00:00 2001
From: Matthew Thomas Anthony <matthew.thomas.anthony@cern.ch>
Date: Fri, 11 Dec 2020 10:45:25 +0000
Subject: [PATCH] add FE eta plot for E>0
---
.../PFlowValidation/PFOHistUtils/PFOHistUtils/PFOPlots.h | 3 +++
.../PFlow/PFlowValidation/PFOHistUtils/src/PFOPlots.cxx | 6 +++++-
2 files changed, 8 insertions(+), 1 deletion(-)
diff --git a/Reconstruction/PFlow/PFlowValidation/PFOHistUtils/PFOHistUtils/PFOPlots.h b/Reconstruction/PFlow/PFlowValidation/PFOHistUtils/PFOHistUtils/PFOPlots.h
index 3bee6b4aacc2..3b08e3b139fd 100644
--- a/Reconstruction/PFlow/PFlowValidation/PFOHistUtils/PFOHistUtils/PFOPlots.h
+++ b/Reconstruction/PFlow/PFlowValidation/PFOHistUtils/PFOHistUtils/PFOPlots.h
@@ -25,6 +25,7 @@ namespace PFO {
TH1* m_PFO_phi;
TH1* m_PFO_m;
TH1* m_PFO_charge;
+
/** Pt Histogram binned in eta */
TH1* m_PFO_pt_etaBinA;
@@ -37,6 +38,8 @@ namespace PFO {
TH1* m_FE_phi;
TH1* m_FE_m;
TH1* m_FE_charge;
+
+ TH1* m_FE_eta_posE; // extra bin for positive energy FE
/** Pt Histogram binned in eta */
TH1* m_FE_pt_etaBinA;
diff --git a/Reconstruction/PFlow/PFlowValidation/PFOHistUtils/src/PFOPlots.cxx b/Reconstruction/PFlow/PFlowValidation/PFOHistUtils/src/PFOPlots.cxx
index a4c450652f8b..d1c941b7ebc6 100644
--- a/Reconstruction/PFlow/PFlowValidation/PFOHistUtils/src/PFOPlots.cxx
+++ b/Reconstruction/PFlow/PFlowValidation/PFOHistUtils/src/PFOPlots.cxx
@@ -38,6 +38,7 @@ namespace PFO {
m_FE_pt_etaBinA = Book1D("_Pt_A",m_sFEContainerName + "_Pt (Entries/1 GeV, |eta| < 1 )",30,-10.0,20.0);
m_FE_pt_etaBinB = Book1D("_Pt_B",m_sFEContainerName + "_Pt (Entries/1 GeV, 1 <= |eta| < 2 )",30,-10.0,20.0);
m_FE_pt_etaBinC = Book1D("_Pt_C",m_sFEContainerName + "_Pt (Entries/1 GeV, |eta| >= 2 )",30,-10.0,20.0);
+ m_FE_eta_posE=Book1D("_EtaPosE",m_sFEContainerName+"_Eta (E>0) (Entries/0.1)",100,-5.0,5.0);
}
}
@@ -64,7 +65,10 @@ namespace PFO {
if (fabs(FE.eta()) < 1) m_FE_pt_etaBinA->Fill(FE.pt()/1000.0);
else if (fabs(FE.eta()) < 2) m_FE_pt_etaBinB->Fill(FE.pt()/1000.0);
else m_FE_pt_etaBinC->Fill(FE.pt()/1000.0);
-
+ // additional debug plot: Eta given FlowElem energy>0
+ if (FE.e()>0){
+ m_FE_eta_posE->Fill(FE.eta());
+ }
}
}
--
GitLab