Ctvmft/interface/dctvmft.h

C========================= Include ctvmft.inc ==========================
C============ include file for the constrained fit module CTVMFT =======

C ----------------------------------------------------------------------
C Single precision stuff

C Counters of the used numbers of Vtx ..

      INTEGER     NVERTX                ! Number of vertices
      INTEGER     NMASSC                ! Number of mass constraints
      INTEGER     NTRACK                ! Number of tracks, this fit

C Vertex track associations

      LOGICAL     TRKVTX(MAXTRK,MAXVTX) ! (Track,Vertex) assoc. table
      INTEGER     VTXPNT(MAXVTX,2)      ! Vertex association information

C Conversion constraints info

      INTEGER     CVTX  (MAXVTX)        ! Conversion constraint specs
      LOGICAL     TRKMCN(MAXTRK,MAXMCN) ! (Track,Mass_Constraint) table
      REAL        CMASS (MAXMCN)        ! Constraint mass

C  Required Input Track Specifications

      CHARACTER*4 TKBANK(MAXTRK)        ! Track bank type (e.g. 'TRKS')
      INTEGER     LIST  (MAXTRK)        ! Track bank number
      REAL        TMASS (MAXTRK)        ! Track mass assignment

C  Fit input, collected using the input specifications

      REAL        PAR0  (5,MAXTRK)      ! Original trk parameter values
      REAL        G     (5,5,MAXTRK)    ! Fitted trk parameter covar.

C  Optional Input, Required Data for Pointing Constraints

      REAL        XYZPV0(3)             ! Primary vertex: XYZ
      REAL        EXYZPV(3,3)           ! Primary vertex: covariance
      REAL        XZSLOPE,YZSLOPE       ! beamline slopes vs z direction
      INTEGER     RUNNUM,TRGNUM         ! Run and Event Number
      INTEGER     ITER                  ! Number of iterative steps
      INTEGER     NTSCUT                ! Total number of cut steps

C  Output; Results of the Fit

      INTEGER     MATDIM                ! Dimension of matrix, this fit
      INTEGER     NDOF                  ! Number of degrees of freedom
      LOGICAL     VTXVTX(MAXVTX,MAXVTX) ! Vertex geneology
      REAL        CHISQR(0:MAXITR)      ! Overall fit Chi Square result
      REAL        CHIT  (MAXTRK)        ! Track contr. to Chi Square
      REAL        CHIV  (0:MAXVTX)      ! Vertex fit Chi Squares
      REAL        CHIM  (MAXMCN)        ! Mass constraint Chi Squares
      REAL        PAR   (5,MAXTRK)      ! Fitted trk helix parameters
      REAL        PARDIF(5,MAXTRK)      ! Parameter diff (fit-input)
      REAL        FMCDIF(MAXMCN)        ! Mass constraint residuals
      REAL        PCON  (MAXVTX,2)      ! Pointing constraint tests
      REAL        SANG  (MAXVTX,2)      ! Pointing constraint tests
      REAL        XYZVRT(3,0:MAXVTX)    ! Primary, Second vertices
      REAL        TRKP4 (MAXTRK,6)      ! Track Px,Py,Pz,E,Pt P
      REAL        VTXP4 (4,MAXVTX)      ! Vertex 4-momentum sum
      REAL        MCNP4 (4,MAXMCN)      ! 4-mom sum for mass constr.
      REAL        DDA   (MAXTRK,8)      ! d(Px,Py)/d(Crv,Phi,Xs,Ys)
      REAL        DXYZPV(3)             ! Primary displacement in fit
      INTEGER     VOFF  (MAXVTX)        ! Vertex      offset pointers
      INTEGER     POFF  (MAXVTX)        ! Point Phi   offset pointers
      INTEGER     COFF  (MAXVTX)        ! Point Cotan offset pointers
      INTEGER     TOFF  (MAXTRK)        ! Track       offset pointers
      INTEGER     MOFF                  ! Mass        offset pointers

C Memory

      INTEGER     TKERR (MAXTRK)        ! Trk error flag, previous fit
      INTEGER     IJKERR(3)             ! Error code reporting
      REAL        XYZPV (3)
      INTEGER     UVWXYZ(UDIM)
      EQUIVALENCE(XYZPV,XYZVRT(1,0))    ! Danger: equiv. statements
      EQUIVALENCE(UVWXYZ,ITER)

C  CTC first approximation and track-vertex step parameters

      REAL    DRMAX,RVMAX,DZMAX,TRNMAX,DSMIN

