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().

EvtAbsLineShape() [2/3]

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

Definition at line 44 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 41 of file EvtAbsLineShape.cc.

                                  {
}

EvtAbsLineShape() [3/3]

EvtAbsLineShape::EvtAbsLineShape

(

const EvtAbsLineShape &

x

)

Definition at line 69 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 71 of file EvtAbsLineShape.hh.

{ _addFactorPn = factor;}

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

clone()

EvtAbsLineShape * EvtAbsLineShape::clone ( )

virtual

Reimplemented in EvtFlatLineShape, EvtManyDeltaFuncLineShape, and EvtRelBreitWignerBarrierFact.

Definition at line 96 of file EvtAbsLineShape.cc.

                                        {
 
  return new EvtAbsLineShape(*this);
}

fixForSP8()

void EvtAbsLineShape::fixForSP8 ( )

inline

Definition at line 87 of file EvtAbsLineShape.hh.

{ _applyFixForSP8=true;}

Referenced by EvtPartProp::fixLSForSP8().

getMass()

double EvtAbsLineShape::getMass ( )

inline

Definition at line 41 of file EvtAbsLineShape.hh.

{return _mass;}

Referenced by EvtPartProp::getMass(), and EvtPartProp::newLineShape().

getMassMax()

double EvtAbsLineShape::getMassMax ( )

inline

Definition at line 43 of file EvtAbsLineShape.hh.

{return _massMax;}

Referenced by EvtPartProp::getMassMax().

getMassMin()

double EvtAbsLineShape::getMassMin ( )

inline

Definition at line 42 of file EvtAbsLineShape.hh.

{return _massMin;} 

Referenced by EvtPartProp::getMassMin().

getMassProb()

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

virtual

Reimplemented in EvtFlatLineShape, EvtManyDeltaFuncLineShape, and EvtRelBreitWignerBarrierFact.

Definition at line 171 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 EvtPartProp::getMassProb(), and EvtRelBreitWignerBarrierFact::getMassProb().

getMaxRange()

double EvtAbsLineShape::getMaxRange ( )

inline

Definition at line 44 of file EvtAbsLineShape.hh.

{return _maxRange;}

Referenced by EvtPartProp::getMaxRange(), and EvtPartProp::newLineShape().

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 119 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(), and EvtPartProp::getRandMass().

getSpinType()

EvtSpinType::spintype EvtAbsLineShape::getSpinType ( )

inline

Definition at line 46 of file EvtAbsLineShape.hh.

{return _spin;}

Referenced by EvtPartProp::newLineShape().

getWidth()

double EvtAbsLineShape::getWidth ( )

inline

Definition at line 45 of file EvtAbsLineShape.hh.

{return _width;} 

Referenced by EvtPartProp::getWidth(), and EvtPartProp::newLineShape().

includeBirthFactor()

void EvtAbsLineShape::includeBirthFactor ( bool  yesno )

inline

Definition at line 70 of file EvtAbsLineShape.hh.

{ _includeBirthFact = yesno; }

Referenced by EvtPartProp::includeBirthFactor().

includeDecayFactor()

void EvtAbsLineShape::includeDecayFactor ( bool  yesno )

inline

Definition at line 72 of file EvtAbsLineShape.hh.

{ _includeDecayFact = yesno; }

Referenced by EvtPartProp::includeDecayFactor().

operator=()

EvtAbsLineShape & EvtAbsLineShape::operator=

(

const EvtAbsLineShape &

x

)

Definition at line 82 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 69 of file EvtAbsLineShape.hh.

{};

Referenced by EvtPartProp::reSetBlatt().

reSetMass()

void EvtAbsLineShape::reSetMass ( double  mass )

inline

Definition at line 50 of file EvtAbsLineShape.hh.

{ _mass=mass;}

Referenced by EvtPartProp::reSetMass().

reSetMassMax()

void EvtAbsLineShape::reSetMassMax ( double  mass )

