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 TAILX,TAILY
      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
     $     ,TAILX(NPLANES),TAILY(NPLANES)              !tails in the 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 DAL
      REAL*8 CHI2D,CHI2DOLD,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/
     $      DAL(5)
     $     ,CHI2D(5),CHI2DOLD(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

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