addnoise_findprim_snanas.fcl

Go to the documentation of this file.

# addnoise_findprim_daqsimanas.fcl

# Read in (presumably noise-free) waveforms, add noise and pedestals,
# then find trigger primitives on the (now noisy) waveforms, and save
# the trigger primitives (as "Hits") to an output TTree.

# Any number of trigger primitive finders can be run: useful for
# running the same finder with multiple parameter settings, for
# example. The output of each trigger primitive finder should be
# connected to a DaqSimAna analyzer

#include "SNAna.fcl"
#include "trigprim.fcl"
BEGIN_PROLOG
trigprim_snana: {
   @table::standard_snana
   SaveNeighbourADCs:    false
   SaveTruth:            false
}
END_PROLOG
#include "opticaldetectormodules_dune.fcl"
#include "standard_detsim_dune10kt_1x2x6.fcl"


process_name: PrimSim

services: @local::dunefd_simulation_services
services.TFileService: { fileName: "SNAna_plus_trigprim_multithreshold.root" }
services.TimeTracker:       {}
services.MemoryTracker:         { } # default is one
services.RandomNumberGenerator: {} #ART native random number generator
// services.FileCatalogMetadata:  @local::art_file_catalog_mc
services.NuRandomService:       @local::per_event_NuRandomService  # seedservice.fcl
services.message:              @local::dune_message_services_prod
services.Geometry: @local::dune10kt_1x2x6_geo
#source is now a root file
source: {
  module_type: RootInput
  maxEvents:  -1        # Number of events to create
}

# Define and configure some modules to do work on each event.
# First modules are defined; they are scheduled later.
# Modules are grouped by type.
physics: {
   producers: {
      simwire: {
         module_type: "SimWireDUNE"
         SimChannelLabel: "largeant"
         NoiseOn: true
         PedestalOn: true
         DistortOn: false
         SuppressOn: false
         KeepEmptyChannels: true
         AdcSimulator:      "adcsim_ideal"
         UseRawDigitInput: true
         RawDigitInputLabel: "daq"
      }

      ophit:              @local::dune35t_ophit 

      # You can add any number of trigger primitive finders here. Make
      # sure each one has a distinct name, and has a corresponding
      # daqsimana analyzer otherwise you won't get a corresponding
      # TTree in the output file. Also make sure to add it to the "simulate" list below
      trigprim12: @local::trigprimpass1
      trigprim13: @local::trigprimpass1
      trigprim14: @local::trigprimpass1
      trigprim15: @local::trigprimpass1
      trigprim16: @local::trigprimpass1
      
      trigprim30sig: @local::trigprimpass2
      trigprim40sig: @local::trigprimpass2
      trigprim45sig: @local::trigprimpass2
      trigprim50sig: @local::trigprimpass2
      trigprim55sig: @local::trigprimpass2
      trigprim60sig: @local::trigprimpass2

      rns:            { module_type: "RandomNumberSaver" }
   }
   
   analyzers: {
      # Create an SNAna for each trigger primitive finders. Make
      # sure to add each one to the "ana" list below
      snanatrigprim12:    @local::trigprim_snana
      snanatrigprim13:    @local::trigprim_snana
      snanatrigprim14:    @local::trigprim_snana
      snanatrigprim15:    @local::trigprim_snana
      snanatrigprim16:    @local::trigprim_snana
      
      snanatrigprim30sig: @local::trigprim_snana
      snanatrigprim40sig: @local::trigprim_snana
      snanatrigprim45sig: @local::trigprim_snana
      snanatrigprim50sig: @local::trigprim_snana
      snanatrigprim55sig: @local::trigprim_snana
      snanatrigprim60sig: @local::trigprim_snana
   }

   simulate: [rns, simwire, ophit,
              trigprim12,         trigprim13,         trigprim14,         trigprim15,         trigprim16,
              trigprim30sig,      trigprim40sig,      trigprim45sig,      trigprim50sig,      trigprim55sig,      trigprim60sig]
   ana     : [snanatrigprim12,    snanatrigprim13,    snanatrigprim14,    snanatrigprim15,    snanatrigprim16,
              snanatrigprim30sig, snanatrigprim40sig, snanatrigprim45sig, snanatrigprim50sig, snanatrigprim55sig, snanatrigprim60sig]
   
   trigger_paths: [ simulate ] 
   end_paths: [ ana ] 

}

