Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Standalone geometry to compare energy deposition at different depths for various materials #578

Open
wants to merge 23 commits into
base: develop
Choose a base branch
from
Open

Standalone geometry to compare energy deposition at different depths for various materials #578

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
23 commits
Select commit Hold shift + click to select a range
2b43762
First commit with simple geometry
rahmans1 Feb 17, 2023
9d44fd4
Remove dependence on main geometry folder
rahmans1 Feb 17, 2023
9342e0d
Proxy parallel world volume
rahmans1 Feb 17, 2023
33013d7
Remove dependence on main geometry
rahmans1 Feb 17, 2023
029ba84
Script to generate geometry and run macro for alternate material and …
rahmans1 Feb 17, 2023
a053233
Job submission script for beluga
rahmans1 Feb 17, 2023
88d0d16
Update script so that root file name has more information
rahmans1 Feb 17, 2023
a0f88d1
Uniquely identify each file
rahmans1 Feb 17, 2023
7947741
Generalize script to produce layered spectrometer straight sections
rahmans1 Mar 5, 2023
d6f3647
Example downstream coil config
rahmans1 Mar 5, 2023
293ebd5
Set root file name based on config name provided by the config file
rahmans1 Mar 5, 2023
9f630cc
Add config name to configuration file
rahmans1 Mar 5, 2023
488f694
Adapt job submission script to use configuration file instead of indi…
rahmans1 Mar 5, 2023
b1ccdc3
Fix typo
rahmans1 Mar 5, 2023
6b6dd43
Example upstream configuration file
rahmans1 Mar 5, 2023
9d3d786
Needs relative path from remoll main folder
rahmans1 Mar 5, 2023
34856f5
First iteration of instructions
rahmans1 Mar 5, 2023
8d71ebc
Include information about detector ids
rahmans1 Mar 5, 2023
02495ff
Add material definition of G10 mixture used in full sim
rahmans1 Mar 6, 2023
11eb118
Adapt code to latest submission config
rahmans1 Mar 6, 2023
3d83d3c
Optimal resource requirements on beluga
rahmans1 Mar 6, 2023
c6bd3f2
Fix typo
rahmans1 Mar 6, 2023
581bd40
Merge branch 'JeffersonLab:develop' into 577-feature-standAloneEnergy…
rahmans1 Mar 8, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
59 changes: 59 additions & 0 deletions geometry_standalone_edep/mollerMother.gdml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,59 @@
<?xml version="1.0" encoding="UTF-8" standalone="no" ?>

<gdml xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://service-spi.web.cern.ch/service-spi/app/releases/GDML/schema/gdml.xsd">

<define>
<constant name="thickness" value="4000"/>
</define>

<materials>
<material name="Epoxy" state="solid">
<D value="1.3" unit="g/cm3"/>
<fraction n="0.5354" ref="C"/>
<fraction n="0.1318" ref="H"/>
<fraction n="0.3328" ref="O"/>
</material>
<material name="G10" state="solid">
<D value="1.3" unit="g/cm3"/>
<fraction n="0.773" ref="G4_SILICON_DIOXIDE"/>
<fraction n="0.147" ref="Epoxy"/>
<fraction n="0.080" ref="G4_Cl"/>
</material>
</materials>

<solids>
<box lunit="mm" name="world_solid" x="1000" y="1000" z="9000"/>

<box lunit="mm" name="testMaterial" x="1000" y="1000" z="thickness"/>

</solids>

<structure>

<volume name="layer1_log">
<materialref ref="G4_Bi"/>
<solidref ref="testMaterial"/>
<auxiliary auxtype="SensDet" auxvalue="planeDet"/>
<auxiliary auxtype="Color" auxvalue="Red"/>
<auxiliary auxtype="DetNo" auxvalue="10"/>
</volume>

<volume name="logicMother">
<materialref ref="G4_AIR"/>
<solidref ref="world_solid"/>

<physvol name="layer1_phys">
<volumeref ref="layer1_log"/>
<position name="layer1_pos" unit="mm" z="thickness/2."/>
</physvol>

<auxiliary auxtype="Alpha" auxvalue="0.1"/>
</volume>

</structure>

<setup name="Default" version="1.0">
<world ref="logicMother"/>
</setup>

</gdml>
75 changes: 75 additions & 0 deletions geometry_standalone_edep/mollerParallel.gdml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,75 @@
<?xml version="1.0" encoding="UTF-8" standalone="no" ?>

<gdml xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://service-spi.web.cern.ch/service-spi/app/releases/GDML/schema/gdml.xsd">

<define>
</define>

<materials>
</materials>

<solids>

<!-- Parallel world must be smaller than the physical world,
private communcation with Makoto Asai -->
<box lunit="mm" name="parallel_solid" x="999" y="999" z="8999"/>


<box
name="vertVirtualPlane_solid"
x="999" y="999" z="1" lunit="mm"/>

<box
name="horizVirtualPlane_solid"
x="999" y="1" z="8999" lunit="mm"/>


