From 79e44280735d6ec0749f994411f7b0df978453c7 Mon Sep 17 00:00:00 2001
From: Tamara Vazquez Schroeder <tamara.vazquez.schroeder@cern.ch>
Date: Fri, 3 Nov 2017 16:31:15 +0000
Subject: [PATCH] Merge branch 'JetTagDQA-btagDevel-21.0-2017-11-10' into
 '21.0'

Update PhysValBTag: b-tagging validation code for new DAOD_PHYSVAL

See merge request !5377

(cherry picked from commit dbe407645461ad236a94fe53adce4f0871e67f35 [formerly 9deb94a12640ec758f2b37a5476ac01250eb0d80])

fffb87d0 Update PhysValBTag: b-tagging validation code for DAOD_PHYSVAL

Former-commit-id: bd756f349b04ba571896c46373ef32e854b56a45
---
 .../JetTagDQA/src/PhysValBTag.cxx             | 29 ++++++++++---------
 .../JetTagDQA/src/PhysValBTag.h               |  4 +--
 2 files changed, 18 insertions(+), 15 deletions(-)

diff --git a/PhysicsAnalysis/JetTagging/JetTagValidation/JetTagDQA/src/PhysValBTag.cxx b/PhysicsAnalysis/JetTagging/JetTagValidation/JetTagDQA/src/PhysValBTag.cxx
index 16b1ac3aa17..185504733ed 100644
--- a/PhysicsAnalysis/JetTagging/JetTagValidation/JetTagDQA/src/PhysValBTag.cxx
+++ b/PhysicsAnalysis/JetTagging/JetTagValidation/JetTagDQA/src/PhysValBTag.cxx
@@ -41,21 +41,21 @@ namespace JetTagDQA {
     //   m_jetPlots(0, "Summary/Jet/", "Jet"),
     //    m_trkvtxPlots(0, "Summary/TrackAndVertex/"),
     m_isData(false),
-    m_antiKt2PV0TrackJetPlots(0, "BTag/AntiKt2PV0TrackJets/", "antiKt2PV0TrackJets"),
-    m_antiKt3PV0TrackJetPlots(0, "BTag/AntiKt3PV0TrackJets/", "antiKt3PV0TrackJets"),
-    m_antiKt4PV0TrackJetPlots(0, "BTag/AntiKt4PV0TrackJets/", "antiKt4PV0TrackJets"),
-    m_antiKt4EMTopoPlots     (0, "BTag/AntiKt4EMTopo/"      , "antiKt4EMTopo"),
-    m_antiKt4LCTopoPlots     (0, "BTag/AntiKt4LCTopo/"      , "antiKt4LCTopo"),
+    m_antiKt2PV0TrackJetPlots             (0, "BTag/AntiKt2PV0TrackJets/"          , "antiKt2PV0TrackJets"),
+    m_antiKt4PV0TrackJetPlots             (0, "BTag/AntiKt4PV0TrackJets/"          , "antiKt4PV0TrackJets"),
+    m_antiKt4EMTopoPlots                  (0, "BTag/AntiKt4EMTopo/"                , "antiKt4EMTopo"),
+    m_antiKtVR30Rmax4Rmin02TrackJetsPlots (0, "BTag/AntiKtVR30Rmax4Rmin02TrackJets", "antiKtVR30Rmax4Rmin02TrackJets"),
+    m_antiKt4EMPFlowJetsPlots             (0, "BTag/AntiKt4EMPFlowJets"            , "antiKtVR30Rmax4Rmin02TrackJets"),
     m_nevents(0)
   {
 
     declareProperty( "isData", m_isData );
 
     declareProperty( "JetContainerName1", m_jetName1 = "AntiKt2PV0TrackJets");
-    declareProperty( "JetContainerName2", m_jetName2 = "AntiKt3PV0TrackJets");
-    declareProperty( "JetContainerName3", m_jetName3 = "AntiKt4PV0TrackJets");
-    declareProperty( "JetContainerName4", m_jetName4 = "AntiKt4EMTopoJets" );
-    declareProperty( "JetContainerName5", m_jetName5 = "AntiKt4LCTopoJets");
+    declareProperty( "JetContainerName2", m_jetName2 = "AntiKt4PV0TrackJets");
+    declareProperty( "JetContainerName3", m_jetName3 = "AntiKt4EMTopoJets" );
+    declareProperty( "JetContainerName4", m_jetName4 = "AntiKtVR30Rmax4Rmin02TrackJets");
+    declareProperty( "JetContainerName5", m_jetName5 = "AntiKt4EMPFlowJets");
 
     declareProperty( "TrackContainerName", m_trackName = "InDetTrackParticles" );
     declareProperty( "VertexContainerName", m_vertexName = "PrimaryVertices" );
@@ -74,10 +74,10 @@ namespace JetTagDQA {
 	ATH_CHECK(ManagedMonitorToolBase::initialize());
 
 	m_btagplots.insert(std::make_pair(m_jetName1, m_antiKt2PV0TrackJetPlots));
-	m_btagplots.insert(std::make_pair(m_jetName2, m_antiKt3PV0TrackJetPlots));
-	m_btagplots.insert(std::make_pair(m_jetName3, m_antiKt4PV0TrackJetPlots));
-	m_btagplots.insert(std::make_pair(m_jetName4, m_antiKt4EMTopoPlots));
-	m_btagplots.insert(std::make_pair(m_jetName5, m_antiKt4LCTopoPlots));
+	m_btagplots.insert(std::make_pair(m_jetName2, m_antiKt4PV0TrackJetPlots));
+	m_btagplots.insert(std::make_pair(m_jetName3, m_antiKt4EMTopoPlots));
+	m_btagplots.insert(std::make_pair(m_jetName4, m_antiKtVR30Rmax4Rmin02TrackJetsPlots));
+	m_btagplots.insert(std::make_pair(m_jetName5, m_antiKt4EMPFlowJetsPlots));
    
         return StatusCode::SUCCESS;
   }
@@ -171,6 +171,9 @@ namespace JetTagDQA {
         int label(1000);
         //double dR(1000);
 
+	//Arnaud: JVT cut to remove horns
+	//if(jet->getAttribute<float>("Jvt") < 0.59 && jet->pt() < 60e3 && std::abs(jet->eta()) < 2.4) continue;
+
         if(jet->pt() > 20000 && std::abs(jet->eta()) < 2.5){
           (plot_i->second).fill(jet);
 
diff --git a/PhysicsAnalysis/JetTagging/JetTagValidation/JetTagDQA/src/PhysValBTag.h b/PhysicsAnalysis/JetTagging/JetTagValidation/JetTagDQA/src/PhysValBTag.h
index 4d5b217679f..1f5027f2a42 100644
--- a/PhysicsAnalysis/JetTagging/JetTagValidation/JetTagDQA/src/PhysValBTag.h
+++ b/PhysicsAnalysis/JetTagging/JetTagValidation/JetTagDQA/src/PhysValBTag.h
@@ -93,10 +93,10 @@ namespace JetTagDQA {
     std::map<std::string, JetTagDQA::BTaggingValidationPlots> m_btagplots;
 
     JetTagDQA::BTaggingValidationPlots m_antiKt2PV0TrackJetPlots;
-    JetTagDQA::BTaggingValidationPlots m_antiKt3PV0TrackJetPlots;
     JetTagDQA::BTaggingValidationPlots m_antiKt4PV0TrackJetPlots;
     JetTagDQA::BTaggingValidationPlots m_antiKt4EMTopoPlots;
-    JetTagDQA::BTaggingValidationPlots m_antiKt4LCTopoPlots;
+    JetTagDQA::BTaggingValidationPlots m_antiKtVR30Rmax4Rmin02TrackJetsPlots;
+    JetTagDQA::BTaggingValidationPlots m_antiKt4EMPFlowJetsPlots;
 
     int m_nevents;
     //int m_nTruthB;
-- 
GitLab