EvtGen Class Reference

#include <EvtGen.hh>

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

Detailed Description

Definition at line 46 of file EvtGen.hh.

Constructor & Destructor Documentation

EvtGen()

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

Definition at line 77 of file EvtGen.cc.

                                    {
 
 
  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::~EvtGen

(

)

Definition at line 64 of file EvtGen.cc.

               {
 
  //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

generateDecay() [1/2]

void EvtGen::generateDecay

(

EvtParticle *

p

)

Definition at line 183 of file EvtGen.cc.

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

generateDecay() [2/2]

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

Definition at line 152 of file EvtGen.cc.

                                            {
 
 
  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();
 
 
}

Referenced by generateDecay(), and generateEvent().

generateEvent() [1/2]

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

Definition at line 240 of file EvtGen.cc.

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

generateEvent() [2/2]

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

Definition at line 218 of file EvtGen.cc.

                                                                            {
 
  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();  
 
}

Referenced by generateEvent().

readUDecay()

void EvtGen::readUDecay

(

const char *const

udecay_name

)

Definition at line 126 of file EvtGen.cc.

                                                   {
 
  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();
    }
  }
  
}

Referenced by EvtDecay::initialize().

---

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

• EvtGen.hh
• EvtGen.cc