DVMissingETFilterTool.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
*/
 
// DVMissingETFilterTool.cxx, (c) ATLAS Detector software
 
#include "LongLivedParticleDPDMaker/DVMissingETFilterTool.h"
#include <vector>
#include <string>
 
// Constructor
DerivationFramework::DVMissingETFilterTool::DVMissingETFilterTool( const std::string& t,
                                                 const std::string& n,
                                                 const IInterface* p ) : 
  AthAlgTool(t,n,p),
  m_ntot(0),
  m_npass(0),
  m_metCut(50000.0),
  m_jetPtCut(50000.0),
  m_applyDeltaPhiCut(false),
  m_deltaPhiCut(2.)
  {
    declareInterface<DerivationFramework::ISkimmingTool>(this);
    declareProperty("METCut", m_metCut);
    declareProperty("JetPtCut", m_jetPtCut);
    declareProperty("applyDeltaPhiCut", m_applyDeltaPhiCut);
    declareProperty("DeltaPhiCut", m_deltaPhiCut);
  }
  
// Destructor
DerivationFramework::DVMissingETFilterTool::~DVMissingETFilterTool() {
}  
 
// Athena initialize and finalize
StatusCode DerivationFramework::DVMissingETFilterTool::initialize()
{
     ATH_MSG_VERBOSE("initialize() ...");
     ATH_CHECK(m_metSGKey.initialize());
     ATH_CHECK(m_jetSGKey.initialize());
     return StatusCode::SUCCESS;
}
StatusCode DerivationFramework::DVMissingETFilterTool::finalize()
{
     ATH_MSG_VERBOSE("finalize() ...");
     ATH_MSG_INFO("Processed "<< m_ntot <<" events, "<< m_npass<<" events passed filter ");
     return StatusCode::SUCCESS;
}
 
// The filter itself
bool DerivationFramework::DVMissingETFilterTool::eventPassesFilter() const
{
     ++m_ntot;
 
     SG::ReadHandle<xAOD::MissingETContainer> metContainer(m_metSGKey);
     if( !metContainer.isValid() ) {
       msg(MSG::WARNING) << "No MET container found, will skip this event" << endmsg;
       return false;
     } 
 
     if (metContainer->size() ==1) {
     
       double MET = metContainer->at(0)->met();
       bool passesEvent= (MET > m_metCut) ;
 
 
       
       if (m_applyDeltaPhiCut) {
     double phiMET=metContainer->at(0)->phi();
         SG::ReadHandle<xAOD::JetContainer> jetContainer(m_jetSGKey);
         if( !jetContainer.isValid() ) {
       msg(MSG::WARNING) << "No jet container found, will skip this event" << endmsg;
       return false;
     }
     if ((jetContainer->size() > 0) && (jetContainer->at(0)->pt() > m_jetPtCut)){
       double phiJet = jetContainer->at(0)->phi();
       double deltaPhi = std::abs(phiMET-phiJet);
       if (deltaPhi > M_PI) deltaPhi = 2.0*M_PI - deltaPhi;
       passesEvent = passesEvent && (deltaPhi > m_deltaPhiCut);
     }
       }
       
       if (passesEvent) {
     ++m_npass;
       }
       return passesEvent;
     }
     return false;
 
}