#include <SkimmingToolHIGG5VBF.h>

Inheritance diagram for DerivationFramework::SkimmingToolHIGG5VBF:

Collaboration diagram for DerivationFramework::SkimmingToolHIGG5VBF:

Public Member Functions
	SkimmingToolHIGG5VBF (const std::string &t, const std::string &n, const IInterface *p)
	Constructor with parameters. More...

	~SkimmingToolHIGG5VBF ()
	Destructor. More...

virtual StatusCode	initialize () override

virtual StatusCode	finalize () override

virtual bool	eventPassesFilter () const override
	Check that the current event passes this filter. More...

Private Member Functions
bool	checkAllJetQuality (const TLorentzVector &jet) const

bool	checkCentralJetQuality (const TLorentzVector &jet) const

TLorentzVector	getCalibedJets (const xAOD::Jet *jet) const

Private Attributes
bool	m_debug

ToolHandle< Trig::TrigDecisionTool >	m_trigDecisionTool

std::atomic< unsigned int >	m_ntot

std::atomic< unsigned int >	m_npass

std::string	m_jetSGKey

std::string	m_calibedJetMomKey

bool	m_reqNAllJets

unsigned int	m_nAllJets

double	m_allJetPtCut

double	m_allJetEtaCut

bool	m_reqNCentralJets

unsigned int	m_nCentralJets

double	m_centralJetPtCut

double	m_centralJetEtaCut

bool	m_reqTrigger

std::vector< std::string >	m_triggers

bool	m_reqVbfMjj

double	m_vbfMjjCut

std::string	m_phSGKey

bool	m_reqPh

double	m_phPtCut

double	m_centralPhEtaCut

Detailed Description

Author: Yasuy.nosp@m.uki..nosp@m.Okumu.nosp@m.ra@c.nosp@m.ern.c.nosp@m.h

Definition at line 34 of file SkimmingToolHIGG5VBF.h.

Constructor & Destructor Documentation

◆ SkimmingToolHIGG5VBF()

DerivationFramework::SkimmingToolHIGG5VBF::SkimmingToolHIGG5VBF	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p
	)

Constructor with parameters.

Definition at line 19 of file SkimmingToolHIGG5VBF.cxx.

                                                      : 
   base_class(t, n, p),
   m_trigDecisionTool("Trig::TrigDecisionTool/TrigDecisionTool"),
   m_ntot(0),
   m_npass(0) 
 {
   
   // for jet multiplicity requirement
   declareProperty("JetContainerKey", m_jetSGKey="AntiKt4EMTopoJets");
   declareProperty("CalibedJetMomentKey", m_calibedJetMomKey="DFCommonJets_Calib");
   
   // for jet multiplicity
   declareProperty("ReqNAllJets", m_reqNAllJets=false);
   declareProperty("NumberOfAllJets", m_nAllJets=4);
   declareProperty("AllJetPtCut", m_allJetPtCut=40.*CLHEP::GeV); 
   declareProperty("AllJetEtaCut", m_allJetEtaCut=4.9); 
   
   declareProperty("ReqNCentralJets", m_reqNCentralJets=false);
   declareProperty("NumberOfCentralJets", m_nCentralJets=2);
   declareProperty("CentralJetPtCut", m_centralJetPtCut=40.*CLHEP::GeV); 
   declareProperty("CentralJetEtaCut", m_centralJetEtaCut=2.6); 
   
   // for trigger requirement
   declareProperty("ReqTrigger", m_reqTrigger=false); 
   declareProperty("Triggers", m_triggers=std::vector<std::string>()); 
   
   // for Mjj requirement ()
   declareProperty("ReqVBFMjj", m_reqVbfMjj=false); // logical "ORed" with ReqPhoton
   declareProperty("MjjCut", m_vbfMjjCut=0.*CLHEP::GeV);
   
   declareProperty("DoDebug", m_debug=false);
  
   // photon requirement (p. rose)
   declareProperty("PhotonContainerKey", m_phSGKey="Photons");
   declareProperty("ReqPhoton", m_reqPh=false); // logical "ORed" with ReqVBFMjj
   declareProperty("PhotonPtCut", m_phPtCut=0.*CLHEP::GeV);
   declareProperty("CentralPhotonEtaCut", m_centralPhEtaCut=2.6);
   
 }

◆ ~SkimmingToolHIGG5VBF()

DerivationFramework::SkimmingToolHIGG5VBF::~SkimmingToolHIGG5VBF ( )

