Skip to content

CAF POT Histogram (develop) #597

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.

Sign up for GitHub

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

Merged
kjplows merged 5 commits into from
Nov 5, 2025
Merged

CAF POT Histogram (develop) #597

Changes from 2 commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
c0b9c7d
add histograms to file even when total events are 0 for POT accounting
lynnt20 Oct 24, 2025
78a102f
remove debug cout
lynnt20 Oct 27, 2025
c7dde1e
add functionality to save ngenevt on a run-level in new histogram
lynnt20 Nov 3, 2025
5ee3fb2
rename NGenEvents to TotalGenEvents to share similar syntax with POT …
lynnt20 Nov 3, 2025
63bd2b4
Merge branch 'develop' into feature/lynnt_cafpot_dev
kjplows Nov 5, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
51 changes: 25 additions & 26 deletions sbncode/CAFMaker/CAFMaker_module.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2802,29 +2802,26 @@ void CAFMaker::endSubRun(art::SubRun& sr) {

//......................................................................
void CAFMaker::endJob() {
if (fTotalEvents == 0) {

std::cerr << "No events processed in this file. Aborting rather than "
"produce an empty CAF."
// Only produce empty recTree/GenieTree since it relies on non-zero art events.
// Still want to keep POT histograms.
if (fTotalEvents == 0) {
std::cerr << "No events processed in this file. Producing empty recTree/GenieTree."
<< std::endl;
// n.b. changed abort() to return so that eny exceptions thrown during startup
// still get printed to the user by art
return;
}



if(fFile){

AddHistogramsToFile(fFile);
fRecTree->SetDirectory(fFile);
if(fGenieTree){
fGenieTree->BuildIndex("SourceFileHash", "GENIEEntry");
fGenieTree->SetDirectory(fFile);
}
if (fParams.CreateBlindedCAF()) {
fRecTreeb->SetDirectory(fFileb);
fRecTreep->SetDirectory(fFilep);
if (fTotalEvents > 0) {
if(fGenieTree){
fGenieTree->BuildIndex("SourceFileHash", "GENIEEntry");
fGenieTree->SetDirectory(fFile);
}
fRecTree->SetDirectory(fFile);
if (fParams.CreateBlindedCAF()) {
fRecTreeb->SetDirectory(fFileb);
fRecTreep->SetDirectory(fFilep);
}
}
fFile->cd();
fFile->Write();
Expand All @@ -2840,17 +2837,19 @@ void CAFMaker::endJob() {
}

if(fFlatFile){

AddHistogramsToFile(fFlatFile);
fFlatTree->SetDirectory(fFlatFile);
if(fFlatGenieTree){
fFlatGenieTree->BuildIndex("SourceFileHash", "GENIEEntry");
fFlatGenieTree->SetDirectory(fFlatFile);
}
if (fParams.CreateBlindedCAF() && fFlatFileb) {
fFlatTreeb->SetDirectory(fFlatFileb);
fFlatTreep->SetDirectory(fFlatFilep);
if (fTotalEvents > 0) {
if(fFlatGenieTree){
fFlatGenieTree->BuildIndex("SourceFileHash", "GENIEEntry");
fFlatGenieTree->SetDirectory(fFlatFile);
}
fFlatTree->SetDirectory(fFlatFile);
if (fParams.CreateBlindedCAF() && fFlatFileb) {
fFlatTreeb->SetDirectory(fFlatFileb);
fFlatTreep->SetDirectory(fFlatFilep);
}
}

fFlatFile->cd();
fFlatFile->Write();
if (fParams.CreateBlindedCAF()) {
Expand Down