QjetMassCondition.cxx

Go to the documentation of this file.

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