Author Topic: Penelope/Penepma/Penfluor/Fanal Monte-Carlo Software and Physics  (Read 22739 times)

Probeman

  • Emeritus
  • *****
  • Posts: 2020
  • Never sleeps...
    • John Donovan
Re: Penelope/Penepma/Penfluor/Fanal Monte-Carlo Software and Physics
« Reply #60 on: February 05, 2020, 12:19:01 pm »
It appears that the NBTH and NBPH parameters are based on the Penepma 2012 "bilayer.in" file which is the default input file for thin film calculations. Apparently the boys in Barcelona must have changed this parameter in Penepma 2016.

So the best solution is to simply edit the bilayer.in file for the new parameter. I've attached a new version below that should work with Penepma 2016 and will upload a new Penepma16.zip file with this modification.

All I did was remove these two lines:

NBTH   45                     [No. of bins for the polar angle THETA]
NBPH   30                   [No. of bins for the azimuthal angle PHI]

And replace them with this line:

NBANGL 45 30              [Nos. of bins for the angles THETA and PHI]

That should work.
« Last Edit: February 05, 2020, 06:43:02 pm by Probeman »
The only stupid question is the one not asked!

John Donovan

  • Administrator
  • Emeritus
  • *****
  • Posts: 2572
  • Other duties as assigned...
    • Probe Software
Re: Penelope/Penepma/Penfluor/Fanal Monte-Carlo Software and Physics
« Reply #61 on: May 06, 2020, 12:06:39 pm »
OK, we've uploaded a ZIP file containing executable and support files for Penepma 2016. You can access the link on this page:

https://probesoftware.com/Technical.html

or you can just grab it with this direct link:

https://probesoftware.com/download/PENEPMA16.ZIP

Here are the short instructions:

Quote
Penepma 2016 is almost identical with Penepma 2012, except that it provides support for emission lines of Li, Be and B. However, Penepma 2016 is not distributed with CalcZAF/Standard because the physics database files for Penepma 2016 are *not* compatible with the Penfluor/Fanal secondary fluorescence application files. Therefore Penepma 2016 is a separate download (here) which can be unzipped and placed in a folder named UserData\Penepma16, if desired for Li, Be and B emission calculations. To utilize the Penepma GUI in the Standard application one should edit the following enties in their Probewin.ini file:

    PENDBASE_Path="C:\Userdata\Penepma16\Pendbase"
    PENEPMA_Path="C:\Userdata\Penepma16\Penepma"
    PENEPMA_Root="C:\Userdata\Penepma16"
    PENEPMA_PAR_Path="C:\Userdata\Penepma12\Penfluor"

The above path to the Penfluor files isn't changed because the Penfluor and Fanal programs won't work with the Penepma 2016 physics database files. 

We probably should re-do the code relating to the Penfluor path in CalcZAF, Standard and PFE, but right now if you want to run the Penfluor/Fanal GUI as described here:

https://probesoftware.com/smf/index.php?topic=58.0

for secondary fluorescence calculations at boundaries, you'll need to re-edit the Penepma paths in the Probewin.ini file back to Penepma12.

In the meantime this Penepma16.ZIP file contains sub folders of output files for pure Li, Be and B at 5, 10, 15, 20 and 25 keV. See the Penepma16\penepma\Batch folder.
« Last Edit: May 06, 2020, 12:11:31 pm by John Donovan »
John J. Donovan, Pres. 
(541) 343-3400

"Not Absolutely Certain, Yet Reliable"

jrminter

  • Professor
  • ****
  • Posts: 57
Re: Penelope/Penepma/Penfluor/Fanal Monte-Carlo Software and Physics
« Reply #62 on: May 07, 2020, 03:33:31 pm »
Many thanks, John. You put a lot of work into this...

Probeman

  • Emeritus
  • *****
  • Posts: 2020
  • Never sleeps...
    • John Donovan
Re: Penelope/Penepma/Penfluor/Fanal Monte-Carlo Software and Physics
« Reply #63 on: May 10, 2020, 10:46:31 am »
Many thanks, John. You put a lot of work into this...

Thank-you John.

I should also mention that in order to calculate intensities down to 50 eV (necessary for Li ka), the number of secondaries generated are too large for the "forcing factor" arrays as specified by default in the Penepma GUI.

For example, here are the default forcing factors generated by selecting the "Optimize Production of Characteristic X-rays" option (from the Cu_Cha.in file):

       >>>>>>>> Interaction forcing.
IFORCE 1 1 4 -10    0.1 1.0           [KB,KPAR,ICOL,FORCER,WLOW,WHIG]
IFORCE 1 1 5 -400   0.1 1.0           [KB,KPAR,ICOL,FORCER,WLOW,WHIG]

It's the FORCER factor in red (-400) that causes the problem, so if you intend to calculate characteristic spectra down to 50 eV or so you should edit this number to -200 as seen here:

       >>>>>>>> Interaction forcing.
IFORCE 1 1 4 -10    0.1 1.0           [KB,KPAR,ICOL,FORCER,WLOW,WHIG]
IFORCE 1 1 5 -200   0.1 1.0           [KB,KPAR,ICOL,FORCER,WLOW,WHIG]


Just FYI.
john
The only stupid question is the one not asked!