<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p><span style="font-size:12pt;" dir="ltr" id="divtagdefaultwrapper"></p>
<div style="margin-top:0;margin-bottom:0;">Dear John, <br>
</div>
<div style="margin-top:0;margin-bottom:0;"><br>
</div>
<div style="margin-top:0;margin-bottom:0;">Thank you for your answer! <br>
</div>
<div style="margin-top:0;margin-bottom:0;"><br>
</div>
<div style="margin-top:0;margin-bottom:0;">I tried to set the step function as follows:
<br>
</div>
<div style="margin-top:0;margin-bottom:0;"> </div>
<div>/gate/physics/addPhysicsList QGSP_BIC_EMZ<br>
/gate/run/initialize <br>
</div>
<div>/process/eLoss/StepFunction <font size="3" face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols" color="black">
<span style="font-size:12pt;" id="divtagdefaultwrapper">0.05 0.1 mm</span></font><br>
/process/eLoss/StepFunctionMuHad <font size="3" face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols" color="black">
<span style="font-size:12pt;" id="divtagdefaultwrapper">0.05 0.1 mm</span></font></div>
<div style="margin-top:0;margin-bottom:0;"> </div>
<div style="margin-top:0;margin-bottom:0;"><br>
</div>
<div style="margin-top:0;margin-bottom:0;">Since the simulation output on the physics did not change I added other commands, which were displayed correctly in the physics output:
<br>
</div>
<div style="margin-top:0;margin-bottom:0;">/process/eLoss/linLossLimit 0.44</div>
<div style="margin-top:0;margin-bottom:0;">/process/msc/RangeFactor 0.18<br>
</div>
<div style="margin-top:0;margin-bottom:0;">/process/msc/RangeFactorMuHad 0.18</div>
<div style="margin-top:0;margin-bottom:0;"><br>
</div>
<div style="margin-top:0;margin-bottom:0;">output:
<div>
<div>msc: for e- SubType= 10<br>
<b>RangeFactor= 0.18</b>, stepLimitType: 3, latDisplacement: 1, skin= 1, geomFactor= 2.5<br>
===== EM models for the G4Region DefaultRegionForTheWorld ======<br>
UrbanMsc : Emin= 0 eV Emax= 100 MeV Table with 120 bins Emin= 100 eV Emax= 100 MeV<br>
WentzelVIUni : Emin= 100 MeV Emax= 10 TeV Table with 100 bins Emin= 100 MeV Emax= 10 TeV<br>
<br>
eIoni: for e- SubType= 2<br>
dE/dx and range tables from 100 eV to 10 TeV in 220 bins<br>
Lambda tables from threshold to 10 TeV, 20 bins per decade, spline: 1<br>
<b>finalRange(mm)= 0.1, dRoverRange= 0.2</b>, integral: 1, fluct: 1, <b>linLossLimit= 0.44</b><br>
===== EM models for the G4Region DefaultRegionForTheWorld ======<br>
PenIoni : Emin= 0 eV Emax= 1 MeV<br>
MollerBhabha : Emin= 1 MeV Emax= 10 TeV deltaVI</div>
<br>
</div>
<div><br>
</div>
<div>Do you have any idea, why this is the case? Did I miss something? Thank you in advance.
<br>
</div>
<div><br>
</div>
<div>Best regards, <br>
</div>
<div>Jana <br>
</div>
<div><br>
</div>
<div>P.S. Sorry to everyone for sending this email twice!<br>
</div>
</div>
</span><br>
<p></p>
</div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>Von:</b> Jana Kretschmer<br>
<b>Gesendet:</b> Mittwoch, 18. Dezember 2019 15:57:04<br>
<b>An:</b> John Apostolakis<br>
<b>Cc:</b> gate-users@lists.opengatecollaboration.org<br>
<b>Betreff:</b> AW: [Gate-users] Change physics list/step function</font>
<div> </div>
</div>
<div><style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p>Dear John, <br>
</p>
<p><br>
</p>
<p>Thank you for your answer! <br>
</p>
<p><br>
</p>
<p>I tried to set the step function as follows: <br>
</p>
<p><br>
</p>
<p></p>
<div>/gate/physics/addPhysicsList QGSP_BIC_EMZ<br>
/gate/run/initialize <br>
</div>
<div>/process/eLoss/StepFunction <font size="3" face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols" color="black">
<span style="font-size:12pt;" id="divtagdefaultwrapper">0.05 0.1 mm</span></font><br>
/process/eLoss/StepFunctionMuHad <font size="3" face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols" color="black">
<span style="font-size:12pt;" id="divtagdefaultwrapper">0.05 0.1 mm</span></font></div>
<p></p>
<p><br>
</p>
<p>Since the simulation output on the physics did not change I added other commands, which were displayed correctly in the physics output:
<br>
</p>
<p><span>/process/eLoss/linLossLimit 0.44</span></p>
<p><span>/process/msc/RangeFactor 0.18</span><br>
</p>
<p><span>/process/msc/RangeFactorMuHad 0.18</span></p>
<p><span><br>
</span></p>
<p><span></span>output: <br>
</p>
</div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>Von:</b> John Apostolakis <John.Apostolakis@cern.ch><br>
<b>Gesendet:</b> Dienstag, 17. Dezember 2019 12:07:43<br>
<b>An:</b> Jana Kretschmer<br>
<b>Cc:</b> gate-users@lists.opengatecollaboration.org<br>
<b>Betreff:</b> Re: [Gate-users] Change physics list/step function</font>
<div> </div>
</div>
<div>Dear Jana,
<div class=""><br class="">
</div>
<div class="">If you are using any physics list with the ending _EMZ, then the electromagnetic processes is instantiated by the G4PhysicsFactory in the following line:</div>
<div class=""><br class="">
</div>
<div class="">
<div class=""> } else if(4 == em_opt) {</div>
<div class=""><span class="Apple-tab-span" style="white-space:pre"></span>p->ReplacePhysics(new G4EmStandardPhysics_option4(verbose)); </div>
</div>
<div class=""> } … </div>
<div class=""><br class="">
</div>
<div class="">The capability of GATE to create each process of a physics list via UI command uses different code, so you cannot mix and match between these options.</div>
<div class=""><br class="">
</div>
<div class="">But after you define the physics list QGSP_BIC_EMZ, it is possible to change the properties of the ionisation process, using the equivalent Geant4 UI command, which I believe is in one of the two directories:</div>
<div class=""><br class="">
</div>
<div class=""> /process/eLoss</div>
<div class=""> /process/em</div>
<div class=""><br class="">
</div>
<div class="">You may need to experiment with the order in which you place these UI commands.</div>
<div class=""><br class="">
</div>
<div class="">Best regards,</div>
<div class="">
<div class=""><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-variant-ligatures: normal; font-variant-east-asian: normal; font-variant-position: normal; line-height: normal; border-spacing: 0px; -webkit-text-decorations-in-effect: none;">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; border-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-stroke-width: 0px;">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; border-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-stroke-width: 0px;">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; border-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-stroke-width: 0px;">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; border-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-stroke-width: 0px;">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; border-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-stroke-width: 0px;">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div class="">
<div class="">John</div>
<div class="">
<div class="">
<div apple-content-edited="true" class=""><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; border-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-stroke-width: 0px;"><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; border-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-stroke-width: 0px;">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div class="">
<div class="">===================================================</div>
<div class="">John Apostolakis, EP Department, CERN</div>
<div class="">SFT (SoFTware for Experiments) Group</div>
<div class="">Office: Building 32/ room R-003 (ground floor), Mail: J27210</div>
</div>
</div>
</span></span></div>
</div>
<div class="">
<div apple-content-edited="true" class=""><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; border-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-stroke-width: 0px;"><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; border-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-stroke-width: 0px;">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div class="">
<div class=""><a href="mailto:john.apostolakis@cern.ch" class="">Email: john.apostolakis@cern.ch</a> Office Tel: +41-22-767-7239</div>
</div>
</div>
</span></span></div>
</div>
<div class="">
<div class="">------------------------------------------------------------------------------------------</div>
</div>
</div>
</div>
</div>
</span></div>
</span></div>
</span></div>
</span></div>
</span></div>
</span></div>
<div><br class="">
</div>
<div>I note that it is possible also to use C++ code (e.g. revised messengers or elsewhere) to do this also:</div>
<div><br class="">
</div>
<div>
<div class=""> #include “G4EmParameters.hh”</div>
<div class=""> ….</div>
<div class=""> G4EmParametters::SetStepFunction(G4double, G4double)</div>
</div>
<div><br class="">
</div>
<div><br class="">
<blockquote type="cite" class="">
<div class="">On 12 Dec 2019, at 14:33, Jana Kretschmer <<a href="mailto:jana.kretschmer@uni-oldenburg.de" class="">jana.kretschmer@uni-oldenburg.de</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div id="divtagdefaultwrapper" dir="ltr" style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-size: 12pt; font-family: Calibri, Helvetica, sans-serif, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols;" class="">
<div style="margin-top: 0px; margin-bottom: 0px;" class="">Hi everyone,<span class="Apple-converted-space"> </span><br class="">
</div>
<div style="margin-top: 0px; margin-bottom: 0px;" class=""><br class="">
</div>
<div style="margin-top: 0px; margin-bottom: 0px;" class="">I am facing a problem concerning the physics settings. I would like to use the physics list<span class="Apple-converted-space"> </span><span class="">QGSP_BIC_EMZ but change the step function for electron<span class="Apple-converted-space"> </span><span class="">ionisation<span class="Apple-converted-space"> </span><span class="">(dRoverRange
and finalRange)</span></span>.<span class="Apple-converted-space"> </span><span class="">Is there anyone who has done this before or knows how I can do this in Gate?<span class="Apple-converted-space"> </span></span><br class="">
</span></div>
<div style="margin-top: 0px; margin-bottom: 0px;" class=""><span class=""><br class="">
</span></div>
<div style="margin-top: 0px; margin-bottom: 0px;" class=""><span class="">I tried to do this via the commands:<span class="Apple-converted-space"> </span><br class="">
</span></div>
<div style="margin-top: 0px; margin-bottom: 0px;" class=""><span class=""><span class="">/gate/physics/addPhysicsList QGSP_BIC_EMZ</span></span></div>
<div style="margin-top: 0px; margin-bottom: 0px;" class=""><span class=""><span class="">/gate/physics/processes/ElectronIonisation/setStepFunction e- 0.05 0.1 mm</span></span></div>
<div style="margin-top: 0px; margin-bottom: 0px;" class=""><span class=""><span class=""></span></span><br class="">
</div>
<div style="margin-top: 0px; margin-bottom: 0px;" class="">Unfortunately this did not work<span class=""></span>, since my simulation output was:<span class="Apple-converted-space"> </span><br class="">
</div>
<p style="margin-top: 0px; margin-bottom: 0px;" class=""></p>
<div class="">eIoni: for e- SubType= 2<br class="">
dE/dx and range tables from 100 eV to 10 TeV in 220 bins<br class="">
Lambda tables from threshold to 10 TeV, 20 bins per decade, spline: 1<br class="">
<span class="Apple-converted-space"> </span><b class="">finalRange(mm)= 0.01, dRoverRange= 0.2</b>, integral: 1, fluct: 1, linLossLimit= 0.01<br class="">
===== EM models for the G4Region DefaultRegionForTheWorld ======<br class="">
PenIoni : Emin= 0 eV Emax= 1 MeV<br class="">
MollerBhabha : Emin= 1 MeV Emax= 10 TeV deltaVI</div>
<br class="">
<p style="margin-top: 0px; margin-bottom: 0px;" class=""></p>
<div style="margin-top: 0px; margin-bottom: 0px;" class="">I also tried to build my own physics list by adding all processes used in<span class="Apple-converted-space"> </span><span class="">QGSP_BIC_EMZ</span><span class="Apple-converted-space"> </span>and
then setting the step function but I was not able to specify the physics models<span class=""><span class=""><span class="Apple-converted-space"> </span>needed</span></span>.</div>
<div style="margin-top: 0px; margin-bottom: 0px;" class=""><br class="">
</div>
<div style="margin-top: 0px; margin-bottom: 0px;" class="">Best regards,<span class="Apple-converted-space"> </span><br class="">
</div>
<div style="margin-top: 0px; margin-bottom: 0px;" class="">Jana<br class="">
</div>
<div style="margin-top: 0px; margin-bottom: 0px;" class=""><br class="">
</div>
</div>
<span style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;" class="">_______________________________________________</span><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<span style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;" class="">Gate-users
mailing list</span><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<a href="mailto:Gate-users@lists.opengatecollaboration.org" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px;" class="">Gate-users@lists.opengatecollaboration.org</a><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px;" class="">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a></div>
</blockquote>
</div>
<br class="">
</div>
</div>
</div>
</body>
</html>