Warning, /include/hbook/hcqcom.inc is written in an unsupported language. File is not indexed.
0001 *
0002 * $Id$
0003 *
0004 * $Log$
0005 * Revision 1.1 1996/01/16 17:07:52 mclareni
0006 * Initial revision
0007 *
0008 *
0009 #ifndef CERNLIB_HBOOK_HCQCOM_INC
0010 #define CERNLIB_HBOOK_HCQCOM_INC
0011 *
0012 *
0013 * hcqcom.inc
0014 *
0015 *CMZ : 4.19/00 13/04/93 16.36.40 by Rene Brun
0016 *-- Author :
0017 * Basic parameters and data for multiquadric fit.
0018 * VPRINT is true if verbose option selected.
0019 LOGICAL VPRINT
0020 * NTUPLE is .TRUE. if ntuple is being processed, otherwise .FALSE.
0021 LOGICAL NTUPLE
0022 * IMQFUN function flag: = 0, no function evaluated yet.
0023 * = 1, linear sum of multiquadrics.
0024 * = 2, exponential of linear sum of multiquadrics.
0025 * (Note: HQF... routines take account of this, HQD... always calculate the
0026 * simple sum of multiquadrics.)
0027 INTEGER IMQFUN
0028 * IDMQ stores current histogram no. when required.
0029 INTEGER IDMQ
0030 * JMQ, KMQ store current multiquadric basis function nos. when required.
0031 INTEGER JMQ, KMQ
0032 * NSIG is number of significant points.
0033 INTEGER NSIG
0034 * NFREE is number of free parameters in the fitting (often equal to NSIG).
0035 INTEGER NFREE
0036 * Limits and recommended values of no. of significant points.
0037 INTEGER NSMIN, NSMIN1, NSMIN2, NSMIN3, NSREC1, NSREC2, NSREC3
0038 PARAMETER (NSMIN1 = 5, NSMIN2 = 10, NSMIN3 = 10)
0039 PARAMETER (NSREC1 = 15, NSREC2 = 50, NSREC3 = 50)
0040 * NSMAX its maximum no. of significant points.
0041 INTEGER NSMAX
0042 PARAMETER (NSMAX = 500)
0043 * NDIM is the number of dimensions - NDMAX is its maximum.
0044 INTEGER NDIM, NDMAX
0045 PARAMETER (NDMAX = 3)
0046 * NDIMC is dimensionality of data currently under consideration (<= NDIM).
0047 * NDIMWH is which of the NDIM1 dimensions are currently being considered.
0048 INTEGER NDIMC, NDIMWH (NDMAX)
0049 * NBINS is no. of bins in each coordinate, NBTOT the total
0050 * no. of bins in current histogram.
0051 INTEGER NBINS (NDMAX)
0052 INTEGER NX, NY, NZ, NBTOT
0053 EQUIVALENCE (NX, NBINS (1)), (NY, NBINS (2)), (NZ, NBINS (3))
0054 * Similarly, NBINSN is no. of bins for which first significant points where
0055 * found when surveying Ntuple data.
0056 INTEGER NBINSN (NDMAX)
0057 INTEGER NXN, NYN, NZN
0058 EQUIVALENCE (NXN, NBINSN (1)), (NYN, NBINSN (2)),
0059 +(NZN, NBINSN (3))
0060 * NNBINn is no. of bins used to bin ntuple data in n-dimensions.
0061 INTEGER NNBIN1, NNBIN2, NNBIN3
0062 PARAMETER (NNBIN1 = 100, NNBIN2 = 40, NNBIN3 = 40)
0063 * NMQEVS is no. of events (see also WMQEVS).
0064 INTEGER NMQEVS
0065 * NMQFUL is no. of histogram bins with non-zero content.
0066 INTEGER NMQFUL
0067 * NINTVS, NQUADS are no. of intervals and no. of gaussian quadrature points
0068 * per interval for integration of multiquadric functions.
0069 INTEGER NINTVS, NQUADS
0070 PARAMETER (NINTVS = 50, NQUADS = 8)
0071 * LAREA is working space for ZEBRA link maintenance.
0072 * LHQUAD is pointer to parent bank in Q in /PAWC/ in SEQ HCBOOK.
0073 * (For bank structure see HQINIT header.)
0074 * LnH and LnV are pointers to current n-D histogram contents and variances.
0075 * LHQCOV is pointer to covariances.
0076 * LHQDER is working bank for double derivatives.
0077 * LHQWK1/2.. are more working banks.
0078 * LHQDJN is pointer to temporary bank holding function values.
0079 * LLAST is a dummy pointer - always last!
0080 INTEGER LAREA (2), LHQUAD, L1H, L1V, L2H, L2V, L3H, L3V,
0081 +LHQCOV, LHQDER, LHQWK1, LHQWK2, LHQWK3,
0082 +LHQDJN, LLAST
0083 *
0084 * SIGGRD contains the gradient of the function with respect to the strengths.
0085 * SIGDEN the (approximate) density of the function at the significant point.
0086 REAL SIGGRD (NSMAX)
0087 REAL SIGDEN (NSMAX)
0088 * SIGV contains significant points' coordinate vectors.
0089 REAL SIGV (NSMAX, NDMAX)
0090 REAL SIGX (NSMAX), SIGY (NSMAX), SIGZ (NSMAX)
0091 EQUIVALENCE (SIGX, SIGV), (SIGY, SIGV (1, 2)), (SIGZ, SIGV (1, 3))
0092 * SIGDEL contains the scale parameter Delta for each coordinate.
0093 REAL SIGDEL (NSMAX)
0094 * SIGVMI, SIGVMA, SIGVBI, SIGVT are lower and upper limits, bin size
0095 * and total range.
0096 REAL SIGVMI (NDMAX), SIGVMA (NDMAX), SIGVBI (NDMAX), SIGVT (NDMAX)
0097 REAL XMI, XMA, DX, DXT, YMI, YMA, DY, DYT, ZMI, ZMA, DZ, DZT
0098 EQUIVALENCE (XMI, SIGVMI (1)), (XMA, SIGVMA (1))
0099 EQUIVALENCE (DX, SIGVBI (1)), (DXT, SIGVT (1))
0100 EQUIVALENCE (YMI, SIGVMI (2)), (YMA, SIGVMA (2))
0101 EQUIVALENCE (DY, SIGVBI (2)), (DYT, SIGVT (2))
0102 EQUIVALENCE (ZMI, SIGVMI (3)), (ZMA, SIGVMA (3))
0103 EQUIVALENCE (DZ, SIGVBI (3)), (DZT, SIGVT (3))
0104 * BINV is used to store current coordinates - usually bin centres, but can
0105 * be edges 0. or 1.
0106 REAL BINV (NDMAX)
0107 REAL XBIN, YBIN, ZBIN
0108 EQUIVALENCE (XBIN, BINV (1)), (YBIN, BINV (2)), (ZBIN, BINV (3))
0109 * VOLTOT is total "volume" and VOLBIN is "volume" of bin in current histogram.
0110 REAL VOLTOT, VOLBIN
0111 * WMQEVS is weighted no. of events (see also NMQEVS).
0112 REAL WMQEVS
0113 * VSCALE scales variances (only if incorrectly assigned errors suspected).
0114 REAL VSCALE
0115 * VMEAN1, VMEAN2 are mean variance per histogram entry and its square.
0116 REAL VMEAN1, VMEAN2
0117 * VCONST is a constant multiplying exponential of sum of multiquadrics for
0118 * Poisson likelihood fitting.
0119 REAL VCONST
0120 * SENSIT is a sensitivity parameter - normally 1 - which controls the
0121 * selection of significant points.
0122 REAL SENSIT
0123 * SLCUT is current value of cut for Laplacian significance.
0124 REAL SLCUT
0125 * SLRMS is rms of Laplacian significance.
0126 REAL SLRMS
0127 * SPREAD is a smoothness parameter - normally 1 - which multiplies the
0128 * scale parameter (radius of curvature).
0129 REAL SPREAD
0130 * SMALLD is a "small" value of the multiquadric parameter DELTA for edges, etc.
0131 * (usually set to 1.E-6 in HQUAD, HQUADN).
0132 REAL SMALLD
0133 * HQMIN, HQMAX, HQFMIN, HQFMAX are current minimum and maximum bin heights and
0134 * function values (at bin centres).
0135 REAL HQMIN, HQMAX, HQFMIN, HQFMAX
0136 *
0137 COMMON /HCQCOM/ SIGGRD, SIGDEN, SIGV, SIGDEL,
0138 + SIGVMI, SIGVMA, SIGVBI, SIGVT,
0139 + BINV, VOLTOT, VOLBIN, WMQEVS, VSCALE, VMEAN1, VMEAN2,
0140 + VCONST, SENSIT, SLCUT, SLRMS, SPREAD, SMALLD,
0141 + HQMIN, HQMAX, HQFMIN, HQFMAX,
0142 + IMQFUN, IDMQ, JMQ, KMQ,
0143 + NSIG, NFREE, NSMIN, NDIM, NDIMC, NDIMWH,
0144 + NBINS, NBTOT, NBINSN, NMQEVS, NMQFUL,
0145 + LAREA, LHQUAD, L1H, L1V, L2H, L2V, L3H, L3V,
0146 + LHQCOV, LHQDER, LHQWK1, LHQWK2, LHQWK3,
0147 + LHQDJN, LLAST,
0148 + VPRINT, NTUPLE
0149
0150
0151 #endif
