Page MenuHomeHEPForge
Differential D90 Diff 388 src/EvtGenModels/EvtLambdacPHH.cpp

Changeset View

Standalone View

View Options

src/EvtGenModels/EvtLambdacPHH.cpp

Show First 20 Lines • Show All 162 Lines • ▼ Show 20 Lines
{ {
// Generate events uniform in the Lambda_c Dalitz plot and find the // Generate events uniform in the Lambda_c Dalitz plot and find the
// normalisation integrals of the Breit-Wigner lineshapes // normalisation integrals of the Breit-Wigner lineshapes
// Lambda_c -> K- pi+ p // Lambda_c -> K- pi+ p
int nDaug( 3 ); int nDaug( 3 );
EvtVector4R p4Daug[3]; EvtVector4R p4Daug[3];
double mDaug[3] = {EvtPDL::getMeanMass( EvtPDL::getId( "K-" ) ), double mDaug[3] = { EvtPDL::getMeanMass( EvtPDL::getId( "K-" ) ),
EvtPDL::getMeanMass( EvtPDL::getId( "pi+" ) ), EvtPDL::getMeanMass( EvtPDL::getId( "pi+" ) ),
EvtPDL::getMeanMass( EvtPDL::getId( "p+" ) )}; EvtPDL::getMeanMass( EvtPDL::getId( "p+" ) ) };
double norm[3] = {0.0, 0.0, 0.0}; double norm[3] = { 0.0, 0.0, 0.0 };
// sample size // sample size
int N( 100000 ); int N( 100000 );
for ( int i = 0; i < N; i++ ) { for ( int i = 0; i < N; i++ ) {
double mParent = EvtPDL::getMass( EvtPDL::getId( "Lambda_c+" ) ); double mParent = EvtPDL::getMass( EvtPDL::getId( "Lambda_c+" ) );
EvtVector4R p0( mParent, 0.0, 0.0, 0.0 ); EvtVector4R p0( mParent, 0.0, 0.0, 0.0 );
// Generate uniform 4 momenta // Generate uniform 4 momenta
Show All 40 Lines
{ {
// Generate events uniform in the Lambda_c Dalitz plot and find the // Generate events uniform in the Lambda_c Dalitz plot and find the
// fit fractions for each resonance // fit fractions for each resonance
// Lambda_c -> K- pi+ p // Lambda_c -> K- pi+ p
int nDaug( 3 ); int nDaug( 3 );
EvtVector4R p4Daug[3]; EvtVector4R p4Daug[3];
double mDaug[3] = {EvtPDL::getMeanMass( EvtPDL::getId( "K-" ) ), double mDaug[3] = { EvtPDL::getMeanMass( EvtPDL::getId( "K-" ) ),
EvtPDL::getMeanMass( EvtPDL::getId( "pi+" ) ), EvtPDL::getMeanMass( EvtPDL::getId( "pi+" ) ),
EvtPDL::getMeanMass( EvtPDL::getId( "p+" ) )}; EvtPDL::getMeanMass( EvtPDL::getId( "p+" ) ) };
double FitFracTop[4] = {0.0, 0.0, 0.0, 0.0}; double FitFracTop[4] = { 0.0, 0.0, 0.0, 0.0 };
double FitFracDenom = 0.0; double FitFracDenom = 0.0;
// sample size // sample size
int N( 100000 ); int N( 100000 );
for ( int i = 0; i < N; i++ ) { for ( int i = 0; i < N; i++ ) {
double mParent = EvtPDL::getMass( EvtPDL::getId( "Lambda_c+" ) ); double mParent = EvtPDL::getMass( EvtPDL::getId( "Lambda_c+" ) );
EvtVector4R p0( mParent, 0.0, 0.0, 0.0 ); EvtVector4R p0( mParent, 0.0, 0.0, 0.0 );
▲ Show 20 Lines • Show All 437 Lines • Show Last 20 Lines