- Table of contents
- Dossier Data Sets, Universes and FCL File Generation
Dossier Data Sets, Universes and FCL File Generation¶
Dossier Lists¶
Combinatorics¶
Between these two experiments there's a rich set of data, but it isn't that uniform in coverage:
Expt | HARP, ITEP |
beam particle ("probe") | piminus, piplus, proton |
beam momentum | 1.4, 2, 3, 4, 5, 6, 6.25, 6.5, 7, 7.5, 8, 8.25, 8.5, 8.9, 9, 12 |
target | Al, Be, C, Cu, Fe, Nb, Pb, Sn, Ta, Ti, U |
Not all possible combinations are represented in the data.
In total there are 171 "experiment setups" in the *_dossier.fcl
files (146 in the _pass1
files that were used), ala: ITEP_proton_on_U_at_8p5GeV
which represent a collection of data. For such a collection there might be multiple histograms representing the energy distribution of different secondaries in different angular ranges.
For the same choice of beam particle, energy, and target ("exptsetup"
, sans experiment name), which is what the Geant4 stage must cover, there are 156 (137 in this pass) combinations. For the cases where both HARP and ITEP-771 ran the same "exptsetup"
one needs to run both analysis modules on the same Geant4 resulting particles; for cases where only one experiment ran that configuration one must run only the appropriate analyzer module.
Dossier information FCL files¶
There are separate scripts for creating FCL include files (attached below) which map Dossier records to which experimental setups for HARP
and ITEP-771
. These are generated and used effectively as "include" files in the art
module configurations.
these need documenting once the Dossier infrastructure for extracting records by experiment (or INSPIRE id) is cleaned up
HARP_dossier.fcl | complete set of HARP experiment setups |
ITEP_dossier.fcl | complete set of ITEP-771 experiment setups |
HARP_dossier_pass1.fcl | used in this initial study (missing entries) |
ITEP_dossier_pass1.fcl | used in this initial study (missing entries) |
These files have contents that look like:
BEGIN_PROLOG
HARP_piminus_on_Be_at_3GeV : [
9212, 9213, 9214, 9215, 9216,
9217, 9218, 9219, 9220, 9221,
9222, 9223, 9224, # piminus
9160, 9161, 9162, 9163, 9164,
9165, 9166, 9167, 9168, 9169,
9170, 9171, 9172 # piplus
]
HARP_piminus_on_Be_at_5GeV : [
9225, 9226, 9227, 9228, 9229,
9230, 9231, 9232, 9233, 9234,
9235, 9236, 9237, # piminus
9173, 9174, 9175, 9176, 9177,
9178, 9179, 9180, 9181, 9182,
9183, 9184, 9185 # piplus
]
...
ITEP_proton_on_U_at_8p5GeV : [
2632, 2633, 2630, 2631 # proton
]
END_PROLOG
These should list dossier entries clustered by produced secondary and in angular bin order, but that's only for convenience.
Random Universes¶
To explore the impact of model parameters we run Geant4 with those parameter modified from the default. Each combination of model parameters constitutes a "universe".
We need for each "universe" to represent the same choice of model parameters independent of the experimental setup.
The same set of "universe" configurations must be used when generating each of the experimental setups for the fitting procedure to make sense. Thus generating this file is an infrequent event (generally) signaling a new set of model parameters to explore.
To do this we construct another FCL file with parameter choices which we'll call: multiverse170208_Bertini.fcl
root -b -q generate_universes.C\(\"multiverse170208\",999\)
# or more generally
# root -b -q generate_universes.C\(\"basename\",999,\"hadronmodel\",true,iseed\)
# basename to make things unique
# 999 = # random universes (always includes the default unmodified model as universe 0 )
# if negative then generate N steps for each parameter over allowed range (one at a time)
# hadron model name ["Bertini"]
# is process level (vs. physics list (which isn't yet supported))
# iseed - random number seed [0]
The 170208
represents the day it was generated (2017-02-08). The script internally is set to explore 4 parameters of the Bertini model.
BEGIN_PROLOG
# config for default Bertini universe
BertiniDefault : {
module_type: ProcLevelMPVaryProducer
Verbosity: 0
HadronicModel: {
DefaultPhysics: true
ModelParameters: {
} # end-of-ModelParameters
} # end-of-HadronicModel
} # end-of-BertiniDefault
# parameters for block of Bertini universes
# Bertini RadiusScale 1 3.5 flat
# Bertini XSecScale 0.1 3 flat
# Bertini FermiScale 0.5 1 flat
# Bertini TrailingRadius 0 5 flat
BertiniRandom4Univ0001 : {
module_type: ProcLevelMPVaryProducer
Verbosity: 0
HadronicModel: {
DefaultPhysics: false
ModelParameters: {
RadiusScale : 2.54237
XSecScale : 1.90091
FermiScale : 0.836028
TrailingRadius : 3.93846
} # end-of-ModelParameters
} # end-of-HadronicModel
} # end-of-BertiniRandom4Univ0001
...
BertiniRandom4Univ0999 : {
module_type: ProcLevelMPVaryProducer
Verbosity: 0
HadronicModel: {
DefaultPhysics: false
ModelParameters: {
RadiusScale : 3.26108
XSecScale : 0.170275
FermiScale : 0.848354
TrailingRadius : 4.7255
} # end-of-ModelParameters
} # end-of-HadronicModel
} # end-of-BertiniRandom4Univ0999
END_PROLOG