[G12] Lepton Simulation
Craig Bookwalter
craigb at jlab.org
Tue Oct 16 01:28:13 EDT 2012
OK Paolone is right -- for each track in the PART sector 0 bank in the
.part file, you need to set its vertex to be the MVRT vertex from the
original reconstructed data.
Code I used for the feedback study is here:
svn co https://jlabsvn.jlab.org/svnroot/clas/users/craigb/gsimfbk
`set_part_vertex.cpp' is the one. You can read the code yourself, but you
give it a BOS file on the command line, and it implicitly expects to find
an ascii ntuple named `vertices.dat' in the cwd where each line is:
<eventno> <proton vertex> <pip vtx> <pim vtx>
where I was too lazy to write x, y, and z for each of the vertices. So 7
entries per line, one line per event.
Hopefully that can get you started.
--cb
> Hey MK -- just read your mail. I'm nosing around to try and figure out
> what exactly I did.
>
>
>> Hey MK,
>>
>> I'm not entirely sure how Craig did it, but here's my 2 cents. You might
>> need to modify your GAMP output to have a z vertex, and then modify
>> gamp2part to propagate that vertex into the PART bank (instead of a
>> randomly assigned z within the target cell). If this is what Craig did,
>> he may still have his modified gamp2part.
>>
>> -Michael
>>
>>
>>> Greetings G12,
>>>
>>> I am attempting to do a study similar to Craig Bookwalter's GSIM study.
>>> http://clasweb.jlab.org/rungroups/g12/wiki/index.php/Simulation
>>>
>>> Instead of exclusive proton pi+ pi- events, I will be using proton e+
>>> e-
>>> events.
>>>
>>> In his study he mentions, "
>>>
>>>
>>> Vertex correction
>>>
>>> It was pointed out that these events should be generated with their
>>> vertices the same as reconstructed by a1c from the real data--in the
>>> first pass the vertex was randomly distributed within the target.
>>>
>>> "
>>>
>>> Can someone point out to me how to feed GSIM with the vertex
>>> information.
>>> The event I am using are data, that have been corrected with beam
>>> corrections and momentum corrections, so the events I pass through are
>>> in the form of a GAMP file.
>>>
>>> Thanks
>>> MK
>>> _______________________________________________
>>> G12 mailing list
>>> G12 at jlab.org
>>> https://mailman.jlab.org/mailman/listinfo/g12
>>>
>>
>>
>> --
>> Michael Paolone
>> Assistant Research Professor
>> Temple University
>> Phone: (215)-204-7240
>> Email: mpaolone at jlab.org
>>
>>
>> _______________________________________________
>> G12 mailing list
>> G12 at jlab.org
>> https://mailman.jlab.org/mailman/listinfo/g12
>>
>
>
> _______________________________________________
> G12 mailing list
> G12 at jlab.org
> https://mailman.jlab.org/mailman/listinfo/g12
>
More information about the G12
mailing list