cone source fix

MonteCarlo/run2.mac

@@ -23,54 +23,55 @@
 #/gps/energy 600 GeV
 
 #/gps/particle gamma
-/gps/particle neutron
+#/gps/particle neutron
+/gps/particle proton
 
 # the beam energy is in gaussian profile
 /gps/ene/type Gauss
-/gps/ene/mono 500 GeV
+/gps/ene/mono 1 TeV
 /gps/ene/sigma 1 MeV
 
 
 # General particle source
+
+# FUNCTIONING CONE
+/gps/ang/type iso
+/gps/pos/centre 0. 0. -99. mm
+/gps/ang/maxtheta 0.003125 deg
+/gps/ang/rot1 0 1 0
+/gps/ang/rot2 1 0 0
+
+
+
+# FUNCTIONING PLANE
 #/gps/pos/type Plane
 #/gps/pos/shape Rectangle
-#/gps/pos/centre 0. 0. -100. mm
+#/gps/pos/centre 0. 0. -99. mm
 #/gps/pos/halfx 20 mm
 #/gps/pos/halfy 20 mm
+#/gps/ang/rot1 0 1 0
+#/gps/ang/rot2 1 0 0
 
+#/gps/ang/type beam1d
+
+# FUNCTIONING BEAM
 #/gps/pos/type Beam
 #/gps/pos/shape Circle
-#/gps/pos/centre 0. 0. -1500. mm
+#/gps/pos/centre 0. 0. -99. mm
 #/gps/pos/radius 1. mm
 #/gps/pos/sigma_r 2. mm
-
-/gps/ang/type iso
-/gps/ang/maxtheta 10 deg
-/gps/ang/rot1 0 0 1
-/gps/ang/rot2 0 1 0
-/gps/pos/centre 0. 0. -100. mm
-
-#
-#
-# the incident surface is in the x-y plane
-/gps/pos/rot1 1 0 0
-/gps/pos/rot2 0 1 0
-#
-#
-# the beam is travelling along the z-axis
-/gps/ang/rot1 0 1 0
-/gps/ang/rot2 1 0 0
-/gps/ang/type beam1d
-# the beam angular dispersion
-/gps/ang/sigma_r 0.00 deg
+#/gps/ang/rot1 0 1 0
+#/gps/ang/rot2 1 0 0
+#/gps/ang/type beam1d
+#/gps/ang/sigma_r 0.00 deg
 
 
 #if using these commands in interactive mode to visualize, use this to allow more events to be displayed (do not use in batch mode)
-#/vis/ogl/set/displayListLimit 1000000000
+/vis/ogl/set/displayListLimit 1000000000
 
 ##############################################################
 ##############################################################
 # number of events
-/run/beamOn 100
+/run/beamOn 1
 ##############################################################
 ##############################################################

MonteCarlo/vis.mac

@@ -58,16 +58,17 @@
 /vis/scene/add/hits
 #
 
-/vis/filtering/trajectories/create/chargeFilter
-/vis/filtering/trajectories/chargeFilter-0/add 1
+
+#/vis/filtering/trajectories/create/chargeFilter
+#/vis/filtering/trajectories/chargeFilter-0/add 1
 #/vis/filtering/trajectories/chargeFilter-0/invert true
 
 # To draw only gammas:
-#/vis/filtering/trajectories/create/particleFilter
+/vis/filtering/trajectories/create/particleFilter
 #/vis/filtering/trajectories/particleFilter-0/add opticalphoton
 #/vis/filtering/trajectories/particleFilter-0/add e-
 #/vis/filtering/trajectories/particleFilter-0/add e+
-#/vis/filtering/trajectories/particleFilter-0/add proton
+/vis/filtering/trajectories/particleFilter-0/add proton
 #/vis/filtering/trajectories/particleFilter-0/add alpha
 #/vis/filtering/trajectories/particleFilter-0/add GenericIon
 #/vis/filtering/trajectories/particleFilter-0/add gamma
@@ -75,16 +76,12 @@
 #
 # To invert the above, drawing all particles except gammas,
 # keep the above two lines but also add:
-#/vis/filtering/trajectories/particleFilter-0/invert true
+/vis/filtering/trajectories/particleFilter-0/invert true
 
 
 #  To increase display limit
 #/vis/ogl/set/displayListLimit 5000000
-# to filter out uncharge particles
 
-#/vis/filtering/trajectories/create/chargeFilter
-#/vis/filtering/trajectories/chargeFilter-0/add 1
-#/vis/filtering/trajectories/chargeFilter-0/add -1
 
 
 # Many other options are available with /vis/modeling and /vis/filtering.