Changeset View
Changeset View
Standalone View
Standalone View
src/EvtGenBase/EvtSemiLeptonicBaryonAmp.cpp
| Show First 20 Lines • Show All 433 Lines • ▼ Show 20 Lines | void EvtSemiLeptonicBaryonAmp::CalcAmp( EvtParticle* parent, EvtAmp& amp, | ||||
| // Handle spin-1/2 daughter baryon Dirac spinor cases | // Handle spin-1/2 daughter baryon Dirac spinor cases | ||||
| if ( EvtPDL::getSpinType( parent->getDaug( 0 )->getId() ) == | if ( EvtPDL::getSpinType( parent->getDaug( 0 )->getId() ) == | ||||
| EvtSpinType::DIRAC ) { | EvtSpinType::DIRAC ) { | ||||
| // Set the form factors | // Set the form factors | ||||
| double f1, f2, f3, g1, g2, g3; | double f1, f2, f3, g1, g2, g3; | ||||
| FormFactors->getdiracff( par_num, bar_num, q2, baryonmass, &f1, &f2, | FormFactors->getdiracff( par_num, bar_num, q2, baryonmass, &f1, &f2, | ||||
| &f3, &g1, &g2, &g3 ); | &f3, &g1, &g2, &g3 ); | ||||
| const double form_fact[6] = {f1, f2, f3, g1, g2, g3}; | const double form_fact[6] = { f1, f2, f3, g1, g2, g3 }; | ||||
| EvtVector4C b11, b12, b21, b22, l1, l2; | EvtVector4C b11, b12, b21, b22, l1, l2; | ||||
| // Lepton Current | // Lepton Current | ||||
| if ( l_num == EM || l_num == MUM || l_num == TAUM ) { | if ( l_num == EM || l_num == MUM || l_num == TAUM ) { | ||||
| l1 = EvtLeptonVACurrent( parent->getDaug( 1 )->spParent( 0 ), | l1 = EvtLeptonVACurrent( parent->getDaug( 1 )->spParent( 0 ), | ||||
| parent->getDaug( 2 )->spParentNeutrino() ); | parent->getDaug( 2 )->spParentNeutrino() ); | ||||
| l2 = EvtLeptonVACurrent( parent->getDaug( 1 )->spParent( 1 ), | l2 = EvtLeptonVACurrent( parent->getDaug( 1 )->spParent( 1 ), | ||||
| ▲ Show 20 Lines • Show All 92 Lines • ▼ Show 20 Lines | void EvtSemiLeptonicBaryonAmp::CalcAmp( EvtParticle* parent, EvtAmp& amp, | ||||
| // Rarita-Schwinger spinor cases | // Rarita-Schwinger spinor cases | ||||
| else if ( EvtPDL::getSpinType( parent->getDaug( 0 )->getId() ) == | else if ( EvtPDL::getSpinType( parent->getDaug( 0 )->getId() ) == | ||||
| EvtSpinType::RARITASCHWINGER ) { | EvtSpinType::RARITASCHWINGER ) { | ||||
| // Set the form factors | // Set the form factors | ||||
| double f1, f2, f3, f4, g1, g2, g3, g4; | double f1, f2, f3, f4, g1, g2, g3, g4; | ||||
| FormFactors->getraritaff( par_num, bar_num, q2, baryonmass, &f1, &f2, | FormFactors->getraritaff( par_num, bar_num, q2, baryonmass, &f1, &f2, | ||||
| &f3, &f4, &g1, &g2, &g3, &g4 ); | &f3, &f4, &g1, &g2, &g3, &g4 ); | ||||
| const double form_fact[8] = {f1, f2, f3, f4, g1, g2, g3, g4}; | const double form_fact[8] = { f1, f2, f3, f4, g1, g2, g3, g4 }; | ||||
| EvtId l_num = parent->getDaug( 1 )->getId(); | EvtId l_num = parent->getDaug( 1 )->getId(); | ||||
| EvtVector4C b11, b12, b21, b22, b13, b23, b14, b24, l1, l2; | EvtVector4C b11, b12, b21, b22, b13, b23, b14, b24, l1, l2; | ||||
| // Lepton Current | // Lepton Current | ||||
| if ( l_num == EM || l_num == MUM || l_num == TAUM ) { | if ( l_num == EM || l_num == MUM || l_num == TAUM ) { | ||||
| // Lepton Current | // Lepton Current | ||||
| ▲ Show 20 Lines • Show All 229 Lines • Show Last 20 Lines | |||||