BOSS 7.0.9
BESIII Offline Software System
Loading...
Searching...
No Matches
Public Member Functions | List of all members
EvtPhokhara_KK Class Reference

#include <EvtPhokhara_KK.hh>

+ Inheritance diagram for EvtPhokhara_KK:

Public Member Functions

 EvtPhokhara_KK ()
 
virtual ~EvtPhokhara_KK ()
 
void getName (std::string &name)
 
EvtDecayBase * clone ()
 
void decay (EvtParticle *p)
 
std::string commandName ()
 
void command (std::string cmd)
 
void init ()
 
void init_mode (EvtParticle *p)
 
void init_evt (EvtParticle *p)
 
void initProbMax ()
 
int getTotalEvt ()
 
void PhokharaInit (int dummy)
 
void ExclusiveDecay (EvtParticle *p)
 
 EvtPhokhara_KK ()
 
virtual ~EvtPhokhara_KK ()
 
void getName (std::string &name)
 
EvtDecayBase * clone ()
 
void decay (EvtParticle *p)
 
std::string commandName ()
 
void command (std::string cmd)
 
void init ()
 
void init_mode (EvtParticle *p)
 
void init_evt (EvtParticle *p)
 
void initProbMax ()
 
int getTotalEvt ()
 
void PhokharaInit (int dummy)
 
void ExclusiveDecay (EvtParticle *p)
 
- Public Member Functions inherited from EvtDecayIncoherent
void makeDecay (EvtParticle *p)
 
virtual ~EvtDecayIncoherent ()
 
void setDaughterSpinDensity (int daughter)
 
int isDaughterSpinDensitySet (int daughter)
 
void makeDecay (EvtParticle *p)
 
virtual ~EvtDecayIncoherent ()
 
void setDaughterSpinDensity (int daughter)
 
int isDaughterSpinDensitySet (int daughter)
 
- Public Member Functions inherited from EvtDecayBase
virtual void getName (std::string &name)=0
 
virtual void decay (EvtParticle *p)=0
 
virtual void makeDecay (EvtParticle *p)=0
 
virtual EvtDecayBase * clone ()=0
 
virtual void init ()
 
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 ()
 
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 ()
 
virtual void getName (std::string &name)=0
 
virtual void decay (EvtParticle *p)=0
 
virtual void makeDecay (EvtParticle *p)=0
 
virtual EvtDecayBase * clone ()=0
 
virtual void init ()
 
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 ()
 
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)
 
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 ()
 
bool daugsDecayedByParentModel ()
 
- Protected Attributes inherited from EvtDecayBase
bool _daugsDecayedByParentModel
 

Detailed Description

Definition at line 33 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtPhokhara_KK.hh.

Constructor & Destructor Documentation

◆ EvtPhokhara_KK() [1/2]

EvtPhokhara_KK::EvtPhokhara_KK ( )

Definition at line 67 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtPhokhara_KK.cc.

67{}

Referenced by clone().

◆ ~EvtPhokhara_KK() [1/2]

EvtPhokhara_KK::~EvtPhokhara_KK ( )
virtual

Definition at line 68 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtPhokhara_KK.cc.

68 {
69 int i;
70 //the deletion of commands is really uggly!
71
72 if (nphokharadecays==0) {
73 delete [] commands;
74 commands=0;
75 return;
76 }
77
78 for(i=0;i<nphokharadecays;i++){
79 if (phokharadecays[i]==this){
80 phokharadecays[i]=phokharadecays[nphokharadecays-1];
81 nphokharadecays--;
82 if (nphokharadecays==0) {
83 delete [] commands;
84 commands=0;
85 }
86 return;
87 }
88 }
89
90 report(ERROR,"EvtGen") << "Error in destroying Phokhara model!"<<endl;
91
92}
report
ostream & report(Severity severity, const char *facility)
Definition: bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtReport.cc:36

◆ EvtPhokhara_KK() [2/2]

EvtPhokhara_KK::EvtPhokhara_KK ( )

◆ ~EvtPhokhara_KK() [2/2]

virtual EvtPhokhara_KK::~EvtPhokhara_KK ( )
virtual

Member Function Documentation

◆ clone() [1/2]

EvtDecayBase * EvtPhokhara_KK::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 101 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtPhokhara_KK.cc.

101 {
102
103 return new EvtPhokhara_KK;
104
105}
EvtPhokhara_KK::EvtPhokhara_KK
EvtPhokhara_KK()
Definition: bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtPhokhara_KK.cc:67

◆ clone() [2/2]

