DerivationFramework::SkimmingToolEXOT14 Class Reference

#include <SkimmingToolEXOT14.h>

Inheritance diagram for DerivationFramework::SkimmingToolEXOT14:

Collaboration diagram for DerivationFramework::SkimmingToolEXOT14:

Public Member Functions
	SkimmingToolEXOT14 (const std::string &t, const std::string &n, const IInterface *p)
	Constructor with parameters. More...
	~SkimmingToolEXOT14 ()
	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 Types
using	LeadingJets_t = std::array< TLorentzVector, 2 >

Private Member Functions
std::optional< LeadingJets_t >	SubcutPreselect () const
bool	SubcutGoodRunList () const
bool	SubcutLArError () const
bool	SubcutTrigger () const
bool	SubcutJetPts (const LeadingJets_t &jets) const
bool	SubcutJetDEta (const LeadingJets_t &jets) const
bool	SubcutDijetMass (const LeadingJets_t &jets) const
bool	SubcutJetDPhi (const LeadingJets_t &jets) const
std::string	TriggerVarName (std::string s) const

Private Attributes
ToolHandle< Trig::TrigDecisionTool >	m_trigDecisionTool
ToolHandle< IJetCalibrationTool >	m_JESTool
std::string	m_jetSGKey
bool	m_reqGRL
bool	m_reqLArError
bool	m_reqTrigger
bool	m_reqPreselection
bool	m_reqJetPts
bool	m_reqJetsDEta
bool	m_reqDiJetMass
bool	m_reqJetsDPhi
std::string	m_goodRunList
std::string	m_defaultTrigger
std::vector< std::string >	m_triggers
double	m_minJetPt
double	m_maxEta
double	m_leadingJetPt
double	m_subleadingJetPt
double	m_etaSeparation
double	m_dijetMass
double	m_jetDPhi
std::atomic< unsigned int >	m_n_tot
std::atomic< unsigned int >	m_n_passGRL
std::atomic< unsigned int >	m_n_passLArError
std::atomic< unsigned int >	m_n_passTrigger
std::atomic< unsigned int >	m_n_passPreselect
std::atomic< unsigned int >	m_n_passJetPts
std::atomic< unsigned int >	m_n_passJetsDEta
std::atomic< unsigned int >	m_n_passDiJetMass
std::atomic< unsigned int >	m_n_passJetsDPhi
std::atomic< unsigned int >	m_n_pass

Static Private Attributes
static const double	s_MZ

Detailed Description

Author

jsaxo.nosp@m.n@ce.nosp@m.rn.ch (adapted from Susumu Oda)

Definition at line 45 of file SkimmingToolEXOT14.h.

Member Typedef Documentation

LeadingJets_t

using DerivationFramework::SkimmingToolEXOT14::LeadingJets_t = std::array<TLorentzVector, 2>

private

Definition at line 103 of file SkimmingToolEXOT14.h.

Constructor & Destructor Documentation

SkimmingToolEXOT14()

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

Constructor with parameters.

Definition at line 23 of file SkimmingToolEXOT14.cxx.

                                                     : 
  base_class(t, n, p),
  m_trigDecisionTool("Trig::TrigDecisionTool/TrigDecisionTool"),
  m_n_tot(0),
  m_n_passGRL(0),
  m_n_passLArError(0),
  m_n_passTrigger(0),
  m_n_passPreselect(0),
  m_n_passJetPts(0),
  m_n_passJetsDEta(0),
  m_n_passDiJetMass(0),
  m_n_passJetsDPhi(0),
  m_n_pass(0)
{
 
 
  declareProperty("JetContainer",          m_jetSGKey = "AntiKt4LCTopoJets");
 
  declareProperty("RequireGRL",            m_reqGRL          = true);
  declareProperty("ReqireLArError",        m_reqLArError     = true);
  declareProperty("RequireTrigger",        m_reqTrigger      = true);
  declareProperty("RequirePreselection",   m_reqPreselection = true);
  declareProperty("RequireJetPts",         m_reqJetPts       = true);
  declareProperty("RequireJetsDEta",       m_reqJetsDEta     = true);
  declareProperty("RequireDiJetMass",      m_reqDiJetMass    = true);
  declareProperty("RequireJetsDPhi",       m_reqJetsDPhi     = true);
 
  declareProperty("GoodRunList",           m_goodRunList = "");
 
  declareProperty("DefaultTrigger",        m_defaultTrigger = "HLT_xe100");
  declareProperty("Triggers",              m_triggers = std::vector<std::string>()); 
 
  declareProperty("MinimumJetPt",          m_minJetPt = 25*CLHEP::GeV);
  declareProperty("MaxEta",                m_maxEta = 4.8);
 
  declareProperty("LeadingJetPtCut",       m_leadingJetPt    = 75*CLHEP::GeV);
  declareProperty("SubleadingJetPtCut",    m_subleadingJetPt = 50*CLHEP::GeV);
 
  declareProperty("EtaSeparation",         m_etaSeparation = 4.8);
 
  declareProperty("DiJetsMass",            m_dijetMass = 1000*CLHEP::GeV);
 
  declareProperty("DiJetDPhi",             m_jetDPhi = 2.5);
 
}

