# A. Tramontana
# tramontana@lns.infn.it
#
# Macro file recommended for the use with laser-driven beam line and built-in physic list.
# This macro uses the physics list already contained in the Geant4 distribution

#######################
# Set of the verboses # 
#######################
#
/control/verbose 1
/tracking/verbose 0
/run/verbose 1
/event/verbose 0

###########################################################################################################
#					ACTIVATED PHYSICS MODELS


/Physics/addPhysics standard_opt4

/geometrySetup/selectGeometry LaserDriven

# Initialize geometry and physic 
/run/initialize

# Set a very high time threshold to allow all decays to happen
/process/had/rdm/thresholdForVeryLongDecayTime 1.0e+60 year

########################################################################################################## 

/gps/particle proton


/gps/pos/type Beam
/gps/pos/shape Circle

/gps/pos/centre -171.3 0.0 0.0 cm   
/gps/pos/radius 0.0 mm   
/gps/pos/sigma_r 0.1 mm  

# the incident surface is in the y-z plane
#
/gps/pos/rot1 0 1 0
/gps/pos/rot2 0 0 1

# the beam is travelling along the x-axis without any angular dispersion (angular despersion set to 0.0)
#
/gps/ang/rot1 0 0 1
/gps/ang/rot2 0 1 0
/gps/ang/type beam1d
/gps/ang/sigma_r 0.028 deg

# ###########the beam energy is in gaussian profile
#/gps/ene/type Gauss
#/gps/ene/mono 7 MeV  
#/gps/ene/sigma 0.0 MeV
#/gps/ang/type iso
#/gps/ang/mintheta 0 deg 
#/gps/ang/maxtheta .1 deg 

################### the beam energy is uniform profile from 4.0 MeV to 12 MeV
#/gps/ene/type User
#/gps/ene/type User
#/gps/hist/type energy
#/gps/hist/point 4.000001 0
#/gps/hist/point 12. 1
#/gps/direction 1 0 0
 
#####################################################################
# TARANIS energy spectrum without angular dipendence: exp(-.44E) for energy between 1. to 12 MeV
/gps/ene/type Exp
/gps/ene/min 1. MeV
/gps/ene/max 12. MeV
/gps/ene/ezero 2.272727273

######################################################################
/vis/open OGLI 600x600-0+0

# Disable auto refresh and quieten vis messages whilst scene and
# trajectories are established:
/vis/viewer/set/autoRefresh false
/vis/verbose errors

/vis/drawVolume
/vis/viewer/set/background white
/vis/viewer/set/viewpointThetaPhi 60 130
/vis/viewer/zoom 1.0
/vis/viewer/panTo -.27 -.0
/vis/scene/add/axes 0 0 0 1 m
/vis/scene/add/trajectories
/vis/scene/endOfRunAction accumulate 
/vis/modeling/trajectories/create/drawByCharge
/vis/scene/endOfEventAction accumulate -1
/vis/viewer/update



##########################################################################################################
############### Define the entrance beam: TARANIS angular distribution
##  beam->1-2
#/gps/verbose 2
#/gps/source/intensity 3.6
#/gps/particle proton
#/gps/pos/type Beam
#/gps/pos/shape Circle
#/gps/pos/centre -171.3 0.0 0.0 cm
#/gps/pos/radius 0.007 mm
#/gps/pos/sigma_r 0. mm
#/gps/ang/rot1 0 0 1
#/gps/ang/rot2 0 1 0
#/gps/pos/rot1 0 0 1
#/gps/pos/rot2 0 1 0
#/gps/ang/type iso
#/gps/ang/mintheta 0 deg
#/gps/ang/maxtheta 26.69 deg
#/gps/ene/type Arb
#/gps/hist/type arb
#/gps/hist/point 1.1 3.5
#/gps/hist/point 3.1399999999999999 1.2
#/gps/hist/inter Exp


################# change first collimator #################
#/LaserDriven/EnergySelector/FirstCollimator/Radius 10.0 mm
##/LaserDriven/EnergySelector/FirstCollimator/Thickness 20.0 mm
##/LaserDriven/EnergySelector/FirstCollimator/zPosizion 42. mm
################# change second collimator #################
#/LaserDriven/EnergySelector/SecondCollimator/Radius 20.0 mm
#/LaserDriven/EnergySelector/SecondCollimator/Thickness 20.0 mm
#/LaserDriven/EnergySelector/SecondCollimator/zPosizion 42. mm
######################## change slit ########################
##/LaserDriven/EnergySelector/Slit/thickness 0.1 mm
/LaserDriven/EnergySelector/Slit/HoleDimensionY 8. mm
/LaserDriven/EnergySelector/Slit/HoleDimensionZ 1. mm
/LaserDriven/EnergySelector/Slit/HolePositionZ 35. mm
#
#/LaserDriven/EnergySelector/Disable
#/LaserDriven/Quadrupoles/DisableQuads
#############################################################################################
# CUT AND STEP MAX 
# Suggested values of cut and step:
#
# Set cuts OUTSIDE the detector
#/run/setCut 1 mm
# Set cuts ONLY inside the detector 
#
#/run/setCutForRegion DetectorLog 0.5 mm

/Step/waterPhantomStepMax 1 mm 


##################################################				
#	    DETECTOR DESCRIPTION
###################################################
# USE THE FOLLOWING FOR A VOXELLIZED DETECTOR
###################################################
#
/changePhantom/position 100 100 0 cm
/changePhantom/size 40 40 40  cm
/changeDetector/size 4 4 4 cm
/changeDetector/voxelSize 1 40 40 mm

# Put the detector in the lower left corner of the phantom
/changeDetector/displacement 0 18 18 cm
/changePhantom/update
/event/printEventNumber 100

/run/beamOn 500