EvtDecayBase * EvtPhokhara_KK::clone ( )
virtual

Implements EvtDecayBase.

◆ command() [1/2]

void EvtPhokhara_KK::command ( std::string  cmd)
virtual

Reimplemented from EvtDecayBase.

Definition at line 159 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtPhokhara_KK.cc.

159 {
160
161 if (ncommand==lcommand){
162
163 lcommand=10+2*lcommand;
164
165 std::string* newcommands=new std::string[lcommand];
166
167 int i;
168
169 for(i=0;i<ncommand;i++){
170 newcommands[i]=commands[i];
171 }
172
173 delete [] commands;
174
175 commands=newcommands;
176
177 }
178
179 commands[ncommand]=cmd;
180
181 ncommand++;
182
183}

◆ command() [2/2]

void EvtPhokhara_KK::command ( std::string  cmd)
virtual

Reimplemented from EvtDecayBase.

◆ commandName() [1/2]

std::string EvtPhokhara_KK::commandName ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 153 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtPhokhara_KK.cc.

153 {
154
155 return std::string("PhokharaPar");
156
157}

◆ commandName() [2/2]

std::string EvtPhokhara_KK::commandName ( )
virtual

Reimplemented from EvtDecayBase.

◆ decay() [1/2]

void EvtPhokhara_KK::decay ( EvtParticle *  p)
virtual

Implements EvtDecayBase.

Definition at line 187 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtPhokhara_KK.cc.

