BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtJ2BB1 Class Reference

#include <EvtJ2BB1.hh>

Inheritance diagram for EvtJ2BB1:

Public Member Functions

 EvtJ2BB1 ()
virtual ~EvtJ2BB1 ()
void getName (std::string &name)
EvtDecayBaseclone ()
void decay (EvtParticle *p)
void init ()
void init (int, double *)
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 initProbMax ()
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 ()
EvtIdgetDaugs ()
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 ()

Public Attributes

double alpha
double u

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 EvtDecayAmp
EvtAmp _amp2
Protected Attributes inherited from EvtDecayBase
bool _daugsDecayedByParentModel

Detailed Description

Definition at line 28 of file EvtJ2BB1.hh.

Constructor & Destructor Documentation

◆ EvtJ2BB1()

EvtJ2BB1::EvtJ2BB1 ( )
inline

Definition at line 31 of file EvtJ2BB1.hh.

31{}

Referenced by clone().

◆ ~EvtJ2BB1()

EvtJ2BB1::~EvtJ2BB1 ( )
virtual

Definition at line 36 of file EvtJ2BB1.cc.

36{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtJ2BB1::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 40 of file EvtJ2BB1.cc.

40{ return new EvtJ2BB1; }

◆ decay()

void EvtJ2BB1::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 56 of file EvtJ2BB1.cc.

56 {
57
59
60 EvtParticle *v, *s1;
61 EvtVector4R pv, ps, ppr;
62
63 v = p->getDaug( 0 ); // Baryon
64 s1 = p->getDaug( 1 ); // Baryon
65 pv = v->getP4();
66 ps = s1->getP4();
67 ppr = p->getP4();
68 // Put phase space results into the daughters.
69
70 EvtHelSys angles( ppr, pv ); // using helicity sys.angles
71 double theta = angles.getHelAng( 1 );
72 double phi = angles.getHelAng( 2 );
73 double gamma = 0;
74 double costheta = cos( theta ); // using helicity angles in parent system
75 // double costheta=ppr.get(3)/ppr.d3mag(); //using Lab system
76
77 if ( getNArg() > 0 ) { alpha = getArg( 0 ); }
78 else
79 {
80 double mass_b1 = EvtPDL::getMass( getDaug( 0 ) );
81 double mass_b2 = EvtPDL::getMass( getDaug( 1 ) );
82 double mass_M = EvtPDL::getMass( getParentId() );
83 double c1 = pow( pv.d3mag(), 2.0 );
84 double R = ( mass_b1 + pv.get( 0 ) ) * ( mass_b2 + ps.get( 0 ) );
85 u = pow( ( R - c1 ) / ( R + c1 ), 2.0 );
86 alpha = ( pow( 4 + u, 2.0 ) - u * pow( 2 + 3 * u, 2.0 ) ) /
87 ( pow( 4 + u, 2.0 ) + u * pow( 2 + 3 * u, 2.0 ) );
88 }
89
90 double F01 = sqrt( ( 1 + alpha ) / 2 );
91 double F00 = sqrt( ( 1 - alpha ) / 4 );
92 double F10 = F01;
93 double F11 = F00;
94
95 vertex( 0, 0, 0, Djmn( 1, 1, 0, phi, theta, gamma ) * F00 ); // J/psi helicity =1
96 // corresponding index=0
97 vertex( 0, 0, 1, Djmn( 1, 1, 1, phi, theta, gamma ) * F01 );
98 vertex( 0, 1, 0, Djmn( 1, 1, -1, phi, theta, gamma ) * F10 );
99 vertex( 0, 1, 1, Djmn( 1, 1, 0, phi, theta, gamma ) * F11 );
100 vertex( 1, 0, 0, Djmn( 1, -1, 0, phi, theta, gamma ) * F00 ); // J/psi helicity =-1
101 // corresponding index=1
102 vertex( 1, 0, 1, Djmn( 1, -1, 1, phi, theta, gamma ) * F01 );
103 vertex( 1, 1, 0, Djmn( 1, -1, -1, phi, theta, gamma ) * F10 );
104 vertex( 1, 1, 1, Djmn( 1, -1, 0, phi, theta, gamma ) * F11 );
105 vertex( 2, 0, 0, Djmn( 1, 0, 0, phi, theta, gamma ) * F00 ); // J/psi helicity =0
106 // corresponding index=2
107 vertex( 2, 0, 1, Djmn( 1, 0, 1, phi, theta, gamma ) * F01 );
108 vertex( 2, 1, 0, Djmn( 1, 0, -1, phi, theta, gamma ) * F10 );
109 vertex( 2, 1, 1, Djmn( 1, 0, 0, phi, theta, gamma ) * F11 );
110
111 return;
112}
EvtComplex Djmn(int j, int m, int n, double phi, double theta, double gamma)
Definition EvtHelSys.cc:165
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_alfinv DOUBLE PRECISION m_Xenph INTEGER m_KeyWtm INTEGER m_idyfs DOUBLE PRECISION m_zini DOUBLE PRECISION m_q2 DOUBLE PRECISION m_Wt_KF DOUBLE PRECISION m_WtCut INTEGER m_KFfin *COMMON c_KarLud $ !Input CMS energy[GeV] $ !CMS energy after beam spread beam strahlung[GeV] $ !Beam energy spread[GeV] $ !z boost due to beam spread $ !electron beam mass *ff pair spectrum $ !minimum v
Definition KarLud.h:35
void vertex(const EvtComplex &amp)
double getArg(int j)
EvtId getParentId()
EvtId * getDaugs()
EvtId getDaug(int i)
double u
Definition EvtJ2BB1.hh:43
double alpha
Definition EvtJ2BB1.hh:42
static double getMass(EvtId i)
Definition EvtPDL.hh:44
const EvtVector4R & getP4() const
EvtParticle * getDaug(int i)
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
double get(int i) const
double d3mag() const
complex_t R(double Q2, double M2, double G, double Mp2, double Mm2)
Definition TUtil.h:22

◆ getName()

void EvtJ2BB1::getName ( std::string & name)
virtual

Implements EvtDecayBase.

Definition at line 38 of file EvtJ2BB1.cc.

38{ model_name = "J2BB1"; }

◆ init() [1/2]

void EvtJ2BB1::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 49 of file EvtJ2BB1.cc.

49 {
50 checkNDaug( 2 );
54}
void checkSpinDaughter(int d1, EvtSpinType::spintype sp)
void checkSpinParent(EvtSpinType::spintype sp)
void checkNDaug(int d1, int d2=-1)

◆ init() [2/2]

void EvtJ2BB1::init ( int ,
double *  )

Member Data Documentation

◆ alpha

double EvtJ2BB1::alpha

Definition at line 42 of file EvtJ2BB1.hh.

Referenced by decay().

◆ u

double EvtJ2BB1::u

Definition at line 43 of file EvtJ2BB1.hh.

Referenced by decay().


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