|
Updates and Results
Talks and Posters
Advice
Ideas
Important Figures
Write-Ups
Outreach
How-To
Funding Opportunities
GENETIS
|
| Important Plots, Tables, and Measurements |
 |
|
|
 Sun Apr 23 14:54:50 2017, Hannah Hasan, Other, Simulation, Plotting ShelfMC Parameter Space, Other 
|
 Fri Oct 6 15:15:53 2017, Hannah Hasan, Other, Simulation, Plotting ShelfMC Parameter Space, Other 6x
|
|
|
Message ID: 27
Entry time: Fri Oct 6 15:15:53 2017
In reply to: 16
|
| Author: |
Hannah Hasan |
| Type: |
Other |
| Category: |
Simulation |
| Subject: |
Plotting ShelfMC Parameter Space |
| Project: |
Other |
|
|
Attached are instructions and scripts for carrying out a parameter space scan with ShelfMC, the simulation package for the ARIANNA detector.
Because some of the plotted outputs looked like colored stripes and did not offer any insight into how effective volume changed with some variables, I made some changes to the simulation and plotting scripts so that different maximum, minimum, and increment values can be chosen for each variable. Now rather than having fixed, hard-coded values for all variables, the parameter space scan and plotting is more flexible for use with variable inputs.
| Quote: |
|
I am trying to write a script that will plot a 2d histogram of effective volume versus two of ShelfMC's parameters.
The script prompts the user for which two parameters (out of five that we vary in our parameter space scan) to plot along the x- and y-axes, as well as what values to hold the other 3 parameters constant at. It then collects the necessary root files from simulation results, generates a plotting script, and runs the plotting script to produce a plot in pdf form.
After many struggles I have the script written to the point where it functions, but the plots don't look right. Some plots look like they could be actual data (like Veff_A_I), and others just look flat-out wrong (like Veff_R_S).
I have yet to pin down the cause of this, but hopefully will be able to sometime in the near future.
|
|
|
|
|
#!/bin/bash
#Jude Rajasekera 3/20/17
tmpShelfmc=$HOME/ShelfMC/git_shelfmc #location of Shelfmc
runName=e18_ParamSpaceScan #name of run
cd $tmpShelfmc #move to home directory
AttenMin=500 # MINIMUM attenuation length for the simulation
AttenMax=1000 # MAXIMUM attenuation length
AttenInc=100 # INCREMENT for attenuation length
RadiusMin=3 # MINIMUM radius
RadiusMax=31 # MAXIMUM radius
RadiusInc=7 # INCREMENT
IceMin=500 # etc...
IceMax=2900
IceInc=400
FirnMin=60
FirnMax=140
FirnInc=20
StDepthMin=0
StDepthMax=200
StDepthInc=50
if [ ! -f ./jobList.txt ]; then #see if there is an existing job file
echo "Creating new job List"
for ((L=$AttenMin;L<=$AttenMax;L+=$AttenInc)) #Attenuation length
do
for ((AR=$RadiusMin;AR<=$RadiusMax;AR+=$RadiusInc)) #Station radius (measured from center of radius to antenna
do
for ((T=$IceMin;T<=$IceMax;T+=$IceInc)) #Ice thickness
do
for ((FT=$FirnMin;FT<=$FirnMax;FT+=$FirnInc)) #Firn thickness
do
for ((SD=$StDepthMin;SD<=$StDepthMax;SD+=$StDepthInc)) #Station depth
do
echo "cd $runName/Atten_Up$L/AntennaRadius$AR/IceThick$T/FirnThick$FT/StationDepth$SD" >> jobList.txt
done
done
done
done
done
else
echo "Picking up from last job"
fi
numbLeft=$(wc -l < ./jobList.txt)
while [ $numbLeft -gt 0 ];
do
jobs=$(showq | grep "cond0092") #change username here
echo '__________Current Running Jobs__________'
echo "$jobs"
echo ''
runningJobs=$(showq | grep "cond0092" | wc -l) #change username here
echo Number of Running Jobs = $runningJobs
echo Number of jobs left = $numbLeft
if [ $runningJobs -le 20 ];then
line=$(head -n 1 jobList.txt)
$line
echo Submit Job && pwd
qsub run_shelfmc_multithread.sh
cd $tmpShelfmc
sed -i 1d jobList.txt
else
echo "Full Capacity"
fi
sleep 1
numbLeft=$(wc -l < ./jobList.txt)
done
|
|
|
|