~SkimmingToolEXOT14()

DerivationFramework::SkimmingToolEXOT14::~SkimmingToolEXOT14

(

)

Destructor.

Definition at line 72 of file SkimmingToolEXOT14.cxx.

                                                           {
}  

Member Function Documentation

eventPassesFilter()

bool DerivationFramework::SkimmingToolEXOT14::eventPassesFilter ( ) const

overridevirtual

Check that the current event passes this filter.

Definition at line 121 of file SkimmingToolEXOT14.cxx.

{
 
  m_n_tot++;
 
  bool writeEvent(false);
 
  // int *leading    = new int(0);
  // if (!evtStore()->contains<int>("leading"))    CHECK(evtStore()->record(leading,    "leading"));
 
  if (!SubcutGoodRunList() && m_reqGRL      ) return false;
  if (!SubcutLArError()    && m_reqLArError ) return false;
  if (!SubcutTrigger()     && m_reqTrigger  ) return false;
 
  const auto jets = SubcutPreselect();
  if (!m_reqPreselection) writeEvent = true;        
 
  // There *must* be two jets for the remaining 
  // pieces, but you can still save the event...
  if (jets) {
 
    bool passDiJets(true);     
    if (!SubcutJetPts(jets.value())        && m_reqJetPts   ) passDiJets = false;
    if (!SubcutJetDEta(jets.value())       && m_reqJetsDEta ) passDiJets = false;
    if (!SubcutDijetMass(jets.value())     && m_reqDiJetMass) passDiJets = false;
    if (!SubcutJetDPhi(jets.value())       && m_reqJetsDPhi ) passDiJets = false;
    if (passDiJets) writeEvent = true; 
  }
 
  if (!writeEvent) return false;
  
  m_n_pass++;
  return true;
 
}

finalize()

StatusCode DerivationFramework::SkimmingToolEXOT14::finalize ( )

overridevirtual

Definition at line 102 of file SkimmingToolEXOT14.cxx.

{
  ATH_MSG_VERBOSE("finalize() ...");
  ATH_MSG_INFO("Processed " << m_n_tot << " events, " << m_n_pass << " events passed filter ");
 
  ATH_MSG_INFO("GRL       :: " << m_n_passGRL);
  ATH_MSG_INFO("LAr Error :: " << m_n_passLArError);
  ATH_MSG_INFO("Trigger   :: " << m_n_passTrigger);
  ATH_MSG_INFO("Preselect :: " << m_n_passPreselect);
  ATH_MSG_INFO("JetPts    :: " << m_n_passJetPts);
  ATH_MSG_INFO("JetsDEta  :: " << m_n_passJetsDEta);
  ATH_MSG_INFO("DijetMass :: " << m_n_passDiJetMass);
  ATH_MSG_INFO("JetsDPhi  :: " << m_n_passJetsDPhi);
 
  return StatusCode::SUCCESS;
 
}

initialize()

StatusCode DerivationFramework::SkimmingToolEXOT14::initialize ( )

overridevirtual

Definition at line 76 of file SkimmingToolEXOT14.cxx.

{
  ATH_MSG_VERBOSE("INITIALIZING VBFINV SELECTOR TOOL");
 
  // trigger decision tool
  if(m_trigDecisionTool.retrieve(DisableTool{!m_reqTrigger}).isFailure()) {
    ATH_MSG_FATAL("Failed to retrieve tool: " << m_trigDecisionTool);
    return StatusCode::FAILURE;
  }
  if (!m_triggers.size()) m_triggers.push_back(m_defaultTrigger);
  ATH_MSG_INFO("Retrieved tool: " << m_trigDecisionTool);
  
 
  // jet energy calibration
  m_JESTool.setTypeAndName("JetCalibrationTool/m_JESTool");
  CHECK( m_JESTool.retrieve() ); //optional, just forces initializing the tool here instead of at first use
  ATH_MSG_INFO("Retrieved tool: " << m_JESTool);
 
 
  return StatusCode::SUCCESS;
}

