Skip to content

Light Calorimetry: add SRLightCalo and SRTrueDeposit classes #181

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

Open
lynnt20 wants to merge 1 commit into
base: develop
Choose a base branch
from
+64 −1
Open

Light Calorimetry: add SRLightCalo and SRTrueDeposit classes #181

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
27 changes: 27 additions & 0 deletions sbnanaobj/StandardRecord/SRLightCalo.h
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,27 @@
/**
* @file sbnanaobj/StandardRecord/SRLightCalo.h
* @brief Defines data structures to store light calorimetry info (header only).
* @author Lynn Tung
* @date November 19, 2025
*/

#ifndef SBNANAOBJ_STANDARDRECORD_SRLIGHTCALO_H
#define SBNANAOBJ_STANDARDRECORD_SRLIGHTCALO_H

#include "sbnanaobj/StandardRecord/SRConstants.h"

namespace caf {

/**
* @brief A struct to store the reconstructed visible energy, light, and charge in a recob::Slice.
*
*/
struct SRLightCalo {
float charge = caf::kSignalingNaN; // the total charge in the slice (on bestplane) [# electrons]
float light = caf::kSignalingNaN; // the reconstructed total number of photons [# photons]
float energy = caf::kSignalingNaN; // (charge+light)*scaling, [GeV]
int bestplane = caf::kUninitializedInt; // plane with most number of hits
};
}

#endif
3 changes: 3 additions & 0 deletions sbnanaobj/StandardRecord/SRSlice.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -19,6 +19,7 @@
#include "sbnanaobj/StandardRecord/SRNuID.h"
#include "sbnanaobj/StandardRecord/SRConstants.h"
#include "sbnanaobj/StandardRecord/SRCVNScore.h"
#include "sbnanaobj/StandardRecord/SRLightCalo.h"

#include <climits>

Expand Down Expand Up @@ -54,6 +55,8 @@ namespace caf

SRCorrectedOpFlash correctedOpFlash; //!< OpFlash corrected with using tpc information matched to this slice

SRLightCalo lightcalo; // !< Reconstructed visible energy, light, and charge

SRFakeReco fake_reco;

bool is_clear_cosmic { false }; //!< Whether pandora marks the slice as a "clear" cosmic
Expand Down
26 changes: 26 additions & 0 deletions sbnanaobj/StandardRecord/SRTrueDeposit.h
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,26 @@
/**
* @file sbnanaobj/StandardRecord/SrTrueDeposit.h
* @brief Defines data structures to store summed SimEnergyDeposits (header only).
* @author Lynn Tung
* @date November 18, 2025
*/

#ifndef SBNANAOBJ_STANDARDRECORD_SRTRUEDEPOSIT_H
#define SBNANAOBJ_STANDARDRECORD_SRTRUEDEPOSIT_H

#include "sbnanaobj/StandardRecord/SRConstants.h"

namespace caf {

/**
* @brief A struct to store true deposited energy, charge, and photons for each MCTruth interaction.
*
*/
struct SRTrueDeposit {
float electrons = caf::kSignalingNaN; ///< # of electrons at the site of energy deposition (after recombination, before drift)
float photons = caf::kSignalingNaN; ///< # of photons at the site of energy deposition (after recombination, before propagation)
float energy = caf::kSignalingNaN; ///< total **deposited energy** [GeV]
};
}

#endif
2 changes: 2 additions & 0 deletions sbnanaobj/StandardRecord/SRTruthBranch.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -5,6 +5,7 @@
#define SRTRUTHBRANCH_H

#include "sbnanaobj/StandardRecord/SRTrueInteraction.h"
#include "sbnanaobj/StandardRecord/SRTrueDeposit.h"
#include "sbnanaobj/StandardRecord/SRMeVPrtl.h"

#include <vector>
Expand All @@ -19,6 +20,7 @@ namespace caf
~SRTruthBranch();

std::vector<SRTrueInteraction> nu; ///< Vector of true nu or cosmic
std::vector<SRTrueDeposit> dep; ///< Vector of true energy deposits (summed), same size as nu
size_t nnu; ///< Number of true nu or cosmic

std::vector<SRMeVPrtl> prtl; ///< If present -- information on decay of MeV "Portal" particle
Expand Down
7 changes: 6 additions & 1 deletion sbnanaobj/StandardRecord/classes_def.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -95,6 +95,7 @@
</class>

<class name="caf::SROpT0Finder" />
<class name="caf::SRLightCalo" />

<class name="caf::SRSliceRecoBranch" ClassVersion="13">
<version ClassVersion="13" checksum="4211357366"/>
Expand Down Expand Up @@ -174,6 +175,8 @@
<version ClassVersion="10" checksum="4056415501"/>
</class>

<class name="caf::SRTrueDeposit"/>

<class name="caf::SRTrueParticlePlaneInfo" ClassVersion="11">
<version ClassVersion="11" checksum="317847253"/>
<version ClassVersion="10" checksum="1689803547"/>
Expand Down Expand Up @@ -215,7 +218,8 @@
<version ClassVersion="10" checksum="1609299342"/>
</class>

<class name="caf::SRTruthBranch" ClassVersion="11">
<class name="caf::SRTruthBranch" ClassVersion="12">
<version ClassVersion="12" checksum="3925767433"/>
<version ClassVersion="11" checksum="1553456992"/>
<version ClassVersion="10" checksum="2608785459"/>
</class>
Expand Down Expand Up @@ -336,6 +340,7 @@
<class name="std::vector<caf::SRTrueInteraction>" />
<class name="std::vector<caf::SRMultiverse>" />
<class name="std::vector<caf::SRTrueParticle>" />
<class name="std::vector<caf::SRTrueDeposit>" />
<class name="std::vector<caf::SRParticleMatch>" />
<class name="std::vector<caf::SRTrueCaloPoint>" />
<class name="std::vector<caf::SRCaloPoint>" />
Expand Down