Hi Chris, Vitaliy,
Like I said some time ago, I did (and do) all my analysis from ep_skim,
which does not have "tml" - MC generated theta angle - but only has the
reconstructed quantities (like real data)...
So I never looked at what actually goes in, just what goes out. I
generated the ntuple by default and never looked at it - good thing I did
not delete it!
Adrian
On Mon, 17 Apr 2006, Christopher Crawford wrote:
> Hi Vitaliy,
> I didn't generate it, I just used Adrian's ntuple; see his previous message
> and the log file in my previous email for the details. Thanks for looking at
> it.
> --Chris
> _______________________________________
>
> TA-53/MPF-1/D111 P-23 MS H803
> LANL, Los Alamos, NM 87545
> 505-665-9804(o) 665-4121(f) 662-0639(h)
> _______________________________________
>
>
> On Apr 17, 2006, at 16:54:38, Vitaliy Ziskin wrote:
>
>> Chris,
>> Can you tell me what channel you used to generate these events. Also, can
>> you send me your kine.init file and the output lines when you ran the code.
>> My guess that all events in the tail are somehow in this peak around 0. I
>> have no idea where this peak around 0 comes from and I'm surprised that
>> Adrian did not see it.
>>
>> Cheers, Vitaliy
>>
>> Christopher Crawford wrote:
>>> I took an initial look at the radiative corrections in MC, and there are
>>> problems. I include two problem histograms: the first one is just raw
>>> theta in the left sector. It has a significant peak with structure around
>>> theta=0, which should not be present. If I had to guess, I would say that
>>> the code is generating radiative kinematics in the frame of the outgoing
>>> electron instead of in the BLAST frame (although that does not make sense
>>> either, since you do see a small tail above v=0.10 GeV^2).
>>> The second is a histogram of the invariant mass W for 25 < th_e < 35 deg.
>>> Out of the 266406 events in the elastic peak (v < .010 GeV^2 or W-M <
>>> 5.6 MeV), there were only 45 events in the radiative tail. The original
>>> code gives delta(.010)=0.804 and delta(.4427)=0.978, implying
>>> 266406*(978-804)/978 = 47397 events in the tail, up to only W=1150 MeV.
>>> Note: v = W^2 - M^2 is the inelasticity (a parameter in MASCARAD).
>>> Due to small statistics, I was unable to verify the shape of the tail.
>>> Anyways, someone must correct these problems before this MC can be trusted
>>> for radiative corrections. I will report on corrections based on the
>>> original MASCARAD code at the next meeting.
>>> --Chris
>>> _______________________________________
>>> TA-53/MPF-1/D111 P-23 MS H803
>>> LANL, Los Alamos, NM 87545
>>> 505-665-9804(o) 665-4121(f) 662-0639(h)
>>> _______________________________________
>>> ------------------------------------------------------------------------
>>
>> --
>> ***********************************************************************
>> Dr. Vitaliy Ziskin 77 Massachusetts Ave
>> Research Scientist NW13-202
>> Nuclear Science and Engineering Cambridge, Ma, 02134
>> MIT (617)253-5466
>> vziskin@mit.edu
>> ***********************************************************************
-------------------------------------------------------------------------- Adrian Sindile Room 205A, DeMeritt Hall Research Scientist phone: (603)862-1217 Xenon Imaging Center & FAX: (603)862-2998 Nuclear Physics Group email: asindile@unh.edu University of New Hampshire http://einstein.unh.edu/~adrian/ ------------------------------------------------------------------------
This archive was generated by hypermail 2.1.2 : Mon Feb 24 2014 - 14:07:33 EST