| 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 | |
|---|
