EvtPhokhara_Lambda Class Reference

#include <EvtPhokhara_Lambda.hh>

Inheritance diagram for EvtPhokhara_Lambda:

Public Member Functions
	EvtPhokhara_Lambda ()
virtual	~EvtPhokhara_Lambda ()
void	getName (std::string &name)
EvtDecayBase *	clone ()
void	decay (EvtParticle *p)
std::string	commandName ()
void	command (std::string cmd)
void	init ()
void	init_mode (EvtParticle *p)
void	init_evt (EvtParticle *p)
void	initProbMax ()
int	getTotalEvt ()
void	PhokharaInit (int dummy)
void	ExclusiveDecay (EvtParticle *p)
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
virtual void	getName (std::string &name)=0
virtual void	decay (EvtParticle *p)=0
virtual void	makeDecay (EvtParticle *p)=0
virtual EvtDecayBase *	clone ()=0
virtual void	init ()
virtual void	initProbMax ()
virtual std::string	commandName ()
virtual void	command (std::string cmd)
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 ()

Additional Inherited Members
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)
Protected Member Functions inherited from EvtDecayBase
bool	daugsDecayedByParentModel ()
Protected Attributes inherited from EvtDecayBase
bool	_daugsDecayedByParentModel

Detailed Description

Definition at line 33 of file EvtPhokhara_Lambda.hh.

Constructor & Destructor Documentation

EvtPhokhara_Lambda()

EvtPhokhara_Lambda::EvtPhokhara_Lambda

(

)

Definition at line 68 of file EvtPhokhara_Lambda.cc.

{}

Referenced by clone().

~EvtPhokhara_Lambda()

EvtPhokhara_Lambda::~EvtPhokhara_Lambda ( )

virtual

Definition at line 69 of file EvtPhokhara_Lambda.cc.

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

Member Function Documentation

clone()

EvtDecayBase * EvtPhokhara_Lambda::clone ( )

virtual

Implements EvtDecayBase.

Definition at line 102 of file EvtPhokhara_Lambda.cc.

                                       {
 
  return new EvtPhokhara_Lambda;
 
}

command()

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

virtual

Reimplemented from EvtDecayBase.

Definition at line 160 of file EvtPhokhara_Lambda.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 EvtPhokhara_Lambda::commandName ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 154 of file EvtPhokhara_Lambda.cc.

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

decay()

void EvtPhokhara_Lambda::decay ( EvtParticle *  p )

virtual

Implements EvtDecayBase.

