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 RESX,RESY
      REAL*8 XGOOD,YGOOD
      REAL*8 DEDXTRK_X          
      REAL*8 DEDXTRK_Y         
      REAL*8 ZINI     ! ***PP***
      REAL*8 PFIXED   ! ***PP***
      REAL*8 CHI2
      REAL*8 XV,YV,ZV
      REAL*8 AXV,AYV
      REAL*8 COV
      REAL*8 FACT
      INTEGER NSTEP
      INTEGER IDCAND            
      INTEGER TRACKMODE
      INTEGER ISTEPMIN
      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) 
     $     ,ZINI
     $     ,PFIXED
     $     ,CHI2
     $     ,XV(NPLANES),YV(NPLANES),ZV(NPLANES)
     $     ,AXV(NPLANES),AYV(NPLANES)
     $     ,COV(5,5)
     $     ,FACT                ! factor scale for tracking precision
     $     ,NSTEP
     $     ,IDCAND              !(2) track-candidate ID (to relate some level1 info)
     $     ,TRACKMODE           ! tracking mode, 0=normal, 1=with bremsstrahlung 
     $     ,ISTEPMIN            ! minimum number of step during tracking
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
*     -------------------------------------------------------------------------
c$$$      REAL*8 ZINI     ! ***PP***
c$$$      REAL*8 PFIXED   ! ***PP***
c$$$      REAL*8 CHI2
c$$$      REAL*8 XV,YV,ZV
c$$$      REAL*8 AXV,AYV
c$$$      REAL*8 COV
c$$$      COMMON/FIT_INFO/
c$$$     $      ZINI
c$$$     $     ,PFIXED
c$$$     $     ,CHI2
c$$$     $     ,XV(NPLANES),YV(NPLANES),ZV(NPLANES)
c$$$     $     ,AXV(NPLANES),AYV(NPLANES)
c$$$     $     ,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             !
***PP***      REAL*8 ZINI               !z coordinate of reference plane
      COMMON/MINICOM/
     $      CHI2D(5),CHI2DD(5,5)
     $     ,DXDAL(NPLANES,5),DYDAL(NPLANES,5)
     $     ,STEPAL(5)
