Tutorial GRB Analisi prompt emission GRB

Transcription

1 Tutorial GRB Analisi prompt emission GRB Sara Cutini

2 Posizione Swift o LAT

3 Trigger GBM

4 Scaricate i dati dal Fermi Science Support Center: Link per LAT: Link per GBM: ftp://legacy.gsfc.nasa.gov/fermi/data/gbm/bursts/

5 Spacecraft: L*SC00.FIT Fotoni: L*PH00.FIT Dati LAT

6 Dati GBM

7 fkeyprint glg_trigdat_all_bn _v01.fit DET_MASK= b0 sta dalla parte dei n0-n5 b1 sta dalla parte dei n6-n11 hanno triggerato il 3,6 e 8 -> il segnale migliore è in n7, n8 + b1 e b0 due file per ogni detector: glg_cspec_n7_bn _v03.rsp glg_tte_n7_bn _v00.fit

8 Curva di luce (GBM) gtbin ~/tutorial_fermi]$ gtbin This is gtbin version ScienceTools-v9r15p2-fssc Type of output file (CCUBE CMAP LC PHA1 PHA2) [] LC Event data file name[] glg_tte_n8_bn _v00.fit Output file name[] glg_tte_n8_bn _v00.lc Spacecraft data file name[] none Algorithm for defining time bins (FILE LIN SNR) LIN Start value for first time bin in MET[] INDEF Stop value for last time bin in MET[] INDEF Width of linearly uniform time bins in MET[] 0.05

9 Curva di luce (LAT) LAT selezione in tempo e spazio: [cutini@sagrantino ~/tutorial_fermi]$ gtselect evclsmin=1 evclsmax=3 Input FT1 file[] L E0D2F37E92_PH00.fits Output FT1 file[] grb fits RA for new search center (degrees) (0:360) [] Dec for new search center (degrees) (-90:90) [] radius of new search region (degrees) (0:180) [] 10 start time (MET in s) (0:) [] end time (MET in s) (0:) [] lower energy limit (MeV) (0:) [] 100 upper energy limit (MeV) (0:) [] maximum zenith angle value (degrees) (0:180) [] 105 Done. curva di luce del LAT cutini@sagrantino ~/tutorial_fermi]$ gtbin This is gtbin version ScienceTools-v9r15p2-fssc Type of output file (CCUBE CMAP LC PHA1 PHA2) [LC] Event data file name[] grb fits Output file name[] grb lc Spacecraft data file name[] L E0D2F37E75_SC00.fits Algorithm for defining time bins (FILE LIN SNR) [] LIN Start value for first time bin in MET[] Stop value for last time bin in MET[] Width of linearly uniform time bins in MET[0.05]0.05

10

11 Mappa conteggi (LAT) gtbin This is gtbin version ScienceTools-v9r15p2-fssc Type of output file (CCUBE CMAP LC PHA1 PHA2) [] CMAP Event data file name[] grb fits Output file name[] grb cmap Spacecraft data file name[] L E0D2F37E75_SC00.fits Size of the X axis in pixels[]120 Size of the Y axis in pixels[]120 Image scale (in degrees/pixel)[] 0.25 Coordinate system (CEL - celestial, GAL -galactic) (CEL GAL) []CEL First coordinate of image center in degrees (RA or galactic l)[] Second coordinate of image center in degrees (DEC or galactic b)[]

12

13 Analisi spettarale Joint 1) Costruire il pha2 per i file del GBM Mytime txt > bkg -21s - -1s > burst 0.5s - 1s cutini@sagrantino ~/tutorial_fermi]$ gtbindef This is gtbindef version ScienceTools-v9r15p2-fssc Type of bins, E -- energy or T -- time[] T File containing the ascii bin definitions[] time.txt Output file name[] time.fit

14 gtbin This is gtbin version ScienceTools-v9r15p2-fssc Type of output file (CCUBE CMAP LC PHA1 PHA2) [] PHA2 Event data file name[] glg_tte_n8_bn _v00.fit Output file name[] glg_tte_n8_bn _v00.pha Spacecraft data file name[] none Algorithm for defining time bins (FILE LIN SNR) [] FILE Name of the file containing the time bin definition[] time.fit 1) Costruire il pha per il file del LAT-> stessa exposure del GBM- ritagliare il file LAT con gtselect gtbin This is gtbin version ScienceTools-v9r15p2-fssc Type of output file (CCUBE CMAP LC PHA1 PHA2) [] PHA1 Event data file name[] grb090510_0.5-1.fits Output file name[] grb090510_0.5-1.pha Spacecraft data file name[] L E0D2F37E75_SC00.fits Algorithm for defining energy bins (FILE LIN LOG) LOG Start value for first energy bin in MeV[] 100 Stop value for last energy bin in MeV[] Number of logarithmically uniform energy bins[] 30

15 3) Costruire la matrice di risposta del LAT ~/tutorial_fermi]$ gtrspgen This is gtrspgen version ScienceTools-v9r15p2-fssc Response calculation method (GRB, PS) (GRB PS) GRB Spectrum file name[]grb090510_0.5-1.pha Spacecraft data file name[]l e0d2f37e75_sc00.fits Output file name[]grb090510_0.5-1.rsp Time of GRB (s)[] Response function to use, Handoff DC2 DC2A DC2FA DC2BA DC2FB etc P6_V3_TRANSIENT Algorithm for defining true energy bins (FILE LIN LOG) LOG Start value for first energy bin in MeV[] 100 Stop value for last energy bin in MeV[] Number of logarithmically uniform energy bins[]100 Run xspec

16

17 Macro per xspec: data 1 glg_tte_b1_bn _v00.pha{2} data 2 glg_tte_n7_bn _v00.pha{2} data 3 glg_tte_n8_bn _v00.pha{2} data 4 grb090510_0.5-1.pha data 5 glg_tte_b0_bn _v00.pha{2} resp 1 glg_cspec_b1_bn _v03.rsp resp 2 glg_cspec_n7_bn _v03.rsp resp 3 glg_cspec_n8_bn _v03.rsp resp 5 glg_cspec_b0_bn _v03.rsp resp 4 grb090510_0.5-1.rsp back 1 glg_tte_b1_bn _v00.pha{1} back 2 glg_tte_n7_bn _v00.pha{1} back 3 glg_tte_n8_bn _v00.pha{1} back 5 glg_tte_b0_bn _v00.pha{1} ign 1:1-5, 100-** ign 2:1-5, 100-** ign 3:1-5, 100-** ign 5:1-5, 100-** ign 4: ** ign bad setpl ene setpl area setpl rebin setpl rebin setpl rebin setpl rebin 5 5 3

18 Afterglow analisi->point Source gtburstfit-> fit delle curve di luce del LAT con la seguente funzione+ bkg costante : A * exp(tau1 / (t - time0) + (t - time0) / tau2) cutini@sagrantino ~/tutorial_fermi]$ gtburstfit This is gtburstfit version N/A Name of file containing events[grb lc] How will initial values for fit be determined (AUTO MAN CALC) [MAN] AUTO After initial guess but before any fitting, reduced chi square is Parameters are: Peak 1: Amplitude = Time0 = Tau1 = Tau2 = Background = After fit, reduced chi square is Parameters are: Peak 1: Amplitude = Time0 = Tau1 = Tau2 = Background = Bayesian Blocks computed for this data set are: Interval Average Counts [ , ] [ , ] [ , ] [ , ]