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