SubcutDijetMass()

bool DerivationFramework::SkimmingToolEXOT14::SubcutDijetMass ( const LeadingJets_t &  jets ) const

private

Definition at line 280 of file SkimmingToolEXOT14.cxx.

                                                                                           {
 
  const double DiJetMass = (jets[0] + jets[1]).M();
  const bool passDiJetMass = DiJetMass > m_dijetMass;
 
  if (passDiJetMass) m_n_passDiJetMass++;
  return passDiJetMass;
 
}

SubcutGoodRunList()

bool DerivationFramework::SkimmingToolEXOT14::SubcutGoodRunList ( ) const

private

Definition at line 157 of file SkimmingToolEXOT14.cxx.

                                                                    {
 
  // Placeholder
  m_n_passGRL++;
  return true;
 
}

SubcutJetDEta()

bool DerivationFramework::SkimmingToolEXOT14::SubcutJetDEta ( const LeadingJets_t &  jets ) const

private

Definition at line 269 of file SkimmingToolEXOT14.cxx.

                                                                                         {
 
  const double JetsDEta = fabs(jets[0].Eta() - jets[1].Eta());
  const bool passJetsDEta = JetsDEta > m_etaSeparation;
 
  if (passJetsDEta) m_n_passJetsDEta++;
  return passJetsDEta;
 
}

SubcutJetDPhi()

bool DerivationFramework::SkimmingToolEXOT14::SubcutJetDPhi ( const LeadingJets_t &  jets ) const

private

Definition at line 290 of file SkimmingToolEXOT14.cxx.

                                                                                         {
 
  const double JetsDPhi = fabs(jets[0].DeltaPhi(jets[1]));
  const bool passJetsDPhi = JetsDPhi < m_jetDPhi;
 
  if (passJetsDPhi) m_n_passJetsDPhi++;
  return passJetsDPhi;
 
}

SubcutJetPts()

bool DerivationFramework::SkimmingToolEXOT14::SubcutJetPts ( const LeadingJets_t &  jets ) const

private

Definition at line 259 of file SkimmingToolEXOT14.cxx.

                                                                                        {
 
  bool passJetPts =  (!m_leadingJetPt    || jets[0].Pt() > m_leadingJetPt);
  passJetPts &= (!m_subleadingJetPt || jets[1].Pt() > m_subleadingJetPt);
 
  if (passJetPts) m_n_passJetPts++;
  return passJetPts;
 
}

SubcutLArError()

bool DerivationFramework::SkimmingToolEXOT14::SubcutLArError ( ) const

private

Definition at line 166 of file SkimmingToolEXOT14.cxx.

                                                                 {
 
  // Retrieve EventInfo
  const xAOD::EventInfo *eventInfo(0);
  ATH_CHECK(evtStore()->retrieve(eventInfo), false);
 
  const bool passLArError = !(eventInfo->errorState(xAOD::EventInfo::LAr) == xAOD::EventInfo::Error);
  
  if (passLArError) m_n_passLArError++;
  return passLArError;
 
}

SubcutPreselect()

std::optional< DerivationFramework::SkimmingToolEXOT14::LeadingJets_t > DerivationFramework::SkimmingToolEXOT14::SubcutPreselect ( ) const

private

