Skip to content

Adding fcl files to produce STM samples at both VD116 and VD101 #521

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

Closed
YongyiBWu wants to merge 1 commit into from
Closed

Adding fcl files to produce STM samples at both VD116 and VD101 #521

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
166 changes: 166 additions & 0 deletions JobConfig/pileup/STM/BeamTo2VD.fcl
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,166 @@
# Re-sample EleBeamCat and MuBeamCat datasets, propagate them to a chosen virtual detector, and write out the events with SimParticles and StepPointMCs travelling through that virtual detector only
# Original author: Yuri Oksuzian, 2019
# Updated for MDC2020 (DetectorSteps): D. Brown
# Updated for STM studies: P. Plesniak

# Offline includes
#include "Offline/CommonMC/fcl/prolog.fcl"
#include "Offline/fcl/standardServices.fcl"
#include "Offline/STMMC/fcl/prolog.fcl"

# Production includes
#include "Production/JobConfig/common/prolog.fcl"
#include "Production/JobConfig/pileup/prolog.fcl"
#include "Production/JobConfig/pileup/STM/prolog.fcl"


process_name: BeamToVD

source : {
module_type : EmptyEvent
maxEvents : @nil
}

services : @local::Services.Sim

physics: {
producers : {
@table::Common.producers
@table::Pileup.producers
extractVirtualDetectorSteps101 : {
# Creates "mu2e::StepPointMCs_extractVirtualDetectorSteps__BeamToVD" products that only store StepPointMCs in the selected virtual detector
# Want to only keep events with the given virtual detector ID and filter out the rest
# Empty collections will be generated if the selected virtual detector has no associated StepPointMCs to allow the compression to take place
module_type : STMResamplingProducer
StepPointMCsTag : @local::ResamplingProducer.StepPointMCsTag
VirtualDetectorID : 101
}
extractVirtualDetectorSteps116 : {
module_type : STMResamplingProducer
StepPointMCsTag : @local::ResamplingProducer.StepPointMCsTag
VirtualDetectorID : 116
}
compressDetStepMCsSTM101 : {
# Creates "mu2e::SimParticlemv_compressDetStepMCsSTM__BeamToVD" and "mu2e::StepPointMCs_compressDetStepMCsSTM__BeamToVD" products from "mu2e::StepPointMCs_extractVirtualDetectorSteps__BeamToVD" products
# Keeps SimParticles and StepPointMCs going through the selected virtual detector only
module_type : CompressDetStepMCs
strawGasStepTag : ""
caloShowerStepTag : ""
surfaceStepTag : ""
crvStepTag : ""
simParticleTags : ["g4run"]
debugLevel : 0
stepPointMCTags : ["extractVirtualDetectorSteps101"]
compressionOptions : {
strawGasStepCompressionLevel: "noCompression"
caloShowerStepCompressionLevel: "noCompression"
crvStepCompressionLevel: "noCompression"
surfaceStepCompressionLevel: "noCompression"
simParticleCompressionLevel : "noCompression"
stepPointMCCompressionLevel : "noCompression"
keepNGenerations : -1
mcTrajectoryCompressionLevel : "noCompression"
}
mcTrajectoryTag : ""
}
compressDetStepMCsSTM116 : {
# Creates "mu2e::SimParticlemv_compressDetStepMCsSTM__BeamToVD" and "mu2e::StepPointMCs_compressDetStepMCsSTM__BeamToVD" products from "mu2e::StepPointMCs_extractVirtualDetectorSteps__BeamToVD" products
# Keeps SimParticles and StepPointMCs going through the selected virtual detector only
module_type : CompressDetStepMCs
strawGasStepTag : ""
caloShowerStepTag : ""
surfaceStepTag : ""
crvStepTag : ""
simParticleTags : ["g4run"]
debugLevel : 0
stepPointMCTags : ["extractVirtualDetectorSteps116"]
compressionOptions : {
strawGasStepCompressionLevel: "noCompression"
caloShowerStepCompressionLevel: "noCompression"
crvStepCompressionLevel: "noCompression"
surfaceStepCompressionLevel: "noCompression"
simParticleCompressionLevel : "noCompression"
stepPointMCCompressionLevel : "noCompression"
keepNGenerations : -1
mcTrajectoryCompressionLevel : "noCompression"
}
mcTrajectoryTag : ""
}
}
filters : {
@table::Common.filters
@table::Pileup.filters
filterVirtualDetectorSteps101: {
# Filters events based on the size of "mu2e::StepPointMCs_compressDetStepMCsSTM__BeamToVD"
# If there are no StepPointMCs in the selected virtual detector, the event is not kept
# Required as extractVirtualDetectorSteps will create a StepPointMCCollection even if the collection is empty
module_type: STMResamplingFilter
StepPointMCsTag : "compressDetStepMCsSTM101"
}
filterVirtualDetectorSteps116: {
# Filters events based on the size of "mu2e::StepPointMCs_compressDetStepMCsSTM__BeamToVD"
# If there are no StepPointMCs in the selected virtual detector, the event is not kept
# Required as extractVirtualDetectorSteps will create a StepPointMCCollection even if the collection is empty
module_type: STMResamplingFilter
StepPointMCsTag : "compressDetStepMCsSTM116"
}
}
analyzers : {
@table::Common.analyzers
countVirtualDetectorHits : {
# Runs a counter over the virtual detectors to keep track of hits
module_type : CountVirtualDetectorHits
stepPointMCsTag : @local::VirtualDetectorCounter.StepPointMCsTag
virtualDetectorIDs : [88, 89, 90, 100, 101, 116]
}
}
# setup paths
STMCompressedPath101 : [@sequence::Pileup.beamResamplerSequence, @sequence::Common.g4Sequence, extractVirtualDetectorSteps101, compressDetStepMCsSTM101, filterVirtualDetectorSteps101]
STMCompressedPath116 : [@sequence::Pileup.beamResamplerSequence, @sequence::Common.g4Sequence, extractVirtualDetectorSteps116, compressDetStepMCsSTM116, filterVirtualDetectorSteps116]
trigger_paths: [STMCompressedPath101, STMCompressedPath116]
outPathCompressed : [genCountLogger, countVirtualDetectorHits, compressedOutput101, compressedOutput116]
end_paths: [outPathCompressed]
}

