Fixes to make it work!

Author: PetrilloAtWork

boosteddmanalysis/Reconstruction/SelectFinalState.cxx

@@ -28,16 +28,14 @@ bdm::SelectFinalState::Selection_t bdm::SelectFinalState::select
   for (std::size_t index = 0U; index < nParticles; ++index) {
     auto const& particle = particles[index];
 
-    std::cout << "Particle #" << index << " mother: " << particle.Mother() << std::endl;
-    
     // TODO
 
     // filter neutrons out
     if (!fKeepNeutrons && (std::abs(particle.PdgCode()) == kAbsNeutron))
       continue;
 
     // exclude the particles that are not "primary" (i.e. not generated)
-    if (particle.Mother() < 0) continue;
+    if (particle.Mother() > 0) continue;
 
     selected.push_back(index);
 

boosteddmanalysis/Reconstruction/bdmparticlereconstructionchain.fcl

@@ -0,0 +1,156 @@
+#
+# File:    bdmparticlereconstructionchain.fcl
+# Brief:   Runs Boosted Dark Matter particle selection and reconstruction algorithms.
+# Author:  Gianluca Petrillo (petrillo@slac.stanford.edu)
+# Date:    20181006
+# Version: 1.0
+# 
+# This is an example of "full chain" where generated particles are selected
+# and "reconstructed".
+# 
+# Input:
+# - collection of simb::MCParticle
+#
+# Output:
+# - collection of bdm::SmearedMCParticle
+# - text file with the content of the original MCParticle list
+# - text file with the content of the reconstructed SmearedMCParticle
+#
+#
+# Changes:
+# 20181006 (petrillo@fnal.gov) [v1.0]
+#   original version
+#
+#
+
+#include "seedservice.fcl"
+#include "particleselectors.fcl"
+#include "smearedreconstructionsettings.fcl"
+
+#-------------------------------------------------------------------------------
+process_name: RecoSmear
+
+#-------------------------------------------------------------------------------
+# no source specified; use -s option from command line
+
+services: {
+  
+  message: {
+  #   debugModules: [ "*" ]
+    destinations: {
+      
+      # grab all the "DumpSmearedMCParticles" messages
+      # and put them in DumpSmearedMCParticles.log
+      DumpSmearedMCParticles: {
+        append: false
+        categories: {
+          DumpSmearedMCParticles: { limit: -1 }
+          default: { limit: 0 }
+        }
+        filename: "DumpSmearedMCParticles.log"
+        threshold: "INFO"
+        type: "file"
+      } # DumpSmearedMCParticles
+      
+      # grab all the "DumpMCParticles" messages
+      # and put them in DumpMCParticles.log
+      DumpMCParticles: {
+        append: false
+        categories: {
+          DumpMCParticles: { limit: -1 }
+          default: { limit: 0 }
+        }
+        filename: "DumpMCParticles.log"
+        threshold: "INFO"
+        type: "file"
+      } # DumpMCParticles
+      
+      LogStandardOut: {
+        categories: {
+          DumpMCParticles:        { limit: 0 }
+          DumpSmearedMCParticles: { limit: 0 }
+          default:                { limit: -1 }
+        }
+        threshold: "INFO"
+        type: "cout"
+      } # LogStandardOut
+      
+    } # destinations
+  } # message
+  
+  RandomNumberGenerator: {}
+  NuRandomService:       @local::per_event_NuRandomService
+  
+} # services
+
+
+#-------------------------------------------------------------------------------
+physics: {
+  
+  producers: {
+  
+    selector: {
+      module_type: ParticleSelector
+      
+      particles: "largeant"
+      
+      selector: @local::standard_select_final_state
+      
+      
+    } # selector
+    
+    smear: {
+      module_type: SmearedReconstruction
+      
+      particles: "selector"
+      
+      reconstruction: @local::smeared_reconstruction_dune10kt
+      
+    } # smear
+    
+  } # producers
+  
+  analyzers: {
+    
+    dumpmcparticles: {
+      
+      module_type:  DumpMCParticles
+      
+      OutputCategory: "DumpMCParticles"
+      InputParticles: "largeant"
+      
+    } # dumpsmearedmcparticles
+    
+    dumpsmearedmcparticles: {
+      
+      module_type:  DumpSmearedMCParticles
+      
+      OutputCategory: "DumpSmearedMCParticles"
+      InputParticles: "smear"
+      
+    } # dumpsmearedmcparticles
+    
+  } # analyzers
+  
+  
+  reco: [ selector, smear ]
+  dumpers: [ output ]
+  stream: [ dumpmcparticles, dumpsmearedmcparticles ]
+  
+  trigger_paths: [ reco ]
+  end_paths: [ stream, dumpers ]
+  
+} # physics
+
+
+#-------------------------------------------------------------------------------
+outputs: {
+  output: {
+    module_type: RootOutput
+    
+    fileName: "%ifb_%tc_%p.root"
+    
+  } # output
+} # outputs
+
+#-------------------------------------------------------------------------------

boosteddmanalysis/Reconstruction/dump_smearedmcparticles.fcl

@@ -69,7 +69,7 @@ physics: {
       # output category ("OutputCategory" by default), useful for filtering (see above)
       OutputCategory: "DumpSmearedMCParticles"
 
-      # specify the label of the simb::MCParticle data product (or producer)
+      # specify the label of the bdm::SmearedMCParticle data product (or producer)
       InputParticles: "smear"
 
     } # dumpsmearedmcparticles

boosteddmanalysis/Reconstruction/particleselectors.fcl

@@ -19,6 +19,8 @@
 #
 #
 
+BEGIN_PROLOG
+
 ################################################################################
 standard_select_final_state: {
 
@@ -28,3 +30,5 @@ standard_select_final_state: {
 
 
 ################################################################################
+
+END_PROLOG