try: if just_prefactor: fit_prefactor(roi, name) else: roi.fit(fit_bg_first=fit_bg_first) except Exception, ex: print 'ERROR spectral fitting pointlike for hypothesis %s:' % hypothesis, ex traceback.print_exc(file=sys.stdout) print_summary() fit(just_prefactor=True) fit(fit_bg_first=True) fit() if localize: paranoid_localize(roi, name, verbosity=4) if fit_extension: roi.fit_extension(which=name) paranoid_localize(roi, name) fit() print 'Making pointlike SED for hypothesis %s' % hypothesis sed = PointlikeSED(roi, name, verbosity=4) sed.save('%s/sed_pointlike_4bpd_%s_%s.yaml' % (pipeline.dirdict['seds'], hypothesis, name)) sed.plot('%s/sed_pointlike_4bpd_%s_%s.png' % (pipeline.dirdict['seds'], hypothesis, name)) print_summary()
try: if just_prefactor: fit_prefactor(roi, name) else: roi.fit(fit_bg_first=fit_bg_first) except Exception, ex: print 'ERROR spectral fitting pointlike for hypothesis %s:' % hypothesis, ex traceback.print_exc(file=sys.stdout) print_summary() fit(just_prefactor=True) fit(fit_bg_first=True) fit() if localize: paranoid_localize(roi, name, verbosity=4) if fit_extension: roi.fit_extension(which=name) paranoid_localize(roi, name) fit() print 'Making pointlike SED for hypothesis %s' % hypothesis sed = PointlikeSED(roi, name, verbosity=4) sed.save('%s/sed_pointlike_4bpd_%s_%s.yaml' % (pipeline.dirdict['seds'],hypothesis,name)) sed.plot('%s/sed_pointlike_4bpd_%s_%s.png' % (pipeline.dirdict['seds'],hypothesis,name)) print_summary() p = source_dict(roi, name)
traceback.print_exc(file=sys.stdout) print_summary() fit(just_prefactor=True) fit(fit_bg_first=True) fit() frozen = freeze_far_away(roi, roi.get_source(name).skydir, max_free) if localize: if override_localization is None: print 'About to Grid localize for hypothesis %s' % hypothesis grid=GridLocalize(roi, name, size=0.5, pixelsize=0.1, verbosity=4) print_summary() ellipse = paranoid_localize(roi, name, verbosity=4) print 'Localization Ellipse:',ellipse else: print 'Override Localization for source %s:', name print 'override_localization=',override_localization roi.modify(which=None, skydir=override_localization['init_position']) assert override_localization['method'] == 'MinuitLocalizer' m=MinuitLocalizer(roi, name, verbosity=4) print 'Localization Ellipse:',m.todict() if fit_extension: roi.fit_extension(which=name) ellipse = paranoid_localize(roi, name) print ellipse