Added GPD subdirectory.

This commit is contained in:
prokscha 2020-12-03 18:01:12 +01:00
parent 6d16ea7f0a
commit 68d00b1591
20 changed files with 4256 additions and 0 deletions

240
run/GPD/40111.mac Normal file
View File

@ -0,0 +1,240 @@
#-----------------------------------------------------------------------
# Macro file for muSR instrument simulations
# Unless specified otherwises, the default units are mm, ns, MeV, MeV/c.
# Lines starting with star "#" are comments.
###################################################################################
############################# G E O M E T R Y ###################################
###################################################################################
# ROTATION MATRIXES:
/musr/command rotation matrix1 1 0 0 90
/musr/command rotation matrix2 0 1 0 28.81
/musr/command rotation matrix3 0 1 0 151.19
/musr/command rotation matrix4 0 0 1 45
/musr/command rotation matrix5 0 1 0 90
/musr/command rotation matrix6 180 90 90
/musr/command rotation matrix7 0 90 90
# --- WORLD VOLUME ---
/musr/command construct box World 2100 2100 11000 G4_AIR 0 0 0 no_logical_volume norot dead -1
# --- BEAMPIPE ---
/musr/command construct tubs beampipeA 0 2005 4500 0 360 G4_Galactic 0 0 -6100 log_World norot dead 230
/musr/command construct tubs beampipe 0 405 500 0 360 G4_Galactic 0 0 -1100 log_World norot dead 231
/musr/command construct tubs beampipeshell 100 105 500 0 360 G4_Al 0 0 0 log_beampipe norot dead 232
/musr/command construct tubs beampipewindow 0 100 0.05 0 360 G4_MYLAR 0 0 499.9 log_beampipe norot dead 233
# PASSIVE COLLIMATOR (for BEAMPIPE vers.1)
# GPDcollimator - thickness in the thinner location is x2+(x9-x4) = 20+59-59
/musr/command construct GPDcollimator GPDcollimator 40 15 60 59 60.01 0 360 0 59 0 2.5 6 G4_Cu 0 0 -59 log_World matrix1 dead 251
/musr/command construct box GPDcollimator2 70 70 15 G4_Pb 0 0 -100 log_World norot dead 255
/musr/command construct box GPDcollimator2hole 2 5 15 G4_Galactic 0 0 0 log_GPDcollimator2 norot dead 256
/musr/command construct box FirstCollimator 45 75 15 G4_Pb 0 0 -250 log_World norot dead 261
/musr/command construct tubs FirstCollHole 0 8 15 0 360 G4_Galactic 0 0 0 log_FirstCollimator norot dead 262
#/musr/command construct box AlumPlate 30 90 5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -80 log_World norot musr/ScintSD 401
#/musr/command construct box AlumPlateHole 2 5 5 G4_Galactic 0 0 0 log_AlumPlate norot dead 267
/musr/command construct box AlumPlate 60 90 5 G4_Al 0 0 -80 log_World norot dead 265
/musr/command construct tubs AlumPlateHole 0 25 5 0 360 G4_Galactic 0 0 0 log_AlumPlate norot dead 267
# --- SAMPLE HOLDER ---
/musr/command construct GPDsampleHolderA GPDsampleHolderA 0 0 0 0 0 G4_Al 0 -111.25 0 log_World matrix1 dead 300
/musr/command construct tubs GPDsampleHolderB 37.51 58 7 0 360 G4_Al 0 57 0 log_World matrix1 dead 301
/musr/command construct tubs GPDsampleHolderC 37.51 58 7 0 360 G4_Al 0 -57 0 log_World matrix1 dead 302
####/musr/command construct tubs GPDsampleHolderD 37 37.5 50 136 88 G4_Al 0 0 0 log_World matrix1 dead 303
####/musr/command construct tubs GPDsampleHolderE 53.5 54 50 136 88 G4_Al 0 0 0 log_World matrix1 dead 304
/musr/command construct box GPDsampleHolderD 12 0.5 50 G4_Al 38 0 28 log_World matrix1 dead 303
/musr/command construct box GPDsampleHolderE 12 0.5 50 G4_Al 38 0 -28 log_World matrix1 dead 304
/musr/command construct box GPDsampleHolderF 12 0.5 50 G4_Al -38 0 28 log_World matrix1 dead 305
/musr/command construct box GPDsampleHolderG 12 0.5 50 G4_Al -38 0 -28 log_World matrix1 dead 306
/musr/command construct tubs GPDsampleHolderH 0 9 50 180 50 G4_Al 48 0 29 log_World matrix1 dead 307
/musr/command construct tubs GPDsampleHolderI 0 9 50 130 50 G4_Al 48 0 -29 log_World matrix1 dead 308
/musr/command construct tubs GPDsampleHolderJ 0 9 50 310 50 G4_Al -48 0 29 log_World matrix1 dead 309
/musr/command construct tubs GPDsampleHolderK 0 9 50 0 50 G4_Al -48 0 -29 log_World matrix1 dead 310
/musr/command construct tubs GPDsampleHolderW 58.01 58.2 64 25 130 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 313
/musr/command construct tubs GPDsampleHolderX 58.01 58.2 64 205 130 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 314
/musr/command construct tubs GPDsampleHolderY 36.79 36.99 64 45 90 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 315
/musr/command construct tubs GPDsampleHolderZ 36.79 36.99 64 225 90 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 316
# --- SAMPLE CELL ---
/musr/command construct tubs GPDsampleCell 0 12 50 0 360 G4_Cu 0 0 0 log_World matrix1 dead 211
# --- SAMPLE ---
/musr/command construct tubs target 0 3.5 7 0 360 G4_Cu 0 0 0 log_GPDsampleCell norot dead 201
#/musr/command construct tubs targetFieldVol 0 3.5 7 0 360 G4_Cu 0 0 0 log_GPDsampleCell norot dead 201
# --- SCINTILLATOR ---
/musr/command construct box ScintB1 13 45 2.5 G4_PLASTIC_SC_VINYLTOLUENE 22 0 -40 log_World matrix2 musr/ScintSD 1
/musr/command construct box ScintB2 13 45 2.5 G4_PLASTIC_SC_VINYLTOLUENE -22 0 -40 log_World matrix3 musr/ScintSD 2
/musr/command construct box ScintF1 13 45 2.5 G4_PLASTIC_SC_VINYLTOLUENE -22 0 40 log_World matrix2 musr/ScintSD 11
/musr/command construct box ScintF2 13 45 2.5 G4_PLASTIC_SC_VINYLTOLUENE 22 0 40 log_World matrix3 musr/ScintSD 12
/musr/command construct box ScintF3 7 45 2.5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 46 log_World norot musr/ScintSD 13
#/musr/command construct box KoincB1 20 90 2.5 G4_PLASTIC_SC_VINYLTOLUENE 33 0 -60 log_World matrix2 musr/ScintSD 21
#/musr/command construct box KoincB2 20 90 2.5 G4_PLASTIC_SC_VINYLTOLUENE -33 0 -60 log_World matrix3 musr/ScintSD 22
# --- M0 ---
/musr/command construct box M0 5 5 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -46 log_World matrix4 musr/ScintSD 102
/musr/command construct GPDmHolder GPDmHolder 5 25 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 -25 -46 log_World norot dead 103
# --- Magnet ---
/musr/command construct tubs Magnet1 180 375 78 0 360 G4_Cu 133 0 0 log_World matrix5 dead 291
/musr/command construct tubs Magnet2 180 375 78 0 360 G4_Cu -133 0 0 log_World matrix5 dead 292
/musr/command construct uprofile Uprof1 25 10 400 5 G4_Al -45 0 200 log_World matrix6 dead 293
/musr/command construct uprofile Uprof2 25 10 400 5 G4_Al -45 0 -200 log_World matrix6 dead 293
/musr/command construct uprofile Uprof3 25 10 400 5 G4_Al 45 0 200 log_World matrix7 dead 293
/musr/command construct uprofile Uprof4 25 10 400 5 G4_Al 45 0 -200 log_World matrix7 dead 293
#
##---------------------------------------------------------
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_1 6
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_2 5
#/musr/command logicalVolumeToBeReweighted mu log_activeCollim 4
##---------------------------------------------------------
/musr/command visattributes log_World invisible
/musr/command visattributes log_beampipewindow red
/musr/command visattributes log_GPDsampleHolderD red
/musr/command visattributes log_GPDsampleHolderE red
/musr/command visattributes log_GPDsampleHolderF red
/musr/command visattributes log_GPDsampleHolderG red
/musr/command visattributes log_GPDcollimator green
/musr/command visattributes log_GPDcollimator2 blue
/musr/command visattributes log_GPDmHolder magenta
/musr/command visattributes log_target blue
/musr/command visattributes log_M0 yellow
#/musr/command visattributes G4_PLASTIC_SC_VINYLTOLUENE red
#/musr/command visattributes G4_Pb blue
####################################################################################
######################### M A G N E T I C F I E L D #########################
###################################################################################
# Set magnetic field (set field intensity in T and sigma in mm)
# syntax for magneticfield: fromfile filename fieldValue
# uniform fieldValue
# gaussian fieldValue sigma
#
/musr/command globalfield centralSolenoidField 0. 0. 0. fromfile 2D GPD_coil_primitive.table log_target 0.03
/musr/command globalfield setparameter SetLargestAcceptableStep 200
/musr/command globalfield setparameter SetMinimumEpsilonStep 5e-5
/musr/command globalfield setparameter SetMaximumEpsilonStep 0.001
/musr/command globalfield setparameter SetDeltaOneStep 0.1
/musr/command globalfield setparameter SetDeltaIntersection 0.01
/musr/command globalfield printparameters
/musr/command globalfield printFieldValueAtPoint 0 0 0
/musr/command globalfield printFieldValueAtPoint 0 0 10
/musr/command globalfield printFieldValueAtPoint 0 0 20
/musr/command globalfield printFieldValueAtPoint 0 0 30
/musr/command globalfield printFieldValueAtPoint 0 0 50
/musr/command globalfield printFieldValueAtPoint 0 0 100
/musr/command globalfield printFieldValueAtPoint 0 0 200
/musr/command globalfield printFieldValueAtPoint 0 0 300
/musr/command globalfield printFieldValueAtPoint 0 0 400
/musr/command globalfield printFieldValueAtPoint 0 0 500
/musr/command globalfield printFieldValueAtPoint 10 0 0
/musr/command globalfield printFieldValueAtPoint 100 0 0
/musr/command globalfield printFieldValueAtPoint 300 0 0
/musr/command globalfield printFieldValueAtPoint 0 10 0
/musr/command globalfield printFieldValueAtPoint 0 100 0
/musr/command globalfield printFieldValueAtPoint 0 300 0
###################################################################################
######################### P H Y S I C S P R O C E S S E S ##################
###################################################################################
# Geant 4.9.3
/musr/command process addDiscreteProcess gamma G4PhotoElectricEffect
/musr/command process addDiscreteProcess gamma G4ComptonScattering
/musr/command process addDiscreteProcess gamma G4GammaConversion
/musr/command process addDiscreteProcess gamma G4RayleighScattering
/musr/command process addProcess e- G4eMultipleScattering -1 1 1
/musr/command process addProcess e- G4eIonisation -1 2 2
/musr/command process addProcess e- G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eMultipleScattering -1 1 1
/musr/command process addProcess e+ G4eIonisation -1 2 2
/musr/command process addProcess e+ G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eplusAnnihilation 0 -1 4
/musr/command process addProcess mu- G4MuMultipleScattering -1 1 1
/musr/command process addProcess mu- G4MuIonisation -1 2 2
/musr/command process addProcess mu- G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu- G4MuPairProduction -1 4 4
/musr/command process addProcess mu+ G4MuMultipleScattering -1 1 1
/musr/command process addProcess mu+ G4MuIonisation -1 2 2
/musr/command process addProcess mu+ G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu+ G4MuPairProduction -1 4 4
#
#/home/install/geant4.9.3/source/processes/electromagnetic/utils/src/G4EnergyLossMessenger.cc
######## /process/msc/StepLimit Minimal | UseDistanceToBoundary | UseSafety
/process/msc/StepLimit UseSafety
#/process/msc/LateralDisplacement
/process/msc/RangeFactor 0.04
/process/msc/GeomFactor 2.5
#/process/msc/FactorForAngleLimit
/process/msc/Skin 3.0
#/process/msc/ThetaLimit 0.2 rad
###################################################################################
################## S O M E O T H E R P A R A M E T E R S ##################
###################################################################################
# Set the overall range cut
#/run/setCut 1 mm
# Set the range cut on a particular volumes (in mm)
#/musr/command SetUserLimits log_target 0.01
#/musr/command SetUserLimits log_targetscint 0.01
#/musr/command SetUserLimits log_M0 0.01
# Store all events into the ROOT tree or just the interesting ones ? (true is default)
/musr/command storeOnlyEventsWithHits false
# Set the minimum time separation between two subsequent signals in the same detector (in ns)
/musr/command signalSeparationTime 10
#
/musr/run/howOftenToPrintEvent 10000
/musr/command maximumRunTimeAllowed 220000
/musr/run/randomOption 2
###################################################################################
######################### R O O T O U T P U T ##############################
###################################################################################
/musr/command rootOutput det_VvvKine off
/musr/command rootOutput det_VvvX off
/musr/command rootOutput det_VvvY off
/musr/command rootOutput det_VvvZ off
/musr/command rootOutput det_VvvVolID off
/musr/command rootOutput det_VvvProcID off
/musr/command rootOutput det_VvvTrackID off
/musr/command rootOutput det_VvvParticleID off
/musr/command rootOutput nOptPhot off
/musr/command rootOutput odet_ID off
/musr/command rootOutput odet_nPhot off
/musr/command rootOutput odet_timeFirst off
/musr/command rootOutput odet_timeA off
/musr/command rootOutput odet_timeB off
/musr/command rootOutput odet_timeC off
/musr/command rootOutput odet_timeD off
/musr/command rootOutput odet_timeE off
/musr/command rootOutput odet_timeLast off
###################################################################################
######################### V I S U A L I S A T I O N ##############################
###################################################################################
#/vis/disable
#/control/execute visFromToni.mac
#/control/execute visDawn1.mac
/control/execute visVRML.mac
###################################################################################
######################### P A R T I C L E G U N #################################
###################################################################################
/gun/vertex 0 0 -1000 mm
# FWHM 10mm ==> sigma = 10/2.354 = 4.2481mm
/gun/vertexsigma 25 25 0 mm
#---/gun/vertexboundary: rMaxAllowed, zMinAllowed, zMaxAllowed
/gun/vertexboundary 100 -999999 999999 mm
/gun/momentum 100.0 MeV
# sigma = 3% ==> sigma 27*0.03 = 0.81
/gun/momentumsmearing 3.00 MeV
#---/gun/momentumboundary: pMinAllowed, pMaxAllowed, dummy
#/gun/momentumboundary 20 40 0 MeV
# TURTLE
#/gun/turtlefilename FOR061_2008_04_22.DAT
#/gun/turtlefilename FOR070_2008_10_17_XXII.DAT
#/gun/turtlefilename FOR077_pie3_HiField_d05_x30.dat
#/gun/turtlefilename FOR077_reggiani_Jan2010_NEW.dat
#/gun/turtleZ0position -900 mm
#/gun/turtlefilename FOR070_GPD_beam_profile_45mm.dat
#/gun/turtleZ0position -500 mm
#---/gun/tilt: xangle, yangle, dummy
#/gun/tilt 0 0.5 0 deg
#---/gun/tiltsigma: xangleSigma, yangleSigma, dummy (1 degree on 1 meter ~ 17mm)
/gun/tiltsigma 0.15 0.15 0 deg
/gun/pitch 0.08 deg
/gun/muonPolarizVector 0 0 -1
#---/gun/decaytimelimits: decayMin, decayMax, decayTime
/gun/decaytimelimits -1 -1 2197.03 ns
###################################################################################
######################## B E A M O N #########################################
###################################################################################
#/run/beamOn 16000000
/run/beamOn 5

260
run/GPD/40505.mac Normal file
View File

@ -0,0 +1,260 @@
#-----------------------------------------------------------------------
# Macro file for muSR instrument simulations
# Unless specified otherwises, the default units are mm, ns, MeV, MeV/c.
# Lines starting with star "#" are comments.
###################################################################################
############################# G E O M E T R Y ###################################
###################################################################################
# ROTATION MATRIXES:
/musr/command rotation matrix1 1 0 0 90
/musr/command rotation matrix2 0 1 0 28.81
/musr/command rotation matrix3 0 1 0 151.19
/musr/command rotation matrix4 0 0 1 45
/musr/command rotation matrix5 0 1 0 90
/musr/command rotation matrix6 180 90 90
/musr/command rotation matrix7 0 90 90
# --- WORLD VOLUME ---
/musr/command construct box World 2100 2100 11000 G4_AIR 0 0 0 no_logical_volume norot dead -1
# --- BEAMPIPE ---
/musr/command construct tubs beampipeA 0 2005 4500 0 360 G4_Galactic 0 0 -6100 log_World norot dead 230
/musr/command construct tubs beampipe 0 405 500 0 360 G4_Galactic 0 0 -1100 log_World norot dead 231
/musr/command construct tubs beampipeshell 100 105 500 0 360 G4_Al 0 0 0 log_beampipe norot dead 232
/musr/command construct tubs beampipewindow 0 100 0.05 0 360 G4_MYLAR 0 0 499.9 log_beampipe norot dead 233
# PASSIVE COLLIMATOR (for BEAMPIPE vers.1)
# GPDcollimator - thickness in the thinner location is x2+(x9-x4) = 20+59-59
##/musr/command construct GPDcollimator GPDcollimator 40 15 60 59 60.01 0 360 0 59 0 2.5 6 G4_Cu 0 0 -59 log_World matrix1 dead 251
/musr/command construct box GPDcollimator2 48 55 15 G4_Pb 0 0 -66 log_World norot dead 255
/musr/command construct box GPDcollimator2hole 2.5 4 15 G4_Galactic 0 0 0 log_GPDcollimator2 norot dead 256
/musr/command construct box FirstCollimator 45 75 15 G4_Pb 0 0 -250 log_World norot dead 261
/musr/command construct tubs FirstCollHole 0 8 15 0 360 G4_Galactic 0 0 0 log_FirstCollimator norot dead 262
###/musr/command construct box AlumPlate 30 90 5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -80 log_World norot musr/ScintSD 401
###/musr/command construct box AlumPlateHole 2 5 5 G4_Galactic 0 0 0 log_AlumPlate norot dead 267
/musr/command construct box Insert 20 20 10 G4_Pb 0 0 -91 log_World norot dead 265
/musr/command construct box InsertHole 2.5 4 10 G4_Galactic 0 0 0 log_Insert norot dead 267
# --- SAMPLE HOLDER ---
##/musr/command construct GPDsampleHolderA GPDsampleHolderA 0 0 0 0 0 G4_Al 0 -111.25 0 log_World matrix1 dead 300
##/musr/command construct tubs GPDsampleHolderB 37.51 58 7 0 360 G4_Al 0 57 0 log_World matrix1 dead 301
##/musr/command construct tubs GPDsampleHolderC 37.51 58 7 0 360 G4_Al 0 -57 0 log_World matrix1 dead 302
######/musr/command construct tubs GPDsampleHolderD 37 37.5 50 136 88 G4_Al 0 0 0 log_World matrix1 dead 303
######/musr/command construct tubs GPDsampleHolderE 53.5 54 50 136 88 G4_Al 0 0 0 log_World matrix1 dead 304
##/musr/command construct box GPDsampleHolderD 12 0.5 50 G4_Al 38 0 28 log_World matrix1 dead 303
##/musr/command construct box GPDsampleHolderE 12 0.5 50 G4_Al 38 0 -28 log_World matrix1 dead 304
##/musr/command construct box GPDsampleHolderF 12 0.5 50 G4_Al -38 0 28 log_World matrix1 dead 305
##/musr/command construct box GPDsampleHolderG 12 0.5 50 G4_Al -38 0 -28 log_World matrix1 dead 306
##/musr/command construct tubs GPDsampleHolderH 0 9 50 180 50 G4_Al 48 0 29 log_World matrix1 dead 307
##/musr/command construct tubs GPDsampleHolderI 0 9 50 130 50 G4_Al 48 0 -29 log_World matrix1 dead 308
##/musr/command construct tubs GPDsampleHolderJ 0 9 50 310 50 G4_Al -48 0 29 log_World matrix1 dead 309
##/musr/command construct tubs GPDsampleHolderK 0 9 50 0 50 G4_Al -48 0 -29 log_World matrix1 dead 310
##/musr/command construct tubs GPDsampleHolderW 58.01 58.2 64 25 130 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 313
##/musr/command construct tubs GPDsampleHolderX 58.01 58.2 64 205 130 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 314
##/musr/command construct tubs GPDsampleHolderY 36.79 36.99 64 45 90 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 315
##/musr/command construct tubs GPDsampleHolderZ 36.79 36.99 64 225 90 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 316
# --- SAMPLE CELL ---
/musr/command construct tubs GPDsampleCell 0 12 50 0 360 G4_Cu 0 0 0 log_World matrix1 dead 211
# --- SAMPLE ---
/musr/command construct tubs target 0 3 5 0 360 G4_Cu 0 0 0 log_GPDsampleCell norot dead 201
#/musr/command construct tubs targetFieldVol 0 3 5 0 360 G4_Cu 0 0 0 log_GPDsampleCell norot dead 201
# --- SCINTILLATOR ---
/musr/command construct box ScintB1 21.5 50 2 G4_PLASTIC_SC_VINYLTOLUENE 26.5 0 -44 log_World norot musr/ScintSD 1
/musr/command construct box ScintB2 21.5 50 2 G4_PLASTIC_SC_VINYLTOLUENE -26.5 0 -44 log_World norot musr/ScintSD 2
/musr/command construct box ScintB4 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE 44 0 -25 log_World norot musr/ScintSD 4
/musr/command construct box ScintB5 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE -44 0 -25 log_World norot musr/ScintSD 5
/musr/command construct box ScintB6 2 50 13 G4_PLASTIC_SC_VINYLTOLUENE 44 0 0 log_World norot musr/ScintSD 6
/musr/command construct box ScintB7 2 50 13 G4_PLASTIC_SC_VINYLTOLUENE -44 0 0 log_World norot musr/ScintSD 7
/musr/command construct box ScintF1 20 50 2 G4_PLASTIC_SC_VINYLTOLUENE 28 0 44 log_World norot musr/ScintSD 11
/musr/command construct box ScintF2 20 50 2 G4_PLASTIC_SC_VINYLTOLUENE -28 0 44 log_World norot musr/ScintSD 12
/musr/command construct box ScintF3 7 50 2 G4_PLASTIC_SC_VINYLTOLUENE 0 0 44 log_World norot musr/ScintSD 13
/musr/command construct box ScintF4 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE 44 0 25 log_World norot musr/ScintSD 14
/musr/command construct box ScintF5 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE -44 0 25 log_World norot musr/ScintSD 15
#/musr/command construct box KoincB1 20 90 2 G4_PLASTIC_SC_VINYLTOLUENE 33 0 -60 log_World matrix2 musr/ScintSD 21
#/musr/command construct box KoincB2 20 90 2 G4_PLASTIC_SC_VINYLTOLUENE -33 0 -60 log_World matrix3 musr/ScintSD 22
# --- M0 ---
/musr/command construct box M0holder 3 4 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -20 log_World norot musr/ScintSD 103
/musr/command construct box M0 4 30 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -44 log_World norot musr/ScintSD 102
#/musr/command construct GPDmHolder GPDmHolder 5 25 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 -25 -43 log_World norot dead 103
# --- Magnet ---
/musr/command construct tubs Magnet1 180 375 78 0 360 G4_Cu 133 0 0 log_World matrix5 dead 291
/musr/command construct tubs Magnet2 180 375 78 0 360 G4_Cu -133 0 0 log_World matrix5 dead 292
/musr/command construct uprofile Uprof1 25 10 400 5 G4_Al -45 0 200 log_World matrix6 dead 293
/musr/command construct uprofile Uprof2 25 10 400 5 G4_Al -45 0 -200 log_World matrix6 dead 293
/musr/command construct uprofile Uprof3 25 10 400 5 G4_Al 45 0 200 log_World matrix7 dead 293
/musr/command construct uprofile Uprof4 25 10 400 5 G4_Al 45 0 -200 log_World matrix7 dead 293
#
##---------------------------------------------------------
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_1 6
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_2 5
#/musr/command logicalVolumeToBeReweighted mu log_activeCollim 4
##---------------------------------------------------------
/musr/command visattributes log_World invisible
/musr/command visattributes log_beampipeA invisible
/musr/command visattributes log_beampipe invisible
/musr/command visattributes log_beampipeshell invisible
/musr/command visattributes log_beampipewindow invisible
/musr/command visattributes log_Magnet1 invisible
/musr/command visattributes log_Magnet2 invisible
/musr/command visattributes log_Uprof1 invisible
/musr/command visattributes log_Uprof2 invisible
/musr/command visattributes log_Uprof3 invisible
/musr/command visattributes log_Uprof4 invisible
/musr/command visattributes log_GPDsampleHolderD red
/musr/command visattributes log_GPDsampleHolderE red
/musr/command visattributes log_GPDsampleHolderF red
/musr/command visattributes log_GPDsampleHolderG red
/musr/command visattributes log_GPDcollimator green
/musr/command visattributes log_GPDcollimator2 blue
/musr/command visattributes log_GPDmHolder magenta
/musr/command visattributes log_target blue
/musr/command visattributes log_M0 yellow
#/musr/command visattributes G4_PLASTIC_SC_VINYLTOLUENE red
#/musr/command visattributes G4_Pb blue
####################################################################################
######################### M A G N E T I C F I E L D #########################
###################################################################################
# Set magnetic field (set field intensity in T and sigma in mm)
# syntax for magneticfield: fromfile filename fieldValue
# uniform fieldValue
# gaussian fieldValue sigma
#
/musr/command globalfield centralSolenoidField 0. 0. 0. fromfile 2D GPD_coil_primitive.table log_target 0.03
/musr/command globalfield setparameter SetLargestAcceptableStep 200
/musr/command globalfield setparameter SetMinimumEpsilonStep 5e-5
/musr/command globalfield setparameter SetMaximumEpsilonStep 0.001
/musr/command globalfield setparameter SetDeltaOneStep 0.1
/musr/command globalfield setparameter SetDeltaIntersection 0.01
/musr/command globalfield printparameters
/musr/command globalfield printFieldValueAtPoint 0 0 0
/musr/command globalfield printFieldValueAtPoint 0 0 10
/musr/command globalfield printFieldValueAtPoint 0 0 20
/musr/command globalfield printFieldValueAtPoint 0 0 30
/musr/command globalfield printFieldValueAtPoint 0 0 50
/musr/command globalfield printFieldValueAtPoint 0 0 100
/musr/command globalfield printFieldValueAtPoint 0 0 200
/musr/command globalfield printFieldValueAtPoint 0 0 300
/musr/command globalfield printFieldValueAtPoint 0 0 400
/musr/command globalfield printFieldValueAtPoint 0 0 500
/musr/command globalfield printFieldValueAtPoint 10 0 0
/musr/command globalfield printFieldValueAtPoint 100 0 0
/musr/command globalfield printFieldValueAtPoint 300 0 0
/musr/command globalfield printFieldValueAtPoint 0 10 0
/musr/command globalfield printFieldValueAtPoint 0 100 0
/musr/command globalfield printFieldValueAtPoint 0 300 0
###################################################################################
######################### P H Y S I C S P R O C E S S E S ##################
###################################################################################
# Geant 4.9.3
/musr/command process addDiscreteProcess gamma G4PhotoElectricEffect
/musr/command process addDiscreteProcess gamma G4ComptonScattering
/musr/command process addDiscreteProcess gamma G4GammaConversion
/musr/command process addDiscreteProcess gamma G4RayleighScattering
/musr/command process addProcess e- G4eMultipleScattering -1 1 1
/musr/command process addProcess e- G4eIonisation -1 2 2
/musr/command process addProcess e- G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eMultipleScattering -1 1 1
/musr/command process addProcess e+ G4eIonisation -1 2 2
/musr/command process addProcess e+ G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eplusAnnihilation 0 -1 4
/musr/command process addProcess mu- G4MuMultipleScattering -1 1 1
/musr/command process addProcess mu- G4MuIonisation -1 2 2
/musr/command process addProcess mu- G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu- G4MuPairProduction -1 4 4
/musr/command process addProcess mu+ G4MuMultipleScattering -1 1 1
/musr/command process addProcess mu+ G4MuIonisation -1 2 2
/musr/command process addProcess mu+ G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu+ G4MuPairProduction -1 4 4
#
#/home/install/geant4.9.3/source/processes/electromagnetic/utils/src/G4EnergyLossMessenger.cc
######## /process/msc/StepLimit Minimal | UseDistanceToBoundary | UseSafety
/process/msc/StepLimit UseSafety
#/process/msc/LateralDisplacement
/process/msc/RangeFactor 0.04
/process/msc/GeomFactor 2.5
#/process/msc/FactorForAngleLimit
/process/msc/Skin 3.0
#/process/msc/ThetaLimit 0.2 rad
###################################################################################
################## S O M E O T H E R P A R A M E T E R S ##################
###################################################################################
# Set the overall range cut
#/run/setCut 1 mm
# Set the range cut on a particular volumes (in mm)
#/musr/command SetUserLimits log_target 0.01
#/musr/command SetUserLimits log_targetscint 0.01
#/musr/command SetUserLimits log_M0 0.01
# Store all events into the ROOT tree or just the interesting ones ? (true is default)
/musr/command storeOnlyEventsWithHits false
# Set the minimum time separation between two subsequent signals in the same detector (in ns)
/musr/command signalSeparationTime 10
#
#/musr/run/howOftenToPrintEvent 10
/musr/command maximumRunTimeAllowed 220000
/musr/run/randomOption 2
###################################################################################
######################### R O O T O U T P U T ##############################
###################################################################################
/musr/command rootOutput det_VvvKine off
/musr/command rootOutput det_VvvX off
/musr/command rootOutput det_VvvY off
/musr/command rootOutput det_VvvZ off
/musr/command rootOutput det_VvvVolID off
/musr/command rootOutput det_VvvProcID off
/musr/command rootOutput det_VvvTrackID off
/musr/command rootOutput det_VvvParticleID off
/musr/command rootOutput nOptPhot off
/musr/command rootOutput odet_ID off
/musr/command rootOutput odet_nPhot off
/musr/command rootOutput odet_timeFirst off
/musr/command rootOutput odet_timeA off
/musr/command rootOutput odet_timeB off
/musr/command rootOutput odet_timeC off
/musr/command rootOutput odet_timeD off
/musr/command rootOutput odet_timeE off
/musr/command rootOutput odet_timeLast off
###################################################################################
######################### V I S U A L I S A T I O N ##############################
###################################################################################
#/vis/disable
#/control/execute visFromToni.mac
#/control/execute visDawn40123.mac
/control/execute visVRML.mac
###################################################################################
######################### P A R T I C L E G U N #################################
###################################################################################
/gun/vertex 0 0 -1000 mm
# FWHM 10mm ==> sigma = 10/2.354 = 4.2481mm
/gun/vertexsigma 25 25 0 mm
#---/gun/vertexboundary: rMaxAllowed, zMinAllowed, zMaxAllowed
/gun/vertexboundary 100 -999999 999999 mm
/gun/momentum 100.0 MeV
# sigma = 3% ==> sigma 27*0.03 = 0.81
/gun/momentumsmearing 3.00 MeV
#---/gun/momentumboundary: pMinAllowed, pMaxAllowed, dummy
#/gun/momentumboundary 20 40 0 MeV
# TURTLE
#/gun/turtlefilename FOR061_2008_04_22.DAT
#/gun/turtlefilename FOR070_2008_10_17_XXII.DAT
#/gun/turtlefilename FOR077_pie3_HiField_d05_x30.dat
#/gun/turtlefilename FOR077_reggiani_Jan2010_NEW.dat
#/gun/turtleZ0position -900 mm
#/gun/turtlefilename FOR070_GPD_beam_profile_45mm.dat
#/gun/turtleZ0position -500 mm
#---/gun/tilt: xangle, yangle, dummy
#/gun/tilt 0 0.5 0 deg
#---/gun/tiltsigma: xangleSigma, yangleSigma, dummy (1 degree on 1 meter ~ 17mm)
/gun/tiltsigma 0.15 0.15 0 deg
/gun/pitch 0.08 deg
/gun/muonPolarizVector 0 0 -1
#---/gun/decaytimelimits: decayMin, decayMax, decayTime
/gun/decaytimelimits -1 -1 2197.03 ns
###################################################################################
######################## B E A M O N #########################################
###################################################################################
#/run/beamOn 16000000
/run/beamOn 1

