Bhlumi Class Reference

#include <Bhlumi.h>

List of all members.

Public Member Functions
	Bhlumi (const std::string &name, ISvcLocator *pSvcLocator)
StatusCode	initialize ()
StatusCode	execute ()
StatusCode	finalize ()
Private Attributes
double	xpar [100]
int	npar [100]
IBesRndmGenSvc *	p_BesRndmGenSvc
int	m_angleMode
double	m_cmEnergy
double	m_minThetaAngle
double	m_maxThetaAngle
double	m_infraredCut
std::vector< int >	m_initSeed

---

Detailed Description

Definition at line 25 of file Bhlumi.h.

---

Constructor & Destructor Documentation

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

---

Member Function Documentation

StatusCode Bhlumi::execute

(

)

Definition at line 221 of file Bhlumi.cxx.

References BHLUMI, DUMPS, calibUtil::ERROR, Bes_Common::INFO, MOMSET, npar, and xpar.

{ 
  MsgStream log(messageService(), name());
  log << MSG::INFO << "Bhlumi executing" << endreq;


  BHLUMI( 0,xpar,npar);

  if( log.level() < MSG::INFO )
  {
    DUMPS(6);
   // dump output to file
   //  DUMPS(16);
  }

  int npart = 0;
  
  // Fill event information
  GenEvent* evt = new GenEvent(1,1);

  GenVertex* prod_vtx = new GenVertex();
//  prod_vtx->add_particle_out( p );
  evt->add_vertex( prod_vtx );
                          
  // incoming beam e+
  GenParticle* p = new GenParticle( CLHEP::HepLorentzVector( MOMSET.p1[0],  MOMSET.p1[1], 
                                                      MOMSET.p1[2],  MOMSET.p1[3]),
                                                      -11, 3); 
  p->suggest_barcode( ++npart );
  prod_vtx->add_particle_in(p);

  // incoming beam e-
  p = new GenParticle( CLHEP::HepLorentzVector( MOMSET.q1[0],  MOMSET.q1[1], MOMSET.q1[2],  MOMSET.q1[3]),
                       11, 3); 
  p->suggest_barcode( ++npart );
  prod_vtx->add_particle_in(p);
  
  // scattered e+
  p = new GenParticle( CLHEP::HepLorentzVector( MOMSET.p2[0],  MOMSET.p2[1], 
                                         MOMSET.p2[2],  MOMSET.p2[3]),
                                         -11, 1); 
  p->suggest_barcode( ++npart );
  prod_vtx->add_particle_out(p);

  // scattered e-
  p = new GenParticle( CLHEP::HepLorentzVector( MOMSET.q2[0],  MOMSET.q2[1], MOMSET.q2[2],  MOMSET.q2[3]),
                       11, 1); 
  p->suggest_barcode( ++npart );
  prod_vtx->add_particle_out(p);

  int iphot=0;
  for (iphot=0; iphot<MOMSET.nphot; iphot++)
  {
    // gamma
    p = new GenParticle( CLHEP::HepLorentzVector( MOMSET.phot[0][iphot],  MOMSET.phot[1][iphot], 
                                           MOMSET.phot[2][iphot],  MOMSET.phot[3][iphot]),
                                           22, 1); 
    p->suggest_barcode( ++npart );
    prod_vtx->add_particle_out(p);
  }
  
  if( log.level() < MSG::INFO )
  {
    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;
    }
    else 
    {
      log << MSG::INFO << "McGenEventCol created and " << npart <<" particles stored in McGenEvent" << endreq;
    }
  }

  // retrieve event from Transient Store (Storegate)
/*  SmartDataPtr<McGenEventCol> McEvtColl(eventSvc(), "/Event/Gen");    
  if ( McEvtColl == 0 ) 
  {
    log << MSG::ERROR << "Could not retrieve McGenEventCollection" << endreq;
    return StatusCode::FAILURE;
  };
  
  McGenEventCol::iterator mcItr;
  for( mcItr = McEvtColl->begin(); mcItr != McEvtColl->end(); mcItr++ )   
  {
     HepMC::GenEvent* hEvt = (*mcItr)->getGenEvt();
    // MeVToGeV( hEvt );
//     callEvtGen( hEvt );
    // GeVToMeV( hEvt );
  };
*/
     
  return StatusCode::SUCCESS; 
}

StatusCode Bhlumi::finalize

(

)

Definition at line 339 of file Bhlumi.cxx.

References BHLUMI, Bes_Common::INFO, npar, and xpar.

{
  MsgStream log(messageService(), name());

  BHLUMI( 2,xpar,npar);

  log << MSG::INFO << "Bhlumi finalized" << endreq;

  return StatusCode::SUCCESS;
}

StatusCode Bhlumi::initialize

(

)

Definition at line 119 of file Bhlumi.cxx.

