[BLAST_ANAWARE] monte-carlo recon woes

From: Vitaliy Ziskin (vziskin@mit.edu)
Date: Thu Nov 13 2003 - 11:48:34 EST


I'm having serious difficulties with the reconstructing the monte-carlo
generated events. Here is what I mean:

Take even #20:

nsed says:
mc 0 -1 : 0.631280 55.494602 1.257000 -0.086000 | 3 electron
0.690000 0.800000 0.000000
mc 1 1 : 0.664320 44.117699 179.764008 -0.086000 | 14 proton
0.690000 0.800000 0.000000
mc 2 0 : 0.061150 75.356102 -166.990005 -0.086000 | 13 neutron
0.690000 0.800000 0.000000

recon says:
************************************************************
* Row * pml * pmr * pwl * pwr *
************************************************************
* 20 * 0.6312800 * 0.0611499 * 0.6412974 * 0.8690107 *

This time gets electron momentum rather well in both "pml" and "pwl".
 However "pmr" assumes momentum of a neutron (???) and "pwr" reconstruct
god knows what

Next event #82:

nsed says:
mc 0 -1 : 0.761120 24.690100 3.544000 -5.241700 | 3 electron
0.690000 0.800000 0.000000
mc 1 0 : 0.050170 106.363701 -17.608000 -5.241700 | 13 neutron
0.690000 0.800000 0.000000
mc 2 1 : 0.402160 64.585701 -179.197006 -5.241700 | 14 proton
0.690000 0.800000 0.000000

recons says:
82 * 0.0501700 * 0.4021600 * 0.4482221 * 0.6218268 *

Here "pml" assumes a value of a neutron (????) and "pwl" is again junk.
 "pmr" correctly takes a value of a proton but "pwr" is again shit.

Most of the time when the "pml" is wrong it asummes the value of a
neutron momentum. However, the reconstruction itself is all over the
place.

This was done for electron on the left, but the results are the same for
the electron on the right. This precludes me from doing any sort of
monte-carlo studies ahead of Dr. Arenhovel's visit. Can some
reconstruction wizards look at it fairly soon?

                                                          Cheers, Vitaliy



This archive was generated by hypermail 2.1.2 : Mon Feb 24 2014 - 14:07:30 EST