genie_dune.fcl

Go to the documentation of this file.

#include "genie.fcl"

BEGIN_PROLOG

# beam neutrino configurations

dunefd_genie:                         @local::standard_genie
dunefd_genie.BeamName:                "lbnf"
dunefd_genie.EventsPerSpill:          1
dunefd_genie.POTPerSpill:             0
dunefd_genie.GlobalTimeOffset:        0
dunefd_genie.RandomTimeOffset:        0
dunefd_genie.EventGeneratorList:      "Default"

dunefd_genie_hist:                    @local::dunefd_genie
dunefd_genie_hist.FluxType:           "histogram"
dunefd_genie_hist.FluxFiles:          ["g4lbne_v3r2p4b_FHC_FD_RIK.root"]
dunefd_genie_hist.FluxCopyMethod:     "DIRECT"
dunefd_genie_hist.FluxSearchPaths:    "/cvmfs/dune.osgstorage.org/pnfs/fnal.gov/usr/dune/persistent/stash/TaskForce_Flux/GenieHistFluxFiles/"
dunefd_genie_hist.BeamCenter:         [0., 0., -1.] 
dunefd_genie_hist.BeamDirection:      [0, 0.101, 0.995] # beam angle is 101mrad
dunefd_genie_hist.BeamRadius:         7. #sqrt(pow(3.6,2)+pow(6,2)
dunefd_genie_hist.TopVolume:          "volCryostat"

dunefddphase_genie_hist:             @local::dunefd_genie_hist
dunefddphase_genie_hist.BeamCenter:         [0., 0., -1.] 
dunefddphase_genie_hist.BeamDirection:      [0.101, 0, 0.995] # beam angle is 101mrad
dunefddphase_genie_hist.BeamRadius:         7. 
dunefddphase_genie_hist.FiducialCut:  "mbox:-610,-160,-10,610,160,610"

dunefd_genie_simple_neutrino:         @local::dunefd_genie
dunefd_genie_simple_neutrino.TopVolume:        "volDetEnclosure"

dune35t_genie_simple_neutrino:                @local::dunefd_genie
dune35t_genie_simple_neutrino.TopVolume:      "volCryostat"

dune35t_genie_hist:                           @local::dunefd_genie
dune35t_genie_hist.BeamCenter:                [100., 50., 0.]
dune35t_genie_hist.BeamDirection:             [0., 0., 0.]

dune_fd_genie_nu_simple_flux_workspace_window: @local::dunefd_genie
dune_fd_genie_nu_simple_flux_workspace_window.FluxType: "simple_flux"
dune_fd_genie_nu_simple_flux_workspace_window.FluxSearchPaths: "/cvmfs/dune.osgstorage.org/pnfs/fnal.gov/usr/dune/persistent/stash/TaskForce_Flux/ref_01/neutrino/gsimple/far/nominal/"
dune_fd_genie_nu_simple_flux_workspace_window.FluxFiles: ["gsimple_LBNE-FarDet_Workspace_g4lbne_v3r3p9_QGSP_BERT_Nominal_neutrino_*root"]
dune_fd_genie_nu_simple_flux_workspace_window.FluxCopyMethod: "DIRECT"
dune_fd_genie_nu_simple_flux_workspace_window.TopVolume: "volCryostat"

dune_fd_genie_nu_simple_flux_workspace_window_rotated_tpc: @local::dune_fd_genie_nu_simple_flux_workspace_window
dune_fd_genie_nu_simple_flux_workspace_window_rotated_tpc.FluxSearchPaths: "/cvmfs/dune.osgstorage.org/pnfs/fnal.gov/usr/dune/persistent/stash/TaskForce_Flux/ref_01/neutrino/gsimple/far/windowRotated10DegreesAroundY/"
dune_fd_genie_nu_simple_flux_workspace_window_rotated_tpc.FluxFiles: ["gsimple_LBNE-FarDet_Workspace_Rotated_Window_g4lbne_v3r3p9_QGSP_BERT_Nominal_neutrino_*.root"]

dune_fd_genie_anu_simple_flux_workspace_window: @local::dune_fd_genie_nu_simple_flux_workspace_window
dune_fd_genie_anu_simple_flux_workspace_window.FluxSearchPaths: "/cvmfs/dune.osgstorage.org/pnfs/fnal.gov/usr/dune/persistent/stash/TaskForce_Flux/ref_01/antineutrino/gsimple/far/nominal/"
dune_fd_genie_anu_simple_flux_workspace_window.FluxFiles: ["gsimple_LBNE-FarDet_Workspace_g4lbne_v3r3p9_QGSP_BERT_Nominal_antineutrino_*.root"]