174
run/GPD/40505f.v1190 Normal file
View File

@ -0,0 +1,174 @@
# TDC V1190 Set up file generated by Deltat : 16-Apr-2008 10:30:27
# bins are always in TDC channel units (195.3125 ps/channel)
INSTRUMENT=GPS
DESCRIPTION=No Veto -- 1port
TYPE=TDCV1190
RESOLUTION=100
MDELAY=0
PDELAY=2000
#REWINDTIMEBINS=10000000
MCOINCIDENCEW=50
PCOINCIDENCEW=50
VCOINCIDENCEW=100
# 30000 mu/s
MUONRATEFACTOR=0.0937185
DATAWINDOWMIN=-5
DATAWINDOWMAX=10.0
PILEUPWINDOWMIN=-10.2
PILEUPWINDOWMAX=10.2
102; "M up"; M; 0.4; 2005; ;
1; "B1"; P; 0.4; 2005; ; B1; 1485; 1515; 50995;
2; "B2"; P; 0.4; 2005; ; B2; 1485; 1515; 50995;
4; "B4"; P; 0.4; 2005; ; B4; 1485; 1515; 50995;
5; "B5"; P; 0.4; 2005; ; B5; 1485; 1515; 50995;
6; "B6"; P; 0.4; 2005; ; B6; 1485; 1515; 50995;
7; "B6"; P; 0.4; 2005; ; B7; 1485; 1515; 50995;
11; "F1"; P; 0.4; 2005; ; F11; 1485; 1515; 50995;
12; "F2"; P; 0.4; 2005; ; F12; 1485; 1515; 50995;
13; "F3"; P; 0.4; 2005; ; F13; 1485; 1515; 50995;
14; "F4"; P; 0.4; 2005; ; F14; 1485; 1515; 50995;
15; "F5"; P; 0.4; 2005; ; F15; 1485; 1515; 50995;
!103; "M embed" V; 0.1; 2005;
! 21; "Coinc B1" K; 0.4; 2005;
! 22; "Coinc B2" K; 0.4; 2005;
!401; "Act Veto" V; 0.1; 2005;
! 3; "F right"; K; 0.4; 1263;
! 4; "B left"; P; 0.2; 1311; 5 -10 -11 -12 -13; Back; 1485; 1515; 50995;
! 5; "B right"; K; 0.2; 1262;
! 6; Up; P; 0.2; 1306; -10 -11 -12 -13; Up; 1490; 1520; 51000;
! 7; Down; P; 0.2; 1264; -10 -11 -12 -13; Down; 1525; 1555; 51035;
! 8; "R int"; P; 0.2; 1309; -10 -11 -12 -13; R_int; 1480; 1510; 50990;
! 9; "R ext"; P; 0.3; 1267; -10 -11 -12 -13; R_ext; 1520; 1550; 51020; 8; Rite
!10; "Bveto up"; V; 0.3; 2085;
!11; "Bveto down"; V; 0.3; 2089;
!12; "Bveto left"; V; 0.3; 2080;
!13; "Bveto right"; V; 0.3; 2079;
!14; "F center"; P; 0.2; 2073; -10 -11 -12 -13; F_cntr; 720; 750; 50230; 2; Forw
!15; Clock; C;
!16; Rejected; R;
!17; Ch17; N;
!18; Ch18; N;
!19; Ch19; N;
!20; Ch20; N;
!21; Ch21; N;
!22; Ch22; N;
!23; Ch23; N;
!24; Ch24; N;
!25; Ch25; N;
!26; Ch26; N;
!27; Ch27; N;
!28; Ch28; N;
!29; Ch29; N;
!30; Ch30; N;
!31; Ch31; N;
#DEBUGEVENT 704 10
#DEBUGEVENT 1427 10
#DEBUGEVENT 3222 10
#DEBUGEVENT 5223 10
#DEBUGEVENT 8351 10
#DEBUGEVENT 8581 10
#DEBUGEVENT 9041 10
#DEBUGEVENT 9645 10
counterPhaseShifts 1 0 2 0 11 180 12 180 13 180
counterGrouping B 1 2
counterGrouping F 11 12 13
counterGrouping L 4 5
counterGrouping R 14 15
musrTH2D hMuDecayMap "Muon decay map;z [mm]; r[mm]" 100 -80. 20. 40 0. 20. muDecayPosZ muDecayPosR
musrTH2D hMuDecayMappileup "Pileup muon decay map;z [mm]; r[mm]" 400 -300. 100. 200 0. 200. pileup_muDecayPosZ pileup_muDecayPosR
musrTH1D hMuDecayPosZ "Penetration of muons into the sample;z[mm];N" 100 -5.0 5. muDecayPosZ
musrTH1D hMuDecayPosZpileup "Decay of the muons contributing to pileup bgr.;z[mm];N" 400 -300 100 pileup_muDecayPosZ
musrTH1D hMuDecayPosX "X of decayed muons;x[mm];N" 100 -5.0 5. muDecayPosX
musrTH1D hdet_m0edep "Energy deposited in the trigger;E[MeV];N" 100 0. 2 det_m0edep
musrTH1D hdet_posEdep "Energy deposited in the positron Counters;E[MeV];N" 100 0. 2 det_posEdep
#musrTH1D hdet_time0 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time0
#musrTH1D hdet_time1 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time1
#musrTH1D hdet_time10 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time10
musrTH1D hgen_time10 "generated muSR spectra;time[#mus];N" 220 -11. 11. gen_time10
musrTH1D hdet_time10_MINUS_gen_time10 "detected - generated ;time_{det}-time_{gen} [ps];N" 100 -500. 500. det_time10_MINUS_gen_time10
#musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 4000 -2000. 2000. muDecayDetID
musrTH1D hdet_time10 "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time10
musrTH1D hdet_time10_Bgr "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time10
musrTH1D hdet_time20 "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time20
musrTH1D hdet_time20_Bgr "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time20
musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 2002 -1001. 1001. muDecayDetID
musrTH1D hMuDecayDetIDpileup "Detector ID where the pileup muons decay;Detector ID; N" 2002 -1001. 1001. pileup_muDecayDetID
humanDecayHistograms hMuDecayDetID hMuDecayDetIDpileup 201 Sample 211 SampleCell 300 Holder 301 Holder 302 Holder 303 Holder 304 Holder 305 Holder 306 Holder 307 Holder 308 Holder 309 Holder 310 Holder 313 Holder 314 Holder 315 Holder 316 Holder 251 Col_Cu 255 Col_1 257 Col_1 261 Col_2 263 Col_2 265 Al_Plate 102 M0 103 M0Holder 1 BC 2 BC 11 FC 12 FC 13 FC 21 BBC 22 BBC 401 Veto -1 World 291 Magnet 292 Magnet 293 Al_Uprof 230 Elsewhere 231 Elsewhere 232 Elsewhere 233 Elsewhere 256 Elsewhere 262 Elsewhere 267 Elsewhere -1000 Escaped
condition 0 alwaysTrue
condition 1 oncePerEvent
condition 2 muonDecayedInSample_gen
condition 3 muonTriggered_gen
condition 4 muonTriggered_det
condition 5 positronHit_det
condition 6 goodEvent_det
condition 7 goodEvent_gen
condition 8 goodEvent_det_AND_goodEvent_gen
condition 9 pileupEvent
condition 10 goodEvent_F_det
condition 11 goodEvent_B_det
condition 12 goodEvent_F_det_AND_pileupEvent
condition 13 goodEvent_B_det_AND_pileupEvent
condition 14 goodEvent_L_det
condition 15 goodEvent_R_det
fit hdet_time20_Bgr pol0 111 -1.9 -0.1 0
#fit hdet_time20 funct1 111 0.2 9.8 0 0 0.3 100
#fit hdet_time20 funct2 111 0.2 9.8 0 0 0.3 100 0
fit hdet_time20 funct4 111 0.2 9.8 0 0 0.3 100 0
#draw hMuDecayPosZ 0
#draw hMuDecayPosZ 1
#draw hMuDecayPosZ 2
#draw hMuDecayPosZ 3
#draw hMuDecayPosZ 4
###draw hMuDecayMap 0
###draw hMuDecayMap 1
#draw hdet_m0edep 0
#draw hdet_m0edep 1
#draw hdet_m0edep 2
#draw hdet_m0edep 3
#draw hdet_m0edep 4
#draw hdet_m0edep 5
#draw hdet_m0edep 6
#draw hdet_time10 6
draw hdet_time20 10
draw hdet_time20 11
draw hdet_time20 14
draw hdet_time20 15
draw hdet_time20_Bgr 10
draw hdet_time20_Bgr 11
draw hdet_time20_Bgr 14
draw hdet_time20_Bgr 15
#draw hgen_time10 7
#draw hdet_time10_MINUS_gen_time10 8
#draw hdet_posEdep 5
#draw hMuDecayDetID 1
#draw hdet_m0edep 9
#draw hdet_time10 9
#draw hMuDecayDetIDpileup 1
#draw hMuDecayDetIDpileup 6
#draw hMuDecayDetIDpileup 9
#draw humanDecayHistograms 1
#draw humanDecayHistograms 6
#draw humanDecayHistograms 9
#draw humanDecayPileupHistograms 1
#draw humanDecayPileupHistograms 6
#draw humanDecayPileupHistograms 9
#draw hdet_time20 4
#draw hdet_time20 6
$

174
run/GPD/40505k.v1190 Normal file
View File

@ -0,0 +1,174 @@
# TDC V1190 Set up file generated by Deltat : 16-Apr-2008 10:30:27
# bins are always in TDC channel units (195.3125 ps/channel)
INSTRUMENT=GPS
DESCRIPTION=No Veto -- 1port
TYPE=TDCV1190
RESOLUTION=100
MDELAY=0
PDELAY=2000
#REWINDTIMEBINS=10000000
MCOINCIDENCEW=50
PCOINCIDENCEW=50
VCOINCIDENCEW=100
# 30000 mu/s
MUONRATEFACTOR=0.0937185
DATAWINDOWMIN=-5
DATAWINDOWMAX=10.0
PILEUPWINDOWMIN=-10.2
PILEUPWINDOWMAX=10.2
102; "M up"; M; 0.4; 2005; 103;
1; "B1"; P; 0.4; 2005; ; B1; 1485; 1515; 50995;
2; "B2"; P; 0.4; 2005; ; B2; 1485; 1515; 50995;
4; "B4"; P; 0.4; 2005; ; B4; 1485; 1515; 50995;
5; "B5"; P; 0.4; 2005; ; B5; 1485; 1515; 50995;
6; "B6"; P; 0.4; 2005; ; B6; 1485; 1515; 50995;
7; "B6"; P; 0.4; 2005; ; B7; 1485; 1515; 50995;
11; "F1"; P; 0.4; 2005; ; F11; 1485; 1515; 50995;
12; "F2"; P; 0.4; 2005; ; F12; 1485; 1515; 50995;
13; "F3"; P; 0.4; 2005; ; F13; 1485; 1515; 50995;
14; "F4"; P; 0.4; 2005; ; F14; 1485; 1515; 50995;
15; "F5"; P; 0.4; 2005; ; F15; 1485; 1515; 50995;
103; "M embed" K; 0.1; 2005;
! 21; "Coinc B1" K; 0.4; 2005;
! 22; "Coinc B2" K; 0.4; 2005;
!401; "Act Veto" V; 0.1; 2005;
! 3; "F right"; K; 0.4; 1263;
! 4; "B left"; P; 0.2; 1311; 5 -10 -11 -12 -13; Back; 1485; 1515; 50995;
! 5; "B right"; K; 0.2; 1262;
! 6; Up; P; 0.2; 1306; -10 -11 -12 -13; Up; 1490; 1520; 51000;
! 7; Down; P; 0.2; 1264; -10 -11 -12 -13; Down; 1525; 1555; 51035;
! 8; "R int"; P; 0.2; 1309; -10 -11 -12 -13; R_int; 1480; 1510; 50990;
! 9; "R ext"; P; 0.3; 1267; -10 -11 -12 -13; R_ext; 1520; 1550; 51020; 8; Rite
!10; "Bveto up"; V; 0.3; 2085;
!11; "Bveto down"; V; 0.3; 2089;
!12; "Bveto left"; V; 0.3; 2080;
!13; "Bveto right"; V; 0.3; 2079;
!14; "F center"; P; 0.2; 2073; -10 -11 -12 -13; F_cntr; 720; 750; 50230; 2; Forw
!15; Clock; C;
!16; Rejected; R;
!17; Ch17; N;
!18; Ch18; N;
!19; Ch19; N;
!20; Ch20; N;
!21; Ch21; N;
!22; Ch22; N;
!23; Ch23; N;
!24; Ch24; N;
!25; Ch25; N;
!26; Ch26; N;
!27; Ch27; N;
!28; Ch28; N;
!29; Ch29; N;
!30; Ch30; N;
!31; Ch31; N;
#DEBUGEVENT 704 10
#DEBUGEVENT 1427 10
#DEBUGEVENT 3222 10
#DEBUGEVENT 5223 10
#DEBUGEVENT 8351 10
#DEBUGEVENT 8581 10
#DEBUGEVENT 9041 10
#DEBUGEVENT 9645 10
counterPhaseShifts 1 0 2 0 11 180 12 180 13 180
counterGrouping B 1 2
counterGrouping F 11 12 13
counterGrouping L 4 5
counterGrouping R 14 15
musrTH2D hMuDecayMap "Muon decay map;z [mm]; r[mm]" 100 -80. 20. 40 0. 20. muDecayPosZ muDecayPosR
musrTH2D hMuDecayMappileup "Pileup muon decay map;z [mm]; r[mm]" 400 -300. 100. 200 0. 200. pileup_muDecayPosZ pileup_muDecayPosR
musrTH1D hMuDecayPosZ "Penetration of muons into the sample;z[mm];N" 100 -5.0 5. muDecayPosZ
musrTH1D hMuDecayPosZpileup "Decay of the muons contributing to pileup bgr.;z[mm];N" 400 -300 100 pileup_muDecayPosZ
musrTH1D hMuDecayPosX "X of decayed muons;x[mm];N" 100 -5.0 5. muDecayPosX
musrTH1D hdet_m0edep "Energy deposited in the trigger;E[MeV];N" 100 0. 2 det_m0edep
musrTH1D hdet_posEdep "Energy deposited in the positron Counters;E[MeV];N" 100 0. 2 det_posEdep
#musrTH1D hdet_time0 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time0
#musrTH1D hdet_time1 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time1
#musrTH1D hdet_time10 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time10
musrTH1D hgen_time10 "generated muSR spectra;time[#mus];N" 220 -11. 11. gen_time10
musrTH1D hdet_time10_MINUS_gen_time10 "detected - generated ;time_{det}-time_{gen} [ps];N" 100 -500. 500. det_time10_MINUS_gen_time10
#musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 4000 -2000. 2000. muDecayDetID
musrTH1D hdet_time10 "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time10
musrTH1D hdet_time10_Bgr "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time10
musrTH1D hdet_time20 "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time20
musrTH1D hdet_time20_Bgr "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time20
musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 2002 -1001. 1001. muDecayDetID
musrTH1D hMuDecayDetIDpileup "Detector ID where the pileup muons decay;Detector ID; N" 2002 -1001. 1001. pileup_muDecayDetID
humanDecayHistograms hMuDecayDetID hMuDecayDetIDpileup 201 Sample 211 SampleCell 300 Holder 301 Holder 302 Holder 303 Holder 304 Holder 305 Holder 306 Holder 307 Holder 308 Holder 309 Holder 310 Holder 313 Holder 314 Holder 315 Holder 316 Holder 251 Col_Cu 255 Col_1 257 Col_1 261 Col_2 263 Col_2 265 Al_Plate 102 M0 103 M0Holder 1 BC 2 BC 11 FC 12 FC 13 FC 21 BBC 22 BBC 401 Veto -1 World 291 Magnet 292 Magnet 293 Al_Uprof 230 Elsewhere 231 Elsewhere 232 Elsewhere 233 Elsewhere 256 Elsewhere 262 Elsewhere 267 Elsewhere -1000 Escaped
condition 0 alwaysTrue
condition 1 oncePerEvent
condition 2 muonDecayedInSample_gen
condition 3 muonTriggered_gen
condition 4 muonTriggered_det
condition 5 positronHit_det
condition 6 goodEvent_det
condition 7 goodEvent_gen
condition 8 goodEvent_det_AND_goodEvent_gen
condition 9 pileupEvent
condition 10 goodEvent_F_det
condition 11 goodEvent_B_det
condition 12 goodEvent_F_det_AND_pileupEvent
condition 13 goodEvent_B_det_AND_pileupEvent
condition 14 goodEvent_L_det
condition 15 goodEvent_R_det
fit hdet_time20_Bgr pol0 111 -1.9 -0.1 0
#fit hdet_time20 funct1 111 0.2 9.8 0 0 0.3 100
#fit hdet_time20 funct2 111 0.2 9.8 0 0 0.3 100 0
fit hdet_time20 funct4 111 0.2 9.8 0 0 0.3 100 0
#draw hMuDecayPosZ 0
draw hMuDecayPosZ 1
#draw hMuDecayPosZ 2
#draw hMuDecayPosZ 3
#draw hMuDecayPosZ 4
###draw hMuDecayMap 0
###draw hMuDecayMap 1
draw hdet_m0edep 0
draw hdet_m0edep 1
draw hdet_m0edep 2
draw hdet_m0edep 3
draw hdet_m0edep 4
draw hdet_m0edep 5
draw hdet_m0edep 6
draw hdet_time10 6
draw hdet_time20 10
draw hdet_time20 11
draw hdet_time20 14
draw hdet_time20 15
draw hdet_time20_Bgr 10
draw hdet_time20_Bgr 11
draw hdet_time20_Bgr 14
draw hdet_time20_Bgr 15
draw hgen_time10 7
draw hdet_time10_MINUS_gen_time10 8
draw hdet_posEdep 5
draw hMuDecayDetID 1
draw hdet_m0edep 9
draw hdet_time10 9
draw hMuDecayDetIDpileup 1
draw hMuDecayDetIDpileup 6
draw hMuDecayDetIDpileup 9
draw humanDecayHistograms 1
draw humanDecayHistograms 6
draw humanDecayHistograms 9
draw humanDecayPileupHistograms 1
draw humanDecayPileupHistograms 6
draw humanDecayPileupHistograms 9
draw hdet_time20 4
draw hdet_time20 6
$

260
run/GPD/40506.mac Normal file
View File

@ -0,0 +1,260 @@
#-----------------------------------------------------------------------
# Macro file for muSR instrument simulations
# Unless specified otherwises, the default units are mm, ns, MeV, MeV/c.
# Lines starting with star "#" are comments.
###################################################################################
############################# G E O M E T R Y ###################################
###################################################################################
# ROTATION MATRIXES:
/musr/command rotation matrix1 1 0 0 90
/musr/command rotation matrix2 0 1 0 28.81
/musr/command rotation matrix3 0 1 0 151.19
/musr/command rotation matrix4 0 0 1 45
/musr/command rotation matrix5 0 1 0 90
/musr/command rotation matrix6 180 90 90
/musr/command rotation matrix7 0 90 90
# --- WORLD VOLUME ---
/musr/command construct box World 2100 2100 11000 G4_AIR 0 0 0 no_logical_volume norot dead -1
# --- BEAMPIPE ---
/musr/command construct tubs beampipeA 0 2005 4500 0 360 G4_Galactic 0 0 -6100 log_World norot dead 230
/musr/command construct tubs beampipe 0 405 500 0 360 G4_Galactic 0 0 -1100 log_World norot dead 231
/musr/command construct tubs beampipeshell 100 105 500 0 360 G4_Al 0 0 0 log_beampipe norot dead 232
/musr/command construct tubs beampipewindow 0 100 0.05 0 360 G4_MYLAR 0 0 499.9 log_beampipe norot dead 233
# PASSIVE COLLIMATOR (for BEAMPIPE vers.1)
# GPDcollimator - thickness in the thinner location is x2+(x9-x4) = 20+59-59
##/musr/command construct GPDcollimator GPDcollimator 40 15 60 59 60.01 0 360 0 59 0 2.5 6 G4_Cu 0 0 -59 log_World matrix1 dead 251
/musr/command construct box GPDcollimator2 48 55 15 G4_Pb 0 0 -66 log_World norot dead 255
/musr/command construct box GPDcollimator2hole 2.5 4 15 G4_Galactic 0 0 0 log_GPDcollimator2 norot dead 256
/musr/command construct box FirstCollimator 45 75 15 G4_Pb 0 0 -250 log_World norot dead 261
/musr/command construct tubs FirstCollHole 0 8 15 0 360 G4_Galactic 0 0 0 log_FirstCollimator norot dead 262
###/musr/command construct box AlumPlate 30 90 5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -80 log_World norot musr/ScintSD 401
###/musr/command construct box AlumPlateHole 2 5 5 G4_Galactic 0 0 0 log_AlumPlate norot dead 267
/musr/command construct box Insert 20 20 10 G4_Pb 0 0 -91 log_World norot dead 265
/musr/command construct box InsertHole 2.5 4 10 G4_Galactic 0 0 0 log_Insert norot dead 267
# --- SAMPLE HOLDER ---
##/musr/command construct GPDsampleHolderA GPDsampleHolderA 0 0 0 0 0 G4_Al 0 -111.25 0 log_World matrix1 dead 300
##/musr/command construct tubs GPDsampleHolderB 37.51 58 7 0 360 G4_Al 0 57 0 log_World matrix1 dead 301
##/musr/command construct tubs GPDsampleHolderC 37.51 58 7 0 360 G4_Al 0 -57 0 log_World matrix1 dead 302
######/musr/command construct tubs GPDsampleHolderD 37 37.5 50 136 88 G4_Al 0 0 0 log_World matrix1 dead 303
######/musr/command construct tubs GPDsampleHolderE 53.5 54 50 136 88 G4_Al 0 0 0 log_World matrix1 dead 304
##/musr/command construct box GPDsampleHolderD 12 0.5 50 G4_Al 38 0 28 log_World matrix1 dead 303
##/musr/command construct box GPDsampleHolderE 12 0.5 50 G4_Al 38 0 -28 log_World matrix1 dead 304
##/musr/command construct box GPDsampleHolderF 12 0.5 50 G4_Al -38 0 28 log_World matrix1 dead 305
##/musr/command construct box GPDsampleHolderG 12 0.5 50 G4_Al -38 0 -28 log_World matrix1 dead 306
##/musr/command construct tubs GPDsampleHolderH 0 9 50 180 50 G4_Al 48 0 29 log_World matrix1 dead 307
##/musr/command construct tubs GPDsampleHolderI 0 9 50 130 50 G4_Al 48 0 -29 log_World matrix1 dead 308
##/musr/command construct tubs GPDsampleHolderJ 0 9 50 310 50 G4_Al -48 0 29 log_World matrix1 dead 309
##/musr/command construct tubs GPDsampleHolderK 0 9 50 0 50 G4_Al -48 0 -29 log_World matrix1 dead 310
##/musr/command construct tubs GPDsampleHolderW 58.01 58.2 64 25 130 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 313
##/musr/command construct tubs GPDsampleHolderX 58.01 58.2 64 205 130 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 314
##/musr/command construct tubs GPDsampleHolderY 36.79 36.99 64 45 90 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 315
##/musr/command construct tubs GPDsampleHolderZ 36.79 36.99 64 225 90 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 316
# --- SAMPLE CELL ---
/musr/command construct tubs GPDsampleCell 0 12 50 0 360 G4_Cu 0 0 0 log_World matrix1 dead 211
# --- SAMPLE ---
/musr/command construct tubs target 0 3 5 0 360 G4_Cu 0 0 0 log_GPDsampleCell norot dead 201
#/musr/command construct tubs targetFieldVol 0 3 5 0 360 G4_Cu 0 0 0 log_GPDsampleCell norot dead 201
# --- SCINTILLATOR ---
/musr/command construct box ScintB1 21.5 50 2 G4_PLASTIC_SC_VINYLTOLUENE 26.5 0 -44 log_World norot musr/ScintSD 1
/musr/command construct box ScintB2 21.5 50 2 G4_PLASTIC_SC_VINYLTOLUENE -26.5 0 -44 log_World norot musr/ScintSD 2
/musr/command construct box ScintB4 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE 44 0 -25 log_World norot musr/ScintSD 4
/musr/command construct box ScintB5 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE -44 0 -25 log_World norot musr/ScintSD 5
/musr/command construct box ScintB6 2 50 13 G4_PLASTIC_SC_VINYLTOLUENE 44 0 0 log_World norot musr/ScintSD 6
/musr/command construct box ScintB7 2 50 13 G4_PLASTIC_SC_VINYLTOLUENE -44 0 0 log_World norot musr/ScintSD 7
/musr/command construct box ScintF1 20 50 2 G4_PLASTIC_SC_VINYLTOLUENE 28 0 44 log_World norot musr/ScintSD 11
/musr/command construct box ScintF2 20 50 2 G4_PLASTIC_SC_VINYLTOLUENE -28 0 44 log_World norot musr/ScintSD 12
/musr/command construct box ScintF3 7 50 2 G4_PLASTIC_SC_VINYLTOLUENE 0 0 44 log_World norot musr/ScintSD 13
/musr/command construct box ScintF4 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE 44 0 25 log_World norot musr/ScintSD 14
/musr/command construct box ScintF5 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE -44 0 25 log_World norot musr/ScintSD 15
#/musr/command construct box KoincB1 20 90 2 G4_PLASTIC_SC_VINYLTOLUENE 33 0 -60 log_World matrix2 musr/ScintSD 21
#/musr/command construct box KoincB2 20 90 2 G4_PLASTIC_SC_VINYLTOLUENE -33 0 -60 log_World matrix3 musr/ScintSD 22
# --- M0 ---
/musr/command construct box M0holder 2.5 4 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -20 log_World norot musr/ScintSD 103
/musr/command construct box M0 4 30 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -44 log_World norot musr/ScintSD 102
#/musr/command construct GPDmHolder GPDmHolder 5 25 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 -25 -43 log_World norot dead 103
# --- Magnet ---
/musr/command construct tubs Magnet1 180 375 78 0 360 G4_Cu 133 0 0 log_World matrix5 dead 291
/musr/command construct tubs Magnet2 180 375 78 0 360 G4_Cu -133 0 0 log_World matrix5 dead 292
/musr/command construct uprofile Uprof1 25 10 400 5 G4_Al -45 0 200 log_World matrix6 dead 293
/musr/command construct uprofile Uprof2 25 10 400 5 G4_Al -45 0 -200 log_World matrix6 dead 293
/musr/command construct uprofile Uprof3 25 10 400 5 G4_Al 45 0 200 log_World matrix7 dead 293
/musr/command construct uprofile Uprof4 25 10 400 5 G4_Al 45 0 -200 log_World matrix7 dead 293
#
##---------------------------------------------------------
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_1 6
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_2 5
#/musr/command logicalVolumeToBeReweighted mu log_activeCollim 4
##---------------------------------------------------------
/musr/command visattributes log_World invisible
/musr/command visattributes log_beampipeA invisible
/musr/command visattributes log_beampipe invisible
/musr/command visattributes log_beampipeshell invisible
/musr/command visattributes log_beampipewindow invisible
/musr/command visattributes log_Magnet1 invisible
/musr/command visattributes log_Magnet2 invisible
/musr/command visattributes log_Uprof1 invisible
/musr/command visattributes log_Uprof2 invisible
/musr/command visattributes log_Uprof3 invisible
/musr/command visattributes log_Uprof4 invisible
/musr/command visattributes log_GPDsampleHolderD red
/musr/command visattributes log_GPDsampleHolderE red
/musr/command visattributes log_GPDsampleHolderF red
/musr/command visattributes log_GPDsampleHolderG red
/musr/command visattributes log_GPDcollimator green
/musr/command visattributes log_GPDcollimator2 blue
/musr/command visattributes log_GPDmHolder magenta
/musr/command visattributes log_target blue
/musr/command visattributes log_M0 yellow
#/musr/command visattributes G4_PLASTIC_SC_VINYLTOLUENE red
#/musr/command visattributes G4_Pb blue
####################################################################################
######################### M A G N E T I C F I E L D #########################
###################################################################################
# Set magnetic field (set field intensity in T and sigma in mm)
# syntax for magneticfield: fromfile filename fieldValue
# uniform fieldValue
# gaussian fieldValue sigma
#
/musr/command globalfield centralSolenoidField 0. 0. 0. fromfile 2D GPD_coil_primitive.table log_target 0.03
/musr/command globalfield setparameter SetLargestAcceptableStep 200
/musr/command globalfield setparameter SetMinimumEpsilonStep 5e-5
/musr/command globalfield setparameter SetMaximumEpsilonStep 0.001
/musr/command globalfield setparameter SetDeltaOneStep 0.1
/musr/command globalfield setparameter SetDeltaIntersection 0.01
/musr/command globalfield printparameters
/musr/command globalfield printFieldValueAtPoint 0 0 0
/musr/command globalfield printFieldValueAtPoint 0 0 10
/musr/command globalfield printFieldValueAtPoint 0 0 20
/musr/command globalfield printFieldValueAtPoint 0 0 30
/musr/command globalfield printFieldValueAtPoint 0 0 50
/musr/command globalfield printFieldValueAtPoint 0 0 100
/musr/command globalfield printFieldValueAtPoint 0 0 200
/musr/command globalfield printFieldValueAtPoint 0 0 300
/musr/command globalfield printFieldValueAtPoint 0 0 400
/musr/command globalfield printFieldValueAtPoint 0 0 500
/musr/command globalfield printFieldValueAtPoint 10 0 0
/musr/command globalfield printFieldValueAtPoint 100 0 0
/musr/command globalfield printFieldValueAtPoint 300 0 0
/musr/command globalfield printFieldValueAtPoint 0 10 0
/musr/command globalfield printFieldValueAtPoint 0 100 0
/musr/command globalfield printFieldValueAtPoint 0 300 0
###################################################################################
######################### P H Y S I C S P R O C E S S E S ##################
###################################################################################
# Geant 4.9.3
/musr/command process addDiscreteProcess gamma G4PhotoElectricEffect
/musr/command process addDiscreteProcess gamma G4ComptonScattering
/musr/command process addDiscreteProcess gamma G4GammaConversion
/musr/command process addDiscreteProcess gamma G4RayleighScattering
/musr/command process addProcess e- G4eMultipleScattering -1 1 1
/musr/command process addProcess e- G4eIonisation -1 2 2
/musr/command process addProcess e- G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eMultipleScattering -1 1 1
/musr/command process addProcess e+ G4eIonisation -1 2 2
/musr/command process addProcess e+ G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eplusAnnihilation 0 -1 4
/musr/command process addProcess mu- G4MuMultipleScattering -1 1 1
/musr/command process addProcess mu- G4MuIonisation -1 2 2
/musr/command process addProcess mu- G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu- G4MuPairProduction -1 4 4
/musr/command process addProcess mu+ G4MuMultipleScattering -1 1 1
/musr/command process addProcess mu+ G4MuIonisation -1 2 2
/musr/command process addProcess mu+ G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu+ G4MuPairProduction -1 4 4
#
#/home/install/geant4.9.3/source/processes/electromagnetic/utils/src/G4EnergyLossMessenger.cc
######## /process/msc/StepLimit Minimal | UseDistanceToBoundary | UseSafety
/process/msc/StepLimit UseSafety
#/process/msc/LateralDisplacement
/process/msc/RangeFactor 0.04
/process/msc/GeomFactor 2.5
#/process/msc/FactorForAngleLimit
/process/msc/Skin 3.0
#/process/msc/ThetaLimit 0.2 rad
###################################################################################
################## S O M E O T H E R P A R A M E T E R S ##################
###################################################################################
# Set the overall range cut
#/run/setCut 1 mm
# Set the range cut on a particular volumes (in mm)
#/musr/command SetUserLimits log_target 0.01
#/musr/command SetUserLimits log_targetscint 0.01
#/musr/command SetUserLimits log_M0 0.01
# Store all events into the ROOT tree or just the interesting ones ? (true is default)
/musr/command storeOnlyEventsWithHits false
# Set the minimum time separation between two subsequent signals in the same detector (in ns)
/musr/command signalSeparationTime 10
#
#/musr/run/howOftenToPrintEvent 10
/musr/command maximumRunTimeAllowed 220000
/musr/run/randomOption 2
###################################################################################
######################### R O O T O U T P U T ##############################
###################################################################################
/musr/command rootOutput det_VvvKine off
/musr/command rootOutput det_VvvX off
/musr/command rootOutput det_VvvY off
/musr/command rootOutput det_VvvZ off
/musr/command rootOutput det_VvvVolID off
/musr/command rootOutput det_VvvProcID off
/musr/command rootOutput det_VvvTrackID off
/musr/command rootOutput det_VvvParticleID off
/musr/command rootOutput nOptPhot off
/musr/command rootOutput odet_ID off
/musr/command rootOutput odet_nPhot off
/musr/command rootOutput odet_timeFirst off
/musr/command rootOutput odet_timeA off
/musr/command rootOutput odet_timeB off
/musr/command rootOutput odet_timeC off
/musr/command rootOutput odet_timeD off
/musr/command rootOutput odet_timeE off
/musr/command rootOutput odet_timeLast off
###################################################################################
######################### V I S U A L I S A T I O N ##############################
###################################################################################
/vis/disable
#/control/execute visFromToni.mac
#/control/execute visDawn40123.mac
#/control/execute visVRML.mac
###################################################################################
######################### P A R T I C L E G U N #################################
###################################################################################
/gun/vertex 0 0 -1000 mm
# FWHM 10mm ==> sigma = 10/2.354 = 4.2481mm
/gun/vertexsigma 25 25 0 mm
#---/gun/vertexboundary: rMaxAllowed, zMinAllowed, zMaxAllowed
/gun/vertexboundary 100 -999999 999999 mm
/gun/momentum 100.0 MeV
# sigma = 3% ==> sigma 27*0.03 = 0.81
/gun/momentumsmearing 3.00 MeV
#---/gun/momentumboundary: pMinAllowed, pMaxAllowed, dummy
#/gun/momentumboundary 20 40 0 MeV
# TURTLE
#/gun/turtlefilename FOR061_2008_04_22.DAT
#/gun/turtlefilename FOR070_2008_10_17_XXII.DAT
#/gun/turtlefilename FOR077_pie3_HiField_d05_x30.dat
#/gun/turtlefilename FOR077_reggiani_Jan2010_NEW.dat
#/gun/turtleZ0position -900 mm
#/gun/turtlefilename FOR070_GPD_beam_profile_45mm.dat
#/gun/turtleZ0position -500 mm
#---/gun/tilt: xangle, yangle, dummy
#/gun/tilt 0 0.5 0 deg
#---/gun/tiltsigma: xangleSigma, yangleSigma, dummy (1 degree on 1 meter ~ 17mm)
/gun/tiltsigma 0.15 0.15 0 deg
/gun/pitch 0.08 deg
/gun/muonPolarizVector 0 0 -1
#---/gun/decaytimelimits: decayMin, decayMax, decayTime
/gun/decaytimelimits -1 -1 2197.03 ns
###################################################################################
######################## B E A M O N #########################################
###################################################################################
/run/beamOn 16000000
#/run/beamOn 1

