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 |
|
44 |
|
45 |
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 |
REAL*8 RESX,RESY |
50 |
REAL*8 TAILX,TAILY |
51 |
REAL*8 XGOOD,YGOOD |
52 |
REAL*8 DEDXTRK_X |
53 |
REAL*8 DEDXTRK_Y |
54 |
REAL*8 ZINI ! ***PP*** |
55 |
REAL*8 PFIXED ! ***PP*** |
56 |
REAL*8 CHI2 |
57 |
REAL*8 XV,YV,ZV |
58 |
REAL*8 AXV,AYV |
59 |
REAL*8 COV |
60 |
REAL*8 FACT |
61 |
INTEGER NSTEP |
62 |
INTEGER IDCAND |
63 |
INTEGER TRACKMODE |
64 |
INTEGER ISTEPMIN |
65 |
COMMON/TRACK/ |
66 |
$ AL(5) !track parameters |
67 |
$ ,XM(NPLANES),YM(NPLANES),ZM(NPLANES) !couple => xyz point |
68 |
$ ,XM_A(NPLANES),YM_A(NPLANES),ZM_A(NPLANES) !single cl => segment AB |
69 |
$ ,XM_B(NPLANES),YM_B(NPLANES),ZM_B(NPLANES) ! |
70 |
$ ,RESX(NPLANES),RESY(NPLANES) !spatial resolution |
71 |
$ ,TAILX(NPLANES),TAILY(NPLANES) !tails in the resolution |
72 |
$ ,XGOOD(NPLANES),YGOOD(NPLANES) !mask of included views |
73 |
c $ ,DEDXTRK(NPLANES) !energy release |
74 |
$ ,DEDXTRK_X(NPLANES),DEDXTRK_Y(NPLANES) |
75 |
$ ,ZINI |
76 |
$ ,PFIXED |
77 |
$ ,CHI2 |
78 |
$ ,XV(NPLANES),YV(NPLANES),ZV(NPLANES) |
79 |
$ ,AXV(NPLANES),AYV(NPLANES) |
80 |
$ ,COV(5,5) |
81 |
$ ,FACT ! factor scale for tracking precision |
82 |
$ ,NSTEP |
83 |
$ ,IDCAND !(2) track-candidate ID (to relate some level1 info) |
84 |
$ ,TRACKMODE ! tracking mode, 0=normal, 1=with bremsstrahlung |
85 |
$ ,ISTEPMIN ! minimum number of step during tracking |
86 |
c****************************************************** |
87 |
cccccc 06/10/2005 modified by elena vannuccini ---> (1) |
88 |
cccccc 17/08/2006 modified by elena vannuccini ---> (2) |
89 |
c****************************************************** |
90 |
|
91 |
|
92 |
* ------------------------------------------------------------------------- |
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 |
c$$$ REAL*8 ZINI ! ***PP*** |
101 |
c$$$ REAL*8 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 |
|
114 |
|
115 |
|
116 |
* ------------------------------------------------------------------------- |
117 |
* Other commons needed by the minimization procedures: |
118 |
* ------------------------------------------------------------------------- |
119 |
* 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 |
REAL*8 DXDAL,DYDAL !first derivatives of x and y on AL-pha |
123 |
REAL*8 STEPAL ! |
124 |
***PP*** REAL*8 ZINI !z coordinate of reference plane |
125 |
COMMON/MINICOM/ |
126 |
$ DAL(5) |
127 |
$ ,CHI2D(5),CHI2DOLD(5),CHI2DD(5,5) |
128 |
$ ,DXDAL(NPLANES,5),DYDAL(NPLANES,5) |
129 |
$ ,STEPAL(5) |
130 |
***PP*** $ ,ZINI |
131 |
|
132 |
* ------------------------------------------------------------------------- |
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 |
|
142 |
|
143 |
* ------------------------------------------------------------------------- |
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 |
* ------------------------------------------------------------------------- |
153 |
* variables for the DeltaB coorection |
154 |
REAL*8 DELTA0,DELTA1,DLT ! DELTA-DEF = DELTA0 + DELTA1 x Y |
155 |
COMMON/DELTAB/DELTA0,DELTA1,DLT |
156 |
|
157 |
|