DUNE / 2x2_sim

Wrappers for ArgonCube 2x2 simulation at NERSC
Apache License 2.0
5 stars 10 forks source link

Limit each hit segment to one contributor #20

Closed mjkramer closed 1 year ago

mjkramer commented 1 year ago

As suggested by Zach, update the edep-sim macro to include something like

/edep/hitSeparation volLArActive -1 mm

This should ensure that each hit segment has only one contributing trajectory, which would be good, since we only save one contributor during the HDF5 conversion, so multiple contributors would lead to a loss of information.

https://github.com/ClarkMcGrew/edep-sim/commit/49e0ac9bd6b4ac58c3800009039335d0f2cfcdfd https://github.com/ClarkMcGrew/edep-sim/blob/0a6975949368dc50c4802f1511f4bbd52d5b7cf9/src/EDepSimSegmentSD.cc#L77 https://github.com/ClarkMcGrew/edep-sim/blob/0a6975949368dc50c4802f1511f4bbd52d5b7cf9/src/EDepSimHitSegment.cc#L96

mjkramer commented 1 year ago

Confirmed (for the first file of MiniRun4) that, with this change, each hit in the LAr is only associated with one trajectory. Prior to this change, there were up to four trajectories associated with each hit in the LAr.

We still can have multiple trajectories per hit in MINERvA, but since they aren't using the HDF5 files (where we truncate the list of contributors), that shouldn't be a problem. If we do want to make this same change for MINERvA, we will either have to painfully repeat this line for every SensDet in MINERvA, or change the default hitSeparation in edep-sim.