# Modify the settings of each of the trigger primitive finders. If you
# want to use a different trigger primitive finder class, set
# "physics.producers.<itemname>.finder.tool_type: MyNewTriggerPrimitiveFinderTool"
physics.producers.trigprim12.finder.Threshold: 1200
physics.producers.trigprim13.finder.Threshold: 1300
physics.producers.trigprim14.finder.Threshold: 1400
physics.producers.trigprim15.finder.Threshold: 1500
physics.producers.trigprim16.finder.Threshold: 1600
physics.producers.trigprim30sig.finder.ThresholdInSigma: 300
physics.producers.trigprim40sig.finder.ThresholdInSigma: 400
physics.producers.trigprim45sig.finder.ThresholdInSigma: 450
physics.producers.trigprim50sig.finder.ThresholdInSigma: 500
physics.producers.trigprim55sig.finder.ThresholdInSigma: 550
physics.producers.trigprim60sig.finder.ThresholdInSigma: 600


# Set the HitLabel to make sure each DaqAnaSim reads hits from its
#  corresponding TriggerPrimitiveFinder
physics.analyzers.snanatrigprim12.HitLabel: "trigprim12"
physics.analyzers.snanatrigprim13.HitLabel: "trigprim13"
physics.analyzers.snanatrigprim14.HitLabel: "trigprim14"
physics.analyzers.snanatrigprim15.HitLabel: "trigprim15"
physics.analyzers.snanatrigprim16.HitLabel: "trigprim16"
physics.analyzers.snanatrigprim30sig.HitLabel: "trigprim30sig"
physics.analyzers.snanatrigprim40sig.HitLabel: "trigprim40sig"
physics.analyzers.snanatrigprim45sig.HitLabel: "trigprim45sig"
physics.analyzers.snanatrigprim50sig.HitLabel: "trigprim50sig"
physics.analyzers.snanatrigprim55sig.HitLabel: "trigprim55sig"
physics.analyzers.snanatrigprim60sig.HitLabel: "trigprim60sig"

# Only save IDEs in one of the output trees, since they're always the same and will only take up space otherwise
# physics.analyzers.snanatrigprim800.SaveIDEs: true

# =====================================================================
# You probably won't need to change anything below here

# Shut up the backtracker
services.message.destinations.LogStandardOut.categories.BackTracker.limit: 0
services.message.destinations.LogStandardError.categories.BackTracker.limit: 0

# Display all the INFO message for SNAna_module
services.message.destinations.LogStandardOut.categories.SNAna_module.limit:-1

# If you want to have the mf::LogDebug() messages, uncomment the following:
# services.message.destinations.LogStandardOut.threshold:"DEBUG"
# services.message.destinations.LogStandardOut.categories.SNAna_module.limit:-1
# services.message.destinations.LogStandardOut.categories.default.limit:0


# DetSim services. These are used by SimWire to add noise and
# pedestals to the incoming (noise-free) waveforms
services.SimChannelExtractService: @local::scxgeneric
services.ChannelNoiseService:      @local::chnoiseold
services.PedestalAdditionService:  @local::padprovided
services.AdcDistortService:        @local::stuckbits
services.AdcSuppressService:       @local::zsnone
services.AdcCompressService:       @local::cmpreplace
services.DetPedestalService:       @local::dune_fixedpeds

#Optical detector configuration
physics.producers.ophit.InputModule: "opdigi"