Changeset View
Changeset View
Standalone View
Standalone View
src/EvtGenModels/EvtRareLbToLllFF.cpp
| Show First 20 Lines • Show All 148 Lines • ▼ Show 20 Lines | void EvtRareLbToLllFF::init() | ||||
| L1520->H5.param( -0.0517, -0.0173, 0.00259, 0.333, 0.308 ); | L1520->H5.param( -0.0517, -0.0173, 0.00259, 0.333, 0.308 ); | ||||
| L1520->H6.param( 0.0206, 0.00679, -0.000220, 0.333, 0.308 ); | L1520->H6.param( 0.0206, 0.00679, -0.000220, 0.333, 0.308 ); | ||||
| FFMap_[EvtPDL::getId( "Lambda0" ).getId()] = L1115.get(); | FFMap_[EvtPDL::getId( "Lambda0" ).getId()] = L1115.get(); | ||||
| FFMap_[EvtPDL::getId( "anti-Lambda0" ).getId()] = L1115.get(); | FFMap_[EvtPDL::getId( "anti-Lambda0" ).getId()] = L1115.get(); | ||||
| FFMap_[EvtPDL::getId( "Lambda(1520)0" ).getId()] = L1520.get(); | FFMap_[EvtPDL::getId( "Lambda(1520)0" ).getId()] = L1520.get(); | ||||
| FFMap_[EvtPDL::getId( "anti-Lambda(1520)0" ).getId()] = L1520.get(); | FFMap_[EvtPDL::getId( "anti-Lambda(1520)0" ).getId()] = L1520.get(); | ||||
| FF_ = {std::move( L1115 ), std::move( L1520 )}; | FF_ = { std::move( L1115 ), std::move( L1520 ) }; | ||||
| EvtGenReport( EVTGEN_INFO, "EvtGen" ) | EvtGenReport( EVTGEN_INFO, "EvtGen" ) | ||||
| << " EvtRareLbToLll is using form factors from arXiv:1108.6129 " | << " EvtRareLbToLll is using form factors from arXiv:1108.6129 " | ||||
| << std::endl; | << std::endl; | ||||
| } | } | ||||
| //============================================================================= | //============================================================================= | ||||
| ▲ Show 20 Lines • Show All 107 Lines • ▼ Show 20 Lines | if ( isNatural( lambda ) ) { | ||||
| FF.GT_[0] = -( MB + M ) * H1 + M * ( 1. + vdotv ) * H2 + | FF.GT_[0] = -( MB + M ) * H1 + M * ( 1. + vdotv ) * H2 + | ||||
| MB * ( 1. + vdotv ) * H3 + MB * H5 + M * H6; | MB * ( 1. + vdotv ) * H3 + MB * H5 + M * H6; | ||||
| FF.GT_[1] = MB * H1 - M * H2 - MB * H3; | FF.GT_[1] = MB * H1 - M * H2 - MB * H3; | ||||
| FF.GT_[2] = M * H1 - M * H2 - MB * H3 - M * H6; | FF.GT_[2] = M * H1 - M * H2 - MB * H3 - M * H6; | ||||
| FF.GT_[3] = -( MB - M ) * H5 - M * ( 1. - vdotv ) * H6; | FF.GT_[3] = -( MB - M ) * H5 - M * ( 1. - vdotv ) * H6; | ||||
| } | } | ||||
| } | } | ||||
| void EvtRareLbToLllFF::getFF( const EvtParticle& parent, const EvtParticle& lambda, | void EvtRareLbToLllFF::getFF( const EvtParticle& parent, | ||||
| const EvtParticle& lambda, | |||||
| EvtRareLbToLllFF::FormFactors& FF ) const | EvtRareLbToLllFF::FormFactors& FF ) const | ||||
| { | { | ||||
| // Find the FF information for this particle, start by setting all to zero | // Find the FF information for this particle, start by setting all to zero | ||||
| FF.areZero(); | FF.areZero(); | ||||
| // Are the FF's for the particle known? | // Are the FF's for the particle known? | ||||
| auto it = FFMap_.find( lambda.getId().getId() ); | auto it = FFMap_.find( lambda.getId().getId() ); | ||||
| Show All 23 Lines | |||||