Destructor.

Definition at line 62 of file SkimmingToolHIGG5VBF.cxx.

62 {

63 }

Member Function Documentation

◆ checkAllJetQuality()

bool DerivationFramework::SkimmingToolHIGG5VBF::checkAllJetQuality ( const TLorentzVector & jet ) const

private

Definition at line 176 of file SkimmingToolHIGG5VBF.cxx.

 {
   if(jet.Pt()<m_allJetPtCut) return false;
   if(fabs(jet.Eta())>m_allJetEtaCut) return false;
   
   return true;
 }

◆ checkCentralJetQuality()

bool DerivationFramework::SkimmingToolHIGG5VBF::checkCentralJetQuality ( const TLorentzVector & jet ) const

private

Definition at line 184 of file SkimmingToolHIGG5VBF.cxx.

 {
   if(jet.Pt()<m_centralJetPtCut) return false;
   if(fabs(jet.Eta())>m_centralJetEtaCut) return false;
   
   return true;
 }

◆ eventPassesFilter()

bool DerivationFramework::SkimmingToolHIGG5VBF::eventPassesFilter ( ) const

overridevirtual

Check that the current event passes this filter.

Definition at line 88 of file SkimmingToolHIGG5VBF.cxx.

 {
   m_ntot++;
   bool acceptEvent(true);
   
   StatusCode sc(StatusCode::SUCCESS);
   
   bool isTriggerFired(m_triggers.empty());
   for(unsigned int i(0); i<m_triggers.size(); i++) {
     if(m_trigDecisionTool->isPassed(m_triggers.at(i))) {
       isTriggerFired = true;
       break;
     }
   }
   
   // (1) Count Jet Multiplicity
   std::vector<TLorentzVector> goodAllJets;
   std::vector<TLorentzVector> goodCentralJets;
  
   const xAOD::JetContainer *jets(nullptr); 
   ATH_CHECK(evtStore()->retrieve(jets, m_jetSGKey), false);
   xAOD::JetContainer::const_iterator jet_itr(jets->begin());
   xAOD::JetContainer::const_iterator jet_end(jets->end());
   for(; jet_itr != jet_end; ++jet_itr) {
     TLorentzVector jetP4 = getCalibedJets( (*jet_itr) );
     if(this->checkAllJetQuality(jetP4)) { goodAllJets.push_back(jetP4); }
     if(this->checkCentralJetQuality(jetP4)) { goodCentralJets.push_back(jetP4); }
   }
   
   // (2) evaluate maximum Mjj in the event
   double maxM = 0.;
   for(unsigned int jet_i = 0; jet_i<goodAllJets.size(); jet_i++) {
     const TLorentzVector& iP4 = goodAllJets.at(jet_i);
     
     for(unsigned int jet_k=jet_i+1; jet_k<goodAllJets.size(); jet_k++) {
       const TLorentzVector& kP4 = goodAllJets.at(jet_k);
       
       const TLorentzVector jjP4 = iP4 + kP4;
       const double jjM = jjP4.M();
       
       if (maxM<jjM) {maxM=jjM;}
     }    
   }
   
   //get max pt of any photon in the event (p. rose)
   double maxPhPt=0.;
   if(m_reqPh){
     const xAOD::PhotonContainer *phots(nullptr);
     ATH_CHECK(evtStore()->retrieve(phots, m_phSGKey), false);
     for (const auto* ph : *phots){
       if(abs(ph->eta())<m_centralPhEtaCut)
     if(ph->pt()>maxPhPt) maxPhPt = ph->pt();
     }//for
   }//if
   
  
  
   const bool passNAllJet     = (goodAllJets.size()>=m_nAllJets);
   const bool passNCentralJet = (goodCentralJets.size()>=m_nCentralJets);
   const bool passMjjCut      = (maxM>m_vbfMjjCut);
   const bool passPhPtCut     = (maxPhPt>m_phPtCut);
  
   if (m_reqNAllJets)     { if (not passNAllJet)     {acceptEvent=false;} }
   if (m_reqNCentralJets) { if (not passNCentralJet) {acceptEvent=false;} }
   //if (m_reqVbfMjj)       { if (not passMjjCut)      {acceptEvent=false;} }
   if (m_reqTrigger)      { if (not isTriggerFired)  {acceptEvent=false;} } 
   //vbf+gamma addition -- logical OR of mjj and phpt cut
   if(m_reqVbfMjj and not m_reqPh)     { if (not passMjjCut)                   {acceptEvent=false;} }
   if(m_reqPh     and not m_reqVbfMjj) { if (not passPhPtCut)                  {acceptEvent=false;} }
   if(m_reqPh     and     m_reqVbfMjj) { if (not (passMjjCut or passPhPtCut) ) {acceptEvent=false;} }
   
   if (acceptEvent) {m_npass++;}
   
   if (m_debug) {
     printf("dbg> L%3d : event accepted [%s] " 
        "(NJets=%2d [%6s], NCentralJets=%2d [%6s] Mjj=%10.1f [%6s], Trigger [%6s]) \n", 
        __LINE__, 
        acceptEvent? "Y" : "N",
        (int)goodAllJets.size(),  passNAllJet ? "PASSED" : "FAILED",
        (int)goodCentralJets.size(),  passNCentralJet ? "PASSED" : "FAILED",
        maxM, passMjjCut ? "PASSED" : "FAILED",
        isTriggerFired ? "PASSED" : "FAILED"
        );
   }
   
   return acceptEvent;
 }

