*-- Author : Stephen Burke
SUBROUTINE FKANAL(JPL,RES,CHISQ,NPS,NRS)
*-----------------------------------------Updates 06/08/92-------
**: FKANAL.......SB. Minor bug fixed.
*-----------------------------------------Updates 24/01/92-------
**: FKANAL 30205.SB. Digi with wrong drift sign counts as a bad hit.
**: FKANAL 30205.SB. Phi histograms added.
*-----------------------------------------Updates----------------
**********************************************************************
* *
* Analyse the KF output *
* *
**********************************************************************
DOUBLE PRECISION RES(2),CHISQ,FKDPHI
*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,FKCNTL.
*KEEP,FKMEAS.
*KEEP,FKTRUE.
*KEEP,FKPROJ.
*KEEP,FKSMTH.
*KEEP,FKRSID.
*KEEP,FKFLAG.
*KEEP,FFDBG.
*KEEP,FKHIST.
*KEEP,FKDBG.
*KEND.
**********************************************************************
CALL HCDIR('//PAWC/'//CKDBG,' ')
CALL HFILL(IBASE+1+5*IRP(JPL),SNGL(RES(1)),0.,1.)
IF (MES(JPL).EQ.2) CALL HFILL(IBASE+12,SNGL(RES(2)),0.,1.)
CALL HFILL(IBASE+5+5*IRP(JPL),SNGL(CHISQ),0.,1.)
* Decode steering
& CALL HFILL(101,SNGL((SSMT(1,JSTART)-TRUE(1,JSTART))/
& CALL HFILL(102,SNGL((SSMT(2,JSTART)-TRUE(2,JSTART))/
& CALL HFILL(103,SNGL((SSMT(3,JSTART)-TRUE(3,JSTART))/
& CALL HFILL(104,SNGL((SSMT(4,JSTART)-TRUE(4,JSTART))/
& CALL HFILL(105,SNGL((SSMT(5,JSTART)-TRUE(5,JSTART))/
C D1OP = DABS(SSMT(3,JSTART)) - DABS(TRUE(3,JSTART))
DPHI = FKDPHI(SSMT(5,JSTART),TRUE(5,JSTART))
CALL HFILL(230,FLOAT(NPS+NRS),D1OP,1.)
CALL HFILL(231,FLOAT(NPS+NRS),D1OP,1.)
CALL HFILL(232,FLOAT(NPS+NRS),D1OP,1.)
CALL HFILL(233,PLOG,D1OP,1.)
CALL HFILL(234,PLOG,DTHETA,1.)
CALL HFILL(235,PLOG,DPHI,1.)
CALL HFILL(243,THETA,D1OP,1.)
CALL HFILL(244,THETA,DTHETA,1.)
CALL HFILL(245,THETA,DPHI,1.)
CALL HFILL(253,PHI,D1OP,1.)
CALL HFILL(254,PHI,DTHETA,1.)
CALL HFILL(255,PHI,DPHI,1.)
CALL HFILL(IBASE+1+5*IRP(KPL),SNGL(RSMT(1,KPL)),0.,1.)
& CALL HFILL(IBASE+12,SNGL(RSMT(2,KPL)),0.,1.)
CALL HFILL(IBASE+5+5*IRP(KPL),SNGL(CHISMT(KPL)),0.,1.)
CALL HFILL(KPL+500,CHPROB,0.,1.)
CALL HFILL(KPL+600,CHPROB,0.,1.)
CALL HFILL(KPL+700,CHPROB,0.,1.)
CALL HFILL(KPL+800,SNGL(RSMT(1,KPL)),0.,1.)
& CALL HFILL(KPL+900,SNGL(RSMT(2,KPL)),0.,1.)
*