C Useful conversion parameter

      REAL    PSCALE                    ! Mom.scale conversion
                                        ! [cm^-1 -> Gev/c]

      COMMON /DCTVMFR/
     >        RUNNUM,TRGNUM,ITER,NTSCUT
     >       ,NVERTX,NMASSC,NTRACK
     >       ,TRKVTX,TRKMCN,VTXPNT,CMASS,CVTX
     >       ,VTXVTX
     >       ,TKBANK,LIST,TMASS
     >       ,MATDIM
     >       ,TKERR
     >       ,NDOF,CHISQR,CHIT,CHIV,CHIM
     >       ,XYZPV0,EXYZPV
     >       ,XZSLOPE,YZSLOPE                       
     >       ,XYZVRT,DXYZPV
     >       ,PAR,G
     >       ,TRKP4,VTXP4,MCNP4,DDA
     >       ,VOFF,TOFF,POFF,COFF,MOFF
     >       ,PAR0,PARDIF
     >       ,FMCDIF,PCON,SANG
     >       ,DRMAX,RVMAX,DZMAX,TRNMAX,DSMIN
     >       ,IJKERR
     >       ,PSCALE

C ----------------------------------------------------------------------
C Double precision stuff

C Covariance matrix, fit parameters
      REAL*8      VMAT(MAXDIM,MAXDIM+1)
      COMMON /DCTVMFD/
     >        VMAT