dune_fd_genie_anu_simple_flux_workspace_window_rotated_tpc: @local::dune_fd_genie_nu_simple_flux_workspace_window
dune_fd_genie_anu_simple_flux_workspace_window_rotated_tpc.FluxSearchPaths: "/cvmfs/dune.osgstorage.org/pnfs/fnal.gov/usr/dune/persistent/stash/TaskForce_Flux/ref_01/antineutrino/gsimple/far/windowRotated10DegreesAroundY/"
dune_fd_genie_anu_simple_flux_workspace_window_rotated_tpc.FluxFiles: ["gsimple_LBNE-FarDet_Workspace_Rotated_Window_g4lbne_v3r3p9_QGSP_BERT_Nominal_antineutrino_*.root"]

#Configurations for 1x2x6 geometry
dune_fd_genie_nu_simple_flux_1x2x6_window: @local::dunefd_genie
dune_fd_genie_nu_simple_flux_1x2x6_window.FiducialCut:  "mbox:-380,-620,-20,380,620,1420"
dune_fd_genie_nu_simple_flux_1x2x6_window.FluxType: "simple_flux"
dune_fd_genie_nu_simple_flux_1x2x6_window.FluxSearchPaths: "/cvmfs/dune.osgstorage.org/pnfs/fnal.gov/usr/dune/persistent/stash/Flux/g4lbne/v3r5p4/QGSP_BERT/OptimizedEngineeredNov2017/neutrino/dune10kt_horizdrift_1x2x6/"
dune_fd_genie_nu_simple_flux_1x2x6_window.FluxFiles: ["gsimple_dune10kt_v1_1x2x6_g4lbne_v3r5p4_QGSP_BERT_OptimizedEngineeredNov2017_neutrino_*root"]
dune_fd_genie_nu_simple_flux_1x2x6_window.FluxCopyMethod: "DIRECT"
dune_fd_genie_nu_simple_flux_1x2x6_window.TopVolume: "volCryostat"

dune_fd_genie_nu_simple_flux_1x2x6_window_rotated_tpc: @local::dune_fd_genie_nu_simple_flux_1x2x6_window
dune_fd_genie_nu_simple_flux_1x2x6_window_rotated_tpc.FluxSearchPaths: "/cvmfs/dune.osgstorage.org/pnfs/fnal.gov/usr/dune/persistent/stash/TaskForce_Flux/opt_03/neutrino/gsimple/far/Rotate90Degrees/"
dune_fd_genie_nu_simple_flux_1x2x6_window_rotated_tpc.FluxFiles: ["gsimple_dune10kt_v1_1x2x6_Rotate90Degrees_g4lbne_v3r4p2_QGSP_BERT_CP_run15_12388_80GeV_neutrino_*root"]

dune_fd_genie_anu_simple_flux_1x2x6_window: @local::dune_fd_genie_nu_simple_flux_1x2x6_window
dune_fd_genie_anu_simple_flux_1x2x6_window.FluxSearchPaths: "/cvmfs/dune.osgstorage.org/pnfs/fnal.gov/usr/dune/persistent/stash/Flux/g4lbne/v3r5p4/QGSP_BERT/OptimizedEngineeredNov2017/antineutrino/dune10kt_horizdrift_1x2x6/"
dune_fd_genie_anu_simple_flux_1x2x6_window.FluxFiles: ["gsimple_dune10kt_v1_1x2x6_g4lbne_v3r5p4_QGSP_BERT_OptimizedEngineeredNov2017_antineutrino*root"]

dune_fd_genie_anu_simple_flux_1x2x6_window_rotated_tpc: @local::dune_fd_genie_nu_simple_flux_1x2x6_window
dune_fd_genie_anu_simple_flux_1x2x6_window_rotated_tpc.FluxSearchPaths: "/cvmfs/dune.osgstorage.org/pnfs/fnal.gov/usr/dune/persistent/stash/TaskForce_Flux/opt_03/antineutrino/gsimple/far/Rotate90Degrees/"
dune_fd_genie_anu_simple_flux_1x2x6_window_rotated_tpc.FluxFiles: ["gsimple_dune10kt_v1_1x2x6_Rotate90Degrees_g4lbne_v3r4p2_QGSP_BERT_CP_run15_12388_80GeV_antineutrino_*root"]

