EvtRhoPi Class Reference

#include <EvtRhoPi.hh>

Inheritance diagram for EvtRhoPi:

Public Member Functions
	EvtRhoPi ()
virtual	~EvtRhoPi ()
void	getName (std::string &name)
EvtDecayBase *	clone ()
void	initProbMax ()
void	init ()
void	decay (EvtParticle *p)
double	AmplitudeSquare (double ResonanceMass, double ResonanceWidth, double r1, double r2, double phase1, double phase2)
EvtVector4R	GetDaugMomLab (int i)
EvtVector4R	GetDaugMomCM (int i)
EvtVector4R	GetDaugMomHel (int i)
int	GetNdaug ()
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 EvtRhoPi.hh.

Constructor & Destructor Documentation

EvtRhoPi()

EvtRhoPi::EvtRhoPi ( )

inline

Definition at line 32 of file EvtRhoPi.hh.

{}

Referenced by clone().

~EvtRhoPi()

EvtRhoPi::~EvtRhoPi ( )

virtual

Definition at line 45 of file EvtRhoPi.cc.

{}

Member Function Documentation

AmplitudeSquare()

double EvtRhoPi::AmplitudeSquare	(	double	ResonanceMass,
		double	ResonanceWidth,
		double	r1,
		double	r2,
		double	phase1,
		double	phase2 )

Definition at line 128 of file EvtRhoPi.cc.

                                                                            {
  EvtVector4R dp1 = GetDaugMomLab( 0 );
  EvtVector4R dp2 = GetDaugMomLab( 1 );
  EvtVector4R dp3 = GetDaugMomLab( 2 );
  VVS Jpsi2rhopi( dp1, dp2, dp3, ResonanceMass, ResonanceWidth, r1, r2, phase1, phase2 );
  return Jpsi2rhopi.amps();
}

Referenced by decay().

clone()

EvtDecayBase * EvtRhoPi::clone ( )

virtual

Implements EvtDecayBase.

Definition at line 49 of file EvtRhoPi.cc.

{ return new EvtRhoPi; }

decay()

void EvtRhoPi::decay ( EvtParticle * p )

virtual

default use rho(770) as the intermediate state, this set parameter will generate the symmetric momentum distribution for pi+ pi-.

Implements EvtDecayBase.

Definition at line 65 of file EvtRhoPi.cc.

                                     {
 
  ///// default use rho(770) as the intermediate state, this set parameter will generate the
  /// symmetric momentum distribution for pi+ pi-.
  double ResonanceMass  = 0.7710;
  double ResonanceWidth = 0.1492;
  double r1             = 0.9;
  double phase1         = 0;
  double r2             = 1.1;
  double phase2         = 1.5;
  if ( getNArg() > 0 )
  {
    ResonanceMass  = getArg( 0 );
    ResonanceWidth = getArg( 1 );
    r1             = getArg( 2 );
    phase1         = getArg( 3 );
    r2             = getArg( 4 );
    phase2         = getArg( 5 );
  }
  double amps, SamAmps, rd1;
  // calculated the max amplitude square in 20000 events, is it enough larger?
  if ( nrun == 1 )
  {
    int ir, nd;
    for ( ir = 0; ir <= 20000; ir++ )
    {
      p->initializePhaseSpace( getNDaug(), getDaugs() );
      int nd = p->getNDaug(), i;
      _nd    = nd;
      for ( i = 0; i <= nd - 1; i++ )
      {
        _p4Lab[i] = p->getDaug( i )->getP4Lab();
        _p4CM[i]  = p->getDaug( i )->getP4();
      }
      amps = AmplitudeSquare( ResonanceMass, ResonanceWidth, r1, r2, phase1, phase2 );
      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"
                              << endl;
    abort();
  }
  // cout<<"max_amp="<<max_amps<<endl;
 
loop:
  p->initializePhaseSpace( getNDaug(), getDaugs() );
  int i;
  for ( i = 0; i <= p->getNDaug() - 1; i++ )
  {
    _p4Lab[i] = p->getDaug( i )->getP4Lab();
    _p4CM[i]  = p->getDaug( i )->getP4();
  }
  //  Put phase space results into the daughters.
  amps    = AmplitudeSquare( ResonanceMass, ResonanceWidth, r1, r2, phase1, phase2 );
  SamAmps = amps / max_amps;
  rd1     = EvtRandom::Flat( 0.0, 1.0 );
  if ( rd1 >= SamAmps ) goto loop;
  return;
}

GetDaugMomCM()

EvtVector4R EvtRhoPi::GetDaugMomCM ( int i )

inline

Definition at line 49 of file EvtRhoPi.hh.

{ return _p4CM[i]; }

GetDaugMomHel()

EvtVector4R EvtRhoPi::GetDaugMomHel ( int i )

inline

Definition at line 50 of file EvtRhoPi.hh.

{ return _p4Hel[i]; }

GetDaugMomLab()

EvtVector4R EvtRhoPi::GetDaugMomLab ( int i )

inline

Definition at line 48 of file EvtRhoPi.hh.

{ return _p4Lab[i]; }

Referenced by AmplitudeSquare().

getName()

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

virtual

Implements EvtDecayBase.

Definition at line 47 of file EvtRhoPi.cc.

{ model_name = "RhoPi"; }

GetNdaug()

int EvtRhoPi::GetNdaug ( )

inline

Definition at line 52 of file EvtRhoPi.hh.

{ return _nd; }

init()

void EvtRhoPi::init ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 51 of file EvtRhoPi.cc.

                    {
 
  //  checkNArg(6);
  EvtSpinType::spintype parenttype = EvtPDL::getSpinType( getParentId() );
  //  if ( parenttype == EvtSpinType::SCALAR ){
  //  report(ERROR,"EvtGen")<<"Scalar decays with flat distribution"<<endl;
  //  ::abort();
  // }
}

initProbMax()

void EvtRhoPi::initProbMax ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 60 of file EvtRhoPi.cc.

{ noProbMax(); }

---

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

• EvtRhoPi.hh
• EvtRhoPi.cc