#include <RandomSignalHypoJetVectorGenerator.h>

Inheritance diagram for RandomSignalHypoJetVectorGenerator:

Collaboration diagram for RandomSignalHypoJetVectorGenerator:

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

	~RandomSignalHypoJetVectorGenerator ()

virtual StatusCode	initialize () override

virtual HypoJetVector	get () override

virtual std::string	toString () const override

Private Attributes
Gaudi::Property< int >	m_n_sig_lo {this, "n_sig_lo", 2, "Lower bound nsig jets to generate "}

Gaudi::Property< int >	m_n_sig_hi {this, "n_sig_hi", 6, "Upper bound nsig jets to generate "}

Gaudi::Property< double >	m_e_sig_lo {this, "e_sig_lo", 50., "Lower bound e sig to generate "}

Gaudi::Property< double >	m_e_sig_hi {this, "e_sig_hi", 100., "Upper bound e sig to generate "}

Gaudi::Property< double >	m_eta_sig_hi {this, "eta_sig_hi", 3.2, "Upper bound eta sig to generate "}

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

RandUint	m_randNSig

RandDouble	m_randEtaSig

RandDouble	m_randEnSig

RandDouble	m_randEtaBkd

RandDouble	m_randEnBkd

RandUint	m_flip

Detailed Description

Definition at line 15 of file RandomSignalHypoJetVectorGenerator.h.

Constructor & Destructor Documentation

◆ RandomSignalHypoJetVectorGenerator()

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

Definition at line 16 of file RandomSignalEventGenerator.cxx.

                                                                  :
   base_class(type, name, parent){
   
 }

◆ ~RandomSignalHypoJetVectorGenerator()

RandomSignalHypoJetVectorGenerator::~RandomSignalHypoJetVectorGenerator ( )

inline

Definition at line 24 of file RandomSignalHypoJetVectorGenerator.h.

24 {}

Member Function Documentation

◆ get()

HypoJetVector RandomSignalHypoJetVectorGenerator::get ( )

overridevirtual

Definition at line 63 of file RandomSignalEventGenerator.cxx.

                                         {
  
   auto factory = std::unique_ptr<ITLorentzVectorFactory>(nullptr);
   if (m_useEtaEtNotEtaE){
     factory.reset(new TLorentzVectorFactory());
   } else {
     factory.reset(new TLorentzVectorFactoryEtaE());
   }
  
   HypoJetVector result;
  
   auto nj =  m_randNSig();
   for(std::size_t ij = 0; ij  < nj;  ++ij){
     const auto& eta = m_randEtaSig();
     const auto& en = m_randEnSig();
     // depending on the factory, en is e or et
     auto tlv = factory->make(eta, en);
     if(m_flip()){
       tlv.SetZ(-tlv.Z()); // randomly flip the sign of the +vly generated tlv.
     }
     result.push_back(std::shared_ptr<const HypoJet::IJet>(new TLorentzVectorAsIJet(tlv)));
   }
  
   if(m_nbkgd){
     for(int i = 0; i < m_nbkgd; ++i){
       const auto& eta = m_randEtaBkd();
       const auto& en = m_randEnBkd();
       // depending on the factory, en is e or et
       auto tlv = factory->make(eta, en);
       if(m_flip()){
     tlv.SetZ(-tlv.Z()); // randomly flip the sign of the +vly generated tlv.
       }  
       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 RandomSignalHypoJetVectorGenerator::initialize ( )

overridevirtual

Definition at line 25 of file RandomSignalEventGenerator.cxx.

                                                {
   if(m_n_sig_lo >= m_n_sig_hi ||
      m_e_sig_lo >= m_e_sig_hi
      ){
     ATH_MSG_ERROR("Config error: lower signal limit exceeds uper limit");
     return StatusCode::FAILURE;
   }
  
   m_randNSig = RandUint(m_n_sig_lo, m_n_sig_hi);
   m_flip = RandUint(0, 1);
  
   m_randEtaSig = RandDouble(-m_eta_sig_hi, m_eta_sig_hi);
   m_randEnSig = RandDouble(m_e_sig_lo, m_e_sig_hi);
  
   m_randEtaBkd =  RandDouble(-m_bkgd_etamax, m_bkgd_etamax);
   m_randEnBkd = RandDouble(5., m_bkgd_etmax);
  
   return StatusCode::SUCCESS;  
 }

◆ toString()

std::string RandomSignalHypoJetVectorGenerator::toString ( ) const

overridevirtual

Definition at line 45 of file RandomSignalEventGenerator.cxx.

                                                              {
   std::stringstream ss;
   ss << name() <<":: "
      << m_n_sig_lo << " "
      << m_n_sig_hi << " "
      << m_e_sig_lo << " "
      << m_e_sig_hi << " "
      << m_eta_sig_hi << " ";
     
   ss << m_nbkgd << " "
      << m_bkgd_etmax << " "
      << m_bkgd_etamax << '\n';
  
   return ss.str();
 }