# Keep only the necessary data
outputs: {
compressedOutput101 : {
module_type: RootOutput
outputCommands: [
"drop *_*_*_*",
@sequence::STMPileup.BeamToVDKeptProducts
]
SelectEvents: [STMCompressedPath101]
fileName : "dts.owner.BeamToVD101.version.sequencer.art"
}
compressedOutput116 : {
module_type: RootOutput
outputCommands: [
"drop *_*_*_*",
@sequence::STMPileup.BeamToVDKeptProducts
]
SelectEvents: [STMCompressedPath116]
fileName : "dts.owner.BeamToVD116.version.sequencer.art"
}
}

# Point Mu2eG4 to the pre-simulated data
physics.producers.g4run.inputs: {
primaryType: "StepPoints"
primaryTag: "beamResampler:Beam"
inputMCTrajectories: ""
simStageOverride: 1
inputPhysVolumeMultiInfo: "beamResampler"
updateEventLevelVolumeInfos: {
input: "beamResampler:eventlevel"
outInstance: "eventlevel"
}
}
physics.producers.g4run.SDConfig.preSimulatedHits: ["beamResampler:virtualdetector"]
physics.producers.g4run.SDConfig.enableSD: [virtualdetector]
physics.producers.g4run.Mu2eG4CommonCut: {}

# Epilog includes
# comment out of MT issue #include "Production/JobConfig/common/MT.fcl"
#include "Production/JobConfig/common/epilog.fcl"
#include "Production/JobConfig/pileup/epilog.fcl"
168 changes: 168 additions & 0 deletions JobConfig/pileup/STM/BeamTo2VD1809.fcl
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,168 @@
# Generate and propagate 1809 keV muon stop signal photons from muon target stops
# Original author: Dave Brown, LBNL
# Updated for STM studies: Pawel Plesniak

# Offline includes
#include "Offline/fcl/standardServices.fcl"
#include "Offline/EventGenerator/fcl/prolog.fcl"
#include "Offline/STMMC/fcl/prolog.fcl"

# Production includes
#include "Production/JobConfig/common/prolog.fcl"
#include "Production/JobConfig/pileup/prolog.fcl"
#include "Production/JobConfig/primary/prolog.fcl"
#include "Production/JobConfig/pileup/STM/prolog.fcl"

process_name: BeamToVD1809

source: {
module_type : EmptyEvent
maxEvents : @nil
}

services : @local::Services.Sim

