EvtIntervalDecayAmp< T > Class Template Reference

#include <EvtIntervalDecayAmp.hh>

Inheritance diagram for EvtIntervalDecayAmp< T >:

Public Member Functions
	EvtIntervalDecayAmp ()
	EvtIntervalDecayAmp (const EvtIntervalDecayAmp< T > &other)
virtual	~EvtIntervalDecayAmp ()
virtual void	init ()
virtual void	initProbMax ()
virtual void	decay (EvtParticle *p)
virtual EvtAmpFactory< T > *	createFactory (const EvtMultiChannelParser &parser)=0
virtual std::vector< EvtVector4R >	initDaughters (const T &p) const =0
const T &	x () const
EvtComplex	amplNonCP (const T &x)
EvtPdfSum< T > *	getPC ()
Public Member Functions inherited from EvtDecayAmp
void	makeDecay (EvtParticle *p)
void	setWeight (double weight)
void	vertex (const EvtComplex &amp)
void	vertex (int i1, const EvtComplex &amp)
void	vertex (int i1, int i2, const EvtComplex &amp)
void	vertex (int i1, int i2, int i3, const EvtComplex &amp)
void	vertex (int *i1, const EvtComplex &amp)
virtual	~EvtDecayAmp ()
Public Member Functions inherited from EvtDecayBase
virtual void	getName (std::string &name)=0
virtual EvtDecayBase *	clone ()=0
virtual std::string	commandName ()
virtual void	command (std::string cmd)
double	getProbMax (double prob)
double	resetProbMax (double prob)
	EvtDecayBase ()
virtual	~EvtDecayBase ()
virtual bool	matchingDecay (const EvtDecayBase &other) const
EvtId	getParentId ()
double	getBranchingFraction ()
void	disableCheckQ ()
void	checkQ ()
int	getNDaug ()
EvtId *	getDaugs ()
EvtId	getDaug (int i)
int	getNArg ()
int	getPHOTOS ()
void	setPHOTOS ()
void	setVerbose ()
void	setSummary ()
double *	getArgs ()
std::string *	getArgsStr ()
double	getArg (int j)
std::string	getArgStr (int j)
std::string	getModelName ()
int	getDSum ()
int	summary ()
int	verbose ()
void	saveDecayInfo (EvtId ipar, int ndaug, EvtId *daug, int narg, std::vector< std::string > &args, std::string name, double brfr)
void	printSummary ()
void	setProbMax (double prbmx)
void	noProbMax ()
void	checkNArg (int a1, int a2=-1, int a3=-1, int a4=-1)
void	checkNDaug (int d1, int d2=-1)
void	checkSpinParent (EvtSpinType::spintype sp)
void	checkSpinDaughter (int d1, EvtSpinType::spintype sp)
virtual int	nRealDaughters ()

Protected Attributes
double	_probMax
int	_nScan
T	_x
EvtAmpFactory< T > *	_fact
Protected Attributes inherited from EvtDecayAmp
EvtAmp	_amp2
Protected Attributes inherited from EvtDecayBase
bool	_daugsDecayedByParentModel

Additional Inherited Members
Static Public Member Functions inherited from EvtDecayBase
static void	findMasses (EvtParticle *p, int ndaugs, EvtId daugs[10], double masses[10])
static void	findMass (EvtParticle *p)
static double	findMaxMass (EvtParticle *p)
Protected Member Functions inherited from EvtDecayBase
bool	daugsDecayedByParentModel ()

Detailed Description

template<class T>
class EvtIntervalDecayAmp< T >

Definition at line 38 of file EvtIntervalDecayAmp.hh.

Constructor & Destructor Documentation

EvtIntervalDecayAmp() [1/2]

template<class T>

EvtIntervalDecayAmp< T >::EvtIntervalDecayAmp ( )

inline

Definition at line 41 of file EvtIntervalDecayAmp.hh.

: _probMax( 0. ), _nScan( 0 ), _fact( 0 ) {}

Referenced by EvtIntervalDecayAmp().

EvtIntervalDecayAmp() [2/2]

template<class T>

EvtIntervalDecayAmp< T >::EvtIntervalDecayAmp ( const EvtIntervalDecayAmp< T > & other )

inline

Definition at line 43 of file EvtIntervalDecayAmp.hh.

      : _probMax( other._probMax ), _nScan( other._nScan ), COPY_PTR( _fact ) {}

~EvtIntervalDecayAmp()

template<class T>

virtual EvtIntervalDecayAmp< T >::~EvtIntervalDecayAmp ( )

inlinevirtual

Definition at line 46 of file EvtIntervalDecayAmp.hh.

{ delete _fact; }

Member Function Documentation

amplNonCP()

template<class T>

EvtComplex EvtIntervalDecayAmp< T >::amplNonCP ( const T & x )

inline

Definition at line 157 of file EvtIntervalDecayAmp.hh.

{ return _fact->getAmp()->evaluate( x ); }

Referenced by EvtBtoKD3P::decay(), and decay().

createFactory()

template<class T>

virtual EvtAmpFactory< T > * EvtIntervalDecayAmp< T >::createFactory ( const EvtMultiChannelParser & parser )

pure virtual

Implemented in EvtPto3P.

Referenced by init().

decay()

template<class T>

virtual void EvtIntervalDecayAmp< T >::decay ( EvtParticle * p )

inlinevirtual

Implements EvtDecayBase.