Definition at line 188 of file EvtPhokhara_Lambda.cc.

                                             {
  bool debug=0;
  EvtId myvpho=EvtPDL::getId("vpho");
  if(p->getId()!=myvpho) {std::cout<<"Parent particle is required to be vpho for Phokhara model"<<std::endl;abort();}
  if(nevtgen==0) {init_mode(p);
  std::cout<<"PHOKHARA :  Lambda anti-Lambda  mode "<<std::endl;
  }
  else{init_evt(p);
  if(debug)cout<<"init_evt works!"<<endl;
  }
  //else{init_mode(p);}
  if(debug) {
      std::cout<<"flags :ph0,nlo,pion,fsr,fsrnlo,ivac,FF_pion,f0_model,FF_kaon,narr_res,FF_pp,chi_sw,chi_pion,FF_Pgg,nlo2 "<<std::endl;
      std::cout<<"= "<<flags_.ph0<<","<<flags_.nlo<<","<<flags_.pion<<","<<flags_.fsr<<","<<flags_.fsrnlo<<","<<flags_.ivac<<","<<flags_.FF_pion<<","<<flags_.f0_model<<","<<flags_.FF_kaon<<","<<flags_.narr_res<<","<<flags_.FF_pp<<","<<flags_.chi_sw<<","<<flags_.be_r<<","<<flags_.FF_Pgg<<","<<flags_.nlo2<<std::endl;
      std::cout<<"ctes: Sp = "<<ctes_.Sp<<std::endl;
      std::cout<<"cuts : w,q2min,q2_min_c,gmin,phot1cut,phot2cut,pi1cut,pi2cut"<<std::endl<<"= "<<cuts_.w<<","<<cuts_.q2min<<","<<cuts_.q2_min_c<<","<<cuts_.gmin<<","<<cuts_.phot1cut<<","<<cuts_.phot2cut<<","<<cuts_.pi1cut<<","<<cuts_.pi2cut<<std::endl;
  }
 
 
 
  int istdheppar=EvtPDL::getStdHep(p->getId());
  int ntrials = 0;
  int tr_old[3];
  tr_old[0] = (int)maxima_.tr[0];
  tr_old[1] = (int)maxima_.tr[1];
  tr_old[2] = (int)maxima_.tr[2];
 
  while( ntrials < 1000000)
    {
      ievent++;
      RANLXDF(Ar_r,1);
      Ar[1] = Ar_r[0];
      
     if (Ar[1] <= (maxima_.Mmax[0]/(maxima_.Mmax[0]+maxima_.Mmax[1]+maxima_.Mmax[2]))) {
        maxima_.count[0] = maxima_.count[0]+1.0;
        GEN_0PH(2,qqmin,ctes_.Sp,cos3min,cos3max);
       }else 
    if (Ar[1] <= ( (maxima_.Mmax[0]+maxima_.Mmax[1])/(maxima_.Mmax[0]+maxima_.Mmax[1]+maxima_.Mmax[2]))) {
        maxima_.count[1] = maxima_.count[1]+1.0;
        GEN_1PH(2,qqmin,qqmax,cos1min,cos1max,cos3min,cos3max);
       }
       else {
         maxima_.count[2] = maxima_.count[2]+1.0;
         GEN_2PH(2,qqmin,cos1min,cos1max,cos2min,cos2max,cos3min,cos3max);
       }
      
      if( ((int)maxima_.tr[0]+(int)maxima_.tr[1]+(int)maxima_.tr[2]) > (tr_old[0]+tr_old[1]+tr_old[2]) ) // event accepted after cuts
        {
          goto storedEvents;
        }
      ntrials ++;
    }
  std::cout <<"FATAL: Could not satisfy cuts after " << ntrials << "trials. Terminate." <<std::endl;
  //----
 storedEvents:
  int more=0;
  int numstable=0;
  int numparton=0;
  EvtId evtnumstable[100];//
  EvtVector4R p4[20];
 
  // except ISR photos, products depending on channel
  if (flags_.pion == 0) { // mu+ mu-
    // mu+
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-13);
    p4[numstable].set(ctes_.momenta[0][5],ctes_.momenta[1][5], ctes_.momenta[2][5], ctes_.momenta[3][5]);     
    numstable++;
    // mu -
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(13);
    p4[numstable].set(ctes_.momenta[0][6],ctes_.momenta[1][6], ctes_.momenta[2][6], ctes_.momenta[3][6]);     
    numstable++;
  }
  if (flags_.pion == 1) { // pi+ pi-
    // pi+
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(211);
    p4[numstable].set(ctes_.momenta[0][5],ctes_.momenta[1][5], ctes_.momenta[2][5], ctes_.momenta[3][5]);     
    numstable++;
    // pi -
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-211);
    p4[numstable].set(ctes_.momenta[0][6],ctes_.momenta[1][6], ctes_.momenta[2][6], ctes_.momenta[3][6]);     
    numstable++;
  }
  if (flags_.pion == 2) { // pi+ pi-2pi0
    // pi0
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(111);
    p4[numstable].set(ctes_.momenta[0][5],ctes_.momenta[1][5], ctes_.momenta[2][5], ctes_.momenta[3][5]);     
    numstable++;
    // pi0
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(111);
    p4[numstable].set(ctes_.momenta[0][6],ctes_.momenta[1][6], ctes_.momenta[2][6], ctes_.momenta[3][6]);     
    numstable++;
    // pi-
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-211);
    p4[numstable].set(ctes_.momenta[0][7],ctes_.momenta[1][7], ctes_.momenta[2][7], ctes_.momenta[3][7]);     
    numstable++;
    // pi +
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(211);
    p4[numstable].set(ctes_.momenta[0][8],ctes_.momenta[1][8], ctes_.momenta[2][8], ctes_.momenta[3][8]);     
    numstable++;
  }   
  if (flags_.pion == 3) { // 2(pi+ pi-)
    // pi+
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(211);
    p4[numstable].set(ctes_.momenta[0][5],ctes_.momenta[1][5], ctes_.momenta[2][5], ctes_.momenta[3][5]);     
    numstable++;
    // pi-
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-211);
    p4[numstable].set(ctes_.momenta[0][6],ctes_.momenta[1][6], ctes_.momenta[2][6], ctes_.momenta[3][6]);     
    numstable++;
    // pi+
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-211);
    p4[numstable].set(ctes_.momenta[0][7],ctes_.momenta[1][7], ctes_.momenta[2][7], ctes_.momenta[3][7]);     
    numstable++;
    // pi -
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(211);
    p4[numstable].set(ctes_.momenta[0][8],ctes_.momenta[1][8], ctes_.momenta[2][8], ctes_.momenta[3][8]);     
    numstable++;
  }   
  if (flags_.pion == 4) { // ppbar
    // pbar
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-2212);
    p4[numstable].set(ctes_.momenta[0][5],ctes_.momenta[1][5], ctes_.momenta[2][5], ctes_.momenta[3][5]);     
    numstable++;
    // p
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(2212);
    p4[numstable].set(ctes_.momenta[0][6],ctes_.momenta[1][6], ctes_.momenta[2][6], ctes_.momenta[3][6]);     
    numstable++;
  }
  if (flags_.pion == 5) { // nnbar
    // pbar
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-2112);
    p4[numstable].set(ctes_.momenta[0][5],ctes_.momenta[1][5], ctes_.momenta[2][5], ctes_.momenta[3][5]);     
    numstable++;
    // p
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(2112);
    p4[numstable].set(ctes_.momenta[0][6],ctes_.momenta[1][6], ctes_.momenta[2][6], ctes_.momenta[3][6]);     
    numstable++;
  }
  if (flags_.pion == 6) { // K+ K-
    // K+
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(321);
    p4[numstable].set(ctes_.momenta[0][5],ctes_.momenta[1][5], ctes_.momenta[2][5], ctes_.momenta[3][5]);     
    numstable++;
    // K -
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-321);
    p4[numstable].set(ctes_.momenta[0][6],ctes_.momenta[1][6], ctes_.momenta[2][6], ctes_.momenta[3][6]);     
    numstable++;
  }
  if (flags_.pion == 7) { // K0K0bar
    // Kbar
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(311);
    p4[numstable].set(ctes_.momenta[0][5],ctes_.momenta[1][5], ctes_.momenta[2][5], ctes_.momenta[3][5]);     
    numstable++;
    // K0
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-311);
    p4[numstable].set(ctes_.momenta[0][6],ctes_.momenta[1][6], ctes_.momenta[2][6], ctes_.momenta[3][6]);     
    numstable++;
  }
  if (flags_.pion == 8) { // pi+ pi-pi0
    // pi+
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(211);
    p4[numstable].set(ctes_.momenta[0][5],ctes_.momenta[1][5], ctes_.momenta[2][5], ctes_.momenta[3][5]);     
    numstable++;
    // pi-
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-211);
    p4[numstable].set(ctes_.momenta[0][6],ctes_.momenta[1][6], ctes_.momenta[2][6], ctes_.momenta[3][6]);     
    numstable++;
    // pi0
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(111);
    p4[numstable].set(ctes_.momenta[0][7],ctes_.momenta[1][7], ctes_.momenta[2][7], ctes_.momenta[3][7]);     
    numstable++;
  } 
  /*if (flags_.pion == 9) { //Lambda Lambdabar-> pi+ pi- ppbar
    // pi+
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(211);
    p4[numstable].set(ctes_.momenta[0][7],ctes_.momenta[1][7], ctes_.momenta[2][7], ctes_.momenta[3][7]);     
    numstable++;
    // pbar
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-2212);
    p4[numstable].set(ctes_.momenta[0][8],ctes_.momenta[1][8], ctes_.momenta[2][8], ctes_.momenta[3][8]);     
    numstable++;
    // pi-
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-211);
    p4[numstable].set(ctes_.momenta[0][9],ctes_.momenta[1][9], ctes_.momenta[2][9], ctes_.momenta[3][9]);     
    numstable++; 
    // p
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(2212);
    p4[numstable].set(ctes_.momenta[0][10],ctes_.momenta[1][10], ctes_.momenta[2][10], ctes_.momenta[3][10]);     
    numstable++;
  }*/
  int LambdaDaus=0;
  EvtId LambdaMode[100];//
  EvtVector4R Lambdap4[20];
  EvtVector4R ZERO(0,0,0,0);
  EvtVector4R tmp(0,0,0,0);
 
  if (flags_.pion == 9) { //Lambda Lambdabar-> pi+ pi- ppbar
    // anti-Lambda0
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-3122);
    p4[numstable].set(ctes_.momenta[0][5],ctes_.momenta[1][5], ctes_.momenta[2][5], ctes_.momenta[3][5]);     
    if (debug)std::cout<<"Phokhara_Lambda: anti-Lambda0 p4[numstable] = "<<p4[numstable]<<std::endl;
    numstable++;
    // Lambda0
    evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(3122);
    p4[numstable].set(ctes_.momenta[0][6],ctes_.momenta[1][6], ctes_.momenta[2][6], ctes_.momenta[3][6]);     
    if (debug) std::cout<<"Phokhara_Lambda: Lambda0 p4[numstable] = "<<p4[numstable]<<std::endl;
    numstable++;
    /////Lambda Decay 
    // pi+
    LambdaMode[LambdaDaus]=EvtPDL::evtIdFromStdHep(211);
    Lambdap4[LambdaDaus].set(ctes_.momenta[0][7],ctes_.momenta[1][7], ctes_.momenta[2][7], ctes_.momenta[3][7]);     
    if (debug) std::cout<<"Phokhara_Lambda: Pi+ p4[numstable] = "<<Lambdap4[LambdaDaus]<<std::endl;
    ////////boost pi+ to antiLambda0 frame
    //tmp = ZERO-p4[0];wrong calculation
    tmp.set(p4[0].get(0),-p4[0].get(1),-p4[0].get(2),-p4[0].get(3));
    Lambdap4[LambdaDaus] = boostTo(Lambdap4[LambdaDaus],tmp);
    if (debug) std::cout<<"Phokhara_Lambda:Boosted Pi+ p4[numstable] = "<<Lambdap4[LambdaDaus]<<std::endl;
    LambdaDaus++;
    // pbar
    LambdaMode[LambdaDaus]=EvtPDL::evtIdFromStdHep(-2212);
    Lambdap4[LambdaDaus].set(ctes_.momenta[0][8],ctes_.momenta[1][8], ctes_.momenta[2][8], ctes_.momenta[3][8]);     
    if (debug) std::cout<<"Phokhara_Lambda: pbar p4[numstable] = "<<Lambdap4[LambdaDaus]<<std::endl;
    ////////boost pbar to antiLambda0 frame
    //tmp = ZERO-p4[0];
    //tmp = p4[0];
    tmp.set(p4[0].get(0),-p4[0].get(1),-p4[0].get(2),-p4[0].get(3));
    Lambdap4[LambdaDaus] = boostTo(Lambdap4[LambdaDaus],tmp);
    if (debug) std::cout<<"Phokhara_Lambda:Boosted pbar p4[numstable] = "<<Lambdap4[LambdaDaus]<<std::endl;
    LambdaDaus++;
    // pi-
    LambdaMode[LambdaDaus]=EvtPDL::evtIdFromStdHep(-211);
    Lambdap4[LambdaDaus].set(ctes_.momenta[0][9],ctes_.momenta[1][9], ctes_.momenta[2][9], ctes_.momenta[3][9]);     
    if (debug) std::cout<<"Phokhara_Lambda: pi- p4[numstable] = "<<Lambdap4[LambdaDaus]<<std::endl;
    ////////boost pi- to Lambda0 frame
    //tmp = ZERO-p4[1];
    //tmp = p4[1];
    tmp.set(p4[1].get(0),-p4[1].get(1),-p4[1].get(2),-p4[1].get(3));
    Lambdap4[LambdaDaus] = boostTo(Lambdap4[LambdaDaus],tmp);
    if (debug) std::cout<<"Phokhara_Lambda:Boosted pi- p4[numstable] = "<<Lambdap4[LambdaDaus]<<std::endl;
    LambdaDaus++; 
    // p
    LambdaMode[LambdaDaus]=EvtPDL::evtIdFromStdHep(2212);
    Lambdap4[LambdaDaus].set(ctes_.momenta[0][10],ctes_.momenta[1][10], ctes_.momenta[2][10], ctes_.momenta[3][10]);     
    if (debug) std::cout<<"Phokhara_Lambda: p p4[numstable] = "<<Lambdap4[LambdaDaus]<<std::endl;
    ////////boost p+ to Lambda0 frame
    //tmp = ZERO-p4[1];
    //tmp = p4[1];
    tmp.set(p4[1].get(0),-p4[1].get(1),-p4[1].get(2),-p4[1].get(3));
    Lambdap4[LambdaDaus] = boostTo(Lambdap4[LambdaDaus],tmp);
    if (debug) std::cout<<"Phokhara_Lambda:Boosted p p4[numstable] = "<<Lambdap4[LambdaDaus]<<std::endl;
    LambdaDaus++;
  }
  
 // ISR gamma      
  evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(22);
  p4[numstable].set(ctes_.momenta[0][2],ctes_.momenta[1][2], ctes_.momenta[2][2], ctes_.momenta[3][2]);     
  if (debug) std::cout<<"Phokhara_Lambda: first gamma p4[numstable] = "<<p4[numstable]<<std::endl;
  numstable++;
  if( ctes_.momenta[0][3] != 0 ) // second photon exists
    {
      evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(22);
      p4[numstable].set(ctes_.momenta[0][3],ctes_.momenta[1][3], ctes_.momenta[2][3], ctes_.momenta[3][3]);     
      if (debug) std::cout<<"Phokhara_Lambda: second gamma p4[numstable] = "<<p4[numstable]<<std::endl;
      numstable++;
    }
  
  int channel=EvtDecayTable::inChannelList(p->getId(),numstable,evtnumstable);
  more=(channel!=-1);
  if(more) {std::cout<<"EvtPhokhara_Lambda:Existence of mode "<<channel<<" in exclusive decay list has the same final state as this one"<<std::endl;abort(); }
  
  p->makeDaughters(numstable,evtnumstable);
  //double weight = p->initializePhaseSpace(getNDaug(),getDaugs());
  if(EvtParticle::_NextLevelDauNum==0){
    //std::cout<<"NextLevelDauNum==0,LambdaDaus = "<<LambdaDaus<<std::endl;
    EvtParticle::_NextLevelDauNum = LambdaDaus;
    for(int i = 0; i< LambdaDaus;i++){
      EvtParticle::_NextLevelId[i] = LambdaMode[i];
      EvtParticle::_NextLevelP4[i] = Lambdap4[i];
    }
  }
  int ndaugFound=0;
  EvtVector4R SUMP4(0,0,0,0);
  for(int i=0;i<numstable;i++){
    p->getDaug(i)->init(evtnumstable[i],p4[i]);
    ndaugFound++;
  }
  if ( ndaugFound == 0 ) {
    report(ERROR,"EvtGen") << "Phokhara has failed to do a decay ";
    report(ERROR,"EvtGen") << EvtPDL::name(p->getId()).c_str() << " " << p->mass()<<endl;
    assert(0);
  }
  if (debug) std::cout<<"EvtPhokhara_Lambda SUMMARY: part p4"<<p->getP4Lab()<<std::endl;
  if (debug) std::cout<<"EvtPhokhara_Lambda SUMMARY: Daug0 p4"<<p->getDaug(0)->getP4Lab()<<std::endl;
  if (debug) std::cout<<"EvtPhokhara_Lambda SUMMARY: Daug1 p4"<<p->getDaug(1)->getP4Lab()<<std::endl;
  if (debug) std::cout<<"EvtPhokhara_Lambda SUMMARY: Daug2 p4"<<p->getDaug(2)->getP4Lab()<<std::endl;
 
  nevtgen++;  
  return ;
  
}

