DummyParticleGun Class Reference

Inheritance diagram for DummyParticleGun:

Public Member Functions
	DummyParticleGun (const std::string &name, ISvcLocator *pSvcLocator)
StatusCode	initialize () override
McGenEventCol	operator() (const EventContext &ctx) const override

Detailed Description

Definition at line 20 of file DummyParticleGun.cxx.

Constructor & Destructor Documentation

DummyParticleGun()

DummyParticleGun::DummyParticleGun ( const std::string & name, ISvcLocator * pSvcLocator )

inline

Definition at line 22 of file DummyParticleGun.cxx.

      : Transformer( name, pSvcLocator, { "McGenEventLocation", "/Event/Gen" } ) {}

Member Function Documentation

initialize()

StatusCode DummyParticleGun::initialize ( )

inlineoverride

Definition at line 25 of file DummyParticleGun.cxx.

                                   {
    info() << "Initializing DummyParticleGun" << endmsg;
 
    StatusCode sc = Transformer::initialize();
    if ( !sc.isSuccess() ) return sc;
 
    // check args
    if ( m_parts.size() == 0 )
    {
      error() << "No PDG codes specified" << endmsg;
      return StatusCode::FAILURE;
    }
 
    if ( m_p_range.size() != 2 )
    {
      error() << "Momentum range must be of size 2" << endmsg;
      return StatusCode::FAILURE;
    }
 
    if ( m_costheta_range.size() != 2 )
    {
      error() << "Cos(theta) range must be of size 2" << endmsg;
      return StatusCode::FAILURE;
    }
 
    if ( m_phi_range.size() != 2 )
    {
      error() << "Phi range must be of size 2" << endmsg;
      return StatusCode::FAILURE;
    }
 
    // set ranges
    m_p1        = m_p_range[0];
    m_p2        = m_p_range[1];
    m_costheta1 = m_costheta_range[0];
    m_costheta2 = m_costheta_range[1];
    m_phi1      = m_phi_range[0];
    m_phi2      = m_phi_range[1];
 
    if ( m_p2 <= m_p1 ) m_p2 = m_p1;
    if ( m_costheta2 <= m_costheta1 ) m_costheta2 = m_costheta1;
    if ( m_phi2 <= m_phi1 ) m_phi2 = m_phi1;
 
    // retrieve random enging
    m_rndm_svc = service<IBesRndmGenSvc>( "BesRndmGenSvc" );
    if ( !m_rndm_svc )
    {
      error() << "Could not get random number service" << endmsg;
      return StatusCode::FAILURE;
    }
 
    m_rndm_engine = m_rndm_svc->GetEngine( "DummyParticleGun" );
    if ( !m_rndm_engine )
    {
      error() << "Could not get random engine" << endmsg;
      return StatusCode::FAILURE;
    }
 
    // retrieve particle data table
    auto part_prop_svc = service<IPartPropSvc>( "PartPropSvc" );
    if ( !part_prop_svc )
    {
      error() << "Could not get Particle Properties Service" << endmsg;
      return StatusCode::FAILURE;
    }
 
    m_particle_table = part_prop_svc->PDT();
    if ( !m_particle_table )
    {
      error() << "Could not get Particle Data Table" << endmsg;
      return StatusCode::FAILURE;
    }
 
    return StatusCode::SUCCESS;
  }

operator()()

McGenEventCol DummyParticleGun::operator() ( const EventContext & ctx ) const

inlineoverride

Definition at line 101 of file DummyParticleGun.cxx.

                                                                     {
    GenEvent* evt = new GenEvent( 1, 1 );
 
    // make vertex (0, 0, 0, 0)
    GenVertex* vtx = new GenVertex();
    evt->add_vertex( vtx );
 
    // sample particles
    for ( auto pdg_code : m_parts.value() )
    {
      const HepPDT::ParticleData* part =
          m_particle_table->particle( HepPDT::ParticleID( abs( pdg_code ) ) );
 
      double mass = part ? part->mass().value() : 0.0;
 
      // sample momentum
      double p         = CLHEP::RandFlat::shoot( m_rndm_engine, m_p1, m_p2 );
      double cos_theta = CLHEP::RandFlat::shoot( m_rndm_engine, m_costheta1, m_costheta2 );
      double phi       = CLHEP::RandFlat::shoot( m_rndm_engine, m_phi1, m_phi2 );
 
      // calculate px, py, pz
      double sin_theta = sqrt( 1.0 - cos_theta * cos_theta );
      double pt        = p * sin_theta;
 
      double px = pt * cos( phi );
      double py = pt * sin( phi );
      double pz = p * cos_theta;
 
      CLHEP::HepLorentzVector p4;
      p4.setVectM( CLHEP::Hep3Vector( px, py, pz ), mass );
 
      vtx->add_particle_out( new GenParticle( p4, pdg_code, 1 ) );
    }
 
    std::vector<long> seeds{ m_rndm_engine->getSeeds()[0], m_rndm_engine->getSeeds()[1] };
    evt->set_random_states( seeds );
 
    // add event to collection
    McGenEventCol mc_col;
    mc_col.push_back( new McGenEvent( evt ) );
    return mc_col;
  };

---

The documentation for this class was generated from the following file:

• DummyParticleGun.cxx