Definition at line 96 of file EvtIntervalDecayAmp.hh.

                                       {
    // Set things up in most general way
 
    static EvtId B0  = EvtPDL::getId( "B0" );
    static EvtId B0B = EvtPDL::getId( "anti-B0" );
    double       t;
    EvtId        other_b;
    EvtComplex   ampl( 0., 0. );
 
    // Sample using pole-compensator pdf
 
    EvtPdfSum<T>* pc = getPC();
    _x               = pc->randomPoint();
 
    if ( _fact->isCPModel() )
    {
 
      EvtCPUtil::OtherB( p, t, other_b );
      EvtComplex A    = _fact->getAmp()->evaluate( _x );
      EvtComplex Abar = _fact->getAmpConj()->evaluate( _x );
      double     dm   = _fact->dm();
 
      if ( other_b == B0B )
        ampl = A * cos( dm * t / ( 2 * EvtConst::c ) ) +
               Abar * sin( dm * t / ( 2 * EvtConst::c ) );
      if ( other_b == B0 )
        ampl = A * cos( dm * t / ( 2 * EvtConst::c ) ) -
               Abar * sin( dm * t / ( 2 * EvtConst::c ) );
    }
    else { ampl = amplNonCP( _x ); }
 
    // Pole-compensate
 
    double comp = sqrt( pc->evaluate( _x ) );
    assert( comp > 0 );
    vertex( ampl / comp );
 
    // Now generate random angles, rotate and setup
    // the daughters
 
    std::vector<EvtVector4R> v = initDaughters( _x );
 
    int N = p->getNDaug();
    if ( v.size() != N )
    {
 
      report( INFO, "EvtGen" ) << "Number of daughters " << N << std::endl;
      report( INFO, "EvtGen" ) << "Momentum vector size " << v.size() << std::endl;
      assert( 0 );
    }
 
    int i;
    for ( i = 0; i < N; i++ ) { p->getDaug( i )->init( getDaugs()[i], v[i] ); }
  }

getPC()

template<class T>

EvtPdfSum< T > * EvtIntervalDecayAmp< T >::getPC ( )

inline

Definition at line 158 of file EvtIntervalDecayAmp.hh.

{ return _fact->getPC(); }

Referenced by EvtBtoKD3P::decay(), and decay().

init()

template<class T>

virtual void EvtIntervalDecayAmp< T >::init ( )

inlinevirtual

Reimplemented from EvtDecayBase.

Definition at line 50 of file EvtIntervalDecayAmp.hh.

                      {
    // Collect model parameters and parse them
 
    vector<std::string> args;
    int                 i;
    for ( i = 0; i < getNArg(); i++ ) args.push_back( getArgStr( i ) );
    EvtMultiChannelParser parser;
    parser.parse( args );
 
    // Create factory and interval
 
    if ( VERBOSE ) report( INFO, "EvtGen" ) << "Create factory and interval" << std::endl;
    _fact = createFactory( parser );
 
    // Maximum PDF value over the Dalitz plot can be specified, or a scan
    // can be performed.
 
    _probMax = parser.pdfMax();
    _nScan   = parser.nScan();
    if ( VERBOSE ) report( INFO, "EvtGen" ) << "Pdf maximum " << _probMax << std::endl;
    if ( VERBOSE ) report( INFO, "EvtGen" ) << "Scan number " << _nScan << std::endl;
  }

initDaughters()

template<class T>

virtual std::vector< EvtVector4R > EvtIntervalDecayAmp< T >::initDaughters ( const T & p ) const

pure virtual

Implemented in EvtPto3P.

Referenced by decay().

initProbMax()

template<class T>

virtual void EvtIntervalDecayAmp< T >::initProbMax ( )

inlinevirtual

Reimplemented from EvtDecayBase.

Definition at line 73 of file EvtIntervalDecayAmp.hh.

                             {
    if ( 0 == _nScan )
    {
 
      if ( _probMax > 0 ) setProbMax( _probMax );
      else assert( 0 );
    }
    else
    {
 
      double        factor = 1.2; // increase maximum probability by 20%
      EvtAmpPdf<T>  pdf( *_fact->getAmp() );
      EvtPdfSum<T>* pc = _fact->getPC();
      EvtPdfDiv<T>  pdfdiv( pdf, *pc );
      printf( "Sampling %d points to find maximum\n", _nScan );
      EvtPdfMax<T> x = pdfdiv.findMax( *pc, _nScan );
      _probMax       = factor * x.value();
      printf( "Found maximum %f\n", x.value() );
      printf( "Increase to   %f\n", _probMax );
      setProbMax( _probMax );
    }
  }

x()

template<class T>

const T & EvtIntervalDecayAmp< T >::x ( ) const

inline

Definition at line 156 of file EvtIntervalDecayAmp.hh.

{ return _x; }

Referenced by amplNonCP(), and initProbMax().

Member Data Documentation

_fact

template<class T>

EvtAmpFactory<T>* EvtIntervalDecayAmp< T >::_fact

protected

Definition at line 165 of file EvtIntervalDecayAmp.hh.

Referenced by amplNonCP(), decay(), EvtIntervalDecayAmp(), EvtIntervalDecayAmp(), getPC(), init(), initProbMax(), and ~EvtIntervalDecayAmp().

_nScan

template<class T>

int EvtIntervalDecayAmp< T >::_nScan

protected

Definition at line 162 of file EvtIntervalDecayAmp.hh.

Referenced by EvtIntervalDecayAmp(), EvtIntervalDecayAmp(), init(), and initProbMax().

_probMax

template<class T>

double EvtIntervalDecayAmp< T >::_probMax

protected

Definition at line 161 of file EvtIntervalDecayAmp.hh.

Referenced by EvtIntervalDecayAmp(), EvtIntervalDecayAmp(), init(), and initProbMax().

_x

template<class T>

T EvtIntervalDecayAmp< T >::_x

protected

Definition at line 163 of file EvtIntervalDecayAmp.hh.

Referenced by decay(), and x().

---

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

• EvtIntervalDecayAmp.hh