+PATCH,$CORR. +DECK,CORR1. Updates version: 3.09/00 to 3.10/00 +REP,*TITLE*,TITLE,1. * CBKFIT 3.10/00 10/09/96 12.20.55 +ADD,CFCOMMON,CFCOMMON,0. *CMZ : 3.10/00 10/09/96 10.58.53 by Unknown +REP,CFCOMMON,CFCOMMON,98. L ,CFCOLL,CFKSKM,CFKSKS,CFCHCV +REP,CFCOMMON,CFCOMMON,108. LOGICAL CFCOLL,CFKSKM,CFKSKS,CFCHCV +ADD,CFCOMMON,CFCOMMON,155. * CFCHCR : take momenta from TVVP * CFCHCV : take also cov.matrix from TVVP (-> error prop.) +REP,CFCOMMON,CFCOMMON,315. DOUBLE PRECISION CHVMOM(3,4),CHVPAR(9,4),CHKMOM(3,2) +ADD,CBKFIT,CFNSQD,0. *CMZ : 3.10/00 10/09/96 11.33.01 by Unknown +REP,CBKFIT,CFNSQD,58-63. DOUBLE PRECISION GGG(3,NCMAX) DOUBLE PRECISION JJ(NCMAX),WORKSP(99),VERT0(3) DOUBLE PRECISION HHH2(NCMAX,NCMAX),ICUTS,DELXN(NVMAX) DOUBLE PRECISION LAM REAL CHISPR INTEGER I,K,MU,NU,IDGT,IER,N,III,NNC,JJJ,NPAR1,IRC,NPAR2 C LAM=1.0D0 +ADD,CBKFIT,CFKINI,0. *CMZ : 3.10/00 10/09/96 11.27.22 by Unknown +REP,CBKFIT,CFKINI,76-82. DOUBLE PRECISION EPSZ DOUBLE PRECISION ZINCMX REAL EMISS,MM LOGICAL EM INTEGER EMR1,EMR2,JRES,M,IG,K,IC,I,II INTEGER J,IIP1,III,INMAX,ICP1,ICP2,IRC INTEGER NPAR1,NPAR2,ICL C EPSZ=0.2D0 ZINCMX=70.D0 +ADD,CBKFIT,CFSQPL,0. *CMZ : 3.10/00 10/09/96 11.34.35 by Unknown +REP,CBKFIT,CFSQPL,27-29. INTEGER IDIAG(3) DOUBLE PRECISION TTT(NCMAX,NCMAX),DVMPY,WORKSP(100) DOUBLE PRECISION UUU(NCMAX,NVMAX),TEMP1,TEMP2,TEMP(NCMAX) C IDIAG(1)=1 IDIAG(2)=3 IDIAG(3)=6 +ADD,CBKFIT,CFTVRS,0. *CMZ : 3.10/00 10/09/96 11.47.54 by Unknown +REP,CBKFIT,CFTVRS,15-34. CHARACTER*60 S(14) C CHARACTER*60 S(14) / C 1 ' ', C 2 ' ***********************', C 3 ' * *', C 4 ' * WELCOME TO *', C 5 ' * *', C 6 ' * C B K F I T *', C 7 ' * *', C 8 ' * *', C 9 ' * VERSION 1.00/00 *', C * ' * *', C 1 ' * *', C 2 ' * *', C 3 ' ***********************', C 4 ' ' / C &SEQ,VERSQQ. S( 1)=' ' S( 2)=' ***********************' S( 3)=' * *' S( 4)=' * WELCOME TO *' S( 5)=' * *' S( 6)=' * C B K F I T *' S( 7)=' * *' S( 8)=' * *' S( 9)=' * VERSION 1.00/00 *' S(10)=' * *' S(11)=' * *' S(12)=' * *' S(13)=' ***********************' S(14)=' ' C ICFVRS= IVERSQ CCFVRS= VERSQQ +DEL,CBKFIT,CFTVRS,41. +ADD,CBKFIT,CFUPDA,0. *CMZ : 3.10/00 10/09/96 11.49.39 by Unknown +REP,CBKFIT,CFUPDA,55. DOUBLE PRECISION FACTOR,PP,PXY,S,C,SSQRT,RHO(3),DEL(2) +REP,CBKFIT,CFUPDA,64-67. LOGICAL RES INTEGER EMR1,EMR2,JRES,M,N,IG,K,IC,LP,IR,NMEM,I,IM,IPF,IV,II,MU INTEGER J,III,JJJ,KTYMIS,NVAR6,NPAR1,INC,ITYP,ICHAR,IVAR,JVAR INTEGER NVAR1,ICC,NPAR2,IRET,IVERK,NPARA,NPARB,NVERK,IVC +ADD,CBKFIT,CFNSEC,0. *CMZ : 3.10/00 10/09/96 11.31.51 by Unknown +REP,CBKFIT,CFNSEC,29-35. DOUBLE PRECISION EPSK DOUBLE PRECISION ZINCMX REAL RMK LOGICAL RES INTEGER K,I,II,ITYP,ICHAR,JJ INTEGER J,INMAX,IRC,IW(3),IFAIL INTEGER NPAR1,NPAR2,ICL C EPSK=0.05D0 ZINCMX=70.D0 +ADD,CBHYPO,CFDOPE,0. *CMZ : 3.10/00 10/09/96 11.18.32 by Unknown +REP,CBHYPO,CFDOPE,51. INTEGER CFPICK +ADD,CBHYPO,CFDOHY,0. *CMZ : 3.10/00 10/09/96 11.17.25 by Unknown +REP,CBHYPO,CFDOHY,83. REAL ARESM,GODNES,PROB,GODBUF(100) +ADD,CBHYPO,CFDOFI,0. *CMZ : 3.10/00 10/09/96 11.16.59 by Unknown +REP,CBHYPO,CFDOFI,70-75. INTEGER IRET,KCP,KCN,KCH,KDE,KGA INTEGER NHM,NDWRES,I,IETYC,NDW,NSLSAV,NSL INTEGER NLI,J,IJUMP,JFI,JFJ,NIJ,INCND,JRET,K,RECRUN,RECEVE INTEGER III,JJJ +ADD,CBHYPO,CFEFFM,0. *CMZ : 3.10/00 10/09/96 11.20.13 by Unknown +REP,CBHYPO,CFEFFM,54. INTEGER JRES,IHY,NMRES,I,JMEAS,JTTKS +ADD,CBHYPO,CFZOUT,0. *CMZ : 3.10/00 10/09/96 11.53.04 by Unknown +REP,CBHYPO,CFZOUT,70-73. REAL PULL(NVMAX),GODNES,ENER,PX,PY,PZ INTEGER IHY,IG,IC,ID,IHYCOD,I,JBIAS,LKSUB,IFP,JKRES,IRO,IMISG,K INTEGER IPCOD,J,IRDP,IDECAY,IPH,M,L,JPCOD,JPAREN,III +ADD,CBHYPO,CFLOCV,0. *CMZ : 3.10/00 10/09/96 10.55.07 by Unknown +REP,CBHYPO,CFLOCV,39-41. INTEGER K,L,M,N REAL RMK C REAL CFT(6),CV(6),DP(6) DOUBLE PRECISION CH,PX,PY,PZ,EN,PP,PXY,KDIRC DOUBLE PRECISION MK,EK,PK,BET,GAM DOUBLE PRECISION COV0(3,3),JAC(3,3),COV1(3,3) +REP,CBHYPO,CFLOCV,96. CALL DLRMEM(CHVPAR,36) +REP,CBHYPO,CFLOCV,139-140. PXY=DSQRT(PX**2+PY**2) EN=Q(JTVVP+7) PP=DSQRT(PX**2+PY**2+PZ**2) C C CV(1)=Q(JTVVP+8) C CV(2)=Q(JTVVP+9) C CV(3)=Q(JTVVP+10) C CV(4)=Q(JTVVP+11) C CV(5)=Q(JTVVP+12) C CV(6)=Q(JTVVP+13) C IF(CFCHCV) THEN COV0(1,1)=Q(JTVVP+8) COV0(2,1)=Q(JTVVP+9) COV0(2,2)=Q(JTVVP+10) COV0(3,1)=Q(JTVVP+11) COV0(3,2)=Q(JTVVP+12) COV0(3,3)=Q(JTVVP+13) COV0(1,2)=COV0(2,1) COV0(1,3)=COV0(3,1) COV0(2,3)=COV0(3,2) C JAC(1,1)=-PY/PXY**2 JAC(1,2)=PX/PXY**2 JAC(1,3)=0.D0 JAC(2,1)=-PX/PXY**3 JAC(2,2)=-PY/PXY**3 JAC(2,3)=0.D0 JAC(3,1)=-PX*PZ/PXY**3 JAC(3,2)=-PY*PZ/PXY**3 JAC(3,3)=1.D0/PXY C DO M=1,3 DO N=1,M COV1(M,N)=0.D0 DO K=1,3 DO L=1,3 COV1(M,N)=COV1(M,N)+JAC(M,K)*JAC(N,L)*COV0(K,L) ENDDO ENDDO ENDDO ENDDO COV1(1,2)=COV1(2,1) COV1(1,3)=COV1(3,1) COV1(2,3)=COV1(3,2) ENDIF +REP,CBHYPO,CFLOCV,156-160. C CHVPAR(1,I)=ATAN2(PY,PX) CHVPAR(2,I)=1.D0/PXY CHVPAR(3,I)=PZ/PXY IF(CFCHCV) THEN CHVPAR(4,I)=COV1(1,1) CHVPAR(5,I)=COV1(2,1) CHVPAR(6,I)=COV1(2,2) CHVPAR(7,I)=COV1(3,1) CHVPAR(8,I)=COV1(3,2) CHVPAR(9,I)=COV1(3,3) ENDIF +REP,CBHYPO,CFLOCV,174-176. C...test of covariances, to be removed from here C DP(1)=Q(JTTKS+21)**2 C DP(2)=0. C DP(3)=Q(JTTKS+22)**2 C DP(4)=0. C DP(5)=0. C DP(6)=Q(JTTKS+23)**2 C DO K=1,6 C CFT(K)=Q(JTTKS+45+K)/CHVPAR(3+K,I) C CALL HFILL(50+K,CFT(K),0.,1.) C CALL HFILL(60+K,CV(K),0.,1.) C CALL HFILL(70+K,DP(K),0.,1.) C ENDDO C...to here Q(JTTKS+40)=CHVPAR(1,I) Q(JTTKS+41)=CHVPAR(2,I) Q(JTTKS+42)=CHVPAR(3,I) IF(CFCHCV) THEN Q(JTTKS+46)=CHVPAR(4,I) Q(JTTKS+47)=CHVPAR(5,I) Q(JTTKS+48)=CHVPAR(6,I) Q(JTTKS+49)=CHVPAR(7,I) Q(JTTKS+50)=CHVPAR(8,I) Q(JTTKS+51)=CHVPAR(9,I) ENDIF +REP,CBHYPO,CFLOCV,187-188. WRITE(LDBG,'(''CHVPAR :''/4G13.5)') & ((CHVPAR(III,JJJ),JJJ=1,4),III=1,9) +ADD,$VERSION,V3_10,*. &DECK,V3_10. *CMZ : 3.10/00 10/09/96 12.20.55 by Pal Hidas *-- Author : Pal Hidas 10/09/96 * *::> VERSION 3.10/00 10/09/96 12.20.55 * * Customized for PC (except for Hollerith constants of Zebra). * Unused variables removed. * * CFCHCV logical option introduced. If it is .TRUE. then the * covariance matrix of TVVP is used instead of TTKS, i.e. the * result of Locater's charged vertex fit is taken into account. * The covariance matrix of (phi,1/pxy,tgl) is calculated from * that of (px,py,pz) by Jacobi-transformation (error propagation). *