Definition at line 204 of file SkimmingToolEXOT14.cxx.

                                                             {
 
  // xAOD::TStore store;
  const xAOD::JetContainer *jets(0); 
  ATH_CHECK(evtStore()->retrieve(jets, m_jetSGKey), {});
  xAOD::JetContainer::const_iterator jet_itr(jets->begin());
  xAOD::JetContainer::const_iterator jet_end(jets->end());
 
  xAOD::JetContainer calibJets;
  calibJets.setStore(new xAOD::JetAuxContainer());
 
  TLorentzVector j1TLV, j2TLV;
 
  // Copy jets into the container to be calibrated
  while(jet_itr != jet_end) {
    xAOD::Jet* jetC = new xAOD::Jet();
    jetC->setJetP4(xAOD::JetFourMom_t((*jet_itr)->pt(), (*jet_itr)->eta(), (*jet_itr)->phi(), (*jet_itr)->m()));
    calibJets.push_back(jetC);
    ++jet_itr;
  }
 
  // Calibrate the jets
  if(m_JESTool->applyCalibration(calibJets).isFailure())
    ATH_MSG_WARNING("Jet calibration returned FAILURE!");
 
  jet_itr = calibJets.begin();
  jet_end = calibJets.end();
  while(jet_itr != jet_end) {
 
    if (abs((*jet_itr)->eta()) > m_maxEta) continue;
 
    if ((*jet_itr)->pt() > j1TLV.Pt()) {
 
      j2TLV = j1TLV;
      j1TLV.SetPtEtaPhiE((*jet_itr)->pt(), (*jet_itr)->eta(), (*jet_itr)->phi(), (*jet_itr)->e());
 
    } else if ((*jet_itr)->pt() > j2TLV.Pt()) {
 
      j2TLV.SetPtEtaPhiE((*jet_itr)->pt(), (*jet_itr)->eta(), (*jet_itr)->phi(), (*jet_itr)->e());
    }
    
    ++jet_itr;
  }
 
  // save this for this code.
  if (j2TLV.Pt() > m_minJetPt) {
    m_n_passPreselect++;
    return LeadingJets_t{j1TLV, j2TLV};
  }
 
  return {};
 
}

SubcutTrigger()

bool DerivationFramework::SkimmingToolEXOT14::SubcutTrigger ( ) const

private

Definition at line 180 of file SkimmingToolEXOT14.cxx.

                                                                {
 
  const xAOD::EventInfo *eventInfo(0);
  ATH_CHECK(evtStore()->retrieve(eventInfo), false);
 
  bool passTrigger = false;
 
  for (unsigned int i = 0; i < m_triggers.size(); i++) {
    bool thisTrig = m_trigDecisionTool->isPassed(m_triggers.at(i));
    SG::Decorator<bool> acc(TriggerVarName(m_triggers.at(i)));
    acc(*eventInfo) = thisTrig;
    // ATH_MSG_INFO("TRIGGER = " << m_triggers.at(i) <<  " -->> " << thisTrig);
    passTrigger |= thisTrig;
  }
  
  //  temporary pass-through of trigger cut for MC
  if (eventInfo->eventType(xAOD::EventInfo::IS_SIMULATION)) passTrigger = true;
 
  if (passTrigger) m_n_passTrigger++;
  return passTrigger;
 
}

TriggerVarName()

std::string DerivationFramework::SkimmingToolEXOT14::TriggerVarName ( std::string  s ) const

private

Definition at line 300 of file SkimmingToolEXOT14.cxx.

                                                                                 {
  std::replace(s.begin(), s.end(), '-', '_'); return s;
}

Member Data Documentation

m_defaultTrigger

std::string DerivationFramework::SkimmingToolEXOT14::m_defaultTrigger

private

Definition at line 88 of file SkimmingToolEXOT14.h.

m_dijetMass

double DerivationFramework::SkimmingToolEXOT14::m_dijetMass

private

Definition at line 96 of file SkimmingToolEXOT14.h.

m_etaSeparation

double DerivationFramework::SkimmingToolEXOT14::m_etaSeparation

private

Definition at line 95 of file SkimmingToolEXOT14.h.

m_goodRunList

std::string DerivationFramework::SkimmingToolEXOT14::m_goodRunList

private

Definition at line 86 of file SkimmingToolEXOT14.h.

m_JESTool

ToolHandle<IJetCalibrationTool> DerivationFramework::SkimmingToolEXOT14::m_JESTool

private

Definition at line 67 of file SkimmingToolEXOT14.h.

m_jetDPhi

double DerivationFramework::SkimmingToolEXOT14::m_jetDPhi

private

Definition at line 97 of file SkimmingToolEXOT14.h.

m_jetSGKey

std::string DerivationFramework::SkimmingToolEXOT14::m_jetSGKey

private

Definition at line 72 of file SkimmingToolEXOT14.h.

m_leadingJetPt

double DerivationFramework::SkimmingToolEXOT14::m_leadingJetPt

private

Definition at line 93 of file SkimmingToolEXOT14.h.

m_maxEta

double DerivationFramework::SkimmingToolEXOT14::m_maxEta

private

Definition at line 92 of file SkimmingToolEXOT14.h.

