Loving ISIS - Confessions of a Former XSPEC User

Summary:

OK, that was a lot of territory to cover in a few pages. But the goal was to cover 98% of the ISIS commands that I need more than once in a blue moon. Plus I wanted to show a worked example of ISIS applied to something that's not Chandra gratings data. ISIS can be used on lots of different kinds of data sets.

Just to demonstrate that it's not nearly as difficult as it might have seemed, here is a condensed version of the analysis script. This version skips some of the slower stuff (e.g., reflection models), and just does a doubly broken power law plus gaussian model fit to the radio/X-ray data. It does use my .isisrc scripts. Including calculation of a complete set of error bars, and calculation of the Gamma_2 vs. Gamma_3 contours (961 grid points), the script runs in about 23 minutes on my 1.5 GHz Powerbook G4. The data for this script can be found here.

   variable radio_data = read_radio("radio.data",50);
   variable radio_id = load_radio(radio_data);

   variable pca_id = load_all("pca.pha");
   set_systematics(pca_id,[0],[50],[0.005]);
   grppha_cut(pca_id,30,3.);
   kev_note(pca_id,3.,22.);

   variable hexte_id = load_hexte("hxt.pha");
   grppha_sn(hexte_id,20.,18.);
   kev_note(hexte_id,18.,200.);

   m_width = 3, r_width=3, d_width=3;
   set_frame_line_width(3); charsize(1.12); point_size(1.5);
   xlog; ylog;
   plot_bsub_counts(2,"xrng",[3.,22],"yrng",[1.01e4,4.3e5]);

   model("constant(Isis_Active_Dataset)*phabs(1)*highecut(1)*(bkn2pow(1)+gaussian(1))");

   newpar(2,0.6,-1);
   newpar(3,40);
   newpar(4,130);
   newpar(5,230);
   newpar(6,0.83);
   newpar(7,4.e-4,0,1.e-6,1.e-2);
   newpar(8,1.68);
   newpar(9,11,0,9,12);
   newpar(10,1.48);
   newpar(11,2.e-3);
   newpar(12,6.4,-1);
   newpar(13,0.3,0,0,1.5);
   newpar(15,0.95,0,0.9,1.1);
   freeze([1,14]);

   () = renorm_counts;
   () = fit;

   variable a = [3,4,5,6,7,8,9,10,11,13,15];
   variable pmin,pmax;
   (pmin,pmax) = error_loop(a,,"bkn2pow_gaussian");
   error_save(a,pmin,pmax,"bkn2pow_gaussian.errors");

   variable px = conf_grid(8,1.67,1.70,31);
   variable py = conf_grid(10,1.43,1.52,31);
   variable cntr = conf_map(px,py);

   save_conf(cntr,"bkn2pow_gauss_g2_vs_g3.contour");

   xlin; ylin; xrange; yrange;

   variable id = open_print("bkn2pow_gauss_g2_vs_g3_contour.ps/vcps");
   resize(21,0.85);

   xlabel("\\fr Photon Index, \\gG\\d2");
   ylabel("\\fr Photon Index, \\gG\\d3");
   plot_conf(cntr);

   close_print(id);

Want more? Periodically check out the Power Tools page.