physics : {
producers : {
@table::Common.producers
@table::Pileup.producers
generate : { # Generates the 1809 keV peaks from muon target stop positions
module_type: Pileup
inputSimParticles: TargetStopResampler
stoppingTargetMaterial : "Al"
verbosity: 0
captureProducts: [{tool_type : "MuCap1809keVGammaGenerator"}]
decayProducts: []
}
extractVirtualDetectorSteps101 : {
# Creates "mu2e::StepPointMCs_extractVirtualDetectorSteps__BeamToVD1809" products that only store StepPointMCs in the selected virtual detector
# Want to only keep events with the given virtual detector ID and filter out the rest
# Empty collections will be generated if the selected virtual detector has no associated StepPointMCs to allow the compression to take place
module_type : STMResamplingProducer
StepPointMCsTag : @local::ResamplingProducer.StepPointMCsTag
VirtualDetectorID : 101
}
extractVirtualDetectorSteps116 : {
module_type : STMResamplingProducer
StepPointMCsTag : @local::ResamplingProducer.StepPointMCsTag
VirtualDetectorID : 116
}
compressDetStepMCsSTM101 : {
# Creates "mu2e::SimParticlemv_compressDetStepMCsSTM__BeamToVD" and "mu2e::StepPointMCs_compressDetStepMCsSTM__BeamToVD" products from "mu2e::StepPointMCs_extractVirtualDetectorSteps__BeamToVD" products
# Keeps SimParticles and StepPointMCs going through the selected virtual detector only
module_type : CompressDetStepMCs
strawGasStepTag : ""
caloShowerStepTag : ""
surfaceStepTag : ""
crvStepTag : ""
simParticleTags : ["g4run"]
debugLevel : 0
stepPointMCTags : ["extractVirtualDetectorSteps101"]
compressionOptions : {
strawGasStepCompressionLevel: "noCompression"
caloShowerStepCompressionLevel: "noCompression"
crvStepCompressionLevel: "noCompression"
surfaceStepCompressionLevel: "noCompression"
simParticleCompressionLevel : "noCompression"
stepPointMCCompressionLevel : "noCompression"
keepNGenerations : -1
mcTrajectoryCompressionLevel : "noCompression"
}
mcTrajectoryTag : ""
}
compressDetStepMCsSTM116 : {
# Creates "mu2e::SimParticlemv_compressDetStepMCsSTM__BeamToVD" and "mu2e::StepPointMCs_compressDetStepMCsSTM__BeamToVD" products from "mu2e::StepPointMCs_extractVirtualDetectorSteps__BeamToVD" products
# Keeps SimParticles and StepPointMCs going through the selected virtual detector only
module_type : CompressDetStepMCs
strawGasStepTag : ""
caloShowerStepTag : ""
surfaceStepTag : ""
crvStepTag : ""
simParticleTags : ["g4run"]
debugLevel : 0
stepPointMCTags : ["extractVirtualDetectorSteps116"]
compressionOptions : {
strawGasStepCompressionLevel: "noCompression"
caloShowerStepCompressionLevel: "noCompression"
crvStepCompressionLevel: "noCompression"
surfaceStepCompressionLevel: "noCompression"
simParticleCompressionLevel : "noCompression"
stepPointMCCompressionLevel : "noCompression"
keepNGenerations : -1
mcTrajectoryCompressionLevel : "noCompression"
}
mcTrajectoryTag : ""
}
}
filters : {
@table::Common.filters
@table::Pileup.filters
TargetStopResampler : @local::Primary.filters.TargetStopResampler
filterVirtualDetectorSteps101: {
# Filters events based on the size of "mu2e::StepPointMCs_compressDetStepMCsSTM__BeamToVD"
# If there are no StepPointMCs in the selected virtual detector, the event is not kept
# Required as extractVirtualDetectorSteps will create a StepPointMCCollection even if the collection is empty
module_type: STMResamplingFilter
StepPointMCsTag : "compressDetStepMCsSTM101"
}
filterVirtualDetectorSteps116: {
# Filters events based on the size of "mu2e::StepPointMCs_compressDetStepMCsSTM__BeamToVD"
# If there are no StepPointMCs in the selected virtual detector, the event is not kept
# Required as extractVirtualDetectorSteps will create a StepPointMCCollection even if the collection is empty
module_type: STMResamplingFilter
StepPointMCsTag : "compressDetStepMCsSTM116"
}
}
analyzers : {
@table::Common.analyzers
countVirtualDetectorHits : {
# Runs a counter over the virtual detectors to keep track of hits
module_type : CountVirtualDetectorHits
stepPointMCsTag : @local::VirtualDetectorCounter.StepPointMCsTag
virtualDetectorIDs : [88, 89, 90, 100, 101, 116]
}
}
PileupPath101 : [TargetStopResampler, @sequence::Common.generateSequence, @sequence::Common.g4Sequence, extractVirtualDetectorSteps101, compressDetStepMCsSTM101, filterVirtualDetectorSteps101]
PileupPath116 : [TargetStopResampler, @sequence::Common.generateSequence, @sequence::Common.g4Sequence, extractVirtualDetectorSteps116, compressDetStepMCsSTM116, filterVirtualDetectorSteps116]
trigger_paths : [PileupPath101, PileupPath116]
outPathCompressed : [genCountLogger, countVirtualDetectorHits, compressedOutput101, compressedOutput116]
end_paths: [outPathCompressed]
}
outputs : {
compressedOutput101 : {
module_type: RootOutput
outputCommands: [
"drop *_*_*_*",
@sequence::STMPileup.BeamToVDKeptProducts
]
SelectEvents: [PileupPath101]
fileName : "dts.owner.BeamToVD101.version.sequencer.art"
}
compressedOutput116 : {
module_type: RootOutput
outputCommands: [
"drop *_*_*_*",
@sequence::STMPileup.BeamToVDKeptProducts
]
SelectEvents: [PileupPath116]
fileName : "dts.owner.BeamToVD116.version.sequencer.art"
}
}

# resampling configuration
physics.producers.g4run.inputs: {
primaryType: StageParticles
primaryTag: "generate"
inputMCTrajectories: ""
simStageOverride: 2
inputPhysVolumeMultiInfo: "TargetStopResampler"
updateEventLevelVolumeInfos: {
input: "TargetStopResampler:eventlevel"
outInstance: "eventlevel"
}
}

# Epilog includes
# comment out of MT issue #include "Production/JobConfig/common/MT.fcl"
#include "Production/JobConfig/common/epilog.fcl"
#include "Production/JobConfig/pileup/epilog.fcl"
Loading