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)


Suzanne Aigrain <visitor2@ast.cam.ac.uk>
Last modified: Mon Feb 19 11:59:23 2001