Rebuilt hdgeant4 pion zapper last night. Submitted jobs around 9:00, and they sat in queue for a while. Just finishing now.
Looking at the resulting vertex files, definitely more pions decaying in the detector, but still probably not enough. Scale the decay length down to 50 and retry...
UPDATE: reran with the max length for each detector element set to 50. This did not produce many more pi- decays. Might have to rethink how I'm handling this...
This is even stranger now that I'm remembering that the lengths in hdgeant4 are in mm...
Resetting max length to 100 mm, but changing max track id from 3 to 6:
if (track.GetTrackID() < 6 &&
Rerunning...
UPDATE: Yep, that did it! Pretty much all of the pions are picked up by this change. Resetting the decay length to 500 mm and rerunning...
UPDATE: Looking good!
NOW that that's all tidied up, will need a new in-depth look at the fastpi MC. This will likely have consequences for the ways that vertexing and NN classification were working out, since the fastpi MC wasn't really and different from the regular ppim MC. Will also have to have Megan histogram the pion decay locations.
No comments:
Post a Comment