1 | |
---|
2 | # --------------------------------------------------------- |
---|
3 | # THIS MACRO PERMIT TO RECONSTRUCT THE SPREAD OUT BRAGG PEAK |
---|
4 | # MODULATING THE PROTON BEAM ENERGY BY MEANS OF A ROTATING |
---|
5 | # MODULATOR WHEEL. |
---|
6 | # HERE THE WHEEL IS ROTATED OF 1 DEGREE AT EACH RUN AND |
---|
7 | # ONE THOUSAND PROTON HYSTORIES ARE SIMULATED IN EACH |
---|
8 | # RUN |
---|
9 | # IN ORDER TO INCREASE THE RESOLUTION MORE EVENT CAN ADDED |
---|
10 | # OR A SMALLER ANGLE CAN BE CHOSEN. |
---|
11 | # |
---|
12 | # NOTE THAT THE MODULATOR MATERIAL IS POLTMETHYLMETHACRILATE |
---|
13 | # (PMMA) FOR DEFAULT. IF ONE WANT CARRY OUT A SIMULATION WITHOUT |
---|
14 | # THE MODULATOR HE/SHE MUST SET "Air" the <<ModMater>> in the |
---|
15 | # <<GetMater>> function of the HadrontherapyModulator.cc class |
---|
16 | # |
---|
17 | # USERS SHOULD GIVE A LOOK TO THE HELP OF THE IDLE TO KNOW |
---|
18 | # THE ACTIVATED MESSSENGERS FOR THE GEOMETRY |
---|
19 | # |
---|
20 | # ADDITIONAL INFORMATIONS ON THE MESSENGER AVAILABLE CAN BE FOUND |
---|
21 | # INSIDE THE HADRONTHERAPY DOCUMENTATION (http://www.ge.infn.it/geant4/examples/). |
---|
22 | # |
---|
23 | # ANYWAY SEND ME AN E-MAIL FOR ANY QUESTION: cirrone@lns.infn.it. |
---|
24 | # -------------------------------------------------------------------------------- |
---|
25 | |
---|
26 | /control/verbose 1 |
---|
27 | /tracking/verbose 0 |
---|
28 | /run/verbose 0 |
---|
29 | /event/verbose 0 |
---|
30 | |
---|
31 | |
---|
32 | # SETTING FOR THE PHYSICS MODELS |
---|
33 | /physics/addPhysics Decay |
---|
34 | /physics/addPhysics EM-Photon-EPDL |
---|
35 | /physics/addPhysics EM-Electron-EEDL |
---|
36 | /physics/addPhysics EM-Positron-Standard |
---|
37 | /physics/addPhysics EM-HadronIon-LowE |
---|
38 | /physics/addPhysics EM-Muon-Standard |
---|
39 | /physics/addPhysics HadronicEl-HadronIon-LElastic |
---|
40 | /physics/addPhysics HadronicInel-ProtonNeutron-LEP |
---|
41 | /physics/addPhysics HadronicInel-Ion-LEP |
---|
42 | /physics/addPhysics HadronicInel-Pion-LEP |
---|
43 | /physics/addPhysics HadronicAtRest-MuonMinus-Capture |
---|
44 | |
---|
45 | /run/initialize |
---|
46 | |
---|
47 | # SET OF THE VISUALISATION CHARACTERISTICS |
---|
48 | /vis/scene/create |
---|
49 | /vis/open OGLIX |
---|
50 | /vis/viewer/zoom 3.5 |
---|
51 | /vis/viewer/set/viewpointThetaPhi 50 130 deg |
---|
52 | /vis/viewer/pan -60 40 cm |
---|
53 | /tracking/storeTrajectory 1 |
---|
54 | /vis/scene/endOfEventAction accumulate |
---|
55 | |
---|
56 | # IF YOU WANT A RANGE SHIFTER PLEASE CHANGE |
---|
57 | # THE FOLLOWIG PARAMETERS |
---|
58 | # SETTING FOR THE RANGE SHIFTER |
---|
59 | |
---|
60 | /beamLine/RangeShifter/RSMat Air |
---|
61 | /beamLine/RangeShifter/thickness 1 cm |
---|
62 | |
---|
63 | |
---|
64 | # START OF THE MODULATOR ROTATION |
---|
65 | # IT CONSISTS IN A COMPLETE ROTATION |
---|
66 | # OF THE WHHEEL (360 DEG) |
---|
67 | |
---|
68 | /run/beamOn 10 |
---|
69 | /vis/viewer/flush |
---|
70 | |
---|
71 | /modulator/angle 1 deg |
---|
72 | /run/beamOn 10 |
---|
73 | /vis/viewer/flush |
---|
74 | |
---|
75 | /modulator/angle 1 deg |
---|
76 | /run/beamOn 10 |
---|
77 | /vis/viewer/flush |
---|
78 | |
---|