***PP***     $     ,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	mocchiut	1.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	pam-fi	1.3	REAL*8 RESX,RESY
48	mocchiut	1.1	REAL*8 XGOOD,YGOOD
49	pam-fi	1.6	REAL*8 DEDXTRK_X
50			REAL*8 DEDXTRK_Y
51	pam-fi	1.3	REAL8 ZINI ! PP*
52			REAL8 PFIXED ! PP*
53			REAL*8 CHI2
54			REAL*8 XV,YV,ZV
55			REAL*8 AXV,AYV
56			REAL*8 COV
57	pam-fi	1.5	REAL*8 FACT
58	pam-fi	1.3	INTEGER NSTEP
59	pam-fi	1.6	INTEGER IDCAND
60	pam-fi	1.4	INTEGER TRACKMODE
61	pam-fi	1.5	INTEGER ISTEPMIN
62	mocchiut	1.1	COMMON/TRACK/
63			$ AL(5) !track parameters
64			$ ,XM(NPLANES),YM(NPLANES),ZM(NPLANES) !couple => xyz point
65			$ ,XM_A(NPLANES),YM_A(NPLANES)!,ZM_A(NPLANES) !single cl => segment AB
66			$ ,XM_B(NPLANES),YM_B(NPLANES)!,ZM_B(NPLANES) !
67			$ ,RESX(NPLANES),RESY(NPLANES) !spatial resolution
68			$ ,XGOOD(NPLANES),YGOOD(NPLANES) !mask of included views
69			c $ ,DEDXTRK(NPLANES) !energy release
70	pam-fi	1.6	$ ,DEDXTRK_X(NPLANES),DEDXTRK_Y(NPLANES)
71	pam-fi	1.3	$ ,ZINI
72			$ ,PFIXED
73			$ ,CHI2
74			$ ,XV(NPLANES),YV(NPLANES),ZV(NPLANES)
75			$ ,AXV(NPLANES),AYV(NPLANES)
76			$ ,COV(5,5)
77	pam-fi	1.5	$ ,FACT ! factor scale for tracking precision
78	pam-fi	1.3	$ ,NSTEP
79	pam-fi	1.2	$ ,IDCAND !(2) track-candidate ID (to relate some level1 info)
80	pam-fi	1.4	$ ,TRACKMODE ! tracking mode, 0=normal, 1=with bremsstrahlung
81	pam-fi	1.5	$ ,ISTEPMIN ! minimum number of step during tracking
82	mocchiut	1.1	c******************************************************
83			cccccc 06/10/2005 modified by elena vannuccini ---> (1)
84	pam-fi	1.2	cccccc 17/08/2006 modified by elena vannuccini ---> (2)
85	mocchiut	1.1	c******************************************************
86
87			* -------------------------------------------------------------------------
88			* The FIT info includes:
89			* - REDUCED chi2
90			* - number of d.o.f.
91			* - evaluated coordinates (track intersection with the silicon planes)
92			* - projected angles
93			* - parameter covariance matrix
94			* -------------------------------------------------------------------------
95	pam-fi	1.3	c$$$ REAL8 ZINI ! PP*
96			c$$$ REAL8 PFIXED ! PP*
97			c$$$ REAL*8 CHI2
98			c$$$ REAL*8 XV,YV,ZV
99			c$$$ REAL*8 AXV,AYV
100			c$$$ REAL*8 COV
101			c$$$ COMMON/FIT_INFO/
102			c$$$ $ ZINI
103			c$$$ $ ,PFIXED
104			c$$$ $ ,CHI2
105			c$$$ $ ,XV(NPLANES),YV(NPLANES),ZV(NPLANES)
106			c$$$ $ ,AXV(NPLANES),AYV(NPLANES)
107			c$$$ $ ,COV(5,5)
108	mocchiut	1.1
109
110
111
112
113			* -------------------------------------------------------------------------
114			* Other commons needed by the minimization procedures:
115			* -------------------------------------------------------------------------
116			* computation of chi^2 and its derivatives
117			REAL*8 CHI2D,CHI2DD !first and second derivatives of chi^2
118			REAL*8 DXDAL,DYDAL !first derivatives of x and y on AL-pha
119			REAL*8 STEPAL !
120	pam-fi	1.3	*PP* REAL*8 ZINI !z coordinate of reference plane
121	mocchiut	1.1	COMMON/MINICOM/
122			$ CHI2D(5),CHI2DD(5,5)
123			$ ,DXDAL(NPLANES,5),DYDAL(NPLANES,5)
124			$ ,STEPAL(5)
125	pam-fi	1.3	*PP* $ ,ZINI
126	mocchiut	1.1	* -------------------------------------------------------------------------
127			* convergence conditions
128			REAL*8 ALTOL !tolerances on AL-pha parameters
129			REAL*8 ALMAX,ALMIN !maximum and minimum value for AL-pha parameters
130			INTEGER ISTEPMAX !maximum number of steps
131			COMMON/MINICON/
132			$ ALTOL(5)
133			$ ,ALMAX(5),ALMIN(5)
134			$ ,ISTEPMAX
135			* -------------------------------------------------------------------------
136			* variables for the tracking (routine posxy)
137			REAL*8 VECT,VECTINI,VOUT !vectors for Runge-Kutta method
138			REAL*8 TOLL !tolerance from plane
139			REAL*8 STEPMAX !maximum number of steps
140			REAL*8 CHARGE !sign of electric charge
141			COMMON/TRACKING/VECT(7),VECTINI(7),VOUT(7),TOLL,STEPMAX,
142			$ CHARGE
143
144