From 9862cebda1b6df5f2cbfc104f10f6a1c8833d93c Mon Sep 17 00:00:00 2001
From: knodlseder <>
Date: Tue, 2 Aug 2011 00:16:33 +0000
Subject: [PATCH] use r80 instead of r68 for background rate estimation

---
 scripts/cta_kb_bkgrates.py | 11 ++++++-----
 1 file changed, 6 insertions(+), 5 deletions(-)

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()
--
libgit2 0.21.4