EvtBtoXsll Class Reference

#include <EvtBtoXsll.hh>

Inheritance diagram for EvtBtoXsll:

Public Member Functions
	EvtBtoXsll ()
virtual	~EvtBtoXsll ()
void	getName (std::string &name)
EvtDecayBase *	clone ()
void	initProbMax ()
void	init ()
void	decay (EvtParticle *p)
Public Member Functions inherited from EvtDecayIncoherent
void	makeDecay (EvtParticle *p)
virtual	~EvtDecayIncoherent ()
void	setDaughterSpinDensity (int daughter)
int	isDaughterSpinDensitySet (int daughter)
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 EvtBtoXsll.hh.

Constructor & Destructor Documentation

EvtBtoXsll()

EvtBtoXsll::EvtBtoXsll ( )

inline

Definition at line 31 of file EvtBtoXsll.hh.

{}

Referenced by clone().

~EvtBtoXsll()

EvtBtoXsll::~EvtBtoXsll ( )

virtual

Definition at line 37 of file EvtBtoXsll.cc.

{}

Member Function Documentation

clone()

EvtDecayBase * EvtBtoXsll::clone ( )

virtual

Implements EvtDecayBase.

Definition at line 41 of file EvtBtoXsll.cc.

{ return new EvtBtoXsll; }

decay()

void EvtBtoXsll::decay ( EvtParticle * p )

virtual

Implements EvtDecayBase.

