DijetMassCondition.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/
 
#include "./DijetMassCondition.h"
#include "./ITrigJetHypoInfoCollector.h"
#include <sstream>
#include <stdexcept>
#include <TLorentzVector.h>
#include <limits>
#include <memory>
 
DijetMassCondition::DijetMassCondition(double massMin,
                       double massMax):
  m_min(massMin),
  m_max(massMax){
  }
  
 
bool
DijetMassCondition::isSatisfied(const HypoJetVector& ips,
                  const std::unique_ptr<ITrigJetHypoInfoCollector>& collector) const{
  
  if(ips.size() != 2){
    std::stringstream ss;
    ss << "DijetMassCondition::isSatisfied must see exactly 2 particles, but received "
       << ips.size()
       << '\n';
    
    throw std::runtime_error(ss.str());
  }
  
  auto j0 = ips[0];
  auto j1 = ips[1];
 
  auto rj0 = j0 -> p4();
  auto rj1 = j1 -> p4();
 
  auto mass = (rj0 + rj1).M();
 
  bool pass = mass >= m_min  and  mass < m_max;
 
   if(collector){
     std::stringstream ss0;
     const void* address = static_cast<const void*>(this);
     ss0 << "DijetMassCondition: (" << address << ") "
     <<  mass << " "
     << std::boolalpha << pass <<  " jet group: \n";
 
     std::stringstream ss1;
 
     for(const auto& ip : ips){
       address = static_cast<const void*>(ip.get());
       ss1 << "    "  << address << " " << ip->eta() << " e " << ip->e() << '\n';
     }
     ss1 << '\n';
     collector -> collect(ss0.str(), ss1.str());
   }
   return pass;
 
}
 
std::string DijetMassCondition::toString() const {
 
 
  std::stringstream ss;
  const void* address = static_cast<const void*>(this);
  ss << "DijetMassCondition: (" << address << ") Capacity: " << s_capacity
    
     << " mass min: " << m_min
     << " mass max: " << m_max <<'\n';
  
  return ss.str();
}