CGEM BOSS 6.6.5.i
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtJscont Class Reference

#include <EvtJscont.hh>

+ Inheritance diagram for EvtJscont:

Public Member Functions

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

Detailed Description

Definition at line 28 of file EvtJscont.hh.

Constructor & Destructor Documentation

◆ EvtJscont()

EvtJscont::EvtJscont ( )
inline

Definition at line 32 of file EvtJscont.hh.

32{}

Referenced by clone().

◆ ~EvtJscont()

EvtJscont::~EvtJscont ( )
virtual

Definition at line 46 of file EvtJscont.cc.

46{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtJscont::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 54 of file EvtJscont.cc.

54 {
55
56 return new EvtJscont;
57
58}

◆ decay()

void EvtJscont::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 76 of file EvtJscont.cc.

76 {
77
79 static int first=1;
80
81 if (first){
82 first=0;
83
84 float val=0.6;
85 if ( getNArg()>1) {
86 val=getArg(1);
87 }
88 char vak[20];
89 sprintf(vak,"PARJ(13)=%f",val);
90 std::string temp(vak);
91 lugive_(temp.c_str(),strlen(temp.c_str()));
92 }
93 EvtVector4R p4[100];
94
95 double energy=p->mass();
96
97 int flavor;
98
99 int i,more;
100 int ndaugjs;
101 int kf[100];
102 EvtId id[100];
103 int type[MAX_DAUG];
104
105 flavor=(int)getArg(0);
106
107 double px[100],py[100],pz[100],e[100];
108
109 if ( p->getNDaug() != 0 ) { return;}
110 do{
111
112 continuum_(&energy,&flavor,&ndaugjs,kf,px,py,pz,e);
113
114 for(i=0;i<ndaugjs;i++){
115
116 id[i]=EvtPDL::evtIdFromStdHep(kf[i]);
117
118 type[i]=EvtPDL::getSpinType(id[i]);
119
120 // have to protect against negative mass^2 for massless particles
121 // i.e. neutrinos and photons.
122 // this is uggly but I need to fix it right now....
123
124 if (px[i]*px[i]+py[i]*py[i]+pz[i]*pz[i]>=e[i]*e[i]){
125
126 e[i]=sqrt(px[i]*px[i]+py[i]*py[i]+pz[i]*pz[i])+0.0000000000001;
127
128 }
129
130 p4[i].set(e[i],px[i],py[i],pz[i]);
131
132 }
133
134 int channel=EvtDecayTable::inChannelList(p->getId(),ndaugjs,id);
135
136 more=((channel!=-1)&&(channel!=p->getChannel()));
137
138
139 }while(more);
140
141 p->makeDaughters(ndaugjs,id);
142
143 for(i=0;i<ndaugjs;i++){
144 p->getDaug(i)->init( id[i], p4[i] );
145 }
146 return ;
147}
void lugive_(const char *cnfgstr, int length)
void continuum_(double *, int *, int *, int *, double *, double *, double *, double *)
const int MAX_DAUG
************Class m_ypar INTEGER m_KeyWgt INTEGER m_nphot INTEGER m_KeyGPS INTEGER m_IsBeamPolarized INTEGER m_EvtGenInterface DOUBLE PRECISION m_Emin DOUBLE PRECISION m_sphot DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_q2 DOUBLE PRECISION m_PolBeam2 DOUBLE PRECISION m_xErrPb *COMMON c_KK2f $ !CMS energy average $ !Spin Polarization vector first beam $ !Spin Polarization vector second beam $ !Beam energy spread[GeV] $ !minimum hadronization energy[GeV] $ !input READ never touch them !$ !debug facility $ !maximum weight $ !inverse alfaQED $ !minimum real photon energy
Definition KK2f.h:50
double getArg(int j)
static int inChannelList(EvtId parent, int ndaug, EvtId *daugs)
Definition EvtId.hh:27
static void jetSetInit()
Definition EvtJetSet.cc:785
static EvtId evtIdFromStdHep(int stdhep)
Definition EvtPDL.cc:244
static EvtSpinType::spintype getSpinType(EvtId i)
Definition EvtPDL.hh:61
void makeDaughters(int ndaug, EvtId *id)
virtual void init(EvtId part_n, const EvtVector4R &p4)=0
EvtId getId() const
int getNDaug() const
EvtParticle * getDaug(int i)
double mass() const
int getChannel() const
void set(int i, double d)
Index first(Pair i)

◆ getName()

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

Implements EvtDecayBase.

Definition at line 48 of file EvtJscont.cc.

48 {
49
50 model_name="JSCONT";
51
52}

◆ init()

void EvtJscont::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 60 of file EvtJscont.cc.

60 {
61
62 // check that there are 1 argument
63
64 checkNArg(1,2);
65
66}
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)

◆ initProbMax()

void EvtJscont::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 69 of file EvtJscont.cc.

69 {
70
71 noProbMax();
72
73}

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