|
1 | 1 | # Overview |
2 | 2 |
|
3 | | -SBNDPDSAnalyzer is an art::analyzer that saves useful information for PDS MC/reco studies at different stages of the simulation: |
4 | | -- Gen: MCTruth |
5 | | -- G4: MCParticles and deposited energy |
6 | | -- Detsim: OpDetWaveforms |
7 | | -- Reco1: deconvolved waveforms, OpHits and OpFlashes |
8 | | -- Reco2: flash matchers scores (SimpleFlash and OpT0Finder) and CRUMBS score |
9 | | - |
10 | | -Two TTrees are created: |
11 | | -- PDSMapTree: Creates a tree with the optical channel IDs, position (X, Y and Z in cm) and the photon detector type (0: CoatedPMT, 1: UncoatedPMT, 2: VUV XARAPUCA and 3: VIS XARAPUCA). |
12 | | -- OpAnaTree: TTree with the analysis variables (described below). Each entry in the TTree corresponds to an event. |
13 | | -- OpAnaPerTrackTree: TTree with the number of SimPhotons produced by each MC particle. Each entry in the TTree corresponds to one MC particle in one event. Requires SimPhotons (not SimPhotonsLite). |
| 3 | +`SBNDPDSAnalyzer` is an `art::analyzer` that saves useful information for PDS MC/reco studies at different stages of the simulation: |
| 4 | +- Gen: `MCTruth` |
| 5 | +- G4: `MCParticles` and deposited energy |
| 6 | +- Detsim: `OpDetWaveforms` |
| 7 | +- Reco1: deconvolved waveforms (`OpDetWaveforms`), `OpHits` and `OpFlashes` |
| 8 | +- Reco2: flash matchers scores (`SimpleFlash` and `OpT0Finder`) and `CRUMBS` score |
| 9 | + |
| 10 | +Three TTrees are created: |
| 11 | +- `PDSMapTree`: Creates a tree with the optical channel IDs, position (X, Y and Z in cm) and the photon detector type (0: CoatedPMT, 1: UncoatedPMT, 2: VUV XARAPUCA and 3: VIS XARAPUCA). |
| 12 | +- `OpAnaTree`: TTree with the analysis variables (described below). Each entry in the TTree corresponds to an event. |
| 13 | +- `OpAnaPerTrackTree`: TTree with the number of SimPhotons produced by each MC particle. Each entry in the TTree corresponds to one MC particle in one event. Requires SimPhotons (not SimPhotonsLite). |
14 | 14 |
|
15 | 15 | Configuration parameters: |
16 | | -- SaveMCTruth, SaveMCParticles, ...: booleans to save different data products at different stages of the simulation/reconstruction chain. |
17 | | -- Verbosity: shows some printous while running the analyzer. |
18 | | -- MakePerTrackTree: option to make a tree with the number of SimPhotons per MCParticle (OpAnaPerTrackTree). |
19 | | -- MakePDSGeoTree: option to dump the PDS mapping in a TTree (PDSMapTree) |
20 | | -- UseSimPhotonsLite: use SimPhotonsLite or SimPhotons (default is true, i.e. use SimPhotonsLite) |
21 | | -- KeepPDGCode: specify what PDG codes will be stored in the TTree from the MCParticle list. If no PDGs are specified, it saves all the MCParticles (default is [], so it saves everything). |
22 | | -- MCTruthOrigin: specify the MCTruth origins that will be saved. Default is [1], i.e. neutrino generated events. |
23 | | -- MCTruthPDG: specify the PDG codes of the candidate vertex particles for the MCTruth. Defaults is [12, -12, 14, -14], i.e. BNB neutrinos. |
24 | | -- PDTypes: vector specifying the PD types for which it will save the SimPhotons (default is ["pmt_coated", "pmt_uncoated"]) |
| 16 | +- `SaveMCTruth`, `SaveMCParticles`, ...: booleans to save different data products at different stages of the simulation/reconstruction chain. |
| 17 | +- `Verbosity`: shows some printous while running the analyzer. |
| 18 | +- `MakePerTrackTree`: option to make a tree with the number of `SimPhotons` per `MCParticle` (`OpAnaPerTrackTree`). |
| 19 | +- `MakePDSGeoTree`: option to dump the PDS mapping in a TTree (`PDSMapTree`) |
| 20 | +- `UseSimPhotonsLite`: use `SimPhotonsLite` or `SimPhotons` (default is true, i.e. use `SimPhotonsLite`) |
| 21 | +- `KeepPDGCode`: specify what PDG codes will be stored in the TTree from the MC particles list. If no PDGs are specified, it saves all the MC particles (default is [], so it saves everything). |
| 22 | +- `MCTruthOrigin`: specify the `MCTruth` origins that will be saved. Default is [1], i.e. neutrino generated events. |
| 23 | +- `MCTruthPDG`: specify the PDG codes of the candidate vertex particles for the `MCTruth`. Defaults is [12, -12, 14, -14], i.e. BNB neutrinos. |
| 24 | +- `PDTypes`: vector specifying the PD types for which it will save the `SimPhotons` (default is ["pmt_coated", "pmt_uncoated"]) |
25 | 25 |
|
26 | | -# Variables in the TTree |
27 | 26 |
|
28 | 27 | # Variables in the TTree |
29 | 28 |
|
@@ -120,6 +119,8 @@ It dumps all the reconstructed OpHits in the event. |
120 | 119 |
|
121 | 120 | It dumps all the reconstructed OpFlash in the event and the associated OpHits associated to each OpFlash. Each entry in the following vectors correspond to one OpFlash. |
122 | 121 |
|
| 122 | +| Branch name | Type | Description | |
| 123 | +|-----------------------------|---------------------------------------------|--------------| |
123 | 124 | | nopflash | int | Total number of reconstructed OpFlash objects | |
124 | 125 | | flash_time| std::vector<double\> | t0 of the reconstructed OpFlashes | |
125 | 126 | | flash_total_pe | std::vector<double\> | Integrated (all optical channels) number of photoelectrons in each OpFlash | |
|
0 commit comments