EvtJ2BB3 Class Reference

#include <EvtJ2BB3.hh>

Inheritance diagram for EvtJ2BB3:

List of all members.

Public Member Functions
void	checkNArg (int a1, int a2=-1, int a3=-1, int a4=-1)
void	checkNDaug (int d1, int d2=-1)
void	checkQ ()
void	checkSpinDaughter (int d1, EvtSpinType::spintype sp)
void	checkSpinParent (EvtSpinType::spintype sp)
EvtDecayBase *	clone ()
virtual void	command (std::string cmd)
virtual std::string	commandName ()
void	decay (EvtParticle *p)
void	disableCheckQ ()
	EvtJ2BB3 ()
double	getArg (int j)
double *	getArgs ()
std::string *	getArgsStr ()
std::string	getArgStr (int j)
double	getBranchingFraction ()
EvtId	getDaug (int i)
EvtId *	getDaugs ()
int	getDSum ()
std::string	getModelName ()
void	getName (std::string &name)
int	getNArg ()
int	getNDaug ()
EvtId	getParentId ()
int	getPHOTOS ()
double	getProbMax (double prob)
void	init (int, double *)
void	init ()
virtual void	initProbMax ()
void	makeDecay (EvtParticle *p)
virtual bool	matchingDecay (const EvtDecayBase &other) const
void	noProbMax ()
virtual int	nRealDaughters ()
void	printSummary ()
double	resetProbMax (double prob)
void	saveDecayInfo (EvtId ipar, int ndaug, EvtId *daug, int narg, std::vector< std::string > &args, std::string name, double brfr)
void	setPHOTOS ()
void	setProbMax (double prbmx)
void	setSummary ()
void	setVerbose ()
void	setWeight (double weight)
int	summary ()
int	verbose ()
void	vertex (int *i1, const EvtComplex &amp)
void	vertex (int i1, int i2, int i3, const EvtComplex &amp)
void	vertex (int i1, int i2, const EvtComplex &amp)
void	vertex (int i1, const EvtComplex &amp)
void	vertex (const EvtComplex &amp)
virtual	~EvtJ2BB3 ()
Static Public Member Functions
void	findMass (EvtParticle *p)
void	findMasses (EvtParticle *p, int ndaugs, EvtId daugs[10], double masses[10])
double	findMaxMass (EvtParticle *p)
Public Attributes
double	alpha
double	u
Protected Member Functions
bool	daugsDecayedByParentModel ()
Protected Attributes
EvtAmp	_amp2
bool	_daugsDecayedByParentModel

---

Constructor & Destructor Documentation

EvtJ2BB3::EvtJ2BB3 (   )  [inline]

{}

EvtJ2BB3::~EvtJ2BB3 (   )  [virtual]

{}

---

Member Function Documentation

void EvtDecayBase::checkNArg (  int  a1, int  a2 = -1, int  a3 = -1, int  a4 = -1 )  [inherited]

{ if ( _narg != a1 && _narg != a2 && _narg != a3 && _narg != a4 ) { report(ERROR,"EvtGen") << _modelname.c_str() << " generator expected "<<endl; report(ERROR,"EvtGen") << a1<<endl;; if ( a2>-1) { report(ERROR,"EvtGen") << " or " << a2<<endl; } if ( a3>-1) { report(ERROR,"EvtGen") << " or " << a3<<endl; } if ( a4>-1) { report(ERROR,"EvtGen") << " or " << a4<<endl; } report(ERROR,"EvtGen") << " arguments but found:"<< _narg << endl; printSummary(); report(ERROR,"EvtGen") << "Will terminate execution!"<<endl; ::abort(); } }

void EvtDecayBase::checkNDaug (  int  d1, int  d2 = -1 )  [inherited]

{ if ( _ndaug != d1 && _ndaug != d2 ) { report(ERROR,"EvtGen") << _modelname.c_str() << " generator expected "; report(ERROR,"EvtGen") << d1; if ( d2>-1) { report(ERROR,"EvtGen") << " or " << d2; } report(ERROR,"EvtGen") << " daughters but found:"<< _ndaug << endl; printSummary(); report(ERROR,"EvtGen") << "Will terminate execution!"<<endl; ::abort(); } }

void EvtDecayBase::checkQ (   )  [inherited]

{ int i; int q=0; int qpar; //If there are no daughters (jetset etc) then we do not //want to do this test. Why? Because sometimes the parent //will have a nonzero charge. if ( _ndaug != 0) { for(i=0; i<_ndaug; i++ ) { q += EvtPDL::chg3(_daug[i]); } qpar = EvtPDL::chg3(_parent); if ( q != qpar ) { report(ERROR,"EvtGen") <<_modelname.c_str()<< " generator expected " << " charge to be conserved, found:"<<endl; report(ERROR,"EvtGen") << "Parent charge of "<<(qpar/3)<<endl; report(ERROR,"EvtGen") << "Sum of daughter charge of "<<(q/3)<<endl; report(ERROR,"EvtGen") << "The parent is "<< EvtPDL::name(_parent).c_str()<<endl; for(i=0; i<_ndaug; i++ ) { report(ERROR,"EvtGen") << "Daughter "<< EvtPDL::name(_daug[i]).c_str()<<endl; } report(ERROR,"EvtGen") << "Will terminate execution!"<<endl; ::abort(); } } }

void EvtDecayBase::checkSpinDaughter (  int  d1, EvtSpinType::spintype  sp )  [inherited]

{ EvtSpinType::spintype parenttype = EvtPDL::getSpinType(getDaug(d1)); if ( parenttype != sp ) { report(ERROR,"EvtGen") << _modelname.c_str() << " did not get the correct daughter spin d=" << d1 << endl; printSummary(); report(ERROR,"EvtGen") << "Will terminate execution!"<<endl; ::abort(); } }

void EvtDecayBase::checkSpinParent (  EvtSpinType::spintype  sp  )  [inherited]

{ EvtSpinType::spintype parenttype = EvtPDL::getSpinType(getParentId()); if ( parenttype != sp ) { report(ERROR,"EvtGen") << _modelname.c_str() << " did not get the correct parent spin\n"; printSummary(); report(ERROR,"EvtGen") << "Will terminate execution!"<<endl; ::abort(); } }

EvtDecayBase * EvtJ2BB3::clone (   )  [virtual]

Implements EvtDecayBase. { return new EvtJ2BB3; }

void EvtDecayBase::command (  std::string  cmd  )  [virtual, inherited]

Reimplemented in EvtJetSet, EvtLunda, EvtLundCharm, EvtPythia, and EvtTauola. { report(ERROR,"EvtGen") << "Should never call EvtDecayBase::command"<<endl; ::abort(); }

std::string EvtDecayBase::commandName (   )  [virtual, inherited]

Reimplemented in EvtJetSet, EvtLunda, EvtLundCharm, EvtPythia, and EvtTauola. { return std::string(""); }

bool EvtDecayBase::daugsDecayedByParentModel (   )  [inline, protected, inherited]

{return _daugsDecayedByParentModel;}

void EvtJ2BB3::decay (  EvtParticle *  p  )  [virtual]

Implements EvtDecayBase. { loop: p->initializePhaseSpace(getNDaug(),getDaugs()); EvtParticle *v,*s1; EvtVector4R pv,ps,ppr; v =p->getDaug(0); s1=p->getDaug(1); pv=v->getP4(); ps=s1->getP4(); ppr=p->getP4(); // Put phase space results into the daughters. EvtHelSys angles(ppr,pv); //using helicity sys. angles double theta =angles.getHelAng(1); double phi =angles.getHelAng(2); double gamma=0; double m_b0 = EvtPDL::getMass(getDaug(0)); //the mass of daughter1 is large than daughter0. double m_b1 = EvtPDL::getMass(getDaug(1)); double m_M = EvtPDL::getMass(getParentId()); double M2=pow(m_M,2.0);double b2_p=pow((m_b0+m_b1),2.0); double b2_m=pow((m_b1-m_b0),2.0); if(M2-b2_p<=0) goto loop; double P_c=sqrt((M2-b2_p)*(M2-b2_m)/(4.0*M2)); double Q=M2-b2_p; double F1=-0.54/(m_b0*b2_p); double F2=-0.54*(m_b1-m_b0)/(m_b0*(m_b0+m_b1)); double F3=-0.89/m_b0; double H1,H0,Hm1; //OK. We have all tested these parameters. if(getNArg()>0) { alpha=getArg(0); H0=1.0;H1=sqrt((1+alpha)/(1-alpha));Hm1=H1; } else{ H1=-1.15*P_c*m_M*(-M2*Q*F1-Q*F2+(M2-m_b0*(m_b0+m_b1))*F3)/(sqrt(Q)*m_b1); H0=-0.82*P_c*M2*(-(m_b1-m_b0)*Q*F1/m_b1-Q*F2/(m_b1*(m_b1-m_b0))+2.0*F3)/sqrt(Q); Hm1=2*P_c*m_M*(m_b0+m_b1)*F3/sqrt(Q); alpha=(H1*H1+Hm1*Hm1-2*H0*H0)/(H1*H1+Hm1*Hm1+2*H0*H0); } //cout<<"------alpha="<<alpha<<";H0="<<H0<<";H1="<<H1<<";Hm1="<<Hm1<<endl; // J/psi helicity =1 corresponding index=0 vertex(0,0,0,Djmn(1, 1, 0,phi,theta,gamma)*H0); vertex(0,0,1,Djmn(1, 1,-1,phi,theta,gamma)*Hm1); vertex(0,0,2,0.0); vertex(0,0,3,Djmn(1, 1, 1,phi,theta,gamma)*H1); vertex(0,1,0,Djmn(1, 1,-1,phi,theta,gamma)*H1); vertex(0,1,1,0.0); vertex(0,1,2,Djmn(1, 1, 1,phi,theta,gamma)*Hm1); vertex(0,1,3,Djmn(1, 1, 0,phi,theta,gamma)*H0); // J/psi helicity =-1 corresponding index=1 vertex(1,0,0,Djmn(1,-1, 0,phi,theta,gamma)*H0); vertex(1,0,1,Djmn(1,-1,-1,phi,theta,gamma)*Hm1); vertex(1,0,2,0.0); vertex(1,0,3,Djmn(1,-1, 1,phi,theta,gamma)*H1); vertex(1,1,0,Djmn(1,-1,-1,phi,theta,gamma)*H1); vertex(1,1,1,0.0); vertex(1,1,2,Djmn(1,-1, 1,phi,theta,gamma)*Hm1); vertex(1,1,3,Djmn(1,-1, 0,phi,theta,gamma)*H0); // J/psi helicity =0 corresponding index=2 vertex(2,0,0,Djmn(1, 0, 0,phi,theta,gamma)*H0); vertex(2,0,1,Djmn(1, 0,-1,phi,theta,gamma)*Hm1); vertex(2,0,2,0.0); vertex(2,0,3,Djmn(1, 0, 1,phi,theta,gamma)*H1); vertex(2,1,0,Djmn(1, 0,-1,phi,theta,gamma)*H1); vertex(2,1,1,0.0); vertex(2,1,2,Djmn(1, 0, 1,phi,theta,gamma)*Hm1); vertex(2,1,3,Djmn(1, 0, 0,phi,theta,gamma)*H0); return ; }

void EvtDecayBase::disableCheckQ (   )  [inline, inherited]

{_chkCharge=0;};

void EvtDecayBase::findMass (  EvtParticle *  p  )  [static, inherited]

{ //Need to also check that this mass does not screw //up the parent //This code assumes that for the ith daughter, 0..i-1 //already have a mass double maxOkMass=findMaxMass(p); int count=0; double mass; bool massOk=false; int i; while (!massOk) { count++; if ( count > 10000 ) { report(INFO,"EvtGen") << "Can not find a valid mass for: " << EvtPDL::name(p->getId()).c_str() <<endl; report(INFO,"EvtGen") << "Now printing parent and/or grandparent tree\n"; if ( p->getParent() ) { if ( p->getParent()->getParent() ) { p->getParent()->getParent()->printTree(); report(INFO,"EvtGen") << p->getParent()->getParent()->mass() <<endl; report(INFO,"EvtGen") << p->getParent()->mass() <<endl; } else{ p->getParent()->printTree(); report(INFO,"EvtGen") << p->getParent()->mass() <<endl; } } else p->printTree(); report(INFO,"EvtGen") << "maxokmass=" << maxOkMass << " " << EvtPDL::getMinMass(p->getId()) << " " << EvtPDL::getMaxMass(p->getId())<<endl; if ( p->getNDaug() ) { for (i=0; i<p->getNDaug(); i++) { report(INFO,"EvtGen") << p->getDaug(i)->mass()<<" "; } report(INFO,"EvtGen") << endl; } if ( maxOkMass >= EvtPDL::getMinMass(p->getId()) ) { report(INFO,"EvtGen") << "taking a default value\n"; p->setMass(maxOkMass); return; } assert(0); } mass = EvtPDL::getMass(p->getId()); //Just need to check that this mass is > than //the mass of all daughters double massSum=0.; if ( p->getNDaug() ) { for (i=0; i<p->getNDaug(); i++) { massSum+= p->getDaug(i)->mass(); } } //some special cases are handled with 0 (stable) or 1 (k0->ks/kl) daughters if (p->getNDaug()<2) massOk=true; if ( p->getParent() ) { if ( p->getParent()->getNDaug()==1 ) massOk=true; } if ( !massOk ) { if (massSum < mass) massOk=true; if ( mass> maxOkMass) massOk=false; } } p->setMass(mass); }

void EvtDecayBase::findMasses (  EvtParticle *  p, int  ndaugs, EvtId  daugs[10], double  masses[10] )  [static, inherited]

{ int i; double mass_sum; int count=0; if (!( p->firstornot() )) { for (i = 0; i < ndaugs; i++ ) { masses[i] = p->getDaug(i)->mass(); } //for } //if else { p->setFirstOrNot(); // if only one daughter do it if (ndaugs==1) { masses[0]=p->mass(); return; } //until we get a combo whose masses are less than _parent mass. do { mass_sum = 0.0; for (i = 0; i < ndaugs; i++ ) { masses[i] = EvtPDL::getMass(daugs[i]); mass_sum = mass_sum + masses[i]; } count++; if(count==10000) { report(ERROR,"EvtGen") <<"Decaying particle:"<< EvtPDL::name(p->getId()).c_str()<<" (m="<<p->mass()<<")"<<endl; report(ERROR,"EvtGen") <<"To the following daugthers"<<endl; for (i = 0; i < ndaugs; i++ ) { report(ERROR,"EvtGen") << EvtPDL::name(daugs[i]).c_str() << endl; } report(ERROR,"EvtGen") << "Has been rejected "<<count << " times, will now take minimal masses " << " of daugthers"<<endl; mass_sum=0.; for (i = 0; i < ndaugs; i++ ) { masses[i] = EvtPDL::getMinMass(daugs[i]); mass_sum = mass_sum + masses[i]; } if (mass_sum > p->mass()){ report(ERROR,"EvtGen") << "Parent mass="<<p->mass() << "to light for daugthers."<<endl << "Will throw the event away."<<endl; //dont terminate - start over on the event. EvtStatus::setRejectFlag(); mass_sum=0.; // ::abort(); } } } while ( mass_sum > p->mass()); } //else return; }

double EvtDecayBase::findMaxMass (  EvtParticle *  p  )  [static, inherited]

{ double maxOkMass=EvtPDL::getMaxMass(p->getId()); //protect against vphotons if ( maxOkMass < 0.0000000001 ) return 10000000.; //and against already determined masses if ( p->hasValidP4() ) maxOkMass=p->mass(); EvtParticle *par=p->getParent(); if ( par ) { double maxParMass=findMaxMass(par); int i; double minDaugMass=0.; for(i=0;i<par->getNDaug();i++){ EvtParticle *dau=par->getDaug(i); if ( dau!=p) { // it might already have a mass if ( dau->isInitialized() || dau->hasValidP4() ) minDaugMass+=dau->mass(); else //give it a bit of phase space minDaugMass+=1.000001*EvtPDL::getMinMass(dau->getId()); } } if ( maxOkMass>(maxParMass-minDaugMass)) maxOkMass=maxParMass-minDaugMass; } return maxOkMass; }

double EvtDecayBase::getArg (  int  j  )  [inherited]

{ // Verify string const char* str = _args[j].c_str(); int i = 0; while(str[i]!=0){ if (isalpha(str[i]) && str[i]!='e') { report(INFO,"EvtGen") << "String " << str << " is not a number" << endl; assert(0); } i++; } char** tc=0; return strtod(_args[j].c_str(),tc); }

double * EvtDecayBase::getArgs (   )  [inherited]

{ if ( _argsD ) return _argsD; //The user has asked for a list of doubles - the arguments //better all be doubles... if ( _narg==0 ) return _argsD; _argsD = new double[_narg]; int i; char * tc; for(i=0;i<_narg;i++) { _argsD[i] = strtod(_args[i].c_str(),&tc); } return _argsD; }

std::string* EvtDecayBase::getArgsStr (   )  [inline, inherited]

{return _args;}

std::string EvtDecayBase::getArgStr (  int  j  )  [inline, inherited]

{return _args[j];}

double EvtDecayBase::getBranchingFraction (   )  [inline, inherited]

{return _brfr;}

EvtId EvtDecayBase::getDaug (  int  i  )  [inline, inherited]

{return _daug[i];}

EvtId* EvtDecayBase::getDaugs (   )  [inline, inherited]

{return _daug;}

int EvtDecayBase::getDSum (   )  [inline, inherited]

{return _dsum; }

std::string EvtDecayBase::getModelName (   )  [inline, inherited]

{return _modelname; }

void EvtJ2BB3::getName (  std::string &  name  )  [virtual]

Implements EvtDecayBase. { model_name="J2BB3"; }

int EvtDecayBase::getNArg (   )  [inline, inherited]

{return _narg;}

int EvtDecayBase::getNDaug (   )  [inline, inherited]

{return _ndaug;}

EvtId EvtDecayBase::getParentId (   )  [inline, inherited]

{return _parent;}

int EvtDecayBase::getPHOTOS (   )  [inline, inherited]

{return _photos;}

double EvtDecayBase::getProbMax (  double  prob  )  [inherited]

{ int i; //diagnostics sum_prob+=prob; if (prob>max_prob) max_prob=prob; if ( defaultprobmax && ntimes_prob<=500 ) { //We are building up probmax with this iteration ntimes_prob += 1; if ( prob > probmax ) { probmax = prob;} if (ntimes_prob==500) { probmax*=1.2; } return 1000000.0*prob; } if ( prob> probmax*1.0001) { report(INFO,"EvtGen") << "prob > probmax:("<<prob<<">"<<probmax<<")"; report(INFO,"") << "("<<_modelname.c_str()<<") "; report(INFO,"") << EvtPDL::name(_parent).c_str()<<" -> "; for(i=0;i<_ndaug;i++){ report(INFO,"") << EvtPDL::name(_daug[i]).c_str() << " "; } report(INFO,"") << endl; if (defaultprobmax) probmax = prob; } ntimes_prob += 1; return probmax; } //getProbMax

void EvtJ2BB3::init (  int  , double *  )

void EvtJ2BB3::init (   )  [virtual]

Reimplemented from EvtDecayBase. { // checkNArg(1); checkNDaug(2); checkSpinParent(EvtSpinType::VECTOR); checkSpinDaughter(0,EvtSpinType::DIRAC); checkSpinDaughter(1,EvtSpinType::RARITASCHWINGER); }

void EvtDecayBase::initProbMax (   )  [virtual, inherited]

Reimplemented in EvtAngH2, EvtAngSam, EvtAngSam3, EvtAV2GV, EvtBody3, EvtBsquark, EvtBto2piCPiso, EvtBTo3piCP, EvtBtoKD3P, EvtBtoKpiCPiso, EvtbTosllAli, EvtbTosllBall, EvtBtoXsEtap, EvtBtoXsgamma, EvtBtoXsll, EvtCBTo3piMPP, EvtCBTo3piP00, EvtD0mixDalitz, EvtDDalitz, EvtDIY, EvtDMix, EvtEtaDalitz, EvtFlatQ2, EvtGoityRoberts, EvtHAngSam3, EvtHelAmp, EvtHQET, EvtHQET2, EvtHypNonLepton, EvtIntervalDecayAmp< T >, EvtISGW2, EvtJetSet, EvtJpipi, EvtJscont, EvtKKLambdaC, EvtKstarstargamma, EvtLambdaP_BarGamma, EvtLNuGamma, EvtLunda, EvtLundCharm, EvtMassH1, EvtMassH2, EvtMBody3, EvtmPhsp, EvtMultibody, EvtOmegaDalitz, EvtP2GC0, EvtP2GC1, EvtP2GC2, EvtPartWave, EvtPBB1, EvtPBB2, EvtPhsp, EvtPi0Dalitz, EvtPycont, EvtPyGaGa, EvtPythia, EvtRhoPi, EvtS2GV, EvtSLBKPole, EvtSLN, EvtSLPole, EvtSPL, EvtSSDCP, EvtSSSCP, EvtSSSCPpng, EvtSTS, EvtSTSCP, EvtSVPCP, EvtSVPHelAmp, EvtSVS, EvtSVSCP, EvtSVSCPiso, EvtSVSCPLH, EvtSVSNONCPEIGEN, EvtSVVCP, EvtSVVCPLH, EvtSVVHelAmp, EvtSVVNONCPEIGEN, EvtT2GV, EvtTauHadnu, EvtTaulnunu, EvtTauola, EvtTauScalarnu, EvtTauVectornu, EvtTSS, EvtTVSPwave, EvtUclaII, EvtVectorIsr, EvtVll, EvtVPHOtoVISR, EvtVPHOtoVISRHi, EvtVSPPwave, EvtVSS, EvtVSSBMixCPT, EvtVSSMix, EvtVub, EvtVubHybrid, EvtVubNLO, EvtVVP, EvtVVpipi, EvtVVPIPI_WEIGHTED, EvtVVSPwave, and EvtIntervalDecayAmp< EvtDalitzPoint >. { //This function is called if the decay does not have a //specialized initialization. //The default is to set the maximum //probability to 0 and the number of times called to 0 //and defaultprobmax to 1 such that the decay will be //generated many many times //in order to generate a reasonable maximum probability //for the decay. defaultprobmax=1; ntimes_prob = 0; probmax = 0.0; } //initProbMax

void EvtDecayAmp::makeDecay (  EvtParticle *  p  )  [virtual, inherited]

Implements EvtDecayBase. { int ntimes=10000; int i,more; EvtSpinDensity rho; double prob,prob_max; _amp2.init(p->getId(),getNDaug(),getDaugs()); // report(INFO,"EvtGen") << "Decaying " << EvtPDL::name(p->getId()) << endl; do{ _daugsDecayedByParentModel=false; _weight = 1.0; decay(p); rho=_amp2.getSpinDensity(); prob=p->getSpinDensityForward().NormalizedProb(rho); if (prob<0.0) { report(ERROR,"EvtGen")<<"Negative prob:"<<p->getId().getId() <<" "<<p->getChannel()<<endl; report(ERROR,"EvtGen") << "rho_forward:"<<endl; report(ERROR,"EvtGen") << p->getSpinDensityForward(); report(ERROR,"EvtGen") << "rho decay:"<<endl; report(ERROR,"EvtGen") << rho <<endl; } if (prob!=prob) { report(DEBUG,"EvtGen") << "Forward density matrix:"<<endl; report(ERROR,"EvtGen") << p->getSpinDensityForward(); report(DEBUG,"EvtGen") << "Decay density matrix:"<<endl; report(ERROR,"EvtGen") << rho; report(DEBUG,"EvtGen") << "prob:"<<prob<<endl; report(DEBUG,"EvtGen") << "Particle:" <<EvtPDL::name(p->getId()).c_str()<<endl; report(DEBUG,"EvtGen") << "channel :"<<p->getChannel()<<endl; report(DEBUG,"EvtGen") << "Momentum:" << p->getP4() << " " << p->mass() << endl; if( p->getParent()!=0){ report(DEBUG,"EvtGen") << "parent:" <<EvtPDL::name( p->getParent()->getId()).c_str()<<endl; report(DEBUG,"EvtGen") << "parent channel :" <<p->getParent()->getChannel()<<endl; int i; report(DEBUG,"EvtGen") << "parent daughters :"; for (i=0;i<p->getParent()->getNDaug();i++){ report(DEBUG,"") << EvtPDL::name( p->getParent()->getDaug(i)->getId()).c_str() << " "; } report(DEBUG,"") << endl; report(DEBUG,"EvtGen") << "daughters :"; for (i=0;i<p->getNDaug();i++){ report(DEBUG,"") << EvtPDL::name( p->getDaug(i)->getId()).c_str() << " "; } report(DEBUG,"") << endl; report(DEBUG,"EvtGen") << "daughter momenta :" << endl;; for (i=0;i<p->getNDaug();i++){ report(DEBUG,"") << p->getDaug(i)->getP4() << " " << p->getDaug(i)->mass(); report(DEBUG,"") << endl; } } } prob/=_weight; prob_max = getProbMax(prob); p->setDecayProb(prob/prob_max); //report(INFO,"EvtGen") << "Prob,prob_max,weight:"<<prob<<" "<<prob_max<<" "<<_weight<<endl; more=prob<EvtRandom::Flat(prob_max); ntimes--; }while(ntimes&&more); //report(INFO,"EvtGen") << "Done\n"; if (ntimes==0){ report(DEBUG,"EvtGen") << "Tried accept/reject:10000" <<" times, and rejected all the times!"<<endl; report(DEBUG,"")<<p->getSpinDensityForward()<<endl; report(DEBUG,"EvtGen") << "Is therefore accepting the last event!"<<endl; report(DEBUG,"EvtGen") << "Decay of particle:"<< EvtPDL::name(p->getId()).c_str()<<"(channel:"<< p->getChannel()<<") with mass "<<p->mass()<<endl; int ii; for(ii=0;ii<p->getNDaug();ii++){ report(DEBUG,"EvtGen") <<"Daughter "<<ii<<":"<< EvtPDL::name(p->getDaug(ii)->getId()).c_str()<<" with mass "<< p->getDaug(ii)->mass()<<endl; } } EvtSpinDensity rho_list[10]; rho_list[0]=p->getSpinDensityForward(); EvtAmp ampcont; if (_amp2._pstates!=1){ ampcont=_amp2.contract(0,p->getSpinDensityForward()); // J2BB2 bugging here } else{ ampcont=_amp2; } // it may be that the parent decay model has already // done the decay - this should be rare and the // model better know what it is doing.. if ( !daugsDecayedByParentModel() ){ // report(INFO,"EvtGen") << "Found " << p->getNDaug() << " daughters\n"; for(i=0;i<p->getNDaug();i++){ rho.SetDim(_amp2.dstates[i]); if (_amp2.dstates[i]==1) { rho.Set(0,0,EvtComplex(1.0,0.0)); } else{ rho=ampcont.contract(_amp2._dnontrivial[i],_amp2); } if (!rho.Check()) { report(ERROR,"EvtGen") << "-------start error-------"<<endl; report(ERROR,"EvtGen")<<"forward rho failed Check:"<< EvtPDL::name(p->getId()).c_str()<<" "<<p->getChannel()<<" "<<i<<endl; report(ERROR,"EvtGen")<<"Parent:"<<EvtPDL::name(p->getParent()->getId()).c_str()<<endl; report(ERROR,"EvtGen")<<"GrandParent:"<<EvtPDL::name(p->getParent()->getParent()->getId()).c_str()<<endl; report(ERROR,"EvtGen")<<"GrandGrandParent:"<<EvtPDL::name(p->getParent()->getParent()->getParent()->getId()).c_str()<<endl; report(ERROR,"EvtGen") << rho; int ii; _amp2.dump(); for(ii=0;ii<i+1;ii++){ report(ERROR,"EvtGen") << rho_list[ii]; } report(ERROR,"EvtGen") << "-------Done with error-------"<<endl; } p->getDaug(i)->setSpinDensityForward(rho); p->getDaug(i)->decay(); rho_list[i+1]= p->getDaug(i)->getSpinDensityBackward(); if (_amp2.dstates[i]!=1){ ampcont=ampcont.contract(_amp2._dnontrivial[i],rho_list[i+1]); } } p->setSpinDensityBackward(_amp2.getBackwardSpinDensity(rho_list)); if (!p->getSpinDensityBackward().Check()) { report(ERROR,"EvtGen")<<"rho_backward failed Check"<< p->getId().getId()<<" "<<p->getChannel()<<endl; report(ERROR,"EvtGen") << p->getSpinDensityBackward(); } } if (getPHOTOS() || EvtRadCorr::alwaysRadCorr()) { EvtRadCorr::doRadCorr(p); } }

bool EvtDecayBase::matchingDecay (  const EvtDecayBase &  other  )  const [virtual, inherited]

{ if ( _ndaug != other._ndaug) return false; if ( _parent != other._parent) return false; std::vector<int> useDs; for ( unsigned int i=0; i<_ndaug; i++) useDs.push_back(0); for ( unsigned int i=0; i<_ndaug; i++) { bool foundIt=false; for ( unsigned int j=0; j<_ndaug; j++) { if ( useDs[j] == 1 ) continue; if ( _daug[i] == other._daug[j] && _daug[i].getAlias() == other._daug[j].getAlias()) { foundIt=true; useDs[j]=1; break; } } if ( foundIt==false) return false; } for ( unsigned int i=0; i<_ndaug; i++) if ( useDs[i]==0) return false; return true; }

void EvtDecayBase::noProbMax (   )  [inherited]

{ defaultprobmax=0; }

virtual int EvtDecayBase::nRealDaughters (   )  [inline, virtual, inherited]

Reimplemented in EvtBtoKD3P, and EvtVSSBMixCPT. { return _ndaug;}

void EvtDecayBase::printSummary (   )  [inherited]

{ int i; if (ntimes_prob>0) { report(INFO,"EvtGen") << "Calls="<<ntimes_prob<<" eff:"<< sum_prob/(probmax*ntimes_prob)<<" frac. max:"<<max_prob/probmax; report(INFO,"") <<" probmax:"<<probmax<<" max:"<<max_prob<<" : "; } report(INFO,"") << EvtPDL::name(_parent).c_str()<<" -> "; for(i=0;i<_ndaug;i++){ report(INFO,"") << EvtPDL::name(_daug[i]).c_str() << " "; } report(INFO,"") << " ("<<_modelname.c_str()<<"):"<< endl; }

double EvtDecayBase::resetProbMax (  double  prob  )  [inherited]

{ report(INFO,"EvtGen") << "Reseting prob max\n"; report(INFO,"EvtGen") << "prob > probmax:("<<prob<<">"<<probmax<<")"; report(INFO,"") << "("<<_modelname.c_str()<<")"; report(INFO,"") << EvtPDL::getStdHep(_parent)<<"->"; for( int i=0;i<_ndaug;i++){ report(INFO,"") << EvtPDL::getStdHep(_daug[i]) << " "; } report(INFO,"") << endl; probmax = 0.0; defaultprobmax = 0; ntimes_prob = 0; return prob; }

void EvtDecayBase::saveDecayInfo (  EvtId  ipar, int  ndaug, EvtId *  daug, int  narg, std::vector< std::string > &  args, std::string  name, double  brfr )  [inherited]

{ int i; _brfr=brfr; _ndaug=ndaug; _narg=narg; _parent=ipar; _dsum=0; if (_ndaug>0) { _daug=new EvtId [_ndaug]; for(i=0;i<_ndaug;i++){ _daug[i]=daug[i]; _dsum+=daug[i].getAlias(); } } else{ _daug=0; } if (_narg>0) { _args=new std::string[_narg+1]; for(i=0;i<_narg;i++){ _args[i]=args[i]; } } else{ _args = 0; } _modelname=name; this->init(); this->initProbMax(); if (_chkCharge){ this->checkQ(); } if (defaultprobmax){ report(INFO,"EvtGen") << "No default probmax for "; report(INFO,"") << "("<<_modelname.c_str()<<") "; report(INFO,"") << EvtPDL::name(_parent).c_str()<<" -> "; for(i=0;i<_ndaug;i++){ report(INFO,"") << EvtPDL::name(_daug[i]).c_str() << " "; } report(INFO,"") << endl; report(INFO,"") << "This is fine for development, but must be provided for production."<<endl; report(INFO,"EvtGen") << "Never fear though - the decay will use the \n"; report(INFO,"EvtGen") << "500 iterations to build up a good probmax \n"; report(INFO,"EvtGen") << "before accepting a decay. "<<endl; } }

void EvtDecayBase::setPHOTOS (   )  [inline, inherited]

{_photos=1;}

void EvtDecayBase::setProbMax (  double  prbmx  )  [inherited]

{ defaultprobmax=0; probmax=prbmx; }

void EvtDecayBase::setSummary (   )  [inline, inherited]

{_summary=1;}

void EvtDecayBase::setVerbose (   )  [inline, inherited]

{_verbose=1;}

void EvtDecayAmp::setWeight (  double  weight  )  [inline, inherited]

{_weight=weight;}

int EvtDecayBase::summary (   )  [inline, inherited]

{return _summary; }

int EvtDecayBase::verbose (   )  [inline, inherited]

{return _verbose; }

void EvtDecayAmp::vertex (  int *  i1, const EvtComplex &  amp )  [inline, inherited]

sets the amplitudes calculated in the decay objects { _amp2.vertex(i1,amp);}

void EvtDecayAmp::vertex (  int  i1, int  i2, int  i3, const EvtComplex &  amp )  [inline, inherited]

sets the amplitudes calculated in the decay objects {_amp2.vertex(i1,i2,i3,amp);}

void EvtDecayAmp::vertex (  int  i1, int  i2, const EvtComplex &  amp )  [inline, inherited]

sets the amplitudes calculated in the decay objects {_amp2.vertex(i1,i2,amp);}

void EvtDecayAmp::vertex (  int  i1, const EvtComplex &  amp )  [inline, inherited]

sets the amplitudes calculated in the decay objects {_amp2.vertex(i1,amp);}

void EvtDecayAmp::vertex (  const EvtComplex &  amp  )  [inline, inherited]

sets the amplitudes calculated in the decay objects {_amp2.vertex(amp);}

---

Member Data Documentation

EvtAmp EvtDecayAmp::_amp2 [protected, inherited]

bool EvtDecayBase::_daugsDecayedByParentModel [protected, inherited]

double EvtJ2BB3::alpha

double EvtJ2BB3::u

---

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/EvtGenModels/EvtJ2BB3.hh
• /data0/mdestefa/bes3soft/Boss_6.5.5/dist/6.5.5/Generator/BesEvtGen/BesEvtGen-00-01-11/src/EvtGen/EvtGenModels/EvtJ2BB3.cc

---