Skip to content
Snippets Groups Projects
Commit dbe40764 authored by Tamara Vazquez Schroeder's avatar Tamara Vazquez Schroeder
Browse files

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

Former-commit-id: 9deb94a12640ec758f2b37a5476ac01250eb0d80
parents 6f17e257 dfe677d1
No related branches found
No related tags found
No related merge requests found
...@@ -41,21 +41,21 @@ namespace JetTagDQA { ...@@ -41,21 +41,21 @@ namespace JetTagDQA {
// m_jetPlots(0, "Summary/Jet/", "Jet"), // m_jetPlots(0, "Summary/Jet/", "Jet"),
// m_trkvtxPlots(0, "Summary/TrackAndVertex/"), // m_trkvtxPlots(0, "Summary/TrackAndVertex/"),
m_isData(false), m_isData(false),
m_antiKt2PV0TrackJetPlots(0, "BTag/AntiKt2PV0TrackJets/", "antiKt2PV0TrackJets"), m_antiKt2PV0TrackJetPlots (0, "BTag/AntiKt2PV0TrackJets/" , "antiKt2PV0TrackJets"),
m_antiKt3PV0TrackJetPlots(0, "BTag/AntiKt3PV0TrackJets/", "antiKt3PV0TrackJets"), m_antiKt4PV0TrackJetPlots (0, "BTag/AntiKt4PV0TrackJets/" , "antiKt4PV0TrackJets"),
m_antiKt4PV0TrackJetPlots(0, "BTag/AntiKt4PV0TrackJets/", "antiKt4PV0TrackJets"), m_antiKt4EMTopoPlots (0, "BTag/AntiKt4EMTopo/" , "antiKt4EMTopo"),
m_antiKt4EMTopoPlots (0, "BTag/AntiKt4EMTopo/" , "antiKt4EMTopo"), m_antiKtVR30Rmax4Rmin02TrackJetsPlots (0, "BTag/AntiKtVR30Rmax4Rmin02TrackJets", "antiKtVR30Rmax4Rmin02TrackJets"),
m_antiKt4LCTopoPlots (0, "BTag/AntiKt4LCTopo/" , "antiKt4LCTopo"), m_antiKt4EMPFlowJetsPlots (0, "BTag/AntiKt4EMPFlowJets" , "antiKtVR30Rmax4Rmin02TrackJets"),
m_nevents(0) m_nevents(0)
{ {
declareProperty( "isData", m_isData ); declareProperty( "isData", m_isData );
declareProperty( "JetContainerName1", m_jetName1 = "AntiKt2PV0TrackJets"); declareProperty( "JetContainerName1", m_jetName1 = "AntiKt2PV0TrackJets");
declareProperty( "JetContainerName2", m_jetName2 = "AntiKt3PV0TrackJets"); declareProperty( "JetContainerName2", m_jetName2 = "AntiKt4PV0TrackJets");
declareProperty( "JetContainerName3", m_jetName3 = "AntiKt4PV0TrackJets"); declareProperty( "JetContainerName3", m_jetName3 = "AntiKt4EMTopoJets" );
declareProperty( "JetContainerName4", m_jetName4 = "AntiKt4EMTopoJets" ); declareProperty( "JetContainerName4", m_jetName4 = "AntiKtVR30Rmax4Rmin02TrackJets");
declareProperty( "JetContainerName5", m_jetName5 = "AntiKt4LCTopoJets"); declareProperty( "JetContainerName5", m_jetName5 = "AntiKt4EMPFlowJets");
declareProperty( "TrackContainerName", m_trackName = "InDetTrackParticles" ); declareProperty( "TrackContainerName", m_trackName = "InDetTrackParticles" );
declareProperty( "VertexContainerName", m_vertexName = "PrimaryVertices" ); declareProperty( "VertexContainerName", m_vertexName = "PrimaryVertices" );
...@@ -74,10 +74,10 @@ namespace JetTagDQA { ...@@ -74,10 +74,10 @@ namespace JetTagDQA {
ATH_CHECK(ManagedMonitorToolBase::initialize()); ATH_CHECK(ManagedMonitorToolBase::initialize());
m_btagplots.insert(std::make_pair(m_jetName1, m_antiKt2PV0TrackJetPlots)); 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_jetName2, m_antiKt4PV0TrackJetPlots));
m_btagplots.insert(std::make_pair(m_jetName3, m_antiKt4PV0TrackJetPlots)); m_btagplots.insert(std::make_pair(m_jetName3, m_antiKt4EMTopoPlots));
m_btagplots.insert(std::make_pair(m_jetName4, m_antiKt4EMTopoPlots)); m_btagplots.insert(std::make_pair(m_jetName4, m_antiKtVR30Rmax4Rmin02TrackJetsPlots));
m_btagplots.insert(std::make_pair(m_jetName5, m_antiKt4LCTopoPlots)); m_btagplots.insert(std::make_pair(m_jetName5, m_antiKt4EMPFlowJetsPlots));
return StatusCode::SUCCESS; return StatusCode::SUCCESS;
} }
...@@ -171,6 +171,9 @@ namespace JetTagDQA { ...@@ -171,6 +171,9 @@ namespace JetTagDQA {
int label(1000); int label(1000);
//double dR(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){ if(jet->pt() > 20000 && std::abs(jet->eta()) < 2.5){
(plot_i->second).fill(jet); (plot_i->second).fill(jet);
......
...@@ -93,10 +93,10 @@ namespace JetTagDQA { ...@@ -93,10 +93,10 @@ namespace JetTagDQA {
std::map<std::string, JetTagDQA::BTaggingValidationPlots> m_btagplots; std::map<std::string, JetTagDQA::BTaggingValidationPlots> m_btagplots;
JetTagDQA::BTaggingValidationPlots m_antiKt2PV0TrackJetPlots; JetTagDQA::BTaggingValidationPlots m_antiKt2PV0TrackJetPlots;
JetTagDQA::BTaggingValidationPlots m_antiKt3PV0TrackJetPlots;
JetTagDQA::BTaggingValidationPlots m_antiKt4PV0TrackJetPlots; JetTagDQA::BTaggingValidationPlots m_antiKt4PV0TrackJetPlots;
JetTagDQA::BTaggingValidationPlots m_antiKt4EMTopoPlots; JetTagDQA::BTaggingValidationPlots m_antiKt4EMTopoPlots;
JetTagDQA::BTaggingValidationPlots m_antiKt4LCTopoPlots; JetTagDQA::BTaggingValidationPlots m_antiKtVR30Rmax4Rmin02TrackJetsPlots;
JetTagDQA::BTaggingValidationPlots m_antiKt4EMPFlowJetsPlots;
int m_nevents; int m_nevents;
//int m_nTruthB; //int m_nTruthB;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment