EvtAbsLineShape Class Reference

#include <EvtAbsLineShape.hh>

Inheritance diagram for EvtAbsLineShape:

Public Member Functions
	EvtAbsLineShape ()
	EvtAbsLineShape (double mass, double width, double maxRange, EvtSpinType::spintype sp)
virtual	~EvtAbsLineShape ()
EvtAbsLineShape &	operator= (const EvtAbsLineShape &x)
	EvtAbsLineShape (const EvtAbsLineShape &x)
double	getMass ()
double	getMassMin ()
double	getMassMax ()
double	getMaxRange ()
double	getWidth ()
EvtSpinType::spintype	getSpinType ()
virtual double	rollMass ()
virtual EvtAbsLineShape *	clone ()
void	reSetMass (double mass)
void	reSetWidth (double width)
void	reSetMassMin (double mass)
void	reSetMassMax (double mass)
virtual void	reSetBlatt (double blatt)
void	includeBirthFactor (bool yesno)
void	addFactorPn (double factor=0.)
void	includeDecayFactor (bool yesno)
void	setPWForDecay (int spin, EvtId d1, EvtId d2)
void	setPWForBirthL (int spin, EvtId par, EvtId othD)
virtual double	getRandMass (EvtId *parId, int nDaug, EvtId *dauId, EvtId *othDaugId, double maxMass, double *dauMasses)
virtual double	getMassProb (double mass, double massPar, int nDaug, double *massDau)
void	fixForSP8 ()

Protected Attributes
bool	_includeDecayFact
bool	_includeBirthFact
double	_addFactorPn
double	_mass
double	_massMin
double	_massMax
double	_width
double	_maxRange
std::vector< EvtId >	_userSetPWD1
std::vector< EvtId >	_userSetPWD2
std::vector< int >	_userSetPW
std::vector< EvtId >	_userSetBirthPar
std::vector< EvtId >	_userSetBirthOthD
std::vector< int >	_userSetBirthPW
EvtSpinType::spintype	_spin
bool	_applyFixForSP8

Detailed Description

Definition at line 31 of file EvtAbsLineShape.hh.

Constructor & Destructor Documentation

EvtAbsLineShape() [1/3]

EvtAbsLineShape::EvtAbsLineShape

(

)

Definition at line 38 of file EvtAbsLineShape.cc.

{}

Referenced by clone(), EvtFlatLineShape::clone(), EvtManyDeltaFuncLineShape::clone(), EvtRelBreitWignerBarrierFact::clone(), EvtAbsLineShape(), EvtRelBreitWignerBarrierFact::EvtRelBreitWignerBarrierFact(), EvtRelBreitWignerBarrierFact::EvtRelBreitWignerBarrierFact(), and operator=().

EvtAbsLineShape() [2/3]

EvtAbsLineShape::EvtAbsLineShape	(	double	mass,
		double	width,
		double	maxRange,
		EvtSpinType::spintype	sp )

Definition at line 42 of file EvtAbsLineShape.cc.

                                                           {
 
  _includeDecayFact = false;
  _includeBirthFact = false;
  _applyFixForSP8   = false;
  _mass             = mass;
  _width            = width;
  _spin             = sp;
  _maxRange         = maxRange;
  double maxdelta   = 15.0 * width;
  // if ( width>0.001 ) {
  //   if ( 5.0*width < 0.6 ) maxdelta = 0.6;
  // }
  if ( maxRange > 0.00001 )
  {
    _massMax = mass + maxdelta;
    _massMin = mass - maxRange;
  }
  else
  {
    _massMax = mass + maxdelta;
    _massMin = mass - 15.0 * width;
  }
  if ( _massMin < 0. ) _massMin = 0.;
  _massMax = mass + maxdelta;
}

~EvtAbsLineShape()

EvtAbsLineShape::~EvtAbsLineShape ( )

virtual

Definition at line 40 of file EvtAbsLineShape.cc.

{}

EvtAbsLineShape() [3/3]

EvtAbsLineShape::EvtAbsLineShape

(

const EvtAbsLineShape &

x

)

Definition at line 70 of file EvtAbsLineShape.cc.

                                                           {
 
  _includeDecayFact = x._includeDecayFact;
  _includeBirthFact = x._includeBirthFact;
  _mass             = x._mass;
  _massMax          = x._massMax;
  _massMin          = x._massMin;
  _width            = x._width;
  _spin             = x._spin;
  _maxRange         = x._maxRange;
  _applyFixForSP8   = x._applyFixForSP8;
}

Member Function Documentation

addFactorPn()

void EvtAbsLineShape::addFactorPn ( double factor = 0. )

inline

Definition at line 72 of file EvtAbsLineShape.hh.

{ _addFactorPn = factor; }

Referenced by EvtRelBreitWignerBarrierFact::getRandMass().

clone()

EvtAbsLineShape * EvtAbsLineShape::clone ( )

virtual

Reimplemented in EvtFlatLineShape, EvtManyDeltaFuncLineShape, and EvtRelBreitWignerBarrierFact.

Definition at line 97 of file EvtAbsLineShape.cc.

{ return new EvtAbsLineShape( *this ); }

Referenced by EvtPartProp::operator=().

fixForSP8()

void EvtAbsLineShape::fixForSP8 ( )

inline

Definition at line 89 of file EvtAbsLineShape.hh.

{ _applyFixForSP8 = true; }

getMass()

double EvtAbsLineShape::getMass ( )

inline

Definition at line 40 of file EvtAbsLineShape.hh.

{ return _mass; }

getMassMax()

double EvtAbsLineShape::getMassMax ( )

inline

Definition at line 42 of file EvtAbsLineShape.hh.

{ return _massMax; }

getMassMin()

double EvtAbsLineShape::getMassMin ( )

inline

Definition at line 41 of file EvtAbsLineShape.hh.

{ return _massMin; }

getMassProb()

double EvtAbsLineShape::getMassProb ( double mass, double massPar, int nDaug, double * massDau )

virtual

Reimplemented in EvtFlatLineShape, EvtManyDeltaFuncLineShape, and EvtRelBreitWignerBarrierFact.

Definition at line 173 of file EvtAbsLineShape.cc.

                                                       {
 
  double dTotMass = 0.;
  if ( nDaug > 1 )
  {
    int i;
    for ( i = 0; i < nDaug; i++ ) { dTotMass += massDau[i]; }
    // report(INFO,"EvtGen") << mass << " " << massPar << " " << dTotMass << " "<< endl;
    //     if ( (mass-dTotMass)<0.0001 ) return 0.;
    if ( ( mass < dTotMass ) ) return 0.;
  }
  if ( _width < 0.0001 ) return 1.;
 
  // no parent - lets not panic
  if ( massPar > 0.0000000001 )
  {
    if ( mass > massPar ) return 0.;
  }
  // Otherwise return the right value.
  // Fortunately we have generated events according to a non-rel BW, so
  // just return..
  // EvtPropBreitWigner bw(_mass,_width);
  // EvtPropFactor<EvtTwoBodyVertex> f(bw);
  // EvtComplex fm=f.eval(mass);
  // EvtComplex fm0=f.eval(_mass);
  // return (abs(fm)*abs(fm))/(abs(fm0)*abs(fm0));
  return 1.0;
}

Referenced by EvtRelBreitWignerBarrierFact::getMassProb().

getMaxRange()

double EvtAbsLineShape::getMaxRange ( )

inline

Definition at line 43 of file EvtAbsLineShape.hh.

{ return _maxRange; }

getRandMass()

double EvtAbsLineShape::getRandMass ( EvtId * parId, int nDaug, EvtId * dauId, EvtId * othDaugId, double maxMass, double * dauMasses )

virtual

Reimplemented in EvtFlatLineShape, EvtManyDeltaFuncLineShape, and EvtRelBreitWignerBarrierFact.

