BOSS 7.0.4
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtDMix Class Reference

#include <EvtDMix.hh>

+ Inheritance diagram for EvtDMix:

Public Member Functions

 EvtDMix ()
 
virtual ~EvtDMix ()
 
void getName (std::string &name)
 
EvtDecayBaseclone ()
 
void initProbMax ()
 
void init ()
 
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 void getName (std::string &name)=0
 
virtual void decay (EvtParticle *p)=0
 
virtual void makeDecay (EvtParticle *p)=0
 
virtual EvtDecayBaseclone ()=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 ()
 
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 30 of file EvtDMix.hh.

Constructor & Destructor Documentation

◆ EvtDMix()

EvtDMix::EvtDMix ( )
inline

Definition at line 34 of file EvtDMix.hh.

34{}

Referenced by clone().

◆ ~EvtDMix()

EvtDMix::~EvtDMix ( )
virtual

Definition at line 31 of file EvtDMix.cc.

31{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtDMix::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 39 of file EvtDMix.cc.

39 {
40
41 return new EvtDMix;
42
43}
EvtDMix()
Definition: EvtDMix.hh:34

◆ decay()

void EvtDMix::decay ( EvtParticle p)
virtual

Implements EvtDecayBase.

Definition at line 61 of file EvtDMix.cc.

61 {
62
63 //unneeded - lange - may13-02
64 //if ( p->getNDaug() != 0 ) {
65 //Will end up here because maxrate multiplies by 1.2
66 // report(DEBUG,"EvtGen") << "In EvtDMix: has "
67 // <<" daugthers should not be here!"<<endl;
68 // return;
69 //}
70
72
73 double ctau=EvtPDL::getctau(p->getId());
74 if ( ctau==0. ) return;
75
76
77 double pdf, random, gt, weight;
78
79 double maxPdf=_rd + sqrt(_rd)*_ypr*50. + 2500.0*(_xpr*_xpr+_ypr*_ypr)/4.0;
80 bool keepGoing=true;
81 while ( keepGoing ) {
82 random=EvtRandom::Flat();
83 gt=-log(random);
84 weight=random;
85 pdf=_rd + sqrt(_rd)*_ypr*gt + gt*gt*(_xpr*_xpr+_ypr*_ypr)/4.0;
86 pdf*=exp(-1.0*gt);
87 pdf/=weight;
88 if ( pdf > maxPdf )
89 std::cout << pdf << " " << weight << " " << maxPdf << " " << gt << std::endl;
90 if ( pdf > maxPdf*EvtRandom::Flat()) keepGoing=false;
91 }
92
93
94 p->setLifetime(gt*ctau);
95
96 return ;
97}
EvtComplex exp(const EvtComplex &c)
Definition: EvtComplex.hh:252
*********Class see also m_nmax DOUBLE PRECISION m_MasPhot DOUBLE PRECISION m_phsu DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_r2 DOUBLE PRECISION m_WtMass INTEGER m_nmax INTEGER m_Nevgen INTEGER m_IsFSR INTEGER m_MarTot *COMMON c_KarFin $ !Output file $ !Event serial number $ !alpha QED at Thomson limit $ !minimum energy at CMS for remooval $ !infrared dimensionless $ !dummy photon IR regulator $ !crude photon multiplicity enhancement factor *EVENT $ !MC crude volume of PhhSpace *Sfactors $ !YFS formfactor IR part only $ !YFS formfactor non IR finite part $ !mass weight
Definition: KarFin.h:34
EvtId * getDaugs()
Definition: EvtDecayBase.hh:65
static double getctau(EvtId i)
Definition: EvtPDL.hh:55
EvtId getId() const
Definition: EvtParticle.cc:113
void setLifetime(double tau)
Definition: EvtParticle.cc:89
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
static double Flat()
Definition: EvtRandom.cc:74

◆ getName()

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

Implements EvtDecayBase.

Definition at line 33 of file EvtDMix.cc.

33 {
34
35 model_name="DMIX";
36
37}

◆ init()

void EvtDMix::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 46 of file EvtDMix.cc.

46 {
47
48 // check that there are 0 arguments
49 checkNArg(3);
50 _rd=getArg(0);
51 _xpr=getArg(1);
52 _ypr=getArg(2);
53}
double getArg(int j)
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)

◆ initProbMax()

void EvtDMix::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 55 of file EvtDMix.cc.

55 {
56
57 noProbMax();
58
59}
void noProbMax()

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