def dRdEi(self, E, mx, i):
R = self.exposure * self.frac[i] * DMU.dRdE(
E, self.N_p[i], self.N_n[i], mx,
[1.0 / self.N_p[i], 0.0, 0.0, 0.0])
if (R < 1e-30):
return 1e-30
else:
return R
def dRdE(self, E, mx, l):
R = 0
for i in range(self.N_iso):
R += self.frac[i] * DMU.dRdE(E, self.N_p[i], self.N_n[i], mx, l)
#print R
if (R < 1e-30):
return 1e-30
else:
return R * self.exposure
