Niu: Created page with "=== Insets as zoom-ins === center

 require("isisscripts");  %%% all necessary data  % get ObsIDs, tstart, and tstop of all used data from so..."

2018-05-02T11:38:39Z

Created page with "=== Insets as zoom-ins === center <pre> require("isisscripts"); %%% all necessary data % get ObsIDs, tstart, and tstop of all used data from so..."

New page

=== Insets as zoom-ins ===

[[File:gx304batlcobs.png|center]]

<pre>
require("isisscripts");

%%% all necessary data

% get ObsIDs, tstart, and tstop of all used data from somewhere
variable data = struct {
obsid = [...],
tstart = [...],
tstop = [...]
}

% load BAT-lightcurve
variable bat = fits_read_table("source.lc.fits");

%%% initialize xfig plot

variable W = 28, H = 10; % width and height of the full-lc plot
variable wt0 = 54850, wt1 = 56800; % xrange in MJD of this plot
variable Wi = 5.2, Hi = 3.5; % width and height of the zoom plots
variable pl = xfig_plot_new(W,H);
% set world; the negative, small pad-values removes the tic-labels
% at the end of each axis (my personal liking)
pl.world(wt0, wt1, 0, 0.8; padx = -1e-6, pady = -1e-6);
% set 2nd world, here the years
pl.world2(yearOfMJD(pl.get_world()[0]), yearOfMJD(pl.get_world()[1]),
__push_array(pl.get_world()[[2,3]]); padx = -1e-6, pady = -1e-6);
pl.x2axis(; major = [2009:2014]+.5-1e-6, minor = [2009:2015],
major_len = 0, minor_len = .2, format= " %.0f");
pl.y2axis(; ticlabels = 0);
pl.y1axis(; format = "%.2f");
pl.xlabel("Time (MJD)");
pl.ylabel("BAT-cts\,s$^{-1}$\,cm$^{-2}$"R);

%%% plot lightcurve

pl.plot(bat.time, bat.rate, bat.error; color = "red", sym = "x",
width = 1, size = .5);

%%% create zoom plots

variable pli = Struct_Type[4]; % an array of structures; for each plot one element
variable t0 = [55260, 55395, 55520, 55655]; % the start (in MJD) for each plot
variable dt = 60; % the length (in days) for every plot (such that the scaling is equal for all)
variable ymx = [0.05, 0.42, 0.27, 0.28]; % the maximum y-value for each plot (in cts/s/cm^2)
variable dx = .24; % the relative spacing between each plot (measured from their centers)

% loop over all plots and create them
variable i;
_for i (0, length(pli)-1, 1) {
%%% initialize the zoom plots
pli[i] = xfig_plot_new(Wi, Hi);
pli[i].world(t0[i], t0[i]+dt, 0, ymx[i]*1.1; pady = -1e-6);
pli[i].xaxis(; major = t0[i]+[.25,.75]*dt, minor = t0[i]+[0:dt:10],
ticlabel_size = "footnotesize");
pli[i].y1axis(; ticlabel_size = "footnotesize", format = "%.2f");
pli[i].x2axis(; ticlabels = 0);

%%% plot the lightcurve
pli[i].plot(bat.time, bat.rate, bat.error; color = "red", sym = "x",
width = 1, size = .5);
% add the observation times
pli[i].plot(.5*(data.tstart+data.tstop), ones(length(data.tstart))*ymx[i];
sym = "darr", color = "blue");

%%% add the object to the full lightcurve plot
variable ix = .5*(1.-dx*length(pli)) + dx*(i+.5); % relative x-coordinate of the zoom plot
pl.add_object(pli[i], ix, .55, 0, -.5; world0);

%%% calculate the xy-positions for the lines connecting
%%% the zoom plots with the full lightcurve

% relative x/y-coordinates of the lower-[left,right]-corner of the zoom plot
variable xi = (pli[i].plot_data.X.x + [0, Wi]) / W;
variable yi = pli[i].plot_data.X[0].y / H * [1,1];
% relative x/y-coordinates of the [beginning,end] of the zoom time range _in_ the full lightcurve plot
variable xl = 1.*(t0[i]+[0,dt]-wt0) / (wt1-wt0);
variable yl = (max(bat.rate[where(t0[i] < bat.time < t0[i]+dt)]) + [0.02, 0.04]) / .8;

%%% plot the connecting lines
pl.plot([xl[0], xl[0], xi[0]], [yl, yi[0]]; world0, color = "gray");
pl.plot([xl[1], xl[1], xi[1]], [yl, yi[1]]; world0, color = "gray");
% shade the area in between
pl.shade_region(
[xl[0], xl[0], xi[0], xi[1], xl[1], xl[1], xl[1]],
[yl[0], yl[1], yi[0], yi[1], yl[1], yl[0], yl[0]];
color = "#DDDDDD", world0
);
}

% render
pl.render("BAT_lc.pdf");
</pre>

[[Category:SLxfig]]

GX304batlcobs (xfig example) - Revision history

Niu: Created page with "=== Insets as zoom-ins === center
require("isisscripts"); %%% all necessary data % get ObsIDs, tstart, and tstop of all used data from so..."