Page MenuHomeHEPForge
Differential D91 Diff 387 src/EvtGenModels/EvtSLBaryonAmp.cpp

Changeset View

Standalone View

View Options

src/EvtGenModels/EvtSLBaryonAmp.cpp

Show First 20 Lines • Show All 421 Lines • ▼ Show 20 Linesvoid EvtSLBaryonAmp::CalcAmp( EvtParticle* parent, EvtAmp& amp,
rho.set( 0, 0, r00 ); rho.set( 0, 0, r00 );
rho.set( 0, 1, r01 ); rho.set( 0, 1, r01 );
rho.set( 1, 0, r10 ); rho.set( 1, 0, r10 );
rho.set( 1, 1, r11 ); rho.set( 1, 1, r11 );
EvtVector4R vector4P = parent->getP4Lab(); EvtVector4R vector4P = parent->getP4Lab();
double pmag = vector4P.d3mag(); double pmag = vector4P.d3mag();
double cosTheta = vector4P.get( 3 ) / pmag; double cosTheta = pmag > 0.0 ? vector4P.get( 3 ) / pmag : 1.0;
tlathamUnsubmitted
Not Done
Inline Actions
double pmag = vector4P.d3mag();
- double cosTheta = pmag > 0.0 ? vector4P.get( 3 ) / pmag : 0.0;
+ double cosTheta = pmag > 0.0 ? vector4P.get( 3 ) / pmag : 1.0;
double theta = acos( cosTheta );

Should maybe be set to 1.0 rather than 0,0?

tlatham: Should maybe be set to 1.0 rather than 0,0?
double theta = acos( cosTheta ); double theta = acos( cosTheta );
double phi = atan2( vector4P.get( 2 ), vector4P.get( 1 ) ); double phi = atan2( vector4P.get( 2 ), vector4P.get( 1 ) );
parent->setSpinDensityForwardHelicityBasis( rho, phi, theta, 0.0 ); parent->setSpinDensityForwardHelicityBasis( rho, phi, theta, 0.0 );
//parent->setSpinDensityForward(rho); //parent->setSpinDensityForward(rho);
// Set the four momentum of the parent baryon in it's rest frame // Set the four momentum of the parent baryon in it's rest frame
Show All 17 Linesvoid EvtSLBaryonAmp::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 106 Lines • ▼ Show 20 Linesvoid EvtSLBaryonAmp::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 254 Lines • Show Last 20 Lines