Changeset View
Changeset View
Standalone View
Standalone View
src/EvtGenModels/EvtBsMuMuKK.cpp
Show First 20 Lines • Show All 320 Lines • ▼ Show 20 Lines | EvtComplex aperp_f2p = AmpTime( q, gplus, gminus, delta_f2p_perp, | ||||
phis_f2p_perp, 1 ); | phis_f2p_perp, 1 ); | ||||
EvtComplex apar_f2p = AmpTime( q, gplus, gminus, delta_f2p_par, | EvtComplex apar_f2p = AmpTime( q, gplus, gminus, delta_f2p_par, | ||||
lambda_f2p_par_abs, A_f2p_par, phis_f2p_par, | lambda_f2p_par_abs, A_f2p_par, phis_f2p_par, | ||||
-1 ); | -1 ); | ||||
// Generate 4-momenta | // Generate 4-momenta | ||||
double mKK = EvtRandom::Flat( kin_lower_limit, kin_upper_limit ); | double mKK = EvtRandom::Flat( kin_lower_limit, kin_upper_limit ); | ||||
double mass[10] = {MJpsi, mKK, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0}; | double mass[10] = { MJpsi, mKK, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 }; | ||||
double Kmass[10] = {MKp, MKm, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0}; | double Kmass[10] = { MKp, MKm, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 }; | ||||
double muMass[10] = {Mmu, Mmu, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0}; | double muMass[10] = { Mmu, Mmu, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 }; | ||||
EvtVector4R mypV[2], mypK[2], mypmu[2]; | EvtVector4R mypV[2], mypK[2], mypmu[2]; | ||||
EvtGenKine::PhaseSpace( 2, mass, mypV, MBs ); | EvtGenKine::PhaseSpace( 2, mass, mypV, MBs ); | ||||
EvtGenKine::PhaseSpace( 2, Kmass, mypK, mKK ); | EvtGenKine::PhaseSpace( 2, Kmass, mypK, mKK ); | ||||
EvtGenKine::PhaseSpace( 2, muMass, mypmu, MJpsi ); | EvtGenKine::PhaseSpace( 2, muMass, mypmu, MJpsi ); | ||||
EvtVector4R p4mup = boostTo( mypmu[0], mypV[0] ); | EvtVector4R p4mup = boostTo( mypmu[0], mypV[0] ); | ||||
EvtVector4R p4mum = boostTo( mypmu[1], mypV[0] ); | EvtVector4R p4mum = boostTo( mypmu[1], mypV[0] ); | ||||
▲ Show 20 Lines • Show All 343 Lines • Show Last 20 Lines |