Project

General

Profile

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