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.
#	User	Rev	Content
1	loizides	1.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