hi chi,
epel does not generate its own vertex, and relies on kivert. i think
you were running and old version of epel, which included a call to
kivert. that call has since been removed for compatibility w/version2,
which was generating it's own vertex (one was used in coda mc, another
in geant).
i think the solution is to move 'kivert' up into 'gukine' like it is
in v2, since most generators use it by default, and yours turns it off
anyways. note: as it stands, the call to 'kivert' doesn't do anything
in 'usergeninter' since it is called after 'gskine', which registers the
vertex.
--chris
zhangchi wrote:
>Hi, I am not sure what is the problem, I RAN epel with v3 and I can show
>you plots.
>
>Make sure in fort file VERT is set appropriately. If you want to generate
>vertex in generator, there is a subroutine in blast_proc.f called GUVERT
>which I recommend, but be sure to set VERT to -1.
>
>Otherwise, set VERT as Aaron indicated in the comment of fort.99.
>
>Sorry, I am in New York right now to get my passport replaced before it
>expires in May.
>
>Chi
>
>On Mon, 10 Mar 2003, Chris Crawford wrote:
>
>
>
>>hi jason, and chi,
>> i figured out what broke the 'epel' generator in blastmc v3. it
>>gives no tracks, because the call to 'kivert' in 'USERGENINTER' comes
>>after it is already needed in 'gukine'.
>> jason, to fix it, you can just move the line 'CALL
>>KIVERT(NVTX,0.0,0)' to the beginning of 'USERGENINTER'.
>> chi, is there any reason not to have v3 the same as v2, where
>>'kivert' is called back in 'gukine'? if not, i'll make the switch and
>>check it in.
>>--chris
>>
>>
>>
>>Jason Seely wrote:
>>
>>
>>
>>>hey guys,
>>>
>>>i know i keep hassling you for help with the monte-carlo and
>>>generator, but i have to bother you again. i am trying to get the ep
>>>elastic generator working. i know that it's generating vertex
>>>information and passing it along to the monte-carlo, and if i write to
>>>an ascii file, it manages to print the vertex information out.
>>>however, that's all it does. there are no error messages or anything,
>>>it's just not putting the events through the detectors. the
>>>monte-carlo doesn't draw any tracks, or output any hits or anything.
>>>so i think the generator itself is fine, but either it's not passing
>>>all the information necessary, or something else with the interface
>>>withe BLAST. i noticed that, for example, the variable 'IPART1' is
>>>always '0', when the generator is actually producing protons ('14')
>>>and electrons ('3'). does anyone have any ideas? i'm kind of stuck
>>>right now.
>>>
>>>thanks in advance,
>>>
>>>jason
>>>
>>>
>>>
>>
>>
This archive was generated by hypermail 2.1.2 : Mon Feb 24 2014 - 14:07:29 EST