# run1.mac # Run in batch: ./exampleB1 run1.mac # Or interactively: /control/execute run1.mac ########## Print commands ########## #/control/verbose 2 #/run/verbose 2 #/event/verbose 0 #/tracking/verbose 1 #/process/had/verbose 0 #/process/em/verbose 0 #/tracking/verbose 2 ########## Electron/positron annihiliation to muons bias ########## #/physics_lists/em/PositronToMuonsFactor 1000 #/physics_lists/em/PositronToMuons true ########## Scintillation ########## /process/optical/processActivation Cerenkov false /process/optical/processActivation OpAbsorption true /process/optical/processActivation OpBoundary false /process/optical/processActivation Scintillation true /process/optical/processActivation OpRayleigh false /process/optical/processActivation OpMieHG false /process/optical/processActivation OpWLS false #/process/optical/scintillation/setEnhancedTimeConstants true #/process/optical/scintillation/setYieldFactor 10 #/process/optical/scintillation/setExcitationRatio .5 #/process/optical/scintillation/setByParticleType true #/process/optical/scintillation/setTrackInfo false #/process/optical/scintillation/setFiniteRiseTime false #/process/optical/scintillation/setStackPhotons true #/process/optical/scintillation/setTrackSecondariesFirst true ########## Initilize ########## /run/initialize ########## Scoring mesh Target 1 ########## #/score/create/cylinderMesh boxMesh_1 #/score/mesh/cylinderSize 0.75 0.45 cm #/score/mesh/translate/xyz 0. 0. 0. cm #/score/mesh/nBin 300 30 1 #/score/quantity/energyDeposit eDep #/score/quantity/nOfStep nOfStepGamma #/score/filter/particle gammaFilter gamma #/score/quantity/nOfStep nOfStepEMinus #/score/filter/particle eMinusFilter e- #/score/filter/particle muMinusFilter mu- #/score/quantity/nOfStep nOfStepEPlus #/score/filter/particle ePlusFilter e+ #/score/filter/particle muPlusFilter mu+ #/score/close #/score/list ########## Scoring mesh Target 2 ########## #/score/create/cylinderMesh boxMesh_2 #/score/mesh/cylinderSize 0.75 0.45 cm #/score/mesh/translate/xyz 0. 0. 2.9 cm #/score/mesh/nBin 300 30 1 #/score/quantity/energyDeposit eDep #/score/quantity/nOfStep nOfStepGamma #/score/filter/particle gammaFilter gamma #/score/quantity/nOfStep nOfStepEMinus #/score/filter/particle eMinusFilter e- #/score/filter/particle muMinusFilter mu- #/score/quantity/nOfStep nOfStepEPlus #/score/filter/particle ePlusFilter e+ #/score/filter/particle muPlusFilter mu+ #/score/close #/score/list ########## Scoring mesh Target 3 ########## #/score/create/cylinderMesh boxMesh_3 #/score/mesh/cylinderSize 0.75 0.075 cm #/score/mesh/translate/xyz 0. 0. 2.15 cm #/score/mesh/nBin 300 5 1 #/score/quantity/energyDeposit eDep #/score/quantity/nOfStep nOfStepGamma #/score/filter/particle gammaFilter gamma #/score/quantity/nOfStep nOfStepEMinus #/score/filter/particle eMinusFilter e- #/score/filter/particle muMinusFilter mu- #/score/quantity/nOfStep nOfStepEPlus #/score/filter/particle ePlusFilter e+ #/score/filter/particle muPlusFilter mu+ #/score/close #/score/list ########## Scoring mesh Target 4 ########## #/score/create/cylinderMesh boxMesh_4 #/score/mesh/cylinderSize 0.75 0.075 cm #/score/mesh/translate/xyz 0. 0. 2.15 cm #/score/mesh/nBin 300 5 1 #/score/quantity/energyDeposit eDep #/score/quantity/nOfStep nOfStepGamma #/score/filter/particle gammaFilter gamma #/score/quantity/nOfStep nOfStepEMinus #/score/filter/particle eMinusFilter e- #/score/filter/particle muMinusFilter mu- #/score/quantity/nOfStep nOfStepEPlus #/score/filter/particle ePlusFilter e+ #/score/filter/particle muPlusFilter mu+ #/score/close #/score/list ########## Scoring mesh Target 5 ########## #/score/create/cylinderMesh boxMesh_5 #/score/mesh/cylinderSize 0.75 0.075 cm #/score/mesh/translate/xyz 0. 0. 2.15 cm #/score/mesh/nBin 300 5 1 #/score/quantity/energyDeposit eDep #/score/quantity/nOfStep nOfStepGamma #/score/filter/particle gammaFilter gamma #/score/quantity/nOfStep nOfStepEMinus #/score/filter/particle eMinusFilter e- #/score/filter/particle muMinusFilter mu- #/score/quantity/nOfStep nOfStepEPlus #/score/filter/particle ePlusFilter e+ #/score/filter/particle muPlusFilter mu+ #/score/close #/score/list ########## Scoring mesh Regenerator ########## #/score/create/cylinderMesh boxMesh_reg #/score/mesh/cylinderSize 0.1 1. cm #/score/mesh/translate/xyz 0. 0. 1002.15 cm #/score/mesh/nBin 40 100 1 #/score/quantity/energyDeposit eDep #/score/quantity/nOfStep nOfStepGamma #/score/filter/particle gammaFilter gamma #/score/quantity/nOfStep nOfStepEMinus #/score/filter/particle eMinusFilter e- #/score/filter/particle muMinusFilter mu- #/score/quantity/nOfStep nOfStepEPlus #/score/filter/particle ePlusFilter e+ #/score/filter/particle muPlusFilter mu+ #/score/close #/score/list ########## Scoring mesh Scintillator ########## #/score/create/boxMesh boxMesh_scint #/score/mesh/boxSize 2.5 2.5 0.25 cm #/score/mesh/translate/xyz -11.48 0. -31 cm #/score/mesh/rotate/rotateY -22.5 deg #/score/mesh/nBin 50 50 5 #/score/quantity/energyDeposit eDep #/score/quantity/nOfStep nOfStepGamma #/score/filter/particle gammaFilter gamma #/score/quantity/nOfStep nOfStepEMinus #/score/filter/particle eMinusFilter e- #/score/filter/particle muMinusFilter mu- #/score/quantity/nOfStep nOfStepEPlus #/score/filter/particle ePlusFilter e+ #/score/filter/particle muPlusFilter mu+ #/score/close #/score/list ########## Particle Electrons 1 GeV for Mainz test Flux########## #/gun/particle e- #/gun/energy 1 GeV /gps/particle e- #/gps/energy 1 GeV /gps/energy 800 MeV /gps/ene/type Gauss /gps/ene/sigma 13 keV /gps/pos/type Beam /gps/pos/shape Circle /gps/pos/radius 0. mm #/gps/pos/sigma_r 1 mm /gps/pos/sigma_x 0.01 mm /gps/pos/sigma_y 0.01 mm /gps/ang/sigma_x 0.000001 rad /gps/ang/sigma_y 0.000001 rad /gps/pos/centre 0 0 -100 cm /gps/direction 0 0 1 ########## Particle Positron 45 GeV for LEMMA ########## #/gun/particle e+ #/gun/energy 45 GeV #/gps/particle e+ #/gps/energy 45 GeV #/gps/ene/type Gauss #/gps/ene/sigma 0 GeV #/gps/pos/type Beam #/gps/pos/shape Circle #/gps/pos/radius 0. mm #/gps/pos/sigma_r 1 mm #/gps/pos/sigma_x 0.3 mm #/gps/pos/sigma_y 0.3 mm #/gps/pos/centre 0 0 -100 cm #/gps/direction 0 0 1 ########## Particle Proton 10 GeV for MAP ########## #/gun/particle e+ #/gun/energy 45 GeV #/gps/particle proton #/gps/energy 10 GeV #/gps/ene/type Gauss #/gps/ene/sigma 0 GeV #/gps/pos/type Beam #/gps/pos/shape Circle #/gps/pos/radius 0. mm #/gps/pos/sigma_r 1 mm #/gps/pos/sigma_x 0.3 mm #/gps/pos/sigma_y 0.3 mm #/gps/pos/centre 0 0 -100 cm #/gps/direction 0 0 1 ########## Magnetic Field ########## /B1/field/value 0. tesla ########## Number of particles ########## /run/printProgress 100000 #/run/beamOn 300000000 /run/beamOn 1000000 ########## Scoring Mesh dump ########## #/score/dump #/score/dumpAllQuantitiesToFile boxMesh_1 boxMesh_1.dump #/score/dumpAllQuantitiesToFile boxMesh_2 boxMesh_2.dump #/score/dumpAllQuantitiesToFile boxMesh_3 boxMesh_3.dump #/score/dumpAllQuantitiesToFile boxMesh_4 boxMesh_4.dump #/score/dumpAllQuantitiesToFile boxMesh_5 boxMesh_5.dump #/score/dumpAllQuantitiesToFile boxMesh_reg boxMesh_reg.dump #/score/dumpAllQuantitiesToFile boxMesh_scint boxMesh_scint.dump