+PATCH,$CORR. +DECK,CORR1. Updates version: 1.27/00 to 1.28/00 +REP,*TITLE*,TITLE,1. CBOFF 1.28/00 10/06/95 01.26.54 SE +REP,$VERSION,V1_27,1-2. *CMZ : 1.27/06 09/06/95 19.34.38 by Christian Voelcker *CMZ : 1.27/05 08/05/95 16.49.59 by Christian Voelcker *CMZ : 1.27/04 26/04/95 20.01.32 by Christian Voelcker *CMZ : 1.27/03 26/04/95 11.31.30 by Christian Voelcker *CMZ : 1.27/02 18/04/95 16.02.11 by Christian Voelcker *CMZ : 1.27/01 08/04/95 12.12.32 by Christian Voelcker *CMZ : 1.27/00 07/04/95 19.39.20 by Christian Voelcker *-- Author : Christian Voelcker 07/04/95 * *::> VERSION 1.27/06 09/06/95 19.34.38 * *> new common in TC_MACRO to hold primary * pattern recognition logicals. TCFFUZ becomes * default now, 'FFUZ' and 'RADO' keywords * moved to the jdc card PREC. (see locater) * *::> VERSION 1.27/05 08/05/95 16.49.59 * *> new key for TC_FUZZ, therefore TC_MACRO modified: * new variable KY30TC in /TCFZZZ/. * *::> VERSION 1.27/04 26/04/95 20.01.32 * *> some modifications in dec CFCOMMON. reorganized the * /SQUIDC/ common, one variable got appended to the * /CFVARS/ common. * *::> VERSION 1.27/03 26/04/95 11.31.30 * *> Add two more key variables KY47TC,KY48TC to /TCFZZZ/ * in deck tc_macro for Fast-FUZzy. * *::> VERSION 1.27/02 18/04/95 16.02.10 * *> Add Fast-Fuzzy common blocks in TC_MACRO (C.A.Meyer) * *::> VERSION 1.27/01 08/04/95 12.12.32 * *> more modifications to the /CFVARS/ common in * deck CFCOMMON. C-programmmers will jubilate !! +ADD,$VERSION,V1_28,*. &DECK,V1_28. *CMZ : 1.28/00 10/06/95 01.26.56 by Christian Voelcker *-- Author : Christian Voelcker 10/06/95 * *::> VERSION 1.28/00 10/06/95 01.26.54 * * implement the silicon vertex - first * iteration... * be aware that some of the common blocks * have changed! (CBLINK, CBXDIV, ...) * +ADD,COMMCB,CBCOMMON,0. *CMZ : 1.28/00 10/06/95 00.21.02 by Christian Voelcker +ADD,COMMCB,CBCOMMON,65. * 06-Jun-95 Chr. Voelcker * * 06-Jun-95 Chr. Voelcker * modified four common blocks: * /CBLINK/,/CBXDIV/,/CBDBCO/,/CBHARD/ * as we introduce three more zebra and a new calibration * type for the silicon vertex detector. * Only for the last of these commons we break the * rule to add new words always at the end of the common. * C-programmers won't need /CBDBCO/.... * one new common block CBTBSV for silicon vertex * tables. * +REP,COMMCB,CBCOMMON,127. & LLAST, & LRVTX,LRSVX,LTVXT +REP,COMMCB,CBCOMMON,147. & IXSTOR, & IXRVTX,IXRSVX,IXTVXT +REP,COMMCB,CBCOMMON,165. & LLAST, & LRVTX,LRSVX,LTVXT +REP,COMMCB,CBCOMMON,185. & IXSTOR, & IXRVTX,IXRSVX,IXTVXT +REP,COMMCB,CBCOMMON,331. * GTPWCC,GTVTCC,GTLGCC,GTLPCC,GTLXCC, +REP,COMMCB,CBCOMMON,337. * TBPWCC,TBVTCC,TBLGCC,TBLPCC,TBLXCC, +REP,COMMCB,CBCOMMON,343. * IMPWCC,IMVTCC,IMLGCC,IMLPCC,IMLXCC, +REP,COMMCB,CBCOMMON,350. * GTPWCC,GTVTCC,GTLGCC,GTLPCC,GTLXCC, +REP,COMMCB,CBCOMMON,356. * TBPWCC,TBVTCC,TBLGCC,TBLPCC,TBLXCC, +REP,COMMCB,CBCOMMON,362. * IMPWCC,IMVTCC,IMLGCC,IMLPCC,IMLXCC, +REP,COMMCB,CBCOMMON,496-498. > RDJDC, RDPWC, RDSVTX LOGICAL RDRSV1, RDSLOW, RDLGHT, RDSOFT, RDBC, > RDJDC, RDPWC, RDSVTX +ADD,COMMCB,CBCOMMON,508. * RDSVTX : silicon vertex detector +ADD,COMMCB,CBCOMMON,606. &KEEP,CBTBSV. INTEGER ISPEDS(0:1919),ISTHRE(0:1919) COMMON /CBTBSV/ ISPEDS,ISTHRE SAVE /CBTBSV/ +REP,COMMCB,CFCOMMON,1-5. *CMZ : 1.27/04 26/04/95 20.00.57 by Christian Voelcker ! * copy the latest cfcommon to the cboff.cmz file ! * /SQUIDC/ got re-arranged, some other things changed as well! *CMZU: 3.08/00 21/04/95 11.21.51 by Pal Hidas *CMZU: 3.05/01 29/03/95 15.35.29 by Pal Hidas *CMZ : 3.05/00 25/03/95 19.03.01 by Pal Hidas *CMZ : 3.04/02 07/03/95 17.23.04 by Pal Hidas *CMZ : 3.04/00 04/03/95 16.36.26 by Pal Hidas *CMZ : 3.03/00 01/03/95 18.20.12 by Pal Hidas *CMZU: 3.00/07 11/10/94 17.51.04 by Pal Hidas +REP,COMMCB,CFCOMMON,63. PARAMETER (NCMAX=20,NPMAX=25,NRMAX=10,NVMAX=81) +REP,COMMCB,CFCOMMON,89-102. I ,NPOS,NNEG,NNE,IEVEF,CFTAKH,BADTRA,BADNEU,NPI0CF I ,MERPI0,BADKL,IDRVCF,MNVXCF,IMNNEU,KSKBCF,BADXCF R ,CHISQ,CUTCL,CFBEAM,CFMISM,CFVTER,CFMIER,CFMISX L ,MISGAM,MISSKL,MISNEU,MISPRO,CFTRHY,CFTRSU,CFDBUG L ,FPUTZ,CFVERT,CFVRTZ,CFSUPC,CFDEUT,CFORCE,CFGCOV L ,CFIGNR,CFVERZ,CFVERA,CFKS00,CFCHCR,CFKSCC L ,CFCOLL,CFKSKM,CFKSKS SAVE /CFVARS/ INTEGER NXTCMB(NPMAX),MEASCO(NPMAX),EVTYPE(6),CFCODE(NCDMAX) INTEGER NPOS,NNEG,NNE,IEVEF,CFTAKH,BADTRA,BADNEU,NPI0CF INTEGER MERPI0(NRMAX),BADKL,IDRVCF(10),MNVXCF,IMNNEU(NPMAX) INTEGER KSKBCF,BADXCF REAL CHISQ,CUTCL,CFBEAM,CFMISM,CFVTER(3),CFMIER(6),CFMISX LOGICAL MISGAM,MISSKL,MISNEU,MISPRO,CFTRHY,CFTRSU,CFDBUG LOGICAL FPUTZ,CFVERT,CFVRTZ,CFSUPC,CFDEUT,CFORCE,CFGCOV LOGICAL CFIGNR,CFVERZ,CFVERA,CFKS00,CFCHCR,CFKSCC LOGICAL CFCOLL,CFKSKM,CFKSKS +ADD,COMMCB,CFCOMMON,121. * CFVTER : ERRORSQ ON VERTEX "MEASURMENT" (X,Y,Z) * CFMIER : "COVARIANCE MATRIX ON MISSING PARTICLE" +ADD,COMMCB,CFCOMMON,133. * CFKSCC : LOGICAL FLAG TO FIT XKSKS,KS -> PI+,PI- BOTH +ADD,COMMCB,CFCOMMON,141. * BADKL : TTKS NUMBER OF KLONG WITH BAD ENERGY +REP,COMMCB,CFCOMMON,278-286. C COMMON /SQUIDC/ X,VAR,ERR,VVAR,DPDX,DFDX,F,GF,H,PFIT,DELX, C * X0,VERT,BQRS,DFDVER,DD,FF,BQRS0,KDIRS,DFDY,GH, C * SSS,PULLD,HHH,KSVT,CHVT,CHVMOM,CHVPAR,CHKMOM, C * CFLY,CHCTAU,XM,ETOT,MASS,MRES,CFVENE,CFVMOM, C * NC,NVAR,KTYPE,MISSID,CASE,NRES,MISSI0, C * NPART,MEM,PMEM,RC,ISPE, C * NCH,NSUM,ITER, C * ENCONS COMMON /SQUIDC/ X,VAR,ERR,VVAR,DFDY,DPDX,DFDX,DELX,DD,FF, D F,GF,PFIT,GH,H,X0,VERT,BQRS,DFDVER,BQRS0, D KSVT,CHVT,KDIRS,CHVMOM,CHVPAR,CHKMOM,CFLY, D CHCTAU,SSS,PULLD,HHH,EN4,PPP4, R XM,ETOT,MASS,MRES,CFVENE,CFVMOM, I NC,NVAR,KTYPE,MISSID,CASE,NRES,MISSI0,NPART, I MEM,PMEM,RC,ISPE,NCH,NSUM,ITER, L ENCONS C +REP,COMMCB,CFCOMMON,296-299. DOUBLE PRECISION BQRS0(3,NPMAX),KSVT(3,2),CHVT(3),KDIRS DOUBLE PRECISION CHVMOM(3,4),CHVPAR(3,4),CHKMOM(3) DOUBLE PRECISION CFLY,CHCTAU DOUBLE PRECISION SSS(3,3),PULLD(3),HHH(NCMAX,NCMAX) DOUBLE PRECISION EN4(4),PPP4(4) +REP,COMMCB,CFCOMMON,335-336. C 2 => NEUTRAL C 3 => MISSING CHARGED C 4 => MISSING NEUTRAL C 5 => MISSING (OLD METHOD,CASE=2,3) C 6 => CHARGED WITH MISSING ENERGY C 7 => NEUTRAL WITH MISSING ENERGY +REP,COMMCB,CFCOMMON,347-351. C MISSING (SQRT(E)) C 6 => Z-VERTEX FIT, ALL NEUTRALS, NOTHING C MISSING (PHI,SQRT(E),LAMBDA) C 7 => VERTEX FIT, ALL NEUTRALS, NOTHING C MISSING (PHI,SQRT(E),LAMBDA) C 8 => KSHORT VERTEX FIT, FIXED MAIN VERTEX C NOTHING MISSING C 9 => KSHORT VERTEX FIT, FIXED MAIN VERTEX C MISSING KLONG C 10 => KSHORT VERTEX FIT, FIXED MAIN VERTEX C BAD KLONG (ENERGY MISSING) C 11 => MISSING 3-MOMENTA, MISSING MASS IS C KNOWN C 12 => MISSING ENERGY, MISSING MASS C AND DIRECTION KNOWN +ADD,COMMCB,CFCOMMON,359. C (OLD METHOD, CASE=3) +REP,COMMCB,TC_MACRO,1. *CMZ : 1.28/00 09/06/95 21.38.04 by Christian Voelcker *CMZ : 1.27/06 08/06/95 18.53.59 by Christian Voelcker *CMZ : 1.27/05 08/05/95 15.46.11 by Christian Voelcker *CMZ : 1.27/03 26/04/95 11.30.14 by Christian Voelcker *CMZU: 1.27/02 10/04/95 15.36.58 by Curtis A. Meyer *CMZU: 1.27/00 06/04/95 17.02.28 by Christian Voelcker +REP,COMMCB,TC_MACRO,30. * Add OLDJDC flag to TCFLAG, and OLDJLG to LGHOLD. +REP,COMMCB,TC_MACRO,83-87. * 06/04/95 Chr. Voelcker * Prepare already for the next release, which will include the * Fast Fuzzy pattern recognition: new variable FFUZTC in the * /TCFLAG/ common. * * 10/04/95 C.A.Meyer * Add Fast-Fuzzy common blocks. * * 26/04/95 Chr. Voelcker * Add two more key variables KY47TC,KY48TC to /TCFZZZ/ * * 08/06/95 Chr. Voelcker * fast fuzzy shall become default now. * move the logicals RADOTC and FFUZTC from common /TCFLAG/ * to a new common /TCPREC/. * keep logical RASHTC at the same place as before * (C-programmers..). Introducing GVTXTC as logical for * silivertex reconstruction. * * 09/06/95 R.Ouared/Chr. Voelcker * new arrays added in TPPRMS common to take into account * the silicon vertex detector as well as in TCLIFT,TRKPRM * TCFLAG. * MTVXT(5) and MRSVX(5) added to /TCLIFT/ * create new commons /TCSVXY/, /TCSVTH/, /TCSVTB/, /TCSVOF/, * /TCSVRA/ and /TCSVCL/ to steer svtx parameters via input * cards. * +REP,COMMCB,TC_MACRO,167. INTEGER MTJSL(5),MTJSR(5),MTJSZ(5),MTVXT(5),MRSVX(5) +REP,COMMCB,TC_MACRO,174. & ,MTJSL,MTJSR,MTJSZ,MTVXT,MRSVX +REP,COMMCB,TC_MACRO,180-182. INTEGER LENSV,LENCLS,LENBP PARAMETER (LENTJ=16,LENHT=23,LENTK=30,LENTR=33,LENHP=6) PARAMETER (LENVX=14,LENVP=13,LENJCP=20) PARAMETER (LENPW=10,LENCL=12,LENVH=7) PARAMETER (LENSV=2, LENCLS=16,LENBP=2) +REP,COMMCB,TC_MACRO,205-208. * ASVXTP contains pitch,half depth and slope of each SVTX module. * OSVXTP contains x-coordinate and y-coordinate of first strip in * first module. * XSVXTP contains the x-ccordinate of each STRIP in the SVTX. * YSVXTP contains the y-coordinate of each STRIP in the SVTX. * RSVXTP contains the radius of the STRIPS in the FIRST SVTX. * FSVXTP contains phi for all the strips * REAL RPWCTP(2),APWCTP(2),OPWCTP(2) REAL XPWCTP(0:319),YPWCTP(0:319) REAL RSVXTP(0:127),ASVXTP(3),OSVXTP(2) REAL XSVXTP(0:1919),YSVXTP(0:1919),FSVXTP(0:1919) COMMON /TPPRMS/ RPWCTP,APWCTP,OPWCTP,XPWCTP,YPWCTP, & RSVXTP,ASVXTP,OSVXTP,XSVXTP,YSVXTP,FSVXTP +REP,COMMCB,TC_MACRO,329-338. * GVTXTC Activate routines for the silicon vertex detector. * RASHTC Enables the special fuzzy Radon patt. rec. for short tracks * LOGICAL GPWCTC,TRAKTC,RAWSTC,PATTTC,CIRCTC LOGICAL HELXTC,VERTTC,VER2TC,OLDJDC,RTRKTC LOGICAL GVTXTC,RASHTC COMMON /TCFLAG/ GPWCTC,TRAKTC,RAWSTC,PATTTC,CIRCTC, & HELXTC,VERTTC,VER2TC,OLDJDC,RTRKTC, & GVTXTC,RASHTC +ADD,COMMCB,TC_MACRO,753. &KEEP,TCCOMM. * INTEGER IAAA(1200,30),NPOINTS(2) COMMON /TCCOMM/ IAAA,NPOINTS SAVE /TCCOMM/ * &KEEP,TCFZZZ. * * Control cards for the Fast-Fuzzy package. * REAL KY20TC,KY22TC,KY35TC,KY46TC,KY47TC,KY48TC INTEGER KY25TC,KY30TC COMMON /TCFZZZ/ KY20TC,KY22TC,KY35TC,KY46TC,KY47TC,KY48TC, & KY25TC,KY30TC SAVE /TCFZZZ/ * &KEEP,TCPREC. * * PRECTC controls, which primary Pattern RECognition * should be activated: * FFUZTC Fast FUZzy pattern recognition. (default) * RADOTC fuzzy RADOn pattern recognition. * OLPTTC OLd PaTtern recognition. * The variables can be steered using the PREC card * in the jdc input file (logical unit 81). * The secundary pattern recognition (RASH) will remain as * option to be selected in cboff card CHAM * * INTEGER PRECTC LOGICAL FFUZTC,RADOTC,OLPTTC COMMON /TCPREC/ PRECTC, & FFUZTC,RADOTC,OLPTTC SAVE /TCPREC/ * &KEEP,TCSVXY. * * POSITIONNING STRIP 128 IN FIRST MODULE . * REAL FSTRIP(2) COMMON /TCSVXY/ FSTRIP SAVE /TCSVXY/ * &KEEP,TCSVTH. * * pedestal thresholds (sigma above noise) * REAL SVSIGM(2) COMMON /TCSVTH/ SVSIGM SAVE /TCSVTH/ * &KEEP,TCSVTB. * * external pedestal tables * LOGICAL LSVTBX COMMON /TCSVTB/ LSVTBX SAVE /TCSVTB/ * * &KEEP,TCSVOF. * * RVTX and RTVX word offsets. * INTEGER IOFSET(2) COMMON /TCSVOF/ IOFSET SAVE /TCSVOF/ * &KEEP,TCSVRA. * * has to see with SVTX signal (raw,ped-sub,comon-mode corrected) * REAL AMPLS(1920),PAMPL(1920),CAMPL(1920),CMOD(15) INTEGER NSHORT,NFORM LOGICAL LSUBPED COMMON/TCSVRA/ AMPLS,PAMPL,CAMPL,CMOD,NSHORT,NFORM,LSUBPED SAVE /TCSVRA/ * &KEEP,TCSVCL. * * defines minimum gap between clusters. * INTEGER ISGAP COMMON /TCSVCL/ ISGAP SAVE /TCSVCL/ * +ADD,CBPHYS,CBFFGO,0. *CMZ : 1.28/00 09/06/95 23.43.07 by Christian Voelcker *CMZ : 1.27/06 08/06/95 18.38.41 by Christian Voelcker +ADD,CBPHYS,CBFFGO,92. * June 6th, 1995 Chr. Voelcker * *> keywords FFUZ and RADO migrate to * jdc card PREC. * *> new keyword in CHAM card: GVTX * * * June 6th, 1995 R.Ouared * calibration: GTVTCC is false by default * +DEL,CBPHYS,CBFFGO,266. +REP,CBPHYS,CBFFGO,330-335. GVTXTC = .FALSE. RASHTC = .FALSE. ELSE TRAKTC = IUCOMP(4HTRAK,ITRAK,NTRAK) .NE. 0 GPWCTC = IUCOMP(4HGPWC,ITRAK,NTRAK) .NE. 0 GVTXTC = IUCOMP(4HGVTX,ITRAK,NTRAK) .NE. 0 +REP,CBPHYS,CBFFGO,343-345. RASHTC = IUCOMP(4HRASH,ITRAK,NTRAK) .NE. 0 +REP,CBPHYS,CBFFGO,357-361. IF(GPWCTC.AND.GVTXTC) THEN WRITE(LLOG,1120) 1120 FORMAT(' cannot select both PWC and SVTX detector.', > /,' Will select PWC by default. ' > ,'This may change in USRUN') GVTXTC = .FALSE. +ADD,CBPHYS,CBFFGO,532. * SVTX calibration ... not yet used ... JX= ICLOC('SV',2,CCAL,1,80) IF ( JX .GT.0 ) THEN IF ( (ICLOC('-',1,CCAL,JX+2,JX+4)-JX) .GT. 0 ) THEN GTVTCC = .FALSE. ENDIF ENDIF +ADD,CBPHYS,CBINIT,0. *CMZ : 1.28/00 09/06/95 20.44.19 by Christian Voelcker +ADD,CBPHYS,CBINIT,65. * Modified 6 June, 1995 Rafik R. * new variables for svtx to initialize. +ADD,CBPHYS,CBINIT,206. GTVTCC = .FALSE. ! set .TRUE. if calibration is availiable +ADD,CBPHYS,CBINIT,233. TBVTCC = .TRUE. +ADD,CBPHYS,CBINIT,263. IMVTCC = 0 +ADD,CBPHYS,CBINIT,445. IXTVXT = IXEVNT +ADD,CBPHYS,CBINIT,457. IXTVXT = IXEVNT +ADD,CBPHYS,CBSRUN,0. *CMZ : 1.28/00 10/06/95 01.19.28 by Christian Voelcker +REP,CBPHYS,CBSRUN,371. RDRSV1 = JBIT(NEWREA, 8).NE.0 RDSVTX = JBIT(NEWREA, 7).NE.0 +REP,CBPHYS,CBSRUN,380-404. READOUT(LEN:LEN+2) = '?? ' LEN=LEN+3 ENDIF IF(RDSVTX) THEN READOUT(LEN:LEN+2) = 'SV ' LEN=LEN+3 ENDIF IF(RDSLOW) THEN READOUT(LEN:LEN+2) = 'SC ' LEN=LEN+3 ENDIF IF(RDLGHT) THEN READOUT(LEN:LEN+2) = 'LP ' LEN=LEN+3 ENDIF IF(RDSOFT) THEN READOUT(LEN:LEN+2) = 'SF ' LEN=LEN+3 ENDIF IF(RDBC) THEN READOUT(LEN:LEN+2) = 'BC ' LEN=LEN+3 ENDIF IF(RDJDC) THEN READOUT(LEN:LEN+2) = 'JD ' LEN=LEN+3 ENDIF IF(RDPWC) THEN READOUT(LEN:LEN+1) = 'PW ' +ADD,CBPHYS,CBTABS,0. *CMZ : 1.28/00 09/06/95 20.56.47 by Christian Voelcker +ADD,CBPHYS,CBTABS,61. * Jun-95 CV/RO call to CBCASV to decode svtx tables * +REP,CBPHYS,CBTABS,71. &SEQ,CBTBSV. &SEQ,TCSVTB. * INTEGER JRTVX COMMON /ADRTVX/JRTVX SAVE /ADRTVX/ +REP,CBPHYS,CBTABS,455. * * locate table bank RTVX for silicon vertex detector * IF(.NOT.LSVTBX) THEN JRTVX= LZFIDH(IXHRAW,'RTVX',0) IF(JRTVX.GT.0) THEN CALL CBCASV(JRTVX) ENDIF ENDIF * +ADD,CBPHYS,GETEVT,0. *CMZ : 1.28/00 09/06/95 21.03.25 by Christian Voelcker +ADD,CBPHYS,GETEVT,62. * * Updated: 26 April, 1995 R.Ouared * addresses of new svtx banks have been allocated: * raw bank: LRVTX=LQ(LHRAW-15) * zero suppressed, pedestal substracted, * common mode rejected raw data: * LRSVX=LQ(LHRAW-16) * reconstruction bank: LTVXT=LQ(LHTJD-8) * +ADD,CBPHYS,GETEVT,258. LRVTX= LQ(LHRAW-15) LRSVX= LQ(LHRAW-16) +ADD,CBPHYS,GETEVT,314. LTVXT= LQ(LHTJD-8) +ADD,CBPHYS,GET_EVEN,0. *CMZ : 1.28/00 09/06/95 20.59.21 by Christian Voelcker +REP,CBPHYS,GET_EVEN,47-51. * 'RXDF','RSFT','RLPR','RLAT','RVTX','RSVX', * 3*' '/ DATA HEAD/'EVHD'/,NBANK/15/ DATA ILOC/-20,-1,-2,-3,-4,-4,-5,-6,-7,-8,-9,-10,-11,-12,-13, * -14,-15,-16,3*-20/ +ADD,CBPHYS,GET_EVEN,342. LRVTX= LQ(LHRAW-15) LRSVX= LQ(LHRAW-16) +ADD,CBPHYS,ZBINIT,0. *CMZ : 1.28/00 09/06/95 23.52.33 by Christian Voelcker +REP,CBPHYS,ZBINIT,46-48. * Modified : 9 June, 1995 R.Ouared * new variables to implement svtx: * 'RVTX', 'RSVX' in CH4 array * IXRVTX,IXTVXT,IXRSVX * * References: ZEBRA Manuals. * * Description: This routine was stripped out of the original +REP,CBPHYS,ZBINIT,102-107. * 'RENO','RVTX','RSVX', 2*' '/ DATA HEAD/'EVHD'/,IBANK/16/ DATA JLOC/-20, -1, -2, -3, -4, -5, * -6, -7, -8, -9, -3, * -4,-11,-12,-13,-19, * -18,-15,-16, 2*-20/ +ADD,CBPHYS,ZBINIT,219. IXRVTX = IXHRAW IXRSVX = IXHRAW +ADD,CBPHYS,ZBINIT,227. IXTVXT = IXHTJD +ADD,CBPHYS,CBCASV,*. &DECK,CBCASV. *CMZ : 1.28/00 09/06/95 21.08.37 by Christian Voelcker *CMZ : 1.30/00 25/05/95 16.28.15 by RAFIK OUARED *CMZU: 1.27/06 20/05/95 19.25.37 by RAFIK OUARED *-- Author : RAFIK OUARED 15/05/95 SUBROUTINE CBCASV(LINKTS) * *********************************************************************************** * * Creation date: 15 May, 1995 * * References : * * Description: This routine take the svtx table bank RTVX and fill up common * with pedestal and threshold information. * * Called by : CBTABS * * Subroutines referenced: JBYT,JBIT * * *********************************************************************************** &SELF,IF=-F77,IF=-SUN. IMPLICIT NONE &SELF,IF=-F77,IF=SUN. IMPLICIT UNDEFINED (A-Z) &SELF. * * PASSED VARIABLES * LINKTS : LINK TO THE TABLE BANK RTVX. * INTEGER LINKTS * INTEGER I,J INTEGER NFORM,NSHORT,IBITMA,NDET,NSTR,IPED1,IPED2,IOFF,IOFB2 * INTEGER JBYT,JBIT EXTERNAL JBYT,JBIT * &SEQ,CBLINK. &SEQ,TCLIFT. &SEQ,TPPRMS. &SEQ,PI2PI. &SEQ,TRKPRM. &SEQ,CBUNIT. &SEQ,CBTBSV. &SEQ,USKEYS. &SEQ,TCSVOF. * * INITIALIZE CBTBSV. * DO I=0,1919 ISPEDS(I)=0 ISTHRE(I)=0 ENDDO * IOFF=IOFSET(2) IOFB2=IOFF/2 * WRITE(LLOG,*) 'CBCASV: LINKTS = ',LINKTS NFORM=JBYT(IQ(LINKTS+1),17,16) NSHORT=JBYT(IQ(LINKTS+1),1,16) IBITMA=IQ(LINKTS+2) NDET=JBYT(IQ(LINKTS+3),17,16) NSTR=JBYT(IQ(LINKTS+3),1,16) IPED1=JBYT(IQ(LINKTS+4+IOFB2),17,16)/16 IPED2=JBYT(IQ(LINKTS+4+IOFB2),1,16)/16 &SELF,IF=CBCASV. WRITE(LLOG,*) 'CBCASV: FORMAT =',NFORM WRITE(LLOG,*) 'CBCASV: NSHORT =',NSHORT WRITE(LLOG,*) 'CBCASV: BITMASK =',IBITMA WRITE(LLOG,*) 'CBCASV: NUMBER OF DETECTORS =',NDET WRITE(LLOG,*) 'CBCASV: NUNMBER OF STRIP PER DET =',NSTR WRITE(LLOG,*) 'CBCASV: FIRST PED =',IPED1 WRITE(LLOG,*) 'CBCASV: SECOND PED =',IPED2 &SELF. * FILLUP COMMON DO I=0,1919 J=I+7+IOFF IF(2*(J/2).NE.J) THEN ISPEDS(I)=JBYT(IQ(LINKTS+(J+1)/2),17,16)/16 ISTHRE(I)=JBYT(IQ(LINKTS+1024+(J+1)/2),17,16)/16 ISTHRE(I)=ISTHRE(I)-ISPEDS(I) ISTHRE(I)=NINT(FLOAT(ISTHRE(I))/5.) ELSE ISPEDS(I)=JBYT(IQ(LINKTS+(J+1)/2),1,16)/16 ISTHRE(I)=JBYT(IQ(LINKTS+1024+(J+1)/2),1,16)/16 ISTHRE(I)=ISTHRE(I)-ISPEDS(I) ISTHRE(I)=NINT(FLOAT(ISTHRE(I))/5.) ENDIF &SELF,IF=CBCASV. WRITE(LLOG,*) 'CBCASV: STRIP,PED,THRES...' WRITE(LLOG,*) I,ISPEDS(I),ISTHRE(I) &SELF ENDDO * RETURN END