Definition at line 115 of file EvtAbsLineShape.cc.

                                                                         {
 
  if ( _width < 0.0001 ) return _mass;
  // its not flat - but generated according to a BW
 
  if ( maxMass > 0 && maxMass < _massMin )
  {
    report( ERROR, "EvtGen" ) << "In EvtAbsLineShape::getRandMass" << endl;
    report( ERROR, "EvtGen" ) << "Decaying particle " << EvtPDL::name( *parId )
                              << " with mass " << maxMass << endl;
    report( ERROR, "EvtGen" ) << " to particle"
                              << " with a minimal mass of " << _massMin << endl;
  }
 
  double mMin = _massMin;
  double mMax = _massMax;
  if ( maxMass > -0.5 && maxMass < mMax ) mMax = maxMass;
  double ymin = atan( 2.0 * ( mMin - _mass ) / _width );
  double ymax = atan( 2.0 * ( mMax - _mass ) / _width );
 
loop:
  double themass = ( _mass + ( ( _width / 2.0 ) * tan( EvtRandom::Flat( ymin, ymax ) ) ) );
 
  if ( fabs( _addFactorPn ) > 0.00000001 )
  { // addFactorPn, pingrg-2010-1-10
    double massOthD   = EvtPDL::getMeanMass( *othDaugId );
    double massParent = EvtPDL::getMeanMass( *parId );
    double phsp, maxp, maxp1, maxp2;
    if ( themass + massOthD < massParent )
    {
      EvtTwoBodyVertex vb( themass, massOthD, massParent, 1.0 );
      phsp = vb.pD();
    }
    else { return themass; }
 
    if ( ( massOthD + mMax ) < massParent )
    {
      EvtTwoBodyVertex vb1( massOthD, mMax, massParent, 1 );
      EvtTwoBodyVertex vb2( massOthD, mMin, massParent, 1 );
      maxp  = vb1.pD();
      maxp1 = pow( maxp, _addFactorPn * 2.0 );
      maxp  = vb2.pD();
      maxp2 = pow( maxp, _addFactorPn * 2.0 );
      maxp  = max( maxp1, maxp2 );
    }
    else { return themass; }
 
    double wt  = pow( phsp, _addFactorPn * 2.0 ) / maxp;
    double rdm = EvtRandom::Flat( 0.0, 1.0 );
 
    if ( rdm > wt ) goto loop;
  }
 
  return themass;
  //  return EvtRandom::Flat(_massMin,_massMax);
}

Referenced by EvtRelBreitWignerBarrierFact::getRandMass().

getSpinType()

EvtSpinType::spintype EvtAbsLineShape::getSpinType ( )

inline

Definition at line 45 of file EvtAbsLineShape.hh.

{ return _spin; }

getWidth()

double EvtAbsLineShape::getWidth ( )

inline

Definition at line 44 of file EvtAbsLineShape.hh.

{ return _width; }

includeBirthFactor()

void EvtAbsLineShape::includeBirthFactor ( bool yesno )

inline

Definition at line 71 of file EvtAbsLineShape.hh.

{ _includeBirthFact = yesno; }

includeDecayFactor()

void EvtAbsLineShape::includeDecayFactor ( bool yesno )

inline

Definition at line 73 of file EvtAbsLineShape.hh.

{ _includeDecayFact = yesno; }

operator=()

EvtAbsLineShape & EvtAbsLineShape::operator=

(

const EvtAbsLineShape &

x

)

Definition at line 83 of file EvtAbsLineShape.cc.

                                                                    {
 
  _includeDecayFact = x._includeDecayFact;
  _includeBirthFact = x._includeBirthFact;
  _mass             = x._mass;
  _massMax          = x._massMax;
  _massMin          = x._massMin;
  _width            = x._width;
  _spin             = x._spin;
  _maxRange         = x._maxRange;
  _applyFixForSP8   = x._applyFixForSP8;
  return *this;
}

reSetBlatt()

virtual void EvtAbsLineShape::reSetBlatt ( double blatt )

inlinevirtual