260
run/GPD/40507.mac Normal file
View File

@ -0,0 +1,260 @@
#-----------------------------------------------------------------------
# Macro file for muSR instrument simulations
# Unless specified otherwises, the default units are mm, ns, MeV, MeV/c.
# Lines starting with star "#" are comments.
###################################################################################
############################# G E O M E T R Y ###################################
###################################################################################
# ROTATION MATRIXES:
/musr/command rotation matrix1 1 0 0 90
/musr/command rotation matrix2 0 1 0 28.81
/musr/command rotation matrix3 0 1 0 151.19
/musr/command rotation matrix4 0 0 1 45
/musr/command rotation matrix5 0 1 0 90
/musr/command rotation matrix6 180 90 90
/musr/command rotation matrix7 0 90 90
# --- WORLD VOLUME ---
/musr/command construct box World 2100 2100 11000 G4_AIR 0 0 0 no_logical_volume norot dead -1
# --- BEAMPIPE ---
/musr/command construct tubs beampipeA 0 2005 4500 0 360 G4_Galactic 0 0 -6100 log_World norot dead 230
/musr/command construct tubs beampipe 0 405 500 0 360 G4_Galactic 0 0 -1100 log_World norot dead 231
/musr/command construct tubs beampipeshell 100 105 500 0 360 G4_Al 0 0 0 log_beampipe norot dead 232
/musr/command construct tubs beampipewindow 0 100 0.05 0 360 G4_MYLAR 0 0 499.9 log_beampipe norot dead 233
# PASSIVE COLLIMATOR (for BEAMPIPE vers.1)
# GPDcollimator - thickness in the thinner location is x2+(x9-x4) = 20+59-59
##/musr/command construct GPDcollimator GPDcollimator 40 15 60 59 60.01 0 360 0 59 0 2.5 6 G4_Cu 0 0 -59 log_World matrix1 dead 251
/musr/command construct box GPDcollimator2 48 55 15 G4_Pb 0 0 -66 log_World norot dead 255
/musr/command construct box GPDcollimator2hole 2.5 4 15 G4_Galactic 0 0 0 log_GPDcollimator2 norot dead 256
/musr/command construct box FirstCollimator 45 75 15 G4_Pb 0 0 -250 log_World norot dead 261
/musr/command construct tubs FirstCollHole 0 8 15 0 360 G4_Galactic 0 0 0 log_FirstCollimator norot dead 262
###/musr/command construct box AlumPlate 30 90 5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -80 log_World norot musr/ScintSD 401
###/musr/command construct box AlumPlateHole 2 5 5 G4_Galactic 0 0 0 log_AlumPlate norot dead 267
/musr/command construct box Insert 20 20 10 G4_Pb 0 0 -91 log_World norot dead 265
/musr/command construct box InsertHole 2.5 4 10 G4_Galactic 0 0 0 log_Insert norot dead 267
# --- SAMPLE HOLDER ---
##/musr/command construct GPDsampleHolderA GPDsampleHolderA 0 0 0 0 0 G4_Al 0 -111.25 0 log_World matrix1 dead 300
##/musr/command construct tubs GPDsampleHolderB 37.51 58 7 0 360 G4_Al 0 57 0 log_World matrix1 dead 301
##/musr/command construct tubs GPDsampleHolderC 37.51 58 7 0 360 G4_Al 0 -57 0 log_World matrix1 dead 302
######/musr/command construct tubs GPDsampleHolderD 37 37.5 50 136 88 G4_Al 0 0 0 log_World matrix1 dead 303
######/musr/command construct tubs GPDsampleHolderE 53.5 54 50 136 88 G4_Al 0 0 0 log_World matrix1 dead 304
##/musr/command construct box GPDsampleHolderD 12 0.5 50 G4_Al 38 0 28 log_World matrix1 dead 303
##/musr/command construct box GPDsampleHolderE 12 0.5 50 G4_Al 38 0 -28 log_World matrix1 dead 304
##/musr/command construct box GPDsampleHolderF 12 0.5 50 G4_Al -38 0 28 log_World matrix1 dead 305
##/musr/command construct box GPDsampleHolderG 12 0.5 50 G4_Al -38 0 -28 log_World matrix1 dead 306
##/musr/command construct tubs GPDsampleHolderH 0 9 50 180 50 G4_Al 48 0 29 log_World matrix1 dead 307
##/musr/command construct tubs GPDsampleHolderI 0 9 50 130 50 G4_Al 48 0 -29 log_World matrix1 dead 308
##/musr/command construct tubs GPDsampleHolderJ 0 9 50 310 50 G4_Al -48 0 29 log_World matrix1 dead 309
##/musr/command construct tubs GPDsampleHolderK 0 9 50 0 50 G4_Al -48 0 -29 log_World matrix1 dead 310
##/musr/command construct tubs GPDsampleHolderW 58.01 58.2 64 25 130 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 313
##/musr/command construct tubs GPDsampleHolderX 58.01 58.2 64 205 130 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 314
##/musr/command construct tubs GPDsampleHolderY 36.79 36.99 64 45 90 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 315
##/musr/command construct tubs GPDsampleHolderZ 36.79 36.99 64 225 90 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 316
# --- SAMPLE CELL ---
/musr/command construct tubs GPDsampleCell 0 12 50 0 360 G4_Cu 0 0 0 log_World matrix1 dead 211
# --- SAMPLE ---
/musr/command construct tubs target 0 3 5 0 360 G4_Cu 0 0 0 log_GPDsampleCell norot dead 201
#/musr/command construct tubs targetFieldVol 0 3 5 0 360 G4_Cu 0 0 0 log_GPDsampleCell norot dead 201
# --- SCINTILLATOR ---
/musr/command construct box ScintB1 21.5 50 2 G4_PLASTIC_SC_VINYLTOLUENE 26.5 0 -44 log_World norot musr/ScintSD 1
/musr/command construct box ScintB2 21.5 50 2 G4_PLASTIC_SC_VINYLTOLUENE -26.5 0 -44 log_World norot musr/ScintSD 2
/musr/command construct box ScintB4 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE 44 0 -25 log_World norot musr/ScintSD 4
/musr/command construct box ScintB5 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE -44 0 -25 log_World norot musr/ScintSD 5
/musr/command construct box ScintB6 2 50 13 G4_PLASTIC_SC_VINYLTOLUENE 44 0 0 log_World norot musr/ScintSD 6
/musr/command construct box ScintB7 2 50 13 G4_PLASTIC_SC_VINYLTOLUENE -44 0 0 log_World norot musr/ScintSD 7
/musr/command construct box ScintF1 20 50 2 G4_PLASTIC_SC_VINYLTOLUENE 28 0 44 log_World norot musr/ScintSD 11
/musr/command construct box ScintF2 20 50 2 G4_PLASTIC_SC_VINYLTOLUENE -28 0 44 log_World norot musr/ScintSD 12
/musr/command construct box ScintF3 7 50 2 G4_PLASTIC_SC_VINYLTOLUENE 0 0 44 log_World norot musr/ScintSD 13
/musr/command construct box ScintF4 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE 44 0 25 log_World norot musr/ScintSD 14
/musr/command construct box ScintF5 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE -44 0 25 log_World norot musr/ScintSD 15
#/musr/command construct box KoincB1 20 90 2 G4_PLASTIC_SC_VINYLTOLUENE 33 0 -60 log_World matrix2 musr/ScintSD 21
#/musr/command construct box KoincB2 20 90 2 G4_PLASTIC_SC_VINYLTOLUENE -33 0 -60 log_World matrix3 musr/ScintSD 22
# --- M0 ---
/musr/command construct box M0holder 2.5 4 0.5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -20 log_World norot musr/ScintSD 103
/musr/command construct box M0 4 30 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -44 log_World norot musr/ScintSD 102
#/musr/command construct GPDmHolder GPDmHolder 5 25 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 -25 -43 log_World norot dead 103
# --- Magnet ---
/musr/command construct tubs Magnet1 180 375 78 0 360 G4_Cu 133 0 0 log_World matrix5 dead 291
/musr/command construct tubs Magnet2 180 375 78 0 360 G4_Cu -133 0 0 log_World matrix5 dead 292
/musr/command construct uprofile Uprof1 25 10 400 5 G4_Al -45 0 200 log_World matrix6 dead 293
/musr/command construct uprofile Uprof2 25 10 400 5 G4_Al -45 0 -200 log_World matrix6 dead 293
/musr/command construct uprofile Uprof3 25 10 400 5 G4_Al 45 0 200 log_World matrix7 dead 293
/musr/command construct uprofile Uprof4 25 10 400 5 G4_Al 45 0 -200 log_World matrix7 dead 293
#
##---------------------------------------------------------
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_1 6
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_2 5
#/musr/command logicalVolumeToBeReweighted mu log_activeCollim 4
##---------------------------------------------------------
/musr/command visattributes log_World invisible
/musr/command visattributes log_beampipeA invisible
/musr/command visattributes log_beampipe invisible
/musr/command visattributes log_beampipeshell invisible
/musr/command visattributes log_beampipewindow invisible
/musr/command visattributes log_Magnet1 invisible
/musr/command visattributes log_Magnet2 invisible
/musr/command visattributes log_Uprof1 invisible
/musr/command visattributes log_Uprof2 invisible
/musr/command visattributes log_Uprof3 invisible
/musr/command visattributes log_Uprof4 invisible
/musr/command visattributes log_GPDsampleHolderD red
/musr/command visattributes log_GPDsampleHolderE red
/musr/command visattributes log_GPDsampleHolderF red
/musr/command visattributes log_GPDsampleHolderG red
/musr/command visattributes log_GPDcollimator green
/musr/command visattributes log_GPDcollimator2 blue
/musr/command visattributes log_GPDmHolder magenta
/musr/command visattributes log_target blue
/musr/command visattributes log_M0 yellow
#/musr/command visattributes G4_PLASTIC_SC_VINYLTOLUENE red
#/musr/command visattributes G4_Pb blue
####################################################################################
######################### M A G N E T I C F I E L D #########################
###################################################################################
# Set magnetic field (set field intensity in T and sigma in mm)
# syntax for magneticfield: fromfile filename fieldValue
# uniform fieldValue
# gaussian fieldValue sigma
#
/musr/command globalfield centralSolenoidField 0. 0. 0. fromfile 2D GPD_coil_primitive.table log_target 0.03
/musr/command globalfield setparameter SetLargestAcceptableStep 200
/musr/command globalfield setparameter SetMinimumEpsilonStep 5e-5
/musr/command globalfield setparameter SetMaximumEpsilonStep 0.001
/musr/command globalfield setparameter SetDeltaOneStep 0.1
/musr/command globalfield setparameter SetDeltaIntersection 0.01
/musr/command globalfield printparameters
/musr/command globalfield printFieldValueAtPoint 0 0 0
/musr/command globalfield printFieldValueAtPoint 0 0 10
/musr/command globalfield printFieldValueAtPoint 0 0 20
/musr/command globalfield printFieldValueAtPoint 0 0 30
/musr/command globalfield printFieldValueAtPoint 0 0 50
/musr/command globalfield printFieldValueAtPoint 0 0 100
/musr/command globalfield printFieldValueAtPoint 0 0 200
/musr/command globalfield printFieldValueAtPoint 0 0 300
/musr/command globalfield printFieldValueAtPoint 0 0 400
/musr/command globalfield printFieldValueAtPoint 0 0 500
/musr/command globalfield printFieldValueAtPoint 10 0 0
/musr/command globalfield printFieldValueAtPoint 100 0 0
/musr/command globalfield printFieldValueAtPoint 300 0 0
/musr/command globalfield printFieldValueAtPoint 0 10 0
/musr/command globalfield printFieldValueAtPoint 0 100 0
/musr/command globalfield printFieldValueAtPoint 0 300 0
###################################################################################
######################### P H Y S I C S P R O C E S S E S ##################
###################################################################################
# Geant 4.9.3
/musr/command process addDiscreteProcess gamma G4PhotoElectricEffect
/musr/command process addDiscreteProcess gamma G4ComptonScattering
/musr/command process addDiscreteProcess gamma G4GammaConversion
/musr/command process addDiscreteProcess gamma G4RayleighScattering
/musr/command process addProcess e- G4eMultipleScattering -1 1 1
/musr/command process addProcess e- G4eIonisation -1 2 2
/musr/command process addProcess e- G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eMultipleScattering -1 1 1
/musr/command process addProcess e+ G4eIonisation -1 2 2
/musr/command process addProcess e+ G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eplusAnnihilation 0 -1 4
/musr/command process addProcess mu- G4MuMultipleScattering -1 1 1
/musr/command process addProcess mu- G4MuIonisation -1 2 2
/musr/command process addProcess mu- G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu- G4MuPairProduction -1 4 4
/musr/command process addProcess mu+ G4MuMultipleScattering -1 1 1
/musr/command process addProcess mu+ G4MuIonisation -1 2 2
/musr/command process addProcess mu+ G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu+ G4MuPairProduction -1 4 4
#
#/home/install/geant4.9.3/source/processes/electromagnetic/utils/src/G4EnergyLossMessenger.cc
######## /process/msc/StepLimit Minimal | UseDistanceToBoundary | UseSafety
/process/msc/StepLimit UseSafety
#/process/msc/LateralDisplacement
/process/msc/RangeFactor 0.04
/process/msc/GeomFactor 2.5
#/process/msc/FactorForAngleLimit
/process/msc/Skin 3.0
#/process/msc/ThetaLimit 0.2 rad
###################################################################################
################## S O M E O T H E R P A R A M E T E R S ##################
###################################################################################
# Set the overall range cut
#/run/setCut 1 mm
# Set the range cut on a particular volumes (in mm)
#/musr/command SetUserLimits log_target 0.01
#/musr/command SetUserLimits log_targetscint 0.01
#/musr/command SetUserLimits log_M0 0.01
# Store all events into the ROOT tree or just the interesting ones ? (true is default)
/musr/command storeOnlyEventsWithHits false
# Set the minimum time separation between two subsequent signals in the same detector (in ns)
/musr/command signalSeparationTime 10
#
#/musr/run/howOftenToPrintEvent 10
/musr/command maximumRunTimeAllowed 220000
/musr/run/randomOption 2
###################################################################################
######################### R O O T O U T P U T ##############################
###################################################################################
/musr/command rootOutput det_VvvKine off
/musr/command rootOutput det_VvvX off
/musr/command rootOutput det_VvvY off
/musr/command rootOutput det_VvvZ off
/musr/command rootOutput det_VvvVolID off
/musr/command rootOutput det_VvvProcID off
/musr/command rootOutput det_VvvTrackID off
/musr/command rootOutput det_VvvParticleID off
/musr/command rootOutput nOptPhot off
/musr/command rootOutput odet_ID off
/musr/command rootOutput odet_nPhot off
/musr/command rootOutput odet_timeFirst off
/musr/command rootOutput odet_timeA off
/musr/command rootOutput odet_timeB off
/musr/command rootOutput odet_timeC off
/musr/command rootOutput odet_timeD off
/musr/command rootOutput odet_timeE off
/musr/command rootOutput odet_timeLast off
###################################################################################
######################### V I S U A L I S A T I O N ##############################
###################################################################################
/vis/disable
#/control/execute visFromToni.mac
#/control/execute visDawn40123.mac
#/control/execute visVRML.mac
###################################################################################
######################### P A R T I C L E G U N #################################
###################################################################################
/gun/vertex 0 0 -1000 mm
# FWHM 10mm ==> sigma = 10/2.354 = 4.2481mm
/gun/vertexsigma 25 25 0 mm
#---/gun/vertexboundary: rMaxAllowed, zMinAllowed, zMaxAllowed
/gun/vertexboundary 100 -999999 999999 mm
/gun/momentum 100.0 MeV
# sigma = 3% ==> sigma 27*0.03 = 0.81
/gun/momentumsmearing 3.00 MeV
#---/gun/momentumboundary: pMinAllowed, pMaxAllowed, dummy
#/gun/momentumboundary 20 40 0 MeV
# TURTLE
#/gun/turtlefilename FOR061_2008_04_22.DAT
#/gun/turtlefilename FOR070_2008_10_17_XXII.DAT
#/gun/turtlefilename FOR077_pie3_HiField_d05_x30.dat
#/gun/turtlefilename FOR077_reggiani_Jan2010_NEW.dat
#/gun/turtleZ0position -900 mm
#/gun/turtlefilename FOR070_GPD_beam_profile_45mm.dat
#/gun/turtleZ0position -500 mm
#---/gun/tilt: xangle, yangle, dummy
#/gun/tilt 0 0.5 0 deg
#---/gun/tiltsigma: xangleSigma, yangleSigma, dummy (1 degree on 1 meter ~ 17mm)
/gun/tiltsigma 0.15 0.15 0 deg
/gun/pitch 0.08 deg
/gun/muonPolarizVector 0 0 -1
#---/gun/decaytimelimits: decayMin, decayMax, decayTime
/gun/decaytimelimits -1 -1 2197.03 ns
###################################################################################
######################## B E A M O N #########################################
###################################################################################
/run/beamOn 16000000
#/run/beamOn 1

175
run/GPD/40508g.v1190 Normal file
View File

@ -0,0 +1,175 @@
# TDC V1190 Set up file generated by Deltat : 16-Apr-2008 10:30:27
# bins are always in TDC channel units (195.3125 ps/channel)
INSTRUMENT=GPS
DESCRIPTION=No Veto -- 1port
TYPE=TDCV1190
RESOLUTION=100
MDELAY=0
PDELAY=2000
#REWINDTIMEBINS=10000000
MCOINCIDENCEW=50
PCOINCIDENCEW=50
VCOINCIDENCEW=100
# 30000 mu/s
MUONRATEFACTOR=0.0937185
DATAWINDOWMIN=-5
DATAWINDOWMAX=10.0
PILEUPWINDOWMIN=-10.2
PILEUPWINDOWMAX=10.2
102; "M up"; M; 0.4; 2005; 103 -104;
1; "B1"; P; 0.4; 2005; -104; B1; 1485; 1515; 50995;
2; "B2"; P; 0.4; 2005; -104; B2; 1485; 1515; 50995;
4; "B4"; P; 0.4; 2005; -104; B4; 1485; 1515; 50995;
5; "B5"; P; 0.4; 2005; -104; B5; 1485; 1515; 50995;
6; "B6"; P; 0.4; 2005; -104; B6; 1485; 1515; 50995;
7; "B6"; P; 0.4; 2005; -104; B7; 1485; 1515; 50995;
11; "F1"; P; 0.4; 2005; -104; F11; 1485; 1515; 50995;
12; "F2"; P; 0.4; 2005; -104; F12; 1485; 1515; 50995;
13; "F3"; P; 0.4; 2005; -104; F13; 1485; 1515; 50995;
14; "F4"; P; 0.4; 2005; -104; F14; 1485; 1515; 50995;
15; "F5"; P; 0.4; 2005; -104; F15; 1485; 1515; 50995;
103; "M embed" K; 0.1; 2005;
104; "M2 Veto" K; 0.1; 2005;
! 21; "Coinc B1" K; 0.4; 2005;
! 22; "Coinc B2" K; 0.4; 2005;
!401; "Act Veto" V; 0.1; 2005;
! 3; "F right"; K; 0.4; 1263;
! 4; "B left"; P; 0.2; 1311; 5 -10 -11 -12 -13; Back; 1485; 1515; 50995;
! 5; "B right"; K; 0.2; 1262;
! 6; Up; P; 0.2; 1306; -10 -11 -12 -13; Up; 1490; 1520; 51000;
! 7; Down; P; 0.2; 1264; -10 -11 -12 -13; Down; 1525; 1555; 51035;
! 8; "R int"; P; 0.2; 1309; -10 -11 -12 -13; R_int; 1480; 1510; 50990;
! 9; "R ext"; P; 0.3; 1267; -10 -11 -12 -13; R_ext; 1520; 1550; 51020; 8; Rite
!10; "Bveto up"; V; 0.3; 2085;
!11; "Bveto down"; V; 0.3; 2089;
!12; "Bveto left"; V; 0.3; 2080;
!13; "Bveto right"; V; 0.3; 2079;
!14; "F center"; P; 0.2; 2073; -10 -11 -12 -13; F_cntr; 720; 750; 50230; 2; Forw
!15; Clock; C;
!16; Rejected; R;
!17; Ch17; N;
!18; Ch18; N;
!19; Ch19; N;
!20; Ch20; N;
!21; Ch21; N;
!22; Ch22; N;
!23; Ch23; N;
!24; Ch24; N;
!25; Ch25; N;
!26; Ch26; N;
!27; Ch27; N;
!28; Ch28; N;
!29; Ch29; N;
!30; Ch30; N;
!31; Ch31; N;
#DEBUGEVENT 704 10
#DEBUGEVENT 1427 10
#DEBUGEVENT 3222 10
#DEBUGEVENT 5223 10
#DEBUGEVENT 8351 10
#DEBUGEVENT 8581 10
#DEBUGEVENT 9041 10
#DEBUGEVENT 9645 10
counterPhaseShifts 1 0 2 0 11 180 12 180 13 180
counterGrouping B 1 2
counterGrouping F 11 12 13
counterGrouping L 4 5
counterGrouping R 14 15
musrTH2D hMuDecayMap "Muon decay map;z [mm]; r[mm]" 100 -80. 20. 40 0. 20. muDecayPosZ muDecayPosR
musrTH2D hMuDecayMappileup "Pileup muon decay map;z [mm]; r[mm]" 400 -300. 100. 200 0. 200. pileup_muDecayPosZ pileup_muDecayPosR
musrTH1D hMuDecayPosZ "Penetration of muons into the sample;z[mm];N" 100 -5.0 5. muDecayPosZ
musrTH1D hMuDecayPosZpileup "Decay of the muons contributing to pileup bgr.;z[mm];N" 400 -300 100 pileup_muDecayPosZ
musrTH1D hMuDecayPosX "X of decayed muons;x[mm];N" 100 -5.0 5. muDecayPosX
musrTH1D hdet_m0edep "Energy deposited in the trigger;E[MeV];N" 100 0. 2 det_m0edep
musrTH1D hdet_posEdep "Energy deposited in the positron Counters;E[MeV];N" 100 0. 2 det_posEdep
#musrTH1D hdet_time0 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time0
#musrTH1D hdet_time1 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time1
#musrTH1D hdet_time10 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time10
musrTH1D hgen_time10 "generated muSR spectra;time[#mus];N" 220 -11. 11. gen_time10
musrTH1D hdet_time10_MINUS_gen_time10 "detected - generated ;time_{det}-time_{gen} [ps];N" 100 -500. 500. det_time10_MINUS_gen_time10
#musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 4000 -2000. 2000. muDecayDetID
musrTH1D hdet_time10 "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time10
musrTH1D hdet_time10_Bgr "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time10
musrTH1D hdet_time20 "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time20
musrTH1D hdet_time20_Bgr "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time20
musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 2002 -1001. 1001. muDecayDetID
musrTH1D hMuDecayDetIDpileup "Detector ID where the pileup muons decay;Detector ID; N" 2002 -1001. 1001. pileup_muDecayDetID
humanDecayHistograms hMuDecayDetID hMuDecayDetIDpileup 201 Sample 211 SampleCell 300 Holder 301 Holder 302 Holder 303 Holder 304 Holder 305 Holder 306 Holder 307 Holder 308 Holder 309 Holder 310 Holder 313 Holder 314 Holder 315 Holder 316 Holder 251 Col_Cu 255 Col_1 257 Col_1 261 Col_2 263 Col_2 265 Al_Plate 102 M0 103 M0Holder 1 BC 2 BC 11 FC 12 FC 13 FC 21 BBC 22 BBC 401 Veto -1 World 291 Magnet 292 Magnet 293 Al_Uprof 230 Elsewhere 231 Elsewhere 232 Elsewhere 233 Elsewhere 256 Elsewhere 262 Elsewhere 267 Elsewhere -1000 Escaped
condition 0 alwaysTrue
condition 1 oncePerEvent
condition 2 muonDecayedInSample_gen
condition 3 muonTriggered_gen
condition 4 muonTriggered_det
condition 5 positronHit_det
condition 6 goodEvent_det
condition 7 goodEvent_gen
condition 8 goodEvent_det_AND_goodEvent_gen
condition 9 pileupEvent
condition 10 goodEvent_F_det
condition 11 goodEvent_B_det
condition 12 goodEvent_F_det_AND_pileupEvent
condition 13 goodEvent_B_det_AND_pileupEvent
condition 14 goodEvent_L_det
condition 15 goodEvent_R_det
fit hdet_time20_Bgr pol0 111 -1.9 -0.1 0
#fit hdet_time20 funct1 111 0.2 9.8 0 0 0.3 100
#fit hdet_time20 funct2 111 0.2 9.8 0 0 0.3 100 0
fit hdet_time20 funct4 111 0.2 9.8 0 0 0.3 100 0
#draw hMuDecayPosZ 0
draw hMuDecayPosZ 1
#draw hMuDecayPosZ 2
#draw hMuDecayPosZ 3
#draw hMuDecayPosZ 4
###draw hMuDecayMap 0
###draw hMuDecayMap 1
draw hdet_m0edep 0
draw hdet_m0edep 1
draw hdet_m0edep 2
draw hdet_m0edep 3
draw hdet_m0edep 4
draw hdet_m0edep 5
draw hdet_m0edep 6
draw hdet_time10 6
draw hdet_time20 10
draw hdet_time20 11
draw hdet_time20 14
draw hdet_time20 15
draw hdet_time20_Bgr 10
draw hdet_time20_Bgr 11
draw hdet_time20_Bgr 14
draw hdet_time20_Bgr 15
draw hgen_time10 7
draw hdet_time10_MINUS_gen_time10 8
draw hdet_posEdep 5
draw hMuDecayDetID 1
draw hdet_m0edep 9
draw hdet_time10 9
draw hMuDecayDetIDpileup 1
draw hMuDecayDetIDpileup 6
draw hMuDecayDetIDpileup 9
draw humanDecayHistograms 1
draw humanDecayHistograms 6
draw humanDecayHistograms 9
draw humanDecayPileupHistograms 1
draw humanDecayPileupHistograms 6
draw humanDecayPileupHistograms 9
draw hdet_time20 4
draw hdet_time20 6
$

178
run/GPD/40508h.v1190 Normal file
View File

