Started working on raw event generator. Functionality to generate files with physical ratios of the charged decay modes is working, but I somehow broke the part of the code that assigns each particle's decay vertex index and parent particle index....
Update... seems to be working now.
Also working is addition of decay mechanism flag into the data stream. This required changing the raw file structure slightly AND adding this feature to the sld_raw_2_hddm.py code.
This may also need to be added later on so that the mech info stays with the recon files. Will check it out.
BUT! At this point, I think that we're ready to gen some MC!
Generating 500 files (5M events) of pim_gam MC...
And also 500 files worth of charged mixed MC (#yolo). This should produce an expected 1.5 sl_mu events per file for roughly 750 events total (yikes, that's grim) in the raw. Detector effects will knock this down by a factor of 10. So, it's likely that I'll need to generate roughly 100 times more at some point to actually do this.
No comments:
Post a Comment