-
Notifications
You must be signed in to change notification settings - Fork 16
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Feature/jlarkin numu selection #84
Conversation
Jacob, this looks like exactly what we talked about to me. I'll approve, but maybe give Chris a chance to look at the code before you merge |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks good to me. Couple of very minor comments.
const SpillVar dummy_var([&out](const caf::SRSpillProxy* spill){ | ||
for(const auto& slc: spill->slc) { | ||
if(kIcarus202208QELikeContainedMuon(&slc)) { | ||
out << "Run: " << spill->hdr.run << "\tSubrun: " << spill->hdr.subrun | ||
<< "\tEvent: " << spill->hdr.evt << "\tVertex: (" | ||
<< slc.vertex.x << ", " << slc.vertex.y << ", " << slc.vertex.z << ")\n"; | ||
return 1; | ||
} | ||
} | ||
return 0; | ||
}); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This works, and I guess I can't fault you for not using EventList since it doesn't exist in SBNAna yet ;)
See issue #56
#include "sbnana/SBNAna/Cuts/NumuCutsIcarus202208.h" | ||
#include "sbnana/SBNAna/Vars/NumuVarsIcarus202208.h" | ||
|
||
#include "sbnanaobj/StandardRecord/Proxy/FwdDeclare.h" |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
No need for the fwd declare when you're actually including the whole thing on the next line anyway.
|
||
namespace ana { | ||
|
||
static bool contained(const caf::SRTrackProxy& trk) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Put Icarus202208 here too before it clashes with something.
return PTrackInd; | ||
}); | ||
|
||
static bool proton_cut(const caf::SRTrackProxy& trk) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Ditto
Add new selection for Icarus Numu
I tried to add Icarus202208 in front of all names to avoid naming conflicts with other Vars and Cuts.
I also added a CAFAna macro which prints the run, subrun, and event number for events with selected slices. This macro is intended to find events for the scanning group to look at.