BOSS 7.1.2
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtVVPIPI_WEIGHTED Class Reference

#include <EvtVVPIPI_WEIGHTED.hh>

+ Inheritance diagram for EvtVVPIPI_WEIGHTED:

Public Member Functions

 EvtVVPIPI_WEIGHTED ()
 
virtual ~EvtVVPIPI_WEIGHTED ()
 
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 30 of file EvtVVPIPI_WEIGHTED.hh.

Constructor & Destructor Documentation

◆ EvtVVPIPI_WEIGHTED()

EvtVVPIPI_WEIGHTED::EvtVVPIPI_WEIGHTED ( )
inline

Definition at line 34 of file EvtVVPIPI_WEIGHTED.hh.

34{}

Referenced by clone().

◆ ~EvtVVPIPI_WEIGHTED()

EvtVVPIPI_WEIGHTED::~EvtVVPIPI_WEIGHTED ( )
virtual

Definition at line 34 of file EvtVVPIPI_WEIGHTED.cc.

34{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtVVPIPI_WEIGHTED::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 43 of file EvtVVPIPI_WEIGHTED.cc.

43 {
44
45 return new EvtVVPIPI_WEIGHTED;
46
47}

◆ decay()

void EvtVVPIPI_WEIGHTED::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 90 of file EvtVVPIPI_WEIGHTED.cc.

90 {
91
92 psi_prime->initializePhaseSpace(getNDaug(),getDaugs());
93
94 EvtParticle *jpsi,*pi1,*pi2;
95
96 jpsi=psi_prime->getDaug(0);
97 pi1=psi_prime->getDaug(1);
98 pi2=psi_prime->getDaug(2);
99
100// Put phase space results into the daughters.
101
102 EvtVector4C ep0,ep1,ep2;
103
104 ep0=psi_prime->eps(0);
105 ep1=psi_prime->eps(1);
106 ep2=psi_prime->eps(2);
107
108 EvtVector4C e0,e1,e2;
109
110 e0 = jpsi->epsParent(0);
111 e1 = jpsi->epsParent(1);
112 e2 = jpsi->epsParent(2);
113
114 double mass2 = (pi1->getP4()+pi2->getP4()).mass2();
115
116 double fac = mass2-4*pi1->mass()*pi2->mass();
117
118 fac *= reweight_event(sqrt(mass2));
119
120 vertex(0,0,fac*(ep0*e0.conj()));
121 vertex(0,1,fac*(ep0*e1.conj()));
122 vertex(0,2,fac*(ep0*e2.conj()));
123
124 vertex(1,0,fac*(ep1*e0.conj()));
125 vertex(1,1,fac*(ep1*e1.conj()));
126 vertex(1,2,fac*(ep1*e2.conj()));
127
128 vertex(2,0,fac*(ep2*e0.conj()));
129 vertex(2,1,fac*(ep2*e1.conj()));
130 vertex(2,2,fac*(ep2*e2.conj()));
131
132 return ;
133
134}
Double_t e1
Double_t e2
character *LEPTONflag integer iresonances real pi2
double reweight_event(double pipi_mass)
void vertex(const EvtComplex &amp)
EvtId * getDaugs()
virtual EvtVector4C epsParent(int i) const
const EvtVector4R & getP4() const
EvtParticle * getDaug(int i)
double mass() const
EvtVector4C conj() const

◆ getName()

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

Implements EvtDecayBase.

Definition at line 36 of file EvtVVPIPI_WEIGHTED.cc.

36 {
37
38 model_name="VVPIPI_WEIGHTED";
39
40}

◆ init()

void EvtVVPIPI_WEIGHTED::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 49 of file EvtVVPIPI_WEIGHTED.cc.

49 {
50
51 static EvtId PIP=EvtPDL::getId("pi+");
52 static EvtId PIM=EvtPDL::getId("pi-");
53 static EvtId PI0=EvtPDL::getId("pi0");
54
55 // check that there are 0 arguments
56 checkNArg(0);
57 checkNDaug(3);
58
61
62
63
64 if ((!(getDaug(1)==PIP&&getDaug(2)==PIM))&&
65 (!(getDaug(1)==PI0&&getDaug(2)==PI0))) {
66 report(ERROR,"EvtGen") << "EvtVVPIPI_WEIGHTED generator expected "
67 << " pi+ and pi- (or pi0 and pi0) "
68 << "as 2nd and 3rd daughter. "<<endl;
69 report(ERROR,"EvtGen") << "Will terminate execution!"<<endl;
70 ::abort();
71 }
72
73}
ostream & report(Severity severity, const char *facility)
Definition EvtReport.cc:36
@ 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)
Definition EvtId.hh:27
static EvtId getId(const std::string &name)
Definition EvtPDL.cc:287

◆ initProbMax()

void EvtVVPIPI_WEIGHTED::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 75 of file EvtVVPIPI_WEIGHTED.cc.

75 {
76 //Hard coded... should not be hard to calculate...
77 setProbMax(0.08* 1.13 );
78}
void setProbMax(double prbmx)

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