[Halld-offline] Clarification of "primary"

[Matthew Shepherd]

Hi Beni,

This seems to help some, but still doesn't completely solve the problem.  If I plot the generated energy for photons tagged primary, I still have a bit of a peak at low energies (even though the photons were generated uniform in energy from 0-2 GeV).

Richard said I can revert to the old behavior by commenting out some lines, but what lines do I need to comment out in order to do this?

Thanks,

Matt


On Mar 30, 2011, at 8:54 AM, Beni Zihlmann wrote:

> 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
> 
> _______________________________________________
> Halld-offline mailing list
> Halld-offline at jlab.org
> https://mailman.jlab.org/mailman/listinfo/halld-offline