inline

Definition at line 68 of file EvtAbsLineShape.hh.

{ _massMax=mass;}

Referenced by EvtPartProp::reSetMassMax().

reSetMassMin()

void EvtAbsLineShape::reSetMassMin ( double  mass )

inline

Definition at line 67 of file EvtAbsLineShape.hh.

{ _massMin=mass;}

Referenced by EvtPartProp::reSetMassMin().

reSetWidth()

void EvtAbsLineShape::reSetWidth ( double  width )

inline

Definition at line 51 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;
      // ---> //
  }

Referenced by EvtPartProp::reSetWidth().

rollMass()

double EvtAbsLineShape::rollMass ( )

virtual

Definition at line 102 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;
  }
}

Referenced by EvtPartProp::rollMass().

setPWForBirthL()

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

inline

Definition at line 78 of file EvtAbsLineShape.hh.

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

Referenced by EvtPartProp::setPWForBirthL().

setPWForDecay()

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

inline

Definition at line 73 of file EvtAbsLineShape.hh.

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

Referenced by EvtPartProp::setPWForDecay().

Member Data Documentation

_addFactorPn

double EvtAbsLineShape::_addFactorPn

protected

Definition at line 92 of file EvtAbsLineShape.hh.

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

_applyFixForSP8

bool EvtAbsLineShape::_applyFixForSP8

protected

Definition at line 112 of file EvtAbsLineShape.hh.

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

_includeBirthFact

bool EvtAbsLineShape::_includeBirthFact

protected

Definition at line 91 of file EvtAbsLineShape.hh.

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

_includeDecayFact

bool EvtAbsLineShape::_includeDecayFact

protected

Definition at line 90 of file EvtAbsLineShape.hh.

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

_mass

double EvtAbsLineShape::_mass

protected

Definition at line 93 of file EvtAbsLineShape.hh.

Referenced by EvtAbsLineShape(), EvtFlatLineShape::EvtFlatLineShape(), 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 95 of file EvtAbsLineShape.hh.

Referenced by EvtAbsLineShape(), EvtFlatLineShape::EvtFlatLineShape(), EvtManyDeltaFuncLineShape::EvtManyDeltaFuncLineShape(), 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 94 of file EvtAbsLineShape.hh.

Referenced by EvtAbsLineShape(), EvtFlatLineShape::EvtFlatLineShape(), EvtManyDeltaFuncLineShape::EvtManyDeltaFuncLineShape(), 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 97 of file EvtAbsLineShape.hh.

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

_spin

EvtSpinType::spintype EvtAbsLineShape::_spin

protected

Definition at line 110 of file EvtAbsLineShape.hh.

Referenced by EvtAbsLineShape(), EvtFlatLineShape::EvtFlatLineShape(), 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 107 of file EvtAbsLineShape.hh.

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

_userSetBirthPar

std::vector<EvtId> EvtAbsLineShape::_userSetBirthPar

protected

Definition at line 107 of file EvtAbsLineShape.hh.

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

_userSetBirthPW

std::vector<int> EvtAbsLineShape::_userSetBirthPW

protected

Definition at line 108 of file EvtAbsLineShape.hh.

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

_userSetPW

std::vector<int> EvtAbsLineShape::_userSetPW

protected

Definition at line 104 of file EvtAbsLineShape.hh.

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

_userSetPWD1

std::vector<EvtId> EvtAbsLineShape::_userSetPWD1

protected

Definition at line 103 of file EvtAbsLineShape.hh.

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

_userSetPWD2

std::vector<EvtId> EvtAbsLineShape::_userSetPWD2

protected

Definition at line 103 of file EvtAbsLineShape.hh.

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

_width

double EvtAbsLineShape::_width

protected

Definition at line 96 of file EvtAbsLineShape.hh.

Referenced by EvtAbsLineShape(), EvtFlatLineShape::EvtFlatLineShape(), 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