[Halld-offline] Clarification of "primary"

[Beni Zihlmann]

Hi Matt,
here is my short term solution for your problem:

declare after the extenal variable:

#extern s_HDDM_t* thisInputEvent;


then you can access the number of initially thrown particles like this:

int a =  
thisInputEvent->physicsEvents->in[0].reactions->in[0].vertices->in[0].products->mult;

and you modify the line
showers->in[0].primary = (stack == 0);

to
showers->in[0].primary = (track <= a);

that will restore the value "primary" to 1 if the particle was initially 
thrown and 0 otherwise.

cheers,
Beni


> Hi all,
>
> While we debate how/if such a change can be implemented, can someone provide quick guidance for how to restore primary tag back to its proper functionality so I can continue some BCAL reconstruction studies?
>
> Is it as simple is commenting out some change that Beni made to gustep?  Or how do I gain access to s_Reactions_t* in the framework?
>
> I can try to fix this on my own, but it is not part of the code I'm familiar with -- I'd just prefer to have some expert guidance.
>
> Thanks!
>
> Matt
>
> On Mar 29, 2011, at 4:31 PM, Richard Jones wrote:
>
>> Beni and all,
>>
>> The meaning of "primary particle" is, I think, pretty well defined.  Why do we need to change language that everyone understands?
>>
>> This change confuses me, and I think screws up the particle stacking sematics of geant, not to speak of the clear separation between primaries and secondaries that exists in our minds -- and in geant as originally designed.  Instead of rewriting the dictionary, can you propose what additional features you would like to see in the hdgeant output, and let's have a discussion on how best to accomplish that?
>>
>> As it is, the changes you have introduced break the familiar pattern of how geant is supposed to work.
>>
>> -Richard J.
>>
>> _______________________________________________
>> Halld-offline mailing list
>> Halld-offline at jlab.org
>> https://mailman.jlab.org/mailman/listinfo/halld-offline
>
> _______________________________________________
> Halld-offline mailing list
> Halld-offline at jlab.org
> https://mailman.jlab.org/mailman/listinfo/halld-offline