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 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 28 of file EvtBsquark.hh.

Constructor & Destructor Documentation

EvtBsquark()

EvtBsquark::EvtBsquark ( )

inline

Definition at line 32 of file EvtBsquark.hh.

{}

Referenced by clone().

~EvtBsquark()

EvtBsquark::~EvtBsquark ( )

virtual

Definition at line 36 of file EvtBsquark.cc.

{}

Member Function Documentation

clone()

EvtDecayBase * EvtBsquark::clone ( )

virtual

Implements EvtDecayBase.

Definition at line 45 of file EvtBsquark.cc.

                               {
 
  return new EvtBsquark;
 
}

decay()

void EvtBsquark::decay ( EvtParticle *  p )

virtual

Implements EvtDecayBase.

Definition at line 65 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 38 of file EvtBsquark.cc.

                                             {
 
  model_name="BSQUARK";     
 
}

init()

void EvtBsquark::init ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 51 of file EvtBsquark.cc.

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

initProbMax()

void EvtBsquark::initProbMax ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 57 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