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 129 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 192 of file BesBdkRc.cxx.

{ 
  MsgStream log(messageService(), name());
  log << MSG::INFO << "BesBdkRc executing" << endreq;
  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(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;
    }
  }
  // string s;
  // getline(cin,s);
 
  return StatusCode::SUCCESS; 
}

finalize()

StatusCode BesBdkRc::finalize

(

)

Definition at line 244 of file BesBdkRc.cxx.

{
  MsgStream log(messageService(), name());
  log << MSG::INFO << "BesBdkRc finalized" << endreq;
  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 152 of file BesBdkRc.cxx.

                               {
  MsgStream log(messageService(), name());
  log << MSG::WARNING << "BesBdkRc initialize" << endreq;
 
  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" << endreq;                      
    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