-
Notifications
You must be signed in to change notification settings - Fork 6
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Merge pull request #133 from rest-for-physics/cosmics-benchmark
Comparison of cosmic generator with disk
- Loading branch information
Showing
15 changed files
with
588 additions
and
14 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,4 @@ | ||
|
|
||
| This is a comparison of the new Cosmic Generator vs the old approach of a disk above the detector. | ||
|
|
||
| In the old approach most primaries (depending on the size of the disk) never reach the detector, leading to a huge inefficiency. |
69 changes: 69 additions & 0 deletions
69
examples/12.Generators/CosmicGeneratorBenchmarks/cosmicGeneratorIAXO.rml
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,69 @@ | ||
| <?xml version="1.0" encoding="UTF-8" standalone="no" ?> | ||
|
|
||
| <restG4> | ||
| <globals> | ||
| <parameter name="mainDataPath" value=""/> | ||
| <parameter name="verboseLevel" value="essential"/> | ||
| </globals> | ||
|
|
||
| <TRestRun> | ||
| <parameter name="experimentName" value="CosmicGenerator"/> | ||
| <parameter name="runType" value="simulation"/> | ||
| <parameter name="runNumber" value="1"/> | ||
| <parameter name="runTag" value="CosmicGenerator"/> | ||
| <parameter name="outputFileName" value="CosmicGenerator.root"/> | ||
| <parameter name="runDescription" value=""/> | ||
| <parameter name="user" value="${USER}"/> | ||
| <parameter name="overwrite" value="off"/> | ||
| <parameter name="readOnly" value="false"/> | ||
| </TRestRun> | ||
|
|
||
| <TRestGeant4Metadata> | ||
| <parameter name="gdmlFile" value="../../13.IAXO/geometry/setup.gdml"/> | ||
| <parameter name="seed" value="137"/> | ||
| <parameter name="nRequestedEntries" value="1000000"/> | ||
| <parameter name="saveAllEvents" value="false"/> | ||
|
|
||
| <generator type="cosmic"> | ||
| <source use="cosmics" direction="(0,-1,0)" filename="${REST_PATH}/data/distributions/CosmicsCry.root" | ||
| particles="muon_minus"/> | ||
| </generator> | ||
|
|
||
| <detector> | ||
| <volume name="gasVolume" sensitive="true" maxStepSize="1mm"/> | ||
| </detector> | ||
|
|
||
| </TRestGeant4Metadata> | ||
|
|
||
| <TRestGeant4PhysicsLists> | ||
| <parameter name="cutForGamma" value="1" units="mm"/> | ||
| <parameter name="cutForElectron" value="1" 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"/> | ||
|
|
||
| <physicsList name="G4EmLivermorePhysics"> | ||
| <option name="pixe" value="false"/> | ||
| </physicsList> | ||
|
|
||
| <physicsList name="G4HadronElasticPhysicsHP"/> | ||
| <physicsList name="G4IonBinaryCascadePhysics"/> | ||
| <physicsList name="G4HadronPhysicsQGSP_BIC_HP"/> | ||
| <physicsList name="G4EmExtraPhysics"/> | ||
|
|
||
| <physicsList name="G4DecayPhysics"/> | ||
| <physicsList name="G4RadioactiveDecayPhysics"/> | ||
| <physicsList name="G4RadioactiveDecay"> | ||
| <option name="ICM" value="true"/> | ||
| <option name="ARM" value="true"/> | ||
| </physicsList> | ||
|
|
||
| <!-- | ||
| <physicsList name="G4NeutronTrackingCut"/> | ||
| --> | ||
|
|
||
| </TRestGeant4PhysicsLists> | ||
| </restG4> |
69 changes: 69 additions & 0 deletions
69
examples/12.Generators/CosmicGeneratorBenchmarks/cosmicGeneratorSimple.rml
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,69 @@ | ||
| <?xml version="1.0" encoding="UTF-8" standalone="no" ?> | ||
|
|
||
| <restG4> | ||
| <globals> | ||
| <parameter name="mainDataPath" value=""/> | ||
| <parameter name="verboseLevel" value="essential"/> | ||
| </globals> | ||
|
|
||
| <TRestRun> | ||
| <parameter name="experimentName" value="CosmicGenerator"/> | ||
| <parameter name="runType" value="simulation"/> | ||
| <parameter name="runNumber" value="1"/> | ||
| <parameter name="runTag" value="CosmicGenerator"/> | ||
| <parameter name="outputFileName" value="CosmicGenerator.root"/> | ||
| <parameter name="runDescription" value=""/> | ||
| <parameter name="user" value="${USER}"/> | ||
| <parameter name="overwrite" value="off"/> | ||
| <parameter name="readOnly" value="false"/> | ||
| </TRestRun> | ||
|
|
||
| <TRestGeant4Metadata> | ||
| <parameter name="gdmlFile" value="../geometry/setup.gdml"/> | ||
| <parameter name="seed" value="137"/> | ||
| <parameter name="nRequestedEntries" value="1000000"/> | ||
| <parameter name="saveAllEvents" value="false"/> | ||
|
|
||
| <generator type="cosmic"> | ||
| <source use="cosmics" direction="(0,-1,0)" filename="${REST_PATH}/data/distributions/CosmicsCry.root" | ||
| particles="muon_minus"/> | ||
| </generator> | ||
|
|
||
| <detector> | ||
| <volume name="gasVolume" sensitive="true" maxStepSize="1mm"/> | ||
| </detector> | ||
|
|
||
| </TRestGeant4Metadata> | ||
|
|
||
| <TRestGeant4PhysicsLists> | ||
| <parameter name="cutForGamma" value="1" units="mm"/> | ||
| <parameter name="cutForElectron" value="1" 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"/> | ||
|
|
||
| <physicsList name="G4EmLivermorePhysics"> | ||
| <option name="pixe" value="false"/> | ||
| </physicsList> | ||
|
|
||
| <physicsList name="G4HadronElasticPhysicsHP"/> | ||
| <physicsList name="G4IonBinaryCascadePhysics"/> | ||
| <physicsList name="G4HadronPhysicsQGSP_BIC_HP"/> | ||
| <physicsList name="G4EmExtraPhysics"/> | ||
|
|
||
| <physicsList name="G4DecayPhysics"/> | ||
| <physicsList name="G4RadioactiveDecayPhysics"/> | ||
| <physicsList name="G4RadioactiveDecay"> | ||
| <option name="ICM" value="true"/> | ||
| <option name="ARM" value="true"/> | ||
| </physicsList> | ||
|
|
||
| <!-- | ||
| <physicsList name="G4NeutronTrackingCut"/> | ||
| --> | ||
|
|
||
| </TRestGeant4PhysicsLists> | ||
| </restG4> |
69 changes: 69 additions & 0 deletions
69
examples/12.Generators/CosmicGeneratorBenchmarks/cosmicGeneratorSphere.rml
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,69 @@ | ||
| <?xml version="1.0" encoding="UTF-8" standalone="no" ?> | ||
|
|
||
| <restG4> | ||
| <globals> | ||
| <parameter name="mainDataPath" value=""/> | ||
| <parameter name="verboseLevel" value="essential"/> | ||
| </globals> | ||
|
|
||
| <TRestRun> | ||
| <parameter name="experimentName" value="CosmicGenerator"/> | ||
| <parameter name="runType" value="simulation"/> | ||
| <parameter name="runNumber" value="1"/> | ||
| <parameter name="runTag" value="CosmicGenerator"/> | ||
| <parameter name="outputFileName" value="CosmicGenerator.root"/> | ||
| <parameter name="runDescription" value=""/> | ||
| <parameter name="user" value="${USER}"/> | ||
| <parameter name="overwrite" value="off"/> | ||
| <parameter name="readOnly" value="false"/> | ||
| </TRestRun> | ||
|
|
||
| <TRestGeant4Metadata> | ||
| <parameter name="gdmlFile" value="../geometry-sphere/setup.gdml"/> | ||
| <parameter name="seed" value="137"/> | ||
| <parameter name="nRequestedEntries" value="1000000"/> | ||
| <parameter name="saveAllEvents" value="false"/> | ||
|
|
||
| <generator type="cosmic"> | ||
| <source use="cosmics" direction="(0,-1,0)" filename="${REST_PATH}/data/distributions/CosmicsCry.root" | ||
| particles="muon_minus"/> | ||
| </generator> | ||
|
|
||
| <detector> | ||
| <volume name="gasVolume" sensitive="true" maxStepSize="1mm"/> | ||
| </detector> | ||
|
|
||
| </TRestGeant4Metadata> | ||
|
|
||
| <TRestGeant4PhysicsLists> | ||
| <parameter name="cutForGamma" value="1" units="mm"/> | ||
| <parameter name="cutForElectron" value="1" 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"/> | ||
|
|
||
| <physicsList name="G4EmLivermorePhysics"> | ||
| <option name="pixe" value="false"/> | ||
| </physicsList> | ||
|
|
||
| <physicsList name="G4HadronElasticPhysicsHP"/> | ||
| <physicsList name="G4IonBinaryCascadePhysics"/> | ||
| <physicsList name="G4HadronPhysicsQGSP_BIC_HP"/> | ||
| <physicsList name="G4EmExtraPhysics"/> | ||
|
|
||
| <physicsList name="G4DecayPhysics"/> | ||
| <physicsList name="G4RadioactiveDecayPhysics"/> | ||
| <physicsList name="G4RadioactiveDecay"> | ||
| <option name="ICM" value="true"/> | ||
| <option name="ARM" value="true"/> | ||
| </physicsList> | ||
|
|
||
| <!-- | ||
| <physicsList name="G4NeutronTrackingCut"/> | ||
| --> | ||
|
|
||
| </TRestGeant4PhysicsLists> | ||
| </restG4> |
75 changes: 75 additions & 0 deletions
75
examples/12.Generators/CosmicGeneratorBenchmarks/diskIAXO.rml
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,75 @@ | ||
| <?xml version="1.0" encoding="UTF-8" standalone="no" ?> | ||
|
|
||
| <restG4> | ||
| <globals> | ||
| <parameter name="mainDataPath" value=""/> | ||
| <parameter name="verboseLevel" value="essential"/> | ||
| </globals> | ||
|
|
||
| <TRestRun> | ||
| <parameter name="experimentName" value="CosmicGenerator"/> | ||
| <parameter name="runType" value="simulation"/> | ||
| <parameter name="runNumber" value="1"/> | ||
| <parameter name="runTag" value="CosmicGenerator"/> | ||
| <parameter name="outputFileName" value="CosmicGenerator.root"/> | ||
| <parameter name="runDescription" value=""/> | ||
| <parameter name="user" value="${USER}"/> | ||
| <parameter name="overwrite" value="off"/> | ||
| <parameter name="readOnly" value="false"/> | ||
| </TRestRun> | ||
|
|
||
| <TRestGeant4Metadata> | ||
| <parameter name="gdmlFile" value="../../13.IAXO/geometry/setup.gdml"/> | ||
| <parameter name="seed" value="137"/> | ||
| <parameter name="nRequestedEntries" value="1000000"/> | ||
| <parameter name="saveAllEvents" value="false"/> | ||
|
|
||
| <generator type="surface" shape="circle" | ||
| position="(0,177.79904,0)cm" size="(1777.9904,0,0)cm" | ||
| rotationAngle="90deg" rotationAxis="(1,0,0)"> | ||
| <source particle="mu-"> | ||
| <energy type="TH2D" name="muon_minus_energy_zenith" | ||
| file="${REST_PATH}/data/distributions/CosmicsCry.root"/> | ||
| <angular type="TH2D" name="muon_minus_energy_zenith" direction="(0,-1,0)" | ||
| file="${REST_PATH}/data/distributions/CosmicsCry.root" /> | ||
| </source> | ||
| </generator> | ||
|
|
||
| <detector> | ||
| <volume name="gasVolume" sensitive="true" maxStepSize="1mm"/> | ||
| </detector> | ||
|
|
||
| </TRestGeant4Metadata> | ||
|
|
||
| <TRestGeant4PhysicsLists> | ||
| <parameter name="cutForGamma" value="1" units="mm"/> | ||
| <parameter name="cutForElectron" value="1" 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"/> | ||
|
|
||
| <physicsList name="G4EmLivermorePhysics"> | ||
| <option name="pixe" value="false"/> | ||
| </physicsList> | ||
|
|
||
| <physicsList name="G4HadronElasticPhysicsHP"/> | ||
| <physicsList name="G4IonBinaryCascadePhysics"/> | ||
| <physicsList name="G4HadronPhysicsQGSP_BIC_HP"/> | ||
| <physicsList name="G4EmExtraPhysics"/> | ||
|
|
||
| <physicsList name="G4DecayPhysics"/> | ||
| <physicsList name="G4RadioactiveDecayPhysics"/> | ||
| <physicsList name="G4RadioactiveDecay"> | ||
| <option name="ICM" value="true"/> | ||
| <option name="ARM" value="true"/> | ||
| </physicsList> | ||
|
|
||
| <!-- | ||
| <physicsList name="G4NeutronTrackingCut"/> | ||
| --> | ||
|
|
||
| </TRestGeant4PhysicsLists> | ||
| </restG4> |
76 changes: 76 additions & 0 deletions
76
examples/12.Generators/CosmicGeneratorBenchmarks/diskSimple.rml
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,76 @@ | ||
| <?xml version="1.0" encoding="UTF-8" standalone="no" ?> | ||
|
|
||
| <restG4> | ||
| <globals> | ||
| <parameter name="mainDataPath" value=""/> | ||
| <parameter name="verboseLevel" value="essential"/> | ||
| </globals> | ||
|
|
||
| <TRestRun> | ||
| <parameter name="experimentName" value="CosmicGenerator"/> | ||
| <parameter name="runType" value="simulation"/> | ||
| <parameter name="runNumber" value="1"/> | ||
| <parameter name="runTag" value="CosmicGenerator"/> | ||
| <parameter name="outputFileName" value="CosmicGenerator.root"/> | ||
| <parameter name="runDescription" value=""/> | ||
| <parameter name="user" value="${USER}"/> | ||
| <parameter name="overwrite" value="off"/> | ||
| <parameter name="readOnly" value="false"/> | ||
| </TRestRun> | ||
|
|
||
| <TRestGeant4Metadata> | ||
| <parameter name="gdmlFile" value="../geometry/setup.gdml"/> | ||
| <parameter name="seed" value="137"/> | ||
| <parameter name="nRequestedEntries" value="1000000"/> | ||
| <parameter name="saveAllEvents" value="false"/> | ||
|
|
||
| <!-- Use r=radius of the cosmic generator --> | ||
| <generator type="surface" shape="circle" | ||
| position="(0,173.20508,0)mm" size="(1732.0508,0,0)mm" | ||
| rotationAngle="90deg" rotationAxis="(1,0,0)"> | ||
| <source particle="mu-"> | ||
| <energy type="TH2D" name="muon_minus_energy_zenith" | ||
| file="${REST_PATH}/data/distributions/CosmicsCry.root"/> | ||
| <angular type="TH2D" name="muon_minus_energy_zenith" direction="(0,-1,0)" | ||
| file="${REST_PATH}/data/distributions/CosmicsCry.root" /> | ||
| </source> | ||
| </generator> | ||
|
|
||
| <detector> | ||
| <volume name="gasVolume" sensitive="true" maxStepSize="1mm"/> | ||
| </detector> | ||
|
|
||
| </TRestGeant4Metadata> | ||
|
|
||
| <TRestGeant4PhysicsLists> | ||
| <parameter name="cutForGamma" value="1" units="mm"/> | ||
| <parameter name="cutForElectron" value="1" 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"/> | ||
|
|
||
| <physicsList name="G4EmLivermorePhysics"> | ||
| <option name="pixe" value="false"/> | ||
| </physicsList> | ||
|
|
||
| <physicsList name="G4HadronElasticPhysicsHP"/> | ||
| <physicsList name="G4IonBinaryCascadePhysics"/> | ||
| <physicsList name="G4HadronPhysicsQGSP_BIC_HP"/> | ||
| <physicsList name="G4EmExtraPhysics"/> | ||
|
|
||
| <physicsList name="G4DecayPhysics"/> | ||
| <physicsList name="G4RadioactiveDecayPhysics"/> | ||
| <physicsList name="G4RadioactiveDecay"> | ||
| <option name="ICM" value="true"/> | ||
| <option name="ARM" value="true"/> | ||
| </physicsList> | ||
|
|
||
| <!-- | ||
| <physicsList name="G4NeutronTrackingCut"/> | ||
| --> | ||
|
|
||
| </TRestGeant4PhysicsLists> | ||
| </restG4> |
Oops, something went wrong.