Commit 6a82a26a authored by Jean-Gregoire Ducoin's avatar Jean-Gregoire Ducoin

adding the calculation for the X_pipeline distance exclusion

parent 24e6aed1
Pipeline #50446 failed with stages
in 7 minutes and 27 seconds
......@@ -244,7 +244,7 @@ def compose_emcoinc(gracedb_id, authors=(),
'/home/grb.exttrig/Online/O2/processor/runs/X/GRBE298784/'
'auto_web/E298784_online_openbox_summary.txt'
) # Michal example
path_output_CBC = 'need to be defined'
#path_output_CBC = 'need to be defined'
if os.path.exists(path_output_X):
output_X = open(path_output_X, 'r')
......@@ -253,6 +253,16 @@ def compose_emcoinc(gracedb_id, authors=(),
if infos_X[2] == 'sgc150q9':
UL50 = infos_X[3]
# need to build a function to convert UL50 to exclusion distance ?
#definition of physical constantes
Msolar = 1.9891e30 #solar mass in kg
G = 6.673*1e-11 #gravitational constant in m^3/(kg*s^2)
c = 299792458 #speed of light in m/s
pc = 3.0856775807*1e16 #parsec in m
Mpc = pc*1e6 #mega-parsec in m
refEnergy = 1e-2*Msolar*c^2 #energy reference
w = 150 #pulsation
X_exclusion_distance = np.sqrt(10*G*refEnergy/(c**3*UL50**2*(2*pi*w)**2))/Mpc
if client is None:
client = rest.GraceDb(service)
......@@ -280,7 +290,7 @@ def compose_emcoinc(gracedb_id, authors=(),
CITE_GCN='(Klebesadel, et al., GCN Circ. XXXXX)',
NSNS_Mpc='XXX',
NSBH_Mpc='XXX',
Burst_Mpc=12,
Burst_Mpc=X_exclusion_distance,
CBC_FAR=CBC_FAR,
Burst_FAR=Burst_FAR)
subject = (
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment