add FE eta plot for E>0

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;

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());
+    }
   }
 
 }