</solids>

<structure>

<volume name="detVirtualPlane_log">
<materialref ref="G4_Galactic"/>
<solidref ref="vertVirtualPlane_solid"/>
<auxiliary auxtype="Visibility" auxvalue="true"/>
<auxiliary auxtype="Color" auxvalue="yellow"/>
<auxiliary auxtype="Alpha" auxvalue="1"/>
<auxiliary auxtype="SensDet" auxvalue="planeDet"/>
<auxiliary auxtype="DetNo" auxvalue="100"/>
</volume>

<volume name="floorVirtualPlane_log">
<materialref ref="G4_Galactic"/>
<solidref ref="horizVirtualPlane_solid"/>
<auxiliary auxtype="Visibility" auxvalue="true"/>
<auxiliary auxtype="Color" auxvalue="yellow"/>
<auxiliary auxtype="Alpha" auxvalue="1"/>
<auxiliary auxtype="SensDet" auxvalue="planeDet"/>
<auxiliary auxtype="DetNo" auxvalue="100"/>
</volume>


<volume name="ParallelWorld_logical">
<materialref ref="G4_Galactic"/>
<solidref ref="parallel_solid"/>

<physvol name="detVirtualPlane_phys">
<volumeref ref="detVirtualPlane_log"/>
<position name="detVirtualPlane_pos" unit="mm" z="4498"/>
</physvol>

<physvol name="floorVirtualPlane_phys">
<volumeref ref="floorVirtualPlane_log"/>
<position name="floorVirtualPlane_pos" unit="mm" x="0" y="-498" z="0"/>
</physvol>

<auxiliary auxtype="Alpha" auxvalue="0.1"/>
</volume>

</structure>

<setup name="Default" version="1.0">
<world ref="ParallelWorld_logical"/>
</setup>

</gdml>
29 changes: 29 additions & 0 deletions scripts/standalone_edep_studies/DSconfig.list
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,29 @@
config_name, DSconfig
single_slab, false
nevents, 2000
particle, e+
energy, 100
n_conductors, 6
single_conductor_height, 9.7
single_conductor_width, 24.0
sideplate_width, 0.2
first_insulation_width, 0.4
second_insulation_width, 0.4
topbot_insulation_width, 1.0
stack_height, 61.75
conductor_hole_radius, 3.05
insulation_mat, mix6337
sideplate_mat, mix6337
conductor_hole_mat, G4_WATER
conductor_mat, G4_Cu
x_rotation, 90
y_rotation, 0
z_rotation, 0
beam_x_pos, 0
beam_y_pos, 0
beam_z_pos, -40
beam_th_dir, 10
beam_ph_dir, 0
max_step_limit_charged_particles, 0.1
run_macro, macros/run.mac
output_geometry, geometry_standalone_edep/mollerMother.gdml
33 changes: 33 additions & 0 deletions scripts/standalone_edep_studies/README.md
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,33 @@
# README

To generate a new config:
```python generateStandalone.py -c <config_file>```

Important variables:

Position and Direction:
For upstream coils, leave x_rotation, y_rotation and z_rotation as 0. To change direction of beam hitting a spot on the coil, you need to modify the beam pos and beam dir variables.

For downstream coils, set x_rotation to 90. Leave the other two as 0. To change direction of beam hitting a spot on the coil, you need to modify the beam pos and beam dir variables.


Charged particle Step Length:
You can check the effect of varying the charged particle step length. This is set using the max_step_limit_charged_particles variable in the config file. Note that this will not set limits on gammas.

Running simulations:
Short jobs (nEvents<2000) can be run from command line:
./build/remoll macros/run.mac

For long jobs (nEvents), adapt and use the given job submission script.

Detector IDs:
1, 6 = Outer most insulation layer
2, 5 = Side Plates
3, 4 = Inner most insulation layer
7, 8 = Top and bottom insulation layer

50-59: Conductor hole
60-69: Conductor
70-79: Gaps between conductor layers filled with insulation


29 changes: 29 additions & 0 deletions scripts/standalone_edep_studies/USconfig.list
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,29 @@
config_name, USconfig1
single_slab, false
nevents, 100
particle, e+
energy, 100
n_conductors, 5
single_conductor_height, 9.0
single_conductor_width, 9.0
sideplate_width, 0.5
first_insulation_width, 0.5
second_insulation_width, 0.5
topbot_insulation_width, 1.0
stack_height, 47.845
conductor_hole_radius, 2.85
insulation_mat, mix6337
sideplate_mat, G4_W
conductor_hole_mat, G4_WATER
conductor_mat, G4_Cu
x_rotation, 0
y_rotation, 0
z_rotation, 0
beam_x_pos, -220.81
beam_y_pos, 59.17
beam_z_pos, -20
beam_th_dir, 85
beam_ph_dir, -15
run_macro, macros/run.mac
output_geometry, geometry_standalone_edep/mollerMother.gdml
max_step_limit_charged_particles, 0.1
Loading