#include <SimpleHypoJetVectorGenerator.h>

Inheritance diagram for SimpleHypoJetVectorGenerator:

Collaboration diagram for SimpleHypoJetVectorGenerator:

Public Member Functions
	SimpleHypoJetVectorGenerator (const std::string &type, const std::string &name, const IInterface *parent)

	~SimpleHypoJetVectorGenerator ()

virtual StatusCode	initialize () override

virtual HypoJetVector	get () override

virtual std::string	toString () const override

Private Attributes
Gaudi::Property< std::vector< double > >	m_es {this, "ets", {}, "jet E values"}

Gaudi::Property< std::vector< double > >	m_etas {this, "etas", {}, "jet eta values"}

Gaudi::Property< int >	m_nbkgd {this, "n_bkgd", 0, "No of background jets to generate "}

Gaudi::Property< double >	m_bkgd_etmax {this, "bkgd_etmax", 0., "Max Et background jets"}

Gaudi::Property< double >	m_bkgd_etamax {this, "bkgd_etamax", 5., "Lower limit Eta range background "}

Gaudi::Property< bool >	m_useEtaEtNotEtaE {this, "useEtaEtNotEtaE", true, "e values or et values "}

std::default_random_engine	m_shuffle_re

Detailed Description

Definition at line 13 of file SimpleHypoJetVectorGenerator.h.

Constructor & Destructor Documentation

◆ SimpleHypoJetVectorGenerator()

SimpleHypoJetVectorGenerator::SimpleHypoJetVectorGenerator	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent
	)

Definition at line 17 of file SimpleHypoJetVectorGenerator.cxx.

                                                          :
   base_class(type, name, parent){
   
 }

◆ ~SimpleHypoJetVectorGenerator()

SimpleHypoJetVectorGenerator::~SimpleHypoJetVectorGenerator ( )

inline

Definition at line 22 of file SimpleHypoJetVectorGenerator.h.

22 {}

Member Function Documentation

◆ get()

HypoJetVector SimpleHypoJetVectorGenerator::get ( )

overridevirtual

Definition at line 51 of file SimpleHypoJetVectorGenerator.cxx.

                                   {
  
   auto factory = std::unique_ptr<ITLorentzVectorFactory>(nullptr);
   if (m_useEtaEtNotEtaE){
     factory.reset(new TLorentzVectorFactory());
   } else {
     factory.reset(new TLorentzVectorFactoryEtaE());
   }
  
   HypoJetVector result;
   
   for(std::size_t i = 0; i < m_es.size(); ++i){
     const auto& eta = m_etas[i];
     const auto& et = m_es[i];
     auto tlv = factory->make(eta, et);
     result.push_back(std::shared_ptr<const HypoJet::IJet>(new TLorentzVectorAsIJet(tlv)));
   }
   
   if(m_nbkgd){
     std::default_random_engine generator;
     std::uniform_real_distribution<double> eta_distribution(-m_bkgd_etamax,
                                 m_bkgd_etamax);
     
     std::uniform_real_distribution<double> en_distribution(5.,
                                m_bkgd_etmax);
  
     
     
     for(int i = 0; i < m_nbkgd; ++i){
       const auto& eta = eta_distribution(generator);
       const auto& en = en_distribution(generator);
       // depending on the factory, en is e or et
       auto tlv = factory->make(eta, en);
       result.push_back(std::shared_ptr<const HypoJet::IJet>(new TLorentzVectorAsIJet(tlv)));
     }
   }
  
   std::shuffle(result.begin(), result.end(), m_shuffle_re);
   return result;
 }

◆ initialize()

StatusCode SimpleHypoJetVectorGenerator::initialize ( )

overridevirtual

Definition at line 26 of file SimpleHypoJetVectorGenerator.cxx.

                                          {
   if(m_es.size() != m_etas.size()){
     ATH_MSG_ERROR("Config error: no of Et values != no of Eta values"
           << m_es.size() << " " << m_etas.size());
     return StatusCode::FAILURE;
   }
  
   return StatusCode::SUCCESS;
   
 }

◆ toString()

std::string SimpleHypoJetVectorGenerator::toString ( ) const

overridevirtual

Definition at line 37 of file SimpleHypoJetVectorGenerator.cxx.

                                                        {
   std::stringstream ss;
   ss << name() <<":: e: ";
   for(const auto& e : m_es){ss << e <<  ", ";}
   ss << " eta:  ";
   for(const auto& e : m_etas){ss << e << ", ";}
   ss<<" bkgd: "<< m_nbkgd << " max E" << m_bkgd_etmax << " max abs eta: "
     << m_bkgd_etamax << '\n';
  
   return ss.str();
 }