*-- Author : I.O. Skillicorn SUBROUTINE FTLSEG(NSS,IM) **: FTLSEG 40000 SM. New monitoring histos. **---------------------------------------------------------------------- C AUTHOR: I.O.SKILLICORN *KEEP,FRDIMS. *KEEP,FH1WORK. * Planar geometry * * Radial geometry * * Radial data... * * Planar Data * * Pointers into DIGI bank for IOS labelled hits * * Track segment data * * Fit data * * *KEEP,FSGPAR. *KEEP,FPTVTX. **the common/VERTEX/ becomes /VERTVV/ (in analogy to /VERTFF/) on the ** 17/6/91, since it is in conflict with the VERTEX module (g.bernardi) ** (note that all these common names should start by F in this deck...) *KEND. * LOCAL ARRAYS * C MODS 13/5/91 FOLLOW TO SAVE SPACE VVVVVVVVVVVVVVVVVVVVVVVVV C ADD C CEND ADDITION ******************************************** * * > Chard * Hard-wired cuts if wanted, else take from COMMON/FSGPAR/ filled * by FPTINT from the FRCP bank. * C MINIMUM SIZE OF CLUSTER FOR STARTING TRIPLE FINDING Chard MINHTS=3 C MINIMUM NUMBER OF POINTS/TRACK SEGMENT Chard MINPTS=4 * but keep this to avoid array size problems!! C MAX SIZE OF CLUSTER FOR ANALYSIS C Chard NWIRES= 8 C C C C*********************************************************************** C IPLOT=0 FOR NO DIAGNOSTIC T0, RMS TO LINE SEG, PLOTS: MAX SPEED Chard IPLOT=0 C*********************************************************************** C C WRITE(*,*)' NDP ',NDP C DO 4999 IP=1,36 C DO 5000 I=1,NDP(IP) C WRITE(*,*)' DRI,RM,NW',DRI(I,IP),RM(I,IP),NW(I,IP),I,IP C000 CONTINUE C999 CONTINUE C ENSURE ONLY 1 HIT/WIRE - TAKE MIN DRIFT HIT C SAME WIRE - KILL LONGER DRIFT C C SUM WIRE HITS - ALL PLANES IN Z C CLUSTER SIZE C I=NW(J,IP)+1 C CHANGE 9/12/92 C NOW HAVE STARTING CLUSTERS C IWC= # HITS/WIRE C WRITE(*,*)' CLUSTERS ' C PRINT 1000,IWC C CALL HFILL(70001,FLOAT(KK)+0.1,0.,1.) C IRN = HIT # CDEB PRINT 1001,IM,I,(IRN(K,I),K=1,KK) C IPN = PLANE # C PRINT 1002,IM,I,(IPN(K,I),K=1,KK) C SEARCH FOR LARGEST CLUSTER C IF(IC.NE.0)CALL HFILL(70002,FLOAT(MAXC)+0.1,0.,1.) C WRITE(*,*)' CLUSTER FOUND IC,MAXC ',IC,MAXC C CHECK FOR SUFFICIENT POINTS FOR LARGEST CLUSTER C FINISHED FOR THIS MODULE CDEB WRITE(*,*)' ANALYSE FOUND CLUSTER IC,MAXC ',IC,MAXC C CHECK MAX SIZE OF CLUSTER * * * DIAGNOSTIC PLOTS FOR LOW MULTIPLICITY CLUSTERS ONLY C CALL FTLISA(IRN(1,1),IPN(1,1),IWC(1),IC,MAXC,NSS,IM,IDIAG) C * * * Make Tzero, Resolution and twin-peaks checksums. C PLANE # C POINT # C STORE FIRST WIRE NUMBER AND DRIFT SIGN C KEEP TO SAME WIRE FOR DRIFT RESIDUAL C SIGN CHANGE - STORE WIRE # AT FIRST CROSS C T0 DETERMINATION - 8 WIRES C CHECK 4 WIRES C +++- AND ---+ C SAME WIRE * T-zero... C +--- AND -+++ * T-zero... C ++-- AND --++ C CALL HFILL(90005,TZV,0.,1.) C IF(TZV.GT.0.0)CALL HFILL(90015,TZV,0.,1.) C ++++ AND ---- * * Twin peaks... * Resolution... C T0 DETERMINATION - 8 WIRES C CHECK 5 WIRES TO ENSURE SAFE SIGN CHANGE 13/8/91 C +++-- AND ---++ C SAME WIRE C CALL HFILL(90004,TZ,0.,1.) C ++--- AND --+++ C CALL HFILL(90004,TZ,0.,1.) C NO SIGN CHANGE C IF(LL.GE.4.AND.NCR.EQ.0)CALLHFILL(90002,0.01,0.,1.) C SIGN CHANGE C IF(LL.GE.4.AND.NCR.NE.0)CALLHFILL(90002,WCROSS,0.,1.) C FIT STR. LINE TO SEGMENT IF 4 OR MORE POINTS CALL FTLFT(XX,YY,LL,0,AT,BT,EE) CALL FTLFT(XX,YN,LL,0,AR,BR,ER) C WRITE(*,*)' EE,ER',EE,ER C RESIDUAL PLOT 18/7/92 C CALL HFILL(90000,EE,0.,1.) C CALL HFILL(90001,ER,0.,1.) C RMS TO LINE SEGMENT DRIFT,RADIUS C WRITE(*,*)' AT,XX(LL/2),BT',AT,XX(LL/2),BT,EE,ER C DRIFT AT CENTRE OF LINE SEG * * * Small farm mod by ? * add SEQ BOSMDL so above might work * initialise RR (else DDMIN*RR test might fail? * C 21/01/93 301211908 MEMBER NAME FTLISAA (FILE46) FVS C 06/01/93 301081132 MEMBER NAME FTLISA5 (FILE46) FVS C 08/12/92 301061503 MEMBER NAME FTLISA2 (FILE46) FVS C 06/12/92 212081044 MEMBER NAME FTLISA (GRAPHICS) FVS C 04/12/92 212061354 MEMBER NAME LISALIB (FILE46) FVS C 28/11/92 212011534 MEMBER NAME FTLISA (FILE46) FVS C 28/11/92 211281327 MEMBER NAME FTLISA (GRAPHICS) FVS C 13/11/92 211281108 MEMBER NAME FTLSEGO (GRAPHICS) FVS C 13/11/92 MEMBER NAME FTLSEGO (FILE46) FVS C 04/12/91 MEMBER NAME NEWSEG (FILE46) FVS C FROM FILE46(GRFTRAC2) C C 13/11/92 C ADD SLOPE CHECK ON TRIPLES - SELECT MINIMUM SLOPE C REMOVE EXTRAPOLATION IN PHI C RCUT ---> 20 CMS C SLCUT---> 0.1 C TSCUT---> 0.08 (WAS 0.1) C DMINX---> 0.15 (WAS 0.2) C C PROJECTION CUT WITH DMINX C VERSION FOR TRACK FINDING IN CLUSTER BY TRIPLES