-
Notifications
You must be signed in to change notification settings - Fork 107
Open
Description
In lines 475-488 of the /src/Physics/HEDIS/XSection/HEDISStrucFunc.cxx there is a check for whether either LHAPDF5 or LHAPDF6 are enabled while configuring, if they are it defines a double fPDF.
If neither are enabled, it throws an error because it then goes on to try and access that double the next few lines down.
// Extract PDF requiring then to be higher than zero
#ifdef __GENIE_LHAPDF6_ENABLED__
double fPDF = fmax( pdf->xfxQ2(qpdf1, xPDF, Q2PDF)/z , 0.);
if (qpdf2!= -999) fPDF -= fmax( pdf->xfxQ2(qpdf2, xPDF, Q2PDF)/z , 0.);
#endif
#ifdef __GENIE_LHAPDF5_ENABLED__
double fPDF = fmax( LHAPDF::xfx(xPDF, TMath::Sqrt(Q2PDF), qpdf1)/z , 0.);
if (qpdf2!= -999) fPDF -= fmax( LHAPDF::xfx(xPDF, TMath::Sqrt(Q2PDF), qpdf2)/z , 0.);
#endif
// Compute SF
double tmp = -999;
if ( sf==1 ) tmp = fPDF*Cp2/2;
else if ( sf==2 ) tmp = fPDF*Cp2*z;
else if ( sf==3 ) tmp = fPDF*Cp3*sign3;
Reactions are currently unavailable
Metadata
Metadata
Assignees
Labels
No labels