BesTwogam Class Reference

#include <BesTwogam.h>

List of all members.

Public Member Functions
	BesTwogam (const string &name, ISvcLocator *pSvcLocator)
StatusCode	initialize ()
StatusCode	execute ()
StatusCode	finalize ()
Private Attributes
int	m_numberEvent
int	m_seed
int	m_unw
float	m_cmEnergy
float	m_m2min
std::string	m_fcard
std::vector< std::string >	m_pypars

---

Detailed Description

Definition at line 24 of file BesTwogam.h.

---

Constructor & Destructor Documentation

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

Definition at line 152 of file BesTwogam.cxx.

References m_cmEnergy, m_fcard, m_m2min, m_pypars, m_seed, and m_unw.

                                                                :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);
}

---

Member Function Documentation

StatusCode BesTwogam::execute

(

)

Definition at line 196 of file BesTwogam.cxx.

References CFLAGS, calibUtil::ERROR, GEN1EVT, HEPEVT_CLEAN, Bes_Common::INFO, LUHEPC, LULIST, and m_numberEvent.

{ 
  MsgStream log(messageService(), name());
  log << MSG::INFO << "BesTwogam executing" << endreq;
  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(messageService(), name());
    log << MSG::INFO << "Add McGenEvent to existing collection" << endreq;
    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" << endreq;
      delete mcColl;
      delete evt;
      delete mcEvent;
      return StatusCode::FAILURE;
    }
  }
  return StatusCode::SUCCESS; 
}

StatusCode BesTwogam::finalize

(

)

Definition at line 242 of file BesTwogam.cxx.

References FINISH, and Bes_Common::INFO.

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

StatusCode BesTwogam::initialize

(

)

Definition at line 164 of file BesTwogam.cxx.

References EvtCyclic3::c_str(), DOCUTS, calibUtil::ERROR, ESTMAXCS, EXCUTS, FNAMES, genRecEmupikp::i, INPARC, ISEEDC, m_cmEnergy, m_fcard, m_m2min, m_numberEvent, m_pypars, m_seed, m_unw, PRTCUT, pygive_(), QCDSTA, RDCUTS, SETDEF, and Bes_Common::WARNING.

                                {
  MsgStream log(messageService(), name());
  log << MSG::WARNING << "BesTwogam initialize" << endreq;

  if(m_fcard.size()>254){
  log << MSG::ERROR << "Too big name of BesTwogam card file" << endreq;
  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;
}

---

Member Data Documentation

float BesTwogam::m_cmEnergy [private]

Definition at line 36 of file BesTwogam.h.

Referenced by BesTwogam(), and initialize().

std::string BesTwogam::m_fcard [private]

Definition at line 38 of file BesTwogam.h.

Referenced by BesTwogam(), and initialize().

float BesTwogam::m_m2min [private]

Definition at line 37 of file BesTwogam.h.

Referenced by BesTwogam(), and initialize().

int BesTwogam::m_numberEvent [private]

Definition at line 33 of file BesTwogam.h.

Referenced by execute(), and initialize().

std::vector<std::string> BesTwogam::m_pypars [private]

Definition at line 39 of file BesTwogam.h.

Referenced by BesTwogam(), and initialize().

int BesTwogam::m_seed [private]

Definition at line 34 of file BesTwogam.h.

Referenced by BesTwogam(), and initialize().

int BesTwogam::m_unw [private]

Definition at line 35 of file BesTwogam.h.

Referenced by BesTwogam(), and initialize().

---