@ -0,0 +1,178 @@
# TDC V1190 Set up file generated by Deltat : 16-Apr-2008 10:30:27
# bins are always in TDC channel units (195.3125 ps/channel)
INSTRUMENT=GPS
DESCRIPTION=No Veto -- 1port
TYPE=TDCV1190
RESOLUTION=100
MDELAY=0
PDELAY=2000
#REWINDTIMEBINS=10000000
MCOINCIDENCEW=50
PCOINCIDENCEW=50
VCOINCIDENCEW=100
# 30000 mu/s
MUONRATEFACTOR=0.0937185
DATAWINDOWMIN=-5
DATAWINDOWMAX=10.0
PILEUPWINDOWMIN=-10.2
PILEUPWINDOWMAX=10.2
102; "M up"; M; 0.4; 2005; 103 -104;
1; "B1"; P; 0.4; 2005; ; B1; 1485; 1515; 50995;
2; "B2"; P; 0.4; 2005; ; B2; 1485; 1515; 50995;
4; "B4"; P; 0.4; 2005; ; B4; 1485; 1515; 50995;
5; "B5"; P; 0.4; 2005; ; B5; 1485; 1515; 50995;
6; "B6"; P; 0.4; 2005; ; B6; 1485; 1515; 50995;
7; "B6"; P; 0.4; 2005; ; B7; 1485; 1515; 50995;
11; "F1"; P; 0.4; 2005; ; F11; 1485; 1515; 50995;
12; "F2"; P; 0.4; 2005; ; F12; 1485; 1515; 50995;
13; "F3"; P; 0.4; 2005; ; F13; 1485; 1515; 50995;
14; "F4"; P; 0.4; 2005; ; F14; 1485; 1515; 50995;
15; "F5"; P; 0.4; 2005; ; F15; 1485; 1515; 50995;
103; "M embed" K; 0.1; 2005;
104; "M2 Veto" V; 0.1; 2005;
! 21; "Coinc B1" K; 0.4; 2005;
! 22; "Coinc B2" K; 0.4; 2005;
!401; "Act Veto" V; 0.1; 2005;
! 3; "F right"; K; 0.4; 1263;
! 4; "B left"; P; 0.2; 1311; 5 -10 -11 -12 -13; Back; 1485; 1515; 50995;
! 5; "B right"; K; 0.2; 1262;
! 6; Up; P; 0.2; 1306; -10 -11 -12 -13; Up; 1490; 1520; 51000;
! 7; Down; P; 0.2; 1264; -10 -11 -12 -13; Down; 1525; 1555; 51035;
! 8; "R int"; P; 0.2; 1309; -10 -11 -12 -13; R_int; 1480; 1510; 50990;
! 9; "R ext"; P; 0.3; 1267; -10 -11 -12 -13; R_ext; 1520; 1550; 51020; 8; Rite
!10; "Bveto up"; V; 0.3; 2085;
!11; "Bveto down"; V; 0.3; 2089;
!12; "Bveto left"; V; 0.3; 2080;
!13; "Bveto right"; V; 0.3; 2079;
!14; "F center"; P; 0.2; 2073; -10 -11 -12 -13; F_cntr; 720; 750; 50230; 2; Forw
!15; Clock; C;
!16; Rejected; R;
!17; Ch17; N;
!18; Ch18; N;
!19; Ch19; N;
!20; Ch20; N;
!21; Ch21; N;
!22; Ch22; N;
!23; Ch23; N;
!24; Ch24; N;
!25; Ch25; N;
!26; Ch26; N;
!27; Ch27; N;
!28; Ch28; N;
!29; Ch29; N;
!30; Ch30; N;
!31; Ch31; N;
setSpecialAnticoincidenceTimeWindow 104 -10.2 10.2 mus
#DEBUGEVENT 704 10
#DEBUGEVENT 1427 10
#DEBUGEVENT 3222 10
#DEBUGEVENT 5223 10
#DEBUGEVENT 8351 10
#DEBUGEVENT 8581 10
#DEBUGEVENT 9041 10
#DEBUGEVENT 9645 10
counterPhaseShifts 1 0 2 0 11 180 12 180 13 180
counterGrouping B 1 2
counterGrouping F 11 12 13
counterGrouping L 4 5
counterGrouping R 14 15
musrTH2D hMuDecayMap "Muon decay map;z [mm]; r[mm]" 100 -80. 20. 40 0. 20. muDecayPosZ muDecayPosR
musrTH2D hMuDecayMappileup "Pileup muon decay map;z [mm]; r[mm]" 400 -300. 100. 200 0. 200. pileup_muDecayPosZ pileup_muDecayPosR
musrTH1D hMuDecayPosZ "Penetration of muons into the sample;z[mm];N" 100 -5.0 5. muDecayPosZ
musrTH1D hMuDecayPosZpileup "Decay of the muons contributing to pileup bgr.;z[mm];N" 400 -300 100 pileup_muDecayPosZ
musrTH1D hMuDecayPosX "X of decayed muons;x[mm];N" 100 -5.0 5. muDecayPosX
musrTH1D hdet_m0edep "Energy deposited in the trigger;E[MeV];N" 100 0. 2 det_m0edep
musrTH1D hdet_posEdep "Energy deposited in the positron Counters;E[MeV];N" 100 0. 2 det_posEdep
#musrTH1D hdet_time0 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time0
#musrTH1D hdet_time1 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time1
#musrTH1D hdet_time10 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time10
musrTH1D hgen_time10 "generated muSR spectra;time[#mus];N" 220 -11. 11. gen_time10
musrTH1D hdet_time10_MINUS_gen_time10 "detected - generated ;time_{det}-time_{gen} [ps];N" 100 -500. 500. det_time10_MINUS_gen_time10
#musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 4000 -2000. 2000. muDecayDetID
musrTH1D hdet_time10 "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time10
musrTH1D hdet_time10_Bgr "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time10
musrTH1D hdet_time20 "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time20
musrTH1D hdet_time20_Bgr "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time20
musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 2002 -1001. 1001. muDecayDetID
musrTH1D hMuDecayDetIDpileup "Detector ID where the pileup muons decay;Detector ID; N" 2002 -1001. 1001. pileup_muDecayDetID
humanDecayHistograms hMuDecayDetID hMuDecayDetIDpileup 201 Sample 211 SampleCell 300 Holder 301 Holder 302 Holder 303 Holder 304 Holder 305 Holder 306 Holder 307 Holder 308 Holder 309 Holder 310 Holder 313 Holder 314 Holder 315 Holder 316 Holder 251 Col_Cu 255 Col_1 257 Col_1 261 Col_2 263 Col_2 265 Al_Plate 102 M0 103 M0Holder 1 BC 2 BC 11 FC 12 FC 13 FC 21 BBC 22 BBC 401 Veto -1 World 291 Magnet 292 Magnet 293 Al_Uprof 230 Elsewhere 231 Elsewhere 232 Elsewhere 233 Elsewhere 256 Elsewhere 262 Elsewhere 267 Elsewhere -1000 Escaped
condition 0 alwaysTrue
condition 1 oncePerEvent
condition 2 muonDecayedInSample_gen
condition 3 muonTriggered_gen
condition 4 muonTriggered_det
condition 5 positronHit_det
condition 6 goodEvent_det
condition 7 goodEvent_gen
condition 8 goodEvent_det_AND_goodEvent_gen
condition 9 pileupEvent
condition 10 goodEvent_F_det
condition 11 goodEvent_B_det
condition 12 goodEvent_F_det_AND_pileupEvent
condition 13 goodEvent_B_det_AND_pileupEvent
condition 14 goodEvent_L_det
condition 15 goodEvent_R_det
fit hdet_time20_Bgr pol0 111 -1.9 -0.1 0
#fit hdet_time20 funct1 111 0.2 9.8 0 0 0.3 100
#fit hdet_time20 funct2 111 0.2 9.8 0 0 0.3 100 0
fit hdet_time20 funct4 111 0.2 9.8 0 0 0.3 100 0
#draw hMuDecayPosZ 0
draw hMuDecayPosZ 1
#draw hMuDecayPosZ 2
#draw hMuDecayPosZ 3
#draw hMuDecayPosZ 4
###draw hMuDecayMap 0
###draw hMuDecayMap 1
draw hdet_m0edep 0
draw hdet_m0edep 1
draw hdet_m0edep 2
draw hdet_m0edep 3
draw hdet_m0edep 4
draw hdet_m0edep 5
draw hdet_m0edep 6
draw hdet_time10 6
draw hdet_time20 10
draw hdet_time20 11
draw hdet_time20 14
draw hdet_time20 15
draw hdet_time20_Bgr 10
draw hdet_time20_Bgr 11
draw hdet_time20_Bgr 14
draw hdet_time20_Bgr 15
draw hgen_time10 7
draw hdet_time10_MINUS_gen_time10 8
draw hdet_posEdep 5
draw hMuDecayDetID 1
draw hdet_m0edep 9
draw hdet_time10 9
draw hMuDecayDetIDpileup 1
draw hMuDecayDetIDpileup 6
draw hMuDecayDetIDpileup 9
draw humanDecayHistograms 1
draw humanDecayHistograms 6
draw humanDecayHistograms 9
draw humanDecayPileupHistograms 1
draw humanDecayPileupHistograms 6
draw humanDecayPileupHistograms 9
draw hdet_time20 4
draw hdet_time20 6
$

263
run/GPD/40601.mac Normal file
View File

@ -0,0 +1,263 @@
#-----------------------------------------------------------------------
# Macro file for muSR instrument simulations
# Unless specified otherwises, the default units are mm, ns, MeV, MeV/c.
# Lines starting with star "#" are comments.
###################################################################################
############################# G E O M E T R Y ###################################
###################################################################################
# ROTATION MATRIXES:
/musr/command rotation matrix1 1 0 0 90
/musr/command rotation matrix2 0 1 0 28.81
/musr/command rotation matrix3 0 1 0 151.19
/musr/command rotation matrix4 0 0 1 45
/musr/command rotation matrix5 0 1 0 90
/musr/command rotation matrix6 180 90 90
/musr/command rotation matrix7 0 90 90
# --- WORLD VOLUME ---
/musr/command construct box World 2100 2100 11000 G4_AIR 0 0 0 no_logical_volume norot dead -1
# --- BEAMPIPE ---
/musr/command construct tubs beampipeA 0 2005 4500 0 360 G4_Galactic 0 0 -6100 log_World norot dead 230
/musr/command construct tubs beampipe 0 405 500 0 360 G4_Galactic 0 0 -1100 log_World norot dead 231
/musr/command construct tubs beampipeshell 100 105 500 0 360 G4_Al 0 0 0 log_beampipe norot dead 232
/musr/command construct tubs beampipewindow 0 100 0.05 0 360 G4_MYLAR 0 0 499.9 log_beampipe norot dead 233
# PASSIVE COLLIMATORS (for BEAMPIPE vers.1)
# GPDcollimator - thickness in the thinner location is x2+(x9-x4) = 20+59-59
##/musr/command construct GPDcollimator GPDcollimator 40 15 60 59 60.01 0 360 0 59 0 2.5 6 G4_Cu 0 0 -59 log_World matrix1 dead 251
/musr/command construct box GPDcollimator2 48 55 15 G4_Pb 0 0 -66 log_World norot dead 255
/musr/command construct box GPDcollimator2hole 2.5 4 15 G4_Galactic 0 0 0 log_GPDcollimator2 norot dead 256
/musr/command construct box FirstCollimator 45 75 15 G4_Pb 0 0 -250 log_World norot dead 261
/musr/command construct tubs FirstCollHole 0 8 15 0 360 G4_Galactic 0 0 0 log_FirstCollimator norot dead 262
/musr/command construct box SplitCollimator1 19 75 25 G4_Cu 23 0 -158 log_World norot dead 261
/musr/command construct box SplitCollimator2 19 75 25 G4_Cu -23 0 -158 log_World norot dead 261
###/musr/command construct box AlumPlate 30 90 5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -80 log_World norot musr/ScintSD 401
###/musr/command construct box AlumPlateHole 2 5 5 G4_Galactic 0 0 0 log_AlumPlate norot dead 267
/musr/command construct box Insert 20 20 10 G4_Pb 0 0 -91 log_World norot dead 265
/musr/command construct box InsertHole 2.5 4 10 G4_Galactic 0 0 0 log_Insert norot dead 267
# --- SAMPLE HOLDER ---
##/musr/command construct GPDsampleHolderA GPDsampleHolderA 0 0 0 0 0 G4_Al 0 -111.25 0 log_World matrix1 dead 300
##/musr/command construct tubs GPDsampleHolderB 37.51 58 7 0 360 G4_Al 0 57 0 log_World matrix1 dead 301
##/musr/command construct tubs GPDsampleHolderC 37.51 58 7 0 360 G4_Al 0 -57 0 log_World matrix1 dead 302
######/musr/command construct tubs GPDsampleHolderD 37 37.5 50 136 88 G4_Al 0 0 0 log_World matrix1 dead 303
######/musr/command construct tubs GPDsampleHolderE 53.5 54 50 136 88 G4_Al 0 0 0 log_World matrix1 dead 304
##/musr/command construct box GPDsampleHolderD 12 0.5 50 G4_Al 38 0 28 log_World matrix1 dead 303
##/musr/command construct box GPDsampleHolderE 12 0.5 50 G4_Al 38 0 -28 log_World matrix1 dead 304
##/musr/command construct box GPDsampleHolderF 12 0.5 50 G4_Al -38 0 28 log_World matrix1 dead 305
##/musr/command construct box GPDsampleHolderG 12 0.5 50 G4_Al -38 0 -28 log_World matrix1 dead 306
##/musr/command construct tubs GPDsampleHolderH 0 9 50 180 50 G4_Al 48 0 29 log_World matrix1 dead 307
##/musr/command construct tubs GPDsampleHolderI 0 9 50 130 50 G4_Al 48 0 -29 log_World matrix1 dead 308
##/musr/command construct tubs GPDsampleHolderJ 0 9 50 310 50 G4_Al -48 0 29 log_World matrix1 dead 309
##/musr/command construct tubs GPDsampleHolderK 0 9 50 0 50 G4_Al -48 0 -29 log_World matrix1 dead 310
##/musr/command construct tubs GPDsampleHolderW 58.01 58.2 64 25 130 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 313
##/musr/command construct tubs GPDsampleHolderX 58.01 58.2 64 205 130 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 314
##/musr/command construct tubs GPDsampleHolderY 36.79 36.99 64 45 90 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 315
##/musr/command construct tubs GPDsampleHolderZ 36.79 36.99 64 225 90 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 316
# --- SAMPLE CELL ---
/musr/command construct tubs GPDsampleCell 0 12 50 0 360 G4_Cu 0 0 0 log_World matrix1 dead 211
# --- SAMPLE ---
/musr/command construct tubs target 0 3 5 0 360 G4_Cu 0 0 0 log_GPDsampleCell norot dead 201
#/musr/command construct tubs targetFieldVol 0 3 5 0 360 G4_Cu 0 0 0 log_GPDsampleCell norot dead 201
# --- SCINTILLATOR ---
/musr/command construct box ScintB1 21.5 50 2 G4_PLASTIC_SC_VINYLTOLUENE 26.5 0 -44 log_World norot musr/ScintSD 1
/musr/command construct box ScintB2 21.5 50 2 G4_PLASTIC_SC_VINYLTOLUENE -26.5 0 -44 log_World norot musr/ScintSD 2
/musr/command construct box ScintB4 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE 44 0 -25 log_World norot musr/ScintSD 4
/musr/command construct box ScintB5 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE -44 0 -25 log_World norot musr/ScintSD 5
/musr/command construct box ScintB6 2 50 13 G4_PLASTIC_SC_VINYLTOLUENE 44 0 0 log_World norot musr/ScintSD 6
/musr/command construct box ScintB7 2 50 13 G4_PLASTIC_SC_VINYLTOLUENE -44 0 0 log_World norot musr/ScintSD 7
/musr/command construct box ScintF1 20 50 2 G4_PLASTIC_SC_VINYLTOLUENE 28 0 44 log_World norot musr/ScintSD 11
/musr/command construct box ScintF2 20 50 2 G4_PLASTIC_SC_VINYLTOLUENE -28 0 44 log_World norot musr/ScintSD 12
/musr/command construct box ScintF3 7 50 2 G4_PLASTIC_SC_VINYLTOLUENE 0 0 44 log_World norot musr/ScintSD 13
/musr/command construct box ScintF4 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE 44 0 25 log_World norot musr/ScintSD 14
/musr/command construct box ScintF5 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE -44 0 25 log_World norot musr/ScintSD 15
#/musr/command construct box KoincB1 20 90 2 G4_PLASTIC_SC_VINYLTOLUENE 33 0 -60 log_World matrix2 musr/ScintSD 21
#/musr/command construct box KoincB2 20 90 2 G4_PLASTIC_SC_VINYLTOLUENE -33 0 -60 log_World matrix3 musr/ScintSD 22
# --- M0 ---
/musr/command construct box M0holder 3 4 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -20 log_World norot musr/ScintSD 103
/musr/command construct box M0 4 30 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -44 log_World norot musr/ScintSD 102
#/musr/command construct GPDmHolder GPDmHolder 5 25 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 -25 -43 log_World norot dead 103
# --- Magnet ---
/musr/command construct tubs Magnet1 180 375 78 0 360 G4_Cu 133 0 0 log_World matrix5 dead 291
/musr/command construct tubs Magnet2 180 375 78 0 360 G4_Cu -133 0 0 log_World matrix5 dead 292
/musr/command construct uprofile Uprof1 25 10 400 5 G4_Al -45 0 200 log_World matrix6 dead 293
/musr/command construct uprofile Uprof2 25 10 400 5 G4_Al -45 0 -200 log_World matrix6 dead 293
/musr/command construct uprofile Uprof3 25 10 400 5 G4_Al 45 0 200 log_World matrix7 dead 293
/musr/command construct uprofile Uprof4 25 10 400 5 G4_Al 45 0 -200 log_World matrix7 dead 293
#
##---------------------------------------------------------
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_1 6
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_2 5
#/musr/command logicalVolumeToBeReweighted mu log_activeCollim 4
##---------------------------------------------------------
/musr/command visattributes log_World invisible
/musr/command visattributes log_beampipeA invisible
/musr/command visattributes log_beampipe invisible
/musr/command visattributes log_beampipeshell invisible
/musr/command visattributes log_beampipewindow invisible
/musr/command visattributes log_Magnet1 invisible
/musr/command visattributes log_Magnet2 invisible
/musr/command visattributes log_Uprof1 invisible
/musr/command visattributes log_Uprof2 invisible
/musr/command visattributes log_Uprof3 invisible
/musr/command visattributes log_Uprof4 invisible
/musr/command visattributes log_GPDsampleHolderD red
/musr/command visattributes log_GPDsampleHolderE red
/musr/command visattributes log_GPDsampleHolderF red
/musr/command visattributes log_GPDsampleHolderG red
/musr/command visattributes log_GPDcollimator green
/musr/command visattributes log_GPDcollimator2 blue
/musr/command visattributes log_GPDmHolder magenta
/musr/command visattributes log_target blue
/musr/command visattributes log_M0 yellow
#/musr/command visattributes G4_PLASTIC_SC_VINYLTOLUENE red
#/musr/command visattributes G4_Pb blue
####################################################################################
######################### M A G N E T I C F I E L D #########################
###################################################################################
# Set magnetic field (set field intensity in T and sigma in mm)
# syntax for magneticfield: fromfile filename fieldValue
# uniform fieldValue
# gaussian fieldValue sigma
#
/musr/command globalfield centralSolenoidField 0. 0. 0. fromfile 2D GPD_coil_primitive.table log_target 0.03
/musr/command globalfield setparameter SetLargestAcceptableStep 200
/musr/command globalfield setparameter SetMinimumEpsilonStep 5e-5
/musr/command globalfield setparameter SetMaximumEpsilonStep 0.001
/musr/command globalfield setparameter SetDeltaOneStep 0.1
/musr/command globalfield setparameter SetDeltaIntersection 0.01
/musr/command globalfield printparameters
/musr/command globalfield printFieldValueAtPoint 0 0 0
/musr/command globalfield printFieldValueAtPoint 0 0 10
/musr/command globalfield printFieldValueAtPoint 0 0 20
/musr/command globalfield printFieldValueAtPoint 0 0 30
/musr/command globalfield printFieldValueAtPoint 0 0 50
/musr/command globalfield printFieldValueAtPoint 0 0 100
/musr/command globalfield printFieldValueAtPoint 0 0 200
/musr/command globalfield printFieldValueAtPoint 0 0 300
/musr/command globalfield printFieldValueAtPoint 0 0 400
/musr/command globalfield printFieldValueAtPoint 0 0 500
/musr/command globalfield printFieldValueAtPoint 10 0 0
/musr/command globalfield printFieldValueAtPoint 100 0 0
/musr/command globalfield printFieldValueAtPoint 300 0 0
/musr/command globalfield printFieldValueAtPoint 0 10 0
/musr/command globalfield printFieldValueAtPoint 0 100 0
/musr/command globalfield printFieldValueAtPoint 0 300 0
###################################################################################
######################### P H Y S I C S P R O C E S S E S ##################
###################################################################################
# Geant 4.9.3
/musr/command process addDiscreteProcess gamma G4PhotoElectricEffect
/musr/command process addDiscreteProcess gamma G4ComptonScattering
/musr/command process addDiscreteProcess gamma G4GammaConversion
/musr/command process addDiscreteProcess gamma G4RayleighScattering
/musr/command process addProcess e- G4eMultipleScattering -1 1 1
/musr/command process addProcess e- G4eIonisation -1 2 2
/musr/command process addProcess e- G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eMultipleScattering -1 1 1
/musr/command process addProcess e+ G4eIonisation -1 2 2
/musr/command process addProcess e+ G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eplusAnnihilation 0 -1 4
/musr/command process addProcess mu- G4MuMultipleScattering -1 1 1
/musr/command process addProcess mu- G4MuIonisation -1 2 2
/musr/command process addProcess mu- G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu- G4MuPairProduction -1 4 4
/musr/command process addProcess mu+ G4MuMultipleScattering -1 1 1
/musr/command process addProcess mu+ G4MuIonisation -1 2 2
/musr/command process addProcess mu+ G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu+ G4MuPairProduction -1 4 4
#
#/home/install/geant4.9.3/source/processes/electromagnetic/utils/src/G4EnergyLossMessenger.cc
######## /process/msc/StepLimit Minimal | UseDistanceToBoundary | UseSafety
/process/msc/StepLimit UseSafety
#/process/msc/LateralDisplacement
/process/msc/RangeFactor 0.04
/process/msc/GeomFactor 2.5
#/process/msc/FactorForAngleLimit
/process/msc/Skin 3.0
#/process/msc/ThetaLimit 0.2 rad
###################################################################################
################## S O M E O T H E R P A R A M E T E R S ##################
###################################################################################
# Set the overall range cut
#/run/setCut 1 mm
# Set the range cut on a particular volumes (in mm)
#/musr/command SetUserLimits log_target 0.01
#/musr/command SetUserLimits log_targetscint 0.01
#/musr/command SetUserLimits log_M0 0.01
# Store all events into the ROOT tree or just the interesting ones ? (true is default)
/musr/command storeOnlyEventsWithHits false
# Set the minimum time separation between two subsequent signals in the same detector (in ns)
/musr/command signalSeparationTime 10
#
#/musr/run/howOftenToPrintEvent 10
/musr/command maximumRunTimeAllowed 220000
/musr/run/randomOption 2
###################################################################################
######################### R O O T O U T P U T ##############################
###################################################################################
/musr/command rootOutput det_VvvKine off
/musr/command rootOutput det_VvvX off
/musr/command rootOutput det_VvvY off
/musr/command rootOutput det_VvvZ off
/musr/command rootOutput det_VvvVolID off
/musr/command rootOutput det_VvvProcID off
/musr/command rootOutput det_VvvTrackID off
/musr/command rootOutput det_VvvParticleID off
/musr/command rootOutput nOptPhot off
/musr/command rootOutput odet_ID off
/musr/command rootOutput odet_nPhot off
/musr/command rootOutput odet_timeFirst off
/musr/command rootOutput odet_timeA off
/musr/command rootOutput odet_timeB off
/musr/command rootOutput odet_timeC off
/musr/command rootOutput odet_timeD off
/musr/command rootOutput odet_timeE off
/musr/command rootOutput odet_timeLast off
###################################################################################
######################### V I S U A L I S A T I O N ##############################
###################################################################################
#/vis/disable
#/control/execute visFromToni.mac
#/control/execute visDawn101.mac
#/control/execute visDawn40123.mac
/control/execute visVRML.mac
###################################################################################
######################### P A R T I C L E G U N #################################
###################################################################################
/gun/vertex 0 0 -1000 mm
# FWHM 10mm ==> sigma = 10/2.354 = 4.2481mm
/gun/vertexsigma 25 25 0 mm
#---/gun/vertexboundary: rMaxAllowed, zMinAllowed, zMaxAllowed
/gun/vertexboundary 100 -999999 999999 mm
/gun/momentum 100.0 MeV
# sigma = 3% ==> sigma 27*0.03 = 0.81
/gun/momentumsmearing 3.00 MeV
#---/gun/momentumboundary: pMinAllowed, pMaxAllowed, dummy
#/gun/momentumboundary 20 40 0 MeV
# TURTLE
#/gun/turtlefilename FOR061_2008_04_22.DAT
#/gun/turtlefilename FOR070_2008_10_17_XXII.DAT
#/gun/turtlefilename FOR077_pie3_HiField_d05_x30.dat
#/gun/turtlefilename FOR077_reggiani_Jan2010_NEW.dat
#/gun/turtleZ0position -900 mm
#/gun/turtlefilename FOR070_GPD_beam_profile_45mm.dat
#/gun/turtleZ0position -500 mm
#---/gun/tilt: xangle, yangle, dummy
#/gun/tilt 0 0.5 0 deg
#---/gun/tiltsigma: xangleSigma, yangleSigma, dummy (1 degree on 1 meter ~ 17mm)
/gun/tiltsigma 0.15 0.15 0 deg
/gun/pitch 0.08 deg
/gun/muonPolarizVector 0 0 -1
#---/gun/decaytimelimits: decayMin, decayMax, decayTime
/gun/decaytimelimits -1 -1 2197.03 ns
###################################################################################
######################## B E A M O N #########################################
###################################################################################
#/run/beamOn 16000000
/run/beamOn 1

263
run/GPD/40602.mac Normal file
View File

@ -0,0 +1,263 @@
#-----------------------------------------------------------------------
# Macro file for muSR instrument simulations
# Unless specified otherwises, the default units are mm, ns, MeV, MeV/c.
# Lines starting with star "#" are comments.
###################################################################################
############################# G E O M E T R Y ###################################
###################################################################################
# ROTATION MATRIXES:
/musr/command rotation matrix1 1 0 0 90
/musr/command rotation matrix2 0 1 0 28.81
/musr/command rotation matrix3 0 1 0 151.19
/musr/command rotation matrix4 0 0 1 45
/musr/command rotation matrix5 0 1 0 90
/musr/command rotation matrix6 180 90 90
/musr/command rotation matrix7 0 90 90
# --- WORLD VOLUME ---
/musr/command construct box World 2100 2100 11000 G4_AIR 0 0 0 no_logical_volume norot dead -1
# --- BEAMPIPE ---
/musr/command construct tubs beampipeA 0 2005 4500 0 360 G4_Galactic 0 0 -6100 log_World norot dead 230
/musr/command construct tubs beampipe 0 405 500 0 360 G4_Galactic 0 0 -1100 log_World norot dead 231
/musr/command construct tubs beampipeshell 100 105 500 0 360 G4_Al 0 0 0 log_beampipe norot dead 232
/musr/command construct tubs beampipewindow 0 100 0.05 0 360 G4_MYLAR 0 0 499.9 log_beampipe norot dead 233
# PASSIVE COLLIMATORS (for BEAMPIPE vers.1)
# GPDcollimator - thickness in the thinner location is x2+(x9-x4) = 20+59-59
##/musr/command construct GPDcollimator GPDcollimator 40 15 60 59 60.01 0 360 0 59 0 2.5 6 G4_Cu 0 0 -59 log_World matrix1 dead 251
/musr/command construct box GPDcollimator2 48 55 15 G4_Pb 0 0 -66 log_World norot dead 255
/musr/command construct box GPDcollimator2hole 2.5 4 15 G4_Galactic 0 0 0 log_GPDcollimator2 norot dead 256
/musr/command construct box FirstCollimator 45 75 15 G4_Pb 0 0 -250 log_World norot dead 261
/musr/command construct tubs FirstCollHole 0 8 15 0 360 G4_Galactic 0 0 0 log_FirstCollimator norot dead 262
/musr/command construct box SplitCollimator1 19 75 25 G4_Cu 22 0 -158 log_World norot dead 261
/musr/command construct box SplitCollimator2 19 75 25 G4_Cu -22 0 -158 log_World norot dead 261
###/musr/command construct box AlumPlate 30 90 5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -80 log_World norot musr/ScintSD 401
###/musr/command construct box AlumPlateHole 2 5 5 G4_Galactic 0 0 0 log_AlumPlate norot dead 267
/musr/command construct box Insert 20 20 10 G4_Pb 0 0 -91 log_World norot dead 265
/musr/command construct box InsertHole 2.5 4 10 G4_Galactic 0 0 0 log_Insert norot dead 267
# --- SAMPLE HOLDER ---
##/musr/command construct GPDsampleHolderA GPDsampleHolderA 0 0 0 0 0 G4_Al 0 -111.25 0 log_World matrix1 dead 300
##/musr/command construct tubs GPDsampleHolderB 37.51 58 7 0 360 G4_Al 0 57 0 log_World matrix1 dead 301
##/musr/command construct tubs GPDsampleHolderC 37.51 58 7 0 360 G4_Al 0 -57 0 log_World matrix1 dead 302
######/musr/command construct tubs GPDsampleHolderD 37 37.5 50 136 88 G4_Al 0 0 0 log_World matrix1 dead 303
######/musr/command construct tubs GPDsampleHolderE 53.5 54 50 136 88 G4_Al 0 0 0 log_World matrix1 dead 304
##/musr/command construct box GPDsampleHolderD 12 0.5 50 G4_Al 38 0 28 log_World matrix1 dead 303
##/musr/command construct box GPDsampleHolderE 12 0.5 50 G4_Al 38 0 -28 log_World matrix1 dead 304
##/musr/command construct box GPDsampleHolderF 12 0.5 50 G4_Al -38 0 28 log_World matrix1 dead 305
##/musr/command construct box GPDsampleHolderG 12 0.5 50 G4_Al -38 0 -28 log_World matrix1 dead 306
##/musr/command construct tubs GPDsampleHolderH 0 9 50 180 50 G4_Al 48 0 29 log_World matrix1 dead 307
##/musr/command construct tubs GPDsampleHolderI 0 9 50 130 50 G4_Al 48 0 -29 log_World matrix1 dead 308
##/musr/command construct tubs GPDsampleHolderJ 0 9 50 310 50 G4_Al -48 0 29 log_World matrix1 dead 309
##/musr/command construct tubs GPDsampleHolderK 0 9 50 0 50 G4_Al -48 0 -29 log_World matrix1 dead 310
##/musr/command construct tubs GPDsampleHolderW 58.01 58.2 64 25 130 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 313
##/musr/command construct tubs GPDsampleHolderX 58.01 58.2 64 205 130 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 314
##/musr/command construct tubs GPDsampleHolderY 36.79 36.99 64 45 90 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 315
##/musr/command construct tubs GPDsampleHolderZ 36.79 36.99 64 225 90 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 316
# --- SAMPLE CELL ---
/musr/command construct tubs GPDsampleCell 0 12 50 0 360 G4_Cu 0 0 0 log_World matrix1 dead 211
# --- SAMPLE ---
/musr/command construct tubs target 0 3 5 0 360 G4_Cu 0 0 0 log_GPDsampleCell norot dead 201
#/musr/command construct tubs targetFieldVol 0 3 5 0 360 G4_Cu 0 0 0 log_GPDsampleCell norot dead 201
# --- SCINTILLATOR ---
/musr/command construct box ScintB1 21.5 50 2 G4_PLASTIC_SC_VINYLTOLUENE 26.5 0 -44 log_World norot musr/ScintSD 1
/musr/command construct box ScintB2 21.5 50 2 G4_PLASTIC_SC_VINYLTOLUENE -26.5 0 -44 log_World norot musr/ScintSD 2
/musr/command construct box ScintB4 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE 44 0 -25 log_World norot musr/ScintSD 4
/musr/command construct box ScintB5 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE -44 0 -25 log_World norot musr/ScintSD 5
/musr/command construct box ScintB6 2 50 13 G4_PLASTIC_SC_VINYLTOLUENE 44 0 0 log_World norot musr/ScintSD 6
/musr/command construct box ScintB7 2 50 13 G4_PLASTIC_SC_VINYLTOLUENE -44 0 0 log_World norot musr/ScintSD 7
/musr/command construct box ScintF1 20 50 2 G4_PLASTIC_SC_VINYLTOLUENE 28 0 44 log_World norot musr/ScintSD 11
/musr/command construct box ScintF2 20 50 2 G4_PLASTIC_SC_VINYLTOLUENE -28 0 44 log_World norot musr/ScintSD 12
/musr/command construct box ScintF3 7 50 2 G4_PLASTIC_SC_VINYLTOLUENE 0 0 44 log_World norot musr/ScintSD 13
/musr/command construct box ScintF4 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE 44 0 25 log_World norot musr/ScintSD 14
/musr/command construct box ScintF5 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE -44 0 25 log_World norot musr/ScintSD 15
#/musr/command construct box KoincB1 20 90 2 G4_PLASTIC_SC_VINYLTOLUENE 33 0 -60 log_World matrix2 musr/ScintSD 21
#/musr/command construct box KoincB2 20 90 2 G4_PLASTIC_SC_VINYLTOLUENE -33 0 -60 log_World matrix3 musr/ScintSD 22
# --- M0 ---
/musr/command construct box M0holder 3 4 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -20 log_World norot musr/ScintSD 103
/musr/command construct box M0 4 30 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -44 log_World norot musr/ScintSD 102
#/musr/command construct GPDmHolder GPDmHolder 5 25 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 -25 -43 log_World norot dead 103
# --- Magnet ---
/musr/command construct tubs Magnet1 180 375 78 0 360 G4_Cu 133 0 0 log_World matrix5 dead 291
/musr/command construct tubs Magnet2 180 375 78 0 360 G4_Cu -133 0 0 log_World matrix5 dead 292
/musr/command construct uprofile Uprof1 25 10 400 5 G4_Al -45 0 200 log_World matrix6 dead 293
/musr/command construct uprofile Uprof2 25 10 400 5 G4_Al -45 0 -200 log_World matrix6 dead 293
/musr/command construct uprofile Uprof3 25 10 400 5 G4_Al 45 0 200 log_World matrix7 dead 293
/musr/command construct uprofile Uprof4 25 10 400 5 G4_Al 45 0 -200 log_World matrix7 dead 293
#
##---------------------------------------------------------
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_1 6
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_2 5
#/musr/command logicalVolumeToBeReweighted mu log_activeCollim 4
##---------------------------------------------------------
/musr/command visattributes log_World invisible
/musr/command visattributes log_beampipeA invisible
/musr/command visattributes log_beampipe invisible
/musr/command visattributes log_beampipeshell invisible
/musr/command visattributes log_beampipewindow invisible
/musr/command visattributes log_Magnet1 invisible
/musr/command visattributes log_Magnet2 invisible
/musr/command visattributes log_Uprof1 invisible
/musr/command visattributes log_Uprof2 invisible
/musr/command visattributes log_Uprof3 invisible
/musr/command visattributes log_Uprof4 invisible
/musr/command visattributes log_GPDsampleHolderD red
/musr/command visattributes log_GPDsampleHolderE red
/musr/command visattributes log_GPDsampleHolderF red
/musr/command visattributes log_GPDsampleHolderG red
/musr/command visattributes log_GPDcollimator green
/musr/command visattributes log_GPDcollimator2 blue
/musr/command visattributes log_GPDmHolder magenta
/musr/command visattributes log_target blue
/musr/command visattributes log_M0 yellow
#/musr/command visattributes G4_PLASTIC_SC_VINYLTOLUENE red
#/musr/command visattributes G4_Pb blue
####################################################################################
######################### M A G N E T I C F I E L D #########################
###################################################################################
# Set magnetic field (set field intensity in T and sigma in mm)
# syntax for magneticfield: fromfile filename fieldValue
# uniform fieldValue
# gaussian fieldValue sigma
#
/musr/command globalfield centralSolenoidField 0. 0. 0. fromfile 2D GPD_coil_primitive.table log_target 0.03
/musr/command globalfield setparameter SetLargestAcceptableStep 200
/musr/command globalfield setparameter SetMinimumEpsilonStep 5e-5
/musr/command globalfield setparameter SetMaximumEpsilonStep 0.001
/musr/command globalfield setparameter SetDeltaOneStep 0.1
/musr/command globalfield setparameter SetDeltaIntersection 0.01
/musr/command globalfield printparameters
/musr/command globalfield printFieldValueAtPoint 0 0 0
/musr/command globalfield printFieldValueAtPoint 0 0 10
/musr/command globalfield printFieldValueAtPoint 0 0 20
/musr/command globalfield printFieldValueAtPoint 0 0 30
/musr/command globalfield printFieldValueAtPoint 0 0 50
/musr/command globalfield printFieldValueAtPoint 0 0 100
/musr/command globalfield printFieldValueAtPoint 0 0 200
/musr/command globalfield printFieldValueAtPoint 0 0 300
/musr/command globalfield printFieldValueAtPoint 0 0 400
/musr/command globalfield printFieldValueAtPoint 0 0 500
/musr/command globalfield printFieldValueAtPoint 10 0 0
/musr/command globalfield printFieldValueAtPoint 100 0 0
/musr/command globalfield printFieldValueAtPoint 300 0 0
/musr/command globalfield printFieldValueAtPoint 0 10 0
/musr/command globalfield printFieldValueAtPoint 0 100 0
/musr/command globalfield printFieldValueAtPoint 0 300 0
###################################################################################
######################### P H Y S I C S P R O C E S S E S ##################
###################################################################################
# Geant 4.9.3
/musr/command process addDiscreteProcess gamma G4PhotoElectricEffect
/musr/command process addDiscreteProcess gamma G4ComptonScattering
/musr/command process addDiscreteProcess gamma G4GammaConversion
/musr/command process addDiscreteProcess gamma G4RayleighScattering
/musr/command process addProcess e- G4eMultipleScattering -1 1 1
/musr/command process addProcess e- G4eIonisation -1 2 2
/musr/command process addProcess e- G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eMultipleScattering -1 1 1
/musr/command process addProcess e+ G4eIonisation -1 2 2
/musr/command process addProcess e+ G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eplusAnnihilation 0 -1 4
/musr/command process addProcess mu- G4MuMultipleScattering -1 1 1
/musr/command process addProcess mu- G4MuIonisation -1 2 2
/musr/command process addProcess mu- G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu- G4MuPairProduction -1 4 4
/musr/command process addProcess mu+ G4MuMultipleScattering -1 1 1
/musr/command process addProcess mu+ G4MuIonisation -1 2 2
/musr/command process addProcess mu+ G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu+ G4MuPairProduction -1 4 4
#
#/home/install/geant4.9.3/source/processes/electromagnetic/utils/src/G4EnergyLossMessenger.cc
######## /process/msc/StepLimit Minimal | UseDistanceToBoundary | UseSafety
/process/msc/StepLimit UseSafety
#/process/msc/LateralDisplacement
/process/msc/RangeFactor 0.04
/process/msc/GeomFactor 2.5
#/process/msc/FactorForAngleLimit
/process/msc/Skin 3.0
#/process/msc/ThetaLimit 0.2 rad
###################################################################################
################## S O M E O T H E R P A R A M E T E R S ##################
###################################################################################
# Set the overall range cut
#/run/setCut 1 mm
# Set the range cut on a particular volumes (in mm)
#/musr/command SetUserLimits log_target 0.01
#/musr/command SetUserLimits log_targetscint 0.01
#/musr/command SetUserLimits log_M0 0.01
# Store all events into the ROOT tree or just the interesting ones ? (true is default)
/musr/command storeOnlyEventsWithHits false
# Set the minimum time separation between two subsequent signals in the same detector (in ns)
/musr/command signalSeparationTime 10
#
#/musr/run/howOftenToPrintEvent 10
/musr/command maximumRunTimeAllowed 220000
/musr/run/randomOption 2
###################################################################################
######################### R O O T O U T P U T ##############################
###################################################################################
/musr/command rootOutput det_VvvKine off
/musr/command rootOutput det_VvvX off
/musr/command rootOutput det_VvvY off
/musr/command rootOutput det_VvvZ off
/musr/command rootOutput det_VvvVolID off
/musr/command rootOutput det_VvvProcID off
/musr/command rootOutput det_VvvTrackID off
/musr/command rootOutput det_VvvParticleID off
/musr/command rootOutput nOptPhot off
/musr/command rootOutput odet_ID off
/musr/command rootOutput odet_nPhot off
/musr/command rootOutput odet_timeFirst off
/musr/command rootOutput odet_timeA off
/musr/command rootOutput odet_timeB off
/musr/command rootOutput odet_timeC off
/musr/command rootOutput odet_timeD off
/musr/command rootOutput odet_timeE off
/musr/command rootOutput odet_timeLast off
###################################################################################
######################### V I S U A L I S A T I O N ##############################
###################################################################################
/vis/disable
#/control/execute visFromToni.mac
#/control/execute visDawn101.mac
#/control/execute visDawn40123.mac
#/control/execute visVRML.mac
###################################################################################
######################### P A R T I C L E G U N #################################
###################################################################################
/gun/vertex 0 0 -1000 mm
# FWHM 10mm ==> sigma = 10/2.354 = 4.2481mm
/gun/vertexsigma 25 25 0 mm
#---/gun/vertexboundary: rMaxAllowed, zMinAllowed, zMaxAllowed
/gun/vertexboundary 100 -999999 999999 mm
/gun/momentum 100.0 MeV
# sigma = 3% ==> sigma 27*0.03 = 0.81
/gun/momentumsmearing 3.00 MeV
#---/gun/momentumboundary: pMinAllowed, pMaxAllowed, dummy
#/gun/momentumboundary 20 40 0 MeV
# TURTLE
#/gun/turtlefilename FOR061_2008_04_22.DAT
#/gun/turtlefilename FOR070_2008_10_17_XXII.DAT
#/gun/turtlefilename FOR077_pie3_HiField_d05_x30.dat
#/gun/turtlefilename FOR077_reggiani_Jan2010_NEW.dat
#/gun/turtleZ0position -900 mm
#/gun/turtlefilename FOR070_GPD_beam_profile_45mm.dat
#/gun/turtleZ0position -500 mm
#---/gun/tilt: xangle, yangle, dummy
#/gun/tilt 0 0.5 0 deg
#---/gun/tiltsigma: xangleSigma, yangleSigma, dummy (1 degree on 1 meter ~ 17mm)
/gun/tiltsigma 0.15 0.15 0 deg
/gun/pitch 0.08 deg
/gun/muonPolarizVector 0 0 -1
#---/gun/decaytimelimits: decayMin, decayMax, decayTime
/gun/decaytimelimits -1 -1 2197.03 ns
###################################################################################
######################## B E A M O N #########################################
###################################################################################
/run/beamOn 16000000
#/run/beamOn 1

