EvtMassAmp Class Reference

#include <EvtMassAmp.hh>

Inheritance diagram for EvtMassAmp:

List of all members.

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)
EvtComplex	evaluate (const EvtPoint1D &p) const
Private Attributes
EvtPropBreitWignerRel	_prop
EvtTwoBodyVertex	_vd
EvtTwoBodyVertex *	_vb
double	_addfactor
bool	_useBirthFact
bool	_useDeathFact
bool	_useBirthFactFF
bool	_useDeathFactFF
bool	_fixUpMassForMax

---

Detailed Description

Definition at line 29 of file EvtMassAmp.hh.

---

Constructor & Destructor Documentation

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

Definition at line 20 of file EvtMassAmp.cc.

Referenced by clone().

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

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

(

)

[virtual]

Definition at line 41 of file EvtMassAmp.cc.

References _vb.

{
  if(_vb) delete _vb;
}

---

Member Function Documentation

void EvtMassAmp::addBirthFact

(

)

[inline]

Definition at line 46 of file EvtMassAmp.hh.

References _useBirthFact.

Referenced by EvtRelBreitWignerBarrierFact::getRandMass().

{ _useBirthFact = true; }

void EvtMassAmp::addBirthFactFF

(

)

[inline]

Definition at line 48 of file EvtMassAmp.hh.

References _useBirthFactFF.

Referenced by EvtRelBreitWignerBarrierFact::getRandMass().

{ _useBirthFactFF = true; }

void EvtMassAmp::addDeathFact

(

)

[inline]

Definition at line 47 of file EvtMassAmp.hh.

References _useDeathFact.

Referenced by EvtRelBreitWignerBarrierFact::getRandMass().

{ _useDeathFact = true; }

void EvtMassAmp::addDeathFactFF

(

)

[inline]

Definition at line 49 of file EvtMassAmp.hh.

References _useDeathFactFF.

Referenced by EvtRelBreitWignerBarrierFact::getRandMass().

{ _useDeathFactFF = true; }

void EvtMassAmp::addFactorPn

(

double

factor

)

[inline]

Definition at line 51 of file EvtMassAmp.hh.

References _addfactor.

Referenced by EvtRelBreitWignerBarrierFact::getRandMass().

{_addfactor = factor;}

EvtComplex EvtMassAmp::amplitude

(

const EvtPoint1D &

p

)

const [virtual]

Implements EvtAmplitude< EvtPoint1D >.

Definition at line 47 of file EvtMassAmp.cc.

References _addfactor, _fixUpMassForMax, _prop, _useBirthFact, _useBirthFactFF, _useDeathFact, _useDeathFactFF, _vb, _vd, EvtTwoBodyKine::AB, EvtTwoBodyVertex::formFactor(), EvtPropagator::g0(), EvtTwoBodyVertex::L(), EvtPropagator::m0(), EvtTwoBodyVertex::mA(), EvtTwoBodyVertex::mAB(), EvtTwoBodyVertex::mB(), EvtTwoBodyVertex::pD(), EvtTwoBodyVertex::phaseSpaceFactor(), EvtPoint1D::value(), and EvtTwoBodyVertex::widthFactor().

{
  // 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;
}

virtual EvtAmplitude<EvtPoint1D>* EvtMassAmp::clone

(

)

const [inline, virtual]

Implements EvtAmplitude< EvtPoint1D >.

Definition at line 38 of file EvtMassAmp.hh.

References EvtMassAmp().

  { return new EvtMassAmp(*this); }

EvtComplex EvtAmplitude< EvtPoint1D >::evaluate

(

const EvtPoint1D &

p

)

const [inline, inherited]

Definition at line 26 of file EvtAmplitude.hh.

References EvtAmplitude< T >::amplitude().

  {
    EvtComplex ret(0.,0.);
    if(p.isValid()) ret = amplitude(p);
    return ret;
  }

void EvtMassAmp::fixUpMassForMax

(

)

[inline]

Definition at line 50 of file EvtMassAmp.hh.

References _fixUpMassForMax.

{ _fixUpMassForMax = true;}

void EvtMassAmp::setBirthVtx

(

const EvtTwoBodyVertex &

vb

)

[inline]

Definition at line 41 of file EvtMassAmp.hh.

References _vb.

Referenced by EvtRelBreitWignerBarrierFact::getRandMass().

  {
    _vb = new EvtTwoBodyVertex(vb);
  }

---

Member Data Documentation

double EvtMassAmp::_addfactor [private]

Definition at line 58 of file EvtMassAmp.hh.

Referenced by addFactorPn(), and amplitude().

bool EvtMassAmp::_fixUpMassForMax [private]

Definition at line 64 of file EvtMassAmp.hh.

Referenced by amplitude(), and fixUpMassForMax().

EvtPropBreitWignerRel EvtMassAmp::_prop [private]

Definition at line 55 of file EvtMassAmp.hh.

Referenced by amplitude().

bool EvtMassAmp::_useBirthFact [private]

Definition at line 60 of file EvtMassAmp.hh.

Referenced by addBirthFact(), and amplitude().

bool EvtMassAmp::_useBirthFactFF [private]

Definition at line 62 of file EvtMassAmp.hh.

Referenced by addBirthFactFF(), and amplitude().

bool EvtMassAmp::_useDeathFact [private]

Definition at line 61 of file EvtMassAmp.hh.

Referenced by addDeathFact(), and amplitude().

bool EvtMassAmp::_useDeathFactFF [private]

Definition at line 63 of file EvtMassAmp.hh.

Referenced by addDeathFactFF(), and amplitude().

EvtTwoBodyVertex* EvtMassAmp::_vb [private]

Definition at line 57 of file EvtMassAmp.hh.

Referenced by amplitude(), setBirthVtx(), and ~EvtMassAmp().

EvtTwoBodyVertex EvtMassAmp::_vd [private]

Definition at line 56 of file EvtMassAmp.hh.

Referenced by amplitude().

---