Definition at line 184 of file EvtBtoXsll.cc.

                                       {
 
  p->makeDaughters( getNDaug(), getDaugs() );
 
  EvtParticle* xhadron = p->getDaug( 0 );
  EvtParticle* leptonp = p->getDaug( 1 );
  EvtParticle* leptonn = p->getDaug( 2 );
 
  double mass[3];
 
  findMasses( p, getNDaug(), getDaugs(), mass );
 
  double mB = p->mass();
  double ml = mass[1];
  double pb;
 
  int        im          = 0;
  static int nmsg        = 0;
  double     xhadronMass = -999.0;
 
  EvtVector4R p4xhadron;
  EvtVector4R p4leptonp;
  EvtVector4R p4leptonn;
 
  // require the hadronic system has mass greater than that of a Kaon pion pair
 
  //  while (xhadronMass < 0.6333)
  // the above minimum value of K+pi mass appears to be too close
  // to threshold as far as JETSET is concerned
  // (JETSET gets caught in an infinite loop)
  // so we choose a lightly larger value for the threshold
  while ( xhadronMass < _mxmin )
  {
    im++;
 
    // Apply Fermi motion and determine effective b-quark mass
 
    // Old BaBar MC parameters
    //    double pf = 0.25;
    //    double ms = 0.2;
    //    double mq = 0.3;
 
    double mb = 0.0;
 
    double xbox, ybox;
 
    while ( mb <= 0.0 )
    {
      pb = _calcprob->FermiMomentum( _pf );
 
      // effective b-quark mass
      mb = mB * mB + _mq * _mq - 2.0 * mB * sqrt( pb * pb + _mq * _mq );
    }
    mb = sqrt( mb );
 
    //    cout << "b-quark momentum = " << pb << " mass = " <<  mb << endl;
 
    // generate a dilepton invariant mass
 
    double s    = 0.0;
    double smin = 4.0 * ml * ml;
    double smax = ( mb - _ms ) * ( mb - _ms );
 
    while ( s == 0.0 )
    {
      xbox = EvtRandom::Flat( smin, smax );
      ybox = EvtRandom::Flat( _dGdsProbMax );
      if ( ybox < _calcprob->dGdsProb( mb, _ms, ml, xbox ) ) { s = xbox; }
    }
 
    //    cout << "dGdsProb(s) = " << _calcprob->dGdsProb(mb, _ms, ml, s)
    //         << " for s = " << s << endl;
 
    // two-body decay of b quark at rest into s quark and dilepton pair:
    // b -> s (ll)
 
    EvtVector4R p4sdilep[2];
 
    double msdilep[2];
    msdilep[0] = _ms;
    msdilep[1] = sqrt( s );
 
    EvtGenKine::PhaseSpace( 2, msdilep, p4sdilep, mb );
 
    // generate dilepton decay with the expected asymmetry: (ll) -> l+ l-
 
    EvtVector4R p4ll[2];
 
    double mll[2];
    mll[0] = ml;
    mll[1] = ml;
 
    double tmp = 0.0;
 
    while ( tmp == 0.0 )
    {
      // (ll) -> l+ l- decay in dilepton rest frame
 
      EvtGenKine::PhaseSpace( 2, mll, p4ll, msdilep[1] );
 
      // boost to b-quark rest frame
 
      p4ll[0] = boostTo( p4ll[0], p4sdilep[1] );
      p4ll[1] = boostTo( p4ll[1], p4sdilep[1] );
 
      // compute kinematical variable u
 
      EvtVector4R p4slp = p4sdilep[0] + p4ll[0];
      EvtVector4R p4sln = p4sdilep[0] + p4ll[1];
 
      double u = p4slp.mass2() - p4sln.mass2();
 
      ybox = EvtRandom::Flat( _dGdsdupProbMax );
 
      double prob = _calcprob->dGdsdupProb( mb, _ms, ml, s, u );
      if ( prob > _dGdsdupProbMax && nmsg < 20 )
      {
        report( INFO, "EvtGen" )
            << "d2gdsdup GT d2gdsdup_max:" << prob << " " << _dGdsdupProbMax
            << " for s = " << s << " u = " << u << " mb = " << mb << endl;
        nmsg++;
      }
      if ( ybox < prob )
      {
        tmp = 1.0;
        //        cout << "dGdsdupProb(s) = " << prob
        //             << " for u = " << u << endl;
      }
    }
 
    // assign 4-momenta to valence quarks inside B meson in B rest frame
 
    double phi   = EvtRandom::Flat( EvtConst::twoPi );
    double costh = EvtRandom::Flat( -1.0, 1.0 );
    double sinth = sqrt( 1.0 - costh * costh );
 
    // b-quark four-momentum in B meson rest frame
 
    EvtVector4R p4b( sqrt( mb * mb + pb * pb ), pb * sinth * sin( phi ),
                     pb * sinth * cos( phi ), pb * costh );
 
    // B meson in its rest frame
    //
    //    EvtVector4R p4B(mB, 0.0, 0.0, 0.0);
    //
    // boost B meson to b-quark rest frame
    //
    //    p4B = boostTo(p4B, p4b);
    //
    //    cout << " B meson mass in b-quark rest frame = " << p4B.mass() << endl;
 
    // boost s, l+ and l- to B meson rest frame
 
    //    EvtVector4R p4s = boostTo(p4sdilep[0], p4B);
    //    p4leptonp       = boostTo(p4ll[0],     p4B);
    //    p4leptonn       = boostTo(p4ll[1],     p4B);
 
    EvtVector4R p4s = boostTo( p4sdilep[0], p4b );
    p4leptonp       = boostTo( p4ll[0], p4b );
    p4leptonn       = boostTo( p4ll[1], p4b );
 
    // spectator quark in B meson rest frame
 
    EvtVector4R p4q( sqrt( pb * pb + _mq * _mq ), -p4b.get( 1 ), -p4b.get( 2 ),
                     -p4b.get( 3 ) );
 
    // hadron system in B meson rest frame
 
    p4xhadron   = p4s + p4q;
    xhadronMass = p4xhadron.mass();
 
    //    cout << "Xs mass = " << xhadronMass << " trial " << im << endl;
  }
 
  // initialize the decay products
 
  xhadron->init( getDaug( 0 ), p4xhadron );
 
  // For B-bar mesons (i.e. containing a b quark) we have the normal
  // order of leptons
  if ( p->getId() == EvtPDL::getId( "anti-B0" ) || p->getId() == EvtPDL::getId( "B-" ) )
  {
    leptonp->init( getDaug( 1 ), p4leptonp );
    leptonn->init( getDaug( 2 ), p4leptonn );
  }
  // For B mesons (i.e. containing a b-bar quark) we need to flip the
  // role of the positive and negative leptons in order to produce the
  // correct forward-backward asymmetry between the two leptons
  else
  {
    leptonp->init( getDaug( 1 ), p4leptonn );
    leptonn->init( getDaug( 2 ), p4leptonp );
  }
 
  return;
}

getName()

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

virtual

Implements EvtDecayBase.

Definition at line 39 of file EvtBtoXsll.cc.

{ model_name = "BTOXSLL"; }

init()