dunefdvd_1x6x6_genie_lbnf_dk2nu: {
    @table::dune_fd_genie_nu_simple_flux_1x2x6_window
    FiducialCut: "mbox: -350, -530, -20, 350, 530, 920"
    FluxType: "dk2nu"
    DetectorLocation: "dunevd10kt_v1_1x6x6"
    FluxSearchPaths: "/cvmfs/dune.osgstorage.org/pnfs/fnal.gov/usr/dune/persistent/stash/Flux/g4lbne/v3r5p4/QGSP_BERT/OptimizedEngineeredNov2017/neutrino/flux/"
    FluxFiles: ["g4lbne_v3r5p4_QGSP_BERT_OptimizedEngineeredNov2017_neutrino*.dk2nu.root"]
}

dunefdvd_1x8x14_genie_lbnf_dk2nu: {
    @table::dunefdvd_1x6x6_genie_lbnf_dk2nu
    FiducialCut: "mbox: -370, -750, -30, 370, 750, 2150"
    DetectorLocation: "dunevd10kt_v1_1x8x14"
}

dunefdvd_1x8x6_genie_lbnf_dk2nu: {
    @table::dunefdvd_1x6x6_genie_lbnf_dk2nu
    FiducialCut: "mbox: -370, -750, -30, 370, 750, 920"
    DetectorLocation: "dunevd10kt_v1_1x8x6"
}

dunefdvd_1x8x6_genie_anu_lbnf_dk2nu_fluxpaths: {
    FluxSearchPaths: "/cvmfs/dune.osgstorage.org/pnfs/fnal.gov/usr/dune/persistent/stash/Flux/g4lbne/v3r5p4/QGSP_BERT/OptimizedEngineeredNov2017/antineutrino/flux/"
    FluxFiles: ["g4lbne_v3r5p4_QGSP_BERT_OptimizedEngineeredNov2017_antineutrino*.dk2nu.root"]
}

#atmospheric neutrino configurations