m_minJetPt

double DerivationFramework::SkimmingToolEXOT14::m_minJetPt

private

Definition at line 91 of file SkimmingToolEXOT14.h.

m_n_pass

std::atomic<unsigned int> DerivationFramework::SkimmingToolEXOT14::m_n_pass

mutableprivate

Definition at line 128 of file SkimmingToolEXOT14.h.

m_n_passDiJetMass

std::atomic<unsigned int> DerivationFramework::SkimmingToolEXOT14::m_n_passDiJetMass

mutableprivate

Definition at line 126 of file SkimmingToolEXOT14.h.

m_n_passGRL

std::atomic<unsigned int> DerivationFramework::SkimmingToolEXOT14::m_n_passGRL

mutableprivate

Definition at line 120 of file SkimmingToolEXOT14.h.

m_n_passJetPts

std::atomic<unsigned int> DerivationFramework::SkimmingToolEXOT14::m_n_passJetPts

mutableprivate

Definition at line 124 of file SkimmingToolEXOT14.h.

m_n_passJetsDEta

std::atomic<unsigned int> DerivationFramework::SkimmingToolEXOT14::m_n_passJetsDEta

mutableprivate

Definition at line 125 of file SkimmingToolEXOT14.h.

m_n_passJetsDPhi

std::atomic<unsigned int> DerivationFramework::SkimmingToolEXOT14::m_n_passJetsDPhi

mutableprivate

Definition at line 127 of file SkimmingToolEXOT14.h.

m_n_passLArError

std::atomic<unsigned int> DerivationFramework::SkimmingToolEXOT14::m_n_passLArError

mutableprivate

Definition at line 121 of file SkimmingToolEXOT14.h.

m_n_passPreselect

std::atomic<unsigned int> DerivationFramework::SkimmingToolEXOT14::m_n_passPreselect

mutableprivate

Definition at line 123 of file SkimmingToolEXOT14.h.

m_n_passTrigger

std::atomic<unsigned int> DerivationFramework::SkimmingToolEXOT14::m_n_passTrigger

mutableprivate

Definition at line 122 of file SkimmingToolEXOT14.h.

m_n_tot

std::atomic<unsigned int> DerivationFramework::SkimmingToolEXOT14::m_n_tot

mutableprivate

Definition at line 119 of file SkimmingToolEXOT14.h.

m_reqDiJetMass

bool DerivationFramework::SkimmingToolEXOT14::m_reqDiJetMass

private

Definition at line 81 of file SkimmingToolEXOT14.h.

m_reqGRL

bool DerivationFramework::SkimmingToolEXOT14::m_reqGRL

private

Definition at line 75 of file SkimmingToolEXOT14.h.

m_reqJetPts

bool DerivationFramework::SkimmingToolEXOT14::m_reqJetPts

private

Definition at line 79 of file SkimmingToolEXOT14.h.

m_reqJetsDEta

bool DerivationFramework::SkimmingToolEXOT14::m_reqJetsDEta

private

Definition at line 80 of file SkimmingToolEXOT14.h.

m_reqJetsDPhi

bool DerivationFramework::SkimmingToolEXOT14::m_reqJetsDPhi

private

Definition at line 82 of file SkimmingToolEXOT14.h.

m_reqLArError

bool DerivationFramework::SkimmingToolEXOT14::m_reqLArError

private

Definition at line 76 of file SkimmingToolEXOT14.h.

m_reqPreselection

bool DerivationFramework::SkimmingToolEXOT14::m_reqPreselection

private

Definition at line 78 of file SkimmingToolEXOT14.h.

m_reqTrigger

bool DerivationFramework::SkimmingToolEXOT14::m_reqTrigger

private

Definition at line 77 of file SkimmingToolEXOT14.h.

m_subleadingJetPt

double DerivationFramework::SkimmingToolEXOT14::m_subleadingJetPt

private

Definition at line 94 of file SkimmingToolEXOT14.h.

m_trigDecisionTool

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

private

Definition at line 66 of file SkimmingToolEXOT14.h.

m_triggers

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

private

Definition at line 89 of file SkimmingToolEXOT14.h.

s_MZ

const double DerivationFramework::SkimmingToolEXOT14::s_MZ

staticprivate

Definition at line 130 of file SkimmingToolEXOT14.h.

---

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

• SkimmingToolEXOT14.h
• SkimmingToolEXOT14.cxx