20:01:03 発信元: *Tutor* Tatsuhiko OGAWA : 1.Video/Chat files of lectures and guidance are available from here http://phits.jaea.go.jp/tutorial/rec.html 2.If you are Mac or Linux user and take the optional course “Medical data treatment by RT-PHITS” on Friday, install Python and its libraries following https://phits.jaea.go.jp/tutorial/Install_Python.pdf 3.Please answer a survey by Friday from here https://forms.gle/htPwzn7eSbFMKoAr5 20:18:41 発信元: 140.Mahdi Bakhtiari : Hello, phits tutors. I have a question (not related to this lecture though!) but about the weight window. Can we use it for very low photon energy in the range of ~1keV to 200 keV safely? 20:20:42 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "Hello, phits tutors...." Yes. Absolutely no problems. 20:20:48 発信元: *Tutor* OHNISHI Seiki : Replying to "Hello, phits tutors...." No problem. However, low energy photons make little contribution to tallies, and too many low energy photons make the calculation efficiency worse. 20:21:24 発信元: 140.Mahdi Bakhtiari : Reacted to No problem. However,... with "👍" 20:21:26 発信元: 144 Margarita Tzivaki : is it possible to use the dump file as a backup of all interactions for long simulations? i.e. produce one every couple thousand histories in case the calculation fails somewhere along the way? Or is there an other method for that? 20:21:27 発信元: 140.Mahdi Bakhtiari : Reacted to Yes. Absolutely no p... with "👍" 20:21:31 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "No problem. However,..." with 👍 20:21:42 発信元: 140.Mahdi Bakhtiari : Removed a 👍 reaction from "No problem. However,..." 20:21:47 発信元: 140.Mahdi Bakhtiari : Reacted to No problem. However,... with "👍" 20:22:33 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "is it possible to us..." I agree it is one of the reasonable ways to use dump. 20:23:30 発信元: 144 Margarita Tzivaki : Hat auf "I agree it is one ..." mit 👍 reagiert 20:24:00 発信元: 140.Mahdi Bakhtiari : And is it recommended to use egs in that case? Or the phits itself? I mean for low energy photons. 20:24:48 発信元: 144 Margarita Tzivaki : Antworten an "is it possible to ..." so in that case the dump region would be the whole universe? 20:25:01 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "is it possible to us..." To use as a backup, however, I would rather use standard tally outputs. PHITS can use the energy distribution of tally outputs as the source Energy distribution in the second step. This method is advantageous in the sense that you Don't have to make heavy dump files. 20:27:00 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "is it possible to us..." @144 Margarita Tzivaki >dump region would be the whole universe Not quite because such dump double counts particles. You should dump particles somewhere between the source and tallies. 20:27:30 発信元: *Tutor* OHNISHI Seiki : Replying to "Hello, phits tutors...." The use of EGS is almost always recommended. 20:27:52 発信元: 140.Mahdi Bakhtiari : Reacted to The use of EGS is al... with "👍" 20:28:08 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "Hello, phits tutors...." For low Energy photons, EGS is recommended regardless of variance reduction. Variance reduction algorithm is built to be independent from reaction phosics algorithms like EGS. 20:28:16 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "The use of EGS is al..." with 👍 20:28:33 発信元: 144 Margarita Tzivaki : Antworten an "is it possible to ..." If I use the tallies as a backup I may loose interactions between histories though, no? 20:28:42 発信元: 140.Mahdi Bakhtiari : Reacted to For low Energy photo... with "👍" 20:29:24 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "is it possible to us..." Exactly. So it is matter of compromise between particle corrélations and data size. 20:29:28 発信元: 144 Margarita Tzivaki : Hat auf "If I use the talli..." mit 👍 reagiert 20:30:09 発信元: 144 Margarita Tzivaki : Antworten an "is it possible to ..." Yes, ok that makes sense 20:30:14 発信元: 144 Margarita Tzivaki : Hat auf "Exactly. So it is ..." mit 👍 reagiert 20:30:44 発信元: *Tutor* Hunter Ratliff : Replying to "is it possible to us..." Another thing to keep in mind is that the dump file is written to continuously while the standard tally outputs are only updated at the end of each batch. 20:31:24 発信元: 217 Sukarna Banik : Hello, I want to know how can I design a input file for radioisotope production for example, 89-Zr from the output file of DDNY. how can I put the output data as input data? I'm sharing a demo output file 20:31:24 発信元: 144 Margarita Tzivaki : Antworten an "is it possible to ..." I don't think I completely understand why the dump file double counts particles. Could you please elaborate? 20:31:48 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "is it possible to us..." Actually, PHITS has an option for [parameters] section dumpall=1 This requests to dump all the details of particle transport to an external file. If you later run the same input with icntl=12 PHITS reads the dumpall file. This option is probably the closest to the answer to your initial question. 20:32:06 発信元: 144 Margarita Tzivaki : Hat auf "Another thing to k..." mit 👍 reagiert 20:32:48 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "Hello, I want to kno..." Excuse me but could you ask this question in a breakout session? 20:33:13 発信元: 144 Margarita Tzivaki : Antworten an "is it possible to ..." ooooh I see, yes this is what I was wondering about. So like that I can continuously write a dump file and then automatically delete it if the simulation runs smoothly 20:33:18 発信元: 208 S.Laca : Hello, which tally or what should I use if I want to see how many particles reach some geometry depending of the angle I choose? 20:33:28 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "Hello, I want to kno..." It is not easy to answer your question in detail in chat. 20:33:49 発信元: *Tutor* Hunter Ratliff : Replying to "is it possible to us..." If you had a dump file for a [T-Cross] tally with multiple surfaces a single particle could cross (like a neutron traveling through an xyz mesh), you'll write that neutron to the dump file every time it crosses from one voxel into another. So, you're counting it many times in that circumstance. 20:33:52 発信元: 217 Sukarna Banik : Replying to "Hello, I want to kno..." I'm actually at a conference that's why can't connect through audio 20:34:03 発信元: 060 Dewa Ngurah Yudhi Prasada : 2nd t-cross doesn't need dump parameters? 20:34:13 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "Hello, which tally o..." To calculate angular dépendent flux, T-cross with a-type (angular binning) is the best. 20:34:56 発信元: 208 S.Laca : Respondiendo a "Hello, which tally..." perfect thank you 20:34:57 発信元: *Tutor* Hunter Ratliff : Replying to "is it possible to us..." If you want to use a dump file as a source, you should be careful to avoid such situations (unless it's actually desired behavior for your application). 20:35:13 発信元: 264 + Kebede Shogile : What are the numbers 1 101 101 314.15? 20:35:23 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "2nd t-cross doesn't ..." Right. The second t-cross is intended to calculate the normalization factor between 1st and 2nd régions. 20:35:35 発信元: 140.Mahdi Bakhtiari : Replying to "is it possible to us..." Actually, PHITS has an option for [parameters] section dumpall=1 This requests to dump all the details of particle transport to an external file. If you later run the same input with icntl=12 PHITS reads the dumpall file. This option is probably the closest to the answer to your initial question. What is the format of that file? Can we check the particle transportation for every history number? 20:35:55 発信元: *Tutor* Hunter Ratliff : Replying to "is it possible to us..." But yes, dumpall also exists and is this kind of comprehensive "runtape" style file, but be warned that it gets huge very quickly, as it is truly writing out everything per history. 20:36:16 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "Hello, I want to kno..." Your microphone or speaker are not working? 20:36:44 発信元: *Tutor* Hunter Ratliff : Replying to "is it possible to us..." The manual's section 5.2.26 details the format of the dumpall file 20:36:56 発信元: 217 Sukarna Banik : Replying to "Hello, I want to kno..." What I mean was I can't talk right now 20:36:57 発信元: 144 Margarita Tzivaki : Antworten an "is it possible to ..." ah ok, that also explains the void in the dump section. So dumping the whole geometry is a terrible idea for using it as a source, but may be a reasonable strategy if I want to use it as a backup (accepting that I may get angry calls from the cluster admins because of file size) 20:37:20 発信元: 140.Mahdi Bakhtiari : Reacted to The manual's section... with "👍" 20:37:27 発信元: 144 Margarita Tzivaki : Hat auf "The manual's secti..." mit 👍 reagiert 20:37:30 発信元: *Tutor* Hunter Ratliff : Replying to "is it possible to us..." The dumpall file is, admittedly, quite difficult to parse yourself (as a human) given the huge variety of things it logs. 20:37:37 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "perfect thank you" with 👍 20:37:57 発信元: 217 Sukarna Banik : Replying to "Hello, I want to kno..." Do you have any sample data for this type work in PHITS? I need a sample input file 20:38:07 発信元: 144 Margarita Tzivaki : Antworten an "is it possible to ..." Yes the runtpe functionality is what I was after, thank you! 20:38:26 発信元: 144 Margarita Tzivaki : Hat auf "If you had a dump ..." mit 👍 reagiert 20:38:48 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "What are the numbers..." 1 : Sequencial ID number 101 : cell number of the upstream cell 102 : cell number of the downstream cell 314.15: Area of this tally 20:39:27 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "Yes the runtpe funct..." with 👍 20:39:52 発信元: *Tutor* Hunter Ratliff : Replying to "is it possible to us..." One thing to note is that, unlike MCNP, PHITS can resume a calculation from where it left off without a runtape file. It will just load the random seed (in binary) from batch.out and load in the tally results and statistics from their file. 20:40:10 発信元: *Tutor* Hunter Ratliff : Replying to "is it possible to us..." *PHITS will pick up from the last completed batch 20:40:17 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "One thing to note is..." with 👍 20:40:32 発信元: *Tutor* Hunter Ratliff : Replying to "is it possible to us..." "Restart mode" is accessed with istdev < 0 20:41:25 発信元: 144 Margarita Tzivaki : Antworten an "is it possible to ..." OH! right! I had forgotten about the batches. Well, then the question is somewhat obsolete. This is excellent, thank you 20:41:44 発信元: *Tutor* Hunter Ratliff : Replying to "is it possible to us..." I don't exactly remember the behavior if you had dump files being written to. If you crash a PHITS simulation, the dump files will of course have info up to the very last history simulated pre-crash but the tally outputs will only go up to the last history of the previously completed batch. 20:42:24 発信元: *Tutor* Tatsuhiko OGAWA : @140.Mahdi Bakhtiari "Dumpall = 1" dumps binary file which you can never read. "Dumpall = -1" dumps ascii file which you can harldy read. 20:42:39 発信元: *Tutor* Hunter Ratliff : Reacted to "@140.Mahdi Bakhtiari..." with 😆 20:42:45 発信元: 144 Margarita Tzivaki : Hat auf "@140.Mahdi Bakhtia..." mit 😂 reagiert 20:43:08 発信元: 028 Amiegbereta Edwin Ehis : Good day. please am just joining, whats the subject of discuss and which folder can I find it? Thanks 20:44:00 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "Good day. please am ..." We are at \phits\lecture\advanced\sourceB to learn to use dump files as sources. 20:44:14 発信元: William Asztalos : Is there any reason we set dump = -12 instead of dump = -11 like in the earlier slides? not sure if it matters or is just a logical unit for fortran 20:44:17 発信元: *Tutor* Hunter Ratliff : Replying to "is it possible to us..." I'm not sure if there's some graceful way PHITS can only resume writing to the dump file once it "caught back up" to where it left off (presumably only possible if you were writing batch and history numbers to the dump), but the safest way would probably be to just rename the existing dump file, have the restart calculation make new dump file, and just merge their contents in post processing (with your own filters to prevent "double-counting" of the histories of the last batch that had failed to complete). 20:44:19 発信元: 171 N. Stoffle : What is the difference between -11 and -12 for the dump parameter? 20:44:29 発信元: 100 Herlan Setiawan : Why in [t-cross] : output = flux, instead using output = current 20:44:31 発信元: 214 S. Schmidt : I would like to know, how can I set the y axis of the cross figure to linear scale? 20:44:42 発信元: 144 Margarita Tzivaki : Hat auf "I'm not sure if th..." mit 😮 reagiert 20:45:12 発信元: 028 Amiegbereta Edwin Ehis : Reacted to "We are at \phits\le..." with 👍 20:45:19 発信元: *Tutor* Hunter Ratliff : Replying to "is it possible to us..." (This is only applicable if you're actually trying to use the dump file in your own analyses, not for just using it to re-insert into PHITS as a source.) 20:45:24 発信元: 144 Margarita Tzivaki : Antworten an "is it possible to ..." yeah that is totally fine, I can work with strategically selecting batch numbers in a way that I will not lose too much information, that's actually way more convenient than the runtpe. So the dump file here works more like the source write and read function in mcnp 20:45:32 発信元: 214 S. Schmidt : It is the number of physical quantities being dumped (and stated in the row below) 20:45:42 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "It stands for the nu..." with 👍 20:45:48 発信元: *Tutor* Hunter Ratliff : Replying to "is it possible to us..." Exactly 20:46:05 発信元: 171 N. Stoffle : Reacted to "It is the number of ..." with 👍 20:46:08 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "What is the differen..." Yes. Stefan is right. Therefore in this file it must be dump = -12 1 2 3 4 5 6 7 8 9 10 18 19 20:46:27 発信元: 144 Margarita Tzivaki : Antworten an "is it possible to ..." Sorry for all the mcnp questions :P 20:46:54 発信元: *Tutor* Hunter Ratliff : Replying to "is it possible to us..." I'd learned MCNP before PHITS too, so I understand where you're coming from here haha. 20:46:55 発信元: 144 Margarita Tzivaki : Hat auf "(This is only appl..." mit 👍 reagiert 20:47:02 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "Is there any reason ..." Stefan just answered. 'It is the number of physical quantities being dumped (and stated in the row below)' 20:47:08 発信元: 264 + Kebede Shogile : Replying to "What are the numbers..." 🙏 20:47:21 発信元: 144 Margarita Tzivaki : Hat auf "I'd learned MCNP b..." mit 😂 reagiert 20:47:39 発信元: 144 Margarita Tzivaki : Antworten an "is it possible to ..." phits is currently winning this contest I think :D 20:48:06 発信元: 052 Dali Georgobiani : Replying to "is it possible to us..." What about FLUKA? 🙂 20:48:20 発信元: 217 Sukarna Banik : Replying to "Hello, I want to kno..." Hello! 20:48:39 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "I would like to know..." Please insert angel =p :ylin to the tally 20:49:09 発信元: 144 Margarita Tzivaki : Antworten an "is it possible to ..." haven't gotten to that one yet! can you recommend it?? 20:49:19 発信元: 052 Dali Georgobiani : Replying to "is it possible to us..." Very much so 🙂 20:49:46 発信元: 144 Margarita Tzivaki : Hat auf "Very much so 🙂" mit ❤️ reagiert 20:49:57 発信元: 052 Dali Georgobiani : Replying to "is it possible to us..." Don’t get me wrong, I absolutely love PHITS! It would be just another tool in your pocket 🙂 20:50:00 発信元: *Tutor* Hunter Ratliff : Replying to "is it possible to us..." I am a biased source, but unless MCNP has made PTRAC far more friendly than it was ~8 years ago, PHITS is definitely winning on this event-by-event dumping front. 20:50:17 発信元: William Asztalos : Replying to "Is there any reason ..." Thank you, I did end up seing that someone asked after me and got the answer, but I appreciate that you went back to my question just in case! 20:50:45 発信元: 052 Dali Georgobiani : Reacted to "I am a biased source..." with ❤️ 20:50:51 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "Why in [t-cross] : o..." You are probably right. In theory, this tally must score the number of particles regardless of the angle. If the particles enter the surface normally, this does not impact the result a lot but I will discuss this issue with my colleagues to make sure if we need to correct this material. 20:50:55 発信元: 144 Margarita Tzivaki : Antworten an "is it possible to ..." sounds like the set-up of a code-off 🥳 10/10 would watch 20:51:03 発信元: 115 J. Santos : where are you seeing the relative error of the spectrum? 20:51:07 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "Thank you, I did end..." with 👍 20:51:50 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "Is there any reason ..." The relative errors are written in cross_eng.out 20:51:57 発信元: 144 Margarita Tzivaki : Antworten an "is it possible to ..." I'm not too deep into the ptrac file, but it's still pretty terrible 20:51:57 発信元: 214 S. Schmidt : Replying to "I would like to know..." With this command "angel =p :ylin", I get the following error: ***** Error Message in ANGEL ***** cross_eng.out 74 : ERROR = P: Parameter Description is Wrong in [ p :y ] 20:52:05 発信元: *Tutor* Hunter Ratliff : Reacted to "I'm not too deep int..." with 😆 20:52:08 発信元: *Tutor* Hunter Ratliff : Replying to "is it possible to us..." And combined with [Counter], the dump files in PHITS are really powerful for event-by-event outputs, enabling you to make your own custom tallies/binning in post processing, without resorting to writing your own [T-Userdefined] Fortran tally and recompiling PHITS. 20:52:10 発信元: 140 Mahdi Bakhtiari Main : I am a biased source, but unless MCNP has made PTRAC far more friendly than it was ~8 years ago, PHITS is definitely winning on this event-by-event dumping front. I have also tried to use something like DUMP option in FLUKA, but PHITS is far better customized for this purpose 🙂 20:52:11 発信元: 115 J. Santos : Reacted to "The relative errors ..." with 👍 20:52:33 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "I would like to know..." Sorry. What if angel =p:ylin without extra space 20:52:51 発信元: *Tutor* Hunter Ratliff : Replying to "is it possible to us..." I'll also add that PHITS Tools has full support for parsing dump files, both the binary and text formats. It does not have support for reading dumpall output though, as that's a whole separate beast. 20:53:20 発信元: 144 Margarita Tzivaki : Antworten an "is it possible to ..." only thing I have ever tried to do with it, is make it into a source term for srim and I quickly abandoned that idea. But maybe I could warm up this idea using phits now that I think about it 20:53:20 発信元: wakai : How to know the note value "Note that the relative errors are 10 -20%)" in page 14? 20:53:23 発信元: *Tutor* Takuya Furuta : Replying to "I would like to know..." Probably p: is not needed. Only angel = ylin. 20:53:45 発信元: 144 Margarita Tzivaki : Hat auf "I'll also add that..." mit ❤️ reagiert 20:53:47 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "Probably p: is not n..." with 👍 20:53:52 発信元: 144 Margarita Tzivaki : Hat auf "Don’t get me wron..." mit 👍 reagiert 20:54:08 発信元: 214 S. Schmidt : Replying to "I would like to know..." Thank you, now its working! 20:54:21 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "How to know the note..." Please see cross_eng.out The statistical uncerainty is written in there 20:54:23 発信元: 129 - K. Butalag : Is it possible to have a second step simulation with higher maxcas and maxbch compared to the first step? 20:54:40 発信元: 144 Margarita Tzivaki : Antworten an "is it possible to ..." I'll also add that PHITS Tools has full support for parsing dump files, both the binary and text formats. It does not have support for reading dumpall output though, as that's a whole separate beast. ok this is giving me all sorts of ideas 20:54:50 発信元: *Tutor* Hunter Ratliff : Back in my PhD work, I actually had given up on PTRAC and instead used a code called MCPL (Monte Carlo Particle Lists) to parse MCNP's SSW output file haha 20:55:36 発信元: 144 Margarita Tzivaki : Antworten an "is it possible to ..." one more code to look into.... 20:55:42 発信元: *Tutor* Hunter Ratliff : But yeah, I've been using PHITS dump files as actual outputs to be studied (not just sources to be reinserted) for a lot of simulation work in my own research work. 20:56:11 発信元: 060 Dewa Ngurah Yudhi Prasada : if we write idmpmode = 0, we can assign new maxcas and maxbch? 20:56:34 発信元: 214 S. Schmidt : Reacted to "Probably p: is not n..." with 👍 20:56:43 発信元: *Tutor* Tatsuhiko OGAWA : Yes. If you set idmpmode=1 and dmpmulti = 2 You can use your dump file for 2 cycles (every particle is used twice). If you set idmpmode=0 you will use dumped source until source particle is read by maxcas*maxbch times. Takuya is right. I have corrected this reply. 20:56:59 発信元: William Asztalos : Do we set totfact = 1.0 now? Or is that no longer needed? 20:57:13 発信元: *Tutor* Takuya Furuta : Replying to "Is it possible to ha..." dmpmulti is a option for idmpmode=1. 20:57:18 発信元: *Tutor* Tatsu Sato : Replying to "Hello, I want to kno..." I am not sure what you exactly would like to calculate, but for estimating isotope production yield, recommendation/neutronsource could be a good starting point for your simulation. 20:57:42 発信元: William Asztalos : Replying to "Do we set totfact = ..." Nevermind! The command line tells me 20:57:50 発信元: wakai : Replying to "How to know the note..." Thank you. 20:58:51 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "dmpmulti is a option..." with 🙏 20:58:58 発信元: 217 Sukarna Banik : Replying to "Hello, I want to kno..." okay 20:59:14 発信元: *Tutor* Takuya Furuta : Replying to "Is it possible to ha..." When idmpmode=1 specified, the number of dump data obtained in 1st run will be used. If you use dmpmulti=2, you used doubling the source. When idmpmode=0, you can specify maxcas*maxbch more than dumped data. In this case, the dumped data will be reused. 20:59:52 発信元: 115 J. Santos : what does dump=0 mean? 20:59:55 発信元: 144 Margarita Tzivaki : Antworten an "is it possible to ..." so following up on phitstools, how complicated is it to get exact impulse and location data on a specific type of particle in a region from a dump file? 21:01:07 発信元: 214 S. Schmidt : You need to set the source section to off 21:01:19 発信元: 140 Mahdi Bakhtiari Main : Replying to "Is it possible to ha..." if we set idmpmode=0, i think we need to take care of the ratio of the generated_particle/primary_particle as a correction factor? or am I missing something? 21:01:24 発信元: *Tutor* Takuya Furuta : recommendation/NuclearReaction may be the one you want. 21:02:00 発信元: *Tutor* Tatsu Sato : Replying to "what does dump=0 mea..." When dump = 0, PHITS reads the information of each particle (e.g. 1 2 3 4 5 6 7 8 9 18 19) from the tally output file. 21:02:05 発信元: *Tutor* Tatsuhiko OGAWA : dump is the number of data entries in the dumped source. dump =0 means to read all data automatically in the file (in this case, 12 entries). 21:02:34 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "You need to set the ..." with 👍 21:04:01 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "Is it possible to ha..." Yes. You need to manually calculate the ratio and use it as totfact 21:04:22 発信元: 140 Mahdi Bakhtiari Main : Reacted to "Yes. You need to man..." with 👍 21:04:24 発信元: 144 Margarita Tzivaki : Hat auf "But yeah, I've bee..." mit 👍 reagiert 21:04:48 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "Thank you." with 👍 21:06:04 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "if we write idmpmode..." Yes. When idmpmode = 0, PHITS samples source particles for maxcas*maxbch times irrespective of the maxcas and maxbch in the first run. 21:06:54 発信元: *Tutor* Hunter Ratliff : Replying to "is it possible to us..." I'm not exactly sure of your use case here, but you have a few different levers to pull: 1) The tally settings to determine what actually gets written to the dump file. (Since the dump file does not have a region number entry, only (x,y,z) coordinates, it may be convenient to separate different regions into different tally dumps, or use a [Counter] to distinguish them, since the counter values can also be written to the dump file). 2) The available columns in the dump file (see manual section 6.7.22). So, presuming you've set up the tally and its dump parameters to give you all the information you need to identify your particles/events of interest (and to not be mixed with unwanted events), this is completely doable. 21:07:05 発信元: 144 Margarita Tzivaki : Hat auf "I am a biased sour..." mit ❤️ reagiert 21:07:14 発信元: 144 Margarita Tzivaki : Antworten an "I am a biased sour..." good to know! 21:09:40 発信元: *Tutor* Hunter Ratliff : Replying to "is it possible to us..." You can write dump files for [T-Product], [T-Cross], and [T-Time] tallies, which should cover most of the cases of interest 21:12:02 発信元: *Tutor* Hunter Ratliff : Reacted to "I am a biased source..." with ❤️ 21:12:10 発信元: 144 Margarita Tzivaki : Antworten an "is it possible to ..." Sweet! That sounds excellent. Assuming phits tools is not as frustrating to get working as mcnp tools, the idea to use charged particle information for material shennanigans in SRIM could maybe be resurrected. 21:12:39 発信元: 144 Margarita Tzivaki : Antworten an "is it possible to ..." thank you again for the super detailled explanations! 21:13:33 発信元: 153 M.Onegin : error in geometry 21:14:01 発信元: 214 S. Schmidt : No geometry definition for cell 103 21:14:02 発信元: 052 Dali Georgobiani : Cell 103 is incomplete 21:14:04 発信元: 056 David Ferreira : Cell 103 21:14:32 発信元: 114 Gordeev Main : Maybe chage paramets in [ T - C r o s s ] reg = 3 21:14:57 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "error in geometry" with 👍 21:15:08 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "No geometry definiti..." with 🙏 21:15:09 発信元: *Tutor* Tatsuhiko OGAWA : Removed a 👍 reaction from "error in geometry" 21:15:11 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "error in geometry" with 🙏 21:15:14 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "Cell 103 is incomple..." with 🙏 21:15:16 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "Cell 103" with 🙏 21:15:19 発信元: *Tutor* Hunter Ratliff : Replying to "is it possible to us..." You're welcome! I've tried to make PHITS Tools as user-friendly as possible, but it is still a young software. If you're experienced with Python, both standard tally outputs and dump files are outputted to Pandas DataFrames (in addition to "more traditional" array formats), so hopefully it'll be straightforward for you to use. 😊 21:15:22 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "Maybe chage paramets..." with 🙏 21:16:08 発信元: 144 Margarita Tzivaki : Hat auf "You're welcome! I'..." mit 👍 reagiert 21:16:20 発信元: 144 Margarita Tzivaki : Antworten an "is it possible to ..." pandas ❤️ 21:16:44 発信元: 264 + Kebede Shogile : Why I am getting such output 21:17:23 発信元: 144 Margarita Tzivaki : Antworten an "is it possible to ..." that's all I need 21:17:53 発信元: 214 S. Schmidt : Replying to "72f279dde5ffb8b82364a0cd1fb2859.png" Did you set the first source to off? 21:18:25 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "Why I am getting suc..." It is a right output up to exercise 3 because the cell 102 is outer void 21:19:03 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "Why I am getting suc..." Outer void is a blackhole in other words. All the particles were absorbed when they reached cell 102. 21:19:22 発信元: *Tutor* Hunter Ratliff : The example on GitHub (the script and notebook have the same contents) hopefully are a good start in just getting a hold of the syntax. It shows example DFs for both standard tally output and a dump file (at the very end). https://github.com/Lindt8/PHITS-Tools/blob/main/example/example_PHITS_tools.ipynb 21:19:40 発信元: 264 + Kebede Shogile : Replying to "72f279dde5ffb8b82364a0cd1fb2859.png" Yes 21:19:52 発信元: 144 Margarita Tzivaki : Hat auf "The example on Git..." mit 👍 reagiert 21:20:15 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "72f279dde5ffb8b82364a0cd1fb2859.png" Your run crashed and EPS file was not updated? 21:20:30 発信元: *Tutor* Tatsu Sato : Replying to "72f279dde5ffb8b82364a0cd1fb2859.png" It seems that the material of cell 102 is outer void (-1) 21:21:51 発信元: 144 Margarita Tzivaki : MCPL : MCPL converters for PHITS 21:22:06 発信元: 052 Dali Georgobiani : Reacted to "MCPL : MCPL converte..." with ❤️ 21:22:22 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "MCPL : MCPL converte..." with ❤️ 21:22:28 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "The example on GitHu..." with 👍 21:29:46 発信元: 052 Dali Georgobiani : Thank you very much, Matsuya-san! 21:30:17 発信元: 043 A.Bekouale : Thanks Mr 21:31:47 発信元: 056 David Ferreira : Session 2 and 3 will be recorded? 21:32:40 発信元: 072 E. Musacchio Gonzalez : Where I can find later the other lessons? 21:32:58 発信元: 150 Melih Tahsin Değer Main : Thanks to tutors tutorial and students for great questions 21:33:30 発信元: *Tutor* Tatsuhiko OGAWA : https://phits.jaea.go.jp/tutorial/rec.html Recording files are and will be here 21:33:39 発信元: 083 Firliyani : Reacted to "https://phits.jaea.g..." with 👍 21:33:50 発信元: 072 E. Musacchio Gonzalez : Thank you! 21:34:13 発信元: 214 S. Schmidt : How long will the videos be available? 21:34:30 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "How long will the vi..." About 4 months 21:34:33 発信元: 222 Theo Victor Cavite : Reacted to "https://phits.jaea.g…" with 👍 21:34:34 発信元: 214 S. Schmidt : Reacted to "About 4 months" with 👍 21:34:49 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "Thank you!" with 👍 21:34:53 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "Thanks to tutors tut..." with 👏 21:35:00 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "Thank you very much,..." with 👍 21:35:03 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "Thanks Mr" with 👍 21:42:40 発信元: 085_Fungki Iqlima Nasyidiah : Sensei, I completely sorry. I can't join to break out room because suddenly i got an assignment that had to be done tonight too 21:44:01 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "Sensei, I completely..." Which session do you go? 21:44:49 発信元: 214 S. Schmidt : Replying to "Sensei, I completely..." I think she just wanted to say that she cannot stay in the afternoon session. But I may be wrong 21:44:51 発信元: 083 Firliyani : Replying to "Sensei, I completely..." i think she wants to leave the zoom 21:44:54 発信元: 060 Dewa Ngurah Yudhi Prasada : Replying to "Sensei, I completely..." @*Tutor* Tatsuhiko OGAWA I think she wants to leave 21:45:42 発信元: 085_Fungki Iqlima Nasyidiah : honestly i want go to room the shielding, but i think i should leave right now. Sorry Sensei 21:45:55 発信元: 085_Fungki Iqlima Nasyidiah : thank you sensei 21:58:02 発信元: 083 Firliyani : How about if we want to replace Ir-192 with F-18 or another RI (such as Co-60) as a radiation source? 21:58:25 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "honestly i want go t..." with 👍 21:58:28 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "thank you sensei" with 👍 22:00:36 発信元: *Tutor* Hunter Ratliff : Replying to "How about if we want..." You would simply replace the text "Ir-192" with "F-18" or the other nuclide you want as the source, making sure you also update it's activity next to it (in Bq). 22:01:32 発信元: 014 ahid nurmanjaya : Reacted to "You would simply rep..." with 👍 22:01:59 発信元: *Tutor* Hunter Ratliff : And if you want more than one source isotope at the same time, you'd change ni=1 to ni=X, where X is the number of nuclides you want in your source and the "ni=X" line is followed by X lines each specifying an isotope and its activity. 22:03:02 発信元: *Tutor* Hunter Ratliff : Replying to "How about if we want..." See \phits\lecture\advanced\sourceA for a review on how to use the RI source section. 22:03:17 発信元: 110 I. Casalbore : Hi, sorry I am joining again right now, can you please tell me how to join the shielding lesson? 22:03:43 発信元: 214 S. Schmidt : Replying to "Hi, sorry I am joini..." \phits\lecture\advanced\shielding 22:03:44 発信元: 083 Firliyani : Replying to "How about if we want..." @*Tutor* Hunter Ratliff Sorry if I missed any information about the "ni" and e-type codes for an RI, is there any connection with JENDL? 22:04:15 発信元: 083 Firliyani : Replying to "How about if we want..." @*Tutor* Hunter Ratliff oke thank you 22:04:24 発信元: *Tutor* Hunter Ratliff : Replying to "Hi, sorry I am joini..." This session is the shielding session, so you're in the right place if you're referring to the breakout rooms. 22:05:00 発信元: 110 I. Casalbore : Replying to "Hi, sorry I am joini..." Yes thank you! 22:05:36 発信元: *Tutor* Hunter Ratliff : There are no specific codes you need to enter. The "ni=" is just "number of isotopes" and it tells PHITS how many isotopes it should be looking for in the provided source section (inputted by you) 22:06:29 発信元: *Tutor* Hunter Ratliff : Replying to "How about if we want..." RI source uses e-type=28 in PHITS 22:07:21 発信元: *Tutor* Hunter Ratliff : Replying to "How about if we want..." I believe the decay data for RI sources is actually coming from ICRP 107 (not JENDL) 22:08:03 発信元: 083 Firliyani : Replying to "How about if we want..." @*Tutor* Hunter Ratliff ok, i will practice more after class. thank you for your information 22:08:10 発信元: *Tutor* Hunter Ratliff : Reacted to "@*Tutor* Hunter Ratl..." with 👍 22:08:17 発信元: 088 G.A. Atregenio : how can we set an energy range for the radiation source? like setting the energy source from 10 keV to 100 GeV? what appropriate energy distribution function can I used to it? 22:09:55 発信元: 060 Dewa Ngurah Yudhi Prasada : how will be the result of t-point if we use negs=1? 22:10:08 発信元: wakai : Which folder is t-point? 22:10:51 発信元: *Tutor* Hunter Ratliff : Replying to "how can we set an en..." For RI source, the energies are pre-defined. But, for a more generic source, you have a lot of different options for how to specify the energy distribution with the various e-type values (see manual section 5.3.19), depending on if you already have the distribution as tabulated data points, as a function, etc. 22:11:05 発信元: *Tutor* Shintaro Hashimoto (sub, chat only) : We should not set negs=1 when calculating t-point. 22:11:50 発信元: 088 G.A. Atregenio : Reacted to "For RI source, the e..." with ❤️ 22:11:57 発信元: *Tutor* Shintaro Hashimoto (sub, chat only) : Replying to "how will be the resu..." We can only estimate neutron and photon doses by the t-point tally. 22:12:07 発信元: *Tutor* Hunter Ratliff : Replying to "Which folder is t-po..." The same input file already has the [T-Point] tally in it, at the very bottom of maze.inp, found in \phits\lecture\advanced\shielding\maze 22:12:18 発信元: 088 G.A. Atregenio : Replying to "how can we set an en..." @*Tutor* Hunter Ratliff I see, Thank you sensei 22:14:01 発信元: 214 S. Schmidt : I still did not fully understand the difference in how effective dose is calculated by T-Cross and T-point. Why does T-Point converge much faster? 22:14:36 発信元: wakai : Replying to "Which folder is t-po..." thank you. 22:14:40 発信元: 083 Firliyani : Replying to "how can we set an en..." @*Tutor* Hunter Ratliff But what if we only want to calculate the particle flux at a specific energy, for example in Co-60, only in the energy range 1 - 1.33 MeV. Is it possible? 22:15:32 発信元: 103 A.HUGUES : does it mean that t-point is similar to a F5 tally in MCNP ? 22:16:24 発信元: *Tutor* Shintaro Hashimoto (sub, chat only) : Replying to "does it mean that t-..." Yes! 22:16:44 発信元: 144 Margarita Tzivaki : Antworten an "does it mean that ..." this is what I am wondering as well. In mcnp they can be somewhat statistically unstable and require some fiddling with the exclusion volume. How is that handled in phits? Or is this not an issue? 22:16:45 発信元: 103 A.HUGUES : Reacted to "Yes!" with 🙏 22:17:28 発信元: wakai : Could you teach the meaning of (-200-150) in the "3 200-c7 0 -200-150 1.0" in page 15? 22:17:40 発信元: *Tutor* Hunter Ratliff : Replying to "I still did not full..." T-point is a fairly special tally, and the math behind it is a bit too much to write out here. If you'd like to understand what it's doing in more detail, please see the \phits\utility\tpoint\readme-en.docx file. 22:18:51 発信元: 240 Youngran Kwon : I can't activate t point. 22:19:11 発信元: 240 Youngran Kwon : The result is not different. 22:19:17 発信元: 214 S. Schmidt : Replying to "I still did not full..." I see, thank you! And what is this mat parameter about? mat mset1 22:20:13 発信元: *Tutor* Shintaro Hashimoto (sub, chat only) : Replying to "does it mean that t-..." I’m not sure if it’s the same issue as the exclusion volume problem, but the use of t-point tally requires caution, especially in shielding calculations with complex geometries. 22:21:32 発信元: *Tutor* Hunter Ratliff : Replying to "how can we set an en..." You can use the energy mesh in your tallies to limit the scope of what is scored, if that is what you mean. 22:21:32 発信元: 214 S. Schmidt : Reacted to "T-point is a fairly ..." with 👍 22:22:01 発信元: *Tutor* Shintaro Hashimoto (sub, chat only) : This simply means -350. 22:22:29 発信元: *Tutor* Hunter Ratliff : Replying to "how can we set an en..." If you only wanted to simulate specific emissions of an isotope (like the 1.333 MeV emission from Co-60 but not the 1.173 MeV line), then in that case, I think you should just manually specify the discrete emission energies and probabilities of interest with e-type=8, as had been done in slide 25 of the source lecture. 22:23:57 発信元: 144 Margarita Tzivaki : Antworten an "does it mean that ..." Yes, that sounds like it is a similar issue 22:24:10 発信元: *Tutor* Shintaro Hashimoto (sub, chat only) : Replying to "The result is not di..." Do you want to move to a breakout room? 22:24:41 発信元: 264 + Kebede Shogile : Where we are using this excel file? 22:25:05 発信元: wakai : Replying to "Could you teach the ..." thank you. 22:25:06 発信元: *Tutor* Shintaro Hashimoto (sub, chat only) : Or please check whether you deleted "off" in the [t-point] line or not. 22:25:16 発信元: *Tutor* Hunter Ratliff : Replying to "I still did not full..." This is part of the multiplier section and is what is used to rescale the results to the units we want (3600*1e-6) and with the dose coefficient set we want (-202) 22:25:57 発信元: 214 S. Schmidt : Reacted to "This is part of the ..." with 👍 22:26:57 発信元: *Tutor* Shintaro Hashimoto (sub, chat only) : Replying to "Where we are using t..." \phits\lecture\advanced\shielding 22:27:18 発信元: *Tutor* Shintaro Hashimoto (sub, chat only) : Replying to "Where we are using t..." The file name is AnalyticalMethod.xlsx 22:27:49 発信元: 264 + Kebede Shogile : Replying to "Where we are using t..." I mean, are we using this excel file in the accelerator.inp? Or why we are editing? 22:28:21 発信元: 144 Margarita Tzivaki : Antworten an "does it mean that ..." is there an option to set an exclusion radius in terms of free path lengths? in phits? 22:28:50 発信元: *Tutor* Hunter Ratliff : The idea is just to compare the "rule of thumb" analytical/empirical formulas commonly used for speed (hand calculations) against the results from PHITS calculations 22:31:16 発信元: 083 Firliyani : if protons hit O-18 target and produce photons, do we use the tally t-point to calculate numbers of photon flux? 22:31:39 発信元: *Tutor* Hunter Ratliff : \phits\XS\jendl5_setup_* 22:31:48 発信元: 083 Firliyani : Reacted to "\phits\XS\jendl5_set..." with 👍 22:32:12 発信元: *Tutor* Shintaro Hashimoto (sub, chat only) : Replying to "does it mean that t-..." @144 Margarita Tzivaki No. We have to set r0 directly. 22:32:36 発信元: 264 + Kebede Shogile : Replying to "Where we are using t..." Okay, got it 22:32:47 発信元: 240 Youngran Kwon : I setuped Jendle5. 22:33:15 発信元: 083 Firliyani : Reacted to "Yes!" with 👍 22:33:38 発信元: 144 Margarita Tzivaki : Hat auf "@144 Margarita Tzi..." mit 👍 reagiert 22:34:22 発信元: *Tutor* Hunter Ratliff : Replying to "if protons hit O-18 ..." I think this should be possible as long as particle energies are within the available data library files for the relevant material isotopes. 22:34:24 発信元: wakai : In the other calculation of materials dpa, how do we check/know the threshold energy of material? 22:36:00 発信元: 083 Firliyani : Reacted to "I think this should ..." with 👍 22:36:04 発信元: *Tutor* Hunter Ratliff : Replying to "if protons hit O-18 ..." If this doesn't work, you could do a 2-part PHITS simulation where you (1) run one simulation to score produced photons to a dump file and then (2) a second simulation that loads that dump file as source and is only responsible for transporting those photons, using the [T-Point] tally 22:39:46 発信元: *Tutor* Hunter Ratliff : Replying to "In the other calcula..." I'm not sure if there's a super convenient resource for seeing what energy ranges are available in nuclear data files for individual nuclides. 22:40:22 発信元: *Tutor* Hunter Ratliff : Replying to "In the other calcula..." On the JENDL-5 webpage, you can get some of this information if you do some digging (and of course presuming that you're using JENDL-5). https://wwwndc.jaea.go.jp/ftpnd/jendl/jendl-5-n.html 22:41:30 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "In the other calcula..." If you mean displacement threshold energies, please see phits/src/dpa.f subroutine dpath(z2,ttd,ieth) you can find displacement threshold f various elements. 22:42:07 発信元: *Tutor* Hunter Ratliff : Reacted to "If you mean displace..." with 🙏 22:43:18 発信元: 240 Youngran Kwon : Replying to "The result is not di..." I deleted "off". but result is not different. 22:44:07 発信元: wakai : In the other calculation of materials dpa, is there arc-dpa and dpa as input.inp file? I want to know/calculate the result's difference. 22:44:39 発信元: 214 S. Schmidt : The name of the output file for the T-Point is different I believe, have you checked that? 22:46:56 発信元: *Tutor* Hunter Ratliff : Replying to "In the other calcula..." I'm not expert in materials damage, but the [T-DPA] setting "idpa" is used to toggle between NRT-dpa and arc-dpa. 22:46:56 発信元: 240 Youngran Kwon : Replying to "The result is not di..." OK thanks 22:47:24 発信元: 240 Youngran Kwon : Reacted to "The name of the outp..." with 😂 22:47:36 発信元: *Tutor* Hunter Ratliff : Replying to "In the other calcula..." (idpa=0, dpa without athermal recombination correction (NRT-dpa), is the default setting) 22:47:57 発信元: 083 Firliyani : Is it possible to change the "reg 11" cell material from air to detector crystal material (NaITl), and still use tally t track, or we have to use t deposit? 22:48:14 発信元: *Tutor* Hunter Ratliff : Replying to "In the other calcula..." \phits\sample\tally\t-dpa contains some sample input using [T-DPA] 22:49:29 発信元: *Tutor* Shintaro Hashimoto (sub, chat only) : please also check the paper [Y. Iwamoto, S. Meigo, and S. Hashimoto, J. Nucl. Mater. 538 (2020) 152261. https://doi.org/10.1016/j.jnucmat.2020.152261 ] to confirm the difference between the NRT- and arc-dpa models. 22:50:48 発信元: 214 S. Schmidt : Reacted to "OK thanks" with 👍 22:54:33 発信元: *Tutor* Hunter Ratliff : Replying to "Is it possible to ch..." The [T-Track] tally is using a fluence-to-dose conversion coefficient set (see the multiplier section) to calculate dose. If you want to score the energy deposition in the detector, [T-Deposit] is the right tally to use. 22:56:53 発信元: *Tutor* Hunter Ratliff : Replying to "Is it possible to ch..." Tomorrow's detector response lecture will go into deeper detail for modelling detectors. (If you want to look in advance, see \phits\lecture\advanced\detector) 22:57:30 発信元: wakai : Replying to "In the other calcula..." Thank you many thing information. 23:01:49 発信元: 060 Dewa Ngurah Yudhi Prasada : t-cross and t-track actually score the same unit (fluence) but t-cross is in surface and t-track is in volume. Am I correct? 23:03:10 発信元: 264 + Kebede Shogile : I don't why I am not getting cross_eng.eps 23:03:27 発信元: *Tutor* Hunter Ratliff : That's right. In technical terms, [T-Track] is a track length estimator for fluence while [T-Cross] calculates surface current (divided by cos(theta) for fluence). 23:03:47 発信元: 060 Dewa Ngurah Yudhi Prasada : Replying to "t-cross and t-track ..." thanks!! 23:04:44 発信元: *Tutor* Hunter Ratliff : Replying to "I don't why I am not..." Do you have set epsout = 1 in the corresponding [T-Cross] tally? 23:05:10 発信元: wakai : Replying to "In the other calcula..." Thank you. I will check them. 23:09:05 発信元: 264 + Kebede Shogile : Replying to "I don't why I am not..." Yeah. However, I couldn't get 23:09:26 発信元: 083 Firliyani : sorry for exercise 3, does it mean we don't need to change dump value for any case 23:09:38 発信元: 264 + Kebede Shogile : Replying to "I don't why I am not..." Now solved! 23:10:53 発信元: 171 N. Stoffle : With the dump = -11 in source matching the t-cross generation dump definition here, is this then writing the same dmp.out data over the previous dmp.out file as suggested in the lecture discussion on dump source filtering? And how does that behavior change relative to using dump = 0? 23:11:08 発信元: *Tutor* Hunter Ratliff : Replying to "sorry for exercise 3..." Correct, we are not changing the listed dump parameter values. I think they were in red text on that slide just to highlight them for talking about them 23:11:39 発信元: *Tutor* Hunter Ratliff : Replying to "I don't why I am not..." great! 23:11:44 発信元: William Asztalos : The fact that maxcas and maxbch are fixed in the second part of a two step simulation is making me wonder. In practice, do we always have the dumpfile read more than once? Because, if not, aren't we essentially just repeating the first step of the simulation with a shorter calculation time? The proton interaction is skipped, but ultimately we do not have more particles making it to the shielding, and while the results will be different from the random number generator, it's not as though we are sampling more times. 23:12:16 発信元: William Asztalos : Replying to "The fact that maxcas..." Wordy question, I know, but does that make sense? I guess I'm just asking if we should always read the dumpfile multiple times to increase statistics. 23:13:52 発信元: 083 Firliyani : Reacted to "Correct, we are not ..." with 👍 23:17:48 発信元: *Tutor* Hunter Ratliff : Replying to "The fact that maxcas..." This is a common application of the dump source. The physics for modelling the beam interactions and secondary particle production are far more expensive than the transport of those secondaries in many cases. So, if you have some scenario where the beam+target is kept the same but you're wanting to test other variations (like moving a detector around, changing shielding, etc.), it can be more convenient just to run the slow production physics once then just reload the produced particles in the various detector/shielding/etc. permutations you want to simulate. 23:19:10 発信元: *Tutor* Hunter Ratliff : Replying to "The fact that maxcas..." In general, I'd make sure your initial simulation generated enough secondaries to be used as a future source to at least fairly represent the distribution of the produced secondaries. 23:19:54 発信元: William Asztalos : Replying to "The fact that maxcas..." Ah okay, makes sense, thank you! If you have a sufficient number of histories and batches in the first step, does it make sense to also read the dump file multiple times? I imagine it's important not to do that too much though. 23:20:15 発信元: William Asztalos : Replying to "The fact that maxcas..." But now the usecase is clear! 23:21:39 発信元: *Tutor* Hunter Ratliff : Replying to "The fact that maxcas..." Yes, that's exactly the right way to think about it. The ideal case perhaps is to only read through it once (in fact, this is the default behavior), but you can use the dmpmulti parameter to make it recycle the dump file some extra number of times. 23:27:50 発信元: *Tutor* Shintaro Hashimoto (sub, chat only) : Replying to "With the dump = -11 ..." When you keep the [t-cross] tally activated to create the dmp.out file, this file will be overwritten during the second-step calculation. As for dump=0, this option can be used in the latest PHITS version. The lecture materials are a bit old - they were made using an older version of PHITS. 23:28:35 発信元: 264 + Kebede Shogile : How we download cross section libraries and add to PHITS? 23:29:55 発信元: *Tutor* Hunter Ratliff : Replying to "With the dump = -11 ..." See the contents of \phits\XS There is an installer present for downloading JENDL-5 data 23:30:15 発信元: 140 Mahdi Bakhtiari : What if we do not use 2 step calculations. And use 1 step calculations with WWG? Wouldn't it be efficient than two step case? 23:30:46 発信元: *Tutor* Hunter Ratliff : Replying to "With the dump = -11 ..." Regarding dump file overwriting, if you are running a "restart calculation" with istdev < 0 (we aren't doing this here), the dump file is appended to; otherwise (as is the case here), it is overwritten. 23:32:11 発信元: 171 N. Stoffle : Reacted to "Regarding dump file ..." with 👍 23:32:23 発信元: 171 N. Stoffle : Reacted to "When you keep the [t..." with 👍 23:33:46 発信元: *Tutor* Hunter Ratliff : Replying to "What if we do not us..." I think which is preferable is a bit dependent on your use case. If you intend to re-use the dump source for many different calculations, then the 2-step calculation is probably preferable. They are both just two different methods to help achieve results with less compute/time cost. 23:34:19 発信元: 140 Mahdi Bakhtiari : Reacted to I think which is pre... with "👍" 23:34:56 発信元: *Tutor* Hunter Ratliff : The 2-step calculation saves cost in repeating complicated reaction physics behind secondary radiation production, and WW reduces cost spent transporting particles in uninteresting directions. 23:35:56 発信元: 140 Mahdi Bakhtiari : Replying to "What if we do not us..." Yes that's correct 23:45:13 発信元: 052 Dali Georgobiani : Thank you very much, Ogawa-san!!! 23:45:32 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "Thank you very much,..." with 👍 23:45:37 発信元: 107 Iaroslav Bielov : Thank you!! 23:45:40 発信元: 140 Mahdi Bakhtiari : Thank you so much 23:45:54 発信元: 214 S. Schmidt : Great lecture, thank you! 23:45:55 発信元: 103 A.HUGUES : Thank you! See you tomorrow! 🎃 23:46:05 発信元: 211_Somchai_Tan : Thank you very much 23:46:21 発信元: 129 - K. Butalag : I would like to know if I can use Jendel 5 to asses the dose due to neutron and photon , studying the shielding for a proton accelerator for isotope production. Proton energies between 7 and 70 MeV 23:46:23 発信元: 083 Firliyani : Thank you sensei, i think all of my questions have been answered 23:46:26 発信元: 060 Dewa Ngurah Yudhi Prasada : Thank you Ogawa-san, I hope we meet again in Japan!! 23:46:32 発信元: 150 Melih Tahsin Değer Main : It was Great session, thank you very much 23:46:35 発信元: 256SergioCapdevila : I am trying to execute DCHAIN but this error happen "findstr" no se reconoce como un comando interno o externo, 23:46:35 発信元: 114 Gordeev Main : Thank you very much, goodbye 23:46:51 発信元: 256SergioCapdevila : You know the path of findstr? 23:47:08 発信元: 088 G.A. Atregenio : Thank you very much sensei! 23:47:15 発信元: wakai : Thank you very much, Ogawa-san 23:47:21 発信元: *Tutor* Hunter Ratliff : Replying to "I would like to know..." Yes, JENDL-5 proton data libraries are up to 200 MeV. 23:47:22 発信元: 160 M. Adithya P : thank you 23:47:38 発信元: *Tutor* Hunter Ratliff : Replying to "I would like to know..." https://wwwndc.jaea.go.jp/ftpnd/jendl/jendl-5-p.html