#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();
}