ExclusiveDecay()

void EvtPhokhara_Lambda::ExclusiveDecay

(

EvtParticle *

p

)

getName()

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

virtual

Implements EvtDecayBase.

Definition at line 96 of file EvtPhokhara_Lambda.cc.

                                                     {
 
  model_name="PHOKHARA_LAMBDA";     
 
}

getTotalEvt()

int EvtPhokhara_Lambda::getTotalEvt ( )

inline

Definition at line 52 of file EvtPhokhara_Lambda.hh.

{return nevt;}

init()

void EvtPhokhara_Lambda::init ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 127 of file EvtPhokhara_Lambda.cc.

                             {
  checkNArg(0);
 
  std::string locvp=getenv("BESEVTGENROOT");
  system("cat $BESEVTGENROOT/share/phokhara_10.0.param>phokhara_10.0.param");
  system("cat $BESEVTGENROOT/share/phokhara_10.0.fferr>phokhara_10.0.fferr");
  system("cat $BESEVTGENROOT/share/phokhara_10.0.ffwarn>phokhara_10.0.ffwarn");
 
  
  if (getParentId().isAlias()){
    
    report(ERROR,"EvtGen") << "EvtPhokhara finds that you are decaying the"<<endl
               << " aliased particle "
               << EvtPDL::name(getParentId()).c_str()
               << " with the Phokhara model"<<endl
               << " this does not work, please modify decay table."
               << endl;
    report(ERROR,"EvtGen") << "Will terminate execution!"<<endl;
    ::abort();
 
  }
 
  store(this);
 
}

