@@ -44,7 +44,8 @@ double AtELossCATIMA::GetRange(double energyIni, double energyFin) const
4444 double remainingEnergy{energyIni};
4545 double range{0 };
4646 while (remainingEnergy > energyFin) {
47- catima::Result result = catima::calculate (*fProjectile , *fMaterial , remainingEnergy / fProjectileMassAmu , fConfig );
47+ catima::Result result =
48+ catima::calculate (*fProjectile , *fMaterial , remainingEnergy / fProjectileMassAmu , fConfig );
4849 double dEdx = result.dEdxi * fDensity ;
4950 double DE = dEdx * fRangeStepSize / 10 .;
5051
@@ -65,7 +66,8 @@ double AtELossCATIMA::GetEnergy(double energyIni, double distance) const
6566 double remainingEnergy{energyIni};
6667 double range{0 };
6768 while (range < distance) {
68- catima::Result result = catima::calculate (*fProjectile , *fMaterial , remainingEnergy / fProjectileMassAmu , fConfig );
69+ catima::Result result =
70+ catima::calculate (*fProjectile , *fMaterial , remainingEnergy / fProjectileMassAmu , fConfig );
6971 double dEdx = result.dEdxi * fDensity ;
7072 double DE{};
7173
@@ -142,8 +144,9 @@ AtELossCATIMA::GetBraggCurve(double energy, double rangeStepSize, double totalFr
142144 double range{};
143145 while (remainingEnergy / energy > totalFractionELoss) {
144146
145- catima::Result result = catima::calculate (*fProjectile , *fMaterial , remainingEnergy / fProjectileMassAmu , fConfig );
146- double dEdx = result.dEdxi * fDensity ; // MeV/cm
147+ catima::Result result =
148+ catima::calculate (*fProjectile , *fMaterial , remainingEnergy / fProjectileMassAmu , fConfig );
149+ double dEdx = result.dEdxi * fDensity ; // MeV/cm
147150 braggCurve.push_back (std::make_pair (dEdx / 10.0 , range)); // store in MeV/mm
148151
149152 double DE = dEdx * rangeStepSize / 10 .;
0 commit comments