NAME: lhr_us PURPOSE: This routine calculates and returns (2d array) the likelyhood ratio for a set of identifications using a method similar to Wolstencroft et Al. 1986 with the product Q1*Q2 replaced by Q=1 It also calculates individual reliabilities according to the Sutherland & Saunders (1992) method and returns this and other vectors on demand It then goes onto produce various histograms (except if noplot set) CALLING SEQUENCE: result=lhr_us(xmmdb,wfcdb[,dmax,radec_err=radec_err,$ radec_fromdb=radec_fromdb,rel=rel,]top=top,bot=bot,dr=dr,field=field,$ magmax=magmax,noplot=noplot,one_per_page=one_per_page) INPUTS: xmmdb: xmm database wfcdb: wfc dbatabase for extracting magnitudes of candidates OPTIONAL INPUT PARAMETERS: dmax: max. offset for nearest matches radec_err: no (an error to be assigned to every xmm position) radec_fromdb: if set, position error extracted for each object from xmmdb magmax: if set to a magnitude value, then only objects with mag > max are included. This can only be used if wfcdb has an item labelled 'mag' OPTIONAL OUTPUT PARAMETERS: rel: if set to a named variable, gets filled with 2d array of individual reliability (lhr/(1+lhr)) values rss: if set to a named variable, gets filled with 2d array of individual reliability (S&S) values rrut: if set to a named variable, gets filled with 2d array of individual reliability (Rutledge et al) values last element has p(no-id) top: if set to a named variable, gets filled with 2d array of numerator of lhr values bot: if set to a named variable, gets filled with 2d array of denominator of lhr values dr: if set to a named variable, gets filled with 2d array of offsets noplot: if set, no histogram plots are produced one_per_page: if set as well as ps, only one plot per page RESTRICTIONS: UNTESTED CALLED PROCEDURES: MODIFICATION HISTORY: Original by Ignas Snellen 08-Dec-2000 S. Aigrain: adapted as a separate function 16-Jan-2001 S. Aigrain: this method separately 18-Feb-2001 S. Aigrain: Rutledge et al. reliability included, one more page of plots.
(See /home/xmmssc/soft/idl/pro/lhr/lhr_us.pro)