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 39 of file EvtIntervalDecayAmp.hh.

Constructor & Destructor Documentation

EvtIntervalDecayAmp() [1/2]

template<class T >

EvtIntervalDecayAmp< T >::EvtIntervalDecayAmp ( )

inline

Definition at line 43 of file EvtIntervalDecayAmp.hh.

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

EvtIntervalDecayAmp() [2/2]

template<class T >

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

inline

Definition at line 47 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 52 of file EvtIntervalDecayAmp.hh.

  {
    delete _fact;
  }

Member Function Documentation

amplNonCP()

template<class T >

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

inline

Definition at line 170 of file EvtIntervalDecayAmp.hh.

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

Referenced by EvtBtoKD3P::decay(), and EvtIntervalDecayAmp< T >::decay().

createFactory()

template<class T >

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

pure virtual

Implemented in EvtPto3P.

Referenced by EvtIntervalDecayAmp< T >::init().

decay()

template<class T >

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

inlinevirtual

Implements EvtDecayBase.

Definition at line 108 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 171 of file EvtIntervalDecayAmp.hh.

{return _fact->getPC();}

Referenced by EvtBtoKD3P::decay(), and EvtIntervalDecayAmp< T >::decay().

init()

template<class T >

virtual void EvtIntervalDecayAmp< T >::init ( )

inlinevirtual

Reimplemented from EvtDecayBase.

Definition at line 60 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 EvtIntervalDecayAmp< T >::decay().

initProbMax()

template<class T >

virtual void EvtIntervalDecayAmp< T >::initProbMax ( )

inlinevirtual

Reimplemented from EvtDecayBase.

Definition at line 86 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 169 of file EvtIntervalDecayAmp.hh.

{return _x;}

Referenced by EvtIntervalDecayAmp< T >::amplNonCP(), and EvtIntervalDecayAmp< T >::initProbMax().

Member Data Documentation

_fact

template<class T >

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

protected

Definition at line 178 of file EvtIntervalDecayAmp.hh.

Referenced by EvtIntervalDecayAmp< T >::amplNonCP(), EvtIntervalDecayAmp< T >::decay(), EvtIntervalDecayAmp< T >::getPC(), EvtIntervalDecayAmp< T >::init(), EvtIntervalDecayAmp< T >::initProbMax(), and EvtIntervalDecayAmp< T >::~EvtIntervalDecayAmp().

_nScan

template<class T >

int EvtIntervalDecayAmp< T >::_nScan

protected

Definition at line 175 of file EvtIntervalDecayAmp.hh.

Referenced by EvtIntervalDecayAmp< T >::init(), and EvtIntervalDecayAmp< T >::initProbMax().

_probMax

template<class T >

double EvtIntervalDecayAmp< T >::_probMax

protected

Definition at line 174 of file EvtIntervalDecayAmp.hh.

Referenced by EvtIntervalDecayAmp< T >::init(), and EvtIntervalDecayAmp< T >::initProbMax().

_x

template<class T >

T EvtIntervalDecayAmp< T >::_x

protected

Definition at line 176 of file EvtIntervalDecayAmp.hh.

Referenced by EvtIntervalDecayAmp< T >::decay(), and EvtIntervalDecayAmp< T >::x().

---

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

• EvtIntervalDecayAmp.hh