Revision:	1.1
Committed:	Wed Sep 17 04:01:48 2008 UTC (16 years, 7 months ago) by loizides
Content type:	text/plain
Branch:	MAIN
CVS Tags:	Mit_032, Mit_031, Mit_025c_branch2, Mit_025c_branch1, Mit_030, Mit_029c, Mit_030_pre1, Mit_029a, Mit_029, Mit_029_pre1, Mit_028a, Mit_025c_branch0, Mit_028, Mit_027a, Mit_027, Mit_026, Mit_025e, Mit_025d, Mit_025c, Mit_025b, Mit_025a, Mit_025, Mit_025pre2, Mit_024b, Mit_025pre1, Mit_024a, Mit_024, Mit_023, Mit_022a, Mit_022, Mit_020d, TMit_020d, Mit_020c, Mit_021, Mit_021pre2, Mit_021pre1, Mit_020b, Mit_020a, Mit_020, Mit_020pre1, Mit_018, Mit_017, Mit_017pre3, Mit_017pre2, Mit_017pre1, V07-05-00, Mit_016, Mit_015b, Mit_015a, Mit_015, Mit_014e, Mit_014d, Mit_014c, Mit_014b, ConvRejection-10-06-09, Mit_014a, Mit_014, Mit_014pre3, Mit_014pre2, Mit_014pre1, Mit_013d, Mit_013c, Mit_013b, Mit_013a, Mit_013, Mit_013pre1, Mit_012i, Mit_012g, Mit_012f, Mit_012e, Mit_012d, Mit_012c, Mit_012b, Mit_012a, Mit_012, Mit_011a, Mit_011, Mit_010a, Mit_010, Mit_009c, Mit_009b, Mit_009a, Mit_009, Mit_008, Mit_008pre2, Mit_008pre1, Mit_006b, Mit_006a, Mit_006, Mit_005, Mit_004, HEAD
Branch point for:	Mit_025c_branch
Log Message:	Moved MitVertex contents to MitCommon. MitVertex therefore is obsolute and should not be touched anymore.
#	Content
1	C========================= Include ctvmft.inc ==========================
2	C============ include file for the constrained fit module CTVMFT =======
3
4	C ----------------------------------------------------------------------
5	C Single precision stuff
6
7	C Counters of the used numbers of Vtx ..
8
9	INTEGER NVERTX ! Number of vertices
10	INTEGER NMASSC ! Number of mass constraints
11	INTEGER NTRACK ! Number of tracks, this fit
12
13	C Vertex track associations
14
15	LOGICAL TRKVTX(MAXTRK,MAXVTX) ! (Track,Vertex) assoc. table
16	INTEGER VTXPNT(MAXVTX,2) ! Vertex association information
17
18	C Conversion constraints info
19
20	INTEGER CVTX (MAXVTX) ! Conversion constraint specs
21	LOGICAL TRKMCN(MAXTRK,MAXMCN) ! (Track,Mass_Constraint) table
22	REAL CMASS (MAXMCN) ! Constraint mass
23
24	C Required Input Track Specifications
25
26	CHARACTER*4 TKBANK(MAXTRK) ! Track bank type (e.g. 'TRKS')
27	INTEGER LIST (MAXTRK) ! Track bank number
28	REAL TMASS (MAXTRK) ! Track mass assignment
29
30	C Fit input, collected using the input specifications
31
32	REAL PAR0 (5,MAXTRK) ! Original trk parameter values
33	REAL G (5,5,MAXTRK) ! Fitted trk parameter covar.
34
35	C Optional Input, Required Data for Pointing Constraints
36
37	REAL XYZPV0(3) ! Primary vertex: XYZ
38	REAL EXYZPV(3,3) ! Primary vertex: covariance
39	REAL XZSLOPE,YZSLOPE ! beamline slopes vs z direction
40	INTEGER RUNNUM,TRGNUM ! Run and Event Number
41	INTEGER ITER ! Number of iterative steps
42	INTEGER NTSCUT ! Total number of cut steps
43
44	C Output; Results of the Fit
45
46	INTEGER MATDIM ! Dimension of matrix, this fit
47	INTEGER NDOF ! Number of degrees of freedom
48	LOGICAL VTXVTX(MAXVTX,MAXVTX) ! Vertex geneology
49	REAL CHISQR(0:MAXITR) ! Overall fit Chi Square result
50	REAL CHIT (MAXTRK) ! Track contr. to Chi Square
51	REAL CHIV (0:MAXVTX) ! Vertex fit Chi Squares
52	REAL CHIM (MAXMCN) ! Mass constraint Chi Squares
53	REAL PAR (5,MAXTRK) ! Fitted trk helix parameters
54	REAL PARDIF(5,MAXTRK) ! Parameter diff (fit-input)
55	REAL FMCDIF(MAXMCN) ! Mass constraint residuals
56	REAL PCON (MAXVTX,2) ! Pointing constraint tests
57	REAL SANG (MAXVTX,2) ! Pointing constraint tests
58	REAL XYZVRT(3,0:MAXVTX) ! Primary, Second vertices
59	REAL TRKP4 (MAXTRK,6) ! Track Px,Py,Pz,E,Pt P
60	REAL VTXP4 (4,MAXVTX) ! Vertex 4-momentum sum
61	REAL MCNP4 (4,MAXMCN) ! 4-mom sum for mass constr.
62	REAL DDA (MAXTRK,8) ! d(Px,Py)/d(Crv,Phi,Xs,Ys)
63	REAL DXYZPV(3) ! Primary displacement in fit
64	INTEGER VOFF (MAXVTX) ! Vertex offset pointers
65	INTEGER POFF (MAXVTX) ! Point Phi offset pointers
66	INTEGER COFF (MAXVTX) ! Point Cotan offset pointers
67	INTEGER TOFF (MAXTRK) ! Track offset pointers
68	INTEGER MOFF ! Mass offset pointers
69
70	C Memory
71
72	INTEGER TKERR (MAXTRK) ! Trk error flag, previous fit
73	INTEGER IJKERR(3) ! Error code reporting
74	REAL XYZPV (3)
75	INTEGER UVWXYZ(UDIM)
76	EQUIVALENCE(XYZPV,XYZVRT(1,0)) ! Danger: equiv. statements
77	EQUIVALENCE(UVWXYZ,ITER)
78
79	C CTC first approximation and track-vertex step parameters
80
81	REAL DRMAX,RVMAX,DZMAX,TRNMAX,DSMIN
82
83	C Useful conversion parameter
84
85	REAL PSCALE ! Mom.scale conversion
86	! [cm^-1 -> Gev/c]
87
88	COMMON /DCTVMFR/
89	> RUNNUM,TRGNUM,ITER,NTSCUT
90	> ,NVERTX,NMASSC,NTRACK
91	> ,TRKVTX,TRKMCN,VTXPNT,CMASS,CVTX
92	> ,VTXVTX
93	> ,TKBANK,LIST,TMASS
94	> ,MATDIM
95	> ,TKERR
96	> ,NDOF,CHISQR,CHIT,CHIV,CHIM
97	> ,XYZPV0,EXYZPV
98	> ,XZSLOPE,YZSLOPE
99	> ,XYZVRT,DXYZPV
100	> ,PAR,G
101	> ,TRKP4,VTXP4,MCNP4,DDA
102	> ,VOFF,TOFF,POFF,COFF,MOFF
103	> ,PAR0,PARDIF
104	> ,FMCDIF,PCON,SANG
105	> ,DRMAX,RVMAX,DZMAX,TRNMAX,DSMIN
106	> ,IJKERR
107	> ,PSCALE
108
109	C ----------------------------------------------------------------------
110	C Double precision stuff
111
112	C Covariance matrix, fit parameters
113	REAL*8 VMAT(MAXDIM,MAXDIM+1)
114	COMMON /DCTVMFD/
115	> VMAT