All current external libraries are now available on RHEL7.
Testing ScienceTools-11-03-00
test_evtbin
This is test_evtbin version ScienceTools-11-03-00 Expected: failed to create a BinConfig before prototypes were loaded: BinConfig::create was unable to find a configuration for mission "GLAST", instrument "LAT" while processing file "/nfs/farm/g/glast/u52/rhel7-test/ScienceTools-scons/data/evtbin/ft1tiny.fits" Error encountered at /afs/slac.stanford.edu/g/glast/ground/GLAST_EXT/redhat7-x86_64-64bit-gcc48/Healpix_3.30/src/cxx/Healpix_cxx/healpix_base.cc, line 701 (function void T_Healpix_Base<I>::Set(int, Healpix_Ordering_Scheme) [with I = int]) bad order Expected: failed to create a Healpix Binner with order -1 : bad order Expected: failed to create a Healpix Binner with order 13 : Order needs to be positive and <=12 test_evtbin: WARNING: No spacecraft file: EXPOSURE keyword will be set equal to ontime. test_evtbin: WARNING: No spacecraft file: EXPOSURE keyword will be set equal to ontime.
test_dataSubselector
test_Likelihood
test_astro
test_tip
test_st_facilities
test_sane
bash-4.2$ test_sane time -p gtobssim infile=xmlFiles.txt srclist=source_names.txt scfile=none sctable="SC_DATA" evroot="test" evtable="EVENTS" simtime=86400.0 ltfrac=0.9 tstart=86400.0 nevents=no maxtime=315500000.0 startdate="2001-01-01 00:00:00" offset=0 rockangle="INDEF" use_ac=yes ra=90.0 dec=20.0 radius=30.0 emin=100.0 emax=300000.0 edisp=yes irfs="DC1A" evtype="none" area=1.0 maxrows=1000000 seed=479153 chatter=2 clobber=yes debug=no gui=no mode="ql" added source "anticenter-32mev" added source "Galactic_diffuse" Generating events for a simulation time of 86400 seconds.... Done. real 52.48 user 46.53 sys 0.18 time -p gtmktime scfile=orbSim_scData_0000.fits sctable="SC_DATA" filter="IN_SAA!=T" roicut=yes evfile=test_events_0000.fits evtable="EVENTS" outfile="test_events.fits" apply_filter=yes overwrite=no header_obstimes=yes tstart=0.0 tstop=0.0 gtifile="default" chatter=2 clobber=yes debug=no gui=no mode="ql" DSS keywords required for ROI-based zenith angle cut are not present in the FT1 file. real 1.36 user 0.15 sys 0.06 time -p gtvcut infile=test_events.fits table="EVENTS" suppress_gtis=yes chatter=2 debug=no gui=no mode="ql" DSTYP1: ENERGY DSUNI1: MeV DSVAL1: 100:300000 DSTYP2: POS(RA,DEC) DSUNI2: deg DSVAL2: CIRCLE(90,20,30) DSTYP3: TIME DSUNI3: s DSVAL3: TABLE DSREF3: :GTI GTIs: (suppressed) real 0.72 user 0.11 sys 0.04 time -p gtselect infile=test_events.fits outfile=filtered_events_0000.fits ra=90.0 dec=20.0 rad=20.0 tmin=0.0 tmax=0.0 emin=32.0 emax=200000.0 zmin=0.0 zmax=180.0 evclass="INDEF" evtype="INDEF" convtype=-1 phasemin=0.0 phasemax=1.0 evtable="EVENTS" chatter=2 clobber=yes debug=no gui=no mode="ql" Done. real 1.03 user 0.16 sys 0.05 time -p gtvcut infile=filtered_events_0000.fits table="EVENTS" suppress_gtis=yes chatter=2 debug=no gui=no mode="ql" DSTYP1: POS(RA,DEC) DSUNI1: deg DSVAL1: CIRCLE(90,20,20) DSTYP2: TIME DSUNI2: s DSVAL2: TABLE DSREF2: :GTI GTIs: (suppressed) DSTYP3: ENERGY DSUNI3: MeV DSVAL3: 100:200000 real 0.61 user 0.12 sys 0.04 time -p gtselect infile=test_events.fits outfile=filtered1.fits ra=90.0 dec=20.0 rad=20.0 tmin=86400.0 tmax=129600.0 emin=30.0 emax=300000.0 zmin=0.0 zmax=180.0 evclass="INDEF" evtype="INDEF" convtype=-1 phasemin=0.0 phasemax=1.0 evtable="EVENTS" chatter=2 clobber=yes debug=no gui=no mode="ql" Done. real 1.01 user 0.17 sys 0.05 time -p gtvcut infile=filtered1.fits table="EVENTS" suppress_gtis=yes chatter=2 debug=no gui=no mode="ql" DSTYP1: POS(RA,DEC) DSUNI1: deg DSVAL1: CIRCLE(90,20,20) DSTYP2: TIME DSUNI2: s DSVAL2: TABLE DSREF2: :GTI GTIs: (suppressed) DSTYP3: ENERGY DSUNI3: MeV DSVAL3: 100:300000 real 0.61 user 0.12 sys 0.03 time -p gtltcube evfile="filtered1.fits" evtable="EVENTS" scfile=orbSim_scData_0000.fits sctable="SC_DATA" outfile=expcube1.fits dcostheta=0.05 binsz=1.0 phibins=0 tmin=0.0 tmax=0.0 file_version="1" zmin=0.0 zmax=180.0 chatter=4 clobber=yes debug=no gui=no mode="ql" This is gtltcube version ScienceTools-11-03-00 applying filter: (START >= 86340) && (STOP <= 129450) Working on file orbSim_scData_0000.fits read 1435 rows .....................! real 10.66 user 8.58 sys 0.77 time -p gtselect infile=test_events.fits outfile=filtered2.fits ra=90.0 dec=20.0 rad=20.0 tmin=129600.0 tmax=172800.0 emin=30.0 emax=300000.0 zmin=0.0 zmax=180.0 evclass="INDEF" evtype="INDEF" convtype=-1 phasemin=0.0 phasemax=1.0 evtable="EVENTS" chatter=2 clobber=yes debug=no gui=no mode="ql" Done. real 0.98 user 0.16 sys 0.05 time -p gtvcut infile=filtered2.fits table="EVENTS" suppress_gtis=yes chatter=2 debug=no gui=no mode="ql" DSTYP1: POS(RA,DEC) DSUNI1: deg DSVAL1: CIRCLE(90,20,20) DSTYP2: TIME DSUNI2: s DSVAL2: TABLE DSREF2: :GTI GTIs: (suppressed) DSTYP3: ENERGY DSUNI3: MeV DSVAL3: 100:300000 real 0.59 user 0.11 sys 0.03 time -p gtltcube evfile="filtered2.fits" evtable="EVENTS" scfile=orbSim_scData_0000.fits sctable="SC_DATA" outfile=expcube2.fits dcostheta=0.05 binsz=1.0 phibins=0 tmin=0.0 tmax=0.0 file_version="1" zmin=0.0 zmax=180.0 chatter=4 clobber=yes debug=no gui=no mode="ql" This is gtltcube version ScienceTools-11-03-00 applying filter: (START >= 130590) && (STOP <= 172860) Working on file orbSim_scData_0000.fits read 1408 rows .....................! real 10.41 user 8.52 sys 0.76 time -p gtltsum infile1=expcube1.fits infile2=expcube2.fits table="Exposure" table2="WEIGHTED_EXPOSURE" outfile=expcube_1_day.fits chatter=4 clobber=yes debug=no gui=no mode="ql" This is gtltsum version ScienceTools-11-03-00 real 1.90 user 0.55 sys 0.53 time -p gtexpmap evfile=filtered_events_0000.fits evtable="EVENTS" scfile=orbSim_scData_0000.fits sctable="SC_DATA" expcube=expcube_1_day.fits outfile=expMap.fits irfs="DC1A" evtype="INDEF" srcrad=30.0 nlong=120 nlat=120 nenergies=20 submap=no nlongmin=0 nlongmax=0 nlatmin=0 nlatmax=0 chatter=4 clobber=yes debug=no gui=no mode="ql" This is gtexpmap version ScienceTools-11-03-00 The exposure maps generated by this tool are meant to be used for *unbinned* likelihood analysis only. Do not use them for binned analyses. ResponseFunctions::load: IRF used: DC1A event_types: 0 1 Computing the ExposureMap using expcube_1_day.fits ....................! real 57.30 user 56.26 sys 0.17 time -p gtvcut infile=expMap.fits table="PRIMARY" suppress_gtis=yes chatter=2 debug=no gui=no mode="ql" DSTYP1: POS(RA,DEC) DSUNI1: deg DSVAL1: CIRCLE(90,20,20) DSTYP2: ENERGY DSUNI2: MeV DSVAL2: 100:200000 DSTYP3: TIME DSUNI3: s DSVAL3: TABLE DSREF3: :GTI GTIs: (suppressed) DSTYP4: IRF_VERSION DSUNI4: DIMENSIONLESS DSVAL4: DC1A real 0.61 user 0.11 sys 0.04 time -p gtdiffrsp evfile=filtered_events_0000.fits evtable="EVENTS" scfile=orbSim_scData_0000.fits sctable="SC_DATA" srcmdl=/nfs/farm/g/glast/u52/rhel7-test/ScienceTools-scons/sane/data/srcModel_egretdiffuse.xml irfs="DC1A" evclsmin=0 evclass="INDEF" evtype="INDEF" convert=no chatter=4 clobber=yes debug=no gui=no mode="ql" This is gtdiffrsp version ScienceTools-11-03-00 ResponseFunctions::load: IRF used: DC1A event_types: 0 1 Creating source named EGRET Diffuse MapBase::readFitsFile: creating WcsMap2 object Creating source named Crab Creating source named PKS 0528+134 Creating source named Geminga adding source EGRET Diffuse Working on... filtered_events_0000.fitsEventContainer::getEvents: Out of 3500 events in file filtered_events_0000.fits, 3500 were accepted, and 0 were rejected. ....................! real 234.80 user 233.63 sys 0.10 time -p gtvcut infile=filtered_events_0000.fits table="EVENTS" suppress_gtis=yes chatter=2 debug=no gui=no mode="ql" DSTYP1: POS(RA,DEC) DSUNI1: deg DSVAL1: CIRCLE(90,20,20) DSTYP2: TIME DSUNI2: s DSVAL2: TABLE DSREF2: :GTI GTIs: (suppressed) DSTYP3: ENERGY DSUNI3: MeV DSVAL3: 100:200000 real 0.60 user 0.12 sys 0.03 time -p gtlike irfs="DC1A" expcube=expcube_1_day.fits srcmdl=/nfs/farm/g/glast/u52/rhel7-test/ScienceTools-scons/sane/data/srcModel_egretdiffuse.xml sfile=none check_fit=yes results=results.dat specfile=counts_spectra.fits statistic="UNBINNED" optimizer="MINUIT" ftol=0.0001 toltype="ABS" tsmin=no save=yes refit=no evfile=filtered_events_0000.fits evtable="EVENTS" scfile=orbSim_scData_0000.fits sctable="SC_DATA" expmap=expMap.fits plot=no cmap=none bexpmap=none wmap=none psfcorr=yes phased_expmap=none chatter=3 clobber=yes debug=no gui=no mode="ql" This is gtlike version ScienceTools-11-03-00 ResponseFunctions::load: IRF used: DC1A event_types: 0 1 EventContainer::getEvents: Out of 3500 events in file filtered_events_0000.fits, 3500 were accepted, and 0 were rejected. Creating source named EGRET Diffuse Creating source named Crab Creating source named PKS 0528+134 Creating source named Geminga ....................! ********** ** 1 **SET PRINT 1.000 ********** ********** ** 2 **SET NOWARN ********** PARAMETER DEFINITIONS: NO. NAME VALUE STEP SIZE LIMITS 1 'Integral ' 1.0000 1.0000 .10000E-04 1000.0 2 'Index ' -2.0000 1.0000 -5.0000 -1.0000 3 'Prefactor ' 11.000 1.0000 .10000E-02 1000.0 4 'Integral ' 1.0000 1.0000 .10000E-04 1000.0 5 'Index ' -2.0000 1.0000 -5.0000 -1.0000 6 'Integral ' 1.0000 1.0000 .10000E-04 1000.0 7 'Index ' -2.0000 1.0000 -5.0000 -1.0000 ********** ** 3 **SET ERR .5000 ********** ********** ** 4 **SET GRAD 1.000 ********** ********** ** 5 **MINIMIZE 1500. .2000 ********** FIRST CALL TO USER FUNCTION AT NEW START POINT, WITH IFLAG=4. START MIGRAD MINIMIZATION. STRATEGY 1. CONVERGENCE WHEN EDM .LT. .20E-03 FCN= 32812.48 FROM MIGRAD STATUS=INITIATE 1 CALLS 2 TOTAL EDM= unknown STRATEGY= 1 NO ERROR MATRIX EXT PARAMETER CURRENT GUESS STEP FIRST NO. NAME VALUE ERROR SIZE DERIVATIVE 1 Integral 1.00000 1.0000 .43760E-01 -1781.5 2 Index -2.0000 1.0000 .78442 -241.58 3 Prefactor 11.000 1.0000 .95978E-02 -4574.3 4 Integral 1.00000 1.0000 .43760E-01 -4885.3 5 Index -2.0000 1.0000 .78442 -844.89 6 Integral 1.00000 1.0000 .43760E-01 -559.93 7 Index -2.0000 1.0000 .78442 -45.940 ERR DEF= .500 MIGRAD MINIMIZATION HAS CONVERGED. MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX. FCN= 32156.41 FROM MIGRAD STATUS=CONVERGED 161 CALLS 162 TOTAL EDM= .28E-05 STRATEGY= 1 ERR MATRIX NOT POS-DEF EXT PARAMETER APPROXIMATE STEP FIRST NO. NAME VALUE ERROR SIZE DERIVATIVE 1 Integral 12.155 15.229 .37110E-01 .19677E-01 2 Index -2.0840 .52073 .78442E-01 .12582E-01 3 Prefactor 11.100 .26082 .84096E-02 .24172 4 Integral 11.390 .81855 .24345E-01 -.58289E-01 5 Index -1.7038 .27825E-01 .57956E-01 .54878E-01 6 Integral 10.813 4.0665 .43760E-01 .20105E-01 7 Index -2.5288 .17312 .14234 -.91183E-02 ERR DEF= .500 EXTERNAL ERROR MATRIX. NDIM= 100 NPAR= 7 ERR DEF= .500 .233E+03 -.806E+01 -.244E+01 .329E+00 -.538E-03 .143E+02 .471E+00 -.806E+01 .279E+00 .835E-01 -.112E-01 .197E-04 -.493E+00 -.162E-01 -.244E+01 .835E-01 .680E-01 -.186E-01 -.155E-03 -.721E-01 -.302E-02 .329E+00 -.112E-01 -.186E-01 .670E+00 .155E-01 .524E-02 .350E-03 -.538E-03 .197E-04 -.155E-03 .155E-01 .774E-03 -.766E-04 -.612E-06 .143E+02 -.493E+00 -.721E-01 .524E-02 -.766E-04 .165E+02 .656E+00 .471E+00 -.162E-01 -.302E-02 .350E-03 -.612E-06 .656E+00 .301E-01 ERR MATRIX NOT POS-DEF PARAMETER CORRELATION COEFFICIENTS NO. GLOBAL 1 2 3 4 5 6 7 1 .99812 1.000 -.998 -.611 .026 -.001 .230 .178 2 .99810 -.998 1.000 .606 -.026 .001 -.229 -.177 3 .62795 -.611 .606 1.000 -.087 -.021 -.068 -.067 4 .68679 .026 -.026 -.087 1.000 .683 .002 .002 5 .68367 -.001 .001 -.021 .683 1.000 -.001 .000 6 .93314 .230 -.229 -.068 .002 -.001 1.000 .930 7 .93116 .178 -.177 -.067 .002 .000 .930 1.000 ERR MATRIX NOT POS-DEF Final values: Integral = 12.155 Index = -2.08405 Prefactor = 11.0996 Integral = 11.3898 Index = -1.70379 Integral = 10.8132 Index = -2.52883 ********** ** 6 **HESSE ********** COVARIANCE MATRIX CALCULATED SUCCESSFULLY FCN= 32156.41 FROM HESSE STATUS=OK 50 CALLS 212 TOTAL EDM= .14E-05 STRATEGY= 1 ERROR MATRIX ACCURATE EXT PARAMETER INTERNAL INTERNAL NO. NAME VALUE ERROR STEP SIZE VALUE 1 Integral 12.155 2.2395 .16407E-03 -1.3498 2 Index -2.0840 .77225E-01 .35170E-03 6.7589 3 Prefactor 11.100 .23162 .37161E-04 -1.3597 4 Integral 11.390 1.1859 .10754E-03 -1.3569 5 Index -1.7038 .40110E-01 .25571E-03 15.003 6 Integral 10.813 4.1094 .27118E-03 -1.7791 7 Index -2.5288 .17933 .62369E-03 -6.0454 ERR DEF= .500 EXTERNAL ERROR MATRIX. NDIM= 100 NPAR= 7 ERR DEF= .500 .502E+01 -.153E+00 -.800E-01 .370E-01 .710E-03 .330E+00 .114E-01 -.153E+00 .597E-02 .161E-02 -.639E-03 -.114E-04 -.933E-02 -.323E-03 -.800E-01 .161E-02 .536E-01 -.401E-01 -.864E-03 .902E-01 .184E-02 .370E-01 -.639E-03 -.401E-01 .141E+01 .396E-01 -.376E-01 -.476E-03 .710E-03 -.114E-04 -.864E-03 .396E-01 .161E-02 -.673E-03 -.556E-05 .330E+00 -.933E-02 .902E-01 -.376E-01 -.673E-03 .169E+02 .679E+00 .114E-01 -.323E-03 .184E-02 -.476E-03 -.556E-05 .679E+00 .322E-01 PARAMETER CORRELATION COEFFICIENTS NO. GLOBAL 1 2 3 4 5 6 7 1 .88885 1.000 -.885 -.154 .014 .008 .036 .028 2 .88674 -.885 1.000 .090 -.007 -.004 -.029 -.023 3 .28028 -.154 .090 1.000 -.146 -.093 .095 .044 4 .83532 .014 -.007 -.146 1.000 .832 -.008 -.002 5 .83295 .008 -.004 -.093 .832 1.000 -.004 -.001 6 .92140 .036 -.029 .095 -.008 -.004 1.000 .920 7 .92069 .028 -.023 .044 -.002 -.001 .920 1.000 Minuit fit quality: 3 estimated distance: 1.40233e-06 Minuit parameter uncertainties: 1 2.23966 2 0.0772497 3 0.231623 4 1.18588 5 0.0401144 6 4.11047 7 0.179581 ....! Computing TS values for each source (4 total) Photon fluxes are computed for the energy range 100 to 200000 MeV Crab: Integral: 12.155 +/- 2.23966 Index: -2.08405 +/- 0.0772497 LowerLimit: 20 UpperLimit: 200000 Npred: 242.222 ROI distance: 6.3295 TS value: 468.148 Flux: 2.12421e-06 +/- 1.99726e-07 photons/cm^2/s EGRET Diffuse: Prefactor: 11.0996 +/- 0.231623 Index: -2.1 Scale: 100 Npred: 2745.74 Flux: 0.000456278 +/- 9.51586e-06 photons/cm^2/s Geminga: Integral: 11.3898 +/- 1.18588 Index: -1.70379 +/- 0.0401144 LowerLimit: 20 UpperLimit: 200000 Npred: 414.101 ROI distance: 8.3097 TS value: 1723.39 Flux: 3.65839e-06 +/- 5.88984e-07 photons/cm^2/s PKS 0528+134: Integral: 10.8132 +/- 4.11047 Index: -2.52883 +/- 0.179581 LowerLimit: 20 UpperLimit: 200000 Npred: 97.9565 ROI distance: 9.59745 TS value: 99.5038 Flux: 9.24358e-07 +/- 6.05517e-07 photons/cm^2/s WARNING: Fit may be bad in range [213.847, 312.719] (MeV) Total number of observed counts: 3500 Total number of model events: 3500.02 -log(Likelihood): 32156.40502 Elapsed CPU time: 12.63 real 13.74 user 12.68 sys 0.14 Crab Spectrum: PowerLaw2 0 Integral: 1.218e+01 2.241e+00 1.000e-05 1.000e+03 ( 1.000e-06) 1 Index: -2.085e+00 7.720e-02 -5.000e+00 -1.000e+00 ( 1.000e+00) 2 LowerLimit: 2.000e+01 0.000e+00 2.000e+01 2.000e+05 ( 1.000e+00) fixed 3 UpperLimit: 2.000e+05 0.000e+00 2.000e+01 2.000e+05 ( 1.000e+00) fixed EGRET Diffuse Spectrum: PowerLaw 4 Prefactor: 1.110e+01 2.316e-01 1.000e-03 1.000e+03 ( 1.000e-03) 5 Index: -2.100e+00 0.000e+00 -3.500e+00 -1.000e+00 ( 1.000e+00) fixed 6 Scale: 1.000e+02 0.000e+00 5.000e+01 2.000e+02 ( 1.000e+00) fixed Geminga Spectrum: PowerLaw2 7 Integral: 1.139e+01 1.186e+00 1.000e-05 1.000e+03 ( 1.000e-06) 8 Index: -1.704e+00 4.011e-02 -5.000e+00 -1.000e+00 ( 1.000e+00) 9 LowerLimit: 2.000e+01 0.000e+00 2.000e+01 2.000e+05 ( 1.000e+00) fixed 10 UpperLimit: 2.000e+05 0.000e+00 2.000e+01 2.000e+05 ( 1.000e+00) fixed PKS 0528+134 Spectrum: PowerLaw2 11 Integral: 1.081e+01 4.155e+00 1.000e-05 1.000e+03 ( 1.000e-06) 12 Index: -2.529e+00 1.817e-01 -5.000e+00 -1.000e+00 ( 1.000e+00) 13 LowerLimit: 2.000e+01 0.000e+00 2.000e+01 2.000e+05 ( 1.000e+00) fixed 14 UpperLimit: 2.000e+05 0.000e+00 2.000e+01 2.000e+05 ( 1.000e+00) fixed Ts values: Crab 468.141922683 EGRET Diffuse 4.6e+31 Geminga 1723.394982 PKS 0528+134 99.5015141918 Exercise UpperLimits.py Info in <Minuit2>: VariableMetricBuilder: Tolerance is not sufficient, continue the minimization Info in <Minuit2>: Current Edm is : edm = 0.00243187 Info in <Minuit2>: Required Edm is : edmval = 0.002 Info in <Minuit2>: VariableMetricBuilder: Tolerance is not sufficient, continue the minimization Info in <Minuit2>: Current Edm is : edm = 0.00202712 Info in <Minuit2>: Required Edm is : edmval = 0.002 0 12.1759138504 -3.89065280615e-05 2.12580546324e-06 1 12.8534138504 0.0464672806229 2.17250373333e-06 2 13.5309138504 0.171753884413 2.21674873776e-06 3 14.2084138504 0.364490870554 2.25885725792e-06 4 14.8859138504 0.615180435878 2.29892830279e-06 5 15.5634138504 0.915826330674 2.33694677599e-06 6 16.2409138504 1.25964807016 2.37360424633e-06 7 16.6117672717 1.46288143897 2.39282944572e-06 Crab 2.38262422112e-06 Exercise IntegralUpperLimit.py Info in <Minuit2>: VariableMetricBuilder: Tolerance is not sufficient, continue the minimization Info in <Minuit2>: Current Edm is : edm = 0.00230245 Info in <Minuit2>: Required Edm is : edmval = 0.002 Info in <Minuit2>: VariableMetricBuilder: Tolerance is not sufficient, continue the minimization Info in <Minuit2>: Current Edm is : edm = 0.00384981 Info in <Minuit2>: Required Edm is : edmval = 0.002 Info in <Minuit2>: VariableMetricBuilder: Tolerance is not sufficient, continue the minimization Info in <Minuit2>: Current Edm is : edm = 0.0030095 Info in <Minuit2>: Required Edm is : edmval = 0.002 Info in <Minuit2>: VariableMetricBuilder: Tolerance is not sufficient, continue the minimization Info in <Minuit2>: Current Edm is : edm = 0.00204942 Info in <Minuit2>: Required Edm is : edmval = 0.002 Crab 2.18749862044e-06 time -p gtselect infile=test_events_0000.fits outfile=filtered_events.fits ra="INDEF" dec="INDEF" rad="INDEF" tmin="INDEF" tmax="INDEF" emin=30.0 emax=300000.0 zmin=0.0 zmax=180.0 evclass="INDEF" evtype="INDEF" convtype=-1 phasemin=0.0 phasemax=1.0 evtable="EVENTS" chatter=2 clobber=yes debug=no gui=no mode="ql" Done. real 1.03 user 0.16 sys 0.07 time -p gtbin evfile=filtered_events.fits scfile=orbSim_scData_0000.fits outfile=countsMap.fits algorithm="CCUBE" ebinalg="LOG" emin=100.0 emax=200000.0 enumbins=30 ebinfile=NONE tbinalg="LIN" tbinfile=NONE nxpix=160 nypix=160 binsz=0.25 coordsys="CEL" xref=90.0 yref=20.0 axisrot=0.0 rafield="RA" decfield="DEC" proj="STG" hpx_ordering_scheme="RING" hpx_order=3 hpx_ebin=yes hpx_region="" evtable="EVENTS" sctable="SC_DATA" efield="ENERGY" tfield="TIME" chatter=2 clobber=yes debug=no gui=no mode="ql" This is gtbin version ScienceTools-11-03-00 real 1.91 user 0.35 sys 0.22 time -p gtltcube evfile="filtered_events.fits" evtable="EVENTS" scfile=orbSim_scData_0000.fits sctable="SC_DATA" outfile=ltcube.fits dcostheta=0.025 binsz=1.0 phibins=0 tmin=0.0 tmax=0.0 file_version="1" zmin=0.0 zmax=180.0 chatter=2 clobber=yes debug=no gui=no mode="ql" Working on file orbSim_scData_0000.fits .....................! real 19.96 user 17.12 sys 1.41 time -p gtexpcube2 infile=ltcube.fits cmap=countsMap.fits outfile=bexpmap.fits irfs="DC1A" evtype="INDEF" nxpix="INDEF" nypix="INDEF" binsz="INDEF" coordsys="GAL" xref="INDEF" yref="INDEF" axisrot=0.0 proj="CAR" ebinalg="LOG" emin="INDEF" emax="INDEF" enumbins="INDEF" ebinfile="NONE" bincalc="EDGE" ignorephi=no thmax=180.0 thmin=0.0 table="EXPOSURE" chatter=2 clobber=yes debug=no mode="ql" Computing binned exposure map....................! real 26.00 user 24.90 sys 0.30 time -p gtexpcube2 infile=ltcube.fits cmap=countsMap.fits outfile=bexpmap_allsky.fits irfs="DC1A" evtype="INDEF" nxpix=360 nypix=180 binsz=1.0 coordsys="GAL" xref=0.0 yref=0.0 axisrot=0.0 proj="CAR" ebinalg="LOG" emin="INDEF" emax="INDEF" enumbins="INDEF" ebinfile="NONE" bincalc="EDGE" ignorephi=no thmax=180.0 thmin=0.0 table="EXPOSURE" chatter=2 clobber=yes debug=no mode="ql" Computing binned exposure map....................! real 61.98 user 60.61 sys 0.56 time -p gtsrcmaps scfile=orbSim_scData_0000.fits sctable="SC_DATA" expcube=ltcube.fits cmap=countsMap.fits srcmdl=/nfs/farm/g/glast/u52/rhel7-test/ScienceTools-scons/sane/data/srcModel_egretdiffuse.xml bexpmap=bexpmap.fits wmap=none outfile=sourceMaps.fits irfs="DC1A" evtype="INDEF" convol=yes resample=yes rfactor=2 minbinsz=0.1 ptsrc=yes psfcorr=yes emapbnds=yes copyall=no chatter=2 clobber=yes debug=no gui=no mode="ql" Generating SourceMap for Crab....................! Generating SourceMap for EGRET DiffuseCaught St13runtime_error at the top level: time -p gtsrcmaps scfile=orbSim_scData_0000.fits sctable="SC_DATA" expcube=ltcube.fits cmap=countsMap.fits srcmdl=/nfs/farm/g/glast/u52/rhel7-test/ScienceTools-scons/sane/data/srcModel_egretdiffuse.xml bexpmap=bexpmap.fits wmap=none outfile=sourceMaps.fits irfs="DC1A" evtype="INDEF" convol=yes resample=yes rfactor=2 minbinsz=0.1 ptsrc=yes psfcorr=yes emapbnds=no copyall=no chatter=2 clobber=yes debug=no gui=no mode="ql" Generating SourceMap for Crab....................! Generating SourceMap for EGRET Diffuse....................! Generating SourceMap for Geminga....................! Generating SourceMap for PKS 0528+134....................! real 244.39 user 238.31 sys 3.86 time -p gtsrcmaps scfile=orbSim_scData_0000.fits sctable="SC_DATA" expcube=ltcube.fits cmap=countsMap.fits srcmdl=/nfs/farm/g/glast/u52/rhel7-test/ScienceTools-scons/sane/data/srcModel_egretdiffuse.xml bexpmap=bexpmap_allsky.fits wmap=none outfile=sourceMaps_allsky.fits irfs="DC1A" evtype="INDEF" convol=yes resample=yes rfactor=2 minbinsz=0.1 ptsrc=yes psfcorr=yes emapbnds=no copyall=no chatter=2 clobber=yes debug=no gui=no mode="ql" Generating SourceMap for Crab....................! Generating SourceMap for EGRET Diffuse....................! Generating SourceMap for Geminga....................! Generating SourceMap for PKS 0528+134....................! real 243.39 user 237.10 sys 4.07 time -p gtlike irfs="DC1A" expcube=ltcube.fits srcmdl=/nfs/farm/g/glast/u52/rhel7-test/ScienceTools-scons/sane/data/srcModel_egretdiffuse.xml sfile=binned_fit_model.xml check_fit=yes results=results.dat specfile=counts_spectra.fits statistic="BINNED" optimizer="MINUIT" ftol=0.001 toltype="ABS" tsmin=no save=yes refit=no evfile=none evtable="EVENTS" scfile=orbSim_scData_0000.fits sctable="SC_DATA" expmap=none plot=no cmap=sourceMaps.fits bexpmap=bexpmap.fits wmap=none psfcorr=yes phased_expmap=none chatter=2 clobber=yes debug=no gui=no mode="ql" ********** ** 1 **SET PRINT .000 ********** ********** ** 2 **SET NOWARN ********** PARAMETER DEFINITIONS: NO. NAME VALUE STEP SIZE LIMITS 1 'Integral ' 1.0000 1.0000 .10000E-04 1000.0 2 'Index ' -2.0000 1.0000 -5.0000 -1.0000 3 'Prefactor ' 11.000 1.0000 .10000E-02 1000.0 4 'Integral ' 1.0000 1.0000 .10000E-04 1000.0 5 'Index ' -2.0000 1.0000 -5.0000 -1.0000 6 'Integral ' 1.0000 1.0000 .10000E-04 1000.0 7 'Index ' -2.0000 1.0000 -5.0000 -1.0000 ********** ** 3 **SET ERR .5000 ********** ********** ** 4 **SET GRAD 1.000 ********** ********** ** 5 **MINIMIZE 1500. 2.000 ********** MIGRAD MINIMIZATION HAS CONVERGED. MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX. FCN= 19816.90 FROM MIGRAD STATUS=CONVERGED 141 CALLS 142 TOTAL EDM= .72E-05 STRATEGY= 1 ERR MATRIX NOT POS-DEF EXT PARAMETER APPROXIMATE STEP FIRST NO. NAME VALUE ERROR SIZE DERIVATIVE 1 Integral 12.140 1.5834 .29222E-01 -.24924 2 Index -2.0836 .54515E-01 .78442E-01 -.10655 3 Prefactor 11.110 .18083 .60504E-02 -.21283 4 Integral 11.323 .68077 .19226E-01 -.16073 5 Index -1.7040 .23207E-01 .45737E-01 .29923E-01 6 Integral 11.961 25.038 .43760E-01 -.19713E-02 7 Index -2.5709 .96024 .78442E-01 -.33778E-02 ERR DEF= .500 Final values: Integral = 12.1396 Index = -2.08361 Prefactor = 11.11 Integral = 11.3231 Index = -1.70403 Integral = 11.9609 Index = -2.57088 ********** ** 6 **HESSE ********** FCN= 19816.90 FROM HESSE STATUS=OK 52 CALLS 194 TOTAL EDM= .71E-05 STRATEGY= 1 ERROR MATRIX ACCURATE EXT PARAMETER INTERNAL INTERNAL NO. NAME VALUE ERROR STEP SIZE VALUE 1 Integral 12.140 2.2474 .12918E-03 -1.3500 2 Index -2.0836 .77597E-01 .27692E-03 .47596 3 Prefactor 11.110 .21165 .26737E-04 -1.3596 4 Integral 11.323 1.1914 .84942E-04 -1.3576 5 Index -1.7040 .40413E-01 .20165E-03 2.4367 6 Integral 11.961 4.7317 .22326E-03 -1.3516 7 Index -2.5709 .18859 .49452E-03 .21624 ERR DEF= .500 Minuit fit quality: 3 estimated distance: 7.06209e-06 Minuit parameter uncertainties: 1 2.24752 2 0.0776215 3 0.211654 4 1.19147 5 0.0404182 6 4.73315 7 0.188886 ....! Computing TS values for each source (4 total) Photon fluxes are computed for the energy range 100 to 200000 MeV Crab: Integral: 12.1396 +/- 2.24752 Index: -2.08361 +/- 0.0776215 LowerLimit: 20 UpperLimit: 200000 TS value: 442.947 Flux: 2.123e-06 +/- 2.00225e-07 photons/cm^2/s EGRET Diffuse: Prefactor: 11.11 +/- 0.211654 Index: -2.1 Scale: 100 Flux: 0.000456706 +/- 8.69547e-06 photons/cm^2/s Geminga: Integral: 11.3231 +/- 1.19147 Index: -1.70403 +/- 0.0404182 LowerLimit: 20 UpperLimit: 200000 TS value: 1624.76 Flux: 3.6356e-06 +/- 5.91005e-07 photons/cm^2/s PKS 0528+134: Integral: 11.9609 +/- 4.73315 Index: -2.57088 +/- 0.188886 LowerLimit: 20 UpperLimit: 200000 TS value: 95.2913 Flux: 9.55637e-07 +/- 1.54084e-07 photons/cm^2/s WARNING: Fit may be bad in range [213.847, 275.51] (MeV) WARNING: Fit may be bad in range [977.933, 1259.92] (MeV) Total number of observed counts: 3991 Total number of model events: 3990.93 -log(Likelihood): 19816.90052 Writing fitted model to binned_fit_model.xml Elapsed CPU time: 8.23 real 9.25 user 8.11 sys 0.31 time -p gtlike irfs="DC1A" expcube=ltcube.fits srcmdl=/nfs/farm/g/glast/u52/rhel7-test/ScienceTools-scons/sane/data/srcModel_egretdiffuse.xml sfile=binned_fit_model.xml check_fit=yes results=results.dat specfile=counts_spectra.fits statistic="BINNED" optimizer="MINUIT" ftol=0.001 toltype="ABS" tsmin=no save=yes refit=no evfile=none evtable="EVENTS" scfile=orbSim_scData_0000.fits sctable="SC_DATA" expmap=none plot=no cmap=sourceMaps_allsky.fits bexpmap=bexpmap_allsky.fits wmap=none psfcorr=yes phased_expmap=none chatter=2 clobber=yes debug=no gui=no mode="ql" ********** ** 1 **SET PRINT .000 ********** ********** ** 2 **SET NOWARN ********** PARAMETER DEFINITIONS: NO. NAME VALUE STEP SIZE LIMITS 1 'Integral ' 1.0000 1.0000 .10000E-04 1000.0 2 'Index ' -2.0000 1.0000 -5.0000 -1.0000 3 'Prefactor ' 11.000 1.0000 .10000E-02 1000.0 4 'Integral ' 1.0000 1.0000 .10000E-04 1000.0 5 'Index ' -2.0000 1.0000 -5.0000 -1.0000 6 'Integral ' 1.0000 1.0000 .10000E-04 1000.0 7 'Index ' -2.0000 1.0000 -5.0000 -1.0000 ********** ** 3 **SET ERR .5000 ********** ********** ** 4 **SET GRAD 1.000 ********** ********** ** 5 **MINIMIZE 1500. 2.000 ********** MIGRAD MINIMIZATION HAS CONVERGED. MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX. FCN= 19812.67 FROM MIGRAD STATUS=CONVERGED 142 CALLS 143 TOTAL EDM= .84E-04 STRATEGY= 1 ERR MATRIX NOT POS-DEF EXT PARAMETER APPROXIMATE STEP FIRST NO. NAME VALUE ERROR SIZE DERIVATIVE 1 Integral 12.235 1.5922 .29247E-01 -.33081 2 Index -2.0853 .54441E-01 .78442E-01 -.11353 3 Prefactor 10.991 .17898 .60191E-02 1.6623 4 Integral 11.369 .68301 .19237E-01 .71522 5 Index -1.7048 .23199E-01 .45669E-01 -.43278 6 Integral 11.805 24.699 .43760E-01 -.17697E-01 7 Index -2.5641 .95870 .78442E-01 .37065E-02 ERR DEF= .500 Final values: Integral = 12.2349 Index = -2.08531 Prefactor = 10.9909 Integral = 11.3693 Index = -1.7048 Integral = 11.8053 Index = -2.56408 ********** ** 6 **HESSE ********** FCN= 19812.67 FROM HESSE STATUS=OK 52 CALLS 195 TOTAL EDM= .69E-04 STRATEGY= 1 ERROR MATRIX ACCURATE EXT PARAMETER INTERNAL INTERNAL NO. NAME VALUE ERROR STEP SIZE VALUE 1 Integral 12.235 2.2592 .12929E-03 -1.3491 2 Index -2.0853 .77478E-01 .27620E-03 .47501 3 Prefactor 10.991 .20944 .26599E-04 -1.3607 4 Integral 11.369 1.1953 .84991E-04 -1.3571 5 Index -1.7048 .40398E-01 .20135E-03 2.4372 6 Integral 11.805 4.6449 .22157E-03 -1.3531 7 Index -2.5641 .18738 .49393E-03 .21972 ERR DEF= .500 Minuit fit quality: 3 estimated distance: 6.90909e-05 Minuit parameter uncertainties: 1 2.25938 2 0.0775027 3 0.20944 4 1.19537 5 0.040403 6 4.64629 7 0.187673 ....! Computing TS values for each source (4 total) Photon fluxes are computed for the energy range 100 to 200000 MeV Crab: Integral: 12.2349 +/- 2.25938 Index: -2.08531 +/- 0.0775027 LowerLimit: 20 UpperLimit: 200000 TS value: 446.097 Flux: 2.13383e-06 +/- 2.00566e-07 photons/cm^2/s EGRET Diffuse: Prefactor: 10.9909 +/- 0.20944 Index: -2.1 Scale: 100 Flux: 0.000451808 +/- 8.60452e-06 photons/cm^2/s Geminga: Integral: 11.3693 +/- 1.19537 Index: -1.7048 +/- 0.040403 LowerLimit: 20 UpperLimit: 200000 TS value: 1629.22 Flux: 3.646e-06 +/- 5.92359e-07 photons/cm^2/s PKS 0528+134: Integral: 11.8053 +/- 4.64629 Index: -2.56408 +/- 0.187673 LowerLimit: 20 UpperLimit: 200000 TS value: 95.787 Flux: 9.53561e-07 +/- 1.53515e-07 photons/cm^2/s WARNING: Fit may be bad in range [213.847, 275.51] (MeV) WARNING: Fit may be bad in range [977.933, 1259.92] (MeV) Total number of observed counts: 3991 Total number of model events: 3991.21 -log(Likelihood): 19812.66808 Writing fitted model to binned_fit_model.xml Elapsed CPU time: 8.51 real 9.55 user 8.40 sys 0.30 time -p gtmodel srcmaps=sourceMaps_allsky.fits srcmdl=binned_fit_model.xml outfile=model_map.fits irfs="DC1A" evtype="INDEF" expcube=ltcube.fits bexpmap=bexpmap_allsky.fits convol=yes resample=yes rfactor=2 outtype="CMAP" psfcorr=yes edisp=no chatter=2 clobber=yes debug=no gui=no mode="ql" real 17.07 user 8.67 sys 5.95 time -p gtmodel srcmaps=sourceMaps_allsky.fits srcmdl=binned_fit_model.xml outfile=model_cube.fits irfs="DC1A" evtype="INDEF" expcube=ltcube.fits bexpmap=bexpmap_allsky.fits convol=yes resample=yes rfactor=2 outtype="ccube" psfcorr=yes edisp=no chatter=2 clobber=yes debug=no gui=no mode="ql" real 14.80 user 8.83 sys 4.02 time -p gtobssim infile=xmlFiles.txt srclist=source_names.txt scfile=none sctable="SC_DATA" evroot="Crab" evtable="EVENTS" simtime=86400.0 ltfrac=0.9 tstart="INDEF" nevents=no maxtime=315500000.0 startdate="2001-01-01 00:00:00" offset=0 rockangle="INDEF" use_ac=no ra=0.0 dec=0.0 radius=20.0 emin=1.0 emax=1000000.0 edisp=yes irfs="DC1A" evtype="none" area=1.0 maxrows=1000000 seed=62587 chatter=2 clobber=yes debug=no gui=no mode="ql" added source "_3EG_J0534p2200-32mev" Generating events for a simulation time of 86400 seconds.... Done. real 1.53 user 0.66 sys 0.06 time -p gtselect infile=Crab_events_0000.fits outfile=Crab_events_filtered.fits ra=83.57 dec=22.01 rad=20.0 tmin=0.0 tmax=0.0 emin=30.0 emax=300000.0 zmin=0.0 zmax=180.0 evclass="INDEF" evtype="INDEF" convtype=-1 phasemin=0.0 phasemax=1.0 evtable="EVENTS" chatter=2 clobber=yes debug=no gui=no mode="ql" Done. real 0.91 user 0.15 sys 0.04 time -p gtbin evfile=Crab_events_filtered.fits scfile=Crab_scData_0000.fits outfile=Crab.fits algorithm="CMAP" ebinalg="LOG" emin=30.0 emax=200000.0 enumbins=20 denergy=0.0 ebinfile=NONE tbinalg="LIN" tstart=0.0 tstop=86400.0 dtime=1000.0 tbinfile=numxpix=100 snratio=1.0 lcemin=0.0 lcemax=0.0 nxpix=100 nypix=100 binsz=0.5 coordsys="CEL" xref=83.0 yref=22.0 axisrot=0.0 rafield="RA" decfield="DEC" proj="CAR" hpx_ordering_scheme="RING" hpx_order=3 hpx_ebin=yes hpx_region="" evtable="EVENTS" sctable="SC_DATA" efield="ENERGY" tfield="TIME" chatter=2 clobber=yes debug=no gui=no mode="ql" This is gtbin version ScienceTools-11-03-00 real 1.19 user 0.20 sys 0.04 time -p gtbin evfile=Crab_events_filtered.fits scfile=Crab_scData_0000.fits outfile=Crab.pha algorithm="PHA1" ebinalg="LOG" emin=30.0 emax=200000.0 enumbins=20 denergy=0.0 ebinfile=NONE tbinalg="LIN" tstart=0.0 tstop=86400.0 dtime=1000.0 tbinfile=numxpix=100 snratio=1.0 lcemin=0.0 lcemax=0.0 nxpix=100 nypix=100 binsz=0.5 coordsys="CEL" xref=83.0 yref=22.0 axisrot=0.0 rafield="RA" decfield="DEC" proj="CAR" hpx_ordering_scheme="RING" hpx_order=3 hpx_ebin=yes hpx_region="" evtable="EVENTS" sctable="SC_DATA" efield="ENERGY" tfield="TIME" chatter=2 clobber=yes debug=no gui=no mode="ql" This is gtbin version ScienceTools-11-03-00 real 1.01 user 0.19 sys 0.03 time -p gtbin evfile=Crab_events_filtered.fits scfile=Crab_scData_0000.fits outfile=Crab.lc algorithm="LC" ebinalg="LOG" emin=30.0 emax=200000.0 enumbins=20 denergy=0.0 ebinfile=NONE tbinalg="LIN" tstart=0.0 tstop=86400.0 dtime=1000.0 tbinfile=numxpix=100 snratio=1.0 lcemin=0.0 lcemax=0.0 nxpix=100 nypix=100 binsz=0.5 coordsys="CEL" xref=83.0 yref=22.0 axisrot=0.0 rafield="RA" decfield="DEC" proj="CAR" hpx_ordering_scheme="RING" hpx_order=3 hpx_ebin=yes hpx_region="" evtable="EVENTS" sctable="SC_DATA" efield="ENERGY" tfield="TIME" chatter=2 clobber=yes debug=no gui=no mode="ql" This is gtbin version ScienceTools-11-03-00 real 0.98 user 0.18 sys 0.04 time -p gtrspgen respalg="PS" specfile=Crab.pha scfile=Crab_scData_0000.fits outfile=Crab.rsp irfs="DC1AF" sctable="SC_DATA" resptpl=DEFAULT chatter=2 clobber=yes debug=no gui=no mode="ql" time=1000.0 thetacut=70.0 dcostheta=0.025 phinumbins=1 ebinalg="LOG" efield="ENERGY" emin=30.0 emax=200000.0 enumbins=20 denergy=0.0 ebinfile=NONE This is gtrspgen version ScienceTools-11-03-00 real 1.29 user 0.27 sys 0.05