174
run/GPD/40602f.v1190 Normal file
View File

@ -0,0 +1,174 @@
# TDC V1190 Set up file generated by Deltat : 16-Apr-2008 10:30:27
# bins are always in TDC channel units (195.3125 ps/channel)
INSTRUMENT=GPS
DESCRIPTION=No Veto -- 1port
TYPE=TDCV1190
RESOLUTION=100
MDELAY=0
PDELAY=2000
#REWINDTIMEBINS=10000000
MCOINCIDENCEW=50
PCOINCIDENCEW=50
VCOINCIDENCEW=100
# 30000 mu/s
MUONRATEFACTOR=0.0909316
DATAWINDOWMIN=-5
DATAWINDOWMAX=10.0
PILEUPWINDOWMIN=-10.2
PILEUPWINDOWMAX=10.2
102; "M up"; M; 0.4; 2005; ;
1; "B1"; P; 0.4; 2005; ; B1; 1485; 1515; 50995;
2; "B2"; P; 0.4; 2005; ; B2; 1485; 1515; 50995;
4; "B4"; P; 0.4; 2005; ; B4; 1485; 1515; 50995;
5; "B5"; P; 0.4; 2005; ; B5; 1485; 1515; 50995;
6; "B6"; P; 0.4; 2005; ; B6; 1485; 1515; 50995;
7; "B6"; P; 0.4; 2005; ; B7; 1485; 1515; 50995;
11; "F1"; P; 0.4; 2005; ; F11; 1485; 1515; 50995;
12; "F2"; P; 0.4; 2005; ; F12; 1485; 1515; 50995;
13; "F3"; P; 0.4; 2005; ; F13; 1485; 1515; 50995;
14; "F4"; P; 0.4; 2005; ; F14; 1485; 1515; 50995;
15; "F5"; P; 0.4; 2005; ; F15; 1485; 1515; 50995;
!103; "M embed" V; 0.1; 2005;
! 21; "Coinc B1" K; 0.4; 2005;
! 22; "Coinc B2" K; 0.4; 2005;
!401; "Act Veto" V; 0.1; 2005;
! 3; "F right"; K; 0.4; 1263;
! 4; "B left"; P; 0.2; 1311; 5 -10 -11 -12 -13; Back; 1485; 1515; 50995;
! 5; "B right"; K; 0.2; 1262;
! 6; Up; P; 0.2; 1306; -10 -11 -12 -13; Up; 1490; 1520; 51000;
! 7; Down; P; 0.2; 1264; -10 -11 -12 -13; Down; 1525; 1555; 51035;
! 8; "R int"; P; 0.2; 1309; -10 -11 -12 -13; R_int; 1480; 1510; 50990;
! 9; "R ext"; P; 0.3; 1267; -10 -11 -12 -13; R_ext; 1520; 1550; 51020; 8; Rite
!10; "Bveto up"; V; 0.3; 2085;
!11; "Bveto down"; V; 0.3; 2089;
!12; "Bveto left"; V; 0.3; 2080;
!13; "Bveto right"; V; 0.3; 2079;
!14; "F center"; P; 0.2; 2073; -10 -11 -12 -13; F_cntr; 720; 750; 50230; 2; Forw
!15; Clock; C;
!16; Rejected; R;
!17; Ch17; N;
!18; Ch18; N;
!19; Ch19; N;
!20; Ch20; N;
!21; Ch21; N;
!22; Ch22; N;
!23; Ch23; N;
!24; Ch24; N;
!25; Ch25; N;
!26; Ch26; N;
!27; Ch27; N;
!28; Ch28; N;
!29; Ch29; N;
!30; Ch30; N;
!31; Ch31; N;
#DEBUGEVENT 704 10
#DEBUGEVENT 1427 10
#DEBUGEVENT 3222 10
#DEBUGEVENT 5223 10
#DEBUGEVENT 8351 10
#DEBUGEVENT 8581 10
#DEBUGEVENT 9041 10
#DEBUGEVENT 9645 10
counterPhaseShifts 1 0 2 0 11 180 12 180 13 180
counterGrouping B 1 2
counterGrouping F 11 12 13
counterGrouping L 4 5
counterGrouping R 14 15
musrTH2D hMuDecayMap "Muon decay map;z [mm]; r[mm]" 100 -80. 20. 40 0. 20. muDecayPosZ muDecayPosR
musrTH2D hMuDecayMappileup "Pileup muon decay map;z [mm]; r[mm]" 400 -300. 100. 200 0. 200. pileup_muDecayPosZ pileup_muDecayPosR
musrTH1D hMuDecayPosZ "Penetration of muons into the sample;z[mm];N" 100 -5.0 5. muDecayPosZ
musrTH1D hMuDecayPosZpileup "Decay of the muons contributing to pileup bgr.;z[mm];N" 400 -300 100 pileup_muDecayPosZ
musrTH1D hMuDecayPosX "X of decayed muons;x[mm];N" 100 -5.0 5. muDecayPosX
musrTH1D hdet_m0edep "Energy deposited in the trigger;E[MeV];N" 100 0. 2 det_m0edep
musrTH1D hdet_posEdep "Energy deposited in the positron Counters;E[MeV];N" 100 0. 2 det_posEdep
#musrTH1D hdet_time0 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time0
#musrTH1D hdet_time1 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time1
#musrTH1D hdet_time10 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time10
musrTH1D hgen_time10 "generated muSR spectra;time[#mus];N" 220 -11. 11. gen_time10
musrTH1D hdet_time10_MINUS_gen_time10 "detected - generated ;time_{det}-time_{gen} [ps];N" 100 -500. 500. det_time10_MINUS_gen_time10
#musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 4000 -2000. 2000. muDecayDetID
musrTH1D hdet_time10 "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time10
musrTH1D hdet_time10_Bgr "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time10
musrTH1D hdet_time20 "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time20
musrTH1D hdet_time20_Bgr "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time20
musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 2002 -1001. 1001. muDecayDetID
musrTH1D hMuDecayDetIDpileup "Detector ID where the pileup muons decay;Detector ID; N" 2002 -1001. 1001. pileup_muDecayDetID
humanDecayHistograms hMuDecayDetID hMuDecayDetIDpileup 201 Sample 211 SampleCell 300 Holder 301 Holder 302 Holder 303 Holder 304 Holder 305 Holder 306 Holder 307 Holder 308 Holder 309 Holder 310 Holder 313 Holder 314 Holder 315 Holder 316 Holder 251 Col_Cu 255 Col_1 257 Col_1 261 Col_2 263 Col_2 265 Al_Plate 102 M0 103 M0Holder 1 BC 2 BC 11 FC 12 FC 13 FC 21 BBC 22 BBC 401 Veto -1 World 291 Magnet 292 Magnet 293 Al_Uprof 230 Elsewhere 231 Elsewhere 232 Elsewhere 233 Elsewhere 256 Elsewhere 262 Elsewhere 267 Elsewhere -1000 Escaped
condition 0 alwaysTrue
condition 1 oncePerEvent
condition 2 muonDecayedInSample_gen
condition 3 muonTriggered_gen
condition 4 muonTriggered_det
condition 5 positronHit_det
condition 6 goodEvent_det
condition 7 goodEvent_gen
condition 8 goodEvent_det_AND_goodEvent_gen
condition 9 pileupEvent
condition 10 goodEvent_F_det
condition 11 goodEvent_B_det
condition 12 goodEvent_F_det_AND_pileupEvent
condition 13 goodEvent_B_det_AND_pileupEvent
condition 14 goodEvent_L_det
condition 15 goodEvent_R_det
fit hdet_time20_Bgr pol0 111 -1.9 -0.1 0
#fit hdet_time20 funct1 111 0.2 9.8 0 0 0.3 100
#fit hdet_time20 funct2 111 0.2 9.8 0 0 0.3 100 0
fit hdet_time20 funct4 111 0.2 9.8 0 0 0.3 100 0
#draw hMuDecayPosZ 0
#draw hMuDecayPosZ 1
#draw hMuDecayPosZ 2
#draw hMuDecayPosZ 3
#draw hMuDecayPosZ 4
###draw hMuDecayMap 0
###draw hMuDecayMap 1
#draw hdet_m0edep 0
#draw hdet_m0edep 1
#draw hdet_m0edep 2
#draw hdet_m0edep 3
#draw hdet_m0edep 4
#draw hdet_m0edep 5
#draw hdet_m0edep 6
#draw hdet_time10 6
draw hdet_time20 10
draw hdet_time20 11
draw hdet_time20 14
draw hdet_time20 15
draw hdet_time20_Bgr 10
draw hdet_time20_Bgr 11
draw hdet_time20_Bgr 14
draw hdet_time20_Bgr 15
#draw hgen_time10 7
#draw hdet_time10_MINUS_gen_time10 8
#draw hdet_posEdep 5
#draw hMuDecayDetID 1
#draw hdet_m0edep 9
#draw hdet_time10 9
#draw hMuDecayDetIDpileup 1
#draw hMuDecayDetIDpileup 6
#draw hMuDecayDetIDpileup 9
#draw humanDecayHistograms 1
#draw humanDecayHistograms 6
#draw humanDecayHistograms 9
#draw humanDecayPileupHistograms 1
#draw humanDecayPileupHistograms 6
#draw humanDecayPileupHistograms 9
#draw hdet_time20 4
#draw hdet_time20 6
$

263
run/GPD/40603.mac Normal file
View File

@ -0,0 +1,263 @@
#-----------------------------------------------------------------------
# Macro file for muSR instrument simulations
# Unless specified otherwises, the default units are mm, ns, MeV, MeV/c.
# Lines starting with star "#" are comments.
###################################################################################
############################# G E O M E T R Y ###################################
###################################################################################
# ROTATION MATRIXES:
/musr/command rotation matrix1 1 0 0 90
/musr/command rotation matrix2 0 1 0 28.81
/musr/command rotation matrix3 0 1 0 151.19
/musr/command rotation matrix4 0 0 1 45
/musr/command rotation matrix5 0 1 0 90
/musr/command rotation matrix6 180 90 90
/musr/command rotation matrix7 0 90 90
# --- WORLD VOLUME ---
/musr/command construct box World 2100 2100 11000 G4_AIR 0 0 0 no_logical_volume norot dead -1
# --- BEAMPIPE ---
/musr/command construct tubs beampipeA 0 2005 4500 0 360 G4_Galactic 0 0 -6100 log_World norot dead 230
/musr/command construct tubs beampipe 0 405 500 0 360 G4_Galactic 0 0 -1100 log_World norot dead 231
/musr/command construct tubs beampipeshell 100 105 500 0 360 G4_Al 0 0 0 log_beampipe norot dead 232
/musr/command construct tubs beampipewindow 0 100 0.05 0 360 G4_MYLAR 0 0 499.9 log_beampipe norot dead 233
# PASSIVE COLLIMATORS (for BEAMPIPE vers.1)
# GPDcollimator - thickness in the thinner location is x2+(x9-x4) = 20+59-59
##/musr/command construct GPDcollimator GPDcollimator 40 15 60 59 60.01 0 360 0 59 0 2.5 6 G4_Cu 0 0 -59 log_World matrix1 dead 251
/musr/command construct box GPDcollimator2 48 55 15 G4_Pb 0 0 -66 log_World norot dead 255
/musr/command construct box GPDcollimator2hole 2.5 4 15 G4_Galactic 0 0 0 log_GPDcollimator2 norot dead 256
/musr/command construct box FirstCollimator 45 75 15 G4_Pb 0 0 -250 log_World norot dead 261
/musr/command construct tubs FirstCollHole 0 8 15 0 360 G4_Galactic 0 0 0 log_FirstCollimator norot dead 262
/musr/command construct box SplitCollimator1 19 75 30 G4_Cu 23 0 -163 log_World norot dead 261
/musr/command construct box SplitCollimator2 19 75 30 G4_Cu -23 0 -163 log_World norot dead 261
###/musr/command construct box AlumPlate 30 90 5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -80 log_World norot musr/ScintSD 401
###/musr/command construct box AlumPlateHole 2 5 5 G4_Galactic 0 0 0 log_AlumPlate norot dead 267
/musr/command construct box Insert 20 20 10 G4_Pb 0 0 -91 log_World norot dead 265
/musr/command construct box InsertHole 2.5 4 10 G4_Galactic 0 0 0 log_Insert norot dead 267
# --- SAMPLE HOLDER ---
##/musr/command construct GPDsampleHolderA GPDsampleHolderA 0 0 0 0 0 G4_Al 0 -111.25 0 log_World matrix1 dead 300
##/musr/command construct tubs GPDsampleHolderB 37.51 58 7 0 360 G4_Al 0 57 0 log_World matrix1 dead 301
##/musr/command construct tubs GPDsampleHolderC 37.51 58 7 0 360 G4_Al 0 -57 0 log_World matrix1 dead 302
######/musr/command construct tubs GPDsampleHolderD 37 37.5 50 136 88 G4_Al 0 0 0 log_World matrix1 dead 303
######/musr/command construct tubs GPDsampleHolderE 53.5 54 50 136 88 G4_Al 0 0 0 log_World matrix1 dead 304
##/musr/command construct box GPDsampleHolderD 12 0.5 50 G4_Al 38 0 28 log_World matrix1 dead 303
##/musr/command construct box GPDsampleHolderE 12 0.5 50 G4_Al 38 0 -28 log_World matrix1 dead 304
##/musr/command construct box GPDsampleHolderF 12 0.5 50 G4_Al -38 0 28 log_World matrix1 dead 305
##/musr/command construct box GPDsampleHolderG 12 0.5 50 G4_Al -38 0 -28 log_World matrix1 dead 306
##/musr/command construct tubs GPDsampleHolderH 0 9 50 180 50 G4_Al 48 0 29 log_World matrix1 dead 307
##/musr/command construct tubs GPDsampleHolderI 0 9 50 130 50 G4_Al 48 0 -29 log_World matrix1 dead 308
##/musr/command construct tubs GPDsampleHolderJ 0 9 50 310 50 G4_Al -48 0 29 log_World matrix1 dead 309
##/musr/command construct tubs GPDsampleHolderK 0 9 50 0 50 G4_Al -48 0 -29 log_World matrix1 dead 310
##/musr/command construct tubs GPDsampleHolderW 58.01 58.2 64 25 130 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 313
##/musr/command construct tubs GPDsampleHolderX 58.01 58.2 64 205 130 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 314
##/musr/command construct tubs GPDsampleHolderY 36.79 36.99 64 45 90 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 315
##/musr/command construct tubs GPDsampleHolderZ 36.79 36.99 64 225 90 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 316
# --- SAMPLE CELL ---
/musr/command construct tubs GPDsampleCell 0 12 50 0 360 G4_Cu 0 0 0 log_World matrix1 dead 211
# --- SAMPLE ---
/musr/command construct tubs target 0 3 5 0 360 G4_Cu 0 0 0 log_GPDsampleCell norot dead 201
#/musr/command construct tubs targetFieldVol 0 3 5 0 360 G4_Cu 0 0 0 log_GPDsampleCell norot dead 201
# --- SCINTILLATOR ---
/musr/command construct box ScintB1 21.5 50 2 G4_PLASTIC_SC_VINYLTOLUENE 26.5 0 -44 log_World norot musr/ScintSD 1
/musr/command construct box ScintB2 21.5 50 2 G4_PLASTIC_SC_VINYLTOLUENE -26.5 0 -44 log_World norot musr/ScintSD 2
/musr/command construct box ScintB4 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE 44 0 -25 log_World norot musr/ScintSD 4
/musr/command construct box ScintB5 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE -44 0 -25 log_World norot musr/ScintSD 5
/musr/command construct box ScintB6 2 50 13 G4_PLASTIC_SC_VINYLTOLUENE 44 0 0 log_World norot musr/ScintSD 6
/musr/command construct box ScintB7 2 50 13 G4_PLASTIC_SC_VINYLTOLUENE -44 0 0 log_World norot musr/ScintSD 7
/musr/command construct box ScintF1 20 50 2 G4_PLASTIC_SC_VINYLTOLUENE 28 0 44 log_World norot musr/ScintSD 11
/musr/command construct box ScintF2 20 50 2 G4_PLASTIC_SC_VINYLTOLUENE -28 0 44 log_World norot musr/ScintSD 12
/musr/command construct box ScintF3 7 50 2 G4_PLASTIC_SC_VINYLTOLUENE 0 0 44 log_World norot musr/ScintSD 13
/musr/command construct box ScintF4 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE 44 0 25 log_World norot musr/ScintSD 14
/musr/command construct box ScintF5 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE -44 0 25 log_World norot musr/ScintSD 15
#/musr/command construct box KoincB1 20 90 2 G4_PLASTIC_SC_VINYLTOLUENE 33 0 -60 log_World matrix2 musr/ScintSD 21
#/musr/command construct box KoincB2 20 90 2 G4_PLASTIC_SC_VINYLTOLUENE -33 0 -60 log_World matrix3 musr/ScintSD 22
# --- M0 ---
/musr/command construct box M0holder 3 4 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -20 log_World norot musr/ScintSD 103
/musr/command construct box M0 4 30 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -44 log_World norot musr/ScintSD 102
#/musr/command construct GPDmHolder GPDmHolder 5 25 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 -25 -43 log_World norot dead 103
# --- Magnet ---
/musr/command construct tubs Magnet1 180 375 78 0 360 G4_Cu 133 0 0 log_World matrix5 dead 291
/musr/command construct tubs Magnet2 180 375 78 0 360 G4_Cu -133 0 0 log_World matrix5 dead 292
/musr/command construct uprofile Uprof1 25 10 400 5 G4_Al -45 0 200 log_World matrix6 dead 293
/musr/command construct uprofile Uprof2 25 10 400 5 G4_Al -45 0 -200 log_World matrix6 dead 293
/musr/command construct uprofile Uprof3 25 10 400 5 G4_Al 45 0 200 log_World matrix7 dead 293
/musr/command construct uprofile Uprof4 25 10 400 5 G4_Al 45 0 -200 log_World matrix7 dead 293
#
##---------------------------------------------------------
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_1 6
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_2 5
#/musr/command logicalVolumeToBeReweighted mu log_activeCollim 4
##---------------------------------------------------------
/musr/command visattributes log_World invisible
/musr/command visattributes log_beampipeA invisible
/musr/command visattributes log_beampipe invisible
/musr/command visattributes log_beampipeshell invisible
/musr/command visattributes log_beampipewindow invisible
/musr/command visattributes log_Magnet1 invisible
/musr/command visattributes log_Magnet2 invisible
/musr/command visattributes log_Uprof1 invisible
/musr/command visattributes log_Uprof2 invisible
/musr/command visattributes log_Uprof3 invisible
/musr/command visattributes log_Uprof4 invisible
/musr/command visattributes log_GPDsampleHolderD red
/musr/command visattributes log_GPDsampleHolderE red
/musr/command visattributes log_GPDsampleHolderF red
/musr/command visattributes log_GPDsampleHolderG red
/musr/command visattributes log_GPDcollimator green
/musr/command visattributes log_GPDcollimator2 blue
/musr/command visattributes log_GPDmHolder magenta
/musr/command visattributes log_target blue
/musr/command visattributes log_M0 yellow
#/musr/command visattributes G4_PLASTIC_SC_VINYLTOLUENE red
#/musr/command visattributes G4_Pb blue
####################################################################################
######################### M A G N E T I C F I E L D #########################
###################################################################################
# Set magnetic field (set field intensity in T and sigma in mm)
# syntax for magneticfield: fromfile filename fieldValue
# uniform fieldValue
# gaussian fieldValue sigma
#
/musr/command globalfield centralSolenoidField 0. 0. 0. fromfile 2D GPD_coil_primitive.table log_target 0.03
/musr/command globalfield setparameter SetLargestAcceptableStep 200
/musr/command globalfield setparameter SetMinimumEpsilonStep 5e-5
/musr/command globalfield setparameter SetMaximumEpsilonStep 0.001
/musr/command globalfield setparameter SetDeltaOneStep 0.1
/musr/command globalfield setparameter SetDeltaIntersection 0.01
/musr/command globalfield printparameters
/musr/command globalfield printFieldValueAtPoint 0 0 0
/musr/command globalfield printFieldValueAtPoint 0 0 10
/musr/command globalfield printFieldValueAtPoint 0 0 20
/musr/command globalfield printFieldValueAtPoint 0 0 30
/musr/command globalfield printFieldValueAtPoint 0 0 50
/musr/command globalfield printFieldValueAtPoint 0 0 100
/musr/command globalfield printFieldValueAtPoint 0 0 200
/musr/command globalfield printFieldValueAtPoint 0 0 300
/musr/command globalfield printFieldValueAtPoint 0 0 400
/musr/command globalfield printFieldValueAtPoint 0 0 500
/musr/command globalfield printFieldValueAtPoint 10 0 0
/musr/command globalfield printFieldValueAtPoint 100 0 0
/musr/command globalfield printFieldValueAtPoint 300 0 0
/musr/command globalfield printFieldValueAtPoint 0 10 0
/musr/command globalfield printFieldValueAtPoint 0 100 0
/musr/command globalfield printFieldValueAtPoint 0 300 0
###################################################################################
######################### P H Y S I C S P R O C E S S E S ##################
###################################################################################
# Geant 4.9.3
/musr/command process addDiscreteProcess gamma G4PhotoElectricEffect
/musr/command process addDiscreteProcess gamma G4ComptonScattering
/musr/command process addDiscreteProcess gamma G4GammaConversion
/musr/command process addDiscreteProcess gamma G4RayleighScattering
/musr/command process addProcess e- G4eMultipleScattering -1 1 1
/musr/command process addProcess e- G4eIonisation -1 2 2
/musr/command process addProcess e- G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eMultipleScattering -1 1 1
/musr/command process addProcess e+ G4eIonisation -1 2 2
/musr/command process addProcess e+ G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eplusAnnihilation 0 -1 4
/musr/command process addProcess mu- G4MuMultipleScattering -1 1 1
/musr/command process addProcess mu- G4MuIonisation -1 2 2
/musr/command process addProcess mu- G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu- G4MuPairProduction -1 4 4
/musr/command process addProcess mu+ G4MuMultipleScattering -1 1 1
/musr/command process addProcess mu+ G4MuIonisation -1 2 2
/musr/command process addProcess mu+ G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu+ G4MuPairProduction -1 4 4
#
#/home/install/geant4.9.3/source/processes/electromagnetic/utils/src/G4EnergyLossMessenger.cc
######## /process/msc/StepLimit Minimal | UseDistanceToBoundary | UseSafety
/process/msc/StepLimit UseSafety
#/process/msc/LateralDisplacement
/process/msc/RangeFactor 0.04
/process/msc/GeomFactor 2.5
#/process/msc/FactorForAngleLimit
/process/msc/Skin 3.0
#/process/msc/ThetaLimit 0.2 rad
###################################################################################
################## S O M E O T H E R P A R A M E T E R S ##################
###################################################################################
# Set the overall range cut
#/run/setCut 1 mm
# Set the range cut on a particular volumes (in mm)
#/musr/command SetUserLimits log_target 0.01
#/musr/command SetUserLimits log_targetscint 0.01
#/musr/command SetUserLimits log_M0 0.01
# Store all events into the ROOT tree or just the interesting ones ? (true is default)
/musr/command storeOnlyEventsWithHits false
# Set the minimum time separation between two subsequent signals in the same detector (in ns)
/musr/command signalSeparationTime 10
#
#/musr/run/howOftenToPrintEvent 10
/musr/command maximumRunTimeAllowed 220000
/musr/run/randomOption 2
###################################################################################
######################### R O O T O U T P U T ##############################
###################################################################################
/musr/command rootOutput det_VvvKine off
/musr/command rootOutput det_VvvX off
/musr/command rootOutput det_VvvY off
/musr/command rootOutput det_VvvZ off
/musr/command rootOutput det_VvvVolID off
/musr/command rootOutput det_VvvProcID off
/musr/command rootOutput det_VvvTrackID off
/musr/command rootOutput det_VvvParticleID off
/musr/command rootOutput nOptPhot off
/musr/command rootOutput odet_ID off
/musr/command rootOutput odet_nPhot off
/musr/command rootOutput odet_timeFirst off
/musr/command rootOutput odet_timeA off
/musr/command rootOutput odet_timeB off
/musr/command rootOutput odet_timeC off
/musr/command rootOutput odet_timeD off
/musr/command rootOutput odet_timeE off
/musr/command rootOutput odet_timeLast off
###################################################################################
######################### V I S U A L I S A T I O N ##############################
###################################################################################
/vis/disable
#/control/execute visFromToni.mac
#/control/execute visDawn101.mac
#/control/execute visDawn40123.mac
#/control/execute visVRML.mac
###################################################################################
######################### P A R T I C L E G U N #################################
###################################################################################
/gun/vertex 0 0 -1000 mm
# FWHM 10mm ==> sigma = 10/2.354 = 4.2481mm
/gun/vertexsigma 25 25 0 mm
#---/gun/vertexboundary: rMaxAllowed, zMinAllowed, zMaxAllowed
/gun/vertexboundary 100 -999999 999999 mm
/gun/momentum 100.0 MeV
# sigma = 3% ==> sigma 27*0.03 = 0.81
/gun/momentumsmearing 3.00 MeV
#---/gun/momentumboundary: pMinAllowed, pMaxAllowed, dummy
#/gun/momentumboundary 20 40 0 MeV
# TURTLE
#/gun/turtlefilename FOR061_2008_04_22.DAT
#/gun/turtlefilename FOR070_2008_10_17_XXII.DAT
#/gun/turtlefilename FOR077_pie3_HiField_d05_x30.dat
#/gun/turtlefilename FOR077_reggiani_Jan2010_NEW.dat
#/gun/turtleZ0position -900 mm
#/gun/turtlefilename FOR070_GPD_beam_profile_45mm.dat
#/gun/turtleZ0position -500 mm
#---/gun/tilt: xangle, yangle, dummy
#/gun/tilt 0 0.5 0 deg
#---/gun/tiltsigma: xangleSigma, yangleSigma, dummy (1 degree on 1 meter ~ 17mm)
/gun/tiltsigma 0.15 0.15 0 deg
/gun/pitch 0.08 deg
/gun/muonPolarizVector 0 0 -1
#---/gun/decaytimelimits: decayMin, decayMax, decayTime
/gun/decaytimelimits -1 -1 2197.03 ns
###################################################################################
######################## B E A M O N #########################################
###################################################################################
/run/beamOn 16000000
#/run/beamOn 1

