20:01:28 発信元: *Tutor* Tatsuhiko Ogawa : http://phits.jaea.go.jp/tutorial/rec.html 20:03:54 発信元: 005 A. MAITRALLAIN : Hello, I cannot access the ;http://phits.jaea.go.jp/tutorial/rec.html link it shows a 404 error Thank you 20:04:48 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "Hello, I cannot acce..." Try to access here https://phits.jaea.go.jp/tutorial/rec.html 20:05:03 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "Hello, I cannot acce..." Not http but https 20:05:24 発信元: 005 A. MAITRALLAIN : Reacted to "Try to access here h..." with 👍 20:11:45 発信元: 003M.Tzivaki : how would one go about multiplying the t-deposit tally with energy dependent dose conversion factors to get e.g. effective dose? 20:11:55 発信元: 079 Sujoy Chatterjee : Will therebe any lecture on basics of ANGEL ? 20:12:58 発信元: 207 Poonamjot : Can we calculate attenuation factor linear or mass attenuation directly ??? 20:13:24 発信元: 002 Worathat Paenthong Main : Is it possible to calculate the energy deposition in target that place in different position? 20:14:10 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "Will therebe any lec..." There is no lecture material specific to ANGEL. Please see phits/manual/manualE-angel4.pdf. 20:14:36 発信元: *Tutor* Yusuke Matsuya : Replying to "how would one go abo..." To get the results of effective dose, you can use [t-track] tally with multiplier function. 20:15:02 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "how would one go abo..." Here is an example of multiplier to get effective dose. [ T - T r a c k ] title = Track Detection u mesh = xyz # x-type = 1 # nx = 1 # -5.0 5.0 y-type = 1 # ny = 1 # -5.0 5.0 z-type = 2 # zmin = 30.00000 # zmax = 40.00000 # nz = 50 # e-type = 1 # ne = 1 # 0.0 5000.0 unit = 1 # axis = z # file = track_z.out # part = all gshow = 3 # epsout = 2 # y-txt = Effective dose [p multiplier = all mat mset1 all ( 1 -202 ) 20:15:31 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "Can we calculate att..." Attenuation cofficient cannot be calculated directly. Please calculate it from the change in fluence obtained by t-track or t-cross tally. 20:16:00 発信元: *Tutor* Takuya Furuta : Replying to "Is it possible to ca..." You can do that by using mesh=reg. 20:16:20 発信元: 207 Poonamjot : Replying to "Can we calculate att..." Okay sir, thankyou 20:16:24 発信元: 003M.Tzivaki : Antworten an "how would one go a..." To get the results of effective dose, you can use [t-track] tally with multiplier function. ah yes that makes sense, I would need the flux for this and not the deposited energy. And thank you for the example! 20:16:56 発信元: *Tutor* Takuya Furuta : Replying to "Is it possible to ca..." Or for arbitrary location, you can create more than one tally so that you can get the dose at different location. 20:16:56 発信元: 079 Sujoy Chatterjee : Replying to "Will therebe any lec..." Thank you 20:17:04 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "Thank you " with 👍 20:17:11 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "To get the results o..." with 👍 20:17:16 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "Okay sir, thankyou " with 👍 20:20:06 発信元: 119 A.F.Jundi Main : Does the numerical value of the deposit energy is in the file deposit.out? If so, how to read it? 20:20:28 発信元: 144 Rhakamerta Hijazi : what is the meaning of Gray/source? should we multiply it by our charge and batch? 20:20:57 発信元: *Tutor* Seiki Ohnishi : Replying to "Does the numerical v..." Yes it is. You can read it by some text editors as you like. 20:21:50 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "what is the meaning ..." No. If you want Gy, multiply (Gy/source ) by the number of source particles in your irradiation. 20:23:02 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "what is the meaning ..." This is discussed in exercise 6 and 7 of basic lecture 3 (next lecture). 20:23:06 発信元: 119 A.F.Jundi Main : Replying to "Does the numerical v..." I mean, how to interpret five columns consists of numbers in the deposit.out file? 20:23:11 発信元: *Tutor* Yusuke Matsuya : Replying to "what is the meaning ..." Just in case, “Gy/source” means that the absorbed dose per generated source you defined in the [source] section. 20:23:44 発信元: 164 P.Davey : Reacted to "Just in case, “Gy/so..." with 👍 20:23:55 発信元: 144 Rhakamerta Hijazi : Replying to "what is the meaning ..." oh alright, thank you so much sensei 20:24:02 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "oh alright, thank yo..." with 👍 20:27:26 発信元: *Tutor* Yusuke Matsuya : Replying to "Does the numerical v..." In the case of 2D doses like excercise 11 (or 12?), the reading the specific target is very complicated… Therefore, when outputting in two dimensions, I feel it is better to evaluate the dose as an image in EPS format rather than simply reading the numbers. Of course, it is possible to find the dose you want. 20:28:00 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "Does the numerical v..." If you are referring to line 125 or later of deposit.out, please take a look at the 2d-type parameters in the manual, where the data format is written. 20:28:15 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "Does the numerical v..." The 4 column format h: n x y(all ),hh0l n # z-lower z-upper all r.err -2.0000E+01 -1.9500E+01 0.0000E+00 0.0000 -1.9500E+01 -1.9000E+01 0.0000E+00 0.0000 -1.9000E+01 -1.8500E+01 0.0000E+00 0.0000 -1.8500E+01 -1.8000E+01 0.0000E+00 0.0000 -1.8000E+01 -1.7500E+01 0.0000E+00 0.0000 -1.7500E+01 -1.7000E+01 0.0000E+00 0.0000 -1.7000E+01 -1.6500E+01 0.0000E+00 0.0000 -1.6500E+01 -1.6000E+01 0.0000E+00 0.0000 -1.6000E+01 -1.5500E+01 0.0000E+00 0.0000 -1.5500E+01 -1.5000E+01 0.0000E+00 0.0000 -1.5000E+01 -1.4500E+01 0.0000E+00 0.0000 -1.4500E+01 -1.4000E+01 0.0000E+00 0.0000 -1.4000E+01 -1.3500E+01 0.0000E+00 0.0000 -1.3500E+01 -1.3000E+01 0.0000E+00 0.0000 -1.3000E+01 -1.2500E+01 0.0000E+00 0.0000 1st column is z coordinate lower boundary 2nd column is z coordinate upper boundary 3rd column is dose. 4th is statistical uncertainty 20:28:41 発信元: 003M.Tzivaki : Antworten an "Does the numerical..." I have found a number of python tools available on the internet to help processing of phits mesh files (for example GitHub - Lindt8/PHITS-Tools: Python module for parsing and post-processing output from the PHITS code) 20:28:54 発信元: 003M.Tzivaki : Antworten an "Does the numerical..." is there any one you would recommend? 20:30:17 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "Does the numerical v..." It is a powerful tool to post-process tally outputs. Some details are written here https://meteor.nucl.kyushu-u.ac.jp/phitswiki/phits_and_dchain_tools 20:30:32 発信元: 003M.Tzivaki : Hat auf "It is a powerful t..." mit 👍 reagiert 20:31:53 発信元: 119 A.F.Jundi Main : Replying to "Does the numerical v..." Thank you so much all 20:32:02 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "Thank you so much al..." with 👍 20:37:56 発信元: 068 ahid nurmanjaya : Reacted to "It is a powerful too..." with 👍 20:42:08 発信元: 003M.Tzivaki : looking at the .out file at the very and there are two parameters listed that are called "seed". Could you please explain briefly how random number generation and seeding is done in phits? can the seed be changed? 20:45:28 発信元: *Tutor* Takuya Furuta : Replying to "looking at the .out ..." It is pseudo-random number created by xorshit64 algorithm. In default setting 6.647299061401E+12 is used for the random seed. 20:46:36 発信元: 132_K.Skerratt-Love_Main : struggling to load the deposit.out file in ANGEL 20:46:36 発信元: *Tutor* Takuya Furuta : Replying to "looking at the .out ..." But you can change using "rseed" or "bitsrseed" parameters in [parameter] section. Please see the PHITS manual for details. We will also come to this in the next lecture lec03. 20:46:48 発信元: 003M.Tzivaki : Hat auf "But you can change..." mit 👍 reagiert 20:46:54 発信元: 003M.Tzivaki : Antworten an "looking at the .ou..." thank you! 20:47:15 発信元: *Tutor* Takuya Furuta : Replying to "struggling to load t..." Do you need a breakout room ? 20:47:22 発信元: 132_K.Skerratt-Love_Main : yes please 20:48:06 発信元: 097 Evan Lambert : Yuho's microphone is muted 20:48:14 発信元: 005 A. MAITRALLAIN : Reacted to "Yuho's microphone is..." with 👍 20:48:39 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "Yuho's microphone is..." with 👍 20:51:32 発信元: 185.W Asztalos : This animation is the best feature ever 20:51:41 発信元: 003M.Tzivaki : Hat auf "This animation is ..." mit 👍 reagiert 20:51:54 発信元: 144 Rhakamerta Hijazi : how do you know tend = 4, and frame = 10? Are there any parameters setup about Tend and frame in our tally? 20:52:52 発信元: 003M.Tzivaki : When I tried loading the tracks for the homework in PHIG-3D the picture looked not at all how I expected. Why would the cylindrical source be rectangular in the xy plane? I also tried to make a xymesh and that looked as expected. Does it have with the way tallies are defined? 20:53:02 発信元: 119 A.F.Jundi Main : Reacted to "This animation is th..." with 👍 20:53:43 発信元: 003M.Tzivaki : Antworten an "When I tried loadi..." unfortunately I can't upload the picture in the chat 20:55:33 発信元: 097 Evan Lambert : I am missing the XYZ-mesh Tally button in PHIG-3D. Is there something missing in my input file? 20:56:19 発信元: *Tutor* Seiki Ohnishi : Replying to "how do you know tend..." You can check the time info by the information button at the left of the "clear" button. The "clear" button lies on the almost top of the particle track pane. 20:57:02 発信元: 003M.Tzivaki : Antworten an "When I tried loadi..." oh thank you so much! just got my answer 20:57:53 発信元: *Tutor* Yusuke Matsuya : Replying to "When I tried loading..." I just noticed this. Resolved already? 20:58:01 発信元: 003M.Tzivaki : Hat auf "I just noticed thi..." mit 👍 reagiert 20:58:37 発信元: *Tutor* Yusuke Matsuya : Replying to "When I tried loading..." OK. Enjoy today’s tutorials! 20:59:27 発信元: *Tutor* Seiki Ohnishi : Replying to "I am missing the XYZ..." On the left side of the PHIG-3D window, there are stacks of panes "Cell", "Setting", "XYZ-mesh", and "Particle tracks". If you click the title of them, the contents will be switched. 20:59:54 発信元: 003M.Tzivaki : Antworten an "When I tried loadi..." I am very much, thank you! 21:00:06 発信元: 097 Evan Lambert : Replying to "I am missing the XYZ..." I see the "Cell", "Setting" and "Particle tracks" but I do not see "XYZ-mesh" 21:00:27 発信元: *Tutor* Yusuke Matsuya : Replying to "I am missing the XYZ..." Maybe you use the previous version of phig-3d? 21:00:53 発信元: 097 Evan Lambert : Replying to "I am missing the XYZ..." Oh, possibly, what is the latest version of phig-3d? 21:01:17 発信元: *Tutor* Seiki Ohnishi : Replying to "I am missing the XYZ..." ver 2025.03.18 is the latest. 21:01:51 発信元: 155 Ahadi Damar Prasetya : If I want to simulate a Faraday cup in a chamber for calculating current of a proton beam, then T-Cross is suitable isn't it? 21:01:55 発信元: *Tutor* Yusuke Matsuya : Replying to "I am missing the XYZ..." You can use the “xyz-mesh tally” of phig3d included in phits335. 21:02:14 発信元: 097 Evan Lambert : Replying to "I am missing the XYZ..." I will check to see if everything is updated after today's lecture. Thank you for the information 21:02:53 発信元: 221 David Ferreira : To calculate necessary shielding for a rad source, e.g. dose-rate after a given material, do you recommend T-Point or T-Track? 21:03:22 発信元: 003M.Tzivaki : why is it that t-point tally is so much better than t-track tally for photons and neutrons? I would have thought that the statistics for a single point are worse than for a t-track mesh with fairly large voxels? Or am I misunderstanding something? 21:04:38 発信元: *Tutor* Seiki Ohnishi : Replying to "To calculate necessa..." T-point is NOT recommended. It shows puseudo convergence easily. T-track is stable and has good efficiency. 21:04:43 発信元: 098 Maria Sultana : which tally I should use for double differential calculation? 21:04:50 発信元: *Tutor* Yusuke Matsuya : Replying to "To calculate necessa..." It depends on what physical quantity you want to calculate. For dose, it's t-deposit, and for effective dose, it's t-track (with multiplier). 21:05:15 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "If I want to simulat..." The T-Cross Tally is good for that purpose. 21:05:26 発信元: 132_K.Skerratt-Love_Main : If i have a 60Co source, glass sample, and NaI detector housed in lead, would I need to state these all as materials? And would the geometry for the NaI detector housed in Lead need to be the AND operator? 21:05:45 発信元: *Tutor* Takuya Furuta : Replying to "To calculate necessa..." If you want to compute value at a certain location only, T-point is efficient. On the other hand, if you are interested in distribution, T-track should be used. 21:06:16 発信元: 221 David Ferreira : Replying to "To calculate necessa..." I am calculating dose rate 21:06:26 発信元: 221 David Ferreira : Replying to "To calculate necessa..." after material shielding 21:06:32 発信元: 096 A. Kirov : if the geometry is small ~ mm, can the .eps files show it at mm scale ? 21:06:44 発信元: 144 Rhakamerta Hijazi : i have several questions about tally result, can i go to breakout room please? 21:06:44 発信元: 155 Ahadi Damar Prasetya : Replying to "If I want to simulat..." @*Tutor* SHIGYO Nobuhiro Thanks!😀 21:07:01 発信元: *Tutor* Takuya Furuta : Replying to "To calculate necessa..." If you use dose conversion coefficient (multiplier), you can get the dose from the fluence. 21:07:14 発信元: 221 David Ferreira : Reacted to "It depends on what p..." with 👍🏻 21:07:17 発信元: 221 David Ferreira : Reacted to "T-point is NOT recom..." with 👍🏻 21:07:18 発信元: *Tutor* Takuya Furuta : Replying to "To calculate necessa..." So T-point and T-track can be used. 21:07:21 発信元: 221 David Ferreira : Reacted to "If you want to compu..." with 👍🏻 21:07:23 発信元: *Tutor* Yusuke Matsuya : Replying to "why is it that t-poi..." It is impractical to calculate particle fluence at a specific point or line using the [t-track] tally. To address this, the point estimator tally [t-point] was introduced to estimate these quantities within a short computational time. 21:07:28 発信元: 155 Ahadi Damar Prasetya : Reacted to "i have several quest..." with 👍 21:07:36 発信元: 119 A.F.Jundi Main : Replying to "If I want to simulat..." @155 Ahadi Damar Prasetya what's a Faraday cup, sir? 21:07:54 発信元: 144 Rhakamerta Hijazi : Replying to "If I want to simulat..." you are so good sir 21:07:58 発信元: 221 David Ferreira : Replying to "To calculate necessa..." Both show me different values. Which one could be more accurate? 21:08:20 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "which tally I should..." For this purpose, T-Cross tally is fine. See also phits/sample/benchmark/Iwamoto-JNST2017//NeutronDDX-proton.inp. 21:08:41 発信元: 155 Ahadi Damar Prasetya : Replying to "If I want to simulat..." @119 A.F.Jundi Main it's basically a thing for measuring current of a beam CMIIW 21:09:04 発信元: 155 Ahadi Damar Prasetya : Replying to "If I want to simulat..." @144 Rhakamerta Hijazi You are good also sir 21:09:24 発信元: 119 A.F.Jundi Main : Reacted to "@119 A.F.Jundi Main ..." with 👍 21:10:09 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "if the geometry is s..." Yes. Even nm scale, EPS files are generated and scales are automatically adjusted. 21:10:31 発信元: 226 David Posthuma de Boer - STFC UKRI : In the T-track tally examples we have defined only a single mesh cell (y-type=1, ny=1). Does this mean that all tracks are projected onto the xz-plane? The same for the t-deposit tally with r-type=1. 21:10:38 発信元: 003M.Tzivaki : ah so the t-point tally is an estimator that doesn't just track the actual particles going through, but will estimate the fluence based at a point from more particle histories? 21:11:11 発信元: 098 Maria Sultana : Replying to "which tally I should..." @*Tutor* SHIGYO Nobuhiro Do I have to change all the cell and surface parameter? I'm working with composite materials. 21:11:18 発信元: *Tutor* Yusuke Matsuya : Replying to "If i have a 60Co sou..." It depends on how detailed you want to be. You probably should better to take lead shielding into account, because the spectrum might contain characteristic x-rays from lead. 21:11:24 発信元: 003M.Tzivaki : Antworten an "why is it that t-p..." does it work like the point tally in mcnp with an exclusion radius? 21:12:59 発信元: 132_K.Skerratt-Love_Main : Replying to "If i have a 60Co s..." Just as a basic model for now, as i learn to use the software. 21:13:08 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "which tally I should..." Change the material section if you are using a compound material for the target, and modify the target cell and surface to get the desired target size. 21:13:21 発信元: 098 Maria Sultana : Replying to "which tally I should..." [ S o u r c e ] s-type = 1 # mono-energetic axial source proj = deuteron # kind of incident particle dir = 1.0 # z-direction of beam [cosine] r0 = 0.5 # radius [cm] z0 = 0. # minimum position of z-axis [cm] z1 = 0. # maximum position of z-axis [cm] e0 = 30 # energy of beam [MeV] [ M a t e r i a l ] mat[1] 9Be 1.0 mat[2] 12C 1.0 [ Mat Name Color ] mat name color 1 Beryllium blue 2 Carbon red [ S u r f a c e ] 10 so 100. 11 cz 0.5 12 pz 0 13 pz 0.3216 15 cz 0.5 16 pz 0.3216 17 pz 0.5726 [ C e l l ] 100 0 -10 #102 #104 101 -1 10 102 1 -1.845 -11 12 -13 104 2 -2.26 -15 16 -17 21:13:52 発信元: 098 Maria Sultana : Replying to "which tally I should..." this is my file. how to fix this in NeutronDDx file? 21:14:00 発信元: 003M.Tzivaki : The dchain tally sounds amazing! I see now where time progression can be useful. How many time steps can this handle? And does it also deal with complex decay chains? 21:14:15 発信元: *Tutor* Yusuke Matsuya : Replying to "In the T-track tally..." You can set ny=2 or 3 or more. In that case, for each y plane, 2d results are generated, so please try it. Meanwhile, r-type is for rz-mesh. 21:15:23 発信元: *Tutor* Seiki Ohnishi : Replying to "If i have a 60Co sou..." At least, you should define glass, and NaI materials in the Material section. I don't think Co is needed because the amount of it is very little. How detailed model you need depends on your simulation purpose. Genrally speaking, it is convenient to make detailed model first, then disable part of them later. 21:15:32 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "which tally I should..." mat[1] and mat[2] are not compounds. They are elements. What kind of element do you want to compose? 21:16:50 発信元: 132_K.Skerratt-Love_Main : Replying to "If i have a 60Co s..." OK, thank you, I will try and plan a model. Are there any example scripts for radiation shielding or simialr to this kind of experiment? 21:17:08 発信元: 098 Maria Sultana : Replying to "which tally I should..." @*Tutor* SHIGYO Nobuhiro composite neutron converter with these elements 21:17:50 発信元: *Tutor* Tatsu Sato : Replying to "The dchain tally sou..." DCHAIN can handle as many time step as you want, but larger time step makes computational time longer. The complex decay chain is fully considered in DCHAIN 21:18:09 発信元: *Tutor* Yusuke Matsuya : Replying to "If i have a 60Co sou..." maybe phits/recommendation/DetectorResponse is useful for you. 21:18:51 発信元: 132_K.Skerratt-Love_Main : Reacted to "maybe phits/recomm..." with 👍 21:18:52 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "which tally I should..." So rather than a compound, it involves gluing a plate of Be and a plate of C together? 21:19:00 発信元: *Tutor* Tatsu Sato : Replying to "why is it that t-poi..." Yes! 21:20:01 発信元: 098 Maria Sultana : Replying to "which tally I should..." yes,kind of that 21:20:06 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "which tally I should..." In that case, material, surface, and cell are fine. 21:20:17 発信元: 003M.Tzivaki : Hat auf "Yes!" mit 👍 reagiert 21:21:07 発信元: *Tutor* Tatsu Sato : Replying to "To calculate necessa..." [t-track] is more accurate and reliable, but you need to set a large history number to obtain good statistical uncertainty data. 21:21:13 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "which tally I should..." In the source section, shift z0 and z1 slightly to make z0 = -1 and z1 = -1 to prevent problems. 21:22:22 発信元: 221 David Ferreira : Reacted to "[t-track] is more ac..." with 👍🏻 21:22:36 発信元: 221 David Ferreira : Replying to "To calculate necessa..." Thank you 21:22:49 発信元: 003M.Tzivaki : Antworten an "The dchain tally s..." phantastic! 21:24:33 発信元: 098 Maria Sultana : Replying to "which tally I should..." @*Tutor* SHIGYO Nobuhiro In the NeutronDDX file there has cell section kind of different, how to modify them to calculate double differential cross section? 21:26:42 発信元: 119 A.F.Jundi Main : when using infl command, does the geometry input file must be in the same directory with the main input file? 21:27:42 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "which tally I should..." Cell 1 in that file corresponds to cells 102 and 104 in your file. Similarly, cell 2 corresponds to cell 100, cell 10 corresponds to cell 101, etc. 21:27:59 発信元: *Tutor* Takuya Furuta : It does not have to be. You can specify the file with absolute path or relative path. 21:28:29 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "which tally I should..." It may be easier to understand if you compare both input files in PHIG3D. 21:28:31 発信元: 147 Audric Uy : Can the infl be added into any section? for example after [T i t l e] 21:28:51 発信元: 119 A.F.Jundi Main : Reacted to "It does not have to ..." with 👍 21:29:13 発信元: *Tutor* Takuya Furuta : Replying to "Can the infl be adde..." Yes. But be careful not to specify in the section with OFF. It will be ignored. 21:29:23 発信元: 147 Audric Uy : Replying to "Can the infl be adde..." Ok thank you! 21:30:48 発信元: 098 Maria Sultana : Replying to "which tally I should..." set: c20[10.0] $ radius of sphere set: c30[0.001] $ depth of sample set: c40[0.001] $ radius of sample set: c50[c91/c92*6.02e23*c30] $ target density particles/cm^2 [ S u r f a c e ] 1 cz c40 2 pz 0.0 3 pz c30 99 so c20 21:31:04 発信元: 098 Maria Sultana : Replying to "which tally I should..." how to modify this with my file? 21:33:03 発信元: 185.W Asztalos : for mathematical expressions, are there automatic typing considerations like in C? For example, does having 2 as the exponent versus 2. make a difference in the output? If I recall correctly this isn't a problem in Fortran, but I thought I'd check 21:33:30 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "which tally I should..." That is the surface of the virtual detector. Make sure that c20 is much larger than the thickness of the target. 21:34:31 発信元: *Tutor* Takuya Furuta : Replying to "for mathematical exp..." Are you saying the difference between 2 and 2.0 ? 21:35:06 発信元: 189 Firliyani Rahmatia : how to add milli ampere second (mAs) parameters to the phits input, for example in the case of xray or accelerator? 21:35:36 発信元: 185.W Asztalos : Replying to "for mathematical exp..." Yes. For example, in C language, if you write the expression (5/6) times a floating point number, the integer division is done to evaluate 5/6, so the answer is zero no matter the value of the floating point 21:35:39 発信元: 098 Maria Sultana : Replying to "which tally I should..." @*Tutor* SHIGYO Nobuhiro for double differential cross section calculation for my input file, this virtual surface detector is necessary? 21:35:50 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "how to add milli amp..." Set this in the totfact parameter in the source section. 21:35:53 発信元: 185.W Asztalos : Replying to "for mathematical exp..." So you must use 5.0/6.0 to get the right answer 21:35:54 発信元: *Tutor* Takuya Furuta : Replying to "how to add milli amp..." This is about the normalization. It will come to this topic soon after in this lecture. 21:38:41 発信元: *Tutor* Tatsu Sato : Replying to "for mathematical exp..." In PHITS, it is not necessary to write 5.0/6.0. 5/6 should work. However, in fortran, you have to write 5.0/6.0. This part is slightly different between PHITS and fortran. 21:39:20 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "which tally I should..." The sample file counts the number of particles passing through the virtual detector and calculates a double differential cross section, so I think it will be necessary in your case as well. 21:41:38 発信元: 098 Maria Sultana : Replying to "which tally I should..." @*Tutor* SHIGYO Nobuhiro so will I need two cell and two surface section for the calculation? one is mine and other is the sample file. 21:42:23 発信元: 185.W Asztalos : Replying to "for mathematical exp..." Okay, thank you! 21:45:14 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "which tally I should..." Your input file should have one cell section and one surface section. A cell section describes the virtual detector information for your target and sample, and a surface section describes it as well. 21:47:38 発信元: 098 Maria Sultana : Replying to "which tally I should..." @*Tutor* SHIGYO Nobuhiro that mean I'll have to modify the NeutronDDX cell and surface according to my elements and their diameters. 21:48:04 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "which tally I should..." Yes, that's right. 21:48:54 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "which tally I should..." Sorry, the thing about virtual detectors is incorrect. 21:49:43 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "which tally I should..." In the sample file, the neutrons coming out of the target are counted using t-product tally without using a virtual detector. 21:50:52 発信元: 003M.Tzivaki : what was non again? 21:51:17 発信元: 098 Maria Sultana : Replying to "which tally I should..." set: c20[10.0] $ radius of sphere set: c30[0.001] $ depth of sample set: c40[0.001] $ radius of sample set: c50[c91/c92*6.02e23*c30] $ target density particles/cm^2 , this is fpr one material but for two material how to write this? 21:52:01 発信元: 098 Maria Sultana : Replying to "which tally I should..." *for 21:52:24 発信元: *Tutor* Takuya Furuta : Replying to "what was non again?" It is the command for ANGEL to ignore ("Non") the column to draw the EPS file. 21:52:39 発信元: 003M.Tzivaki : Hat auf "It is the command ..." mit 👍 reagiert 21:53:05 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "which tally I should..." If neutrons are only produced in the Be target, use the information for the Be target only. 21:53:40 発信元: 003M.Tzivaki : Antworten an "what was non again?" ah but the numbers in the non comun? are those the relative error of the volume calculation? 21:53:42 発信元: 119 A.F.Jundi Main : I am designing geometry for input file using SimpleGeo. But until now, I am not sure what length scale the SimpleGeo uses. Is it cm or mm? 21:53:45 発信元: 207 Poonamjot : Sir, i had written the onion geometry code , same as shown here earlier but after execution output is large sphere enclosing other spheres and similar in 2d , or i am not getting same 3d figure..why so?? 21:53:57 発信元: 155 Ahadi Damar Prasetya : Reacted to "I am designing geome..." with 👍 21:54:20 発信元: 155 Ahadi Damar Prasetya : Replying to "I am designing geome..." I was wondering the same, nice question. 21:54:32 発信元: 119 A.F.Jundi Main : Reacted to "I was wondering the ..." with 👍 21:54:40 発信元: 098 Maria Sultana : Replying to "which tally I should..." no, neutrons are finally produce in carbon. should I use C target then? 21:55:28 発信元: 185.W Asztalos : Are statistical errors calculated per batch, or calculated over the total number of histories? 21:58:07 発信元: 185.W Asztalos : Replying to "Are statistical erro..." Specifically, are there two steps to uncertainty calculation: per batch, then averaged over batches? Or, is it done all at the end? For example, in FLUKA, statistical uncertainty is reduced more effectively by increasing the number of cycles (the FLUKA analogue of batches) as opposed to just increasing histories per batch. 21:58:19 発信元: *Tutor* Takuya Furuta : Replying to "what was non again?" With this file, it does not have so much meaning. For normal PHITS results, you get histogram figure where relative error is not shown in that figure. By the way, there is also option to have a histogram with error bars (epsout=2). 21:58:59 発信元: 155 Ahadi Damar Prasetya : Reacted to "Are statistical erro..." with 👍 21:59:38 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "which tally I should..." In this case, we are bonding Be and C plates together, so why not tally the amount of neutrons generated from each target and add them up later? It would be easier to understand if we made two t-product tallies. 22:00:07 発信元: *Tutor* Takuya Furuta : Replying to "Sir, i had written t..." Shall we talk in breakout room ? 22:00:34 発信元: 207 Poonamjot : Replying to "Sir, i had written t..." Okay sir 22:05:03 発信元: 205 Zahidul Islam Apu Main : does maxcas * maxbatch have any effect of number of source?? Also does increasing maxbatch decrease statistical error? 22:05:28 発信元: 096 A. Kirov : Replying to "if the geometry is s..." Thanks, How do you set eps files in mm scale? It did not do it automatically 22:05:35 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "Are statistical erro..." The istdev parameter allows you to change whether the variance is calculated per batch or over the entire history. 22:06:40 発信元: 003M.Tzivaki : from what Tatsuhiko is saying, that column seems to represent the uncertainty. What type of uncertainty measure is shown? Ìs it the standard error 22:06:50 発信元: 185.W Asztalos : Replying to "Are statistical erro..." Okay, awesome 22:09:10 発信元: 156 Desi Listianti main : I want to ask about Tally and History. Can I have a breakout room please? 22:09:16 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "does maxcas * maxbat..." The product of maxcas and maxbch is the number of histories. Increasing maxbch increases the number of histories, which reduces statistical errors. 22:09:19 発信元: 098 Maria Sultana : Replying to "which tally I should..." @*Tutor* SHIGYO Nobuhiro how to do this? I'm not understanding it clearly. 22:11:50 発信元: 003M.Tzivaki : Antworten an "what was non again?" I guess my underlyuing question would be how is this error calculated? 22:14:54 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "which tally I should..." The first t-product tally calculates the neutron DDX emitted from the Be target. In this case, c91 and c92 use Be information. Next, write the second t-product tally to calculate the neutron DDX emitted from the C target. After running PHITS, add the results of the two tallies to obtain the DDX of the composite target. 22:17:03 発信元: 205 Zahidul Islam Apu Main : Replying to "does maxcas * maxbat..." Thank you 22:18:05 発信元: 098 Maria Sultana : Replying to "which tally I should..." how to add the results? 22:18:52 発信元: 098 Maria Sultana : Replying to "which tally I should..." I mean for two material I'll get two PHITS.out file, then how to add them for final result? 22:26:16 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "which tally I should..." Let σ_Be and σ_C be the t-product tally results, n_Be and n_C be the number densities, and t_Be and t_C be the target thickness. The total cross-sectional area is: (n_Be*t_Be*σ_Be + n_C*t_C*σ_C)/(n_Be*t_Be + n_C*t_C) 22:28:43 発信元: 003M.Tzivaki : Antworten an "I am designing geo..." it has been a while since I last used it but I am reasonably sure it is in cm. But you can verify this by making a simple shape in simple geo (sphere) and exporting it to phits and comparing the dimensions (we know phits is in cm so if it's the same number, simplegeo is also in cm). 22:29:34 発信元: *Tutor* SHIGYO Nobuhiro : Reacted to "it has been a while ..." with 👍 22:29:44 発信元: 119 A.F.Jundi Main : If i rerun the input file using istdev=-1 after being terminated by zeroing the remaining batch, will it continued from the last batch? 22:30:11 発信元: 213 Taskin Sultana : stdev = 1 means what? 22:30:29 発信元: *Tutor* Yusuke Matsuya : Replying to "If i rerun the input..." yes 22:30:34 発信元: 155 Ahadi Damar Prasetya : Reacted to "it has been a while ..." with 👍 22:31:53 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "stdev = 1 means what..." A new calculation is begun. Statistical uncertainties are estimated from the standard error of each batch result. 22:32:00 発信元: 119 A.F.Jundi Main : Replying to "If i rerun the input..." thank you 22:32:25 発信元: 119 A.F.Jundi Main : Reacted to "it has been a while ..." with 👍 22:32:30 発信元: 119 A.F.Jundi Main : Replying to "I am designing geome..." ah you have a point. thank you! 22:32:32 発信元: 003M.Tzivaki : I am assuming that the stdcut means the standard deviation expressed as a relative error is that correct? 22:33:23 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "I am assuming that t..." yes 22:33:24 発信元: *Tutor* Takuya Furuta : Replying to "I am assuming that t..." Standard deviation cut-off. 22:33:38 発信元: 213 Taskin Sultana : Replying to "stdev = 1 means what..." Thank you Tutor 22:33:57 発信元: 155 Ahadi Damar Prasetya : Removed a 👍 reaction from "it has been a while ..." 22:33:57 発信元: 003M.Tzivaki : Hat auf "Standard deviation..." mit 👍 reagiert 22:34:06 発信元: 003M.Tzivaki : Hat auf "yes" mit 👍 reagiert 22:39:35 発信元: 098 Maria Sultana : Replying to "which tally I should..." here t-product tally result would be the sum up which we found in output file? 22:40:26 発信元: 098 Maria Sultana : Replying to "which tally I should..." sum over 1.4123E-03 like this? 22:41:46 発信元: *Tutor* SHIGYO Nobuhiro : The sum over value is the result of integrating with respect to energy, so it is the angle differential cross section rather than DDX. 22:43:56 発信元: *Tutor* SHIGYO Nobuhiro : DDX is written for each energy interval above the sum over. 22:47:16 発信元: 003M.Tzivaki : Antworten an "I am assuming that..." I am still very confused about the errors calculated in phits. Is it the standard error or the standard deviation? 22:50:07 発信元: 064 Youngran Kwon : network is unstable. 22:50:10 発信元: 003M.Tzivaki : Antworten an "I am assuming that..." and they are in every case expressed as a relative error yes? so the % value is given in relation to the mean? 22:50:30 発信元: 005 A. MAITRALLAIN : Is totfact necessary a flux or could it be a number of particle per second? 22:50:32 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "I am assuming that t..." The uncertainty output by PHITS is the standard deviation. 22:51:12 発信元: *Tutor* SHIGYO Nobuhiro : It is relative to the average value. 22:54:21 発信元: *Tutor* Seiki Ohnishi : Replying to "Is totfact necessary..." Usually PHITS output a fluix per source particle. So you should multiply the result by the source intensity you have. If the source intensity, ex photons/sec of RI capsule, are known before running PHITS simulation, you can input it as the totfact, then PHITS do multiplication. 22:54:23 発信元: 221 David Ferreira : what is the difference between maxcas and totfact? 22:54:26 発信元: 003M.Tzivaki : Hat auf "It is relative to ..." mit 👍 reagiert 22:54:48 発信元: *Tutor* Takuya Furuta : As it is given in this example, totfact is a factor to multiply the PHITS results. But there is a special case in RI source option such that the results are automatically normalized to per second. You will come to this in the lecture tomorrow (Source A). 22:55:22 発信元: *Tutor* Yusuke Matsuya : Replying to "Is totfact necessary..." Totfact means normalization parameter. 22:55:33 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "what is the differen..." maxcas (and maxbch) are necessary to determine the statistical accuracy. totfact is a parameter for normalizing the results. 22:55:49 発信元: 005 A. MAITRALLAIN : Replying to "Is totfact necessary..." Thank you 22:57:23 発信元: 125 Maheraj Khan : Dr.@*Tutor* Yusuke Matsuya sensei, Would please give me some time today after finishing this session? I've some inquiry about PHITS chem module with you. 22:57:47 発信元: 003M.Tzivaki : Antworten an "I am assuming that..." ah ok! we are assuming a normal distribution yes? 22:59:20 発信元: *Tutor* Yusuke Matsuya : Replying to "Dr.@*Tutor* Yusuke M..." OK. After finishing this session, I will take you to a breakout room. 22:59:33 発信元: 125 Maheraj Khan : Reacted to "OK. After finishing ..." with ❤️ 22:59:52 発信元: 125 Maheraj Khan : Replying to "Dr.@*Tutor* Yusuke M..." @*Tutor* Yusuke Matsuya Thank you Sensei. 23:00:00 発信元: *Tutor* Takuya Furuta : Replying to "I am assuming that t..." Yes 23:02:25 発信元: 131 A.HUGUES : For the exercice 7, we have large relative errors for inners spheres. what is the minimum relative error in order to assume that we are converged? In MCNP it's 10% for F4 et 5% for F5. Do we have the same thinks for PHITS ? 23:02:50 発信元: 097 Evan Lambert : Thank you for the session today! 23:02:52 発信元: 155 Ahadi Damar Prasetya : Thank you! 23:03:02 発信元: 100-Le Huu Loi-Main : Thank you and see you tomorrow! 23:03:03 発信元: 149 M.Rosero : Thank you! 23:03:07 発信元: 147 Audric Uy : Thank you very much! 23:03:09 発信元: *Tutor* Takuya Furuta : Replying to "For the exercice 7, ..." You can think in the same way. 23:03:14 発信元: 003M.Tzivaki : Antworten an "I am assuming that..." So if we have a tally result for example: 23:03:15 発信元: 132_K.Skerratt-Love_Main : is there a way to download the tutorials, since they are only available for upto 4 months? 23:03:17 発信元: *Tutor* Tatsuhiko Ogawa : http://phits.jaea.go.jp/tutorial/rec.html 23:03:25 発信元: 131 A.HUGUES : Reacted to "You can think in the..." with 👍 23:03:27 発信元: 132_K.Skerratt-Love_Main : the video of the tutorials 23:03:31 発信元: *Tutor* Tatsuhiko Ogawa : 1.https://forms.gle/bQjAzdJjXEiNbJhr6 23:03:35 発信元: *Tutor* Tatsu Sato : Replying to "For the exercice 7, ..." We do not recommend any specific value for the appropriate statistical uncertainties because it totally depends on the simulation condition. 23:03:40 発信元: 220 Yusuf Efe Zorlu : Thank you very much sensei. Have a nice day 23:03:46 発信元: 068 ahid nurmanjaya : お疲れ様です 23:03:47 発信元: 125 Maheraj Khan : Thank you very much for this session. 23:03:49 発信元: 066 Serge Takoukam : Thanks, USEFUL 23:03:51 発信元: 032 I.Gordeev main : Thank you very much! 23:03:53 発信元: 207 Poonamjot : Thankyou sir 23:03:57 発信元: *Tutor* Tatsu Sato : Replying to "For the exercice 7, ..." But as Furuta-san wrote, you can think the same way as MCNP 23:03:59 発信元: 042 Melih Tahsin Değer Main : thanks a lot, have a nice day 23:04:01 発信元: 092 Bal Vikram Khatri : Thank you very much 23:04:07 発信元: 178 A.J.Puspitasari : Thank you very much! 23:04:07 発信元: 047 Nina Widiawati : Thank you very much, see you tomorrow 23:04:09 発信元: 016-Chalit Muanglay ("Sub") : Thank you very much. 23:04:11 発信元: 029 Herlan Setiawan : Thank you.. 23:04:11 発信元: 119 A.F.Jundi Main : Thank you very much sensei! 23:04:12 発信元: 216 Chadchayada Samarak : Thank you 23:04:13 発信元: 131 A.HUGUES : Replying to "For the exercice 7, ..." Thanks 23:04:13 発信元: 124 S.Tancharakorn : Thank you very much, see you tomorrow. 23:04:15 発信元: 221 David Ferreira : Replying to "what is the differen..." If I have a source activity, I should set it on maxcas and maxbch? Or I have to use totfact too? 23:04:17 発信元: 181 Muhammad Mada Ramadan Purnama : Thank you very much! 23:04:18 発信元: 022 Devendra Raj Upadhyay : Replying to "Sir, i had written t..." Thank you so much 23:04:23 発信元: 066 Serge Takoukam : 👋 23:04:28 発信元: 030 Iaroslav Bielov : Thank you! See you tomorrow!! 23:04:32 発信元: 003M.Tzivaki : Antworten an "I am assuming that..." in units of Flux: 100 0.01, that means 100 with a stdev of 1% Is that correct? 23:04:33 発信元: 048 Nur Aida : Thank you very much, 23:04:59 発信元: 131 A.HUGUES : Thank you and see you tomorrow!!