◆ finalize()

StatusCode DerivationFramework::SkimmingToolHIGG5VBF::finalize ( )

overridevirtual

Definition at line 80 of file SkimmingToolHIGG5VBF.cxx.

 {
   ATH_MSG_VERBOSE("finalize() ...");
   ATH_MSG_INFO("Processed " << m_ntot << " events, " << m_npass << " events passed filter ");
   return StatusCode::SUCCESS;
 }

◆ getCalibedJets()

TLorentzVector DerivationFramework::SkimmingToolHIGG5VBF::getCalibedJets ( const xAOD::Jet * jet ) const

private

Definition at line 193 of file SkimmingToolHIGG5VBF.cxx.

 {
   TLorentzVector rc;
   
   if(!jet) return rc;
  
   SG::ConstAccessor<float> ptAcc(m_calibedJetMomKey+"_pt");
   SG::ConstAccessor<float> etaAcc(m_calibedJetMomKey+"_eta");
   SG::ConstAccessor<float> phiAcc(m_calibedJetMomKey+"_phi");
   SG::ConstAccessor<float> mAcc(m_calibedJetMomKey+"_m");
   const float& pt =ptAcc(*jet);
   const float& eta=etaAcc(*jet);
   const float& phi=phiAcc(*jet);
   const float& m  =mAcc(*jet);
   
   rc.SetPtEtaPhiM(pt, eta, phi, m);
   
   return rc;
 }

◆ initialize()

StatusCode DerivationFramework::SkimmingToolHIGG5VBF::initialize ( )

overridevirtual

Definition at line 66 of file SkimmingToolHIGG5VBF.cxx.

 {
   ATH_MSG_VERBOSE("initialize() ...");
  
   // trigger decision tool
   if(m_trigDecisionTool.retrieve(DisableTool{!m_reqTrigger}).isFailure()) {
     ATH_MSG_FATAL("Failed to retrieve tool: " << m_trigDecisionTool);
     return StatusCode::FAILURE;
   }
   ATH_MSG_INFO("Retrieved tool: " << m_trigDecisionTool);
  
   return StatusCode::SUCCESS;
 }

Member Data Documentation

◆ m_allJetEtaCut

double DerivationFramework::SkimmingToolHIGG5VBF::m_allJetEtaCut

private

Definition at line 65 of file SkimmingToolHIGG5VBF.h.

◆ m_allJetPtCut

double DerivationFramework::SkimmingToolHIGG5VBF::m_allJetPtCut

private

Definition at line 64 of file SkimmingToolHIGG5VBF.h.

◆ m_calibedJetMomKey

std::string DerivationFramework::SkimmingToolHIGG5VBF::m_calibedJetMomKey

private

Definition at line 59 of file SkimmingToolHIGG5VBF.h.

◆ m_centralJetEtaCut

double DerivationFramework::SkimmingToolHIGG5VBF::m_centralJetEtaCut

private

Definition at line 70 of file SkimmingToolHIGG5VBF.h.

◆ m_centralJetPtCut

double DerivationFramework::SkimmingToolHIGG5VBF::m_centralJetPtCut

private

Definition at line 69 of file SkimmingToolHIGG5VBF.h.

◆ m_centralPhEtaCut

double DerivationFramework::SkimmingToolHIGG5VBF::m_centralPhEtaCut

