addcomp

To find and download HEASARC data in the cloud, you can use astroquery.heasarc or download our new tool, hark.

The HEASARC and NuSTAR teams are greatly saddened by the sudden passing of Katja Pottschmidt. Most recently Katja was the lead scientist for the NuSTAR Guest Observer Facility (GOF), a role she had supported for many years. During her science career she worked on many other high energy astrophysics missions and played an integral role in advancing our knowledge of the universe. She was a wonderful colleague and friend and will be keenly missed by all who knew her.

addcomp

add component to a model

Add a component to the model.

Syntax:	addcomp	[<modeName>:]<n> <comp>

where <n> is the component number before which the new component is to be inserted, and <comp> is the name of the new component. Components are numbered in sequence in order of appearance in the expression entered. The new component is regarded as an operator on the component added if it is not additive.

The optional <modelName> qualifier allows the user to address a named model.

The user is prompted for parameter values for the component. If there are m components in the current model, then acceptable values for the component number added are 1 to m+1.

XSPEC detects the type of the model (additive, multiplicative etc), checks the correctness of the syntax of the output model, and adds the component if the resulting models obeys the syntax rules documented in the model command.

Thus,

XSPEC12> mo  wa(po)

Followed by

XSPEC12> addcomp 2 bb

Yields the model achieved by

XSPEC12> mo  wa(bb + po)

See also delcomp (delete component by number).

Other examples will serve to clarify addcomp's behavior. Suppose that the current model specification is

ga+po

which using the show command would yield the description

model  = gaussian[1] + powerlaw[2]

The comments give the model expression following the entry of addcomp and delcomp commands:

XSPEC12> addcomp 2 wab 
// gaussian[1] + wabs[2](powerlaw[3])
XSPEC12> addcomp 4 pha 
// (gaussian[1] + wabs[2](powerlaw[3]))phabs[4]
XSPEC12> delcomp 1	
// (wabs[1](powerlaw[2]))phabs[3]}
XSPEC12> addcomp 2 zg 
// (wabs[1](zgauss[2] + powerlaw[3]))phabs[4]
XSPEC12> delcomp 3 
// (wabs[1](zgauss[2]))phabs[3]

XSPEC12> mo wa(po)
XSPEC12> addcomp 1 ga
// gauss[1] + wabs[2]*powerlaw[3]
XSPEC12> delcomp 1
XSPEC12> addcomp 1 pha
// phabs[1]*wabs[2]*powerlaw[3]

XSPEC12>mo wabs(po)
XSPEC12> addcomp 3 bb
// wabs[1]*powerlaw[2] + bbody[3]
XSPEC12> delcomp 1
XSPEC12> addcomp 3 pha
// wabs[1]*powerlaw[2]*pha[3]
XSPEC12> addcomp 3 po
// ERROR: po (additive) is interpreted as being added to the 
// multiplicative model pha[3], which is a context error.

For multiply nested models?
XSPEC12> mo wa(po + pha(bb + ga))
XSPEC12> addcomp 6 po
// wabs[1](powerlaw[2] + phabs[3](bbody[4] + ga[5]) + powerlaw[6])
XSPEC12> addcomp 5 peg
// wabs[1](powerlaw[2] + phabs[3](bbody[4] + pegpwlw[5] ga[6]) + powerlaw[7])
XSPEC12> addcomp 7 wa
// wabs[1](powerlaw[2] + phabs[3](bbody[4] + pegpwlw[5] ga[6]) + wabs[7]*powerlaw[8])

NASA | GSFC | Sciences and Exploration

HEASARC

High Energy Astrophysics Science Archive Research Center

addcomp