init_evt()

void EvtPhokhara_Lambda::init_evt

(

EvtParticle *

p

)

Definition at line 525 of file EvtPhokhara_Lambda.cc.

                                               {
 m_pion=9; 
 // mu+mu-(0),pi+pi-(1),2pi0pi+pi-(2),
 // 2pi+2pi-(3),ppbar(4),nnbar(5),
 // K+K-(6),K0K0bar(7),pi+pi-pi0(8), 
 // Lamb Lambbar->pi-pi+ppbar(9) 
#include "Phokhara_init_evt.txt"
}

Referenced by decay().

init_mode()

void EvtPhokhara_Lambda::init_mode

(

EvtParticle *

p

)

Definition at line 116 of file EvtPhokhara_Lambda.cc.

                                                {
 m_pion=9; 
 // mu+mu-(0),pi+pi-(1),2pi0pi+pi-(2),
 // 2pi+2pi-(3),ppbar(4),nnbar(5),
 // K+K-(6),K0K0bar(7),pi+pi-pi0(8), 
 // Lamb Lambbar->pi-pi+ppbar(9) 
#include "Phokhara_init_mode.txt"
}

Referenced by decay().

initProbMax()

void EvtPhokhara_Lambda::initProbMax ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 109 of file EvtPhokhara_Lambda.cc.

                                    {
 
  noProbMax();
 
}

PhokharaInit()

void EvtPhokhara_Lambda::PhokharaInit

(

int

dummy

)

Definition at line 512 of file EvtPhokhara_Lambda.cc.

                                              {
  static int first=1;
  if (first){
    
    first=0;
    //for(int i=0;i<ncommand;i++)
    // lugive0_(commands[i].c_str(),strlen(commands[i].c_str()));
  }
  
}

---

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

• EvtPhokhara_Lambda.hh
• EvtPhokhara_Lambda.cc