94 static double ctauLong =
EvtConst::c / ( Gamma - fabs( deltaGamma ) / 2 );
126 double mt =
exp( -std::max( 0., deltaGamma ) * t / ( 2 *
EvtConst::c ) );
127 double pt =
exp( +std::min( 0., deltaGamma ) * t / ( 2 *
EvtConst::c ) );
143 if ( other_b == BSB ) {
149 cG0P = G0P * ( gplus + lambda_km * gminus );
150 cG1P = G1P * ( gplus + lambda_km * gminus );
151 cG1M = G1M * ( gplus - lambda_km * gminus );
152 }
else if ( other_b == BS0 ) {
156 cG0P = G0P * ( gplus + ( 1.0 / lambda_km ) * gminus );
157 cG1P = G1P * ( gplus + ( 1.0 / lambda_km ) * gminus );
158 cG1M = -G1M * ( gplus - ( 1.0 / lambda_km ) * gminus );
162 <<
"other_b was not BSB or BS0!" << std::endl;
171 AP = ( cG1P + cG1M ) / sqrt( 2.0 );
172 AM = ( cG1P - cG1M ) / sqrt( 2.0 );
187 if ( ( p->
getId() != BS0 ) && ( p->
getId() != BSB ) )
EvtParticle * getParent() const
double getArg(unsigned int j)
static void SVVHel(EvtParticle *parent, EvtAmp &, EvtId n_v1, EvtId n_v2, const EvtComplex &hp, const EvtComplex &h0, const EvtComplex &hm)
std::ostream & EvtGenReport(EvtGenSeverity severity, const char *facility=0)
double getDeltaGamma(const EvtId id)
void initProbMax() override
static double getctau(EvtId i)
EvtDecayBase * clone() override
bool isBsMixed(EvtParticle *p)
void OtherB(EvtParticle *p, double &t, EvtId &otherb)
void setLifetime(double tau)
void setProbMax(double prbmx)
void decay(EvtParticle *p) override
void checkNDaug(int d1, int d2=-1)
void checkSpinParent(EvtSpinType::spintype sp)
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
static EvtId getId(const std::string &name)
void checkSpinDaughter(int d1, EvtSpinType::spintype sp)
EvtComplex exp(const EvtComplex &c)
static EvtCPUtil * getInstance()
double getDeltaM(const EvtId id)
std::string getName() override
EvtId getDaug(int i) const