20:00:34 発信元: Devendra Raj Upadhyay : Greetings from Nepal 20:01:14 発信元: *Tutor* Tatsuhiko Ogawa : Hello everyone. Here are links http://phits.jaea.go.jp/tutorial/rec.html https://forms.gle/bQjAzdJjXEiNbJhr6 20:01:42 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "Greetings from Nepal" with 👍 20:02:08 発信元: 132_K.Skerratt-Love_Main : Question: I have been practicing geometry and I am unsre of how you move an object to a set distance in tthe void. e.g., sphere is at the origin but I have a square and want to move it 30cm to the right. I can't decipher the [cell] nomenclature to work this out. 20:03:51 発信元: *Tutor* SHIGYO Nobuhiro : Try moving the planes that defines the square in the [surface] section by 30cm. 20:04:51 発信元: 132_K.Skerratt-Love_Main : Replying to "Question: I have b..." what would that line of code look like, based on below? 20:04:53 発信元: 132_K.Skerratt-Love_Main : Replying to "Question: I have b..." [Surface] 10 so 2.5 20 rpp -10 10 -10 10 -10 10 999 so 100. [Cell] 100 1 -2.32 -10 120 2 -11.432 10 -20 1000 0 -999 #all 1001 -1 999 20:06:29 発信元: *Tutor* Seiki Ohnishi : Replying to "Question: I have bee..." 20 RPP 10 30 -10 10 -10 10 20:07:59 発信元: *Tutor* Seiki Ohnishi : Replying to "Question: I have bee..." and surface 10 and 20 do not have overlaping any more, so "120 2 -11.432 -20" works well. 20:09:27 発信元: 132_K.Skerratt-Love_Main : Replying to "Question: I have b..." the sampels stills seem to be overlapped. 20:13:28 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "Question: I have bee..." Can you tell where the overlap is occurring? 20:15:25 発信元: 132_K.Skerratt-Love_Main : Replying to "Question: I have b..." [mat 2] is now a rectangle rather than a square and has an x between -10 to above 20 cm, and z between 0 to 10cm, where the sphere intercepts at the origin 20:17:01 発信元: 205 Zahidul Islam Apu Main : Are emin and dmax depends on incoming or particle particle? 20:17:43 発信元: 079 Sujoy Chatterjee : Is JENDL 5 the default library for PHITS 3.35 ? 20:17:49 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "Question: I have bee..." What happens if you change number 120 in the cell section back to "120 2 -11.432 10 -20"? 20:18:52 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "Are emin and dmax de..." emin() and dmax() apply not only to the incident particles but also to the generated secondary particles. 20:19:27 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "Is JENDL 5 the defau..." yes 20:19:37 発信元: 205 Zahidul Islam Apu Main : Replying to "Are emin and dmax de..." okay 20:20:04 発信元: *Tutor* Tatsu Sato : Replying to "Is JENDL 5 the defau..." Only when you install JENDL-5 using its installer 20:20:07 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "Question: I have bee..." Would you like to discuss this in a breakout room? 20:20:25 発信元: 079 Sujoy Chatterjee : Okay 20:20:43 発信元: 064 Youngran Kwon : Replying to "Are emin and dmax de..." ok 20:21:04 発信元: 064 Youngran Kwon : Replying to "Is JENDL 5 the defau..." yes 20:21:38 発信元: 132_K.Skerratt-Love_Main : Replying to "Question: I have b..." I realised i missed the removal of -10 to 10 in surface section. square is now no longer overlapping and has moved in the X-axis. I will wait till the end of the day to follow-on from this so I can keep up with the current exercises. Thank you Shigyo-san 20:21:54 発信元: *Tutor* SHIGYO Nobuhiro : Reacted to "I realised i missed ..." with 👍 20:23:01 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "Is JENDL 5 the defau..." For JENDL5 neutrons cross sections of 31 target nuclei are available by default. To use JENDL5 for other targets or particles other than neutron, you need to install JENDL-5. For details, please see \phits\lecture\advanced\High-E-N-Phys Exercise 6. 20:23:12 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "okay" with 👍 20:29:35 発信元: 064 Youngran Kwon : Replying to "Is JENDL 5 the defau..." thanks 20:30:54 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "thanks" with 👍 20:31:45 発信元: 098 Maria Sultana : [ Forced Collisions ] part = proton reg fcl 1 -1.0 in NeutronDDX sample file what does this section mean? 20:31:54 発信元: 131 A.HUGUES : why for the lec01 and lec02 exercices we did not use the negs parameter ? 20:31:58 発信元: 185.W Asztalos : we can also notice in deposit.out that the energy deposition changed... 20:32:47 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "[ Forced Collisions ..." It can improve tally statistics of secondary particles generated from a certain target in which the collision probability is very low. 20:33:11 発信元: *Tutor* Yusuke Matsuya : Replying to "why for the lec01 an..." The reason is that the calculation time becomes long. 20:33:33 発信元: 013 A.Gillespie : What is the default cutoff energy? 20:33:38 発信元: 131 A.HUGUES : Reacted to "The reason is that t..." with 🙏 20:33:57 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "why for the lec01 an..." I think this is because they did not look at the detailed behavior of electrons or photons. 20:33:58 発信元: 185.W Asztalos : Replying to "we can also notice i..." nevermind, I was thinking of yesterday's value which was for protons 20:34:19 発信元: *Tutor* Tatsu Sato : Replying to "What is the default ..." 100 keV for electrons and positron when EGS mode is used 20:36:06 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "What is the default ..." For other particles, please see page 45. 20:37:02 発信元: 098 Maria Sultana : Replying to "[ Forced Collisions ..." @*Tutor* SHIGYO Nobuhiro should it change value for double differential cross section calculation? 20:38:10 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "[ Forced Collisions ..." Even when calculating DDX, please leave it as -1.0. 20:38:47 発信元: 064 Youngran Kwon : Replying to "[ Forced Collisions ..." great 20:39:05 発信元: 098 Maria Sultana : Replying to "[ Forced Collisions ..." @*Tutor* SHIGYO Nobuhiro okay. 20:39:44 発信元: *Tutor* SHIGYO Nobuhiro : Reacted to "@*Tutor* SHIGYO Nobu..." with 👍 20:39:56 発信元: 007 V. Mendes : Replying to "[ Forced Collisions ..." Hi, Could you explain why Energy and momemtum are not conserved in COnventional Mode ? 20:40:01 発信元: 098 Maria Sultana : Replying to "[ Forced Collisions ..." [ M a t e r i a l ] m1 Al 1.0 set: c91[2.7] $ density of target [g/cm3] set: c92[27] $ molecular weight [g/mol] [ C e l l ] 1 1 -c91 -1 2 -3 $ Target 2 0 -99 #1 10 -1 99 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 how to write these section for two material? Can you send me a demo written for two material? 20:40:26 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "[ Forced Collisions ..." @007 V. Mendes it is now being explained by Takuya 20:42:07 発信元: 098 Maria Sultana : Replying to "[ Forced Collisions ..." also my projectile is deuteron 20:42:16 発信元: Olin Pinto : Just for curiosity: Is there a way in PHITS to edit or draw geometries using a graphical interface, and assign regions to these bodies? For example, similar to how it's done in FLUKA using FLAIR, where you can draw shapes like RPPs and assign regions without manually writing them in the input file? 20:42:25 発信元: 185.W Asztalos : So, in the default mode, secondary pions, for instance, would not be transported? Or is this only for incident neutrons? 20:43:21 発信元: *Tutor* Tatsuhiko Ogawa : This argument is only for incident nuetrons below 20 MeV. For pion production, this argument is not relevant. 20:43:50 発信元: 185.W Asztalos : Replying to "So, in the default m..." Okay, just wanted to make sure. So for a high energy proton beam, this does not apply. 20:44:12 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "So, in the default m..." Exactly. 20:44:22 発信元: 185.W Asztalos : Replying to "So, in the default m..." Thank you! 20:45:12 発信元: *Tutor* Tatsu Sato : Replying to "Just for curiosity: ..." Actually, the latest version of Flair has function to convert FLUKA input file to PHITS input file (and vice versa) 20:45:16 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "[ Forced Collisions ..." @098 Maria Sultana What do you mean by "two material". Why Don't you write 2 independent inputs and run them separately? 20:45:39 発信元: *Tutor* Tatsu Sato : Replying to "Just for curiosity: ..." We are preparing lecture note for using Flair for making PHITS input file 20:45:43 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "Thank you!" with 👍 20:46:57 発信元: Olin Pinto : Replying to "Just for curiosity: ..." That’s really great Thanks for the answer. 20:48:04 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "That’s really great ..." with 👍 20:49:04 発信元: Elizabeth Musacchio Gonzalez : Reacted to "We are preparing lec..." with 👏 20:49:50 発信元: 221 David Ferreira : Reacted to "We are preparing lec..." with 👍🏻 20:50:00 発信元: 221 David Ferreira : Reacted to "Actually, the latest..." with 👍🏻 20:52:57 発信元: 003M.Tzivaki : Hat auf "We are preparing l..." mit 👍 reagiert 20:55:22 発信元: 079 Sujoy Chatterjee : If event mode is turned ON, photon dose will also be zero, just like neutron dose ? 20:55:40 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "[ Forced Collisions ..." May be like this? Please examine carefully. [ M a t e r i a l ] m1 Al 1.0 m2 C 1.0 set: c91[2.7] $ density of target [g/cm3] set: c92[27] $ molecular weight [g/mol] set: c93[1.8] est: c94[12] [ C e l l ] 1 1 -c91 -1 2 -3 $ Target 2 2 -c93 -1 3 -4 2 0 -99 #1 10 -1 99 set: c20[10.0] $ radius of sphere set: c30[0.001] $ depth of sample set: c31[0.001] 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 4 pz c31 99 so c20 20:56:32 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "If event mode is tur..." yes 20:56:47 発信元: 079 Sujoy Chatterjee : Replying to "If event mode is t..." Thank you 20:56:52 発信元: 098 Maria Sultana : Replying to "[ Forced Collisions ..." @*Tutor* Tatsuhiko Ogawa Here the first material output result like the generated neutron and remain deutron will collide to the second material. That's why I need two material input together. I think separate input won't work there for me. 20:58:36 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "If event mode is tur..." "If event mode is turned ON" -> No "If negs=1 is turned ON" -> Yes. negs for photons is analogous to e-mode for neutrons. 20:59:35 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "[ Forced Collisions ..." I see. The input that Nobu proposed should work perfectly. 20:59:45 発信元: 002 Worathat Paenthong Main : How to calculate the recoil particle from material? 20:59:49 発信元: 098 Maria Sultana : Replying to "[ Forced Collisions ..." @*Tutor* SHIGYO Nobuhiro Thank you. I'll try to work on this. 20:59:49 発信元: 079 Sujoy Chatterjee : Replying to "If event mode is t..." Thanks for clarification 21:00:01 発信元: *Tutor* SHIGYO Nobuhiro : Reacted to "@*Tutor* SHIGYO Nobu..." with 👍 21:00:25 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "May be like this? Pl..." with 👍 21:00:39 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "Thanks for clarifica..." with 👍 21:00:41 発信元: 205 Zahidul Islam Apu Main : Does e-mode =2 replace neutron dose value as proton dose value? Does target actually has proton dose? 21:01:29 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "How to calculate the..." For irradiation of neutrons below 20 MeV, please turn on e-mode =2. For the other projectiles, they are calculated by default. 21:01:38 発信元: 079 Sujoy Chatterjee : Replying to "If event mode is t..." Why my proton dose is lower than then the value shown in the tutorial ? 21:01:58 発信元: 079 Sujoy Chatterjee : Replying to "If event mode is t..." When e=mode is set to 2 21:03:21 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "If event mode is tur..." For 2 reasons. 1, Statistics is not enough. 2, Dose is deposited not only by protons but also by recoil nuclei, deutérons, tritons, which is not included to "proton dose" 21:04:33 発信元: 079 Sujoy Chatterjee : Replying to "If event mode is t..." I am using the same maxcas=1000 as shown in tutorial 21:05:13 発信元: *Tutor* Yuho Hirata : Replying to "Does e-mode =2 repla..." Not only protons but also other secondary particles deposit energy. If you want to calculate all the deposited energy, set part=all. 21:05:47 発信元: khelifi : for spectrum simulation, how we can defind the GEB function 21:06:11 発信元: Audric Zachary Uy : Replying to "If event mode is tur..." I think they mentioned yesterday, it could be because they used an older version of phits to obtain the values in the ppt? Is this correct? 21:07:06 発信元: 205 Zahidul Islam Apu Main : Replying to "Does e-mode =2 repla..." Then why we didn't get any proton dose when e-mode off and also didn't get neutron dose when e-mode=2 21:07:13 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "If event mode is tur..." 079 Sujoy Chatterjee Even so, neutron dose calculated using kerma factor is fast. In contrast, e-mode = 2 is much slower because statistics is incremented only when proton is produced and volume of interest is crossed. 21:08:09 発信元: *Tutor* Takuya Furuta : Replying to "If event mode is tur..." The values in PPT are computed by older version of PHITS. So they can be statistically different. 21:08:19 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "If event mode is tur..." @Audric Zachary Uy Ummm, small différences are due to version difference. But whichever version you are using, dose calculated by emode = 2 is always more accurate but slower. 21:09:09 発信元: Audric Zachary Uy : Replying to "If event mode is tur..." I understand, thank you! 21:09:31 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "I understand, thank ..." with 👍 21:10:27 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "for spectrum simulat..." GEB? Do you mean Gaussian Energy Broadening? If so, please use "dresol" and "dfano" in T-deposit tally. 21:10:47 発信元: 098 Maria Sultana : Replying to "[ Forced Collisions ..." @*Tutor* SHIGYO Nobuhiro [ C e l l ] 1 1 -c91 -1 2 -3 $ Target 2 2 -c93 -1 3 -4 2 0 -99 #1 10 -1 99 here cell 2 write twice, one would be 3,right? 21:11:01 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "for spectrum simulat..." See also phits/lecture/advanced/detector. 21:11:37 発信元: 131 A.HUGUES : do i have to change parameters if i want to calculate dose of secondary photon induced by a neutron source ? and if i want to calculate dose of secondary photon induced by an electron source (bremstrhalung effect) ? 21:13:16 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "do i have to change ..." IF you want high spatial resolution or want to trace reaction physics, you need e-mode = 2 for neutrons and negs = 1 for photons. Otherwise, default by kerma factor is enough because it is fairly accurate and anyway fast. 21:13:16 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "[ Forced Collisions ..." This should be fine. [ C e l l ] 1 1 -c91 -1 2 -3 $ Target 2 2 -c93 -1 3 -4 3 0 -99 #1 #2 10 -1 99 21:13:46 発信元: 131 A.HUGUES : Reacted to "IF you want high spa..." with 🙏 21:14:22 発信元: 098 Maria Sultana : Replying to "[ Forced Collisions ..." @*Tutor* SHIGYO Nobuhiro [ Forced Collisions ] part = proton reg fcl 1 -1.0 as my projectile is deuteron, should I change proton to deuteron? 21:14:33 発信元: *Tutor* Yuho Hirata : Replying to "Does e-mode =2 repla..." When e-mode off, no protons are produced. Therefore, PHITS treats the neutrons as depositing energy. When e-mode=2, the neutron energy is transferred to secondary particles (protons), so the neutrons themselves do not deposit energy 21:16:33 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "[ Forced Collisions ..." Of course, change it to deuteron. 21:21:25 発信元: 205 Zahidul Islam Apu Main : Replying to "Does e-mode =2 repla..." Okay now I get it. Thank you. 21:21:56 発信元: *Tutor* Yuho Hirata : Reacted to "Okay now I get it. T..." with 👍 21:22:16 発信元: khelifi : Reacted to "GEB? Do you mean G..." with 👍 21:22:37 発信元: khelifi : Replying to "for spectrum simul..." thnks 21:24:24 発信元: 098 Maria Sultana : Replying to "[ Forced Collisions ..." @*Tutor* SHIGYO Nobuhiro [T-Product] title = [t-product] in region mesh mesh = reg reg = 1 e-type = 3 # e-mesh is log given by emin, emax and ne ne = 300 # number of e-mesh points emin=1 emax=2000 a-type = -1 na = 7 0.0 2.5 22.5 27.5 82.5 87.5 157.5 162.5 unit = 23 axis = eng file = ddx_product.out output = nuclear epsout = 1 part = neutron factor = 1/c50*1.0E27 y-txt = d^2\sigma/dE/d\Omega (mb/sr/MeV) should I change anything in [T-product] tally? As I worked with unit 1 earlier to my other calculation 21:24:35 発信元: *Tutor* Tatsuhiko Ogawa : One reminder to participants. When you would like to dicsuss with us in a breakout session, please "raise hand". 21:25:35 発信元: 221 David Ferreira : Do you have a library of materials? How do you type steel in Material section? 21:25:57 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "[ Forced Collisions ..." To score neutrons from 2 targets reg = 1 -> reg = 1 2 21:26:42 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "Do you have a librar..." Please see \phits\data\material.inp You can find material definition samples. 21:27:13 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "One reminder to part..." you can request a breakout session whenever you want. 21:27:32 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "[ Forced Collisions ..." Since the energy of deuterons is probably not very high, it would be better to reduce emax. 21:27:53 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "Please see \phits\d..." with 👍 21:27:58 発信元: *Tutor* Tatsuhiko Ogawa : Removed a 👍 reaction from "Please see \phits\d..." 21:28:03 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "Since the energy of ..." with 👍 21:28:09 発信元: 119 A.F.Jundi Main : Replying to "Do you have a librar..." You can also try using material compendium report from PNNL (Compendium of Material Composition Data for Radiation Transport Modeling | Report | PNNL) 21:28:50 発信元: 221 David Ferreira : Reacted to "Please see \phits\d..." with 👍🏻 21:28:53 発信元: *Tutor* Yuho Hirata : Reacted to "You can also try usi..." with 👍 21:28:56 発信元: 221 David Ferreira : Reacted to "You can also try usi..." with 👍🏻 21:29:33 発信元: 079 Sujoy Chatterjee : Reacted to "You can also try u..." with 👍 21:31:04 発信元: 098 Maria Sultana : Replying to "[ Forced Collisions ..." @*Tutor* SHIGYO Nobuhiro and the unit should remain 23? 21:31:43 発信元: 132_K.Skerratt-Love_Main : why does ne = 3? 21:31:49 発信元: 132_K.Skerratt-Love_Main : what does it refer to? 21:32:12 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "[ Forced Collisions ..." You must also correct "factor = 1/c50*1.0E27" because this is a conversion factor for Al. For C target result, you have to multiply the result by c50/(c93/c94*6.02e23*c31) 21:33:02 発信元: *Tutor* SHIGYO Nobuhiro : It represents the number of energy bins. There are three bins: [0,50], [50,100], and [100,200]. 21:33:18 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "You can also try usi..." with 👍 21:33:52 発信元: 132_K.Skerratt-Love_Main : Reacted to "It represents the ..." with 👍 21:35:27 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "[ Forced Collisions ..." The unit of unit = 23 is [1/MeV/sr/source]. As Ogawa san mentioned, you can get DDX by multiplying this by a factor. 21:36:01 発信元: 119 A.F.Jundi Main : Why is the last energy does not have intensity input? 21:36:58 発信元: 098 Maria Sultana : Replying to "[ Forced Collisions ..." @*Tutor* Tatsuhiko Ogawa so I should write "factor = c50/(c93/c94*6.02e23*c31)" like this? 21:37:14 発信元: *Tutor* SHIGYO Nobuhiro : Replying to "Why is the last ener..." Since there are three energy bins, three intensities should be enough. 21:37:28 発信元: *Tutor* Tatsuhiko Ogawa : and 200 is the upper bound of the 3rd energy bin. 21:39:49 発信元: 066 Serge Takoukam : I want to know if there is a relationship between energy and intensities 1, 2 and 3. Are the intensities chosen arbitrarily? 21:40:22 発信元: 119 A.F.Jundi Main : Replying to "Why is the last ener..." What if i shifted the intensity ratio so that the lower bound does (0.0) not have the intensity while the upper bound (200.0) have the intensity? Will it change the interpretation? 21:41:09 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "I want to know if th..." Energies and intensities can be defined independently by users. 21:41:48 発信元: 119 A.F.Jundi Main : Reacted to "Energies and intensi..." with 👍 21:42:10 発信元: 066 Serge Takoukam : Replying to "I want to know if th..." Thanks @*Tutor* Tatsuhiko Ogawa 21:42:35 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "Why is the last ener..." If you mean ne = 3 0.0 0 50.0 3 100.0 2 200.0 1 PHITS will crash because you define 4 intensities despite ne = 3. 21:42:36 発信元: *Tutor* Yuho Hirata : If you are using macOS or Linux and sourceA.inp does not work, please check the [parameter] section for the file(1) setting. If a line such as file(1)=c:/phits is present, please delete it. 21:43:42 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "Why is the last ener..." If you wish to use source beyond 200 MeV, it will be ne = 3 50.0 3 100.0 2 200.0 1 300.0 21:44:36 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "Why is the last ener..." for example. 300 MeV is a value that I randomly chose. 21:44:42 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "Thanks @*Tutor* Tats..." with 👍 21:45:05 発信元: *Tutor* Takuya Furuta : Replying to "If you are using mac..." Hirata-san. It should not be for most cases. PHITS installer automatically detect file(1) path and change it. So it is a person specific problem. 21:45:13 発信元: 119 A.F.Jundi Main : Replying to "Why is the last ener..." @*Tutor* Tatsuhiko Ogawa Ooh okay i understand now. Thank you for the examples! 21:46:57 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "@*Tutor* Tatsuhiko O..." with 👍 21:47:28 発信元: 079 Sujoy Chatterjee : How to use the source type when we want to simulate two different particles, say electron and photon with respective intensities 21:48:15 発信元: 160 Benny Lai : Instead of individually inputting photon energy and intensities of nuclides, does PHITS have a function or a library where I can reference the source nuclide and the code will automatically simulate the discrete energies and intensities that are emitted from the nucleus? 21:48:53 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "How to use the sourc..." "Multi source" function does such a job. Please wait for some pages. 21:49:09 発信元: *Tutor* Takuya Furuta : Replying to "Instead of individua..." It is called RI source option, you will see this in this lecture soon. 21:49:20 発信元: 079 Sujoy Chatterjee : Replying to "How to use the sou..." Okay 21:49:29 発信元: 160 Benny Lai : Reacted to "It is called RI sour..." with 👍 21:49:45 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "Okay" with 👍 21:52:46 発信元: 164 Patrick Davey : Reacted to "It is called RI sour..." with 👍 21:52:48 発信元: 164 Patrick Davey : Reacted to "Instead of individua..." with 👍 21:52:53 発信元: 098 Maria Sultana : Replying to "[ Forced Collisions ..." please confirm me. Is it okay? 21:54:58 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "[ Forced Collisions ..." Not quite. As conversion factors are different for target 1 and 2, you cannot use one global factor. I would suggest to use 2 T-product tallies 21:56:15 発信元: 119 A.F.Jundi Main : Does the determination of intensity only help with the visualization of the graphs? Or it has another physical/mathematical meaning? 21:56:42 発信元: 098 Maria Sultana : Replying to "[ Forced Collisions ..." @*Tutor* Tatsuhiko Ogawa how to write this? can you give me a demo please? My first material is Be and second material is C. 21:57:19 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "[ Forced Collisions ..." [T-Product] $ 1st tally title = [t-product] in region mesh mesh = reg reg = 1 e-type = 3 # e-mesh is log given by emin, emax and ne ne = 300 # number of e-mesh points emin=1 emax=2000 a-type = -1 na = 7 0.0 2.5 22.5 27.5 82.5 87.5 157.5 162.5 unit = 23 axis = eng file = ddx_product1.out output = nuclear epsout = 1 part = neutron factor = 1/c50*1.0E27 y-txt = d^2\sigma/dE/d\Omega (mb/sr/MeV) 21:57:31 発信元: *Tutor* Tatsuhiko Ogawa : [T-Product] $ 2nd tally title = [t-product] in region mesh mesh = reg reg = 2 e-type = 3 # e-mesh is log given by emin, emax and ne ne = 300 # number of e-mesh points emin=1 emax=2000 a-type = -1 na = 7 0.0 2.5 22.5 27.5 82.5 87.5 157.5 162.5 unit = 23 axis = eng file = ddx_product2.out output = nuclear epsout = 1 part = neutron set: c51[c93/c94*6.02e23*c31] $ target density particles/cm^2 factor = 1/c51*1.0E27 y-txt = d^2\sigma/dE/d\Omega (mb/sr/MeV) 22:00:04 発信元: *Tutor* Tatsuhiko Ogawa : Intensity scales up all the results (dose, particle flux, etc.). Is it what you are asking? 22:00:34 発信元: 110 A. Yuniarto : If I want to simulate efficiency calibration of Ge detector using mixed point source, how to interpret the PHITS output? 22:00:59 発信元: 221 David Ferreira : On exercise 4 z0 =-10 and z1=40. Is that a point source? 22:01:13 発信元: 066 Serge Takoukam : Reacted to On exercise 4 z0 =-1... with "👍" 22:01:35 発信元: Elizabeth Musacchio Gonzalez : Replying to "On exercise 4 z0 =-1..." I don’t think so, should be a cylindrical source 22:01:55 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "If I want to simulat..." Please try \phits\lecture\advanced\detector This is a lecture includes detector simulation result comprehension. 22:02:17 発信元: 221 David Ferreira : Replying to "On exercise 4 z0 =-1..." But on the exercise says point source 22:02:33 発信元: 119 A.F.Jundi Main : Replying to "Does the determinati..." @*Tutor* Tatsuhiko Ogawa Yes. I tried the exercise 5 with same intensity ratio of 1:1, but the height of the peak is still different 22:03:00 発信元: Elizabeth Musacchio Gonzalez : Replying to "On exercise 4 z0 =-1..." To have point source r, z0 and z1 should be 0, z0=z1 22:03:36 発信元: 119 A.F.Jundi Main : Replying to "Does the determinati..." oh sorry, not the peak, but the figure in track_xz 22:03:55 発信元: 164 Patrick Davey : Reacted to "Please try \phits\le..." with 👍 22:03:57 発信元: 119 A.F.Jundi Main : Replying to "Does the determinati..." the left one is still brighter than the right one 22:04:21 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "Does the determinati..." They must be the same in theory. They might differ due to statistical uncertainty. 22:04:27 発信元: 221 David Ferreira : Reacted to "To have point source..." with 👍🏻 22:05:00 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "To have point source..." with 👍 22:05:18 発信元: 110 A. Yuniarto : Reacted to "Please try \phits\le..." with 👍 22:05:31 発信元: Elizabeth Musacchio Gonzalez : Replying to "On exercise 4 z0 =-1..." To have isotropic then change dir = all 22:07:49 発信元: 221 David Ferreira : Reacted to "To have isotropic th..." with 👍🏻 22:08:02 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "On exercise 4 z0 =-1..." @221 David Ferreira Do you by chance mix up exercise 4 and 5? Exercise 4 does not mention " z1=40". 22:08:16 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "On exercise 4 z0 =-1..." z1=40 appears in exercise 5. 22:12:02 発信元: 221 David Ferreira : No. For some reason it shows up a point source when I run for exercise 4, and I thought I had z1=40. So when I did the exercise 5 and run and get the same point source I got confused. I simulate again with z0=-10 and z1=40 and got a cylindrical source. So I changed z1=-10 on exercise 4 and I didn't realize 😅 22:12:50 発信元: 110 A. Yuniarto : Replying to "If I want to simulat..." Do you have some references for a marinelli beaker source simulation? 22:14:05 発信元: 119 A.F.Jundi Main : Replying to "Does the determinati..." Ah i see. Thank you! 22:15:37 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "If I want to simulat..." Unfortunately not. If the authors of https://www.sciencedirect.com/science/article/pii/S0969806X25004463 are generous, they will give you input files. 22:15:46 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "Ah i see. Thank you!" with 👍 22:16:23 発信元: 110 A. Yuniarto : Replying to "If I want to simulat..." thank you for your information 22:16:28 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "No. For some reason ..." with 👍 22:16:39 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "thank you for your i..." with 👍 22:24:30 発信元: 005 A MAITRALLAIN : Hi, It seems we still need to define a "proj". Does it matter if we use something other than "photon" in this case? As the 60Co will only emit photons anyways? 22:25:30 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "Hi, It seems we stil..." You can use other "proj" For example, proj = electron means to simulate beta-rays from Co-60 22:25:47 発信元: 005 A MAITRALLAIN : Reacted to "You can use other "p..." with 👍 22:25:52 発信元: 005 A MAITRALLAIN : Replying to "Hi, It seems we stil..." Thank you 22:26:43 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "Hi, It seems we stil..." proj = all is much cooler because all possible source particles from the nucleus are emitted, alpha, beta, beta+, Auger electron, gamma, X-ray. 22:26:49 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "Thank you" with 👍 22:27:22 発信元: 005 A MAITRALLAIN : Reacted to "proj = all is much ..." with 👍 22:27:33 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "proj = all is much ..." with 👍 22:27:40 発信元: *Tutor* Tatsuhiko Ogawa : Removed a 👍 reaction from "proj = all is much ..." 22:27:52 発信元: 164 Patrick Davey : Reacted to "proj = all is much ..." with 👍 22:27:57 発信元: 110 A. Yuniarto : Sorry if I missed the explanation on this: what is the difference between ne and ni? 22:28:11 発信元: 119 A.F.Jundi Main : Replying to "Hi, It seems we stil..." @*Tutor* Tatsuhiko Ogawa will using proj=all cause a longer simulation time? 22:28:46 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "Sorry if I missed th..." ne is the number of Energy bins or Energy points. ni is the number of RIs used only for e-type = 28 and 29 22:29:18 発信元: 110 A. Yuniarto : Reacted to "ne is the number of ..." with 👍 22:29:24 発信元: 110 A. Yuniarto : Replying to "Sorry if I missed th..." thank you 22:29:36 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "thank you" with 👍 22:31:33 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "Hi, It seems we stil..." Generally yes. "proj" must be chosen what source particles you wish to calculate. If you are not sure of it, proj = all might be overkill but a safe option. 22:32:20 発信元: 119 A.F.Jundi Main : Reacted to "Generally yes. "pro..." with 👍 22:37:35 発信元: 185.W Asztalos : For a RI source, what is the effect of changing ``proj''? Is it a filter for gamma vs beta vs alpha, or redundant? That is, since we have Cs137, by setting proj = photon, do we ignore the beta particle that normally is emitted as well? 22:38:32 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "For a RI source, wha..." it works as a filter. "by setting proj = photon, do we ignore the beta particle" Totally correct. 22:39:15 発信元: 185.W Asztalos : Replying to "For a RI source, wha..." Thank you! Looks like if I had waited a few minutes I would've gotten my answer, sorry! 22:39:22 発信元: *Tutor* Tatsuhiko Ogawa : Reacted to "Thank you! Looks lik..." with 👍 22:41:00 発信元: 132_K.Skerratt-Love_Main : I am unable to attend tomorrow. I will watch the recording back, but if i have any questions, can email for help, or should I pot to forum? 22:41:11 発信元: 132_K.Skerratt-Love_Main : *post to forum 22:41:32 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "I am unable to atten..." Please use PHITS forum. It is available 24/7. 22:41:40 発信元: 132_K.Skerratt-Love_Main : Reacted to "Please use PHITS f..." with 👍 22:41:43 発信元: *Tutor* Tatsuhiko Ogawa : Replying to "I am unable to atten..." https://meteor.nucl.kyushu-u.ac.jp/phitsforum/ 22:49:08 発信元: 147 Audric Zachary Uy : May I ask if it’s possible to simulate a moving source? 22:49:34 発信元: *Tutor* Tatsuhiko Ogawa : http://phits.jaea.go.jp/tutorial/rec.html https://forms.gle/bQjAzdJjXEiNbJhr6 22:49:50 発信元: 176 K.Kim : Thank you. See you tomorrow 22:49:54 発信元: 042 Melih Tahsin Değer Main : Thanks for today, have a good day... 22:49:56 発信元: 216 Chadchayada Samarak : Thank you 22:49:56 発信元: 119 A.F.Jundi Main : Thank you so much sensei! 22:49:58 発信元: 219 Ismail Hossain : Thank you all 22:50:00 発信元: 124 S.Tancharakorn : Thank you very much, bye bye. 22:50:02 発信元: 029 Herlan Setiawan : thank you 22:50:05 発信元: 034 J.Suello Main : Thank you 22:50:05 発信元: 131 A.HUGUES : Thank you sensei!! 22:50:09 発信元: 207 Poonamjot : Thankyou so much , bye bye 22:50:11 発信元: 048 nur aida : arigao mina sensei 22:50:14 発信元: 064 Youngran Kwon : Thanks, See you tomorrow. 22:50:16 発信元: 032 I.Gordeev main : Thank you very much, bye 22:50:21 発信元: Iaroslav Bielov : Thank you, have a good day! 22:50:22 発信元: 092 Bal Vikram Khatri : Thank you very much 22:50:22 発信元: 022 D. R. Upadhyay : Thank you so much 22:50:30 発信元: 066 Serge Takoukam : Thanks 22:50:32 発信元: 068 ahid nurmanjaya : お疲れ様です