EvtMassAmp Class Reference

#include <EvtMassAmp.hh>

Inheritance diagram for EvtMassAmp:

Public Member Functions
	EvtMassAmp (const EvtPropBreitWignerRel &prop, const EvtTwoBodyVertex &vd)
	EvtMassAmp (const EvtMassAmp &other)
virtual	~EvtMassAmp ()
virtual EvtComplex	amplitude (const EvtPoint1D &p) const
virtual EvtAmplitude< EvtPoint1D > *	clone () const
void	setBirthVtx (const EvtTwoBodyVertex &vb)
void	addBirthFact ()
void	addDeathFact ()
void	addBirthFactFF ()
void	addDeathFactFF ()
void	fixUpMassForMax ()
void	addFactorPn (double factor)
Public Member Functions inherited from EvtAmplitude< EvtPoint1D >
	EvtAmplitude ()
	EvtAmplitude (const EvtAmplitude &)
virtual	~EvtAmplitude ()
EvtComplex	evaluate (const EvtPoint1D &p) const

Additional Inherited Members

Protected Member Functions inherited from EvtAmplitude< EvtPoint1D >

Detailed Description

Definition at line 29 of file EvtMassAmp.hh.

Constructor & Destructor Documentation

EvtMassAmp() [1/2]

EvtMassAmp::EvtMassAmp	(	const EvtPropBreitWignerRel &	prop,
		const EvtTwoBodyVertex &	vd )

Definition at line 20 of file EvtMassAmp.cc.

  : EvtAmplitude<EvtPoint1D>()
  ,_prop(prop), _vd(vd), _vb(0)
  ,_useBirthFact(false), _useDeathFact(false)
  ,_useBirthFactFF(false), _useDeathFactFF(false)
  ,_fixUpMassForMax(false),_addfactor(0.)
{}

Referenced by clone().

EvtMassAmp() [2/2]

EvtMassAmp::EvtMassAmp

(

const EvtMassAmp &

other

)

Definition at line 28 of file EvtMassAmp.cc.

  : EvtAmplitude<EvtPoint1D>(other)
  ,_prop(other._prop), _vd(other._vd)
  ,_vb(other._vb ? new EvtTwoBodyVertex(*other._vb) : 0)
  ,_useBirthFact(other._useBirthFact)
  ,_useDeathFact(other._useDeathFact)
  ,_useBirthFactFF(other._useBirthFactFF)
  ,_useDeathFactFF(other._useDeathFactFF)
  ,_fixUpMassForMax(other._fixUpMassForMax)
  ,_addfactor(other._addfactor)
{}

~EvtMassAmp()

EvtMassAmp::~EvtMassAmp ( )

virtual

Definition at line 41 of file EvtMassAmp.cc.

{
  if(_vb) delete _vb;
}

Member Function Documentation

addBirthFact()

void EvtMassAmp::addBirthFact ( )

inline

Definition at line 46 of file EvtMassAmp.hh.

{ _useBirthFact = true; }

Referenced by EvtRelBreitWignerBarrierFact::getRandMass().

addBirthFactFF()

void EvtMassAmp::addBirthFactFF ( )

inline

Definition at line 48 of file EvtMassAmp.hh.

{ _useBirthFactFF = true; }

Referenced by EvtRelBreitWignerBarrierFact::getRandMass().

addDeathFact()

void EvtMassAmp::addDeathFact ( )

inline

Definition at line 47 of file EvtMassAmp.hh.

{ _useDeathFact = true; }

Referenced by EvtRelBreitWignerBarrierFact::getRandMass().

addDeathFactFF()

void EvtMassAmp::addDeathFactFF ( )

inline

Definition at line 49 of file EvtMassAmp.hh.

{ _useDeathFactFF = true; }

Referenced by EvtRelBreitWignerBarrierFact::getRandMass().

addFactorPn()

void EvtMassAmp::addFactorPn ( double factor )

inline

Definition at line 51 of file EvtMassAmp.hh.

{_addfactor = factor;}

Referenced by EvtRelBreitWignerBarrierFact::getRandMass().

amplitude()

EvtComplex EvtMassAmp::amplitude ( const EvtPoint1D & p ) const

virtual

Implements EvtAmplitude< EvtPoint1D >.

Definition at line 47 of file EvtMassAmp.cc.

{
  // Modified vertex
 
  double m = p.value();
  // keep things from crashing..
 
  if ( m< (_vd.mA()+_vd.mB()) ) return EvtComplex(0.,0.);
 
  EvtTwoBodyKine vd(_vd.mA(),_vd.mB(),m);
  
  // Compute mass-dependent width for relativistic propagator
 
  EvtPropBreitWignerRel bw(_prop.m0(),_prop.g0()*_vd.widthFactor(vd)); 
  EvtComplex amp = bw.evaluate(m);
 
 
  // Birth vertex factors
 
  if(_useBirthFact) {
 
    assert(_vb);
    if ( (m+_vb->mB()) < _vb->mAB() ) {  
      EvtTwoBodyKine vb(m,_vb->mB(),_vb->mAB());
      double phsp= _vb->phaseSpaceFactor(vb,EvtTwoBodyKine::AB);
      amp *= phsp;
 
      if ( _fixUpMassForMax ) {
    amp *= sqrt((vb.p() / _vb->pD()));
      }
      if(_useBirthFactFF) {
    
    assert(_vb);
    amp *= _vb->formFactor(vb);
      }
    }
    else{
      if ( _vb->L() != 0 ) amp=0.;
    }
  }
 
 
  // Decay vertex factors
 
  if(_useDeathFact) {
    amp *= _vd.phaseSpaceFactor(vd,EvtTwoBodyKine::AB);
    if ( _fixUpMassForMax ) {
      amp *= sqrt((vd.p() / _vd.pD()));
    }
  }
  if(_useDeathFactFF) amp *= _vd.formFactor(vd);
 
  // AddFactorPn
 
  if(fabs(_addfactor) > 0.00000001){ //addFactorPn, pingrg-2011-1-6
    assert(_vb);
    double phsp;
    if ( (m+_vb->mB()) < _vb->mAB() ) {  
      EvtTwoBodyKine vb(m,_vb->mB(),_vb->mAB());
      phsp= vb.p(EvtTwoBodyKine::AB);
    } else{
      if ( _vb->L() != 0 ) phsp=0.;
    }
 
    amp *= pow(phsp,_addfactor);
    //   std::cout<<"EvtMassAmp::addFactorPn= "<<pow(phsp,_addfactor)<<std::endl;
   }
 
  return amp;
}

clone()

virtual EvtAmplitude< EvtPoint1D > * EvtMassAmp::clone ( ) const

inlinevirtual

Implements EvtAmplitude< EvtPoint1D >.

Definition at line 38 of file EvtMassAmp.hh.

  { return new EvtMassAmp(*this); }

fixUpMassForMax()

void EvtMassAmp::fixUpMassForMax ( )

inline

Definition at line 50 of file EvtMassAmp.hh.

{ _fixUpMassForMax = true;}

setBirthVtx()

void EvtMassAmp::setBirthVtx ( const EvtTwoBodyVertex & vb )

inline

Definition at line 41 of file EvtMassAmp.hh.

  {
    _vb = new EvtTwoBodyVertex(vb);
  }

Referenced by EvtRelBreitWignerBarrierFact::getRandMass().

---

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

• EvtMassAmp.hh
• EvtMassAmp.cc