DerivationFramework::SkimmingToolEXOT5 Class Reference

#include <SkimmingToolEXOT5.h>

Inheritance diagram for DerivationFramework::SkimmingToolEXOT5:

Collaboration diagram for DerivationFramework::SkimmingToolEXOT5:

Public Member Functions
	SkimmingToolEXOT5 (const std::string &t, const std::string &n, const IInterface *p)
	~SkimmingToolEXOT5 ()
virtual StatusCode	initialize () override
virtual StatusCode	finalize () override
virtual bool	eventPassesFilter () const override

Private Attributes
ToolHandle< IJetCalibrationTool >	m_jetCalibrationTool
std::string	m_jetSGKey
double	m_uncalibMonoJetPt
double	m_monoJetPt
double	m_leadingJetPt
double	m_subleadingJetPt
double	m_Mjj
std::atomic< unsigned int >	m_ntot
std::atomic< unsigned int >	m_npass

Detailed Description

Definition at line 22 of file SkimmingToolEXOT5.h.

Constructor & Destructor Documentation

SkimmingToolEXOT5()

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

Definition at line 27 of file SkimmingToolEXOT5.cxx.

                                                                                 :
  base_class(t, n, p),
  m_ntot(0),
  m_npass(0)
{
 
 
  declareProperty("JetContainer",          m_jetSGKey = "AntiKt4EMTopoJets");
 
  declareProperty("UncalibMonoJetPtCut",   m_uncalibMonoJetPt = 100000.);
  declareProperty("MonoJetPtCut",          m_monoJetPt = 100000.);
  declareProperty("LeadingJetPtCut",       m_leadingJetPt = 40000.);
  declareProperty("SubleadingJetPtCut",    m_subleadingJetPt = 40000.);
  declareProperty("DiJetMassCut",          m_Mjj = 150000.);
}

~SkimmingToolEXOT5()

DerivationFramework::SkimmingToolEXOT5::~SkimmingToolEXOT5

(

)

Definition at line 45 of file SkimmingToolEXOT5.cxx.

                                                         {
}

Member Function Documentation

eventPassesFilter()

bool DerivationFramework::SkimmingToolEXOT5::eventPassesFilter ( ) const

overridevirtual

Definition at line 68 of file SkimmingToolEXOT5.cxx.

{
 
  m_ntot++;
 
  const xAOD::EventInfo* eventInfo = nullptr;
  ATH_CHECK(evtStore()->retrieve(eventInfo), false);
 
  const xAOD::JetContainer* jets = nullptr;
  ATH_CHECK(evtStore()->retrieve(jets, m_jetSGKey), false);
 
  bool passUncalibMonojetCut = false;
  bool passRecoJetCuts = false;
  bool passTruthJetCuts = false;
 
  std::unique_ptr<xAOD::JetContainer> recoJets = make_unique<xAOD::JetContainer>();
  std::unique_ptr<xAOD::JetAuxContainer> recoJetsAux = make_unique<xAOD::JetAuxContainer>();
  recoJets->setStore(recoJetsAux.get());
 
  std::pair<xAOD::JetContainer*, xAOD::ShallowAuxContainer*> jets_shallowCopy = xAOD::shallowCopyContainer(*jets);
  for (const auto jet : *jets_shallowCopy.first) {
    if (jet->pt() > 100000.) passUncalibMonojetCut = true;
  }
  if (m_jetCalibrationTool->applyCalibration(*jets_shallowCopy.first).isFailure()) {
    ATH_MSG_ERROR("Problem applying jet calibration");
    return false;
  }
  for (const auto jet : *jets_shallowCopy.first) {
    xAOD::Jet* newJet = new xAOD::Jet();
    newJet->makePrivateStore(*jet);
    recoJets->push_back(newJet);
  }
  delete jets_shallowCopy.first;
  delete jets_shallowCopy.second;
 
  if (recoJets->size() > 1) std::partial_sort(recoJets->begin(), recoJets->begin()+2, recoJets->end(), DescendingPt());
  float mjj = 0;
  if (recoJets->size() > 1) {
    TLorentzVector jet1 = recoJets->at(0)->p4();
    TLorentzVector jet2 = recoJets->at(1)->p4();
    auto dijet = jet1 + jet2;
    mjj = dijet.M();
  }
  if ((recoJets->size() > 0 && recoJets->at(0)->pt() > m_monoJetPt) || (recoJets->size() > 1 && recoJets->at(0)->pt() > m_leadingJetPt && recoJets->at(1)->pt() > m_subleadingJetPt && mjj > m_Mjj)) passRecoJetCuts = true;
 
  if (eventInfo->eventType(xAOD::EventInfo::IS_SIMULATION)) {
    const xAOD::JetContainer* truthJetContainer = nullptr;
    ATH_CHECK(evtStore()->retrieve(truthJetContainer, "AntiKt4TruthJets"), false);
 
    std::unique_ptr<xAOD::JetContainer> truthJets = make_unique<xAOD::JetContainer>();
    std::unique_ptr<xAOD::JetAuxContainer> truthJetsAux = make_unique<xAOD::JetAuxContainer>();
    truthJets->setStore(truthJetsAux.get());
 
    for (const auto truthJet : *truthJetContainer) {
      xAOD::Jet* newTruthJet = new xAOD::Jet();
      newTruthJet->makePrivateStore(*truthJet);
      truthJets->push_back(newTruthJet);
    }
 
    if (truthJets->size() > 1) std::partial_sort(truthJets->begin(), truthJets->begin()+2, truthJets->end(), DescendingPt());
    float truthMjj = 0;
    if (truthJets->size() > 1) {
      TLorentzVector truthJet1 = truthJets->at(0)->p4();
      TLorentzVector truthJet2 = truthJets->at(1)->p4();
      auto truthDijet = truthJet1 + truthJet2;
      truthMjj = truthDijet.M();
    }
    if ((truthJets->size() > 0 && truthJets->at(0)->pt() > m_monoJetPt) || (truthJets->size() > 1 && truthJets->at(0)->pt() > m_leadingJetPt && truthJets->at(1)->pt() > m_subleadingJetPt && truthMjj > m_Mjj)) passTruthJetCuts = true;
  }
 
  bool acceptEvent = passUncalibMonojetCut || passRecoJetCuts || passTruthJetCuts;
 
  m_npass++;
  return acceptEvent;
 
}