Reimplemented in EvtRelBreitWignerBarrierFact.

Definition at line 70 of file EvtAbsLineShape.hh.

{};

reSetMass()

void EvtAbsLineShape::reSetMass ( double mass )

inline

Definition at line 49 of file EvtAbsLineShape.hh.

{ _mass = mass; }

reSetMassMax()

void EvtAbsLineShape::reSetMassMax ( double mass )

inline

Definition at line 69 of file EvtAbsLineShape.hh.

{ _massMax = mass; }

reSetMassMin()

void EvtAbsLineShape::reSetMassMin ( double mass )

inline

Definition at line 68 of file EvtAbsLineShape.hh.

{ _massMin = mass; }

reSetWidth()

void EvtAbsLineShape::reSetWidth ( double width )

inline

Definition at line 50 of file EvtAbsLineShape.hh.

                                  {
    _width = width;
    // <--- added by L.L. Wang to fix a bug
    double maxdelta = 15.0 * width;
    if ( _maxRange > 0.00001 )
    {
      _massMax = _mass + maxdelta;
      _massMin = _mass - _maxRange;
    }
    else
    {
      _massMax = _mass + maxdelta;
      _massMin = _mass - 15.0 * width;
    }
    if ( _massMin < 0. ) _massMin = 0.;
    _massMax = _mass + maxdelta;
    // ---> //
  }

rollMass()

double EvtAbsLineShape::rollMass ( )

virtual

Definition at line 99 of file EvtAbsLineShape.cc.

                                 {
 
  double ymin, ymax;
  double temp;
 
  if ( _width < 0.0001 ) { return _mass; }
  else
  {
    ymin = atan( 2.0 * ( _massMin - _mass ) / _width );
    ymax = atan( 2.0 * ( _massMax - _mass ) / _width );
 
    temp = ( _mass + ( ( _width / 2.0 ) * tan( EvtRandom::Flat( ymin, ymax ) ) ) );
 
    return temp;
  }
}

setPWForBirthL()

void EvtAbsLineShape::setPWForBirthL ( int spin, EvtId par, EvtId othD )

inline

Definition at line 79 of file EvtAbsLineShape.hh.

                                                         {
    _userSetBirthPW.push_back( spin );
    _userSetBirthOthD.push_back( othD );
    _userSetBirthPar.push_back( par );
  }

setPWForDecay()

void EvtAbsLineShape::setPWForDecay ( int spin, EvtId d1, EvtId d2 )

inline

Definition at line 74 of file EvtAbsLineShape.hh.

                                                             {
    _userSetPW.push_back( spin );
    _userSetPWD1.push_back( d1 );
    _userSetPWD2.push_back( d2 );
  }

Member Data Documentation

_addFactorPn

double EvtAbsLineShape::_addFactorPn

protected

Definition at line 94 of file EvtAbsLineShape.hh.

Referenced by addFactorPn(), getRandMass(), and EvtRelBreitWignerBarrierFact::getRandMass().

_applyFixForSP8

bool EvtAbsLineShape::_applyFixForSP8

protected

Definition at line 114 of file EvtAbsLineShape.hh.

Referenced by EvtAbsLineShape(), EvtAbsLineShape(), fixForSP8(), EvtRelBreitWignerBarrierFact::getRandMass(), and operator=().

_includeBirthFact

bool EvtAbsLineShape::_includeBirthFact

protected

Definition at line 93 of file EvtAbsLineShape.hh.

Referenced by EvtAbsLineShape(), EvtAbsLineShape(), EvtRelBreitWignerBarrierFact::EvtRelBreitWignerBarrierFact(), EvtRelBreitWignerBarrierFact::EvtRelBreitWignerBarrierFact(), EvtRelBreitWignerBarrierFact::getRandMass(), includeBirthFactor(), operator=(), and EvtRelBreitWignerBarrierFact::operator=().

_includeDecayFact

bool EvtAbsLineShape::_includeDecayFact

protected

Definition at line 92 of file EvtAbsLineShape.hh.