# this configuration is tuned on the full geometry (10 kton) of the far detector:
dune_fd_genie_atmo: { 
  @table::standard_genie_atmo_flux    # inherit from the standard configuration
  
  BeamName:            "atm"
  GlobalTimeOffset:    0
  RandomTimeOffset:    0
  EventGeneratorList:  "Default"
  
  #
  # now fix the flux window parameters; from GENIE manual:
  #
  # "For a given direction, determined by the zenith angle and the azimuth angle,
  # the flux generation surface is a circular area, with radius RT,
  # which is tangent to a sphere of radius RL centered at the coordinate system origin."
  #
  # NOTE: in `GENIEHelper` (nutools `v2_26_00`) `BeamCenter` is not used
  # for the atmospheric neutrino flux (probably, because GENIE does not know
  # whet to do with it).
  #
  # Rl is how far from the origin the circular flux window is;
  # in GENIE style, this is meters; and it should be large enough to include
  # all the active volume of the detector (in principle, even more).
  # Now, this radius and center of the circle are pivoting around the origin
  # of the coordinate system, which in DUNE FD is on one of the faces
  # (the upstream face of the neutrino beam).
  # Unfortunately there are like 60 more meters of detector on one side of
  # that origin, which means we need a large radius Rl (more than 60 meters,
  # and in fact closer to the size of the diagonal, 5 (+) 5 (+) 60 m),
  # and the radius of the circle must be large enough that when the circle is
  # facing down (and since the origin is on a side of the detector,
  # half of the circle does not see any detector at all), that is 60 meters
  # of radius. Now, art:
  #
  #       _.     _________________________________________________
  #       /|    /!                                               /|
  #     15m    / !                                              / |
  #     /     /  !                                             /  |
  #   _L     /________________________________________________/   |       y
  #   ^      |   !                                            |   |      ^
  #   |      | * !                                            |   |      |     z
  #   |      |   !............................................|...|      .---->
  #  14m     |  ,                                             |  /      /
  #   |      | ,                                              | /      L x      
  #   v      |,                                               |/ 
  #   _      '________________________________________________/ 
  #         
  #          |<-- - - - - - - - - - 62 m - - - - - - - - - -->|
  #
  #       <========= LBNF is there; beam neutrinos go: -------->
  #
  # The origin of the DUNE coordinate system, which is the pivot GENIE rotates
  # the flux window around, is marked by (*).
  #
  Rl: 70 # meters
  # Rt is the radius of the circular flux window; again, should be large enough
  # that its "shadow" covers the detector
  Rt: 70 # meters
  # given that we picked a window specific to a volume, it makes sense to
  # specify the volume here as well (the dimensions are in the cartoon above):
  TopVolume:  "volCryostat"
  
  
  # now, about the orientation:
  # the flux is defined and interpreted by GENIE in a geographical reference
  # frame where z is orthogonal to earth surface pointing to the zenith,
  # x points toward south and y toward east (right-handed coordinate system).
  # DUNE geometry points z along the LBNF beam, and y in the plane containing z
  # and the vertical (see the axis triad in the cartoon above).
  # The flux needs to be rotated in that system, and GENIEHelper allows us to
  # express a rotation as in: what are the coordinates of the new axes,
  # expressed in the old reference frame.
  # Here we are rotating the geographical local frame ("old") into DUNE frame
  # ("new"). The rotation is complicate, and here we take heavy shortcuts:
  # * DUNE y axis should be the direction of the vertical (which makes the
  #   beam come with an angle); so this is effectively the geographic z axis
  # * DUNE z axis roughly points to Fermilab; since Fermilab is east-south-east
  #   of DUNE, the angle to Fermilab from the geographical x axis as defined
  #   above is positive and close to 90 degrees (in fact, roughly 81 degrees);
  #   but z verse is the opposite, from Fermilab (LBNF) to DUNE FD,
  #   effectively pointing WNW.
  # * DUNE x axis will have enough good taste to adapt and make the frame
  #   right-handed; this would mean pointing SSW.
  # My shortcuts so far: z becomes from ESE becomes just E, and# x axis goes
  # south. Defining phi as the angle from south, measured toward west
  # (counterclockwise), which is the geographical convention ("old"), this is
  # equivalent to Fermilab being at phi = pi/2 (again, z points backward).
  # 
  # To be tested: complete rotation is:
  # [
  #   +sin(phi), -cos(phi), +0.000,
  #     +0.000,    +0.000,  +1.000, # new y is old z
  #   -cos(phi), -sin(phi), +0.000  # new z
  # ]
  #
  # A rough estimation of the actual phi is +81 degrees, or +1.414 rad.
  # Finally, note that the Bartol flux the most standard producton picks
  # (from http://www-pnp.physics.ox.ac.uk/~barr/fluxfiles)
  # have no azimut parametrization, so phi does not really matter yet.
  # Honda flux, on the other hand, does have such parametrization.
  #
  FluxRotCfg: "newxyz"
  FluxRotValues: [
    +1.0, +0.0, +0.0,  # new x axis in old coordinates: be nice and fix things (ends up the same)
    +0.0, +0.0, +1.0,  # new y axis in old coordinates: vertical
    +0.0, -1.0, +0.0   # new z axis in old coordinates: westward
    ]
  # 
  # A more complete rotation is available as
  # `dune_fd_atmo_flux_rotation_precise` below.
  #
  # NOTE: these rotations will need to be removed or fixed if feature request
  # 21185 (https://cdcvs.fnal.gov/redmine/issues/21185) is implemented.
  #
  
} # dune_fd_genie_atmo


#
# the following snipper can be added to an existing configuration
# (including the ones with `dune_fd_genie_atmo` above) to test a
# more precise rotation which correctly aligns the detector to
# the geographic coordinates (see issue #20034).
# Usage: plug it at the end of the generator configuration:
#     
#     physics: {
#       producers: {
#         generator: {
#           ...
#           @table::dune_fd_atmo_flux_rotation_precise
#         }
#       }
#     }
#
# or overwrite the existing configuration like in:
#     
#     physics.producers.generator: {
#       @table::physics.producers.generator
#       @table::dune_fd_atmo_flux_rotation_precise
#     }
#     
# WARNING: this configuration was, at the time it was added
#          (October 29, 2018), UNTESTED.
#
dune_fd_atmo_flux_rotation_precise: {
  # the details of this rotation are explained in `dune_fd_genie_atmo`;
  # the imparted rotation is phi = 1.414 rad (81.016 degrees)
  # from southward direction counterclockwise
  # (and rounded to 4 significant digits).
  FluxRotCfg: "newxyz"
  FluxRotValues: [
    +0.9877, -0.1564, +0.0000,  # new x axis in old coordinates: be nice and fix things (ends up the same)
    +0.0000, +0.0000, +1.0000,  # new y axis in old coordinates: vertical
    -0.1564, -0.9877, +0.0000   # new z axis in old coordinates: away from Batavia, IL
    ]

} # dune_fd_atmo_flux_rotation_precise


