real-event-analysis.md

@@ -14,8 +14,65 @@ _This page is evolving and new changes need to be added_
                  fix-dQuadMonA = 0.0 
 
 - We may also repeat the same run for a BBH case. This is to make sure that the pipeline works well and gives similar results compared to the PE Prod run. 
-- The `webdir` paths is expected to have a given format (see the script given below)
+- The `webdir` path is expected to have a given format (see the script given below)
 - All runs (except for a particular reason) should be using the same installation
+- Once the job submission is done and the results are available we need to do the following steps for post-processing:
+
+  - Step 1: creating a PE Summary page:
+
+             import glob
+              import os
+              events = ['SXXXXX']
+              outdirBase = '/home/YourNAME/public_html/LSC/lalinference/ObervationRUN/SummaryPages/SXXXXX/C01[or any specification]/ks_sym'
+              paths = ['/home/YourNAME/ObervationRUN/SXXXXX/C01[or any specification]/ks_sym']
+
+              print(events)
+              print(paths)
+              
+              for (event,path) in zip(events,paths):
+                  f =  open(glob.glob(os.path.join(path,'lal*sh'))[0])
+                  content = f.read()
+                  webdir =  os.path.join(outdirBase,event,'ks_sym')
+                  samples = glob.glob(os.path.join(path, 'posterior_samples/pos*.hdf5'))[0]
+                  commandline = 'summarypage'
+                  commandline = '{} --webdir {}'.format(commandline, webdir)
+                  commandline = '{} --samples {}'.format(commandline, samples)
+                  commandline = '{} --psd'.format(commandline)
+
+                  if len(content.split('--L1-psd ')) > 1:
+                         psd = content.split('--L1-psd ')[1].split(' ')[0]
+                         commandline = '{} L1:{}'.format(commandline, psd)
+
+                  if len(content.split('--H1-psd ')) > 1:
+                         psd = content.split('--H1-psd ')[1].split(' ')[0]
+                         commandline = '{} H1:{}'.format(commandline, psd)
+
+                  if len(content.split('--V1-psd ')) > 1:
+                         psd = content.split('--V1-psd ')[1].split(' ')[0]
+                  commandline = '{} V1:{}'.format(commandline, psd)
+
+                  commandline = '{} --calibration'.format(commandline)
+
+                  if len(content.split('--L1-spcal-envelope ')) > 1:
+                         cal = content.split('--L1-spcal-envelope ')[1].split(' ')[0]
+                         commandline = '{} L1:{}'.format(commandline, cal)
+
+                  if len(content.split('--H1-spcal-envelope ')) > 1:
+                         cal = content.split('--H1-spcal-envelope ')[1].split(' ')[0]
+                         commandline = '{} H1:{}'.format(commandline, cal)
+
+                  if len(content.split('--V1-spcal-envelope ')) > 1:
+                         cal = content.split('--V1-spcal-envelope ')[1].split(' ')[0]
+                         commandline = '{} V1:{}'.format(commandline, cal)
+
+                 commandline = '{} --config {}'.format(commandline, os.path.join(path, 'config.ini'))
+                 commandline = '{} {}'.format(commandline, '--no_ligo_skymap --labels pv3hm_c01_ks_sym')
+                 posterior = glob.glob('{}/posterior_samples/*.hdf5'.format(path))
+                 if len(posterior) > 0:
+                        print( commandline)
+                        os.system(commandline)
+
+
 
 - Most of the above-mentioned steps would automatically follow for an injection study