inc/F77/common_mini_2.f

***************************************************************************
*
*     Common for the track analysis, performed by the routine mini.f
*
***************************************************************************
c      PARAMETER (PIGR=3.141592654)!???????????????????

*     -------------------------------------------------------------------------
*     A TRACK is defined by: 
*     - the measured points associated to it 
*     - their spatial resolution 
*     - the track status vector AL-pha 
*     
*     This common must be filled when calling mini_2. After the fit, 
*     the initial values of AL-pha (initial guess, for instance from Hough) 
*     are replaced with the minimization result.
*
*     About the measured points:
*     - a COUPLE defines a point in the space; in this case you should set
*
*       XGOOD     = 1
*       YGOOD     = 1
*       XM,YM,ZM  = coordinates of the point in cm
*       RESX      = spatial resolution of cl-x
*       RESY      = spatial resolution of cl-y
*
*     - a SINGLET (single cluster x/y) defines a segment AB, parallel to 
*       the x/y strips; in this case you should set, for example in the case
*       of a cl-x:
*
*       XGOOD     = 1
*       YGOOD     = 0
*       ZM        = average z coordinate of the segment AB in cm
*       XM_A,YM_A = coordinates of extreme point A
*       XM_B,YM_B = coordinates of extreme point B
*       RESX      = spatial resolution of cl-x
*
*     NB!!! ZM must be always initialized, at least with the mechanical positions
*
*     (All the needed quantities for couples and singlets can be evaluated
*     by means of the routine xyz_PAM)
*     -------------------------------------------------------------------------
      REAL*8 AL
      REAL*8 XM,YM,ZM
      REAL*8 XM_A,YM_A!,ZM_A
      REAL*8 XM_B,YM_B!,ZM_B
      REAL*8 XGOOD,YGOOD
      REAL*8 RESX,RESY
c     REAL*8 DEDXTRK            !(1)
      REAL*8 DEDXTRK_X          !(1)
      REAL*8 DEDXTRK_Y          !(1)
      INTEGER IDCAND            !(2)
      COMMON/TRACK/
     $      AL(5)                                      !track parameters
     $     ,XM(NPLANES),YM(NPLANES),ZM(NPLANES)        !couple    => xyz point 
     $     ,XM_A(NPLANES),YM_A(NPLANES)!,ZM_A(NPLANES) !single cl => segment AB
     $     ,XM_B(NPLANES),YM_B(NPLANES)!,ZM_B(NPLANES) !
     $     ,RESX(NPLANES),RESY(NPLANES)                !spatial resolution
     $     ,XGOOD(NPLANES),YGOOD(NPLANES)              !mask of included views
c     $     ,DEDXTRK(NPLANES)                           !energy release
     $     ,DEDXTRK_X(NPLANES),DEDXTRK_Y(NPLANES) !(1)
     $     ,IDCAND              !(2) track-candidate ID (to relate some level1 info)
c******************************************************
cccccc 06/10/2005 modified by elena vannuccini ---> (1)
cccccc 17/08/2006 modified by elena vannuccini ---> (2)
c******************************************************

*     -------------------------------------------------------------------------
*     The FIT info includes:
*     - REDUCED chi2 
*     - number of d.o.f.
*     - evaluated coordinates (track intersection with the silicon planes)
*     - projected angles
*     - parameter covariance matrix
*     -------------------------------------------------------------------------
      REAL*8 CHI2
      REAL*8 XV,YV,ZV
      REAL*8 AXV,AYV
      REAL*8 COV
      COMMON/FIT_INFO/
     $      CHI2
     $     ,XV(NPLANES),YV(NPLANES),ZV(NPLANES)
     $     ,AXV(NPLANES),AYV(NPLANES)
     $     ,COV(5,5)            


*     -------------------------------------------------------------------------
*     Other commons needed by the minimization procedures:
*     -------------------------------------------------------------------------
*     computation of chi^2 and its derivatives 
      REAL*8 CHI2D,CHI2DD       !first and second derivatives of chi^2 
      REAL*8 DXDAL,DYDAL        !first derivatives of x and y on AL-pha 
      REAL*8 STEPAL             !
      REAL*8 ZINI               !z coordinate of reference plane
      COMMON/MINICOM/
     $      CHI2D(5),CHI2DD(5,5)
     $     ,DXDAL(NPLANES,5),DYDAL(NPLANES,5)
     $     ,STEPAL(5)
     $     ,ZINI
*     -------------------------------------------------------------------------
*     convergence conditions 
      REAL*8 ALTOL              !tolerances on AL-pha parameters
      REAL*8 ALMAX,ALMIN        !maximum and minimum value for AL-pha parameters
      INTEGER ISTEPMAX          !maximum number of steps      
      COMMON/MINICON/
     $      ALTOL(5)
     $     ,ALMAX(5),ALMIN(5)
     $     ,ISTEPMAX
*     -------------------------------------------------------------------------
*     variables for the tracking (routine posxy)
      REAL*8 VECT,VECTINI,VOUT  !vectors for Runge-Kutta method
      REAL*8 TOLL               !tolerance from plane
      REAL*8 STEPMAX            !maximum number of steps
      REAL*8 CHARGE             !sign of electric charge
      COMMON/TRACKING/VECT(7),VECTINI(7),VOUT(7),TOLL,STEPMAX,
     $     CHARGE

     
