EvtBsquark Class Reference

#include <EvtBsquark.hh>

Inheritance diagram for EvtBsquark:

Public Member Functions
	EvtBsquark ()
virtual	~EvtBsquark ()
void	getName (std::string &name)
EvtDecayBase *	clone ()
void	initProbMax ()
void	init ()
void	decay (EvtParticle *p)
Public Member Functions inherited from EvtDecayProb
void	makeDecay (EvtParticle *p)
void	setProb (double prob)
double	getProb ()
void	setWeight (double weight)
virtual	~EvtDecayProb ()
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 EvtDecayBase
bool	_daugsDecayedByParentModel

Detailed Description

Definition at line 28 of file EvtBsquark.hh.

Constructor & Destructor Documentation

EvtBsquark()

EvtBsquark::EvtBsquark ( )

inline

Definition at line 31 of file EvtBsquark.hh.

{}

Referenced by clone().

~EvtBsquark()

EvtBsquark::~EvtBsquark ( )

virtual

Definition at line 35 of file EvtBsquark.cc.

{}

Member Function Documentation

clone()

EvtDecayBase * EvtBsquark::clone ( )

virtual

Implements EvtDecayBase.

Definition at line 39 of file EvtBsquark.cc.

{ return new EvtBsquark; }

decay()

void EvtBsquark::decay ( EvtParticle * p )

virtual

Implements EvtDecayBase.

