Cosmic Simulations

From Mu2eWiki
Revision as of 13:00, 15 December 2023 by Sophie (talk | contribs)
Jump to navigation Jump to search

This page gives the details which are specific to simulating cosmic events.

The standard workflow for a given cosmic campaign contains two stages:

  • Stage 1 (S1) - generates cosmics in Mu2e volume using given generator, outputs "sim" .art file
  • Stage 2 (S2) - takes S1 as input and resamples a given number of times to create final sample of "dts" .art files

When generating a given sample the user will probably want a given number of hours of "livetime." This will be the output of S2. The user must configure the jobs with the end result in mind.

This page details the quantities which must be considered.

Livetime

Livetime of a given Cosmic sample can be established with the following command:

mu2e -c Offline/Print/fcl/printCosmicLivetime.fcl /pnfs/mu2e/tape/phy-sim/sim/mu2e/CosmicDSStopsCORSIKA/MDC2020n_10h/art/cb/f0/sim.mu2e.CosmicDSStopsCORSIKA.MDC2020n_10h.103.art

which will print out a list of quantities associated with the provided data-set. These include N Primaries, Area, LowE (lowest energy), HighE (highest energy), FluxConstant and Livetime. For the standard configuration, a single job corresponds to 0.062266 seconds of livetime.

Resampling at S2

The standard assumption is that S2 resamples S1 by a factor of 1000. This means to simulate a total of X hours we need X/1000 hours at S1.

Standard 10hour Configuration

Our standard jobs take "10 hours" of livetime. We configure our 0.062 s jobs in the following way:

  • 500 jobs at Stage 1, this corresponds to 500*0.0622 = 31 s
  • Resampling factor of 1000 = 31*1000 s = 8.6 hours