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 ()
virtual	~EvtAmplitude ()
EvtComplex	evaluate (const EvtPoint1D &p) const

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(), and EvtMassAmp().

EvtMassAmp() [2/2]

EvtMassAmp::EvtMassAmp

(

const EvtMassAmp &

other

)

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

                        {
  if ( _vb ) delete _vb;
}

Member Function Documentation

addBirthFact()

void EvtMassAmp::addBirthFact ( )

inline

Definition at line 41 of file EvtMassAmp.hh.

{ _useBirthFact = true; }

Referenced by EvtRelBreitWignerBarrierFact::getRandMass().

addBirthFactFF()

void EvtMassAmp::addBirthFactFF ( )

inline

Definition at line 43 of file EvtMassAmp.hh.

{ _useBirthFactFF = true; }

Referenced by EvtRelBreitWignerBarrierFact::getRandMass().

addDeathFact()

void EvtMassAmp::addDeathFact ( )

inline

Definition at line 42 of file EvtMassAmp.hh.

{ _useDeathFact = true; }

Referenced by EvtRelBreitWignerBarrierFact::getRandMass().

addDeathFactFF()

void EvtMassAmp::addDeathFactFF ( )

inline

Definition at line 44 of file EvtMassAmp.hh.

{ _useDeathFactFF = true; }

Referenced by EvtRelBreitWignerBarrierFact::getRandMass().

addFactorPn()

void EvtMassAmp::addFactorPn ( double factor )

inline

Definition at line 46 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 48 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 37 of file EvtMassAmp.hh.

{ return new EvtMassAmp( *this ); }

fixUpMassForMax()

void EvtMassAmp::fixUpMassForMax ( )

inline

Definition at line 45 of file EvtMassAmp.hh.

{ _fixUpMassForMax = true; }

setBirthVtx()

void EvtMassAmp::setBirthVtx ( const EvtTwoBodyVertex & vb )

inline

Definition at line 39 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