Re: r4881 - trunk/src/programs/Simulation/HDGeant

[David Lawrence <davidl@jlab.org>]

Hi Ryan,

    No, I don't think there is a better way in this case. Someone else 
may have a different idea (and can chime in if so), but for what you're 
trying to do, I think the path you've chosen is the most practical.

Regards,
-David

Ryan Mitchell wrote:
> Hi David,
>
> Great.  Thanks.  That does look useful for studying FCAL shower 
> profiles and is probably preferable to modifying gustep.F.  Another 
> thing we want to do though is revisit the FCAL radiation damage 
> simulations.  In that case we will have to generate a large number of 
> minimum bias and BEAM events and record energy depositions in small 
> chunks of the FCAL.  Rather than recording event-by-event energy 
> depositions (which would be impractical, I think) I plan to just keep 
> a running sum, with very minor modifications to what I just committed 
> in gustep.F.  Is there a better way to study radiation damage than 
> modifying gustep.F?
>
> Ryan
>
>
> On Feb 20, 2009, at 11:23 AM, David Lawrence wrote:
>
>>
>> Hi Ryan,
>>
>>   FYI: There is another mechanism available to do something similar 
>> to this using trajectory points, in case you find it useful. To use 
>> it, set the TRAJ card in the control.in file to 5 (see the comments 
>> in control.in). Then you will get a huge hdgeant.hddm file, but it 
>> will contain the info needed to make DMCTrajectoryPoint objects in 
>> DANA (see src/libraries/TRACKING/DMCTrajectoryPoint.h). Since this 
>> records every gustep call for every particle tracked in the event, 
>> you can access the full shower development info with just as much 
>> accuracy, but outside of the GEANt program.
>>
>>   It may not provide any advantage over the method you've already 
>> implemented, but it's just something to be aware of.
>>
>> Regards,
>> -David
>>
>> Hall-D.SVN.Repository@jlab.org wrote:
>>> Author: remitche
>>> Date: 2009-02-20 10:58:54 -0500 (Fri, 20 Feb 2009)
>>> New Revision: 4881
>>>
>>> Modified:
>>>   trunk/src/programs/Simulation/HDGeant/gustep.F
>>> Log:
>>> Added an ntuple (FCAL_SHOWER_PROFILE_NTUPLE) that divides the FCAL 
>>> into 1x1x1cm cubes and records energy deposition in each cube.  This 
>>> will give us more detailed information about FCAL shower shapes.
>>>
>>>
>>>
>>
>> -- 
>>
>> ------------------------------------------------------------------------
>> David Lawrence Ph.D.
>> Staff Scientist                 Office: (757)269-5567   [[[  [   [ 
>> [       Jefferson Lab                   Pager:  (757)584-5567   [  [ 
>> [ [ [ [   http://www.jlab.org/~davidl     davidl@jlab.org         
>> [[[  [[ [[ [[[
>> ------------------------------------------------------------------------
>>
>>
>

-- 

------------------------------------------------------------------------
 David Lawrence Ph.D.
 Staff Scientist                 Office: (757)269-5567   [[[  [   [ [       
 Jefferson Lab                   Pager:  (757)584-5567   [  [ [ [ [ [   
 http://www.jlab.org/~davidl     davidl@jlab.org         [[[  [[ [[ [[[
------------------------------------------------------------------------