hi,
there were some questions about how the newton fitter convergence
depended on the initial fast fit (also if it was possible to converge to
two different solutions), so i checked into the matter.
first of all, i found a couple of bugs: the most important one was in
the determination of chi^2 = sum_18 (sim_i - data_i). the effect of the
bug was chi^2 = sum_18 ((sim_i-data_i)*2^(18-i)), so that the inner
hits were entirely discounted in the chi^2. also, we had previously
reported chi^2 ~ 1 (200micron resolution), but now, i am seeing between
30-50, which makes sense, seeing we need a better calibration of the
wire chambers.
i am going to recrunch old data with the new fast fitter and bug
fixes, and see how this affects our resolution.
there was also a bug in the magnetic field, allowing access to points
off the grid, causing segfaults, which has been fixed
the result of my study was that the newton fitter will converge to the
same value starting from on the order of almost any valid track
parameter. the convergence can be fine-tuned with various parameters
labeled FitTrack.*, Newt.* and NewtDer.* in the blastrc file. in
particular, i think NewtDer.dZ0 was too small. that would be a good
project for someone to fine-tune all of these parameters.
anyways, here are plots of convergence in the most dynamic variables,
(p,th) for an electron and proton, taken at random.
--chris
electron:
-----------
// run 1070, event 8, nw 0
float data[18] =
{ 89.3713, 90.2627, 91.1635, 97.8851, 98.8512, 99.7311,
136.753, 138.067, 139.116, 147.394, 148.378, 149.558,
193.022, 194.501, 196.027, 205.61, 207.085, 208.483 };
float v[6] =
{ p=0.768798, theta=27.6792, phi=-169.675, x0=-59.1049, y0=-10.7683,
z0=110.499 };
int pid=-1;
0 theta 20 40 60 80 100
120 140
p
+---------+---------+---------+---------+---------+---------+---------+-----
0.16:
......^.......^.....^^.^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^.^....^^^^^.^..^.:....
0.18:
.....^.^.........^^^...^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^......^^..^^^.^.......
0.20:
..........^...^^.^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^.^.........^...^.^^.^:....
0.22:
...^^^^.....^.^.^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^............^.......^.^..
0.24:
....^..^....^..^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^..........................:
0.26:
...^.^.^....^.^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^.^^................^...^....:^.
0.28:
...^......^^..^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^.......................^.....
0.30:
.^..^^.^.^^..^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^..............................
0.32:
.......^^^..^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^...............................
0.34:
.....^.^..^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^........................^.....:.
0.36:
...^...^^..^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^.................................
0.38:
........^..^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^...............................^..
0.40:
..^...^.....^^^^^^^^^^^^^^^^^^^^^^^^^^^^^...............................:...
0.42:
......^....^^^^^^^^^^^^^^^^^^^^^^^^^^^^^....................................
0.44:
...^^.^^.^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^...................................^.
0.46:
.^......^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^......................................
0.48:
^......^.^^^^^^^^^^^^^^^^^^^^^^^^^^^^.......................................
0.50:
.^^...^.^^^^^^^^^^^^^^^^^^^^^^^^^^^^........................................
0.52:
^^....^..^^^^^^^^^^^^^^^^^^^^^^^^^^.........................................
0.54:
......^.^^^^^^^^^^^^^^^^^^^^^^^^^^.........................................:
0.56:
.....^..^^^^^^^^^^^^^^^^^^^^^^^^^.......................................^...
0.58:
.....^..^^^^^^^^^^^^^^^^^^^^^^^^............................................
0.60:
........^^^^^^^^^^^^^^^^^^^^^.^.............................................
0.62:
.^...^.^.^^^^^^^^^^^^^^^^^^^^^^.............................................
0.64:
.....^.^^^^^^^^^^^^^^^^^^^^^^^..............................................
0.66:
.^...^.^^^^^^^^^^^^^^^^^^^^^^...............................................
0.68:
.....^.^^^^^^^^^^^^^^^^^^^^^................................................
0.70:
.....^.^^^^^^^^^^^^^^^^^^^^^...............................................^
0.72:
.....^^.^^^^^^^^^^^^^^^^^^^.................................................
0.74:
.^...^^^^^^^^^^^^^^^^^^^^^.............................................^....
0.76:
.^..^^^^^^^^^^^X^^^^^^^^^^..................................................
0.78:
^....^^^^^^^^^^^^^^^^^^^^................................................^..
0.80:
^....^^.^^^^^^^^^^^^^^^^^........................^..........................
0.82:
^.......^^^^^^^^^^^^^^^^....................................................
0.84:
.^....^^^^^^^^^^^^^^^^^^....................................................
0.86:
....^.^^^^^^^^^^^^^^^^^....................................................^
0.88:
....^^^^^^^^^^^^^^^^^^^.....................................................
0.90:
......^^^^^^^^^^^^^^^^......................................................
0.92:
....^^^^^^^^^^^^^^^^^^......................................................
0.94:
.^..^.^^^^^^^^^^^^^^^^......................................................
0.96:
^...^^^^^^^^^^^^^^^^^.......................................................
0.98:
^^..^^^^^^^^^^^^^^^^^.......................................................
proton
---------
run 3070, event 2, nw 1
float data[18] =
{ 6.81864, 6.70208, 6.9777, 5.20358, 5.21103, 5.24986,
7.26853, 7.25334, 7.18858, 4.13989, 3.92184, 4.33837,
4.51546, 4.43925, 4.36508, 0.282379, 0.196165, 0.0577528 };
float v[6] =
{p= 0.34338, theta=71.1494, phi=172.919, x0=-82.72, y0=10.2755,
z0=30.5711 };
int pid=1;
0 theta 20 40 60 80 100
120 140
+---------+---------+---------+---------+---------+---------+---------+-----
0.00:
............................................................................
0.02:
.......xxx.x.x...xxxxxxxx....x.....x............x.......xxx..xx.x...........
0.04:
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxx.xx.xxxxxxx.....xxx..xxxx.xx.xxx.x..xx..xx....
0.06:
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx..xxxxx...x..xxx.x..x..xxx...x.x..xx.x
0.08:
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx...x.x..x..x......x....x...x...x...
0.10:
.xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx........xx..x..x.x.xx..x..x...x...
0.12:
.x.xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx........xx.x...x..xx.............
0.14:
.x..xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx...x...x........x..x.x..........
0.16:
..x.x.xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx........x..x..............xx....
0.18:
...xxx.xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx.......xx...x..x...xxx..........
0.20:
....xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx.......x...x.......x...........
0.22:
......xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx..x....x.....x..x..x..x.......x
0.24:
......xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx.......x.......x.......x.x.....
0.26:
.......xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx...........x.............xx...
0.28:
..........xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx.........x.xx...x.............
0.30:
........xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx.......x.x.x..x............x..
0.32:
..........xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx..x...xx......................
0.34:
.........xxxxxxxxxxxxxxxxxxxxxxxxxOxxxxxxxxxxx..............................
0.36:
...........xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx.......x......................
0.38:
..........xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx......x.....x.................
0.40:
...........xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx..........x...............x..x
0.42:
...........xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx........x.x...................
0.44:
............xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx..........x...................
0.46:
............xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx.............x................
0.48:
............xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx.................x........x...
0.50:
............xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx.............................
0.52:
.............xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx.............................
0.54:
..............xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx...........x.................
0.56:
..............x.xxxxxxxxxxxxxxxxxxxxxxxxxxxxxx.......x..................x...
0.58:
...............xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx..............................
0.60:
................xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx.............x...............
0.62:
................xxxxxxxxxxxxxxxxxxxxxxxxxxxxxx..............................
0.64:
..................xxxxxxxxxxxxxxxxxxxxxxxxx.xx..............................
0.66:
...................xxxxxxxxxxxxxxxxxxxxxxxxxxx..............................
0.68:
...................x.xxxxxxxxxxxxxxxxxxxxxxxxx..............................
0.70:
.....................xx.xxxxxxxxxxxxxxxxxxxxxx..............................
0.72:
..........................xx.....x..........................................
0.74:
......................................x.....................................
0.76:
............................................................................
0.78:
........................................................................x...
0.80:
........................................................................x...
0.82:
............................................................................
0.84:
............................................................................
0.86:
............................................................................
0.88:
............................................................................
0.90:
............................................................................
0.92:
............................................................................
0.94:
............................................................................
0.96:
..........................................................................x.
0.98:
............................................................................
1.00:
............................................................................
This archive was generated by hypermail 2.1.2 : Mon Feb 24 2014 - 14:07:29 EST