BesTwogam Class Reference

#include <BesTwogam.h>

Inheritance diagram for BesTwogam:

Public Member Functions
	BesTwogam (const string &name, ISvcLocator *pSvcLocator)
StatusCode	initialize ()
StatusCode	execute ()
StatusCode	finalize ()

Detailed Description

Definition at line 24 of file BesTwogam.h.

Constructor & Destructor Documentation

BesTwogam()

BesTwogam::BesTwogam	(	const string &	name,
		ISvcLocator *	pSvcLocator )

Definition at line 146 of file BesTwogam.cxx.

    : Algorithm( name, pSvcLocator ) {
  declareProperty( "InitialSeed", m_seed = 1001 );
  declareProperty( "CMEnergy", m_cmEnergy = 3.097 ); // 2*Ebeam [GeV]
  declareProperty( "MinimumW2", m_m2min = 0.02 );
  declareProperty( "Unweighted", m_unw = 1 );
  declareProperty( "CardFileName", m_fcard = "BesTwogam.conf" );
  // interface set pythia pars, pingrg
  m_pypars.clear();
  declareProperty( "setPythiaPars", m_pypars );
}

Referenced by BesTwogam().

Member Function Documentation

execute()

StatusCode BesTwogam::execute

(

)

Definition at line 190 of file BesTwogam.cxx.

                              {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "BesTwogam executing" << endmsg;
  HepMC::HEPEVT_Wrapper::set_max_number_entries( 2000 );
  HepMC::HEPEVT_Wrapper::set_sizeof_real( 8 );
  HepMC::IO_HEPEVT HepEvtIO;
  HEPEVT_CLEAN();
  GEN1EVT();
  // PRINTEVT();
  if ( CFLAGS.GOODEV != 1 )
  {
    log << MSG::ERROR << " BesTwogam: fail to generate good event" << endl;
    return StatusCode::FAILURE;
  }
  m_numberEvent++;
  if ( log.level() < MSG::INFO ) LULIST( 1 );
  LUHEPC( 1 );
  //  HEPEVT_PRINT();
  HepMC::GenEvent* evt = HepEvtIO.read_next_event();
  evt->set_event_number( m_numberEvent );
  evt->set_signal_process_id( 1 );
  // Check if the McCollection already exists
  SmartDataPtr<McGenEventCol> anMcCol( eventSvc(), "/Event/Gen" );
  if ( anMcCol != 0 )
  {
    // Add event to existing collection
    MsgStream log( msgSvc(), name() );
    log << MSG::INFO << "Add McGenEvent to existing collection" << endmsg;
    McGenEvent* mcEvent = new McGenEvent( evt );
    anMcCol->push_back( mcEvent );
  }
  else
  {
    // Create Collection and add  to the transient store
    McGenEventCol* mcColl  = new McGenEventCol;
    McGenEvent*    mcEvent = new McGenEvent( evt );
    mcColl->push_back( mcEvent );
    StatusCode sc = eventSvc()->registerObject( "/Event/Gen", mcColl );
    if ( sc != StatusCode::SUCCESS )
    {
      log << MSG::ERROR << "Could not register McGenEvent" << endmsg;
      delete mcColl;
      delete evt;
      delete mcEvent;
      return StatusCode::FAILURE;
    }
  }
  return StatusCode::SUCCESS;
}

finalize()

StatusCode BesTwogam::finalize

(

)

Definition at line 240 of file BesTwogam.cxx.

                               {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "BesTwogam finalized" << endmsg;
  FINISH();
  return StatusCode::SUCCESS;
}

initialize()

StatusCode BesTwogam::initialize

(

)

Definition at line 158 of file BesTwogam.cxx.

                                 {
  MsgStream log( msgSvc(), name() );
  log << MSG::WARNING << "BesTwogam initialize" << endmsg;
 
  if ( m_fcard.size() > 254 )
  {
    log << MSG::ERROR << "Too big name of BesTwogam card file" << endmsg;
    return StatusCode::FAILURE;
  }
  //  FNAMES.dummy=61777;
  strcpy( FNAMES.cardName, m_fcard.c_str() );
  FNAMES.i_cardName = m_fcard.size();
 
  RDCUTS();
  ISEEDC.ISEED  = m_seed;
  EXCUTS.EBEAME = 0.5 * m_cmEnergy;
  EXCUTS.W2MINE = m_m2min;
  INPARC.NTUNW  = m_unw;
 
  SETDEF();
  QCDSTA();
  DOCUTS();
  ESTMAXCS();
  PRTCUT();
  //  USRINI();
  m_numberEvent = 0;
  // for pythia parameter setting,pingrg-2013-12-11
  for ( int i = 0; i < m_pypars.size(); i++ )
  { pygive_( m_pypars[i].c_str(), strlen( m_pypars[i].c_str() ) ); }
  return StatusCode::SUCCESS;
}

---

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

• BesTwogam.h
• BesTwogam.cxx