187 {
188 EvtId myvpho=EvtPDL::getId("vpho");
189 if(p->getId()!=myvpho) {std::cout<<"Parent particle is required to be vpho for Phokhara model"<<std::endl;abort();}
190 if(nevtgen==0) {init_mode(p);}
191 else{init_evt(p);}
192
193 std::cout<<"PHOKHARA : K+K- mode "<<std::endl;
194
195 int istdheppar=EvtPDL::getStdHep(p->getId());
196 int ntrials = 0;
197 int tr_old[3];
198 tr_old[0] = (int)maxima_.tr[0];
199 tr_old[1] = (int)maxima_.tr[1];
200 tr_old[2] = (int)maxima_.tr[2];
201
202 while( ntrials < 1000000)
203 {
204 ievent++;
205 RANLXDF(Ar_r,1);
206 Ar[1] = Ar_r[0];
207
208 if (Ar[1] <= (maxima_.Mmax[0]/(maxima_.Mmax[0]+maxima_.Mmax[1]+maxima_.Mmax[2]))) {
209 maxima_.count[0] = maxima_.count[0]+1.0;
210 GEN_0PH(2,qqmin,ctes_.Sp,cos3min,cos3max);
211 }else
212 if (Ar[1] <= ( (maxima_.Mmax[0]+maxima_.Mmax[1])/(maxima_.Mmax[0]+maxima_.Mmax[1]+maxima_.Mmax[2]))) {
213 maxima_.count[1] = maxima_.count[1]+1.0;
214 GEN_1PH(2,qqmin,qqmax,cos1min,cos1max,cos3min,cos3max);
215 }
216 else {
217 maxima_.count[2] = maxima_.count[2]+1.0;
218 GEN_2PH(2,qqmin,cos1min,cos1max,cos2min,cos2max,cos3min,cos3max);
219 }
220
221 if( ((int)maxima_.tr[0]+(int)maxima_.tr[1]+(int)maxima_.tr[2]) > (tr_old[0]+tr_old[1]+tr_old[2]) ) // event accepted after cuts
222 {
223 goto storedEvents;
224 }
225 ntrials ++;
226 }
227 std::cout <<"FATAL: Could not satisfy cuts after " << ntrials << "trials. Terminate." <<std::endl;
228 //----
229 storedEvents:
230 int more=0;
231 int numstable=0;
232 int numparton=0;
233 EvtId evtnumstable[100];//
234 EvtVector4R p4[20];
235
236 // except ISR photos, products depending on channel
237 if (flags_.pion == 0) { // mu+ mu-
238 // mu+
239 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-13);
240 p4[numstable].set(ctes_.momenta[0][5],ctes_.momenta[1][5], ctes_.momenta[2][5], ctes_.momenta[3][5]);
241 numstable++;
242 // mu -
243 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(13);
244 p4[numstable].set(ctes_.momenta[0][6],ctes_.momenta[1][6], ctes_.momenta[2][6], ctes_.momenta[3][6]);
245 numstable++;
246 }
247 if (flags_.pion == 1) { // pi+ pi-
248 // pi+
249 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(211);
250 p4[numstable].set(ctes_.momenta[0][5],ctes_.momenta[1][5], ctes_.momenta[2][5], ctes_.momenta[3][5]);
251 numstable++;
252 // pi -
253 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-211);
254 p4[numstable].set(ctes_.momenta[0][6],ctes_.momenta[1][6], ctes_.momenta[2][6], ctes_.momenta[3][6]);
255 numstable++;
256 }
257 if (flags_.pion == 2) { // pi+ pi-2pi0
258 // pi0
259 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(111);
260 p4[numstable].set(ctes_.momenta[0][5],ctes_.momenta[1][5], ctes_.momenta[2][5], ctes_.momenta[3][5]);
261 numstable++;
262 // pi0
263 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(111);
264 p4[numstable].set(ctes_.momenta[0][6],ctes_.momenta[1][6], ctes_.momenta[2][6], ctes_.momenta[3][6]);
265 numstable++;
266 // pi-
267 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-211);
268 p4[numstable].set(ctes_.momenta[0][7],ctes_.momenta[1][7], ctes_.momenta[2][7], ctes_.momenta[3][7]);
269 numstable++;
270 // pi +
271 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(211);
272 p4[numstable].set(ctes_.momenta[0][8],ctes_.momenta[1][8], ctes_.momenta[2][8], ctes_.momenta[3][8]);
273 numstable++;
274 }
275 if (flags_.pion == 3) { // 2(pi+ pi-)
276 // pi+
277 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(211);
278 p4[numstable].set(ctes_.momenta[0][5],ctes_.momenta[1][5], ctes_.momenta[2][5], ctes_.momenta[3][5]);
279 numstable++;
280 // pi-
281 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-211);
282 p4[numstable].set(ctes_.momenta[0][6],ctes_.momenta[1][6], ctes_.momenta[2][6], ctes_.momenta[3][6]);
283 numstable++;
284 // pi+
285 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-211);
286 p4[numstable].set(ctes_.momenta[0][7],ctes_.momenta[1][7], ctes_.momenta[2][7], ctes_.momenta[3][7]);
287 numstable++;
288 // pi -
289 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(211);
290 p4[numstable].set(ctes_.momenta[0][8],ctes_.momenta[1][8], ctes_.momenta[2][8], ctes_.momenta[3][8]);
291 numstable++;
292 }
293 if (flags_.pion == 4) { // ppbar
294 // pbar
295 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-2212);
296 p4[numstable].set(ctes_.momenta[0][5],ctes_.momenta[1][5], ctes_.momenta[2][5], ctes_.momenta[3][5]);
297 numstable++;
298 // p
299 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(2212);
300 p4[numstable].set(ctes_.momenta[0][6],ctes_.momenta[1][6], ctes_.momenta[2][6], ctes_.momenta[3][6]);
301 numstable++;
302 }
303 if (flags_.pion == 5) { // nnbar
304 // pbar
305 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-2112);
306 p4[numstable].set(ctes_.momenta[0][5],ctes_.momenta[1][5], ctes_.momenta[2][5], ctes_.momenta[3][5]);
307 numstable++;
308 // p
309 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(2112);
310 p4[numstable].set(ctes_.momenta[0][6],ctes_.momenta[1][6], ctes_.momenta[2][6], ctes_.momenta[3][6]);
311 numstable++;
312 }
313 if (flags_.pion == 6) { // K+ K-
314 // K+
315 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(321);
316 p4[numstable].set(ctes_.momenta[0][5],ctes_.momenta[1][5], ctes_.momenta[2][5], ctes_.momenta[3][5]);
317 numstable++;
318 // K -
319 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-321);
320 p4[numstable].set(ctes_.momenta[0][6],ctes_.momenta[1][6], ctes_.momenta[2][6], ctes_.momenta[3][6]);
321 numstable++;
322 }
323 if (flags_.pion == 7) { // K0K0bar
324 // Kbar
325 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(311);
326 p4[numstable].set(ctes_.momenta[0][5],ctes_.momenta[1][5], ctes_.momenta[2][5], ctes_.momenta[3][5]);
327 numstable++;
328 // K0
329 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-311);
330 p4[numstable].set(ctes_.momenta[0][6],ctes_.momenta[1][6], ctes_.momenta[2][6], ctes_.momenta[3][6]);
331 numstable++;
332 }
333 if (flags_.pion == 8) { // pi+ pi-pi0
334 // pi+
335 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(211);
336 p4[numstable].set(ctes_.momenta[0][5],ctes_.momenta[1][5], ctes_.momenta[2][5], ctes_.momenta[3][5]);
337 numstable++;
338 // pi-
339 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-211);
340 p4[numstable].set(ctes_.momenta[0][6],ctes_.momenta[1][6], ctes_.momenta[2][6], ctes_.momenta[3][6]);
341 numstable++;
342 // pi0
343 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(111);
344 p4[numstable].set(ctes_.momenta[0][7],ctes_.momenta[1][7], ctes_.momenta[2][7], ctes_.momenta[3][7]);
345 numstable++;
346 }
347 if (flags_.pion == 9) { //Lambda Lambdabar-> pi+ pi- ppbar
348 // pi+
349 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(211);
350 p4[numstable].set(ctes_.momenta[0][7],ctes_.momenta[1][7], ctes_.momenta[2][7], ctes_.momenta[3][7]);
351 numstable++;
352 // pbar
353 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-2212);
354 p4[numstable].set(ctes_.momenta[0][8],ctes_.momenta[1][8], ctes_.momenta[2][8], ctes_.momenta[3][8]);
355 numstable++;
356 // pi-
357 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(-211);
358 p4[numstable].set(ctes_.momenta[0][9],ctes_.momenta[1][9], ctes_.momenta[2][9], ctes_.momenta[3][9]);
359 numstable++;
360 // p
361 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(2212);
362 p4[numstable].set(ctes_.momenta[0][10],ctes_.momenta[1][10], ctes_.momenta[2][10], ctes_.momenta[3][10]);
363 numstable++;
364 }
365
366 // ISR gamma
367 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(22);
368 p4[numstable].set(ctes_.momenta[0][2],ctes_.momenta[1][2], ctes_.momenta[2][2], ctes_.momenta[3][2]);
369 numstable++;
370 if( ctes_.momenta[0][3] != 0 ) // second photon exists
371 {
372 evtnumstable[numstable]=EvtPDL::evtIdFromStdHep(22);
373 p4[numstable].set(ctes_.momenta[0][3],ctes_.momenta[1][3], ctes_.momenta[2][3], ctes_.momenta[3][3]);
374 numstable++;
375 }
376
377 int channel=EvtDecayTable::inChannelList(p->getId(),numstable,evtnumstable);
378 more=(channel!=-1);
379 if(more) {std::cout<<"Existence of mode "<<channel<<" in exclusive decay list has the same final state as this one"<<std::endl;abort(); }
380
381 p->makeDaughters(numstable,evtnumstable);
382 //double weight = p->initializePhaseSpace(getNDaug(),getDaugs());
383
384 int ndaugFound=0;
385 EvtVector4R SUMP4(0,0,0,0);
386 for(int i=0;i<numstable;i++){
387 p->getDaug(i)->init(evtnumstable[i],p4[i]);
388 ndaugFound++;
389 }
390 if ( ndaugFound == 0 ) {
391 report(ERROR,"EvtGen") << "Phokhara has failed to do a decay ";
392 report(ERROR,"EvtGen") << EvtPDL::name(p->getId()).c_str() << " " << p->mass()<<endl;
393 assert(0);
394 }
395
396 nevtgen++;
397 return ;
398
399}
flags_
struct @17 flags_
GEN_1PH
#define GEN_1PH(I, QQMIN, QQMAX, COS1MIN, COS1MAX, COS3MIN, COS3MAX)
Definition: bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtPhokharaDef.hh:100
GEN_0PH
#define GEN_0PH(I, QQMIN, SP, COS3MIN, COS3MAX)
Definition: bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtPhokharaDef.hh:106
maxima_
struct @16 maxima_
ctes_
struct @10 ctes_
GEN_2PH
#define GEN_2PH(I, QQMIN, COS1MIN, COS1MAX, COS2MIN, COS2MAX, COS3MIN, COS3MAX)
Definition: bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtPhokharaDef.hh:103
EvtDecayTable::inChannelList
static int inChannelList(EvtId parent, int ndaug, EvtId *daugs)
Definition: bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDecayTable.cc:795
EvtPDL::getId
static EvtId getId(const std::string &name)
Definition: bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtPDL.cc:287
EvtParticle::init
virtual void init(EvtId part_n, const EvtVector4R &p4)=0
EvtPhokhara_KK::init_evt
void init_evt(EvtParticle *p)
Definition: bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtPhokhara_KK.cc:437
EvtPhokhara_KK::init_mode
void init_mode(EvtParticle *p)
Definition: bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtPhokhara_KK.cc:115

