vitaly,
why don't you use lrn ? It should work for both data and montecarlo and it
should be the only one to use.
And what montecarlo events are you trying to generate anyways? It seems an
e'pn final state and the problems you report have to do with neutrons
pointing out of plane, at the edge of our acceptance.
On thing to check if whether geant generates secondaries off your
60 MeV/c neutron. That could be another explanation of why you get such
a strange event. I'd not try to reconstruct that..
then, the recon code may have a preference for n rather than p, but of
course it is because it accepts two tracks only. But what happens if you
use lrn ? What is the reconstructed charge/bending of your "neutron" ??
-- ________________________________________________________________________________ Tancredi Botto, phone: +1-617-253-9204 mobile: +1-978-490-4124 research scientist MIT/Bates, 21 Manning Av Middleton MA, 01949 ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^On Thu, 13 Nov 2003, Vitaliy Ziskin wrote:
> 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