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

#include <EvtP2GC1.hh>

Inheritance diagram for EvtP2GC1:

Public Member Functions

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

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 EvtP2GC1.hh.

Constructor & Destructor Documentation

◆ EvtP2GC1()

EvtP2GC1::EvtP2GC1 ( )
inline

Definition at line 31 of file EvtP2GC1.hh.

31{}

Referenced by clone().

◆ ~EvtP2GC1()

EvtP2GC1::~EvtP2GC1 ( )
virtual

Definition at line 36 of file EvtP2GC1.cc.

36{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtP2GC1::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 40 of file EvtP2GC1.cc.

40{ return new EvtP2GC1; }

◆ decay()

void EvtP2GC1::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 63 of file EvtP2GC1.cc.

63 {
64
66
67 EvtParticle *v, *s1;
68 EvtVector4R pv, ps, ppr;
69
70 v = p->getDaug( 0 ); // gamma
71 s1 = p->getDaug( 1 ); // chi_c1
72 pv = v->getP4();
73 ps = s1->getP4();
74 ppr = p->getP4();
75
76 // Put phase space results into the daughters.
77 EvtHelSys angles( ppr, pv ); // using helicity sys. angles
78 double theta = angles.getHelAng( 1 );
79 double phi = angles.getHelAng( 2 );
80 double gamma = 0;
81 double alpha = -1. / 3.;
82 double A02 = sqrt( ( 1 + alpha ) / 2 );
83 double A00 = sqrt( ( 1 - alpha ) / 4 );
84 double A12 = -A02, A11 = -A00;
85 int j1, j2, j3;
86 for ( j1 = 0; j1 <= 2; )
87 {
88 for ( j2 = 0; j2 <= 1; )
89 {
90 for ( j3 = 0; j3 <= 2; )
91 {
92 vertex( j1, j2, j3, 0.0 );
93 j3++;
94 }
95 j2++;
96 }
97 j1++;
98 }
99
100 vertex( 0, 1, 2, Djmn( 1, 1, -1, phi, theta, gamma ) * A12 ); // helicity for psiprime = 1
101 vertex( 0, 1, 1, Djmn( 1, 1, 0, phi, theta, gamma ) * A11 );
102 vertex( 0, 0, 2, Djmn( 1, 1, 1, phi, theta, gamma ) * A02 );
103 vertex( 0, 0, 0, Djmn( 1, 1, 0, phi, theta, gamma ) * A00 );
104
105 vertex( 1, 1, 2, Djmn( 1, -1, -1, phi, theta, gamma ) * A12 ); // helicity for psiprime = -1
106 vertex( 1, 1, 1, Djmn( 1, -1, 0, phi, theta, gamma ) * A11 );
107 vertex( 1, 0, 2, Djmn( 1, -1, 1, phi, theta, gamma ) * A02 );
108 vertex( 1, 0, 0, Djmn( 1, -1, 0, phi, theta, gamma ) * A00 );
109
110 vertex( 2, 1, 2, Djmn( 1, 0, -1, phi, theta, gamma ) * A12 ); // helicity for psiprime = -1
111 vertex( 2, 1, 1, Djmn( 1, 0, 0, phi, theta, gamma ) * A11 );
112 vertex( 2, 0, 2, Djmn( 1, 0, 1, phi, theta, gamma ) * A02 );
113 vertex( 2, 0, 0, Djmn( 1, 0, 0, phi, theta, gamma ) * A00 );
114
115 return;
116}
EvtComplex Djmn(int j, int m, int n, double phi, double theta, double gamma)
Definition EvtHelSys.cc:165
double alpha
**********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)
EvtId * getDaugs()
const EvtVector4R & getP4() const
EvtParticle * getDaug(int i)
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)

◆ getName()

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

Implements EvtDecayBase.

Definition at line 38 of file EvtP2GC1.cc.

38{ model_name = "P2GC1"; }

◆ init()

void EvtP2GC1::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 42 of file EvtP2GC1.cc.

42 {
43 static EvtId GAM = EvtPDL::getId( "gamma" );
44 // check that there are 0 arguments
45 checkNArg( 0 );
46 checkNDaug( 2 );
50 if ( !( getDaug( 0 ) == GAM ) )
51 {
52 report( ERROR, "EvtGen" ) << "EvtP2GC1 generator radiative decays expected " << endl;
53 ::abort();
54 }
55}
ostream & report(Severity severity, const char *facility)
Definition EvtReport.cc:34
@ ERROR
Definition EvtReport.hh:49
void checkSpinDaughter(int d1, EvtSpinType::spintype sp)
void checkSpinParent(EvtSpinType::spintype sp)
void checkNDaug(int d1, int d2=-1)
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
EvtId getDaug(int i)
static EvtId getId(const std::string &name)
Definition EvtPDL.cc:272

◆ initProbMax()

void EvtP2GC1::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 57 of file EvtP2GC1.cc.

57 {
58
59 // Hard coded... should not be hard to calculate...
60 setProbMax( 1.67 );
61}
void setProbMax(double prbmx)

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