1	***************************************************************************
2	*
3	* Common for the track analysis, performed by the routine mini.f
4	*
5	***************************************************************************
6	c PARAMETER (PIGR=3.141592654)!???????????????????
7
8	* -------------------------------------------------------------------------
9	* A TRACK is defined by:
10	* - the measured points associated to it
11	* - their spatial resolution
12	* - the track status vector AL-pha
13	*
14	* This common must be filled when calling mini_2. After the fit,
15	* the initial values of AL-pha (initial guess, for instance from Hough)
16	* are replaced with the minimization result.
17	*
18	* About the measured points:
19	* - a COUPLE defines a point in the space; in this case you should set
20	*
21	* XGOOD = 1
22	* YGOOD = 1
23	* XM,YM,ZM = coordinates of the point in cm
24	* RESX = spatial resolution of cl-x
25	* RESY = spatial resolution of cl-y
26	*
27	* - a SINGLET (single cluster x/y) defines a segment AB, parallel to
28	* the x/y strips; in this case you should set, for example in the case
29	* of a cl-x:
30	*
31	* XGOOD = 1
32	* YGOOD = 0
33	* ZM = average z coordinate of the segment AB in cm
34	* XM_A,YM_A = coordinates of extreme point A
35	* XM_B,YM_B = coordinates of extreme point B
36	* RESX = spatial resolution of cl-x
37	*
38	* NB!!! ZM must be always initialized, at least with the mechanical positions
39	*
40	* (All the needed quantities for couples and singlets can be evaluated
41	* by means of the routine xyz_PAM)
42	* -------------------------------------------------------------------------
43	REAL*8 AL
44	REAL*8 XM,YM,ZM
45	REAL*8 XM_A,YM_A!,ZM_A
46	REAL*8 XM_B,YM_B!,ZM_B
47	REAL*8 XGOOD,YGOOD
48	REAL*8 RESX,RESY
49	c REAL*8 DEDXTRK !(1)
50	REAL*8 DEDXTRK_X !(1)
51	REAL*8 DEDXTRK_Y !(1)
52	INTEGER IDCAND !(2)
53	COMMON/TRACK/
54	$ AL(5) !track parameters
55	$ ,XM(NPLANES),YM(NPLANES),ZM(NPLANES) !couple => xyz point
56	$ ,XM_A(NPLANES),YM_A(NPLANES)!,ZM_A(NPLANES) !single cl => segment AB
57	$ ,XM_B(NPLANES),YM_B(NPLANES)!,ZM_B(NPLANES) !
58	$ ,RESX(NPLANES),RESY(NPLANES) !spatial resolution
59	$ ,XGOOD(NPLANES),YGOOD(NPLANES) !mask of included views
60	c $ ,DEDXTRK(NPLANES) !energy release
61	$ ,DEDXTRK_X(NPLANES),DEDXTRK_Y(NPLANES) !(1)
62	$ ,IDCAND !(2) track-candidate ID (to relate some level1 info)
63	c******************************************************
64	cccccc 06/10/2005 modified by elena vannuccini ---> (1)
65	cccccc 17/08/2006 modified by elena vannuccini ---> (2)
66	c******************************************************
67
68	* -------------------------------------------------------------------------
69	* The FIT info includes:
70	* - REDUCED chi2
71	* - number of d.o.f.
72	* - evaluated coordinates (track intersection with the silicon planes)
73	* - projected angles
74	* - parameter covariance matrix
75	* -------------------------------------------------------------------------
76	REAL*8 CHI2
77	REAL*8 XV,YV,ZV
78	REAL*8 AXV,AYV
79	REAL*8 COV
80	COMMON/FIT_INFO/
81	$ CHI2
82	$ ,XV(NPLANES),YV(NPLANES),ZV(NPLANES)
83	$ ,AXV(NPLANES),AYV(NPLANES)
84	$ ,COV(5,5)
85
86
87
88
89
90	* -------------------------------------------------------------------------
91	* Other commons needed by the minimization procedures:
92	* -------------------------------------------------------------------------
93	* computation of chi^2 and its derivatives
94	REAL*8 CHI2D,CHI2DD !first and second derivatives of chi^2
95	REAL*8 DXDAL,DYDAL !first derivatives of x and y on AL-pha
96	REAL*8 STEPAL !
97	REAL*8 ZINI !z coordinate of reference plane
98	COMMON/MINICOM/
99	$ CHI2D(5),CHI2DD(5,5)
100	$ ,DXDAL(NPLANES,5),DYDAL(NPLANES,5)
101	$ ,STEPAL(5)
102	$ ,ZINI
103	* -------------------------------------------------------------------------
104	* convergence conditions
105	REAL*8 ALTOL !tolerances on AL-pha parameters
106	REAL*8 ALMAX,ALMIN !maximum and minimum value for AL-pha parameters
107	INTEGER ISTEPMAX !maximum number of steps
108	COMMON/MINICON/
109	$ ALTOL(5)
110	$ ,ALMAX(5),ALMIN(5)
111	$ ,ISTEPMAX
112	* -------------------------------------------------------------------------
113	* variables for the tracking (routine posxy)
114	REAL*8 VECT,VECTINI,VOUT !vectors for Runge-Kutta method
115	REAL*8 TOLL !tolerance from plane
116	REAL*8 STEPMAX !maximum number of steps
117	REAL*8 CHARGE !sign of electric charge
118	COMMON/TRACKING/VECT(7),VECTINI(7),VOUT(7),TOLL,STEPMAX,
119	$ CHARGE
120
121