Changeset View
Changeset View
Standalone View
Standalone View
src/EvtGenModels/EvtBcVHad.cpp
Show First 20 Lines • Show All 55 Lines • ▼ Show 20 Lines | void EvtBcVHad::init() | ||||
idVector = getDaug( 0 ).getId(); | idVector = getDaug( 0 ).getId(); | ||||
whichfit = int( getArg( 0 ) + 0.1 ); | whichfit = int( getArg( 0 ) + 0.1 ); | ||||
ffmodel = std::make_unique<EvtBCVFF2>( idVector, whichfit ); | ffmodel = std::make_unique<EvtBCVFF2>( idVector, whichfit ); | ||||
wcurr = std::make_unique<EvtWHad>(); | wcurr = std::make_unique<EvtWHad>(); | ||||
// determine the code of final hadronic state | // determine the code of final hadronic state | ||||
EvtIdSet thePis( "pi+", "pi-", "pi0" ); | EvtIdSet thePis{ "pi+", "pi-", "pi0" }; | ||||
EvtIdSet theK( "K+", "K-", "K_S0" ); | EvtIdSet theK{ "K+", "K-", "K_S0" }; | ||||
if ( getNDaug() == 2 && thePis.contains( getDaug( 1 ) ) ) { | if ( getNDaug() == 2 && thePis.contains( getDaug( 1 ) ) ) { | ||||
out_code = 1; // pi+ | out_code = 1; // pi+ | ||||
} else if ( getNDaug() == 3 && thePis.contains( getDaug( 1 ) ) && | } else if ( getNDaug() == 3 && thePis.contains( getDaug( 1 ) ) && | ||||
thePis.contains( getDaug( 2 ) ) ) { | thePis.contains( getDaug( 2 ) ) ) { | ||||
out_code = 2; // pi+ pi0 | out_code = 2; // pi+ pi0 | ||||
} else if ( getNDaug() == 4 && thePis.contains( getDaug( 1 ) ) && | } else if ( getNDaug() == 4 && thePis.contains( getDaug( 1 ) ) && | ||||
thePis.contains( getDaug( 2 ) ) && | thePis.contains( getDaug( 2 ) ) && | ||||
thePis.contains( getDaug( 3 ) ) ) { | thePis.contains( getDaug( 3 ) ) ) { | ||||
▲ Show 20 Lines • Show All 303 Lines • Show Last 20 Lines |