/[PAMELA software]/DarthVader/TrackerLevel2/inc/F77/common_mini_2.f
ViewVC logotype

Contents of /DarthVader/TrackerLevel2/inc/F77/common_mini_2.f

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.13 - (show annotations) (download)
Thu Dec 4 14:33:48 2008 UTC (16 years ago) by pam-fi
Branch: MAIN
CVS Tags: v9r00, v9r01, v6r01, v6r00
Changes since 1.12: +2 -6 lines
NEW TRACKING ALGORITHM (for tracks with few hits and for nuclei)

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

  ViewVC Help
Powered by ViewVC 1.1.23