Definition at line 54 of file EvtBsquark.cc.

                                       {
 
  static EvtId cquark     = EvtPDL::getId( "c" );
  static EvtId anticquark = EvtPDL::getId( "anti-c" );
 
  static EvtIdSet leptons( "e-", "mu-", "tau-" );
 
  p->initializePhaseSpace( getNDaug(), getDaugs() );
 
  int charge = 1;
 
  EvtParticle* lepton;
  lepton = p->getDaug( 1 );
  if ( leptons.contains( lepton->getId() ) ) { charge = -1; }
 
  EvtDiracParticle charmquark;
 
  // this is a very crude approximation...
  if ( charge == -1 ) { charmquark.init( cquark, p->getDaug( 0 )->getP4() ); }
  else { charmquark.init( anticquark, p->getDaug( 0 )->getP4() ); }
 
  EvtVector4R p4c = p->getDaug( 0 )->getP4();
 
  EvtVector4R p4sn = p->getDaug( 2 )->getP4();
 
  EvtVector4R p4b( p->mass(), 0.0, 0.0, 0.0 );
 
  EvtComplex M[2][2];
 
  int il, ic;
 
  // project out the right handed current
  EvtGammaMatrix PR = 0.5 * ( EvtGammaMatrix::id() + EvtGammaMatrix::g5() );
 
  double tanbeta = getArg( 1 );
  double cosbeta = cos( atan( tanbeta ) );
  double sinbeta = sin( atan( tanbeta ) );
 
  double mb = 4.9;
  double mc = 1.3;
  double mw = 80.4;
 
  double Mass      = getArg( 2 );
  double mu        = getArg( 3 );
  double mchargino = getArg( 4 );
 
  double tan2phim = 2 * sqrt( 2.0 ) * mw * ( mu * cosbeta + Mass * sinbeta ) /
                    ( Mass * Mass - mu * mu + 2 * mw * mw * cos( 2 * atan( tanbeta ) ) );
 
  double phim = 0.5 * atan( tan2phim );
 
  EvtComplex U11 = cos( phim );
  EvtComplex U12 = sin( phim );
  EvtComplex U21 = -sin( phim );
  EvtComplex U22 = cos( phim );
 
  double tan2phip = 2 * sqrt( 2.0 ) * mw * ( mu * cosbeta + Mass * sinbeta ) /
                    ( Mass * Mass - mu * mu - 2 * mw * mw * cos( 2 * atan( tanbeta ) ) );
 
  double phip = 0.5 * atan( tan2phip );
 
  EvtComplex V11 = cos( phip );
  EvtComplex V12 = sin( phip );
  EvtComplex V21 = -sin( phip );
  EvtComplex V22 = cos( phip );
 
  double theta  = getArg( 0 );
  double ctheta = cos( theta );
  double stheta = sin( theta );
 
  double vcsb  = 0.08;
  double mchi1 = mchargino;
  double mchi2 = mchargino;
 
  // overall scale factor
  double g = 1.0;
 
  EvtComplex a1 =
      mchi1 * ( U11 * ctheta - mb * U12 * stheta / ( sqrt( 2.0 ) * mw * cosbeta ) );
  EvtComplex a2 =
      mchi2 * ( U21 * ctheta - mb * U22 * stheta / ( sqrt( 2.0 ) * mw * cosbeta ) );
 
  EvtComplex b1 = mc * conj( V12 ) * ctheta / ( sqrt( 2.0 ) * mw * sinbeta );
  EvtComplex b2 = mc * conj( V22 ) * ctheta / ( sqrt( 2.0 ) * mw * sinbeta );
 
  EvtComplex f1 = -( g * g * V11 * vcsb ) / ( ( p4b - p4c ).mass2() - mchi1 * mchi1 );
  EvtComplex f2 = -( g * g * V21 * vcsb ) / ( ( p4b - p4c ).mass2() - mchi1 * mchi2 );
 
  // report(INFO,"EvtGen") <<g<<" "<<V11<<" "<<FL<<" "<<vcsb<<" "<<mchi1<<endl;
  // report(INFO,"EvtGen") << "f1:"<<f1<<" "<<(p4b-p4c).mass2()<<endl;
  // report(INFO,"EvtGen") << "f2:"<<f2<<" "<<(p4b-p4c).mass2()<<endl;
 
  // report(INFO,"EvtGen") << "p4sn:"<<p4sn<<endl;
 
  EvtGammaMatrix pslash =
      p4sn.get( 0 ) * EvtGammaMatrix::g0() - p4sn.get( 1 ) * EvtGammaMatrix::g1() -
      p4sn.get( 2 ) * EvtGammaMatrix::g2() - p4sn.get( 3 ) * EvtGammaMatrix::g3();
 
  // report(INFO,"EvtGen") << "pslash:"<<pslash<<endl;
 
  for ( il = 0; il < 2; il++ )
  {
    for ( ic = 0; ic < 2; ic++ )
    {
 
      EvtComplex a = 0.0;
      EvtComplex b = 0.0;
 
      if ( charge == -1 )
      {
        a = charmquark.spParent( ic ) * ( PR * lepton->spParent( il ) );
        b = charmquark.spParent( ic ) * ( ( pslash * PR ) * lepton->spParent( il ) );
      }
      else
      {
        a = lepton->spParent( il ) * ( PR * charmquark.spParent( ic ) );
        b = lepton->spParent( il ) * ( ( pslash * PR ) * charmquark.spParent( ic ) );
      }
 
      // report(INFO,"EvtGen") <<"pslash*PR:"<<pslash*PR<<endl;
      // report(INFO,"EvtGen") <<"sp charm:"<<charmquark.spParent(ic)<<endl;
      // report(INFO,"EvtGen") <<"sp lepton:"<<lepton->spParent(il)<<endl;
 
      M[ic][il] = f1 * ( a1 * a + b1 * b ) + f2 * ( a2 * a + b2 * b );
 
      // report(INFO,"EvtGen") << "Contr1:"<<a1<<" "<<a<<" "<<b1<<" "<<b<<endl;
      // report(INFO,"EvtGen") << "Contr2:"<<a2<<" "<<a<<" "<<b2<<" "<<b<<endl;
 
      // report(INFO,"EvtGen") <<"case1:"<<f1<<" "<<a1<<" "<<b1<<" "<<a<<" "<<b<<endl;
      // report(INFO,"EvtGen") <<"case2:"<<f2<<" "<<a2<<" "<<b2<<" "<<a<<" "<<b<<endl;
    }
  }
 
  double prob = real( M[0][0] * conj( M[0][0] ) + M[1][0] * conj( M[1][0] ) +
                      M[0][1] * conj( M[0][1] ) + M[1][1] * conj( M[1][1] ) );
 
  // report(INFO,"EvtGen") <<"prob:"<<prob<<endl;
 
  setProb( prob );
 
  return;
}

getName()

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

virtual

Implements EvtDecayBase.

Definition at line 37 of file EvtBsquark.cc.

{ model_name = "BSQUARK"; }

init()

void EvtBsquark::init ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 41 of file EvtBsquark.cc.

                      {
 
  // check that there are 5 arguments
  checkNArg( 5 );
}

initProbMax()

void EvtBsquark::initProbMax ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 47 of file EvtBsquark.cc.

                             {
 
  // For now do not set a maximum.
 
  // SetProbMax(0.000000000005);
}

---

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

• EvtBsquark.hh
• EvtBsquark.cc