Skip to content
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

Merged
merged 4 commits into from
Sep 4, 2022

Conversation

jacoblarkin
Copy link
Contributor

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.

@etworcester
Copy link
Contributor

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

Copy link
Contributor

@cjbacchus cjbacchus left a 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.

Comment on lines +22 to +32
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;
});
Copy link
Contributor

@cjbacchus cjbacchus Aug 4, 2022

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"
Copy link
Contributor

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)
Copy link
Contributor

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)
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Ditto

@miquelnebot miquelnebot merged commit 780acca into develop Sep 4, 2022
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

4 participants