◆ decay() [2/2]

void EvtPhokhara_KK::decay ( EvtParticle *  p)
virtual

Implements EvtDecayBase.

◆ ExclusiveDecay() [1/2]

void EvtPhokhara_KK::ExclusiveDecay ( EvtParticle *  p)

◆ ExclusiveDecay() [2/2]

void EvtPhokhara_KK::ExclusiveDecay ( EvtParticle *  p)

◆ getName() [1/2]

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

Implements EvtDecayBase.

Definition at line 95 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtPhokhara_KK.cc.

95 {
96
97 model_name="PHOKHARA_KK";
98
99}

◆ getName() [2/2]

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

Implements EvtDecayBase.

◆ getTotalEvt() [1/2]

int EvtPhokhara_KK::getTotalEvt ( )
inline

Definition at line 52 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtPhokhara_KK.hh.

52{return nevt;}

◆ getTotalEvt() [2/2]

int EvtPhokhara_KK::getTotalEvt ( )
inline

Definition at line 52 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtPhokhara_KK.hh.

52{return nevt;}

◆ init() [1/2]

void EvtPhokhara_KK::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 126 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtPhokhara_KK.cc.

126 {
127 checkNArg(0);
128
129 std::string locvp=getenv("BESEVTGENROOT");
130 system("cat $BESEVTGENROOT/share/phokhara_9.1.param>phokhara_9.1.param");
131 system("cat $BESEVTGENROOT/share/phokhara_9.1.fferr>phokhara_9.1.fferr");
132 system("cat $BESEVTGENROOT/share/phokhara_9.1.ffwarn>phokhara_9.1.ffwarn");
133
134
135 if (getParentId().isAlias()){
136
137 report(ERROR,"EvtGen") << "EvtPhokhara finds that you are decaying the"<<endl
138 << " aliased particle "
139 << EvtPDL::name(getParentId()).c_str()
140 << " with the Phokhara model"<<endl
141 << " this does not work, please modify decay table."
142 << endl;
143 report(ERROR,"EvtGen") << "Will terminate execution!"<<endl;
144 ::abort();
145
146 }
147
148 store(this);
149
150}
EvtDecayBase::checkNArg
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
Definition: bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDecayBase.cc:482