finalize()

StatusCode DerivationFramework::SkimmingToolEXOT5::finalize ( )

overridevirtual

Definition at line 59 of file SkimmingToolEXOT5.cxx.

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

initialize()

StatusCode DerivationFramework::SkimmingToolEXOT5::initialize ( )

overridevirtual

Definition at line 48 of file SkimmingToolEXOT5.cxx.

{
  ATH_MSG_VERBOSE("initialize() ..");
 
  m_jetCalibrationTool.setTypeAndName("JetCalibrationTool/EXOT5JESTool");
  CHECK(m_jetCalibrationTool.retrieve());
  ATH_MSG_INFO("Retrieved tool: " << m_jetCalibrationTool);
 
  return StatusCode::SUCCESS;
}

Member Data Documentation

m_jetCalibrationTool

ToolHandle<IJetCalibrationTool> DerivationFramework::SkimmingToolEXOT5::m_jetCalibrationTool

private

Definition at line 36 of file SkimmingToolEXOT5.h.

m_jetSGKey

std::string DerivationFramework::SkimmingToolEXOT5::m_jetSGKey

private

Definition at line 37 of file SkimmingToolEXOT5.h.

m_leadingJetPt

double DerivationFramework::SkimmingToolEXOT5::m_leadingJetPt

private

Definition at line 41 of file SkimmingToolEXOT5.h.

m_Mjj

double DerivationFramework::SkimmingToolEXOT5::m_Mjj

private

Definition at line 43 of file SkimmingToolEXOT5.h.

m_monoJetPt

double DerivationFramework::SkimmingToolEXOT5::m_monoJetPt

private

Definition at line 40 of file SkimmingToolEXOT5.h.

m_npass

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

mutableprivate

Definition at line 46 of file SkimmingToolEXOT5.h.

m_ntot

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

mutableprivate

Definition at line 45 of file SkimmingToolEXOT5.h.

m_subleadingJetPt

double DerivationFramework::SkimmingToolEXOT5::m_subleadingJetPt

private

Definition at line 42 of file SkimmingToolEXOT5.h.

m_uncalibMonoJetPt

double DerivationFramework::SkimmingToolEXOT5::m_uncalibMonoJetPt

private

Definition at line 39 of file SkimmingToolEXOT5.h.

---

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

• SkimmingToolEXOT5.h
• SkimmingToolEXOT5.cxx