dune_fd_genie_atmo_min_flux:            @local::dune_fd_genie_atmo
dune_fd_genie_atmo_min_flux.TopVolume:  "volCryostat"
dune_fd_genie_atmo_min_flux.FluxType:   "atmo_BARTOL"
dune_fd_genie_atmo_min_flux.FluxSearchPaths: "/cvmfs/dune.osgstorage.org/pnfs/fnal.gov/usr/dune/persistent/stash/TaskForce_Flux/atmos/Bartol/"
dune_fd_genie_atmo_min_flux.FluxFiles: ["fmin20_0401z.sou_nue","fmin20_0401z.sou_num","fmin20_0401z.sou_nbe","fmin20_0401z.sou_nbm"]
#dune_fd_genie_atmo_min_flux.FluxSearchPaths: "/cvmfs/dune.osgstorage.org/pnfs/fnal.gov/usr/dune/persistent/stash/TaskForce_Flux/atmos/Bartol/"
#dune_fd_genie_atmo_min_flux.FluxFiles: ["fmin10_0401z.sou_nue","fmin10_0401z.sou_num","fmin10_0401z.sou_nbe","fmin10_0401z.sou_nbm"]
dune_fd_genie_atmo_min_flux.FluxCopyMethod: "DIRECT"

dune_fd_genie_atmo_max_flux:            @local::dune_fd_genie_atmo
dune_fd_genie_atmo_max_flux.TopVolume:  "volCryostat"
dune_fd_genie_atmo_max_flux.FluxType:   "atmo_BARTOL"
dune_fd_genie_atmo_max_flux.FluxSearchPaths: "/cvmfs/dune.osgstorage.org/pnfs/fnal.gov/usr/dune/persistent/stash/TaskForce_Flux/atmos/Bartol/"
dune_fd_genie_atmo_max_flux.FluxFiles: ["fmax20_0401z.sou_nue","fmax20_0401z.sou_num","fmax20_0401z.sou_nbe","fmax20_0401z.sou_nbm"]
#dune_fd_genie_atmo_max_flux.FluxSearchPaths: "/cvmfs/dune.osgstorage.org/pnfs/fnal.gov/usr/dune/persistent/stash/TaskForce_Flux/atmos/Bartol/"
#dune_fd_genie_atmo_max_flux.FluxFiles: ["fmax10_0401z.sou_nue","fmax10_0401z.sou_num","fmax10_0401z.sou_nbe","fmax10_0401z.sou_nbm"]
dune_fd_genie_atmo_max_flux.FluxCopyMethod: "DIRECT"

dune_fd_genie_atmo_max_hondaflux:            @local::dune_fd_genie_atmo
dune_fd_genie_atmo_max_hondaflux.TopVolume:  "volCryostat"
dune_fd_genie_atmo_max_hondaflux.FluxType:   "atmo_HONDA"
dune_fd_genie_atmo_max_hondaflux.FluxSearchPaths: "/cvmfs/dune.osgstorage.org/pnfs/fnal.gov/usr/dune/persistent/stash/TaskForce_Flux/atmos/Honda/"
dune_fd_genie_atmo_max_hondaflux.FluxFiles: ["hms-ally-20-12-solmax.d","hms-ally-20-12-solmax.d","hms-ally-20-12-solmax.d","hms-ally-20-12-solmax.d"]
dune_fd_genie_atmo_max_hondaflux.FluxCopyMethod: "DIRECT"

dune_fd_genie_atmo_min_hondaflux: @local::dune_fd_genie_atmo_max_hondaflux
dune_fd_genie_atmo_min_hondaflux.FluxFiles: ["hms-ally-20-12-solmin.d","hms-ally-20-12-solmin.d","hms-ally-20-12-solmin.d","hms-ally-20-12-solmin.d"]

END_PROLOG