cone beam

MonteCarlo/config/config.cfg

@@ -20,23 +20,28 @@ ForceDetectorPosition:           YES
 
 EMCALx:							0
 EMCALy:							0
-EMCALz:							0.0
+EMCALz:							1825.0
+#0.0
 
 RPDx:							  0
 RPDy:							  0
-RPDz:							  190.0
+RPDz:							  400.0
+#190 
 
 ZDC1x:							0
 ZDC1y:							0
-ZDC1z:							1525.0
+ZDC1z:							0.0
+#1525
 
 ZDC2x:							0
 ZDC2y:							0
-ZDC2z:							1675.0
+ZDC2z:							150.0
+#1675
 
 ZDC3x:							0
 ZDC3y:							0
-ZDC3z:							1825.0
+ZDC3z:							300.0
+%1825
 
 AbsorberZ:				  0
 

MonteCarlo/run2.mac

@@ -27,7 +27,7 @@
 
 # the beam energy is in gaussian profile
 /gps/ene/type Gauss
-/gps/ene/mono 2.5 TeV
+/gps/ene/mono 500 GeV
 /gps/ene/sigma 1 MeV
 
 
@@ -38,11 +38,17 @@
 #/gps/pos/halfx 20 mm
 #/gps/pos/halfy 20 mm
 
-/gps/pos/type Beam
-/gps/pos/shape Circle
-/gps/pos/centre 0. 0. -1500. mm
-/gps/pos/radius 1. mm
-/gps/pos/sigma_r 2. mm
+#/gps/pos/type Beam
+#/gps/pos/shape Circle
+#/gps/pos/centre 0. 0. -1500. 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
 
 #
 #

MonteCarlo/src/DetectorConstruction.cc

@@ -171,7 +171,7 @@ G4VPhysicalVolume* DetectorConstruction::ConstructDetector()
   ZDC2  = true;
   ZDC3  = true;
   EMCAL = true;
-  RPD   = false;
+  RPD   = true;
 
 	leadblockZ = (AbsorberZ)*mm;
 	lead_in = 	align_run->lead_In;