void EvtBtoXsll::init ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 43 of file EvtBtoXsll.cc.

                      {
 
  // check that there are no arguments
 
  checkNArg( 0, 4, 5 );
 
  checkNDaug( 3 );
 
  // Check that the two leptons are the same type
 
  EvtId lepton1type = getDaug( 1 );
  EvtId lepton2type = getDaug( 2 );
 
  int etyp   = 0;
  int mutyp  = 0;
  int tautyp = 0;
  if ( lepton1type == EvtPDL::getId( "e+" ) || lepton1type == EvtPDL::getId( "e-" ) )
  { etyp++; }
  if ( lepton2type == EvtPDL::getId( "e+" ) || lepton2type == EvtPDL::getId( "e-" ) )
  { etyp++; }
  if ( lepton1type == EvtPDL::getId( "mu+" ) || lepton1type == EvtPDL::getId( "mu-" ) )
  { mutyp++; }
  if ( lepton2type == EvtPDL::getId( "mu+" ) || lepton2type == EvtPDL::getId( "mu-" ) )
  { mutyp++; }
  if ( lepton1type == EvtPDL::getId( "tau+" ) || lepton1type == EvtPDL::getId( "tau-" ) )
  { tautyp++; }
  if ( lepton2type == EvtPDL::getId( "tau+" ) || lepton2type == EvtPDL::getId( "tau-" ) )
  { tautyp++; }
 
  if ( etyp != 2 && mutyp != 2 && tautyp != 2 )
  {
 
    report( ERROR, "EvtGen" ) << "Expect two leptons of the same type in EvtBtoXsll.cc\n";
    ::abort();
  }
 
  // Check that the second and third entries are leptons with positive
  // and negative charge, respectively
 
  int lpos = 0;
  int lneg = 0;
  if ( lepton1type == EvtPDL::getId( "e+" ) || lepton1type == EvtPDL::getId( "mu+" ) ||
       lepton1type == EvtPDL::getId( "tau+" ) )
  { lpos++; }
  if ( lepton2type == EvtPDL::getId( "e-" ) || lepton2type == EvtPDL::getId( "mu-" ) ||
       lepton2type == EvtPDL::getId( "tau-" ) )
  { lneg++; }
 
  if ( lpos != 1 || lneg != 1 )
  {
 
    report( ERROR, "EvtGen" ) << "Expect 2nd and 3rd particles to be positive and negative "
                                 "leptons in EvtBtoXsll.cc\n";
    ::abort();
  }
 
  _mb    = 4.8;
  _ms    = 0.2;
  _mq    = 0.;
  _pf    = 0.41;
  _mxmin = 1.1;
  if ( getNArg() == 4 )
  {
    // b-quark mass
    _mb = getArg( 0 );
    // s-quark mass
    _ms = getArg( 1 );
    // spectator quark mass
    _mq = getArg( 2 );
    // Fermi motion parameter
    _pf = getArg( 3 );
  }
  if ( getNArg() == 5 ) { _mxmin = getArg( 4 ); }
 
  _calcprob = new EvtBtoXsllUtil;
 
  double ml = EvtPDL::getMeanMass( getDaug( 1 ) );
 
  // determine the maximum probability density from dGdsProb
 
  int    i, j;
  int    nsteps   = 100;
  double s        = 0.0;
  double smin     = 4.0 * ml * ml;
  double smax     = ( _mb - _ms ) * ( _mb - _ms );
  double probMax  = -10000.0;
  double sProbMax = -10.0;
  double uProbMax = -10.0;
 
  for ( i = 0; i < nsteps; i++ )
  {
    s           = smin + ( i + 0.0005 ) * ( smax - smin ) / (double)nsteps;
    double prob = _calcprob->dGdsProb( _mb, _ms, ml, s );
    if ( prob > probMax )
    {
      sProbMax = s;
      probMax  = prob;
    }
  }
 
  _dGdsProbMax = probMax;
 
  if ( verbose() )
  {
    report( INFO, "EvtGen" ) << "dGdsProbMax = " << probMax << " for s = " << sProbMax << endl;
  }
 
  // determine the maximum probability density from dGdsdupProb
 
  probMax  = -10000.0;
  sProbMax = -10.0;
 
  for ( i = 0; i < nsteps; i++ )
  {
    s = smin + ( i + 0.0005 ) * ( smax - smin ) / (double)nsteps;
    double umax =
        sqrt( ( s - ( _mb + _ms ) * ( _mb + _ms ) ) * ( s - ( _mb - _ms ) * ( _mb - _ms ) ) );
    for ( j = 0; j < nsteps; j++ )
    {
      double u    = -umax + ( j + 0.0005 ) * ( 2.0 * umax ) / (double)nsteps;
      double prob = _calcprob->dGdsdupProb( _mb, _ms, ml, s, u );
      if ( prob > probMax )
      {
        sProbMax = s;
        uProbMax = u;
        probMax  = prob;
      }
    }
  }
 
  _dGdsdupProbMax = probMax;
 
  if ( verbose() )
  {
    report( INFO, "EvtGen" ) << "dGdsdupProbMax = " << probMax << " for s = " << sProbMax
                             << " and u = " << uProbMax << endl;
  }
}

initProbMax()

void EvtBtoXsll::initProbMax ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 182 of file EvtBtoXsll.cc.

{ noProbMax(); }

---

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

• EvtBtoXsll.hh
• EvtBtoXsll.cc