diff --git a/scripts/cta_kb_bkgrates.py b/scripts/cta_kb_bkgrates.py index cd764df..bea8e07 100755 --- a/scripts/cta_kb_bkgrates.py +++ b/scripts/cta_kb_bkgrates.py @@ -61,17 +61,18 @@ def make_file_function(irfname, filename): # Compute energy (in MeV) energy = math.pow(10.0, logE)*1.0e6 + emin = math.pow(10.0, logE-0.1)*1.0e6 + emax = math.pow(10.0, logE+0.1)*1.0e6 + ewidth = emax-emin - # Compute solid angle of r68 region in sr - omega = 2.0 * math.pi * (1.0 - math.cos(math.radians(r68))) + # Compute solid angle of r80 region in sr + omega = 2.0 * math.pi * (1.0 - math.cos(math.radians(r80))) # Compute background rate per steradian and MeV - bkg_rate = bgd/omega/energy + bkg_rate = bgd/omega/ewidth # Write file.write(str(energy)+" "+str(bkg_rate)+"\n") - - print logE, energy, bgd, omega, bkg_rate # Close files irf.close()