diff --git a/compact/LGC.xml b/compact/LGC.xml
index 15bb04b90d7995c0cc7e21f28eea67fb7a98bd5e..988048f879f63e63f23be4da22f248bc89c33206 100644
--- a/compact/LGC.xml
+++ b/compact/LGC.xml
@@ -190,7 +190,7 @@
     <rotation x="65*degree" y="0" z="0" />
     <comment> TODO: implement PMT material layers </comment>
     <pmt_array dx="16*cm" dy="16*cm" sensitive="yes" vis="LGCPMTVis" />
-    <cone shape="cone" thickness="4*mm" rmin="7.0*cm" rmax="21.5*cm" length="30.*cm" vis="LGCConeVis" />
+    <cone shape="paraboloid" thickness="4*mm" rmin="7.0*cm" rmax="21.5*cm" length="30.*cm" vis="LGCConeVis" />
     <shield material="Lead" thickness="1.*cm" radius="11.53*cm" length="8.*cm"
       shift_z="1.5*cm" vis="LGCShieldVis" />
   </winston_cone>
diff --git a/scripts/subsystems/lgc_sim.py b/scripts/LGC/sim.py
similarity index 90%
rename from scripts/subsystems/lgc_sim.py
rename to scripts/LGC/sim.py
index 39f4e1c423207cae1231f2f02b86b5e04e15990c..7cc1cef09df9dab86395f5e8716e0ca6bc14d7ed 100644
--- a/scripts/subsystems/lgc_sim.py
+++ b/scripts/LGC/sim.py
@@ -5,9 +5,9 @@ Based on M. Frank and F. Gaede runSim.py
 Modified with settings for SoLID simulation
 
 [simulation]:
-    python scripts/subsystems/lgc_sim.py --compactFile solid.xml --runType batch
+    python scripts/LGC/sim.py --compactFile solid.xml --runType batch
 [visualization]:
-    python scripts/subsystems/lgc_sim.py --compactFile solid.xml --runType qt --macro macro/vis.mac
+    python scripts/LGC/sim.py --compactFile solid.xml --runType qt --macro macro/vis.mac
 """
 from __future__ import absolute_import, unicode_literals
 import logging
@@ -65,13 +65,13 @@ if __name__ == "__main__":
     SIM.part.userParticleHandler = ""
 
     # Particle gun settings: electrons with fixed energy and theta, varying phi
-    SIM.numberOfEvents = 100
+    SIM.numberOfEvents = 100000
     SIM.enableGun = True
     SIM.gun.position = (0., 0., "-300*cm")
     SIM.gun.energy = "5*GeV"
     SIM.gun.particle = "e-"
-    SIM.gun.thetaMin = "12.0*deg"
-    SIM.gun.thetaMax = "12.1*deg"
+    SIM.gun.thetaMin = "8.0*deg"
+    SIM.gun.thetaMax = "16.0*deg"
     SIM.gun.distribution = "cos(theta)"
 
     # Output file (assuming CWD)