private

Definition at line 88 of file SkimmingToolHIGG5VBF.h.

◆ m_debug

bool DerivationFramework::SkimmingToolHIGG5VBF::m_debug

private

Definition at line 51 of file SkimmingToolHIGG5VBF.h.

◆ m_jetSGKey

std::string DerivationFramework::SkimmingToolHIGG5VBF::m_jetSGKey

private

Definition at line 58 of file SkimmingToolHIGG5VBF.h.

◆ m_nAllJets

unsigned int DerivationFramework::SkimmingToolHIGG5VBF::m_nAllJets

private

Definition at line 63 of file SkimmingToolHIGG5VBF.h.

◆ m_nCentralJets

unsigned int DerivationFramework::SkimmingToolHIGG5VBF::m_nCentralJets

private

Definition at line 68 of file SkimmingToolHIGG5VBF.h.

◆ m_npass

std::atomic<unsigned int> DerivationFramework::SkimmingToolHIGG5VBF::m_npass

mutableprivate

Definition at line 56 of file SkimmingToolHIGG5VBF.h.

◆ m_ntot

std::atomic<unsigned int> DerivationFramework::SkimmingToolHIGG5VBF::m_ntot

mutableprivate

Definition at line 55 of file SkimmingToolHIGG5VBF.h.

◆ m_phPtCut

double DerivationFramework::SkimmingToolHIGG5VBF::m_phPtCut

private

Definition at line 87 of file SkimmingToolHIGG5VBF.h.

◆ m_phSGKey

std::string DerivationFramework::SkimmingToolHIGG5VBF::m_phSGKey

private

Definition at line 85 of file SkimmingToolHIGG5VBF.h.

◆ m_reqNAllJets

bool DerivationFramework::SkimmingToolHIGG5VBF::m_reqNAllJets

private

Definition at line 62 of file SkimmingToolHIGG5VBF.h.

◆ m_reqNCentralJets

bool DerivationFramework::SkimmingToolHIGG5VBF::m_reqNCentralJets

private

Definition at line 67 of file SkimmingToolHIGG5VBF.h.

◆ m_reqPh

bool DerivationFramework::SkimmingToolHIGG5VBF::m_reqPh

private

Definition at line 86 of file SkimmingToolHIGG5VBF.h.

◆ m_reqTrigger

bool DerivationFramework::SkimmingToolHIGG5VBF::m_reqTrigger

private

Definition at line 73 of file SkimmingToolHIGG5VBF.h.

◆ m_reqVbfMjj

bool DerivationFramework::SkimmingToolHIGG5VBF::m_reqVbfMjj

private

Definition at line 77 of file SkimmingToolHIGG5VBF.h.

◆ m_trigDecisionTool

ToolHandle<Trig::TrigDecisionTool> DerivationFramework::SkimmingToolHIGG5VBF::m_trigDecisionTool

private

Definition at line 53 of file SkimmingToolHIGG5VBF.h.

◆ m_triggers

std::vector<std::string> DerivationFramework::SkimmingToolHIGG5VBF::m_triggers

private

Definition at line 74 of file SkimmingToolHIGG5VBF.h.

◆ m_vbfMjjCut

double DerivationFramework::SkimmingToolHIGG5VBF::m_vbfMjjCut

private

Definition at line 78 of file SkimmingToolHIGG5VBF.h.

The documentation for this class was generated from the following files:

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ SkimmingToolHIGG5VBF()

◆ ~SkimmingToolHIGG5VBF()

Member Function Documentation

◆ checkAllJetQuality()

◆ checkCentralJetQuality()

◆ eventPassesFilter()

◆ finalize()

◆ getCalibedJets()

◆ initialize()

Member Data Documentation

◆ m_allJetEtaCut

◆ m_allJetPtCut

◆ m_calibedJetMomKey

◆ m_centralJetEtaCut

◆ m_centralJetPtCut

◆ m_centralPhEtaCut

◆ m_debug

◆ m_jetSGKey

◆ m_nAllJets

◆ m_nCentralJets

◆ m_npass

◆ m_ntot

◆ m_phPtCut

◆ m_phSGKey

◆ m_reqNAllJets

◆ m_reqNCentralJets

◆ m_reqPh

◆ m_reqTrigger

◆ m_reqVbfMjj

◆ m_trigDecisionTool

◆ m_triggers

◆ m_vbfMjjCut