*-- Author : Stephen Burke
SUBROUTINE FKHBK
*-----------------------------------------Updates 24/01/92-------
**: FKHBK 30205.SB. Phi histograms added.
*-----------------------------------------Updates----------------
**********************************************************************
* *
* Book diagnostic histograms *
* *
**********************************************************************
*KEEP,FKNPL.
*
* Per-track values can go in H1WORK; note that LTRUE and LFIRST must
* be set at least per event.
*
* This is about 36k words long; the remaining common blocks are
* about 3.6k in total. Some of this could be in /H1WORK/, but the
* blocks would have to be reorganised.
*
* /FKPROJ/
* /FKFILT/
* /FKSMTH/
* /FKINT/
* /FKRSID/
* /FKTRUE/
* /FKDBG/
*KEEP,FFSTEE.
*KEEP,FKFLAG.
*KEEP,FKTRUE.
*KEEP,FKHIST.
*KEEP,FKDBG.
*KEND.
**********************************************************************
CALL HCDIR('//PAWC',' ')
CALL HMDIR(CKDBG,'S')
* Decode steering
CALL HBOOK1(101,'PULL ON X$',60,-6.,6.,0.)
CALL HBOOK1(102,'PULL ON Y$',60,-6.,6.,0.)
CALL HBOOK1(103,'PULL ON Q/P$',60,-6.,6.,0.)
CALL HBOOK1(104,'PULL ON TAN(THETA)$',60,-6.,6.,0.)
CALL HBOOK1(105,'PULL ON PHI$',60,-6.,6.,0.)
CALL HBOOK2(230,'DELTA(1/P) VS NSEG (LOW p)$',10,0.,10.,
CALL HBOOK2(231,'DELTA(1/P) VS NSEG (MED P)$',10,0.,10.,
CALL HBOOK2(232,'DELTA(1/P) VS NSEG (HIGH P)$',10,0.,10.,
CALL HBOOK2(233,'DELTA(1/P) VS LOG(P)$',60,-3.,3.,
CALL HBOOK2(234,'DELTA(THETA) VS LOG(P)$',60,-3.,3.,
CALL HBOOK2(235,'DELTA(PHI) VS LOG(P)$',60,-3.,3.,
CALL HBOOK2(243,'DELTA(1/P) VS THETA$',80,0.,0.4,
CALL HBOOK2(244,'DELTA(THETA) VS THETA$',80,0.,0.4,
CALL HBOOK2(245,'DELTA(PHI) VS THETA$',80,0.,0.4,
CALL HBOOK2(253,'DELTA(1/P) VS PHI$',64,0.,6.4,
CALL HBOOK2(254,'DELTA(THETA) VS PHI$',64,0.,6.4,
CALL HBOOK2(255,'DELTA(PHI) VS PHI$',64,0.,6.4,
CALL HBOOK1(301,'DRIFT RESID FOR OLD TRUE POINTS (PLANAR)$'
CALL HBOOK1(305,'CHISQ FOR OLD TRUE POINTS (PLANAR)$'
CALL HBOOK1(306,'DRIFT RESID FOR OLD TRUE POINTS (RADIAL)$'
CALL HBOOK1(307,'RADIUS RESID FOR OLD TRUE POINTS (RADIAL)$'
CALL HBOOK1(310,'CHISQ FOR OLD TRUE POINTS (RADIAL)$'
CALL HBOOK1(311,'DRIFT RESID FOR OLD WRONG POINTS (PLANAR)$'
CALL HBOOK1(315,'CHISQ FOR OLD WRONG POINTS (PLANAR)$'
CALL HBOOK1(316,'DRIFT RESID FOR OLD WRONG POINTS (RADIAL)$'
CALL HBOOK1(317,'RADIUS RESID FOR OLD WRONG POINTS (RADIAL)$'
CALL HBOOK1(320,'CHISQ FOR OLD WRONG POINTS (RADIAL)$'
CALL HBOOK1(321,'DRIFT RESID FOR NEW TRUE POINTS (PLANAR)$'
CALL HBOOK1(325,'CHISQ FOR NEW TRUE POINTS (PLANAR)$'
CALL HBOOK1(326,'DRIFT RESID FOR NEW TRUE POINTS (RADIAL)$'
CALL HBOOK1(327,'RADIUS RESID FOR NEW TRUE POINTS (RADIAL)$'
CALL HBOOK1(330,'CHISQ FOR NEW TRUE POINTS (RADIAL)$'
CALL HBOOK1(331,'DRIFT RESID FOR NEW WRONG POINTS (PLANAR)$'
CALL HBOOK1(335,'CHISQ FOR NEW WRONG POINTS (PLANAR)$'
CALL HBOOK1(336,'DRIFT RESID FOR NEW WRONG POINTS (RADIAL)$'
CALL HBOOK1(337,'RADIUS RESID FOR NEW WRONG POINTS (RADIAL)$'
CALL HBOOK1(340,'CHISQ FOR NEW WRONG POINTS (RADIAL)$'
CALL HBOOK1(IPL+500,'FILTERED PROB$',100,0.,1.,0.)
CALL HBOOK1(IPL+600,'SMOOTHED PROB$',100,0.,1.,0.)
IF (LTRUTH) CALL HBOOK1(IPL+700,'SMOOTHED PROB (TRUE)$'
CALL HBOOK1(IPL+800,'DRIFT RESIDUALS$',100,-0.1,0.1,0.)
CALL HBOOK1(IPL+900,'RADIAL RESIDUALS$',100,-10.,10.,0.)
CALL HBPRO(0,0.)
CALL HMINIM(0,0.)
CALL HIDOPT(0,'INTE')
CALL HCDIR('//PAWC',' ')
*