Referenced by EvtAbsLineShape(), EvtAbsLineShape(), EvtRelBreitWignerBarrierFact::EvtRelBreitWignerBarrierFact(), EvtRelBreitWignerBarrierFact::EvtRelBreitWignerBarrierFact(), EvtRelBreitWignerBarrierFact::getRandMass(), includeDecayFactor(), operator=(), and EvtRelBreitWignerBarrierFact::operator=().

_mass

double EvtAbsLineShape::_mass

protected

Definition at line 95 of file EvtAbsLineShape.hh.

Referenced by EvtAbsLineShape(), EvtAbsLineShape(), EvtFlatLineShape::EvtFlatLineShape(), EvtFlatLineShape::EvtFlatLineShape(), EvtManyDeltaFuncLineShape::EvtManyDeltaFuncLineShape(), EvtManyDeltaFuncLineShape::EvtManyDeltaFuncLineShape(), EvtRelBreitWignerBarrierFact::EvtRelBreitWignerBarrierFact(), getMass(), getRandMass(), EvtRelBreitWignerBarrierFact::getRandMass(), operator=(), EvtFlatLineShape::operator=(), EvtManyDeltaFuncLineShape::operator=(), EvtRelBreitWignerBarrierFact::operator=(), reSetMass(), reSetWidth(), and rollMass().

_massMax

double EvtAbsLineShape::_massMax

protected

Definition at line 97 of file EvtAbsLineShape.hh.

Referenced by EvtAbsLineShape(), EvtAbsLineShape(), EvtFlatLineShape::EvtFlatLineShape(), EvtFlatLineShape::EvtFlatLineShape(), EvtManyDeltaFuncLineShape::EvtManyDeltaFuncLineShape(), EvtManyDeltaFuncLineShape::EvtManyDeltaFuncLineShape(), EvtRelBreitWignerBarrierFact::EvtRelBreitWignerBarrierFact(), EvtRelBreitWignerBarrierFact::EvtRelBreitWignerBarrierFact(), getMassMax(), getRandMass(), EvtFlatLineShape::getRandMass(), EvtManyDeltaFuncLineShape::getRandMass(), EvtRelBreitWignerBarrierFact::getRandMass(), operator=(), EvtFlatLineShape::operator=(), EvtManyDeltaFuncLineShape::operator=(), EvtRelBreitWignerBarrierFact::operator=(), reSetMassMax(), reSetWidth(), and rollMass().

_massMin

double EvtAbsLineShape::_massMin

protected

Definition at line 96 of file EvtAbsLineShape.hh.

Referenced by EvtAbsLineShape(), EvtAbsLineShape(), EvtFlatLineShape::EvtFlatLineShape(), EvtFlatLineShape::EvtFlatLineShape(), EvtManyDeltaFuncLineShape::EvtManyDeltaFuncLineShape(), EvtManyDeltaFuncLineShape::EvtManyDeltaFuncLineShape(), EvtRelBreitWignerBarrierFact::EvtRelBreitWignerBarrierFact(), EvtRelBreitWignerBarrierFact::EvtRelBreitWignerBarrierFact(), getMassMin(), getRandMass(), EvtFlatLineShape::getRandMass(), EvtManyDeltaFuncLineShape::getRandMass(), EvtRelBreitWignerBarrierFact::getRandMass(), operator=(), EvtFlatLineShape::operator=(), EvtManyDeltaFuncLineShape::operator=(), EvtRelBreitWignerBarrierFact::operator=(), reSetMassMin(), reSetWidth(), and rollMass().

_maxRange

double EvtAbsLineShape::_maxRange

protected

Definition at line 99 of file EvtAbsLineShape.hh.