263
run/GPD/40604.mac Normal file
View File

@ -0,0 +1,263 @@
#-----------------------------------------------------------------------
# Macro file for muSR instrument simulations
# Unless specified otherwises, the default units are mm, ns, MeV, MeV/c.
# Lines starting with star "#" are comments.
###################################################################################
############################# G E O M E T R Y ###################################
###################################################################################
# ROTATION MATRIXES:
/musr/command rotation matrix1 1 0 0 90
/musr/command rotation matrix2 0 1 0 28.81
/musr/command rotation matrix3 0 1 0 151.19
/musr/command rotation matrix4 0 0 1 45
/musr/command rotation matrix5 0 1 0 90
/musr/command rotation matrix6 180 90 90
/musr/command rotation matrix7 0 90 90
# --- WORLD VOLUME ---
/musr/command construct box World 2100 2100 11000 G4_AIR 0 0 0 no_logical_volume norot dead -1
# --- BEAMPIPE ---
/musr/command construct tubs beampipeA 0 2005 4500 0 360 G4_Galactic 0 0 -6100 log_World norot dead 230
/musr/command construct tubs beampipe 0 405 500 0 360 G4_Galactic 0 0 -1100 log_World norot dead 231
/musr/command construct tubs beampipeshell 100 105 500 0 360 G4_Al 0 0 0 log_beampipe norot dead 232
/musr/command construct tubs beampipewindow 0 100 0.05 0 360 G4_MYLAR 0 0 499.9 log_beampipe norot dead 233
# PASSIVE COLLIMATORS (for BEAMPIPE vers.1)
# GPDcollimator - thickness in the thinner location is x2+(x9-x4) = 20+59-59
##/musr/command construct GPDcollimator GPDcollimator 40 15 60 59 60.01 0 360 0 59 0 2.5 6 G4_Cu 0 0 -59 log_World matrix1 dead 251
/musr/command construct box GPDcollimator2 48 55 15 G4_Pb 0 0 -66 log_World norot dead 255
/musr/command construct box GPDcollimator2hole 2.5 4 15 G4_Galactic 0 0 0 log_GPDcollimator2 norot dead 256
/musr/command construct box FirstCollimator 45 75 15 G4_Pb 0 0 -250 log_World norot dead 261
/musr/command construct tubs FirstCollHole 0 8 15 0 360 G4_Galactic 0 0 0 log_FirstCollimator norot dead 262
/musr/command construct box SplitCollimator1 19 75 20 G4_Cu 23 0 -153 log_World norot dead 261
/musr/command construct box SplitCollimator2 19 75 20 G4_Cu -23 0 -153 log_World norot dead 261
###/musr/command construct box AlumPlate 30 90 5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -80 log_World norot musr/ScintSD 401
###/musr/command construct box AlumPlateHole 2 5 5 G4_Galactic 0 0 0 log_AlumPlate norot dead 267
/musr/command construct box Insert 20 20 10 G4_Pb 0 0 -91 log_World norot dead 265
/musr/command construct box InsertHole 2.5 4 10 G4_Galactic 0 0 0 log_Insert norot dead 267
# --- SAMPLE HOLDER ---
##/musr/command construct GPDsampleHolderA GPDsampleHolderA 0 0 0 0 0 G4_Al 0 -111.25 0 log_World matrix1 dead 300
##/musr/command construct tubs GPDsampleHolderB 37.51 58 7 0 360 G4_Al 0 57 0 log_World matrix1 dead 301
##/musr/command construct tubs GPDsampleHolderC 37.51 58 7 0 360 G4_Al 0 -57 0 log_World matrix1 dead 302
######/musr/command construct tubs GPDsampleHolderD 37 37.5 50 136 88 G4_Al 0 0 0 log_World matrix1 dead 303
######/musr/command construct tubs GPDsampleHolderE 53.5 54 50 136 88 G4_Al 0 0 0 log_World matrix1 dead 304
##/musr/command construct box GPDsampleHolderD 12 0.5 50 G4_Al 38 0 28 log_World matrix1 dead 303
##/musr/command construct box GPDsampleHolderE 12 0.5 50 G4_Al 38 0 -28 log_World matrix1 dead 304
##/musr/command construct box GPDsampleHolderF 12 0.5 50 G4_Al -38 0 28 log_World matrix1 dead 305
##/musr/command construct box GPDsampleHolderG 12 0.5 50 G4_Al -38 0 -28 log_World matrix1 dead 306
##/musr/command construct tubs GPDsampleHolderH 0 9 50 180 50 G4_Al 48 0 29 log_World matrix1 dead 307
##/musr/command construct tubs GPDsampleHolderI 0 9 50 130 50 G4_Al 48 0 -29 log_World matrix1 dead 308
##/musr/command construct tubs GPDsampleHolderJ 0 9 50 310 50 G4_Al -48 0 29 log_World matrix1 dead 309
##/musr/command construct tubs GPDsampleHolderK 0 9 50 0 50 G4_Al -48 0 -29 log_World matrix1 dead 310
##/musr/command construct tubs GPDsampleHolderW 58.01 58.2 64 25 130 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 313
##/musr/command construct tubs GPDsampleHolderX 58.01 58.2 64 205 130 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 314
##/musr/command construct tubs GPDsampleHolderY 36.79 36.99 64 45 90 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 315
##/musr/command construct tubs GPDsampleHolderZ 36.79 36.99 64 225 90 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 316
# --- SAMPLE CELL ---
/musr/command construct tubs GPDsampleCell 0 12 50 0 360 G4_Cu 0 0 0 log_World matrix1 dead 211
# --- SAMPLE ---
/musr/command construct tubs target 0 3 5 0 360 G4_Cu 0 0 0 log_GPDsampleCell norot dead 201
#/musr/command construct tubs targetFieldVol 0 3 5 0 360 G4_Cu 0 0 0 log_GPDsampleCell norot dead 201
# --- SCINTILLATOR ---
/musr/command construct box ScintB1 21.5 50 2 G4_PLASTIC_SC_VINYLTOLUENE 26.5 0 -44 log_World norot musr/ScintSD 1
/musr/command construct box ScintB2 21.5 50 2 G4_PLASTIC_SC_VINYLTOLUENE -26.5 0 -44 log_World norot musr/ScintSD 2
/musr/command construct box ScintB4 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE 44 0 -25 log_World norot musr/ScintSD 4
/musr/command construct box ScintB5 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE -44 0 -25 log_World norot musr/ScintSD 5
/musr/command construct box ScintB6 2 50 13 G4_PLASTIC_SC_VINYLTOLUENE 44 0 0 log_World norot musr/ScintSD 6
/musr/command construct box ScintB7 2 50 13 G4_PLASTIC_SC_VINYLTOLUENE -44 0 0 log_World norot musr/ScintSD 7
/musr/command construct box ScintF1 20 50 2 G4_PLASTIC_SC_VINYLTOLUENE 28 0 44 log_World norot musr/ScintSD 11
/musr/command construct box ScintF2 20 50 2 G4_PLASTIC_SC_VINYLTOLUENE -28 0 44 log_World norot musr/ScintSD 12
/musr/command construct box ScintF3 7 50 2 G4_PLASTIC_SC_VINYLTOLUENE 0 0 44 log_World norot musr/ScintSD 13
/musr/command construct box ScintF4 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE 44 0 25 log_World norot musr/ScintSD 14
/musr/command construct box ScintF5 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE -44 0 25 log_World norot musr/ScintSD 15
#/musr/command construct box KoincB1 20 90 2 G4_PLASTIC_SC_VINYLTOLUENE 33 0 -60 log_World matrix2 musr/ScintSD 21
#/musr/command construct box KoincB2 20 90 2 G4_PLASTIC_SC_VINYLTOLUENE -33 0 -60 log_World matrix3 musr/ScintSD 22
# --- M0 ---
/musr/command construct box M0holder 3 4 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -20 log_World norot musr/ScintSD 103
/musr/command construct box M0 4 30 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -44 log_World norot musr/ScintSD 102
#/musr/command construct GPDmHolder GPDmHolder 5 25 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 -25 -43 log_World norot dead 103
# --- Magnet ---
/musr/command construct tubs Magnet1 180 375 78 0 360 G4_Cu 133 0 0 log_World matrix5 dead 291
/musr/command construct tubs Magnet2 180 375 78 0 360 G4_Cu -133 0 0 log_World matrix5 dead 292
/musr/command construct uprofile Uprof1 25 10 400 5 G4_Al -45 0 200 log_World matrix6 dead 293
/musr/command construct uprofile Uprof2 25 10 400 5 G4_Al -45 0 -200 log_World matrix6 dead 293
/musr/command construct uprofile Uprof3 25 10 400 5 G4_Al 45 0 200 log_World matrix7 dead 293
/musr/command construct uprofile Uprof4 25 10 400 5 G4_Al 45 0 -200 log_World matrix7 dead 293
#
##---------------------------------------------------------
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_1 6
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_2 5
#/musr/command logicalVolumeToBeReweighted mu log_activeCollim 4
##---------------------------------------------------------
/musr/command visattributes log_World invisible
/musr/command visattributes log_beampipeA invisible
/musr/command visattributes log_beampipe invisible
/musr/command visattributes log_beampipeshell invisible
/musr/command visattributes log_beampipewindow invisible
/musr/command visattributes log_Magnet1 invisible
/musr/command visattributes log_Magnet2 invisible
/musr/command visattributes log_Uprof1 invisible
/musr/command visattributes log_Uprof2 invisible
/musr/command visattributes log_Uprof3 invisible
/musr/command visattributes log_Uprof4 invisible
/musr/command visattributes log_GPDsampleHolderD red
/musr/command visattributes log_GPDsampleHolderE red
/musr/command visattributes log_GPDsampleHolderF red
/musr/command visattributes log_GPDsampleHolderG red
/musr/command visattributes log_GPDcollimator green
/musr/command visattributes log_GPDcollimator2 blue
/musr/command visattributes log_GPDmHolder magenta
/musr/command visattributes log_target blue
/musr/command visattributes log_M0 yellow
#/musr/command visattributes G4_PLASTIC_SC_VINYLTOLUENE red
#/musr/command visattributes G4_Pb blue
####################################################################################
######################### M A G N E T I C F I E L D #########################
###################################################################################
# Set magnetic field (set field intensity in T and sigma in mm)
# syntax for magneticfield: fromfile filename fieldValue
# uniform fieldValue
# gaussian fieldValue sigma
#
/musr/command globalfield centralSolenoidField 0. 0. 0. fromfile 2D GPD_coil_primitive.table log_target 0.03
/musr/command globalfield setparameter SetLargestAcceptableStep 200
/musr/command globalfield setparameter SetMinimumEpsilonStep 5e-5
/musr/command globalfield setparameter SetMaximumEpsilonStep 0.001
/musr/command globalfield setparameter SetDeltaOneStep 0.1
/musr/command globalfield setparameter SetDeltaIntersection 0.01
/musr/command globalfield printparameters
/musr/command globalfield printFieldValueAtPoint 0 0 0
/musr/command globalfield printFieldValueAtPoint 0 0 10
/musr/command globalfield printFieldValueAtPoint 0 0 20
/musr/command globalfield printFieldValueAtPoint 0 0 30
/musr/command globalfield printFieldValueAtPoint 0 0 50
/musr/command globalfield printFieldValueAtPoint 0 0 100
/musr/command globalfield printFieldValueAtPoint 0 0 200
/musr/command globalfield printFieldValueAtPoint 0 0 300
/musr/command globalfield printFieldValueAtPoint 0 0 400
/musr/command globalfield printFieldValueAtPoint 0 0 500
/musr/command globalfield printFieldValueAtPoint 10 0 0
/musr/command globalfield printFieldValueAtPoint 100 0 0
/musr/command globalfield printFieldValueAtPoint 300 0 0
/musr/command globalfield printFieldValueAtPoint 0 10 0
/musr/command globalfield printFieldValueAtPoint 0 100 0
/musr/command globalfield printFieldValueAtPoint 0 300 0
###################################################################################
######################### P H Y S I C S P R O C E S S E S ##################
###################################################################################
# Geant 4.9.3
/musr/command process addDiscreteProcess gamma G4PhotoElectricEffect
/musr/command process addDiscreteProcess gamma G4ComptonScattering
/musr/command process addDiscreteProcess gamma G4GammaConversion
/musr/command process addDiscreteProcess gamma G4RayleighScattering
/musr/command process addProcess e- G4eMultipleScattering -1 1 1
/musr/command process addProcess e- G4eIonisation -1 2 2
/musr/command process addProcess e- G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eMultipleScattering -1 1 1
/musr/command process addProcess e+ G4eIonisation -1 2 2
/musr/command process addProcess e+ G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eplusAnnihilation 0 -1 4
/musr/command process addProcess mu- G4MuMultipleScattering -1 1 1
/musr/command process addProcess mu- G4MuIonisation -1 2 2
/musr/command process addProcess mu- G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu- G4MuPairProduction -1 4 4
/musr/command process addProcess mu+ G4MuMultipleScattering -1 1 1
/musr/command process addProcess mu+ G4MuIonisation -1 2 2
/musr/command process addProcess mu+ G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu+ G4MuPairProduction -1 4 4
#
#/home/install/geant4.9.3/source/processes/electromagnetic/utils/src/G4EnergyLossMessenger.cc
######## /process/msc/StepLimit Minimal | UseDistanceToBoundary | UseSafety
/process/msc/StepLimit UseSafety
#/process/msc/LateralDisplacement
/process/msc/RangeFactor 0.04
/process/msc/GeomFactor 2.5
#/process/msc/FactorForAngleLimit
/process/msc/Skin 3.0
#/process/msc/ThetaLimit 0.2 rad
###################################################################################
################## S O M E O T H E R P A R A M E T E R S ##################
###################################################################################
# Set the overall range cut
#/run/setCut 1 mm
# Set the range cut on a particular volumes (in mm)
#/musr/command SetUserLimits log_target 0.01
#/musr/command SetUserLimits log_targetscint 0.01
#/musr/command SetUserLimits log_M0 0.01
# Store all events into the ROOT tree or just the interesting ones ? (true is default)
/musr/command storeOnlyEventsWithHits false
# Set the minimum time separation between two subsequent signals in the same detector (in ns)
/musr/command signalSeparationTime 10
#
#/musr/run/howOftenToPrintEvent 10
/musr/command maximumRunTimeAllowed 220000
/musr/run/randomOption 2
###################################################################################
######################### R O O T O U T P U T ##############################
###################################################################################
/musr/command rootOutput det_VvvKine off
/musr/command rootOutput det_VvvX off
/musr/command rootOutput det_VvvY off
/musr/command rootOutput det_VvvZ off
/musr/command rootOutput det_VvvVolID off
/musr/command rootOutput det_VvvProcID off
/musr/command rootOutput det_VvvTrackID off
/musr/command rootOutput det_VvvParticleID off
/musr/command rootOutput nOptPhot off
/musr/command rootOutput odet_ID off
/musr/command rootOutput odet_nPhot off
/musr/command rootOutput odet_timeFirst off
/musr/command rootOutput odet_timeA off
/musr/command rootOutput odet_timeB off
/musr/command rootOutput odet_timeC off
/musr/command rootOutput odet_timeD off
/musr/command rootOutput odet_timeE off
/musr/command rootOutput odet_timeLast off
###################################################################################
######################### V I S U A L I S A T I O N ##############################
###################################################################################
/vis/disable
#/control/execute visFromToni.mac
#/control/execute visDawn101.mac
#/control/execute visDawn40123.mac
#/control/execute visVRML.mac
###################################################################################
######################### P A R T I C L E G U N #################################
###################################################################################
/gun/vertex 0 0 -1000 mm
# FWHM 10mm ==> sigma = 10/2.354 = 4.2481mm
/gun/vertexsigma 25 25 0 mm
#---/gun/vertexboundary: rMaxAllowed, zMinAllowed, zMaxAllowed
/gun/vertexboundary 100 -999999 999999 mm
/gun/momentum 100.0 MeV
# sigma = 3% ==> sigma 27*0.03 = 0.81
/gun/momentumsmearing 3.00 MeV
#---/gun/momentumboundary: pMinAllowed, pMaxAllowed, dummy
#/gun/momentumboundary 20 40 0 MeV
# TURTLE
#/gun/turtlefilename FOR061_2008_04_22.DAT
#/gun/turtlefilename FOR070_2008_10_17_XXII.DAT
#/gun/turtlefilename FOR077_pie3_HiField_d05_x30.dat
#/gun/turtlefilename FOR077_reggiani_Jan2010_NEW.dat
#/gun/turtleZ0position -900 mm
#/gun/turtlefilename FOR070_GPD_beam_profile_45mm.dat
#/gun/turtleZ0position -500 mm
#---/gun/tilt: xangle, yangle, dummy
#/gun/tilt 0 0.5 0 deg
#---/gun/tiltsigma: xangleSigma, yangleSigma, dummy (1 degree on 1 meter ~ 17mm)
/gun/tiltsigma 0.15 0.15 0 deg
/gun/pitch 0.08 deg
/gun/muonPolarizVector 0 0 -1
#---/gun/decaytimelimits: decayMin, decayMax, decayTime
/gun/decaytimelimits -1 -1 2197.03 ns
###################################################################################
######################## B E A M O N #########################################
###################################################################################
/run/beamOn 16000000
#/run/beamOn 1

174
run/GPD/40604f.v1190 Normal file
View File

@ -0,0 +1,174 @@
# TDC V1190 Set up file generated by Deltat : 16-Apr-2008 10:30:27
# bins are always in TDC channel units (195.3125 ps/channel)
INSTRUMENT=GPS
DESCRIPTION=No Veto -- 1port
TYPE=TDCV1190
RESOLUTION=100
MDELAY=0
PDELAY=2000
#REWINDTIMEBINS=10000000
MCOINCIDENCEW=50
PCOINCIDENCEW=50
VCOINCIDENCEW=100
# 30000 mu/s
MUONRATEFACTOR=0.0933042
DATAWINDOWMIN=-5
DATAWINDOWMAX=10.0
PILEUPWINDOWMIN=-10.2
PILEUPWINDOWMAX=10.2
102; "M up"; M; 0.4; 2005; ;
1; "B1"; P; 0.4; 2005; ; B1; 1485; 1515; 50995;
2; "B2"; P; 0.4; 2005; ; B2; 1485; 1515; 50995;
4; "B4"; P; 0.4; 2005; ; B4; 1485; 1515; 50995;
5; "B5"; P; 0.4; 2005; ; B5; 1485; 1515; 50995;
6; "B6"; P; 0.4; 2005; ; B6; 1485; 1515; 50995;
7; "B6"; P; 0.4; 2005; ; B7; 1485; 1515; 50995;
11; "F1"; P; 0.4; 2005; ; F11; 1485; 1515; 50995;
12; "F2"; P; 0.4; 2005; ; F12; 1485; 1515; 50995;
13; "F3"; P; 0.4; 2005; ; F13; 1485; 1515; 50995;
14; "F4"; P; 0.4; 2005; ; F14; 1485; 1515; 50995;
15; "F5"; P; 0.4; 2005; ; F15; 1485; 1515; 50995;
!103; "M embed" V; 0.1; 2005;
! 21; "Coinc B1" K; 0.4; 2005;
! 22; "Coinc B2" K; 0.4; 2005;
!401; "Act Veto" V; 0.1; 2005;
! 3; "F right"; K; 0.4; 1263;
! 4; "B left"; P; 0.2; 1311; 5 -10 -11 -12 -13; Back; 1485; 1515; 50995;
! 5; "B right"; K; 0.2; 1262;
! 6; Up; P; 0.2; 1306; -10 -11 -12 -13; Up; 1490; 1520; 51000;
! 7; Down; P; 0.2; 1264; -10 -11 -12 -13; Down; 1525; 1555; 51035;
! 8; "R int"; P; 0.2; 1309; -10 -11 -12 -13; R_int; 1480; 1510; 50990;
! 9; "R ext"; P; 0.3; 1267; -10 -11 -12 -13; R_ext; 1520; 1550; 51020; 8; Rite
!10; "Bveto up"; V; 0.3; 2085;
!11; "Bveto down"; V; 0.3; 2089;
!12; "Bveto left"; V; 0.3; 2080;
!13; "Bveto right"; V; 0.3; 2079;
!14; "F center"; P; 0.2; 2073; -10 -11 -12 -13; F_cntr; 720; 750; 50230; 2; Forw
!15; Clock; C;
!16; Rejected; R;
!17; Ch17; N;
!18; Ch18; N;
!19; Ch19; N;
!20; Ch20; N;
!21; Ch21; N;
!22; Ch22; N;
!23; Ch23; N;
!24; Ch24; N;
!25; Ch25; N;
!26; Ch26; N;
!27; Ch27; N;
!28; Ch28; N;
!29; Ch29; N;
!30; Ch30; N;
!31; Ch31; N;
#DEBUGEVENT 704 10
#DEBUGEVENT 1427 10
#DEBUGEVENT 3222 10
#DEBUGEVENT 5223 10
#DEBUGEVENT 8351 10
#DEBUGEVENT 8581 10
#DEBUGEVENT 9041 10
#DEBUGEVENT 9645 10
counterPhaseShifts 1 0 2 0 11 180 12 180 13 180
counterGrouping B 1 2
counterGrouping F 11 12 13
counterGrouping L 4 5
counterGrouping R 14 15
musrTH2D hMuDecayMap "Muon decay map;z [mm]; r[mm]" 100 -80. 20. 40 0. 20. muDecayPosZ muDecayPosR
musrTH2D hMuDecayMappileup "Pileup muon decay map;z [mm]; r[mm]" 400 -300. 100. 200 0. 200. pileup_muDecayPosZ pileup_muDecayPosR
musrTH1D hMuDecayPosZ "Penetration of muons into the sample;z[mm];N" 100 -5.0 5. muDecayPosZ
musrTH1D hMuDecayPosZpileup "Decay of the muons contributing to pileup bgr.;z[mm];N" 400 -300 100 pileup_muDecayPosZ
musrTH1D hMuDecayPosX "X of decayed muons;x[mm];N" 100 -5.0 5. muDecayPosX
musrTH1D hdet_m0edep "Energy deposited in the trigger;E[MeV];N" 100 0. 2 det_m0edep
musrTH1D hdet_posEdep "Energy deposited in the positron Counters;E[MeV];N" 100 0. 2 det_posEdep
#musrTH1D hdet_time0 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time0
#musrTH1D hdet_time1 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time1
#musrTH1D hdet_time10 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time10
musrTH1D hgen_time10 "generated muSR spectra;time[#mus];N" 220 -11. 11. gen_time10
musrTH1D hdet_time10_MINUS_gen_time10 "detected - generated ;time_{det}-time_{gen} [ps];N" 100 -500. 500. det_time10_MINUS_gen_time10
#musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 4000 -2000. 2000. muDecayDetID
musrTH1D hdet_time10 "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time10
musrTH1D hdet_time10_Bgr "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time10
musrTH1D hdet_time20 "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time20
musrTH1D hdet_time20_Bgr "detected muSR spectra (phase shifted);time[#mus];N" 700 -2.0 12.0 det_time20
musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 2002 -1001. 1001. muDecayDetID
musrTH1D hMuDecayDetIDpileup "Detector ID where the pileup muons decay;Detector ID; N" 2002 -1001. 1001. pileup_muDecayDetID
humanDecayHistograms hMuDecayDetID hMuDecayDetIDpileup 201 Sample 211 SampleCell 300 Holder 301 Holder 302 Holder 303 Holder 304 Holder 305 Holder 306 Holder 307 Holder 308 Holder 309 Holder 310 Holder 313 Holder 314 Holder 315 Holder 316 Holder 251 Col_Cu 255 Col_1 257 Col_1 261 Col_2 263 Col_2 265 Al_Plate 102 M0 103 M0Holder 1 BC 2 BC 11 FC 12 FC 13 FC 21 BBC 22 BBC 401 Veto -1 World 291 Magnet 292 Magnet 293 Al_Uprof 230 Elsewhere 231 Elsewhere 232 Elsewhere 233 Elsewhere 256 Elsewhere 262 Elsewhere 267 Elsewhere -1000 Escaped
condition 0 alwaysTrue
condition 1 oncePerEvent
condition 2 muonDecayedInSample_gen
condition 3 muonTriggered_gen
condition 4 muonTriggered_det
condition 5 positronHit_det
condition 6 goodEvent_det
condition 7 goodEvent_gen
condition 8 goodEvent_det_AND_goodEvent_gen
condition 9 pileupEvent
condition 10 goodEvent_F_det
condition 11 goodEvent_B_det
condition 12 goodEvent_F_det_AND_pileupEvent
condition 13 goodEvent_B_det_AND_pileupEvent
condition 14 goodEvent_L_det
condition 15 goodEvent_R_det
fit hdet_time20_Bgr pol0 111 -1.9 -0.1 0
#fit hdet_time20 funct1 111 0.2 9.8 0 0 0.3 100
#fit hdet_time20 funct2 111 0.2 9.8 0 0 0.3 100 0
fit hdet_time20 funct4 111 0.2 9.8 0 0 0.3 100 0
#draw hMuDecayPosZ 0
#draw hMuDecayPosZ 1
#draw hMuDecayPosZ 2
#draw hMuDecayPosZ 3
#draw hMuDecayPosZ 4
###draw hMuDecayMap 0
###draw hMuDecayMap 1
#draw hdet_m0edep 0
#draw hdet_m0edep 1
#draw hdet_m0edep 2
#draw hdet_m0edep 3
#draw hdet_m0edep 4
#draw hdet_m0edep 5
#draw hdet_m0edep 6
#draw hdet_time10 6
draw hdet_time20 10
draw hdet_time20 11
draw hdet_time20 14
draw hdet_time20 15
draw hdet_time20_Bgr 10
draw hdet_time20_Bgr 11
draw hdet_time20_Bgr 14
draw hdet_time20_Bgr 15
#draw hgen_time10 7
#draw hdet_time10_MINUS_gen_time10 8
#draw hdet_posEdep 5
#draw hMuDecayDetID 1
#draw hdet_m0edep 9
#draw hdet_time10 9
#draw hMuDecayDetIDpileup 1
#draw hMuDecayDetIDpileup 6
#draw hMuDecayDetIDpileup 9
#draw humanDecayHistograms 1
#draw humanDecayHistograms 6
#draw humanDecayHistograms 9
#draw humanDecayPileupHistograms 1
#draw humanDecayPileupHistograms 6
#draw humanDecayPileupHistograms 9
#draw hdet_time20 4
#draw hdet_time20 6
$

263
run/GPD/40605.mac Normal file
View File

@ -0,0 +1,263 @@
#-----------------------------------------------------------------------
# Macro file for muSR instrument simulations
# Unless specified otherwises, the default units are mm, ns, MeV, MeV/c.
# Lines starting with star "#" are comments.
###################################################################################
############################# G E O M E T R Y ###################################
###################################################################################
# ROTATION MATRIXES:
/musr/command rotation matrix1 1 0 0 90
/musr/command rotation matrix2 0 1 0 28.81
/musr/command rotation matrix3 0 1 0 151.19
/musr/command rotation matrix4 0 0 1 45
/musr/command rotation matrix5 0 1 0 90
/musr/command rotation matrix6 180 90 90
/musr/command rotation matrix7 0 90 90
# --- WORLD VOLUME ---
/musr/command construct box World 2100 2100 11000 G4_AIR 0 0 0 no_logical_volume norot dead -1
# --- BEAMPIPE ---
/musr/command construct tubs beampipeA 0 2005 4500 0 360 G4_Galactic 0 0 -6100 log_World norot dead 230
/musr/command construct tubs beampipe 0 405 500 0 360 G4_Galactic 0 0 -1100 log_World norot dead 231
/musr/command construct tubs beampipeshell 100 105 500 0 360 G4_Al 0 0 0 log_beampipe norot dead 232
/musr/command construct tubs beampipewindow 0 100 0.05 0 360 G4_MYLAR 0 0 499.9 log_beampipe norot dead 233
# PASSIVE COLLIMATORS (for BEAMPIPE vers.1)
# GPDcollimator - thickness in the thinner location is x2+(x9-x4) = 20+59-59
##/musr/command construct GPDcollimator GPDcollimator 40 15 60 59 60.01 0 360 0 59 0 2.5 6 G4_Cu 0 0 -59 log_World matrix1 dead 251
/musr/command construct box GPDcollimator2 48 55 15 G4_Pb 0 0 -66 log_World norot dead 255
/musr/command construct box GPDcollimator2hole 2.5 4 15 G4_Galactic 0 0 0 log_GPDcollimator2 norot dead 256
/musr/command construct box FirstCollimator 45 75 15 G4_Pb 0 0 -250 log_World norot dead 261
/musr/command construct tubs FirstCollHole 0 8 15 0 360 G4_Galactic 0 0 0 log_FirstCollimator norot dead 262
/musr/command construct box SplitCollimator1 19 75 20 G4_Pb 23 0 -153 log_World norot dead 261
/musr/command construct box SplitCollimator2 19 75 20 G4_Pb -23 0 -153 log_World norot dead 261
###/musr/command construct box AlumPlate 30 90 5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -80 log_World norot musr/ScintSD 401
###/musr/command construct box AlumPlateHole 2 5 5 G4_Galactic 0 0 0 log_AlumPlate norot dead 267
/musr/command construct box Insert 20 20 10 G4_Pb 0 0 -91 log_World norot dead 265
/musr/command construct box InsertHole 2.5 4 10 G4_Galactic 0 0 0 log_Insert norot dead 267
# --- SAMPLE HOLDER ---
##/musr/command construct GPDsampleHolderA GPDsampleHolderA 0 0 0 0 0 G4_Al 0 -111.25 0 log_World matrix1 dead 300
##/musr/command construct tubs GPDsampleHolderB 37.51 58 7 0 360 G4_Al 0 57 0 log_World matrix1 dead 301
##/musr/command construct tubs GPDsampleHolderC 37.51 58 7 0 360 G4_Al 0 -57 0 log_World matrix1 dead 302
######/musr/command construct tubs GPDsampleHolderD 37 37.5 50 136 88 G4_Al 0 0 0 log_World matrix1 dead 303
######/musr/command construct tubs GPDsampleHolderE 53.5 54 50 136 88 G4_Al 0 0 0 log_World matrix1 dead 304
##/musr/command construct box GPDsampleHolderD 12 0.5 50 G4_Al 38 0 28 log_World matrix1 dead 303
##/musr/command construct box GPDsampleHolderE 12 0.5 50 G4_Al 38 0 -28 log_World matrix1 dead 304
##/musr/command construct box GPDsampleHolderF 12 0.5 50 G4_Al -38 0 28 log_World matrix1 dead 305
##/musr/command construct box GPDsampleHolderG 12 0.5 50 G4_Al -38 0 -28 log_World matrix1 dead 306
##/musr/command construct tubs GPDsampleHolderH 0 9 50 180 50 G4_Al 48 0 29 log_World matrix1 dead 307
##/musr/command construct tubs GPDsampleHolderI 0 9 50 130 50 G4_Al 48 0 -29 log_World matrix1 dead 308
##/musr/command construct tubs GPDsampleHolderJ 0 9 50 310 50 G4_Al -48 0 29 log_World matrix1 dead 309
##/musr/command construct tubs GPDsampleHolderK 0 9 50 0 50 G4_Al -48 0 -29 log_World matrix1 dead 310
##/musr/command construct tubs GPDsampleHolderW 58.01 58.2 64 25 130 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 313
##/musr/command construct tubs GPDsampleHolderX 58.01 58.2 64 205 130 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 314
##/musr/command construct tubs GPDsampleHolderY 36.79 36.99 64 45 90 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 315
##/musr/command construct tubs GPDsampleHolderZ 36.79 36.99 64 225 90 G4_POLYETHYLENE 0 0 0 log_World matrix1 dead 316
# --- SAMPLE CELL ---
/musr/command construct tubs GPDsampleCell 0 12 50 0 360 G4_Cu 0 0 0 log_World matrix1 dead 211
# --- SAMPLE ---
/musr/command construct tubs target 0 3 5 0 360 G4_Cu 0 0 0 log_GPDsampleCell norot dead 201
#/musr/command construct tubs targetFieldVol 0 3 5 0 360 G4_Cu 0 0 0 log_GPDsampleCell norot dead 201
# --- SCINTILLATOR ---
/musr/command construct box ScintB1 21.5 50 2 G4_PLASTIC_SC_VINYLTOLUENE 26.5 0 -44 log_World norot musr/ScintSD 1
/musr/command construct box ScintB2 21.5 50 2 G4_PLASTIC_SC_VINYLTOLUENE -26.5 0 -44 log_World norot musr/ScintSD 2
/musr/command construct box ScintB4 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE 44 0 -25 log_World norot musr/ScintSD 4
/musr/command construct box ScintB5 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE -44 0 -25 log_World norot musr/ScintSD 5
/musr/command construct box ScintB6 2 50 13 G4_PLASTIC_SC_VINYLTOLUENE 44 0 0 log_World norot musr/ScintSD 6
/musr/command construct box ScintB7 2 50 13 G4_PLASTIC_SC_VINYLTOLUENE -44 0 0 log_World norot musr/ScintSD 7
/musr/command construct box ScintF1 20 50 2 G4_PLASTIC_SC_VINYLTOLUENE 28 0 44 log_World norot musr/ScintSD 11
/musr/command construct box ScintF2 20 50 2 G4_PLASTIC_SC_VINYLTOLUENE -28 0 44 log_World norot musr/ScintSD 12
/musr/command construct box ScintF3 7 50 2 G4_PLASTIC_SC_VINYLTOLUENE 0 0 44 log_World norot musr/ScintSD 13
/musr/command construct box ScintF4 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE 44 0 25 log_World norot musr/ScintSD 14
/musr/command construct box ScintF5 2 50 11 G4_PLASTIC_SC_VINYLTOLUENE -44 0 25 log_World norot musr/ScintSD 15
#/musr/command construct box KoincB1 20 90 2 G4_PLASTIC_SC_VINYLTOLUENE 33 0 -60 log_World matrix2 musr/ScintSD 21
#/musr/command construct box KoincB2 20 90 2 G4_PLASTIC_SC_VINYLTOLUENE -33 0 -60 log_World matrix3 musr/ScintSD 22
# --- M0 ---
/musr/command construct box M0holder 3 4 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -20 log_World norot musr/ScintSD 103
/musr/command construct box M0 4 30 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -44 log_World norot musr/ScintSD 102
#/musr/command construct GPDmHolder GPDmHolder 5 25 1.0 G4_PLASTIC_SC_VINYLTOLUENE 0 -25 -43 log_World norot dead 103
# --- Magnet ---
/musr/command construct tubs Magnet1 180 375 78 0 360 G4_Cu 133 0 0 log_World matrix5 dead 291
/musr/command construct tubs Magnet2 180 375 78 0 360 G4_Cu -133 0 0 log_World matrix5 dead 292
/musr/command construct uprofile Uprof1 25 10 400 5 G4_Al -45 0 200 log_World matrix6 dead 293
/musr/command construct uprofile Uprof2 25 10 400 5 G4_Al -45 0 -200 log_World matrix6 dead 293
/musr/command construct uprofile Uprof3 25 10 400 5 G4_Al 45 0 200 log_World matrix7 dead 293
/musr/command construct uprofile Uprof4 25 10 400 5 G4_Al 45 0 -200 log_World matrix7 dead 293
#
##---------------------------------------------------------
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_1 6
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_2 5
#/musr/command logicalVolumeToBeReweighted mu log_activeCollim 4
##---------------------------------------------------------
/musr/command visattributes log_World invisible
/musr/command visattributes log_beampipeA invisible
/musr/command visattributes log_beampipe invisible
/musr/command visattributes log_beampipeshell invisible
/musr/command visattributes log_beampipewindow invisible
/musr/command visattributes log_Magnet1 invisible
/musr/command visattributes log_Magnet2 invisible
/musr/command visattributes log_Uprof1 invisible
/musr/command visattributes log_Uprof2 invisible
/musr/command visattributes log_Uprof3 invisible
/musr/command visattributes log_Uprof4 invisible
/musr/command visattributes log_GPDsampleHolderD red
/musr/command visattributes log_GPDsampleHolderE red
/musr/command visattributes log_GPDsampleHolderF red
/musr/command visattributes log_GPDsampleHolderG red
/musr/command visattributes log_GPDcollimator green
/musr/command visattributes log_GPDcollimator2 blue
/musr/command visattributes log_GPDmHolder magenta
/musr/command visattributes log_target blue
/musr/command visattributes log_M0 yellow
#/musr/command visattributes G4_PLASTIC_SC_VINYLTOLUENE red
#/musr/command visattributes G4_Pb blue
####################################################################################
######################### M A G N E T I C F I E L D #########################
###################################################################################
# Set magnetic field (set field intensity in T and sigma in mm)
# syntax for magneticfield: fromfile filename fieldValue
# uniform fieldValue
# gaussian fieldValue sigma
#
/musr/command globalfield centralSolenoidField 0. 0. 0. fromfile 2D GPD_coil_primitive.table log_target 0.03
/musr/command globalfield setparameter SetLargestAcceptableStep 200
/musr/command globalfield setparameter SetMinimumEpsilonStep 5e-5
/musr/command globalfield setparameter SetMaximumEpsilonStep 0.001
/musr/command globalfield setparameter SetDeltaOneStep 0.1
/musr/command globalfield setparameter SetDeltaIntersection 0.01
/musr/command globalfield printparameters
/musr/command globalfield printFieldValueAtPoint 0 0 0
/musr/command globalfield printFieldValueAtPoint 0 0 10
/musr/command globalfield printFieldValueAtPoint 0 0 20
/musr/command globalfield printFieldValueAtPoint 0 0 30
/musr/command globalfield printFieldValueAtPoint 0 0 50
/musr/command globalfield printFieldValueAtPoint 0 0 100
/musr/command globalfield printFieldValueAtPoint 0 0 200
/musr/command globalfield printFieldValueAtPoint 0 0 300
/musr/command globalfield printFieldValueAtPoint 0 0 400
/musr/command globalfield printFieldValueAtPoint 0 0 500
/musr/command globalfield printFieldValueAtPoint 10 0 0
/musr/command globalfield printFieldValueAtPoint 100 0 0
/musr/command globalfield printFieldValueAtPoint 300 0 0
/musr/command globalfield printFieldValueAtPoint 0 10 0
/musr/command globalfield printFieldValueAtPoint 0 100 0
/musr/command globalfield printFieldValueAtPoint 0 300 0
###################################################################################
######################### P H Y S I C S P R O C E S S E S ##################
###################################################################################
# Geant 4.9.3
/musr/command process addDiscreteProcess gamma G4PhotoElectricEffect
/musr/command process addDiscreteProcess gamma G4ComptonScattering
/musr/command process addDiscreteProcess gamma G4GammaConversion
/musr/command process addDiscreteProcess gamma G4RayleighScattering
/musr/command process addProcess e- G4eMultipleScattering -1 1 1
/musr/command process addProcess e- G4eIonisation -1 2 2
/musr/command process addProcess e- G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eMultipleScattering -1 1 1
/musr/command process addProcess e+ G4eIonisation -1 2 2
/musr/command process addProcess e+ G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eplusAnnihilation 0 -1 4
/musr/command process addProcess mu- G4MuMultipleScattering -1 1 1
/musr/command process addProcess mu- G4MuIonisation -1 2 2
/musr/command process addProcess mu- G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu- G4MuPairProduction -1 4 4
/musr/command process addProcess mu+ G4MuMultipleScattering -1 1 1
/musr/command process addProcess mu+ G4MuIonisation -1 2 2
/musr/command process addProcess mu+ G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu+ G4MuPairProduction -1 4 4
#
#/home/install/geant4.9.3/source/processes/electromagnetic/utils/src/G4EnergyLossMessenger.cc
######## /process/msc/StepLimit Minimal | UseDistanceToBoundary | UseSafety
/process/msc/StepLimit UseSafety
#/process/msc/LateralDisplacement
/process/msc/RangeFactor 0.04
/process/msc/GeomFactor 2.5
#/process/msc/FactorForAngleLimit
/process/msc/Skin 3.0
#/process/msc/ThetaLimit 0.2 rad
###################################################################################
################## S O M E O T H E R P A R A M E T E R S ##################
###################################################################################
# Set the overall range cut
#/run/setCut 1 mm
# Set the range cut on a particular volumes (in mm)
#/musr/command SetUserLimits log_target 0.01
#/musr/command SetUserLimits log_targetscint 0.01
#/musr/command SetUserLimits log_M0 0.01
# Store all events into the ROOT tree or just the interesting ones ? (true is default)
/musr/command storeOnlyEventsWithHits false
# Set the minimum time separation between two subsequent signals in the same detector (in ns)
/musr/command signalSeparationTime 10
#
#/musr/run/howOftenToPrintEvent 10
/musr/command maximumRunTimeAllowed 220000
/musr/run/randomOption 2
###################################################################################
######################### R O O T O U T P U T ##############################
###################################################################################
/musr/command rootOutput det_VvvKine off
/musr/command rootOutput det_VvvX off
/musr/command rootOutput det_VvvY off
/musr/command rootOutput det_VvvZ off
/musr/command rootOutput det_VvvVolID off
/musr/command rootOutput det_VvvProcID off
/musr/command rootOutput det_VvvTrackID off
/musr/command rootOutput det_VvvParticleID off
/musr/command rootOutput nOptPhot off
/musr/command rootOutput odet_ID off
/musr/command rootOutput odet_nPhot off
/musr/command rootOutput odet_timeFirst off
/musr/command rootOutput odet_timeA off
/musr/command rootOutput odet_timeB off
/musr/command rootOutput odet_timeC off
/musr/command rootOutput odet_timeD off
/musr/command rootOutput odet_timeE off
/musr/command rootOutput odet_timeLast off
###################################################################################
######################### V I S U A L I S A T I O N ##############################
###################################################################################
/vis/disable
#/control/execute visFromToni.mac
#/control/execute visDawn101.mac
#/control/execute visDawn40123.mac
#/control/execute visVRML.mac
###################################################################################
######################### P A R T I C L E G U N #################################
###################################################################################
/gun/vertex 0 0 -1000 mm
# FWHM 10mm ==> sigma = 10/2.354 = 4.2481mm
/gun/vertexsigma 25 25 0 mm
#---/gun/vertexboundary: rMaxAllowed, zMinAllowed, zMaxAllowed
/gun/vertexboundary 100 -999999 999999 mm
/gun/momentum 100.0 MeV
# sigma = 3% ==> sigma 27*0.03 = 0.81
/gun/momentumsmearing 3.00 MeV
#---/gun/momentumboundary: pMinAllowed, pMaxAllowed, dummy
#/gun/momentumboundary 20 40 0 MeV
# TURTLE
#/gun/turtlefilename FOR061_2008_04_22.DAT
#/gun/turtlefilename FOR070_2008_10_17_XXII.DAT
#/gun/turtlefilename FOR077_pie3_HiField_d05_x30.dat
#/gun/turtlefilename FOR077_reggiani_Jan2010_NEW.dat
#/gun/turtleZ0position -900 mm
#/gun/turtlefilename FOR070_GPD_beam_profile_45mm.dat
#/gun/turtleZ0position -500 mm
#---/gun/tilt: xangle, yangle, dummy
#/gun/tilt 0 0.5 0 deg
#---/gun/tiltsigma: xangleSigma, yangleSigma, dummy (1 degree on 1 meter ~ 17mm)
/gun/tiltsigma 0.15 0.15 0 deg
/gun/pitch 0.08 deg
/gun/muonPolarizVector 0 0 -1
#---/gun/decaytimelimits: decayMin, decayMax, decayTime
/gun/decaytimelimits -1 -1 2197.03 ns
###################################################################################
######################## B E A M O N #########################################
###################################################################################
/run/beamOn 16000000
#/run/beamOn 1

View File

