Generators

From Mu2eWiki
Revision as of 18:28, 8 May 2017 by Rlc (talk | contribs)
Jump to navigation Jump to search

Construction.jpeg This page is a draft, please help complete it!

Introduction

Generators are art modules which produce the first particle in a simulation job. Some may be configured by fcl only, but some can be configured by specifying a SimpleConfig file.

Protons

This generator creates a proton approaching the proton target. You can set the width of the beam, and offset it for systematic studies.

Since each generated event is a proton on target, normalization can be done by simple counting events, equal to the number of protons on target.

Examples of configurations are under the JobConfig directory of Offline.

It is difficult to generate a billionth of the full luminosity.

Cosmics

This generator creates a cosmic ray approaching the detector. If the generator is the "Daya Bay" generator, then only one particle, a muon, is generated per event. The CRY generator allows multiple correlated particles per event, of several species, simulating an atmospheric shower. Example configurations are under the JobConfig directory of Offline.

The cosmic ray phase space is specified by a horizontal or vertical plane with limited extent. There also may be angle and momentum restrictions.

The Daya Bay generator will print an integral at the start of the job:

Total cosmic rate = 93510.8 Hz

This tells you how much time you have simulated: Ngen/rate=time.

It is not too hard to generate a comics sample corresponding to the full run of mu2e.

Stopped Particles

This generator starts with a list of stopped particle positions, usually in a ntuple made in a earlier job. Each event starts with reading in a stopped particle position. The stopped particle is decayed. Usually the decay mode and daughters are specified. The decay model can be specified. A common choice would be to generate DIO decays of stopped muons. Another common request is to produce a given spectrum of soft photons, or other particles, as would be expected from stopped muon capture. For some studies stopped pions may be used on input.

Here is a an example stanza to configure for DIO decays.

physics: {
   producers: {
      generate: {
         module_type: "StoppedParticleReactionGun"
         muonStops: {
            averageNumRecordsToUse: 500000
            branchName: "stops"
            inputFiles: ["mergedMuonStops/nts.mu2e.cd3-beam-nts3-tgtstops.v533_v543_v563b.001002_00000000.root"]
            treeName: "stoppedMuonDumper/stops"
         }
         physics: {
            elow: 1
            pdgId: 11
            spectrumResolution: 1e-1
            spectrumShape: "Czarnecki"
            spectrumVariable: "totalEnergy"
         }
      }
   }
}

Normalization would have to be based on the jobs that produced the stopped particle list.

Particle guns

Particle guns use the EventGenerator art module and are configured with SimpleConfig. Examples are in the Mu2eG4/test directory of Offline.

GenEventCounter