EvtHypNonLepton Class Reference

#include <EvtHypNonLepton.hh>

Inheritance diagram for EvtHypNonLepton:

Public Member Functions
	EvtHypNonLepton ()
virtual	~EvtHypNonLepton ()
void	getName (std::string &name)
EvtDecayBase *	clone ()
void	decay (EvtParticle *p)
void	init ()
void	initProbMax ()
void	calcAmp (EvtAmp *amp, EvtParticle *parent)
Public Member Functions inherited from EvtDecayAmp
void	makeDecay (EvtParticle *p)
void	setWeight (double weight)
void	vertex (const EvtComplex &amp)
void	vertex (int i1, const EvtComplex &amp)
void	vertex (int i1, int i2, const EvtComplex &amp)
void	vertex (int i1, int i2, int i3, const EvtComplex &amp)
void	vertex (int *i1, const EvtComplex &amp)
virtual	~EvtDecayAmp ()
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 EvtDecayAmp
EvtAmp	_amp2
Protected Attributes inherited from EvtDecayBase
bool	_daugsDecayedByParentModel

Detailed Description

Definition at line 20 of file EvtHypNonLepton.hh.

Constructor & Destructor Documentation

EvtHypNonLepton()

EvtHypNonLepton::EvtHypNonLepton ( )

inline

Definition at line 24 of file EvtHypNonLepton.hh.

{}

Referenced by clone().

~EvtHypNonLepton()

EvtHypNonLepton::~EvtHypNonLepton ( )

virtual

Definition at line 31 of file EvtHypNonLepton.cc.

{}

Member Function Documentation

calcAmp()

void EvtHypNonLepton::calcAmp	(	EvtAmp *	amp,
		EvtParticle *	parent
	)

Definition at line 129 of file EvtHypNonLepton.cc.

                                                            {
 
  static long noTries=0;
  int i;
  EvtComplex Matrix[2][2],B_to_A;
 
  //G_F  = 1.16637e-5;
  //M_pi = 0.13957;
 
  for(i=0;i<4;i++){
    //std::cout << "--------------------------------------------------" << std::endl;
    Matrix[i/2][i%2]  = EvtLeptonSCurrent(parent->sp(i/2),parent->getDaug(0)->spParent(i%2));
    //std::cout << "Matrix = " << Matrix[i/2][i%2] << std::endl;
    Matrix[i/2][i%2] -= m_B_to_A*EvtLeptonPCurrent(parent->sp(i/2),parent->getDaug(0)->spParent(i%2));
    //std::cout << "Matrix = " << Matrix[i/2][i%2] << std::endl;
    //Matrix[i/2][i%2] *= G_F*M_pi*M_pi;
    //std::cout << "Matrix = " << Matrix[i/2][i%2] << std::endl;
    //std::cout << "--------------------------------------------------" << std::endl;
    amp->vertex(i/2,i%2,Matrix[i/2][i%2]);
  }
 
  if(m_noTries>0) if(!((++noTries)%m_noTries)) report(DEBUG,"EvtGen") << " EvtHypNonLepton already finished " << noTries << " matrix element calculations" << std::endl;
}

Referenced by decay().

clone()

EvtDecayBase * EvtHypNonLepton::clone ( )

virtual

Implements EvtDecayBase.

Definition at line 33 of file EvtHypNonLepton.cc.

                                    {
  return new EvtHypNonLepton;
}

decay()

void EvtHypNonLepton::decay ( EvtParticle *  p )

virtual

Implements EvtDecayBase.

Definition at line 122 of file EvtHypNonLepton.cc.

                                              {
 
  parent->initializePhaseSpace(getNDaug(),getDaugs());
  calcAmp(&_amp2,parent);
  
}

getName()

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

virtual

Implements EvtDecayBase.

Definition at line 37 of file EvtHypNonLepton.cc.

                                                  {
  model_name= "HypNonLepton";
}

init()

void EvtHypNonLepton::init ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 41 of file EvtHypNonLepton.cc.

                          {
 
  if(getNArg()<2 || getNArg()>3){ // alpha phi gamma delta
    report(ERROR,"EvtGen") << " ERROR: EvtHypNonLepton generator expected 2 or 3 arguments but found: " << getNArg() << std::endl;
    report(INFO ,"EvtGen") << "  1. Decay asymmetry parameter - alpha" << std::endl;
    report(INFO ,"EvtGen") << "  2. Parameter phi - in degrees (not radians)" << std::endl;
    report(INFO ,"EvtGen") << "  3. Note on every x-th decay" << std::endl;
    ::abort();
  }
 
  if(getNDaug()!=2){ // Check that there are 2 daughters only
    report(ERROR,"EvtGen") << " ERROR: EvtHypNonLepton generator expected 2 daughters but found: " << getNDaug() << std::endl;
    ::abort();
  }
 
  // Check particles spins
  if(EvtSpinType::getSpin2(EvtPDL::getSpinType(getParentId()))!=1){
    report(ERROR,"EvtGen") << " ERROR: EvtHypNonLepton generator expected dirac parent particle, but found " << EvtSpinType::getSpin2(EvtPDL::getSpinType(getParentId())) << " spin degrees of freedom" << std::endl;
    ::abort();
  }
  if(EvtSpinType::getSpin2(EvtPDL::getSpinType(getDaug(0)))!=1){
    report(ERROR,"EvtGen") << " ERROR: EvtHypNonLepton generator expected the first child to be dirac particle, but found " << EvtSpinType::getSpin2(EvtPDL::getSpinType(getDaug(0))) << " spin degrees of freedom" << std::endl;
    ::abort();
  }
  if(EvtSpinType::getSpin2(EvtPDL::getSpinType(getDaug(1)))!=0){
    report(ERROR,"EvtGen") << " ERROR: EvtHypNonLepton generator expected the second child to be scalar particle, but found " << EvtSpinType::getSpin2(EvtPDL::getSpinType(getDaug(1))) << " spin degrees of freedom" << std::endl;
    ::abort();
  }
 
  // Read all parameters
  m_alpha = getArg(0);
  m_phi   = getArg(1)*EvtConst::pi/180;
  if(getNArg()==3) m_noTries = static_cast<long>(getArg(2));
  else             m_noTries = 0;
 
  // calculate additional parameters
  double p,M,m1,m2;
  double p_to_s,beta,delta,gamma;
 
  M  = EvtPDL::getMass(getParentId());
  m1 = EvtPDL::getMass(getDaug(0));
  m2 = EvtPDL::getMass(getDaug(1));
 
  if(m1+m2>=M){
    report(ERROR,"EvtGen") << " ERROR: EvtHypNonLepton found impossible decay: " << M << " --> " << m1 << " + " << m2 << " GeV\n" << std::endl;
    ::abort();
  }
 
  p = sqrt(M*M-(m1+m2)*(m1+m2))*sqrt(M*M-(m1-m2)*(m1-m2))/2./M;
 
  beta = sqrt(1.-m_alpha*m_alpha)*sin(m_phi);
  delta = -atan2(beta,m_alpha);
  gamma = sqrt(1.-m_alpha*m_alpha-beta*beta);
  p_to_s = sqrt((1.-gamma)/(1.+gamma));
 
  m_B_to_A = p_to_s*(m1+sqrt(p*p+m1*m1))/p*EvtComplex(cos(delta),sin(delta));
 
}

initProbMax()

void EvtHypNonLepton::initProbMax ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 100 of file EvtHypNonLepton.cc.

                                 {
 
  double maxProb,m1,m2,M,p;
 
  M=EvtPDL::getMass(getParentId());
  m1=EvtPDL::getMass(getDaug(0));
  m2=EvtPDL::getMass(getDaug(1));
 
  if(m1+m2>=M){
    report(ERROR,"EvtGen") << " ERROR: EvtHypNonLepton found impossible decay: " << M << " --> " << m1 << " + " << m2 << " GeV\n" << std::endl;
    ::abort();
  }
 
  p=sqrt(M*M-(m1+m2)*(m1+m2))*sqrt(M*M-(m1-m2)*(m1-m2))/2/M;
  maxProb=16*M*(sqrt(p*p+m1*m1)+m1+abs(m_B_to_A)*abs(m_B_to_A)*(sqrt(p*p+m1*m1)-m1));
  //maxProb *= G_F*M_pi*M_pi;
 
  setProbMax(maxProb);
  report(INFO,"EvtGen") << " EvtHypNonLepton set up maximum probability to " << maxProb << std::endl;
 
}

---

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

• EvtHypNonLepton.hh
• EvtHypNonLepton.cc