40 |
* (All the needed quantities for couples and singlets can be evaluated |
* (All the needed quantities for couples and singlets can be evaluated |
41 |
* by means of the routine xyz_PAM) |
* by means of the routine xyz_PAM) |
42 |
* ------------------------------------------------------------------------- |
* ------------------------------------------------------------------------- |
43 |
|
|
44 |
|
|
45 |
REAL*8 AL |
REAL*8 AL |
46 |
REAL*8 XM,YM,ZM |
REAL*8 XM,YM,ZM |
47 |
REAL*8 XM_A,YM_A!,ZM_A |
REAL*8 XM_A,YM_A,ZM_A |
48 |
REAL*8 XM_B,YM_B!,ZM_B |
REAL*8 XM_B,YM_B,ZM_B |
49 |
REAL*8 RESX,RESY |
REAL*8 RESX,RESY |
50 |
|
REAL*8 TAILX,TAILY |
51 |
REAL*8 XGOOD,YGOOD |
REAL*8 XGOOD,YGOOD |
52 |
c REAL*8 DEDXTRK !(1) |
REAL*8 DEDXTRK_X |
53 |
REAL*8 DEDXTRK_X !(1) |
REAL*8 DEDXTRK_Y |
|
REAL*8 DEDXTRK_Y !(1) |
|
54 |
REAL*8 ZINI ! ***PP*** |
REAL*8 ZINI ! ***PP*** |
55 |
REAL*8 PFIXED ! ***PP*** |
REAL*8 PFIXED ! ***PP*** |
56 |
REAL*8 CHI2 |
REAL*8 CHI2 |
57 |
REAL*8 XV,YV,ZV |
REAL*8 XV,YV,ZV |
58 |
|
REAL*8 XV_A,YV_A,ZV_A |
59 |
|
REAL*8 XV_B,YV_B,ZV_B |
60 |
REAL*8 AXV,AYV |
REAL*8 AXV,AYV |
61 |
REAL*8 COV |
REAL*8 COV |
62 |
|
REAL*8 FACT |
63 |
INTEGER NSTEP |
INTEGER NSTEP |
64 |
INTEGER IDCAND !(2) |
INTEGER IDCAND |
65 |
|
INTEGER TRACKMODE |
66 |
|
INTEGER ISTEPMIN |
67 |
COMMON/TRACK/ |
COMMON/TRACK/ |
68 |
$ AL(5) !track parameters |
$ AL(5) !track parameters |
69 |
$ ,XM(NPLANES),YM(NPLANES),ZM(NPLANES) !couple => xyz point |
$ ,XM(NPLANES),YM(NPLANES),ZM(NPLANES) !couple => xyz point |
70 |
$ ,XM_A(NPLANES),YM_A(NPLANES)!,ZM_A(NPLANES) !single cl => segment AB |
$ ,XM_A(NPLANES),YM_A(NPLANES),ZM_A(NPLANES) !single cl => segment AB |
71 |
$ ,XM_B(NPLANES),YM_B(NPLANES)!,ZM_B(NPLANES) ! |
$ ,XM_B(NPLANES),YM_B(NPLANES),ZM_B(NPLANES) ! |
72 |
$ ,RESX(NPLANES),RESY(NPLANES) !spatial resolution |
$ ,RESX(NPLANES),RESY(NPLANES) !spatial resolution |
73 |
|
$ ,TAILX(NPLANES),TAILY(NPLANES) !tails in the resolution |
74 |
$ ,XGOOD(NPLANES),YGOOD(NPLANES) !mask of included views |
$ ,XGOOD(NPLANES),YGOOD(NPLANES) !mask of included views |
75 |
c $ ,DEDXTRK(NPLANES) !energy release |
c $ ,DEDXTRK(NPLANES) !energy release |
76 |
$ ,DEDXTRK_X(NPLANES),DEDXTRK_Y(NPLANES) !(1) |
$ ,DEDXTRK_X(NPLANES),DEDXTRK_Y(NPLANES) |
77 |
$ ,ZINI |
$ ,ZINI |
78 |
$ ,PFIXED |
$ ,PFIXED |
79 |
$ ,CHI2 |
$ ,CHI2 |
80 |
$ ,XV(NPLANES),YV(NPLANES),ZV(NPLANES) |
$ ,XV(NPLANES),YV(NPLANES),ZV(NPLANES) |
81 |
|
$ ,XV_A(NPLANES),YV_A(NPLANES),ZV_A(NPLANES) !single cl => segment AB |
82 |
|
$ ,XV_B(NPLANES),YV_B(NPLANES),ZV_B(NPLANES) ! |
83 |
$ ,AXV(NPLANES),AYV(NPLANES) |
$ ,AXV(NPLANES),AYV(NPLANES) |
84 |
$ ,COV(5,5) |
$ ,COV(5,5) |
85 |
|
$ ,FACT ! factor scale for tracking precision |
86 |
$ ,NSTEP |
$ ,NSTEP |
87 |
$ ,IDCAND !(2) track-candidate ID (to relate some level1 info) |
$ ,IDCAND !(2) track-candidate ID (to relate some level1 info) |
88 |
|
$ ,TRACKMODE ! tracking mode, 0=normal, 1=with bremsstrahlung |
89 |
|
$ ,ISTEPMIN ! minimum number of step during tracking |
90 |
c****************************************************** |
c****************************************************** |
91 |
cccccc 06/10/2005 modified by elena vannuccini ---> (1) |
cccccc 06/10/2005 modified by elena vannuccini ---> (1) |
92 |
cccccc 17/08/2006 modified by elena vannuccini ---> (2) |
cccccc 17/08/2006 modified by elena vannuccini ---> (2) |
93 |
c****************************************************** |
c****************************************************** |
94 |
|
|
95 |
|
|
96 |
* ------------------------------------------------------------------------- |
* ------------------------------------------------------------------------- |
97 |
* The FIT info includes: |
* The FIT info includes: |
98 |
* - REDUCED chi2 |
* - REDUCED chi2 |
117 |
|
|
118 |
|
|
119 |
|
|
|
|
|
|
|
|
120 |
* ------------------------------------------------------------------------- |
* ------------------------------------------------------------------------- |
121 |
* Other commons needed by the minimization procedures: |
* Other commons needed by the minimization procedures: |
122 |
* ------------------------------------------------------------------------- |
* ------------------------------------------------------------------------- |
123 |
* computation of chi^2 and its derivatives |
* computation of chi^2 and its derivatives |
124 |
REAL*8 CHI2D,CHI2DD !first and second derivatives of chi^2 |
REAL*8 DAL |
125 |
|
REAL*8 CHI2D,CHI2DOLD,CHI2DD !first and second derivatives of chi^2 |
126 |
REAL*8 DXDAL,DYDAL !first derivatives of x and y on AL-pha |
REAL*8 DXDAL,DYDAL !first derivatives of x and y on AL-pha |
127 |
REAL*8 STEPAL ! |
REAL*8 STEPAL ! |
128 |
***PP*** REAL*8 ZINI !z coordinate of reference plane |
***PP*** REAL*8 ZINI !z coordinate of reference plane |
129 |
COMMON/MINICOM/ |
COMMON/MINICOM/ |
130 |
$ CHI2D(5),CHI2DD(5,5) |
$ DAL(5) |
131 |
|
$ ,CHI2D(5),CHI2DOLD(5),CHI2DD(5,5) |
132 |
$ ,DXDAL(NPLANES,5),DYDAL(NPLANES,5) |
$ ,DXDAL(NPLANES,5),DYDAL(NPLANES,5) |
133 |
$ ,STEPAL(5) |
$ ,STEPAL(5) |
134 |
***PP*** $ ,ZINI |
***PP*** $ ,ZINI |
135 |
|
|
136 |
* ------------------------------------------------------------------------- |
* ------------------------------------------------------------------------- |
137 |
* convergence conditions |
* convergence conditions |
138 |
REAL*8 ALTOL !tolerances on AL-pha parameters |
REAL*8 ALTOL !tolerances on AL-pha parameters |
142 |
$ ALTOL(5) |
$ ALTOL(5) |
143 |
$ ,ALMAX(5),ALMIN(5) |
$ ,ALMAX(5),ALMIN(5) |
144 |
$ ,ISTEPMAX |
$ ,ISTEPMAX |
145 |
|
|
146 |
|
|
147 |
* ------------------------------------------------------------------------- |
* ------------------------------------------------------------------------- |
148 |
* variables for the tracking (routine posxy) |
* variables for the tracking (routine posxy) |
149 |
REAL*8 VECT,VECTINI,VOUT !vectors for Runge-Kutta method |
REAL*8 VECT,VECTINI,VOUT !vectors for Runge-Kutta method |
153 |
COMMON/TRACKING/VECT(7),VECTINI(7),VOUT(7),TOLL,STEPMAX, |
COMMON/TRACKING/VECT(7),VECTINI(7),VOUT(7),TOLL,STEPMAX, |
154 |
$ CHARGE |
$ CHARGE |
155 |
|
|
156 |
|
* ------------------------------------------------------------------------- |
157 |
|
* variables for the DeltaB coorection |
158 |
|
REAL*8 DELTA0,DELTA1,DLT ! DELTA-DEF = DELTA0 + DELTA1 x Y |
159 |
|
COMMON/DELTAB/DELTA0,DELTA1,DLT |
160 |
|
|
161 |
|
|