BesBdkRc Class Reference

#include <BesBdkRc.h>

Inheritance diagram for BesBdkRc:

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

Detailed Description

Definition at line 21 of file BesBdkRc.h.

Constructor & Destructor Documentation

BesBdkRc()

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

Definition at line 126 of file BesBdkRc.cxx.

    : Algorithm( name, pSvcLocator ) {
  //  declareProperty("InitialSeed",m_iseed=1001);
  declareProperty( "CMEnergy", m_CMEnergy = 3.097 ); // 2*Ebeam [GeV]
  declareProperty( "W2min", m_w2min = 0.02 );        // Cut on invariant gamma-gamma mass
  declareProperty( "EstimatedMaxWeight", m_ewe = 2.5 );
  declareProperty( "SoftPhotonMaxEnergy", m_kzero = 0.002 );
  qmass[0] = 0.2;
  qmass[1] = 0.2;
  qmass[2] = 0.5;
  qmass[3] = 1.5;
  qmass[4] = 4.5;
  qmass[5] = 180.;
  declareProperty( "MaxNTry", m_maxNTry );
  declareProperty( "FinalState", m_ifinal );
  declareProperty( "MinTheta", m_tcmin );
  declareProperty( "MinMomentum", m_pcmin );
 
  m_numberEvent = 0;
  toRad         = M_PI / 180.0;
  toDeg         = 180.0 / M_PI;
}

Member Function Documentation

execute()

StatusCode BesBdkRc::execute

(

)

Definition at line 189 of file BesBdkRc.cxx.

                             {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "BesBdkRc executing" << endmsg;
  HepMC::HEPEVT_Wrapper::set_max_number_entries( 2000 );
  HepMC::HEPEVT_Wrapper::set_sizeof_real( 8 );
  HepMC::IO_HEPEVT HepEvtIO;
 
  HEPEVT_CLEAN();
  GEN1EVT();
 
  if ( FLAGS.GOODEVT != 1 )
  {
    log << MSG::ERROR << " BesBdkRc: fail to generate good event" << endl;
    return StatusCode::FAILURE;
  }
 
  m_numberEvent++;
  if ( log.level() < MSG::INFO ) LULIST( 1 );
  //  LULIST(1);
  LUHEPC( 1 );
  HepMC::GenEvent* evt = HepEvtIO.read_next_event();
  evt->set_event_number( m_numberEvent );
  evt->set_signal_process_id( 1 );
  //  evt->print();
  // 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;
    }
  }
  // string s;
  // getline(cin,s);
 
  return StatusCode::SUCCESS;
}

finalize()

StatusCode BesBdkRc::finalize

(

)

Definition at line 245 of file BesBdkRc.cxx.

                              {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "BesBdkRc finalized" << endmsg;
  int itot;
  itot = XSECTC.ITOT;
  double cstot;
  double cserr;
  FINISH( itot, cstot, cserr );
  float effcut   = 0;
  float cscut    = 0;
  float efferr   = 0;
  float cscuterr = 0;
  if ( XSECTC.NEVUNW )
  {
    effcut   = float( DELPHC.NEVCUT ) / float( XSECTC.NEVUNW );
    cscut    = effcut * cstot;
    efferr   = sqrt( effcut * ( 1.0 - effcut ) / float( XSECTC.NEVUNW ) );
    cscuterr = sqrt( cstot * efferr * cstot * efferr + effcut * effcut * cserr * cserr );
  }
  printf( "BDKRC SUMMARY: Cross section after user cuts= %G +- %G nb\n", cscut, cscuterr );
  printf( "               Cut acceptance              = %G +- %G \n", effcut, efferr );
  return StatusCode::SUCCESS;
}

initialize()

StatusCode BesBdkRc::initialize

(

)

Definition at line 149 of file BesBdkRc.cxx.

                                {
  MsgStream log( msgSvc(), name() );
  log << MSG::WARNING << "BesBdkRc initialize" << endmsg;
 
  static const bool CREATEIFNOTTHERE( true );
  StatusCode        RndmStatus = service( "BesRndmGenSvc", p_BesRndmGenSvc, CREATEIFNOTTHERE );
  if ( !RndmStatus.isSuccess() || 0 == p_BesRndmGenSvc )
  {
    log << MSG::ERROR << " Could not initialize Random Number Service" << endmsg;
    return RndmStatus;
  }
  CLHEP::HepRandomEngine* engine = p_BesRndmGenSvc->GetEngine( "BesBdkRc" );
  engine->showStatus();
  BesBdkRcRandom::setRandomEngine( engine );
 
  //"INIT" subroutine action is replaced by c++ code
  // for gen_init  fortran code is called
 
  FINCUT.TCHMIN = m_tcmin * toRad;
  FINCUT.PCHMIN = m_pcmin;
  //  RLUXGO(3,m_iseed,0,0);
 
  DELPHC.MAXNTRY = m_maxNTry;
 
  COLCHC.IFINAL = m_ifinal;
 
  INPUT.EB = 0.5 * m_CMEnergy;
 
  LOCALC.W2MINR = m_w2min;
  LOCALC.EBEAM  = 0.5 * m_CMEnergy;
  LOCALC.MAXESW = m_ewe;
  LOCALC.KZERO  = m_kzero;
  //  LOCALC.ISEED=m_iseed;
  for ( int i = 0; i < 6; i++ ) LOCALC.QMASS[i] = qmass[i];
  GEN_INIT();
  XSECTC.NEVUNW = 0;
 
  return StatusCode::SUCCESS;
}

---

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

• BesBdkRc.h
• BesBdkRc.cxx