Skip to content

Inherit genie config from genie_beam_settings.fcl in sbncode #667

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

Merged
bear-is-asleep merged 3 commits into from
Apr 1, 2025
Merged

Inherit genie config from genie_beam_settings.fcl in sbncode #667

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
b5cda3a
Inherit genie config from genie_beam_settings.fcl in sbncode
marcodeltutto Feb 26, 2025
a0dbf33
Re-add rotated bucket config
marcodeltutto Mar 25, 2025
387e2cf
Accidentally removed flux config in previous commit.
marcodeltutto Mar 27, 2025
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
97 changes: 21 additions & 76 deletions sbndcode/LArSoftConfigurations/gen/genie_sbnd.fcl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -11,8 +11,15 @@
# * GENIE: configuration of the GENIEGen module, to be included in the job
#
#
# Comments:
# Bunch sigma taken from
# https://beamdocs.fnal.gov/AD/DocDB/0050/005000/001/bunchLength_1st_draft.pdf
# Bunch spacing taken from
# https://inspirehep.net/files/610a942fd8632bbbca2c8ad90da86670
# GlobalTimeOffset: 0 # Brailsford 2017/10/09 Simulation currently only reads out one drift frame so having a 1.6ms offset (one drift frame) means almost all events get placed outside the readout window. We COULD make the readout window 2 or 3 drift frames long but we also have a disk space problem (we don't have any free). The solution is to remove the offset and keep one frame readout

#include "genie.fcl"
#include "genie_beam_settings.fcl"

BEGIN_PROLOG

Expand Down Expand Up @@ -255,48 +262,6 @@ sbnd_flux_numi_dirt: {
}


################################################################################
### Bucket structure configuration (Implemented by A. C. Ezeribe)
###
#
# These are FCL tables that setup arranging the neutrino interactions into BNB buckets.
# Use them in its configuration as:
#
# physics: {
# producers: {
# generator: {
# module_type: GenieGen
# # ...
# @table: fluxbucket_bub
# }
# }
# }
#
# or, if needed to override an existing configuration, with the less
# straightforward:
#
# physics.producers.generator: {
# @table::physics.producers.generator
# @table: fluxbucket_bnb
# }
#
#
# Booster Neutrino Beam, bucket
#
# Bunch sigma taken from
# https://beamdocs.fnal.gov/AD/DocDB/0050/005000/001/bunchLength_1st_draft.pdf
# Bunch spacing taken from
# https://inspirehep.net/files/610a942fd8632bbbca2c8ad90da86670

sbnd_fluxbucket_bnb: {
SpillTimeConfig: "evgb::EvtTimeFNALBeam booster, dtbucket=18.936, sigma=1.308"
}

# Booster Neutrino Beam, rotated bucket
#
sbnd_fluxbucket_bnb_rotated: {
SpillTimeConfig: "evgb::EvtTimeFNALBeam booster, dtbucket=18.936, sigma=1.308"
}

################################################################################
### GENIE module configurations
Expand All @@ -319,48 +284,25 @@ sbnd_fluxbucket_bnb_rotated: {
###
###

sbnd_genie: {
@table::standard_genie
@table::sbnd_flux_bnb_nu
@table::sbnd_fluxbucket_bnb
BeamName: "booster"
# GlobalTimeOffset: 1.28e6 #sbnd reads out 1.6ms before the spill
GlobalTimeOffset: 0 #Brailsford 2017/10/09 Simulation currently only reads out one drift frame so having a 1.6ms offset (one drift frame) means almost all events get placed outside the readout window. We COULD make the readout window 2 or 3 drift frames long but we also have a disk space problem (we don't have any free). The solution is to remove the offset and keep one frame readout
#EventGeneratorList: "Default+CCMEC+NCMEC" # Only needed to generate partial samples in GENIE v3
DefinedVtxHistRange: true
VtxPosHistRange: [ -210, 210, -210, 210, -10, 510 ]
AddGenieVtxTime: true #Enable adding time of flight to neutrino vertex
} # sbnd_genie

sbnd_genie_hist: {
@table::sbnd_genie
@table::sbnd_flux_bnb_nu_Bv1_hist
BeamCenter: [0., 0., 0.]
BeamDirection: [0., 0., 1.]
BeamRadius: 4.
}


#
# Basic configuration (BnB beam)
# Basic configurations (BNB beam)
#

sbnd_genie_simple: {
@table::sbnd_genie
@table::sbn_genie_BNB_base
@table::sbnd_flux_bnb_nu
EventsPerSpill: 0
POTPerSpill: 5e12
MaxFluxFileMB: 100
}


sbnd_genie_simple_dirt: {
@table::sbnd_genie
@table::sbn_genie_BNB_base
@table::sbnd_flux_bnb_dirt
EventsPerSpill: 0
POTPerSpill: 5e12
MaxFluxFileMB: 100
TopVolume: "volWorld"
#Start the flux rays at 18m upstream of the TPC frontface. Chosen as this is the distance a muon of 8 GeV (max flux sim. energy) can travel
# Start the flux rays at 18m upstream of the TPC frontface.
# Chosen as this is the distance a muon of 8 GeV (max flux sim. energy) can travel
FluxUpstreamZ: -18
}

Expand All @@ -379,15 +321,18 @@ sbnd_genie_simple_rockbox: {
GeomScan: "file: GENIE/sbnd_rock_maxpathlength_fluxI_gdmlv02_00.xml"
}

# Rotated bucket configuration
sbnd_fluxbucket_bnb_rotated: {
SpillTimeConfig: "evgb::EvtTimeFNALBeam booster, dtbucket=18.936, sigma=1.308"
}

#
# Basic configurations (NuMI beam)
#

sbnd_genie_simple_numi: {
@table::standard_genie
@table::sbn_genie_NuMI_base
@table::sbnd_flux_numi
BeamName: "numi"
GlobalTimeOffset: 1.6e6
POTPerSpill: 3e13
EventsPerSpill: 0
}

sbnd_genie_simple_numi_dirt: {
Expand Down