Skip to content

Commit

Permalink
added new example 08 (#4)
Browse files Browse the repository at this point in the history
  • Loading branch information
lobis committed Apr 27, 2021
1 parent b53be50 commit 8c3a9cf
Show file tree
Hide file tree
Showing 5 changed files with 1,773 additions and 0 deletions.
84 changes: 84 additions & 0 deletions examples/08.RestGenerator/Pb210FromGenerator.rml
Original file line number Diff line number Diff line change
@@ -0,0 +1,84 @@
<?xml version="4.0" encoding="UTF-8" standalone="no" ?>

<restG4>
<globals>
<parameter name="mainDataPath" value=""/>
<parameter name="verboseLevel" value="essential"/>
</globals>

<TRestRun name="Run metadata" title="REST Metadata run info (template)">
<parameter name="experimentName" value="BabyIAXO"/>
<parameter name="runType" value="simulation"/>

// Value "auto" can add simulation run into database ???
<parameter name="runNumber" value="1"/>

// If no runTag is given in a restG4 simulation. The title of TRestGeant4Metadata will be used
<parameter name="runTag" value="BabyIAXO"/>

// The output filename can be built using any member from TRestRun, fRunNumber, fRunTag, ..., etc.
<parameter name="outputFile" value="Run[fRunNumber]_[fRunTag]_[fExperimentName].root"/>

<parameter name="runDescription" value=""/>
<parameter name="user" value="${USER}"/>
<parameter name="verboseLevel" value="3"/>
<parameter name="overwrite" value="off"/>
<parameter name="readOnly" value="false"/>
</TRestRun>

<TRestGeant4Metadata name="Pb210FromGenerator" title="Pb210FromGenerator">

<parameter name="gdml_file" value="geometry/BabyIAXO.gdml"/>
<parameter name="subEventTimeDelay" value="100" units="us"/>

<parameter name="seed" value="17021981"/>

<parameter name="Nevents" value="10"/>

<generator type="volume" from="gas">

<source use="Pb210FromShielding.root"/>

</generator>

<storage sensitiveVolume="gas">
<parameter name="energyRange" value="(0,5)" units="MeV"/>
<activeVolume name="gas" chance="1" maxStepSize="1mm"/>
</storage>
</TRestGeant4Metadata>

<TRestGeant4PhysicsLists name="default" title="First physics list implementation.">

<parameter name="cutForGamma" value="1" units="mm"/>
<parameter name="cutForElectron" value="2" units="mm"/>
<parameter name="cutForPositron" value="1" units="mm"/>

<parameter name="cutForMuon" value="1" units="mm"/>
<parameter name="cutForNeutron" value="1" units="mm"/>
<parameter name="minEnergyRangeProductionCuts" value="1" units="keV"/>
<parameter name="maxEnergyRangeProductionCuts" value="1" units="GeV"/>

// Use only one EM physics list
<!-- EM Physics lists -->
<physicsList name="G4EmLivermorePhysics"></physicsList>
<!-- <physicsList name="G4EmPenelopePhysics"> </physicsList> -->
<!-- <physicsList name="G4EmStandardPhysics_option3"> </physicsList> -->

<!-- Decay physics lists -->
<physicsList name="G4DecayPhysics"></physicsList>
<physicsList name="G4RadioactiveDecayPhysics"></physicsList>
<physicsList name="G4RadioactiveDecay">
<option name="ICM" value="true"/>
<option name="ARM" value="true"/>
</physicsList>

<!-- Hadron physics lists -->
<physicsList name="G4HadronElasticPhysicsHP"></physicsList>
<physicsList name="G4IonBinaryCascadePhysics"></physicsList>
<physicsList name="G4HadronPhysicsQGSP_BIC_HP"></physicsList>
<physicsList name="G4NeutronTrackingCut"></physicsList>
<physicsList name="G4EmExtraPhysics"></physicsList>

</TRestGeant4PhysicsLists>

</restG4>
Binary file added examples/08.RestGenerator/Pb210FromShielding.root
Binary file not shown.
18 changes: 18 additions & 0 deletions examples/08.RestGenerator/README.md
Original file line number Diff line number Diff line change
@@ -0,0 +1,18 @@
**Date:** 27/April/2021

**Author:** Luis Antonio Obis Aparicio ([email protected])

This example reads source data from a `root` file (`Pb210FromShielding.root`) generated via `restGenerator`.
The file in this example comes from the background simulation for the contamination of Pb210 in the shielding of BabyIAXO.

#### Detector geometry

BabyIAXO geometry

#### Event generator

Reading from the `root` file, particles will be generated on a virtual cylinder surrounding the chamber (sensitive volume).

#### Testing the example

TODO.
Loading

0 comments on commit 8c3a9cf

Please sign in to comment.