References BHLUMI, cos(), calibUtil::ERROR, Bes_Common::FATAL, IBesRndmGenSvc::GetEngine(), GLIMIT, Bes_Common::INFO, m_angleMode, m_cmEnergy, m_infraredCut, m_maxThetaAngle, m_minThetaAngle, max, min, npar, p_BesRndmGenSvc, BhlumiRandom::setRandomEngine(), Bes_Common::WARNING, and xpar.

                             {

  MsgStream log(messageService(), name());
  
  log << MSG::INFO << "Bhlumi 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("Bhlumi");                            
  engine->showStatus();
  BhlumiRandom::setRandomEngine(engine);      

  GLIMIT(50000);
                                                                                            




  int KeyGen =   3;  // ! Multiphoton Bhlumi
  int KeyRem =   1;  // ! No remooval of photons below epsCM
  KeyRem =   0;  // ! Remooval of photons below epsCM, Necessary for Z!!!
  int KeyWgt =   2;  // ! weighted events, with t generation down to zero
  KeyWgt =   0;  // ! WT=1 events, for detector simulation
  int KeyRnd =   1;  // ! RANMAR random numbers
  int KeyOpt =1000*KeyGen +100*KeyRem +10*KeyWgt +KeyRnd;
  int KeyPia =   0;  // ! Vacuum polarization OFF
  int KeyMod =   2;  //   ! Matrix element default version 4.x
  KeyPia =   2;  //   ! Vacuum polarization ON
  int  KeyZet =   0;  //   ! Z contribution OFF
  KeyZet =   1;  //   ! Z contribution ON
  int KeyRad =1000*KeyZet +10*KeyMod +KeyPia;
  npar[0]=   KeyOpt;
  npar[1]=   KeyRad;
  double CmsEne =  m_cmEnergy;  
  xpar[0]=   CmsEne;
  double th1,th2,thmin,thmax;
  if(m_angleMode==0){
    th1    =  m_minThetaAngle;  //       ! Detector range ThetaMin [rad]
    th2    =  m_maxThetaAngle;  //       ! Detector range ThetaMax [rad]
    thmin    = 0.7*th1;  //   ! thmin has to be lower  than th1
    thmax    = 2.0*th2;  //   ! thmax has to be higher than th2
    if(thmin<0.||thmax>3.1415926) {
      log << MSG::WARNING << "input angles exceed range (0~pi), so effect will be unprospectable" << endreq;
      return StatusCode::FAILURE;
    }
  }
  else if(m_angleMode==1){
    th1 = m_minThetaAngle*3.1415926/180.;
    th2 = m_maxThetaAngle*3.1415926/180.;
    // not multiply  0.7(2.0) coefficient while large angle
    thmin    = th1;
    thmax    = th2;
  }
  else if(m_angleMode==2){
    th1 = acos(max(m_minThetaAngle,m_maxThetaAngle));
    th2 = acos(min(m_minThetaAngle,m_maxThetaAngle));
    // not multiply  0.7(2.0) coefficient while large angle
    thmin    = th1;
    thmax    = th2;
  }
  else{
    log << MSG::FATAL << "unknown angle mode!" << endreq;
    return StatusCode::FAILURE;
  }
  if(thmin<0.||thmax>3.1415926) {
    log << MSG::FATAL << "input angles exceed range (0~pi), unprospectable" << endreq;
    return StatusCode::FAILURE;
  }
  else if(thmin>thmax) {
    log << MSG::FATAL << "thmin>thmax, unprospectable" << endreq;
    return StatusCode::FAILURE;
  }
  if(KeyWgt == 2) thmin=th1;  //  ! Because generation below th1 is on!!!
  xpar[1]=   CmsEne*CmsEne*(1-cos(thmin))/2;  //  ! TransMin [GeV**2]
  xpar[2]=   CmsEne*CmsEne*(1-cos(thmax))/2;  //  ! TransMax [GeV**2]
  xpar[3]=   m_infraredCut;  //        ! Infrared cut on photon energy

 // MARINI(m_initSeed[0], m_initSeed[1], m_initSeed[2]);

  BHLUMI(-1,xpar,npar);

  return StatusCode::SUCCESS;
}

---

Member Data Documentation

int Bhlumi::m_angleMode [private]

Definition at line 41 of file Bhlumi.h.

Referenced by initialize().

double Bhlumi::m_cmEnergy [private]

Definition at line 43 of file Bhlumi.h.

Referenced by initialize().

double Bhlumi::m_infraredCut [private]

Definition at line 43 of file Bhlumi.h.

Referenced by initialize().

std::vector<int> Bhlumi::m_initSeed [private]

Definition at line 45 of file Bhlumi.h.

double Bhlumi::m_maxThetaAngle [private]

Definition at line 43 of file Bhlumi.h.

Referenced by initialize().

double Bhlumi::m_minThetaAngle [private]

Definition at line 43 of file Bhlumi.h.

Referenced by initialize().

int Bhlumi::npar[100] [private]

Definition at line 36 of file Bhlumi.h.

Referenced by execute(), finalize(), and initialize().

IBesRndmGenSvc* Bhlumi::p_BesRndmGenSvc [private]

Definition at line 38 of file Bhlumi.h.

Referenced by initialize().

double Bhlumi::xpar[100] [private]

Definition at line 35 of file Bhlumi.h.

Referenced by execute(), finalize(), and initialize().

---