Referenced by EvtAbsLineShape(), EvtAbsLineShape(), EvtFlatLineShape::EvtFlatLineShape(), EvtFlatLineShape::EvtFlatLineShape(), EvtManyDeltaFuncLineShape::EvtManyDeltaFuncLineShape(), EvtManyDeltaFuncLineShape::EvtManyDeltaFuncLineShape(), EvtRelBreitWignerBarrierFact::EvtRelBreitWignerBarrierFact(), EvtRelBreitWignerBarrierFact::EvtRelBreitWignerBarrierFact(), getMaxRange(), operator=(), EvtFlatLineShape::operator=(), EvtManyDeltaFuncLineShape::operator=(), EvtRelBreitWignerBarrierFact::operator=(), and reSetWidth().

_spin

EvtSpinType::spintype EvtAbsLineShape::_spin

protected

Definition at line 112 of file EvtAbsLineShape.hh.

Referenced by EvtAbsLineShape(), EvtAbsLineShape(), EvtFlatLineShape::EvtFlatLineShape(), EvtFlatLineShape::EvtFlatLineShape(), EvtManyDeltaFuncLineShape::EvtManyDeltaFuncLineShape(), EvtManyDeltaFuncLineShape::EvtManyDeltaFuncLineShape(), EvtRelBreitWignerBarrierFact::EvtRelBreitWignerBarrierFact(), EvtRelBreitWignerBarrierFact::getRandMass(), getSpinType(), operator=(), EvtFlatLineShape::operator=(), EvtManyDeltaFuncLineShape::operator=(), and EvtRelBreitWignerBarrierFact::operator=().

_userSetBirthOthD

std::vector<EvtId> EvtAbsLineShape::_userSetBirthOthD

protected

Definition at line 109 of file EvtAbsLineShape.hh.

Referenced by EvtRelBreitWignerBarrierFact::getRandMass(), and setPWForBirthL().

_userSetBirthPar

std::vector<EvtId> EvtAbsLineShape::_userSetBirthPar

protected

Definition at line 109 of file EvtAbsLineShape.hh.

Referenced by EvtRelBreitWignerBarrierFact::getRandMass(), and setPWForBirthL().

_userSetBirthPW

std::vector<int> EvtAbsLineShape::_userSetBirthPW

protected

Definition at line 110 of file EvtAbsLineShape.hh.

Referenced by EvtRelBreitWignerBarrierFact::getRandMass(), and setPWForBirthL().

_userSetPW

std::vector<int> EvtAbsLineShape::_userSetPW

protected

Definition at line 106 of file EvtAbsLineShape.hh.

Referenced by EvtRelBreitWignerBarrierFact::getRandMass(), and setPWForDecay().

_userSetPWD1

std::vector<EvtId> EvtAbsLineShape::_userSetPWD1

protected

Definition at line 105 of file EvtAbsLineShape.hh.

Referenced by EvtRelBreitWignerBarrierFact::getRandMass(), and setPWForDecay().

_userSetPWD2

std::vector<EvtId> EvtAbsLineShape::_userSetPWD2

protected

Definition at line 105 of file EvtAbsLineShape.hh.

Referenced by EvtRelBreitWignerBarrierFact::getRandMass(), and setPWForDecay().

_width

double EvtAbsLineShape::_width

protected

Definition at line 98 of file EvtAbsLineShape.hh.

Referenced by EvtAbsLineShape(), EvtAbsLineShape(), EvtFlatLineShape::EvtFlatLineShape(), EvtFlatLineShape::EvtFlatLineShape(), EvtManyDeltaFuncLineShape::EvtManyDeltaFuncLineShape(), EvtManyDeltaFuncLineShape::EvtManyDeltaFuncLineShape(), EvtRelBreitWignerBarrierFact::EvtRelBreitWignerBarrierFact(), getMassProb(), EvtRelBreitWignerBarrierFact::getMassProb(), getRandMass(), EvtManyDeltaFuncLineShape::getRandMass(), EvtRelBreitWignerBarrierFact::getRandMass(), getWidth(), operator=(), EvtFlatLineShape::operator=(), EvtManyDeltaFuncLineShape::operator=(), EvtRelBreitWignerBarrierFact::operator=(), reSetWidth(), and rollMass().

---

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

• EvtAbsLineShape.hh
• EvtAbsLineShape.cc