Chandra X-Ray Observatory
	(CXC)
Skip to the navigation links
Last modified: December 2013

URL: http://cxc.harvard.edu/sherpa/ahelp/jdpileup.html
Jump to: Description · Example · Bugs · See Also


AHELP for CIAO 4.6 Sherpa v2

jdpileup

Context: models

Synopsis

A CCD pileup model developed by John Davis of MIT.

Description

The 1D pileup model can be applied to model the 1D Chandra spectrum obtained in the imaging mode. The pileup model does not work for pileup in dispersed grating spectra or 2D image data. The model was designed for imaging pileup, including pileup by the gratings in zeroth order. It should be used only for energy spectra.

The use of multiple optimization methods is desirable when fitting data with pileup model.

jdpileup Parameters

Number Name Description
1 alpha probability of a good grade when two photons pile together
2 g0 probability of grade 0 assignment
3 f fraction of flux falling into the pileup region
4 n number of detection cells
5 ftime frame time [seconds] (keyword EXPTIME in the event file)
6 fracexp fractional exposure that the point source experienced while dithering on the chip (keyword FRACPROB in the ARF file)
7 nterms maximum number of photons considered for pileup in a single frame.

n, f and g0: The values of n, f and g0 should remain frozen; The full discussion of these parameters is presented in Davis (2001).

alpha: The value of the parameter alpha should be allowed to vary. alpha parameterizes "grade migration" in the detector, and represents the probability, per photon count greater than one, that the piled event is not rejected by the spacecraft software as a "bad event". Specifically, if n photons are piled together in a single frame, the probability of them being retained (as a single photon event with their summed energy) is given by alpha(n-1). In reality, the alpha parameter should be a photon-energy-dependent and detector-chip-dependent matrix; for simplicity, the jdpileup model assumes a constant value.

ftime: The ftime and parameter should be set to the value given in the header keyword EXPTIME of the event file. (Note that EXPTIME is used instead of TIMEDEL because the latter includes the transfer time, which ftime should not.)

fracexp: The fracexp parameter should be set to the value given in the header keyword FRACEXPO of the ARF file.

nterms: This should be left frozen at its maximum value of 30. I.e., the expansion of the model will include terms corresponding to 0, 1, 2, .. 30 photon events landing in the same extraction region during the same frame time.

The integration of models in Sherpa is controlled by an integration flag in each model structure. Refer to "ahelp integrate" for information on integrating model components.

Example

sherpa> set_source(xsphabs.gal*powlaw1d.p1)
sherpa> set_pileup_model(jdpileup.jd)
sherpa> show_model()
Model: 1
apply_rmf(jdpileup.jd(xsphabs.gal * powlaw1d.p1))
   Param        Type          Value          Min          Max      Units
   -----        ----          -----          ---          ---      -----
   jd.alpha     thawed          0.5            0            1
   jd.g0        frozen            1  1.17549e-38            1
   jd.f         thawed         0.95          0.9            1
   jd.n         frozen            1  1.17549e-38          100
   jd.ftime     frozen        3.241  1.17549e-38            5        sec
   jd.fracexp   frozen        0.987            0            1
   jd.nterms    frozen           30            1          100
   gal.nH       frozen        0.014            0       100000 10^22
atoms / cm^2
   p1.gamma     thawed      1.37542          -10           10
   p1.ref       frozen            1 -3.40282e+38  3.40282e+38
   p1.ampl      thawed  0.000238509            0  3.40282e+38

Define an absorbed power-law source model for fitting to data set 1, and add a pileup model component. Print the pileup model parameter values along with the source model parameters.

Bugs

See the bugs pages on the Sherpa website for an up-to-date listing of known bugs.

See Also

models
absorptionedge, absorptiongaussian, absorptionlorentz, absorptionvoigt, accretiondisk, atten, bbody, bbodyfreq, beta1d, beta2d, blackbody, box1d, box2d, bpl1d, bremsstrahlung, brokenpowerlaw, ccm, const1d, const2d, cos, delta1d, delta2d, dered, devaucouleurs2d, disk2d, edge, emissiongaussian, emissionlorentz, emissionvoigt, erf, erfc, exp, exp10, fm, gauss1d, gauss2d, hubblereynolds, linebroad, list_model_components, list_models, lmc, log, log10, logabsorption, logemission, logparabola, lorentz1d, lorentz2d, models, normbeta1d, normgauss1d, normgauss2d, opticalgaussian, poisson, polynom1d, polynom2d, polynomial, powerlaw, powlaw1d, recombination, scale1d, scale2d, schechter, seaton, sersic2d, shell2d, sin, sm, smc, sqrt, stephi1d, steplo1d, tablemodel, tan, xgal, xs, xsabsori, xsacisabs, xsapec, xsbapec, xsbbody, xsbbodyrad, xsbexrav, xsbexriv, xsbkn2pow, xsbknpower, xsbmc, xsbremss, xsbvapec, xsbvvapec, xsc6mekl, xsc6pmekl, xsc6pvmkl, xsc6vmekl, xscabs, xscemekl, xscevmkl, xscflow, xscompbb, xscompls, xscompps, xscompst, xscompth, xscomptt, xsconstant, xscplinear, xscutoffpl, xscyclabs, xsdisk, xsdiskbb, xsdiskir, xsdiskline, xsdiskm, xsdisko, xsdiskpbb, xsdiskpn, xsdust, xsedge, xseplogpar, xseqpair, xseqtherm, xsequil, xsexpabs, xsexpdec, xsexpfac, xsezdiskbb, xsgabs, xsgadem, xsgaussian, xsgnei, xsgrad, xsgrbm, xshighecut, xshrefl, xskerrbb, xskerrd, xskerrdisk, xslaor, xslaor2, xslogpar, xslorentz, xsmeka, xsmekal, xsmkcflow, xsnei, xsnotch, xsnpshock, xsnsa, xsnsagrav, xsnsatmos, xsnsmax, xsnteea, xsnthcomp, xsoptxagn, xsoptxagnf, xspcfabs, xspegpwrlw, xspexmon, xspexrav, xspexriv, xsphabs, xsplabs, xsplcabs, xsposm, xspowerlaw, xspshock, xspwab, xsraymond, xsredden, xsredge, xsrefsch, xssedov, xssirf, xssmedge, xsspexpcut, xsspline, xssrcut, xssresc, xssss_ice, xsstep, xsswind1, xstbabs, xstbgrain, xstbvarabs, xsuvred, xsvapec, xsvarabs, xsvbremss, xsvequil, xsvgadem, xsvgnei, xsvmcflow, xsvmeka, xsvmekal, xsvnei, xsvnpshock, xsvphabs, xsvpshock, xsvraymond, xsvsedov, xsvvapec, xswabs, xswndabs, xsxion, xszbbody, xszbremss, xszdust, xszedge, xszgauss, xszhighect, xszigm, xszpcfabs, xszphabs, xszpowerlw, xszredden, xszsmdust, xsztbabs, xszvarabs, xszvfeabs, xszvphabs, xszwabs, xszwndabs, xszxipcf

Last modified: December 2013
Smithsonian Institute Smithsonian Institute

The Chandra X-Ray Center (CXC) is operated for NASA by the Smithsonian Astrophysical Observatory. 60 Garden Street, Cambridge, MA 02138 USA.   Email:   cxcweb@head.cfa.harvard.edu Smithsonian Institution, Copyright © 1998-2014. All rights reserved.