[807] | 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 | |
---|