@ -0,0 +1,611 @@
301 1 2 1
1 X [LENGU]
2 Y [LENGU]
3 Z [LENGU]
4 HX [FIELU]
5 HZ [FIELU]
6 HMOD [FIELU]
fieldNormalisation 880.169809
0
0.0000000e+00 0 0 0 1.1361444e-03 0
0.0000000e+00 0 20 0 1.1361444e-03 0
1.0000000e+00 0 0 0 1.1362377e-03 0
1.0000000e+00 0 20 0 1.1362377e-03 0
2.0000000e+00 0 0 0 1.1365074e-03 0
2.0000000e+00 0 20 0 1.1365074e-03 0
3.0000000e+00 0 0 0 1.1369226e-03 0
3.0000000e+00 0 20 0 1.1369226e-03 0
4.0000000e+00 0 0 0 1.1374296e-03 0
4.0000000e+00 0 20 0 1.1374296e-03 0
5.0000000e+00 0 0 0 1.1379495e-03 0
5.0000000e+00 0 20 0 1.1379495e-03 0
6.0000000e+00 0 0 0 1.1383740e-03 0
6.0000000e+00 0 20 0 1.1383740e-03 0
7.0000000e+00 0 0 0 1.1385589e-03 0
7.0000000e+00 0 20 0 1.1385589e-03 0
8.0000000e+00 0 0 0 1.1383168e-03 0
8.0000000e+00 0 20 0 1.1383168e-03 0
9.0000000e+00 0 0 0 1.1374066e-03 0
9.0000000e+00 0 20 0 1.1374066e-03 0
1.0000000e+01 0 0 0 1.1355196e-03 0
1.0000000e+01 0 20 0 1.1355196e-03 0
1.1000000e+01 0 0 0 1.1322640e-03 0
1.1000000e+01 0 20 0 1.1322640e-03 0
1.2000000e+01 0 0 0 1.1271451e-03 0
1.2000000e+01 0 20 0 1.1271451e-03 0
1.3000000e+01 0 0 0 1.1195471e-03 0
1.3000000e+01 0 20 0 1.1195471e-03 0
1.4000000e+01 0 0 0 1.1087207e-03 0
1.4000000e+01 0 20 0 1.1087207e-03 0
1.5000000e+01 0 0 0 1.0937932e-03 0
1.5000000e+01 0 20 0 1.0937932e-03 0
1.6000000e+01 0 0 0 1.0738238e-03 0
1.6000000e+01 0 20 0 1.0738238e-03 0
1.7000000e+01 0 0 0 1.0479326e-03 0
1.7000000e+01 0 20 0 1.0479326e-03 0
1.8000000e+01 0 0 0 1.0155087e-03 0
1.8000000e+01 0 20 0 1.0155087e-03 0
1.9000000e+01 0 0 0 9.7642128e-04 0
1.9000000e+01 0 20 0 9.7642128e-04 0
2.0000000e+01 0 0 0 9.3109490e-04 0
2.0000000e+01 0 20 0 9.3109490e-04 0
2.1000000e+01 0 0 0 8.8037814e-04 0
2.1000000e+01 0 20 0 8.8037814e-04 0
2.2000000e+01 0 0 0 8.2530227e-04 0
2.2000000e+01 0 20 0 8.2530227e-04 0
2.3000000e+01 0 0 0 7.6687852e-04 0
2.3000000e+01 0 20 0 7.6687852e-04 0
2.4000000e+01 0 0 0 7.0599560e-04 0
2.4000000e+01 0 20 0 7.0599560e-04 0
2.5000000e+01 0 0 0 6.4339740e-04 0
2.5000000e+01 0 20 0 6.4339740e-04 0
2.6000000e+01 0 0 0 5.7970163e-04 0
2.6000000e+01 0 20 0 5.7970163e-04 0
2.7000000e+01 0 0 0 5.1543274e-04 0
2.7000000e+01 0 20 0 5.1543274e-04 0
2.8000000e+01 0 0 0 4.5105685e-04 0
2.8000000e+01 0 20 0 4.5105685e-04 0
2.9000000e+01 0 0 0 3.8701483e-04 0
2.9000000e+01 0 20 0 3.8701483e-04 0
3.0000000e+01 0 0 0 3.2375422e-04 0
3.0000000e+01 0 20 0 3.2375422e-04 0
3.1000000e+01 0 0 0 2.6176224e-04 0
3.1000000e+01 0 20 0 2.6176224e-04 0
3.2000000e+01 0 0 0 2.0160348e-04 0
3.2000000e+01 0 20 0 2.0160348e-04 0
3.3000000e+01 0 0 0 1.4396412e-04 0
3.3000000e+01 0 20 0 1.4396412e-04 0
3.4000000e+01 0 0 0 8.9697678e-05 0
3.4000000e+01 0 20 0 8.9697678e-05 0
3.5000000e+01 0 0 0 3.9848332e-05 0
3.5000000e+01 0 20 0 3.9848332e-05 0
3.6000000e+01 0 0 0 -4.4041645e-06 0
3.6000000e+01 0 20 0 -4.4041645e-06 0
3.7000000e+01 0 0 0 -4.1948665e-05 0
3.7000000e+01 0 20 0 -4.1948665e-05 0
3.8000000e+01 0 0 0 -7.2089079e-05 0
3.8000000e+01 0 20 0 -7.2089079e-05 0
3.9000000e+01 0 0 0 -9.4844419e-05 0
3.9000000e+01 0 20 0 -9.4844419e-05 0
4.0000000e+01 0 0 0 -1.1095290e-04 0
4.0000000e+01 0 20 0 -1.1095290e-04 0
4.1000000e+01 0 0 0 -1.2157438e-04 0
4.1000000e+01 0 20 0 -1.2157438e-04 0
4.2000000e+01 0 0 0 -1.2793545e-04 0
4.2000000e+01 0 20 0 -1.2793545e-04 0
4.3000000e+01 0 0 0 -1.3111042e-04 0
4.3000000e+01 0 20 0 -1.3111042e-04 0
4.4000000e+01 0 0 0 -1.3195272e-04 0
4.4000000e+01 0 20 0 -1.3195272e-04 0
4.5000000e+01 0 0 0 -1.3110933e-04 0
4.5000000e+01 0 20 0 -1.3110933e-04 0
4.6000000e+01 0 0 0 -1.2906196e-04 0
4.6000000e+01 0 20 0 -1.2906196e-04 0
4.7000000e+01 0 0 0 -1.2616788e-04 0
4.7000000e+01 0 20 0 -1.2616788e-04 0
4.8000000e+01 0 0 0 -1.2269274e-04 0
4.8000000e+01 0 20 0 -1.2269274e-04 0
4.9000000e+01 0 0 0 -1.1883471e-04 0
4.9000000e+01 0 20 0 -1.1883471e-04 0
5.0000000e+01 0 0 0 -1.1474199e-04 0
5.0000000e+01 0 20 0 -1.1474199e-04 0
5.1000000e+01 0 0 0 -1.1052536e-04 0
5.1000000e+01 0 20 0 -1.1052536e-04 0
5.2000000e+01 0 0 0 -1.0626737e-04 0
5.2000000e+01 0 20 0 -1.0626737e-04 0
5.3000000e+01 0 0 0 -1.0202909e-04 0
5.3000000e+01 0 20 0 -1.0202909e-04 0
5.4000000e+01 0 0 0 -9.7855223e-05 0
5.4000000e+01 0 20 0 -9.7855223e-05 0
5.5000000e+01 0 0 0 -9.3777888e-05 0
5.5000000e+01 0 20 0 -9.3777888e-05 0
5.6000000e+01 0 0 0 -8.9819586e-05 0
5.6000000e+01 0 20 0 -8.9819586e-05 0
5.7000000e+01 0 0 0 -8.5995449e-05 0
5.7000000e+01 0 20 0 -8.5995449e-05 0
5.8000000e+01 0 0 0 -8.2314978e-05 0
5.8000000e+01 0 20 0 -8.2314978e-05 0
5.9000000e+01 0 0 0 -7.8783399e-05 0
5.9000000e+01 0 20 0 -7.8783399e-05 0
6.0000000e+01 0 0 0 -7.5402713e-05 0
6.0000000e+01 0 20 0 -7.5402713e-05 0
6.1000000e+01 0 0 0 -7.2172507e-05 0
6.1000000e+01 0 20 0 -7.2172507e-05 0
6.2000000e+01 0 0 0 -6.9090590e-05 0
6.2000000e+01 0 20 0 -6.9090590e-05 0
6.3000000e+01 0 0 0 -6.6153484e-05 0
6.3000000e+01 0 20 0 -6.6153484e-05 0
6.4000000e+01 0 0 0 -6.3356805e-05 0
6.4000000e+01 0 20 0 -6.3356805e-05 0
6.5000000e+01 0 0 0 -6.0695551e-05 0
6.5000000e+01 0 20 0 -6.0695551e-05 0
6.6000000e+01 0 0 0 -5.8164335e-05 0
6.6000000e+01 0 20 0 -5.8164335e-05 0
6.7000000e+01 0 0 0 -5.5757548e-05 0
6.7000000e+01 0 20 0 -5.5757548e-05 0
6.8000000e+01 0 0 0 -5.3469496e-05 0
6.8000000e+01 0 20 0 -5.3469496e-05 0
6.9000000e+01 0 0 0 -5.1294495e-05 0
6.9000000e+01 0 20 0 -5.1294495e-05 0
7.0000000e+01 0 0 0 -4.9226940e-05 0
7.0000000e+01 0 20 0 -4.9226940e-05 0
7.1000000e+01 0 0 0 -4.7261364e-05 0
7.1000000e+01 0 20 0 -4.7261364e-05 0
7.2000000e+01 0 0 0 -4.5392469e-05 0
7.2000000e+01 0 20 0 -4.5392469e-05 0
7.3000000e+01 0 0 0 -4.3615153e-05 0
7.3000000e+01 0 20 0 -4.3615153e-05 0
7.4000000e+01 0 0 0 -4.1924529e-05 0
7.4000000e+01 0 20 0 -4.1924529e-05 0
7.5000000e+01 0 0 0 -4.0315927e-05 0
7.5000000e+01 0 20 0 -4.0315927e-05 0
7.6000000e+01 0 0 0 -3.8784905e-05 0
7.6000000e+01 0 20 0 -3.8784905e-05 0
7.7000000e+01 0 0 0 -3.7327244e-05 0
7.7000000e+01 0 20 0 -3.7327244e-05 0
7.8000000e+01 0 0 0 -3.5938946e-05 0
7.8000000e+01 0 20 0 -3.5938946e-05 0
7.9000000e+01 0 0 0 -3.4616230e-05 0
7.9000000e+01 0 20 0 -3.4616230e-05 0
8.0000000e+01 0 0 0 -3.3355521e-05 0
8.0000000e+01 0 20 0 -3.3355521e-05 0
8.1000000e+01 0 0 0 -3.2153447e-05 0
8.1000000e+01 0 20 0 -3.2153447e-05 0
8.2000000e+01 0 0 0 -3.1006827e-05 0
8.2000000e+01 0 20 0 -3.1006827e-05 0
8.3000000e+01 0 0 0 -2.9912663e-05 0
8.3000000e+01 0 20 0 -2.9912663e-05 0
8.4000000e+01 0 0 0 -2.8868131e-05 0
8.4000000e+01 0 20 0 -2.8868131e-05 0
8.5000000e+01 0 0 0 -2.7870571e-05 0
8.5000000e+01 0 20 0 -2.7870571e-05 0
8.6000000e+01 0 0 0 -2.6917478e-05 0
8.6000000e+01 0 20 0 -2.6917478e-05 0
8.7000000e+01 0 0 0 -2.6006495e-05 0
8.7000000e+01 0 20 0 -2.6006495e-05 0
8.8000000e+01 0 0 0 -2.5135405e-05 0
8.8000000e+01 0 20 0 -2.5135405e-05 0
8.9000000e+01 0 0 0 -2.4302117e-05 0
8.9000000e+01 0 20 0 -2.4302117e-05 0
9.0000000e+01 0 0 0 -2.3504665e-05 0
9.0000000e+01 0 20 0 -2.3504665e-05 0
9.1000000e+01 0 0 0 -2.2741199e-05 0
9.1000000e+01 0 20 0 -2.2741199e-05 0
9.2000000e+01 0 0 0 -2.2009975e-05 0
9.2000000e+01 0 20 0 -2.2009975e-05 0
9.3000000e+01 0 0 0 -2.1309350e-05 0
9.3000000e+01 0 20 0 -2.1309350e-05 0
9.4000000e+01 0 0 0 -2.0637779e-05 0
9.4000000e+01 0 20 0 -2.0637779e-05 0
9.5000000e+01 0 0 0 -1.9993803e-05 0
9.5000000e+01 0 20 0 -1.9993803e-05 0
9.6000000e+01 0 0 0 -1.9376048e-05 0
9.6000000e+01 0 20 0 -1.9376048e-05 0
9.7000000e+01 0 0 0 -1.8783220e-05 0
9.7000000e+01 0 20 0 -1.8783220e-05 0
9.8000000e+01 0 0 0 -1.8214096e-05 0
9.8000000e+01 0 20 0 -1.8214096e-05 0
9.9000000e+01 0 0 0 -1.7667524e-05 0
9.9000000e+01 0 20 0 -1.7667524e-05 0
1.0000000e+02 0 0 0 -1.7142415e-05 0
1.0000000e+02 0 20 0 -1.7142415e-05 0
1.0100000e+02 0 0 0 -1.6637742e-05 0
1.0100000e+02 0 20 0 -1.6637742e-05 0
1.0200000e+02 0 0 0 -1.6152534e-05 0
1.0200000e+02 0 20 0 -1.6152534e-05 0
1.0300000e+02 0 0 0 -1.5685876e-05 0
1.0300000e+02 0 20 0 -1.5685876e-05 0
1.0400000e+02 0 0 0 -1.5236899e-05 0
1.0400000e+02 0 20 0 -1.5236899e-05 0
1.0500000e+02 0 0 0 -1.4804785e-05 0
1.0500000e+02 0 20 0 -1.4804785e-05 0
1.0600000e+02 0 0 0 -1.4388759e-05 0
1.0600000e+02 0 20 0 -1.4388759e-05 0
1.0700000e+02 0 0 0 -1.3988086e-05 0
1.0700000e+02 0 20 0 -1.3988086e-05 0
1.0800000e+02 0 0 0 -1.3602072e-05 0
1.0800000e+02 0 20 0 -1.3602072e-05 0
1.0900000e+02 0 0 0 -1.3230060e-05 0
1.0900000e+02 0 20 0 -1.3230060e-05 0
1.1000000e+02 0 0 0 -1.2871425e-05 0
1.1000000e+02 0 20 0 -1.2871425e-05 0
1.1100000e+02 0 0 0 -1.2525578e-05 0
1.1100000e+02 0 20 0 -1.2525578e-05 0
1.1200000e+02 0 0 0 -1.2191958e-05 0
1.1200000e+02 0 20 0 -1.2191958e-05 0
1.1300000e+02 0 0 0 -1.1870033e-05 0
1.1300000e+02 0 20 0 -1.1870033e-05 0
1.1400000e+02 0 0 0 -1.1559300e-05 0
1.1400000e+02 0 20 0 -1.1559300e-05 0
1.1500000e+02 0 0 0 -1.1259280e-05 0
1.1500000e+02 0 20 0 -1.1259280e-05 0
1.1600000e+02 0 0 0 -1.0969518e-05 0
1.1600000e+02 0 20 0 -1.0969518e-05 0
1.1700000e+02 0 0 0 -1.0689581e-05 0
1.1700000e+02 0 20 0 -1.0689581e-05 0
1.1800000e+02 0 0 0 -1.0419061e-05 0
1.1800000e+02 0 20 0 -1.0419061e-05 0
1.1900000e+02 0 0 0 -1.0157566e-05 0
1.1900000e+02 0 20 0 -1.0157566e-05 0
1.2000000e+02 0 0 0 -9.9047248e-06 0
1.2000000e+02 0 20 0 -9.9047248e-06 0
1.2100000e+02 0 0 0 -9.6601849e-06 0
1.2100000e+02 0 20 0 -9.6601849e-06 0
1.2200000e+02 0 0 0 -9.4236097e-06 0
1.2200000e+02 0 20 0 -9.4236097e-06 0
1.2300000e+02 0 0 0 -9.1946792e-06 0
1.2300000e+02 0 20 0 -9.1946792e-06 0
1.2400000e+02 0 0 0 -8.9730884e-06 0
1.2400000e+02 0 20 0 -8.9730884e-06 0
1.2500000e+02 0 0 0 -8.7585467e-06 0
1.2500000e+02 0 20 0 -8.7585467e-06 0
1.2600000e+02 0 0 0 -8.5507773e-06 0
1.2600000e+02 0 20 0 -8.5507773e-06 0
1.2700000e+02 0 0 0 -8.3495160e-06 0
1.2700000e+02 0 20 0 -8.3495160e-06 0
1.2800000e+02 0 0 0 -8.1545110e-06 0
1.2800000e+02 0 20 0 -8.1545110e-06 0
1.2900000e+02 0 0 0 -7.9655218e-06 0
1.2900000e+02 0 20 0 -7.9655218e-06 0
1.3000000e+02 0 0 0 -7.7823191e-06 0
1.3000000e+02 0 20 0 -7.7823191e-06 0
1.3100000e+02 0 0 0 -7.6046838e-06 0
1.3100000e+02 0 20 0 -7.6046838e-06 0
1.3200000e+02 0 0 0 -7.4324065e-06 0
1.3200000e+02 0 20 0 -7.4324065e-06 0
1.3300000e+02 0 0 0 -7.2652875e-06 0
1.3300000e+02 0 20 0 -7.2652875e-06 0
1.3400000e+02 0 0 0 -7.1031356e-06 0
1.3400000e+02 0 20 0 -7.1031356e-06 0
1.3500000e+02 0 0 0 -6.9457681e-06 0
1.3500000e+02 0 20 0 -6.9457681e-06 0
1.3600000e+02 0 0 0 -6.7930102e-06 0
1.3600000e+02 0 20 0 -6.7930102e-06 0
1.3700000e+02 0 0 0 -6.6446950e-06 0
1.3700000e+02 0 20 0 -6.6446950e-06 0
1.3800000e+02 0 0 0 -6.5006623e-06 0
1.3800000e+02 0 20 0 -6.5006623e-06 0
1.3900000e+02 0 0 0 -6.3607590e-06 0
1.3900000e+02 0 20 0 -6.3607590e-06 0
1.4000000e+02 0 0 0 -6.2248386e-06 0
1.4000000e+02 0 20 0 -6.2248386e-06 0
1.4100000e+02 0 0 0 -6.0927605e-06 0
1.4100000e+02 0 20 0 -6.0927605e-06 0
1.4200000e+02 0 0 0 -5.9643902e-06 0
1.4200000e+02 0 20 0 -5.9643902e-06 0
1.4300000e+02 0 0 0 -5.8395988e-06 0
1.4300000e+02 0 20 0 -5.8395988e-06 0
1.4400000e+02 0 0 0 -5.7182625e-06 0
1.4400000e+02 0 20 0 -5.7182625e-06 0
1.4500000e+02 0 0 0 -5.6002629e-06 0
1.4500000e+02 0 20 0 -5.6002629e-06 0
1.4600000e+02 0 0 0 -5.4854861e-06 0
1.4600000e+02 0 20 0 -5.4854861e-06 0
1.4700000e+02 0 0 0 -5.3738232e-06 0
1.4700000e+02 0 20 0 -5.3738232e-06 0
1.4800000e+02 0 0 0 -5.2651692e-06 0
1.4800000e+02 0 20 0 -5.2651692e-06 0
1.4900000e+02 0 0 0 -5.1594237e-06 0
1.4900000e+02 0 20 0 -5.1594237e-06 0
1.5000000e+02 0 0 0 -5.0564902e-06 0
1.5000000e+02 0 20 0 -5.0564902e-06 0
1.5100000e+02 0 0 0 -4.9562759e-06 0
1.5100000e+02 0 20 0 -4.9562759e-06 0
1.5200000e+02 0 0 0 -4.8586917e-06 0
1.5200000e+02 0 20 0 -4.8586917e-06 0
1.5300000e+02 0 0 0 -4.7636519e-06 0
1.5300000e+02 0 20 0 -4.7636519e-06 0
1.5400000e+02 0 0 0 -4.6710744e-06 0
1.5400000e+02 0 20 0 -4.6710744e-06 0
1.5500000e+02 0 0 0 -4.5808798e-06 0
1.5500000e+02 0 20 0 -4.5808798e-06 0
1.5600000e+02 0 0 0 -4.4929922e-06 0
1.5600000e+02 0 20 0 -4.4929922e-06 0
1.5700000e+02 0 0 0 -4.4073382e-06 0
1.5700000e+02 0 20 0 -4.4073382e-06 0
1.5800000e+02 0 0 0 -4.3238475e-06 0
1.5800000e+02 0 20 0 -4.3238475e-06 0
1.5900000e+02 0 0 0 -4.2424522e-06 0
1.5900000e+02 0 20 0 -4.2424522e-06 0
1.6000000e+02 0 0 0 -4.1630870e-06 0
1.6000000e+02 0 20 0 -4.1630870e-06 0
1.6100000e+02 0 0 0 -4.0856892e-06 0
1.6100000e+02 0 20 0 -4.0856892e-06 0
1.6200000e+02 0 0 0 -4.0101980e-06 0
1.6200000e+02 0 20 0 -4.0101980e-06 0
1.6300000e+02 0 0 0 -3.9365553e-06 0
1.6300000e+02 0 20 0 -3.9365553e-06 0
1.6400000e+02 0 0 0 -3.8647049e-06 0
1.6400000e+02 0 20 0 -3.8647049e-06 0
1.6500000e+02 0 0 0 -3.7945925e-06 0
1.6500000e+02 0 20 0 -3.7945925e-06 0
1.6600000e+02 0 0 0 -3.7261659e-06 0
1.6600000e+02 0 20 0 -3.7261659e-06 0
1.6700000e+02 0 0 0 -3.6593749e-06 0
1.6700000e+02 0 20 0 -3.6593749e-06 0
1.6800000e+02 0 0 0 -3.5941708e-06 0
1.6800000e+02 0 20 0 -3.5941708e-06 0
1.6900000e+02 0 0 0 -3.5305069e-06 0
1.6900000e+02 0 20 0 -3.5305069e-06 0
1.7000000e+02 0 0 0 -3.4683379e-06 0
1.7000000e+02 0 20 0 -3.4683379e-06 0
1.7100000e+02 0 0 0 -3.4076203e-06 0
1.7100000e+02 0 20 0 -3.4076203e-06 0
1.7200000e+02 0 0 0 -3.3483118e-06 0
1.7200000e+02 0 20 0 -3.3483118e-06 0
1.7300000e+02 0 0 0 -3.2903720e-06 0
1.7300000e+02 0 20 0 -3.2903720e-06 0
1.7400000e+02 0 0 0 -3.2337616e-06 0
1.7400000e+02 0 20 0 -3.2337616e-06 0
1.7500000e+02 0 0 0 -3.1784426e-06 0
1.7500000e+02 0 20 0 -3.1784426e-06 0
1.7600000e+02 0 0 0 -3.1243785e-06 0
1.7600000e+02 0 20 0 -3.1243785e-06 0
1.7700000e+02 0 0 0 -3.0715338e-06 0
1.7700000e+02 0 20 0 -3.0715338e-06 0
1.7800000e+02 0 0 0 -3.0198745e-06 0
1.7800000e+02 0 20 0 -3.0198745e-06 0
1.7900000e+02 0 0 0 -2.9693674e-06 0
1.7900000e+02 0 20 0 -2.9693674e-06 0
1.8000000e+02 0 0 0 -2.9199807e-06 0
1.8000000e+02 0 20 0 -2.9199807e-06 0
1.8100000e+02 0 0 0 -2.8716834e-06 0
1.8100000e+02 0 20 0 -2.8716834e-06 0
1.8200000e+02 0 0 0 -2.8244456e-06 0
1.8200000e+02 0 20 0 -2.8244456e-06 0
1.8300000e+02 0 0 0 -2.7782385e-06 0
1.8300000e+02 0 20 0 -2.7782385e-06 0
1.8400000e+02 0 0 0 -2.7330342e-06 0
1.8400000e+02 0 20 0 -2.7330342e-06 0
1.8500000e+02 0 0 0 -2.6888054e-06 0
1.8500000e+02 0 20 0 -2.6888054e-06 0
1.8600000e+02 0 0 0 -2.6455262e-06 0
1.8600000e+02 0 20 0 -2.6455262e-06 0
1.8700000e+02 0 0 0 -2.6031711e-06 0
1.8700000e+02 0 20 0 -2.6031711e-06 0
1.8800000e+02 0 0 0 -2.5617156e-06 0
1.8800000e+02 0 20 0 -2.5617156e-06 0
1.8900000e+02 0 0 0 -2.5211360e-06 0
1.8900000e+02 0 20 0 -2.5211360e-06 0
1.9000000e+02 0 0 0 -2.4814092e-06 0
1.9000000e+02 0 20 0 -2.4814092e-06 0
1.9100000e+02 0 0 0 -2.4425130e-06 0
1.9100000e+02 0 20 0 -2.4425130e-06 0
1.9200000e+02 0 0 0 -2.4044257e-06 0
1.9200000e+02 0 20 0 -2.4044257e-06 0
1.9300000e+02 0 0 0 -2.3671265e-06 0
1.9300000e+02 0 20 0 -2.3671265e-06 0
1.9400000e+02 0 0 0 -2.3305950e-06 0
1.9400000e+02 0 20 0 -2.3305950e-06 0
1.9500000e+02 0 0 0 -2.2948117e-06 0
1.9500000e+02 0 20 0 -2.2948117e-06 0
1.9600000e+02 0 0 0 -2.2597574e-06 0
1.9600000e+02 0 20 0 -2.2597574e-06 0
1.9700000e+02 0 0 0 -2.2254136e-06 0
1.9700000e+02 0 20 0 -2.2254136e-06 0
1.9800000e+02 0 0 0 -2.1917626e-06 0
1.9800000e+02 0 20 0 -2.1917626e-06 0
1.9900000e+02 0 0 0 -2.1587868e-06 0
1.9900000e+02 0 20 0 -2.1587868e-06 0
2.0000000e+02 0 0 0 -2.1264695e-06 0
2.0000000e+02 0 20 0 -2.1264695e-06 0
2.0100000e+02 0 0 0 -2.0947942e-06 0
2.0100000e+02 0 20 0 -2.0947942e-06 0
2.0200000e+02 0 0 0 -2.0637452e-06 0
2.0200000e+02 0 20 0 -2.0637452e-06 0
2.0300000e+02 0 0 0 -2.0333069e-06 0
2.0300000e+02 0 20 0 -2.0333069e-06 0
2.0400000e+02 0 0 0 -2.0034646e-06 0
2.0400000e+02 0 20 0 -2.0034646e-06 0
2.0500000e+02 0 0 0 -1.9742035e-06 0
2.0500000e+02 0 20 0 -1.9742035e-06 0
2.0600000e+02 0 0 0 -1.9455098e-06 0
2.0600000e+02 0 20 0 -1.9455098e-06 0
2.0700000e+02 0 0 0 -1.9173696e-06 0
2.0700000e+02 0 20 0 -1.9173696e-06 0
2.0800000e+02 0 0 0 -1.8897697e-06 0
2.0800000e+02 0 20 0 -1.8897697e-06 0
2.0900000e+02 0 0 0 -1.8626972e-06 0
2.0900000e+02 0 20 0 -1.8626972e-06 0
2.1000000e+02 0 0 0 -1.8361395e-06 0
2.1000000e+02 0 20 0 -1.8361395e-06 0
2.1100000e+02 0 0 0 -1.8100845e-06 0
2.1100000e+02 0 20 0 -1.8100845e-06 0
2.1200000e+02 0 0 0 -1.7845203e-06 0
2.1200000e+02 0 20 0 -1.7845203e-06 0
2.1300000e+02 0 0 0 -1.7594355e-06 0
2.1300000e+02 0 20 0 -1.7594355e-06 0
2.1400000e+02 0 0 0 -1.7348188e-06 0
2.1400000e+02 0 20 0 -1.7348188e-06 0
2.1500000e+02 0 0 0 -1.7106593e-06 0
2.1500000e+02 0 20 0 -1.7106593e-06 0
2.1600000e+02 0 0 0 -1.6869465e-06 0
2.1600000e+02 0 20 0 -1.6869465e-06 0
2.1700000e+02 0 0 0 -1.6636702e-06 0
2.1700000e+02 0 20 0 -1.6636702e-06 0
2.1800000e+02 0 0 0 -1.6408202e-06 0
2.1800000e+02 0 20 0 -1.6408202e-06 0
2.1900000e+02 0 0 0 -1.6183870e-06 0
2.1900000e+02 0 20 0 -1.6183870e-06 0
2.2000000e+02 0 0 0 -1.5963610e-06 0
2.2000000e+02 0 20 0 -1.5963610e-06 0
2.2100000e+02 0 0 0 -1.5747331e-06 0
2.2100000e+02 0 20 0 -1.5747331e-06 0
2.2200000e+02 0 0 0 -1.5534942e-06 0
2.2200000e+02 0 20 0 -1.5534942e-06 0
2.2300000e+02 0 0 0 -1.5326357e-06 0
2.2300000e+02 0 20 0 -1.5326357e-06 0
2.2400000e+02 0 0 0 -1.5121491e-06 0
2.2400000e+02 0 20 0 -1.5121491e-06 0
2.2500000e+02 0 0 0 -1.4920261e-06 0
2.2500000e+02 0 20 0 -1.4920261e-06 0
2.2600000e+02 0 0 0 -1.4722587e-06 0
2.2600000e+02 0 20 0 -1.4722587e-06 0
2.2700000e+02 0 0 0 -1.4528392e-06 0
2.2700000e+02 0 20 0 -1.4528392e-06 0
2.2800000e+02 0 0 0 -1.4337598e-06 0
2.2800000e+02 0 20 0 -1.4337598e-06 0
2.2900000e+02 0 0 0 -1.4150131e-06 0
2.2900000e+02 0 20 0 -1.4150131e-06 0
2.3000000e+02 0 0 0 -1.3965920e-06 0
2.3000000e+02 0 20 0 -1.3965920e-06 0
2.3100000e+02 0 0 0 -1.3784893e-06 0
2.3100000e+02 0 20 0 -1.3784893e-06 0
2.3200000e+02 0 0 0 -1.3606983e-06 0
2.3200000e+02 0 20 0 -1.3606983e-06 0
2.3300000e+02 0 0 0 -1.3432122e-06 0
2.3300000e+02 0 20 0 -1.3432122e-06 0
2.3400000e+02 0 0 0 -1.3260246e-06 0
2.3400000e+02 0 20 0 -1.3260246e-06 0
2.3500000e+02 0 0 0 -1.3091291e-06 0
2.3500000e+02 0 20 0 -1.3091291e-06 0
2.3600000e+02 0 0 0 -1.2925195e-06 0
2.3600000e+02 0 20 0 -1.2925195e-06 0
2.3700000e+02 0 0 0 -1.2761899e-06 0
2.3700000e+02 0 20 0 -1.2761899e-06 0
2.3800000e+02 0 0 0 -1.2601342e-06 0
2.3800000e+02 0 20 0 -1.2601342e-06 0
2.3900000e+02 0 0 0 -1.2443469e-06 0
2.3900000e+02 0 20 0 -1.2443469e-06 0
2.4000000e+02 0 0 0 -1.2288222e-06 0
2.4000000e+02 0 20 0 -1.2288222e-06 0
2.4100000e+02 0 0 0 -1.2135549e-06 0
2.4100000e+02 0 20 0 -1.2135549e-06 0
2.4200000e+02 0 0 0 -1.1985395e-06 0
2.4200000e+02 0 20 0 -1.1985395e-06 0
2.4300000e+02 0 0 0 -1.1837709e-06 0
2.4300000e+02 0 20 0 -1.1837709e-06 0
2.4400000e+02 0 0 0 -1.1692440e-06 0
2.4400000e+02 0 20 0 -1.1692440e-06 0
2.4500000e+02 0 0 0 -1.1549539e-06 0
2.4500000e+02 0 20 0 -1.1549539e-06 0
2.4600000e+02 0 0 0 -1.1408959e-06 0
2.4600000e+02 0 20 0 -1.1408959e-06 0
2.4700000e+02 0 0 0 -1.1270652e-06 0
2.4700000e+02 0 20 0 -1.1270652e-06 0
2.4800000e+02 0 0 0 -1.1134572e-06 0
2.4800000e+02 0 20 0 -1.1134572e-06 0
2.4900000e+02 0 0 0 -1.1000674e-06 0
2.4900000e+02 0 20 0 -1.1000674e-06 0
2.5000000e+02 0 0 0 -1.0868916e-06 0
2.5000000e+02 0 20 0 -1.0868916e-06 0
2.5100000e+02 0 0 0 -1.0739254e-06 0
2.5100000e+02 0 20 0 -1.0739254e-06 0
2.5200000e+02 0 0 0 -1.0611648e-06 0
2.5200000e+02 0 20 0 -1.0611648e-06 0
2.5300000e+02 0 0 0 -1.0486055e-06 0
2.5300000e+02 0 20 0 -1.0486055e-06 0
2.5400000e+02 0 0 0 -1.0362438e-06 0
2.5400000e+02 0 20 0 -1.0362438e-06 0
2.5500000e+02 0 0 0 -1.0240757e-06 0
2.5500000e+02 0 20 0 -1.0240757e-06 0
2.5600000e+02 0 0 0 -1.0120974e-06 0
2.5600000e+02 0 20 0 -1.0120974e-06 0
2.5700000e+02 0 0 0 -1.0003052e-06 0
2.5700000e+02 0 20 0 -1.0003052e-06 0
2.5800000e+02 0 0 0 -9.8869563e-07 0
2.5800000e+02 0 20 0 -9.8869563e-07 0
2.5900000e+02 0 0 0 -9.7726506e-07 0
2.5900000e+02 0 20 0 -9.7726506e-07 0
2.6000000e+02 0 0 0 -9.6601007e-07 0
2.6000000e+02 0 20 0 -9.6601007e-07 0
2.6100000e+02 0 0 0 -9.5492732e-07 0
2.6100000e+02 0 20 0 -9.5492732e-07 0
2.6200000e+02 0 0 0 -9.4401351e-07 0
2.6200000e+02 0 20 0 -9.4401351e-07 0
2.6300000e+02 0 0 0 -9.3326544e-07 0
2.6300000e+02 0 20 0 -9.3326544e-07 0
2.6400000e+02 0 0 0 -9.2267997e-07 0
2.6400000e+02 0 20 0 -9.2267997e-07 0
2.6500000e+02 0 0 0 -9.1225405e-07 0
2.6500000e+02 0 20 0 -9.1225405e-07 0
2.6600000e+02 0 0 0 -9.0198466e-07 0
2.6600000e+02 0 20 0 -9.0198466e-07 0
2.6700000e+02 0 0 0 -8.9186889e-07 0
2.6700000e+02 0 20 0 -8.9186889e-07 0
2.6800000e+02 0 0 0 -8.8190387e-07 0
2.6800000e+02 0 20 0 -8.8190387e-07 0
2.6900000e+02 0 0 0 -8.7208681e-07 0
2.6900000e+02 0 20 0 -8.7208681e-07 0
2.7000000e+02 0 0 0 -8.6241497e-07 0
2.7000000e+02 0 20 0 -8.6241497e-07 0
2.7100000e+02 0 0 0 -8.5288566e-07 0
2.7100000e+02 0 20 0 -8.5288566e-07 0
2.7200000e+02 0 0 0 -8.4349629e-07 0
2.7200000e+02 0 20 0 -8.4349629e-07 0
2.7300000e+02 0 0 0 -8.3424427e-07 0
2.7300000e+02 0 20 0 -8.3424427e-07 0
2.7400000e+02 0 0 0 -8.2512712e-07 0
2.7400000e+02 0 20 0 -8.2512712e-07 0
2.7500000e+02 0 0 0 -8.1614239e-07 0
2.7500000e+02 0 20 0 -8.1614239e-07 0
2.7600000e+02 0 0 0 -8.0728767e-07 0
2.7600000e+02 0 20 0 -8.0728767e-07 0
2.7700000e+02 0 0 0 -7.9856062e-07 0
2.7700000e+02 0 20 0 -7.9856062e-07 0
2.7800000e+02 0 0 0 -7.8995895e-07 0
2.7800000e+02 0 20 0 -7.8995895e-07 0
2.7900000e+02 0 0 0 -7.8148042e-07 0
2.7900000e+02 0 20 0 -7.8148042e-07 0
2.8000000e+02 0 0 0 -7.7312283e-07 0
2.8000000e+02 0 20 0 -7.7312283e-07 0
2.8100000e+02 0 0 0 -7.6488402e-07 0
2.8100000e+02 0 20 0 -7.6488402e-07 0
2.8200000e+02 0 0 0 -7.5676191e-07 0
2.8200000e+02 0 20 0 -7.5676191e-07 0
2.8300000e+02 0 0 0 -7.4875442e-07 0
2.8300000e+02 0 20 0 -7.4875442e-07 0
2.8400000e+02 0 0 0 -7.4085954e-07 0
2.8400000e+02 0 20 0 -7.4085954e-07 0
2.8500000e+02 0 0 0 -7.3307530e-07 0
2.8500000e+02 0 20 0 -7.3307530e-07 0
2.8600000e+02 0 0 0 -7.2539977e-07 0
2.8600000e+02 0 20 0 -7.2539977e-07 0
2.8700000e+02 0 0 0 -7.1783106e-07 0
2.8700000e+02 0 20 0 -7.1783106e-07 0
2.8800000e+02 0 0 0 -7.1036730e-07 0
2.8800000e+02 0 20 0 -7.1036730e-07 0
2.8900000e+02 0 0 0 -7.0300670e-07 0
2.8900000e+02 0 20 0 -7.0300670e-07 0
2.9000000e+02 0 0 0 -6.9574747e-07 0
2.9000000e+02 0 20 0 -6.9574747e-07 0
2.9100000e+02 0 0 0 -6.8858787e-07 0
2.9100000e+02 0 20 0 -6.8858787e-07 0
2.9200000e+02 0 0 0 -6.8152620e-07 0
2.9200000e+02 0 20 0 -6.8152620e-07 0
2.9300000e+02 0 0 0 -6.7456080e-07 0
2.9300000e+02 0 20 0 -6.7456080e-07 0
2.9400000e+02 0 0 0 -6.6769001e-07 0
2.9400000e+02 0 20 0 -6.6769001e-07 0
2.9500000e+02 0 0 0 -6.6091226e-07 0
2.9500000e+02 0 20 0 -6.6091226e-07 0
2.9600000e+02 0 0 0 -6.5422595e-07 0
2.9600000e+02 0 20 0 -6.5422595e-07 0
2.9700000e+02 0 0 0 -6.4762956e-07 0
2.9700000e+02 0 20 0 -6.4762956e-07 0
2.9800000e+02 0 0 0 -6.4112159e-07 0
2.9800000e+02 0 20 0 -6.4112159e-07 0
2.9900000e+02 0 0 0 -6.3470054e-07 0
2.9900000e+02 0 20 0 -6.3470054e-07 0
3.0000000e+02 0 0 0 -6.2836498e-07 0
3.0000000e+02 0 20 0 -6.2836498e-07 0

65
run/GPD/vis.mac Normal file
View File

@ -0,0 +1,65 @@
# This is a macro file for visualizing G4 events.
# It can either be included in another macro or called with /control/exec vis.mac
# Create an OpenGL driver (i.e. a scene handler and viewer)
# Some useful choices: VRML2FILE, OGLSX, OGLIX, DAWNFILE, etc.
#/vis/open VRML2FILE
/vis/open OGLIX 600x600-0+0
#*/vis/open DAWNFILE
# To calculate volumes and masses uncomment the next two lines
#*/vis/open ATree
#*/vis/ASCIITree/verbose 4
# Create a new empty scene and attach it to handler
/vis/scene/create
# Add world volume, trajectories and hits to the scene
/vis/scene/add/volume
/vis/scene/add/trajectories
/vis/scene/add/hits
/vis/sceneHandler/attach
# Configure the viewer (optional)
/vis/viewer/set/viewpointThetaPhi 235 -45
/vis/viewer/set/lightsThetaPhi 120 60
#/vis/viewer/set/hiddenEdge true
#/vis/viewer/set/style surface
#/vis/viewer/zoom 0.5
# Style: s - surface, w - wireframe
# Note: "surface style" and "hiddenEdge true" remove transparency!
# Other viewpoints (25 55) (235 -45) (125 35)
# Store trajectory information for visualisation (set to 0 if too many tracks cause core dump)
/tracking/storeTrajectory 1
#At the end of each event (default behaviour)
#/vis/scene/endOfEventAction refresh
#At the end of run of X events - Data from X events will be superimposed
/vis/scene/endOfEventAction accumulate
#At the end of Y runs - Data from Y runs will be superimposed
/vis/scene/endOfRunAction accumulate
# Coloured trajectories for an easier particle identification:
# PDG IDs and colours: e- 11 red, e+ -11 blue, nu_e 12 yellow,
# mu+ -13 magenta, anti_nu_mu -14 green, gamma 22 grey
#
#/vis/modeling/trajectories/create/drawByCharge
#/vis/modeling/trajectories/drawByCharge-0/set 1 cyan
/vis/modeling/trajectories/create/drawByParticleID
#*/vis/modeling/trajectories/drawByParticleID-0/set gamma grey
/vis/modeling/trajectories/drawByParticleID-0/setRGBA mu+ 1 0 1 1
/vis/modeling/trajectories/drawByParticleID-0/setRGBA e+ 0 0 0.8 0.5
#*/vis/modeling/trajectories/drawByParticleID-0/setRGBA nu_e 0.7 0.7 0 1
#*/vis/modeling/trajectories/drawByParticleID-0/setRGBA anti_nu_mu 0.3 1.0 0 0.5
# Verbosity of hits
#/hits/verbose 2
# Output just the detector geometry
/vis/viewer/flush
/vis/reviewKeptEvents

70
run/GPD/visDawn101.mac Normal file
View File

@ -0,0 +1,70 @@
# This is a macro file for visualizing G4 events.
# It can either be included in another macro or called with /control/exec vis.mac
# Create an OpenGL driver (i.e. a scene handler and viewer)
# Some useful choices: VRML2FILE, OGLSX, OGLIX, DAWNFILE, etc.
#/vis/open VRML2FILE
#*/vis/open OGLIX 600x600-0+0
/vis/open DAWNFILE
# To calculate volumes and masses uncomment the next two lines
#*/vis/open ATree
#*/vis/ASCIITree/verbose 4
# Create a new empty scene and attach it to handler
/vis/scene/create
# Add world volume, trajectories and hits to the scene
/vis/scene/add/volume
/vis/scene/add/trajectories
/vis/scene/add/hits
/vis/sceneHandler/attach
# Configure the viewer (optional)
#/vis/viewer/set/viewpointThetaPhi 235 -45
/vis/viewer/set/viewpointThetaPhi 80 30
#/vis/viewer/set/lightsThetaPhi 120 60
#/vis/viewer/set/hiddenEdge true
/vis/viewer/set/style surface
/vis/viewer/zoom 0.8
# Style: s - surface, w - wireframe
# Note: "surface style" and "hiddenEdge true" remove transparency!
# Other viewpoints (25 55) (235 -45) (125 35)
# Store trajectory information for visualisation (set to 0 if too many tracks cause core dump)
/tracking/storeTrajectory 1
#At the end of each event (default behaviour)
#/vis/scene/endOfEventAction refresh
#At the end of run of X events - Data from X events will be superimposed
#cks
#/vis/scene/endOfEventAction accumulate
#At the end of Y runs - Data from Y runs will be superimposed
#/vis/scene/endOfRunAction accumulate
# Coloured trajectories for an easier particle identification:
# PDG IDs and colours: e- 11 red, e+ -11 blue, nu_e 12 yellow,
# mu+ -13 magenta, anti_nu_mu -14 green, gamma 22 grey
#
#/vis/modeling/trajectories/create/drawByCharge
#/vis/modeling/trajectories/drawByCharge-0/set 1 cyan
/vis/modeling/trajectories/create/drawByParticleID
#*/vis/modeling/trajectories/drawByParticleID-0/set gamma grey
#/vis/modeling/trajectories/drawByParticleID-0/setRGBA gamma 1 1 1 0
/vis/modeling/trajectories/drawByParticleID-0/setRGBA mu+ 1 0 0 1
/vis/modeling/trajectories/drawByParticleID-0/setRGBA e+ 0 0 1 1
/vis/modeling/trajectories/drawByParticleID-0/setRGBA gamma 0 1 0 1
/vis/modeling/trajectories/drawByParticleID-0/setRGBA e- 1 0 1 1
/vis/modeling/trajectories/drawByParticleID-0/setRGBA nu_e 1 1 1 0 1
/vis/modeling/trajectories/drawByParticleID-0/setRGBA anti_nu_mu 1 1 1 0.5
#/vis/modeling/trajectories/drawByParticleID-0/set nu_e white
#/vis/modeling/trajectories/drawByParticleID-0/set anti_nu_mu white
# Verbosity of hits
#/hits/verbose 2
# Output just the detector geometry
/vis/viewer/flush

63
run/GPD/visFromToni.mac Normal file
View File

@ -0,0 +1,63 @@
# This is a macro file for visualizing G4 events.
# It can either be included in another macro or called with /control/exec vis.mac
# Create an OpenGL driver (i.e. a scene handler and viewer)
# Some useful choices: VRML2FILE, OGLSX, OGLIX, DAWNFILE, etc.
#/vis/open VRML2FILE
#*/vis/open OGLIX 600x600-0+0
/vis/open DAWNFILE
# To calculate volumes and masses uncomment the next two lines
#*/vis/open ATree
#*/vis/ASCIITree/verbose 4
# Create a new empty scene and attach it to handler
/vis/scene/create
# Add world volume, trajectories and hits to the scene
/vis/scene/add/volume
/vis/scene/add/trajectories
/vis/scene/add/hits
/vis/sceneHandler/attach
# Configure the viewer (optional)
/vis/viewer/set/viewpointThetaPhi 235 -45
/vis/viewer/set/lightsThetaPhi 120 60
#/vis/viewer/set/hiddenEdge true
#/vis/viewer/set/style surface
#/vis/viewer/zoom 0.5
# Style: s - surface, w - wireframe
# Note: "surface style" and "hiddenEdge true" remove transparency!
# Other viewpoints (25 55) (235 -45) (125 35)
# Store trajectory information for visualisation (set to 0 if too many tracks cause core dump)
/tracking/storeTrajectory 1
#At the end of each event (default behaviour)
#/vis/scene/endOfEventAction refresh
#At the end of run of X events - Data from X events will be superimposed
/vis/scene/endOfEventAction accumulate
#At the end of Y runs - Data from Y runs will be superimposed
#/vis/scene/endOfRunAction accumulate
# Coloured trajectories for an easier particle identification:
# PDG IDs and colours: e- 11 red, e+ -11 blue, nu_e 12 yellow,
# mu+ -13 magenta, anti_nu_mu -14 green, gamma 22 grey
#
#/vis/modeling/trajectories/create/drawByCharge
#/vis/modeling/trajectories/drawByCharge-0/set 1 cyan
/vis/modeling/trajectories/create/drawByParticleID
#*/vis/modeling/trajectories/drawByParticleID-0/set gamma grey
/vis/modeling/trajectories/drawByParticleID-0/setRGBA mu+ 1 0 1 1
/vis/modeling/trajectories/drawByParticleID-0/setRGBA e+ 0 0 0.8 0.5
#*/vis/modeling/trajectories/drawByParticleID-0/setRGBA nu_e 0.7 0.7 0 1
#*/vis/modeling/trajectories/drawByParticleID-0/setRGBA anti_nu_mu 0.3 1.0 0 0.5
# Verbosity of hits
#/hits/verbose 2
# Output just the detector geometry
/vis/viewer/flush

63
run/GPD/visVRML.mac Normal file
View File

@ -0,0 +1,63 @@
# This is a macro file for visualizing G4 events.
# It can either be included in another macro or called with /control/exec vis.mac
# Create an OpenGL driver (i.e. a scene handler and viewer)
# Some useful choices: VRML2FILE, OGLSX, OGLIX, DAWNFILE, etc.
/vis/open VRML2FILE
#*/vis/open OGLIX 600x600-0+0
#/vis/open DAWNFILE
# To calculate volumes and masses uncomment the next two lines
#*/vis/open ATree
#*/vis/ASCIITree/verbose 4
# Create a new empty scene and attach it to handler
/vis/scene/create
# Add world volume, trajectories and hits to the scene
/vis/scene/add/volume
/vis/scene/add/trajectories
/vis/scene/add/hits
/vis/sceneHandler/attach
# Configure the viewer (optional)
/vis/viewer/set/viewpointThetaPhi 235 -45
/vis/viewer/set/lightsThetaPhi 120 60
#/vis/viewer/set/hiddenEdge true
#/vis/viewer/set/style surface
#/vis/viewer/zoom 0.5
# Style: s - surface, w - wireframe
# Note: "surface style" and "hiddenEdge true" remove transparency!
# Other viewpoints (25 55) (235 -45) (125 35)
# Store trajectory information for visualisation (set to 0 if too many tracks cause core dump)
/tracking/storeTrajectory 1
#At the end of each event (default behaviour)
#/vis/scene/endOfEventAction refresh
#At the end of run of X events - Data from X events will be superimposed
/vis/scene/endOfEventAction accumulate
#At the end of Y runs - Data from Y runs will be superimposed
/vis/scene/endOfRunAction accumulate
# Coloured trajectories for an easier particle identification:
# PDG IDs and colours: e- 11 red, e+ -11 blue, nu_e 12 yellow,
# mu+ -13 magenta, anti_nu_mu -14 green, gamma 22 grey
#
#/vis/modeling/trajectories/create/drawByCharge
#/vis/modeling/trajectories/drawByCharge-0/set 1 cyan
/vis/modeling/trajectories/create/drawByParticleID
#*/vis/modeling/trajectories/drawByParticleID-0/set gamma grey
/vis/modeling/trajectories/drawByParticleID-0/setRGBA mu+ 1 0 1 1
/vis/modeling/trajectories/drawByParticleID-0/setRGBA e+ 0 0 0.8 0.5
#*/vis/modeling/trajectories/drawByParticleID-0/setRGBA nu_e 0.7 0.7 0 1
#*/vis/modeling/trajectories/drawByParticleID-0/setRGBA anti_nu_mu 0.3 1.0 0 0.5
# Verbosity of hits
#/hits/verbose 2
# Output just the detector geometry
/vis/viewer/flush