◆ init() [2/2]

void EvtPhokhara_KK::init ( )
virtual

Reimplemented from EvtDecayBase.

◆ init_evt() [1/2]

void EvtPhokhara_KK::init_evt ( EvtParticle *  p)

Definition at line 437 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtPhokhara_KK.cc.

437 {
438 m_pion=6;
439 // mu+mu-(0),pi+pi-(1),2pi0pi+pi-(2),
440 // 2pi+2pi-(3),ppbar(4),nnbar(5),
441 // K+K-(6),K0K0bar(7),pi+pi-pi0(8),
442 // Lamb Lambbar->pi-pi+ppbar(9)
443#include "Phokhara_init_evt.txt"
444}

Referenced by decay().

◆ init_evt() [2/2]

void EvtPhokhara_KK::init_evt ( EvtParticle *  p)

◆ init_mode() [1/2]

void EvtPhokhara_KK::init_mode ( EvtParticle *  p)

Definition at line 115 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtPhokhara_KK.cc.

115 {
116 m_pion=6;
117 // mu+mu-(0),pi+pi-(1),2pi0pi+pi-(2),
118 // 2pi+2pi-(3),ppbar(4),nnbar(5),
119 // K+K-(6),K0K0bar(7),pi+pi-pi0(8),
120 // Lamb Lambbar->pi-pi+ppbar(9)
121#include "Phokhara_init_mode.txt"
122}

Referenced by decay().

◆ init_mode() [2/2]

void EvtPhokhara_KK::init_mode ( EvtParticle *  p)

◆ initProbMax() [1/2]

void EvtPhokhara_KK::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 108 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtPhokhara_KK.cc.

108 {
109
110 noProbMax();
111
112}

◆ initProbMax() [2/2]

void EvtPhokhara_KK::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

◆ PhokharaInit() [1/2]

void EvtPhokhara_KK::PhokharaInit ( int  dummy)

Definition at line 424 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtPhokhara_KK.cc.

424 {
425 static int first=1;
426 if (first){
427
428 first=0;
429 //for(int i=0;i<ncommand;i++)
430 // lugive0_(commands[i].c_str(),strlen(commands[i].c_str()));
431 }
432
433}

◆ PhokharaInit() [2/2]

void EvtPhokhara_KK::PhokharaInit ( int  dummy)
The documentation for this class was generated from the following files: