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 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 21 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 152 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 146 of file EvtHypNonLepton.cc.

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

getName()

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

virtual

Implements EvtDecayBase.

Definition at line 35 of file EvtHypNonLepton.cc.

{ model_name = "HypNonLepton"; }

init()

void EvtHypNonLepton::init ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 37 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 118 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