Good news is that I think that the hddm to vertex file generator is working. Assembling files to send to Megan. Decided against putting the decaying particle names in the vertex files b/c it was too cumbersome to switch for the different topos AND it's inefficient. The particle id numbers in the file refer to the particle ids in the raw files.
In checking these files, it looks like only about 40% of the fast pion events have pions that decay (at all). This could be due to a) not having a short-enough lifetime or b) track.GetTrackID() < 3 not being met for the pions in MyPrimaryPionZapper.cc in hdgeant4.
Looking at the code, it is probably the former, as I currently have the path length set to be up to 5 m! Yikes:
return (rand() % 1000)/1000.0 * 500;
Turning this down to 1 m:
return (rand() % 1000)/1000.0 * 100;
... and rebuilding hdgeant4 and rerunning fastpi mc.
No comments:
Post a Comment