EvtPythia Class Reference

#include <EvtPythia.hh>

Inheritance diagram for EvtPythia:

Public Member Functions
	EvtPythia ()
virtual	~EvtPythia ()
void	getName (std::string &name)
EvtDecayBase *	clone ()
void	decay (EvtParticle *p)
std::string	commandName ()
void	command (std::string cmd)
void	init ()
void	initProbMax ()
Public Member Functions inherited from EvtDecayIncoherent
void	makeDecay (EvtParticle *p)
virtual	~EvtDecayIncoherent ()
void	setDaughterSpinDensity (int daughter)
int	isDaughterSpinDensitySet (int daughter)
Public Member Functions inherited from EvtDecayBase
double	getProbMax (double prob)
double	resetProbMax (double prob)
	EvtDecayBase ()
virtual	~EvtDecayBase ()
virtual bool	matchingDecay (const EvtDecayBase &other) const
EvtId	getParentId ()
double	getBranchingFraction ()
void	disableCheckQ ()
void	checkQ ()
int	getNDaug ()
EvtId *	getDaugs ()
EvtId	getDaug (int i)
int	getNArg ()
int	getPHOTOS ()
void	setPHOTOS ()
void	setVerbose ()
void	setSummary ()
double *	getArgs ()
std::string *	getArgsStr ()
double	getArg (int j)
std::string	getArgStr (int j)
std::string	getModelName ()
int	getDSum ()
int	summary ()
int	verbose ()
void	saveDecayInfo (EvtId ipar, int ndaug, EvtId *daug, int narg, std::vector< std::string > &args, std::string name, double brfr)
void	printSummary ()
void	setProbMax (double prbmx)
void	noProbMax ()
void	checkNArg (int a1, int a2=-1, int a3=-1, int a4=-1)
void	checkNDaug (int d1, int d2=-1)
void	checkSpinParent (EvtSpinType::spintype sp)
void	checkSpinDaughter (int d1, EvtSpinType::spintype sp)
virtual int	nRealDaughters ()

Static Public Member Functions
static void	pythiaInit (int f)
static void	pythiacont (double *, int *, int *, double *, double *, double *, double *)
Static Public Member Functions inherited from EvtDecayBase
static void	findMasses (EvtParticle *p, int ndaugs, EvtId daugs[10], double masses[10])
static void	findMass (EvtParticle *p)
static double	findMaxMass (EvtParticle *p)

Additional Inherited Members
Protected Member Functions inherited from EvtDecayBase
bool	daugsDecayedByParentModel ()
Protected Attributes inherited from EvtDecayBase
bool	_daugsDecayedByParentModel

Detailed Description

Definition at line 32 of file EvtPythia.hh.

Constructor & Destructor Documentation

EvtPythia()

EvtPythia::EvtPythia

(

)

Definition at line 94 of file EvtPythia.cc.

{}

Referenced by clone().

~EvtPythia()

EvtPythia::~EvtPythia ( )

virtual

Definition at line 96 of file EvtPythia.cc.

                     {
 
 
  int i;
 
 
  //the deletion of commands is really uggly!
 
  if (njetsetdecays==0) {
    delete [] commands;
    commands=0;
    return;
  }
 
  for(i=0;i<njetsetdecays;i++){
    if (jetsetdecays[i]==this){
      jetsetdecays[i]=jetsetdecays[njetsetdecays-1];
      njetsetdecays--;
      if (njetsetdecays==0) {
    delete [] commands;
    commands=0;
      }
      return;
    }
  }
 
  report(ERROR,"EvtGen") << "Error in destroying Pythia model!"<<endl;
 
}

Member Function Documentation

clone()

EvtDecayBase * EvtPythia::clone ( )

virtual

Implements EvtDecayBase.

Definition at line 133 of file EvtPythia.cc.

                              {
 
  return new EvtPythia;
 
}

command()

void EvtPythia::command ( std::string cmd )

virtual

Reimplemented from EvtDecayBase.

Definition at line 177 of file EvtPythia.cc.

                                    {
 
  if (ncommand==lcommand){
 
    lcommand=10+2*lcommand;
 
    std::string* newcommands=new std::string[lcommand];
    
    int i;
 
    for(i=0;i<ncommand;i++){
      newcommands[i]=commands[i];
    }
    
    delete [] commands;
 
    commands=newcommands;
 
  }
 
  commands[ncommand]=cmd;
 
  ncommand++;
 
}

commandName()

std::string EvtPythia::commandName ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 170 of file EvtPythia.cc.

                                {
 
  return std::string("JetSetPar");
  
}

