Synrad tutorial: photon simulated desorption
=
[Link to documentation](https://molflow.web.cern.ch/node/139)
# Intro to PSD
Physical process illustration:
![](https://codimd.web.cern.ch/uploads/upload_48f44af47c7d7762115ac247bc1a2b79.png)
![](https://codimd.web.cern.ch/uploads/upload_ddef8987dfe1167f2adbac63eeb5003b.png)
![](https://codimd.web.cern.ch/uploads/upload_0acf1fa41eb63c03efc8e1f6fc440b12.png)
![](https://codimd.web.cern.ch/uploads/upload_a57521b8b1f597d914cf26754e2d74ed.png)
Measurement:
![](https://codimd.web.cern.ch/uploads/upload_63f0bea828c25585d0fc9297a41e84b1.png)
![](https://codimd.web.cern.ch/uploads/upload_37e472ddfa5ffb2a404c31d9b05cf8fd.png)
![](https://codimd.web.cern.ch/uploads/upload_1fae2f0f777378c48d9437e9ec474dec.png)
PSD data:
![](https://codimd.web.cern.ch/uploads/upload_af5c4cae7c9ffe56c6c5ba7e9ecb3225.png)
![](https://codimd.web.cern.ch/uploads/upload_b1abf8327536fcb704f1df6e9885a139.png)
# PSD in Synrad/Molflow
![](https://codimd.web.cern.ch/uploads/upload_00c7182b0e17616b40c0b0c0aa71fb9f.png)
# Extruding the target to create a box
![](https://codimd.web.cern.ch/uploads/upload_a3a78b98de74ae32a25565b9bc9da3f2.png)
![](https://codimd.web.cern.ch/uploads/upload_865837cda77ab96b0c47004a3d36a2c8.png)
The target is already pointing towards the source, so we extrude it by 300cm (it's at Z=200cm)
![](https://codimd.web.cern.ch/uploads/upload_97758385f52484d19fabe5452c59b597.png)
And we ge ta closed box:
![](https://codimd.web.cern.ch/uploads/upload_71f2c3e225063f6c5da88244b0c026f4.png)
Let's run the simulation and save the file.
![](https://codimd.web.cern.ch/uploads/upload_118749de1961776f07192848ece3f722.png)
# Inserting in Molflow
We simply open the original Synrad file:
![](https://codimd.web.cern.ch/uploads/upload_db4bd0f7805fd53799e135ca55ebe1b3.png)
Molflow can open Synrad geometries:
![](https://codimd.web.cern.ch/uploads/upload_defaad38d1d2dfed189cb6e6487686cf.png)
Then we "import" desorption (start the PSD conversion process):
![](https://codimd.web.cern.ch/uploads/upload_36e7690515d17a88914943e87807e942.png)
Since we opened the Synrad file, we'll choose "Use current file":
![](https://codimd.web.cern.ch/uploads/upload_e648424033746a501fc3f082da14b8d7.png)
And we'll imitate a simple curve by a simple equation:
![](https://codimd.web.cern.ch/uploads/upload_b1abf8327536fcb704f1df6e9885a139.png =x300)
![](https://codimd.web.cern.ch/uploads/upload_3ae457a6787959644ab306da7100b71a.png)
Time is 1h (3600 seconds), the longer the system conditions, the better:
![](https://codimd.web.cern.ch/uploads/upload_2de9c3545a5c7ab3c25437c75ac215cb.png)
![](https://codimd.web.cern.ch/uploads/upload_1880aa9d4957cc707ddefa4ed2d90c8d.png)
The facet advanced parameters show that it worked and give some diagnostic info:
![](https://codimd.web.cern.ch/uploads/upload_22734ab2587ca2481182fafdc5522d7e.png)
# Run a simulation
We currently want to diagnose PSD, so we'll set sticking=1 on all facets:
![](https://codimd.web.cern.ch/uploads/upload_1ed651dec78574f9373cb53c244f1504.png)
The simulation looks more or less fine, particles come from where they should:
![](https://codimd.web.cern.ch/uploads/upload_6e811e414fc473898cd893723b664dd4.png)
Turning on Hits view and zooming on the source, we can see that molecules are created at irradiation locations:
![](https://codimd.web.cern.ch/uploads/upload_905df003240e891620e78965d645fac8.png)
Let's, however, create a desorption counting texture to visualize desorption:
![](https://codimd.web.cern.ch/uploads/upload_009997ba05ea3380dc44fe856777e1b2.png)
# Result
On linear texture scaling (Synrad vs Molflow):
![](https://codimd.web.cern.ch/uploads/upload_8220c0c7e569e7e1f27b0614ed579fa2.png)
Flux is "compressed":
![](https://codimd.web.cern.ch/uploads/upload_b1abf8327536fcb704f1df6e9885a139.png)
![](https://codimd.web.cern.ch/uploads/upload_42253858ba0d5ff49f3462c9ea485708.png)
Synrad simulations: several orders of magnitude if possible:
![](https://codimd.web.cern.ch/uploads/upload_da93fca8ade35b614e85cf9cee1368f5.png)