RS Oph Fermi Analysis Mab Bernardos

From my_wiki
Jump to: navigation, search

Overview[edit]

  • Analysis with Fermitools version 2.0.8 and Fermipy version 1.0.1.

Results from this analysis are available here:

   /fefs/aswg/workspace/maria.bernardos/LSTanalysis/full_analysis/results/fermi

Days coincident with LST Observations[edit]

  • Dates: 59434.4 to 59440.4 (MJD)

Analysis details[edit]

  • Data selection:
   selection_dict = {                                                                                                                                              
       'tmin' : 650112485.0,                                                                                                                                      
       'tmax' : 650630885.0,                                                                                                                                      
       'zmax' : 80.0,                                                                                                                                          
       'evclass' : 128,                                                                                                                                        
       'evtype' : 3,                                                                                                                                           
       'ra' : 267.555,                                                                                                                                         
       'dec' : -6.70778,                                                                                                                                       
       'filter'  : 'DATA_QUAL>0 && LAT_CONFIG==1',                                                                                                             
       'roicut' : 'no',                                                                                                                                        
   }
  • Components:
   components:                                                                                                                                                     
    -  selection : { evtype : 48, emin: 50, emax: 100, zmax: 80  }                                                                                                
    -  selection : { evtype: 3, emin: 100, emax: 300000, zmax: 80  }
  • Binning:
  binning:                                                                                                                                                        
    roiwidth   : 20.0                                                                                                                                             
    binsz      : 0.1                                                                                                                                              
    binsperdec : 8
  • Model:
  model:                                                                                                                                                          
    src_roiwidth : 20.0                                                                                                                                           
    galdiff  : '$FERMI_DIR/refdata/fermi/galdiffuse/gll_iem_v07.fits'                                                                                             
    isodiff  : '$FERMI_DIR/refdata/fermi/galdiffuse/iso_P8R3_SOURCE_V3_v1.txt'                                                                                    
    catalogs :                                                                                                                                                    
       - gll_psc_v26.xml                                                                                 
    sources :                                                                                                                                                     
       - { name :  'RS Oph', ra :  267.555, dec :  -6.7078,                                                                                                        
         SpectrumType : 'LogParabola', norm : !!float 1E-11, Scale : 1000, beta : 0.0,                                                                              
         SpatialModel: 'PointSource' }
  • Sources

The settings and analysis are the same as in the MAGIC paper. SEDs and LCs are produced for sources freed within 4 deg radius of RS Oph going down to 100 MeV. Also, different SEDs and LCs allowing the nearby AGN (4FGL J1745.4-0753) normalization, index, and norm and index to be fixed are compared with the previous case. There is no change in the SEDs and LCs for between each case.


SED[edit]

Nova normal rs oph sed.png

Full observation period daily analysis[edit]

  • Dates from 59431.4 to 59461.4

Analysis details[edit]

  • Data selection:
   selection_dict = {                                                                                                                                              
       'tmin' : 649853285,                                                                                                                                      
       'tmax' : 652445285,                                                                                                                                      
       'zmax' : 90.0,                                                                                                                                          
       'evclass' : 128,                                                                                                                                        
       'evtype' : 3,                                                                                                                                           
       'ra' : 267.555,                                                                                                                                         
       'dec' : -6.70778,                                                                                                                                       
       'emin' : 100.0,                                                                                                                                         
       'emax' : 300000,                                                                                                                                        
       'filter'  : 'DATA_QUAL>0 && LAT_CONFIG==1',                                                                                                             
       'roicut' : 'no',                                                                                                                                        
   }
  • Binning
  binning:                                                                                                                                                        
    roiwidth   : 15.0                                                                                                                                             
    binsz      : 0.1                                                                                                                                              
    binsperdec : 8
  • Model:
  model:                                                                                                                                                          
    src_roiwidth : 20.0                                                                                                                                           
    galdiff  : '$FERMI_DIR/refdata/fermi/galdiffuse/gll_iem_v07.fits'                                                                                             
    isodiff  : '$FERMI_DIR/refdata/fermi/galdiffuse/iso_P8R3_SOURCE_V3_v1.txt'                                                                                    
    catalogs :                                                                                                                                                    
       - gll_psc_v26.xml                                                                                 
    sources :                                                                                                                                                     
       - { name :  'RS Oph', ra :  267.555, dec :  -6.7078,                                                                                                        
         SpectrumType : 'LogParabola', norm : !!float 1E-11, Scale : 1000, beta : 0.0,                                                                              
         SpatialModel: 'PointSource' }

Results Daily[edit]

   #DATE (MJD) TS INTEGRAL_FLUX / 1e-7 FLUX_ERR ALPHA ALPHA_ERR BETA BETA_ERR
   59431.4 -0.000 10.988 0.000 2.000 0.000 0.000 0.000
   59432.4 -0.000 10.743 0.000 2.000 0.000 0.000 0.000
   59433.4 -0.000 10.707 0.000 2.000 0.000 0.000 0.000
   59434.4 189.184 18.256 3.085 2.143 0.130 0.064 0.079
   59435.4 1006.853 46.389 3.904 1.970 0.074 0.201 0.052
   59436.4 570.825 39.466 3.897 2.118 0.097 0.170 0.064
   59437.4 444.719 28.093 3.411 1.965 0.096 0.173 0.066
   59438.4 215.113 17.761 3.081 2.087 0.154 0.229 0.111
   59439.4 181.892 12.745 2.539 1.979 0.157 0.219 0.112
   59440.4 94.723 7.325 2.348 1.632 0.210 0.154 0.114
   59441.4 96.776 11.765 3.086 1.973 0.181 0.164 0.128
   59442.4 60.010 4.468 1.633 2.031 0.416 0.930 0.529
   59443.4 81.262 5.491 1.870 1.587 0.282 0.396 0.209
   59444.4 2.209 11.374 0.000 2.000 0.000 0.000 0.000
   59445.4 22.895 8.776 3.856 2.672 0.689 0.242 0.384
   59446.4 15.754 11.820 0.000 2.000 0.000 0.000 0.000
   59447.4 16.549 11.533 0.000 2.000 0.000 0.000 0.000
   59448.4 27.901 11.887 0.000 2.000 0.000 0.000 0.000
   59449.4 2.576 11.429 0.000 2.000 0.000 0.000 0.000
   59450.4 27.113 3.937 2.312 1.931 0.351 0.158 0.229
   59451.4 0.628 11.214 0.000 2.000 0.000 0.000 0.000
   59452.4 2.304 11.073 0.000 2.000 0.000 0.000 0.000
   59453.4 13.250 11.242 0.000 2.000 0.000 0.000 0.000
   59454.4 20.874 5.946 3.020 2.445 0.476 0.163 0.275
   59455.4 4.612 11.530 0.000 2.000 0.000 0.000 0.000
   59456.4 7.357 11.760 0.000 2.000 0.000 0.000 0.000
   59457.4 0.972 10.989 0.000 2.000 0.000 0.000 0.000
   59458.4 6.207 11.656 0.000 2.000 0.000 0.000 0.000
   59459.4 5.886 11.591 0.000 2.000 0.000 0.000 0.000