decay()

void EvtPythia::decay ( EvtParticle * p )

virtual

Implements EvtDecayBase.

Definition at line 211 of file EvtPythia.cc.

                                    {
  
  //added by Lange Jan4,2000
  static EvtId STRNG=EvtPDL::getId("string");
 
  int istdheppar=EvtPDL::getStdHep(p->getId());
 
  if (pycomp_(&istdheppar)==0){
    report(ERROR,"EvtGen") << "Pythia can not decay:"
      <<EvtPDL::name(p->getId()).c_str()<<endl;
     return;  
  }
 
 
  double mp=p->mass();
 
  EvtVector4R p4[20];
  
  int i,more;
  int ip=EvtPDL::getStdHep(p->getId());
 
  int ndaugjs;
  int kf[100];
  EvtId evtnumstable[100],evtnumparton[100];
  int stableindex[100],partonindex[100];
  int numstable;
  int numparton;
  int km[100];
  EvtId type[MAX_DAUG];
  
  cbbeam_.maximum = mp; //pingrg
  if(mp==0){std::cout<<"Particle "<<EvtPDL::name(p->getId())<<" has zero mass"<<std::endl;abort();}
  pythiaInit(0);
  
  double px[100],py[100],pz[100],e[100];
  if ( p->getNDaug() != 0 ) { p->deleteDaughters(true);}
  
  int count=0;
 
  do{
 
    pythiadec_(&ip,&mp,&ndaugjs,kf,km,px,py,pz,e);
 
    
    numstable=0;
    numparton=0;
 
    for(i=0;i<ndaugjs;i++){
      // std::cout<<"EvtPDL::evtIdFromStdHep(kf[i]),i= "<<i<<" "<<EvtPDL::evtIdFromStdHep(kf[i])<<std::endl;
      if (EvtPDL::evtIdFromStdHep(kf[i])==EvtId(-1,-1)) {
    report(ERROR,"EvtGen") << "Pythia returned particle:"<<kf[i]<<endl;
    report(ERROR,"EvtGen") << "This can not be translated to evt number"<<endl;
    report(ERROR,"EvtGen") << "and the decay will be rejected!"<<endl;
    report(ERROR,"EvtGen") << "The decay was of particle:"<<ip<<endl;
    int i=1;
    pylist_(i);
      }
 
      //sort out the partons
      if (abs(kf[i])<=6||kf[i]==21){
    partonindex[numparton]=i;
    evtnumparton[numparton]=EvtPDL::evtIdFromStdHep(kf[i]);
    numparton++;
      }
      else{
    stableindex[numstable]=i;
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(kf[i]); 
    numstable++;
      }
      
      
      // have to protect against negative mass^2 for massless particles
      // i.e. neutrinos and photons.
      // this is uggly but I need to fix it right now....
      
      if (px[i]*px[i]+py[i]*py[i]+pz[i]*pz[i]>=e[i]*e[i]){
    
        e[i]=sqrt(px[i]*px[i]+py[i]*py[i]+pz[i]*pz[i])+0.0000000000001;
    
      }
      
      p4[i].set(e[i],px[i],py[i],pz[i]);
      
      
    }
    
  int channel= EvtDecayTable::inChannelList(p->getId(),numstable,evtnumstable);
   
    more=(channel!=-1);
   
   
   count++;
   
  }while( more && (count<10000) );
 
 
  if (count>9999) {
    report(INFO,"EvtGen") << "Too many loops in EvtPythia!!!"<<endl;
    report(INFO,"EvtGen") << "Parent:"<<EvtPDL::name(getParentId()).c_str()<<endl;
    for(i=0;i<numstable;i++){
      report(INFO,"EvtGen") << "Daug("<<i<<")"<<EvtPDL::name(evtnumstable[i]).c_str()<<endl;
    }
    
  }
  
  
  
  if (numparton==0){
    
    p->makeDaughters(numstable,evtnumstable);
    
    for(i=0;i<numstable;i++){
      p->getDaug(i)->init(evtnumstable[i],p4[stableindex[i]]);
    }
    
    fixPolarizations(p);
    
    return ;
   
  }
  else{
    
    //have partons in JETSET
    
    EvtVector4R p4string(0.0,0.0,0.0,0.0);
    
    for(i=0;i<numparton;i++){
      p4string+=p4[partonindex[i]];
    }
    
    int nprimary=1;
    type[0]=STRNG;
    for(i=0;i<numstable;i++){
      if (km[stableindex[i]]==0){
    type[nprimary++]=evtnumstable[i];
      }
    }
    
    p->makeDaughters(nprimary,type);
    
    p->getDaug(0)->init(STRNG,p4string);
    
    EvtVector4R p4partons[10];
    
    for(i=0;i<numparton;i++){
      p4partons[i]=p4[partonindex[i]];
    }
    
    ((EvtStringParticle*)p->getDaug(0))->initPartons(numparton,p4partons,evtnumparton);
    
    
    
    nprimary=1;
    
    for(i=0;i<numstable;i++){
      
      if (km[stableindex[i]]==0){
    p->getDaug(nprimary++)->init(evtnumstable[i],p4[stableindex[i]]);
      }
    }
    
    
    int nsecond=0;
    for(i=0;i<numstable;i++){
      if (km[stableindex[i]]!=0){
    type[nsecond++]=evtnumstable[i];
      }
    }
    
    
    p->getDaug(0)->makeDaughters(nsecond,type);
    
    nsecond=0;
    for(i=0;i<numstable;i++){
      if (km[stableindex[i]]!=0){
    p4[stableindex[i]]=boostTo(p4[stableindex[i]],p4string);
    p->getDaug(0)->getDaug(nsecond)->init(evtnumstable[i],p4[stableindex[i]]);
    p->getDaug(0)->getDaug(nsecond)->setDiagonalSpinDensity();
    p->getDaug(0)->getDaug(nsecond)->decay();
    nsecond++;
      }
    }
    
    fixPolarizations(p);
    
    return ;
    
  }
 
}

getName()

void EvtPythia::getName ( std::string & name )

virtual

Implements EvtDecayBase.

Definition at line 127 of file EvtPythia.cc.

                                            {
 
  model_name="PYTHIA";     
 
}

init()

void EvtPythia::init ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 147 of file EvtPythia.cc.

                    {
 
  checkNArg(1);
 
 
  if (getParentId().isAlias()){
 
    report(ERROR,"EvtGen") << "EvtPythia finds that you are decaying the"<<endl
                           << " aliased particle "
               << EvtPDL::name(getParentId()).c_str()
               << " with the Pythia model"<<endl
               << " this does not work, please modify decay table."
               << endl;
    report(ERROR,"EvtGen") << "Will terminate execution!"<<endl;
    ::abort();
 
  }
 
  store(this);
 
}

initProbMax()

void EvtPythia::initProbMax ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 140 of file EvtPythia.cc.

                           {
 
  noProbMax();
 
}

pythiacont()

void EvtPythia::pythiacont ( double * energy, int * ndaugjs, int * kf, double * px, double * py, double * pz, double * e )

static

Definition at line 203 of file EvtPythia.cc.

{
  pycontinuum_(energy,ndaugjs,kf,px,py,pz,e);
}

Referenced by EvtPycont::decay(), and EvtPyGaGa::decay().

pythiaInit()

void EvtPythia::pythiaInit ( int f )

static

Definition at line 1051 of file EvtPythia.cc.

                                   {
  
  static int first=1;
  if (first){
    
    first=0;
    
    report(INFO,"EvtGen") << "Will initialize Pythia."<<endl;
    for(int i=0;i<ncommand;i++)
      pygive_(commands[i].c_str(),strlen(commands[i].c_str()));
    
    char fname[200];
    
    char hostBuffer[100];
    
    if ( gethostname( hostBuffer, 100 ) != 0 ){
      report(ERROR,"EvtGen") << " couldn't get hostname." << endl;
      strncpy( hostBuffer, "hostnameNotFound", 100 );
    }
    
    char pid[100];
    
    int thePid=getpid();
    
    if ( sprintf( pid, "%d", thePid ) == 0 ){
      report(ERROR,"EvtGen") << " couldn't get process ID." << endl;
      strncpy( pid, "666", 100 );
    }
    
    strcpy(fname,"jet.d-");
    strcat(fname,hostBuffer);
    strcat(fname,"-");
    strcat(fname,pid);
    
    MakePythiaFile(fname);
    evtpythiainit_(fname,strlen(fname));
    initpythia_(&dummy);
    
    if (0==getenv("EVTSAVEJETD")){
      char delcmd[300];
      strcpy(delcmd,"rm -f ");
      strcat(delcmd,fname); 
      system(delcmd);
    }
    
    report(INFO,"EvtGen") << "Done initializing Pythia."<<endl;
    
  }
 
}

Referenced by EvtPycont::decay(), EvtPyGaGa::decay(), and decay().

---

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

• EvtPythia.hh
• EvtPythia.cc