diff --git a/GeneratorInterface/HydjetInterface/interface/HydjetHadronizer.h b/GeneratorInterface/HydjetInterface/interface/HydjetHadronizer.h index 564bf767b9cf..16677a39dea5 100644 --- a/GeneratorInterface/HydjetInterface/interface/HydjetHadronizer.h +++ b/GeneratorInterface/HydjetInterface/interface/HydjetHadronizer.h @@ -137,14 +137,21 @@ namespace gen { Return the nuclear radius derived from the beam/target atomic mass number. */ - double HydjetHadronizer::nuclear_radius() const { - if (abs(abeamtarget_ - 208.0) < 0.0001) return 6.624; - else if (abs(abeamtarget_ - 197.0) < 0.0001) return 6.38; - else if (abs(abeamtarget_ - 129.0) < 0.0001) return 5.36; - else if (abs(abeamtarget_ - 63.0) < 0.0001) return 4.20; - else if (abs(abeamtarget_ - 20.0) < 0.0001) return 2.992; - else if (abs(abeamtarget_ - 16.0) < 0.0001) return 2.608; - else return 1.15 * pow((double)abeamtarget_, 1. / 3.); + double HydjetHadronizer::nuclear_radius() const { + if (abs(abeamtarget_ - 208.0) < 0.0001) + return 6.624; + else if (abs(abeamtarget_ - 197.0) < 0.0001) + return 6.38; + else if (abs(abeamtarget_ - 129.0) < 0.0001) + return 5.36; + else if (abs(abeamtarget_ - 63.0) < 0.0001) + return 4.20; + else if (abs(abeamtarget_ - 20.0) < 0.0001) + return 2.992; + else if (abs(abeamtarget_ - 16.0) < 0.0001) + return 2.608; + else + return 1.15 * pow((double)abeamtarget_, 1. / 3.); } } // namespace gen #endif