EvtGen Class Reference

#include <EvtGen.hh>

List of all members.

Public Member Functions
	EvtGen (const char *const decayName, const char *const pdtTableName, EvtRandomEngine *randomEngine=0, EvtAbsRadCorr *isrEngine=0)
void	generateDecay (EvtParticle *p)
void	generateDecay (int stdhepid, EvtVector4R P, EvtVector4R D, EvtStdHep *evtStdHep, EvtSpinDensity *spinDensity=0)
void	generateEvent (EvtParticle *p, HepLorentzVector D)
void	generateEvent (int stdhepid, HepLorentzVector P, HepLorentzVector D)
void	readUDecay (const char *const udecay_name)
	~EvtGen ()
Private Attributes
EvtPDL	_pdl

---

Constructor & Destructor Documentation

EvtGen::EvtGen (  const char *const   decayName, const char *const   pdtTableName, EvtRandomEngine *  randomEngine = 0, EvtAbsRadCorr *  isrEngine = 0 )

{ report(INFO,"EvtGen") << "Initializing EvtGen"<<endl; report(INFO,"EvtGen") << "Storing known decay models"<<endl; EvtModelReg dummy; report(INFO,"EvtGen") << "Main decay file name :"<<decayName<<endl; report(INFO,"EvtGen") << "PDT table file name :"<<pdtTableName<<endl; report(INFO,"EvtGen") << "Initializing RadCorr=PHOTOS"<<endl; if (isrEngine==0){ static EvtPHOTOS defaultRadCorrEngine; EvtRadCorr::setRadCorrEngine(&defaultRadCorrEngine); report(INFO,"EvtGen") <<"No RadCorr engine given in " <<"EvtGen::EvtGen constructor, " <<"will use default EvtPHOTOS."<<endl; } else{ EvtRadCorr::setRadCorrEngine(isrEngine); } _pdl.readPDT(pdtTableName); ifstream indec; EvtDecayTable::readDecayFile(decayName); if (randomEngine==0){ static EvtRandomEngine defaultRandomEngine; EvtRandom::setRandomEngine(&defaultRandomEngine); report(INFO,"EvtGen") <<"No random engine given in " <<"EvtGen::EvtGen constructor, " <<"will use default EvtRandomEngine."<<endl; } else{ EvtRandom::setRandomEngine(randomEngine); } report(INFO,"EvtGen") << "Done initializing EvtGen"<<endl; }

EvtGen::~EvtGen (   )

{ //This is a bit uggly, should not do anything //in a destructor. This will fail if EvtGen is made a static //because then this destructor might be called _after_ //the destructoin of objects that it depends on, e.g., EvtPDL. if (getenv("EVTINFO")){ EvtDecayTable::printSummary(); } }

---

Member Function Documentation

void EvtGen::generateDecay (  EvtParticle *  p  )

{ int times=0; do{ times+=1; EvtStatus::initRejectFlag(); p->decay(); //ok then finish. if ( EvtStatus::getRejectFlag()==0 ) { times=0; } else{ int ii; for (ii=0;ii<p->getNDaug();ii++){ EvtParticle *temp=p->getDaug(ii); temp->deleteTree(); } p->resetFirstOrNot(); p->resetNDaug(); } if ( times==10000) { report(ERROR,"EvtGen") << "Your event has been rejected 10000 times!"<<endl; report(ERROR,"EvtGen") << "Will now abort."<<endl; ::abort(); times=0; } } while (times); }

void EvtGen::generateDecay (  int  stdhepid, EvtVector4R  P, EvtVector4R  D, EvtStdHep *  evtStdHep, EvtSpinDensity *  spinDensity = 0 )

{ EvtParticle *p; if(spinDensity==0){ p=EvtParticleFactory::particleFactory(EvtPDL::evtIdFromStdHep(stdhepid),P); } else{ p=EvtParticleFactory::particleFactory(EvtPDL::evtIdFromStdHep(stdhepid), P,*spinDensity); } generateDecay(p); // p->Decay(); evtStdHep->init(); p->makeStdHep(*evtStdHep); evtStdHep->translate(D); p->deleteTree(); }

void EvtGen::generateEvent (  EvtParticle *  p, HepLorentzVector  D )

{ int i; static int nevent=0; nevent++; static EvtStdHep evtstdhep; // static EvtSecondary evtsecondary; int j; int istat; int partnum; double px,py,pz,e,m; double x,y,z,t; EvtVector4R p4,x4; generateDecay(root_part); // root_part->Decay(); int npart=0; EvtId list_of_stable[10]; EvtParticle* stable_parent[10]; list_of_stable[0]=EvtId(-1,-1); stable_parent[0]=0; evtstdhep.init(); // evtsecondary.init(); // root_part->makeStdHep(evtstdhep,evtsecondary,list_of_stable); root_part->makeStdHep(evtstdhep); //report(INFO,"EvtGen") << evtstdhep; //report(INFO,"EvtGen") << evtsecondary; npart=evtstdhep.getNPart(); for(i=0;i<evtstdhep.getNPart();i++){ j=i+1; int jmotherfirst=evtstdhep.getFirstMother(i)+1; int jmotherlast=evtstdhep.getLastMother(i)+1; int jdaugfirst=evtstdhep.getFirstDaughter(i)+1; int jdauglast=evtstdhep.getLastDaughter(i)+1; partnum=evtstdhep.getStdHepID(i); istat=evtstdhep.getIStat(i); p4=evtstdhep.getP4(i); x4=evtstdhep.getX4(i); px=p4.get(1); py=p4.get(2); pz=p4.get(3); e=p4.get(0); x=x4.get(1)+D.x(); y=x4.get(2)+D.y(); z=x4.get(3)+D.z(); t=x4.get(0)+D.t(); m=p4.mass(); begevtgenstore_(&j,&nevent,&npart,&istat, &partnum,&jmotherfirst,&jmotherlast, &jdaugfirst,&jdauglast, &px,&py,&pz,&e, &m,&x,&y,&z,&t); } }

void EvtGen::generateEvent (  int  stdhepid, HepLorentzVector  P, HepLorentzVector  D )

{ EvtParticle *root_part; EvtVectorParticle *vector_part; vector_part=new EvtVectorParticle; EvtVector4R p_init; p_init.set(P.t(),P.x(),P.y(),P.z()); vector_part->init(EvtPDL::evtIdFromStdHep(stdhepid),p_init); root_part=(EvtParticle *)vector_part; root_part->setVectorSpinDensity(); generateEvent(root_part,D); root_part->deleteTree(); }

void EvtGen::readUDecay (  const char *const   udecay_name  )

{ ifstream indec; if ( uDecayName[0] == 0) { report(INFO,"EvtGen") << "Is not reading a user decay file!"<<endl; } else{ indec.open(uDecayName); if (indec) { EvtDecayTable::readDecayFile(uDecayName); report(INFO,"EvtGen") << "Reading "<<uDecayName <<" to override decay table."<<endl; } else{ report(INFO,"EvtGen") << "Can not find UDECAY file '" <<uDecayName<<"'. Stopping"<<endl; ::abort(); } } }

---

Member Data Documentation

EvtPDL EvtGen::_pdl [private]

---

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

• /data0/mdestefa/bes3soft/Boss_6.5.5/dist/6.5.5/Generator/BesEvtGen/BesEvtGen-00-01-11/src/EvtGen/EvtGen.hh
• /data0/mdestefa/bes3soft/Boss_6.5.5/dist/6.5.5/Generator/BesEvtGen/BesEvtGen-00-01-11/src/EvtGen/EvtGen.cc

---