EvtFDC Class Reference

#include <EvtFDC.hh>

Inheritance diagram for EvtFDC:

Public Member Functions
	EvtFDC ()
virtual	~EvtFDC ()
void	getName (std::string &name)
EvtDecayBase *	clone ()
void	initProbMax ()
void	init ()
void	decay (EvtParticle *p)
double	AmplitudeSquare ()
int	GetNdaug ()
void	setEvtMomentum (EvtVector4R *p4)
Public Member Functions inherited from EvtDecayIncoherent
void	makeDecay (EvtParticle *p)
virtual	~EvtDecayIncoherent ()
void	setDaughterSpinDensity (int daughter)
int	isDaughterSpinDensitySet (int daughter)
Public Member Functions inherited from EvtDecayBase
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 ()

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 ()
Protected Attributes inherited from EvtDecayBase
bool	_daugsDecayedByParentModel

Detailed Description

Definition at line 29 of file EvtFDC.hh.

Constructor & Destructor Documentation

EvtFDC()

EvtFDC::EvtFDC ( )

inline

Definition at line 32 of file EvtFDC.hh.

{}

Referenced by clone().

~EvtFDC()

EvtFDC::~EvtFDC ( )

virtual

Definition at line 61 of file EvtFDC.cc.

{}

Member Function Documentation

AmplitudeSquare()

double EvtFDC::AmplitudeSquare

(

)

clone()

EvtDecayBase * EvtFDC::clone ( )

virtual

Implements EvtDecayBase.

Definition at line 65 of file EvtFDC.cc.

{ return new EvtFDC; }

decay()

void EvtFDC::decay ( EvtParticle * p )

virtual

Implements EvtDecayBase.

Definition at line 92 of file EvtFDC.cc.

                                   {
  double amps, SamAmps, rd1;
  // calculated the max amplitude square in 20000 events, is it enough larger?
  if ( nrun == 1 )
  {
    int ir, nd;
    std::cout << "Wait for a moment for calculation of maxiumal amplitude squared"
              << std::endl;
    for ( ir = 0; ir <= 200000; ir++ )
    {
    loop0:
      p->initializePhaseSpace( getNDaug(), getDaugs() );
      int nd = p->getNDaug();
      _nd    = nd;
      for ( int i = 0; i <= nd - 1; i++ ) { _p4[i] = p->getDaug( i )->getP4Lab(); }
      for ( int i = 0; i < _nd; i++ )
      {
        e[i]  = _p4[i].get( 0 );
        px[i] = _p4[i].get( 1 );
        py[i] = _p4[i].get( 2 );
        pz[i] = _p4[i].get( 3 );
      }
      // for debuging
      amps = fdcevtgen_( e, px, py, pz );
      // std::cout<<"nrun "<<nrun<<std::endl;
      if ( amps < 0 ) goto loop0;
      if ( amps > max_amps ) max_amps = amps * 1.01;
      nrun++;
    }
  }
  if ( max_amps == 0.0 )
  {
    report( ERROR, "EvtGen" )
        << "The decay amplitude square should be positive number, but it is " << max_amps
        << endl;
    abort();
  }
  // cout<<"max_amp="<<max_amps<<endl;
 
loop:
  p->initializePhaseSpace( getNDaug(), getDaugs() );
  int i;
  for ( i = 0; i <= p->getNDaug() - 1; i++ ) { _p4[i] = p->getDaug( i )->getP4Lab(); }
  setEvtMomentum( _p4 );
  //  Put phase space results into the daughters.
  amps = fdcevtgen_( e, px, py, pz );
  if ( amps < 0 ) goto loop;
  SamAmps = amps / max_amps;
  rd1     = EvtRandom::Flat( 0.0, 1.0 );
  if ( rd1 >= SamAmps ) goto loop;
  if ( amps > max_amps ) max_amps = amps * 1.01;
  nrun++;
}

getName()

void EvtFDC::getName ( std::string & name )

virtual

Implements EvtDecayBase.

Definition at line 63 of file EvtFDC.cc.

{ model_name = "FDC"; }

GetNdaug()

int EvtFDC::GetNdaug ( )

inline

Definition at line 47 of file EvtFDC.hh.

{ return _nd; }

init()

void EvtFDC::init ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 67 of file EvtFDC.cc.

                  {
  init_fit_();
  initevtgen_();
  int ndiag = -1; // Feyman diagram number, started from 1, if -1, generate all intermediate
                  // states including interference effects
  if ( getNArg() == 1 ) { ndiag = getArg( 0 ); }
  setamp_( &ndiag );
 
  EvtSpinType::spintype parenttype = EvtPDL::getSpinType( getParentId() );
}

initProbMax()

void EvtFDC::initProbMax ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 77 of file EvtFDC.cc.

{ noProbMax(); }

setEvtMomentum()

void EvtFDC::setEvtMomentum

(

EvtVector4R *

p4

)

Definition at line 82 of file EvtFDC.cc.

                                             {
  for ( int i = 0; i < _nd; i++ )
  {
    e[i]  = _p4[i].get( 0 );
    px[i] = _p4[i].get( 1 );
    py[i] = _p4[i].get( 2 );
    pz[i] = _p4[i].get( 3 );
  }
}

Referenced by decay().

---

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

• EvtFDC.hh
• EvtFDC.cc