20:00:40	 発信元: 099 hyoriyaz : good morning!
20:01:06	 発信元: 240 Youngran Kwon : Good evening ^^
20:01:14	 発信元: *Tutor* Tatsuhiko OGAWA : 1.Recording and guidance of this tutorial is available from here http://phits.jaea.go.jp/tutorial/rec.html
2.Please install Python and its libraries to take the optional course “Medical data treatment by RT-PHITS” on Friday
3.Please replace all files in \phits\lecture\advanced\DCHAIN1 with those in https://git.oecd-nea.org/phits/package/phits-lecture/-/tree/main/advanced/DCHAIN1 . They have been drastically updated.
20:14:24	 発信元: 115 J. Santos : What would be the weight of the particle if the Wlow in cell 4 would be for example 1.5?
20:16:35	 発信元: 028 Amiegbereta Edwin Ehis : Good day. just joining in. Please what's the topic of discuss and where can I find it?
20:16:36	 発信元: 264 + Kebede Shogile : How to set the right weight for [t-wwg]?
20:17:02	 発信元: Anton Moroz : Replying to "Good day. just joini..."

phits\lecture\advanced\variance-reduction
20:17:13	 発信元: 028 Amiegbereta Edwin Ehis : Replying to "Good day. just joini..."

thanks
20:17:17	 発信元: 110 I. Casalbore : Is wupn always supposed to be 5 or does it depend on the case?
20:17:28	 発信元: *Tutor* Hunter Ratliff : Replying to "Good day. just joini..."

and "WW-concept" section right now
20:18:00	 発信元: *Tutor* Tatsu Sato : Replying to "What would be the we..."

Chance to survive becomes lower and weight should be over 1.5
20:18:23	 発信元: *Tutor* Hunter Ratliff : Replying to "How to set the right..."

We will get to using [t-wwg] later in this lecture :)
20:19:07	 発信元: *Tutor* Tatsu Sato : Replying to "Is wupn always suppo..."

You can change the value from your input file, but in most cases, it is not necessary to be changed.
20:20:12	 発信元: *Tutor* Yusuke Matsuya : Replying to "How to set the right..."

You can learn about [t-wwg] p38-49.
20:20:29	 発信元: 110 I. Casalbore : Reacted to "You can change the v..." with ❤️
20:51:06	 発信元: 152 Meutia Rahmi Afifa : so we run it three times?
20:51:55	 発信元: *Tutor* Ohnishi Seiki : Replying to "so we run it three t..."

Yes we do.
20:52:49	 発信元: 208 S.Laca : it says wwg.out doesn't exist. what do I do?
20:53:02	 発信元: Anton Moroz : can we use different meshes for particles with different range?
20:53:15	 発信元: 110 I. Casalbore : What is the purpose of wwg-xz.out? If I understand correctly wwg.out is showing the weight window section, but I do not understand wwg-xz
20:53:39	 発信元: *Tutor* Ohnishi Seiki : Replying to "it says wwg.out does..."

After the first run, it is created in the same folder as the input file.
20:55:14	 発信元: *Tutor* Ohnishi Seiki : Replying to "can we use different..."

Yes we can, but the mesh should cover the important regions.
20:55:44	 発信元: *Tutor* Takuya Furuta : Replying to "What is the purpose ..."

For the file = wwg-xz.out, axis = xz is specified above. So the flux results are given to the wwg-xz.out file while weight window values are given in wwg.out.
20:56:06	 発信元: Anton Moroz : Replying to "can we use different..."

for this do I need to create two t-wwg sections? 
I have protons that produce alphas, and I want to transport alphas (10 um range)
20:56:39	 発信元: *Tutor* Takuya Furuta : Replying to "What is the purpose ..."

The purpose of the wwg-xz.out file is to check the flux in EPS figure.
20:58:03	 発信元: *Tutor* Takuya Furuta : Replying to "can we use different..."

Yes. Only one mesh can be specified in one [t-wwg]. So you need two t-wwg.
20:58:12	 発信元: Anton Moroz : Reacted to "Yes. Only one mesh c..." with 👏🏼
20:58:13	 発信元: 110 I. Casalbore : Replying to "What is the purpose ..."

So this is because t-wwg calculates the flux in order to create the weight window?
20:58:41	 発信元: *Tutor* Yusuke Matsuya : Replying to "What is the purpose ..."

In additon, now we focus on the xz meshes. You are checking the flux in xz meshes. That is why we use the www-xz.out.
20:59:07	 発信元: *Tutor* Takuya Furuta : Replying to "What is the purpose ..."

Yes. That's right. So with axis=xz, you can check the flux results as well.
20:59:27	 発信元: 110 I. Casalbore : Replying to "What is the purpose ..."

Thank you very much :)
21:00:17	 発信元: *Tutor* Yusuke Matsuya : Replying to "it says wwg.out does..."

Yes. We need to rform 2 step calculations!
21:00:26	 発信元: Anton Moroz : does the mesh in t-track have to match in some way to the mesh in t-wwg?
21:00:39	 発信元: *Tutor* Yusuke Matsuya : Replying to "it says wwg.out does..."

Yes. We need to perform 2 step calculations!
21:01:05	 発信元: *Tutor* Ohnishi Seiki : Replying to "does the mesh in t-t..."

No. You can use different meshes in t-track and t-wwg.
21:01:32	 発信元: Anton Moroz : Replying to "does the mesh in t-t..."

e.g. I have a 10um resolution for t-wwg. will I see the particles in t-track with 10 mm resolution?
21:03:50	 発信元: 144 Margarita Tzivaki : I am not sure I understand how the T-wwg section setup translates to the actually produced ww-mesh that we see in the wwg.out file. Would you mind a quick recap?
21:08:03	 発信元: *Tutor* Ohnishi Seiki : Replying to "does the mesh in t-t..."

You can use finer mesh as long as enough number of particles enters. However, too small mesh may fail because of numerical errors.
21:08:18	 発信元: Anton Moroz : Replying to "does the mesh in t-t..."

thank you, this is very helpful
21:09:54	 発信元: 264 + Kebede Shogile : When we need to use weight window in addition to shielding?
21:10:07	 発信元: Anton Moroz : can a mesh be nonuniform? for example it becomes finer with increasing depth
21:10:52	 発信元: *Tutor* Ohnishi Seiki : Replying to "can a mesh be nonuni..."

Of cause you can. Non-equi-interval meshes can be used.
21:11:02	 発信元: Anton Moroz : Reacted to "Of cause you can. No..." with 👍🏼
21:12:34	 発信元: *Tutor* Takuya Furuta : The purpose of the Weight Window Generator is to produce weight window factors to achieve uniform statistics everywhere. The statistics is determined by the number of particles calculating the quantities. Therefore what we need to do is to have weight window factor such that particle number coming to the region is constant every where. In principle, this can be achieved when we know the flux values at each location. Namely weight=flux. This is the idea of the T-wwg.
21:15:19	 発信元: *Tutor* Tatsu Sato : Replying to "When we need to use ..."

In the case that you feel computational time is too long
21:15:24	 発信元: 144 Margarita Tzivaki : thank you!
21:16:22	 発信元: 264 + Kebede Shogile : Please, would you explain more the relation between wwg and t-wwg?
21:18:31	 発信元: *Tutor* Ohnishi Seiki : Replying to "Please, would you ex..."

t-wwg is Tally function for Weight Window Generation. It generates appropriate Weight Window lower values set.
21:20:07	 発信元: 115 J. Santos : in how many iterations do you accomplish certain convergence in the values of the wwg? I guess it could depend on the geometry
21:21:54	 発信元: 144 Margarita Tzivaki : do you have any guidance how to set up the t-wwg section? as in how large to set the mesh (only a region of the geometry or the whole geometry?) and how to choose nx, ny, nz?
21:22:23	 発信元: *Tutor* Ohnishi Seiki : Replying to "in how many iteratio..."

Your guess is right. It depends on the geometry, particle, source energy and so on. That's somewhat empirical.
21:22:48	 発信元: 144 Margarita Tzivaki : Hat auf "The purpose of the..." mit 👍 reagiert
21:22:50	 発信元: 115 J. Santos : Reacted to "Your guess is right...." with 👍
21:25:26	 発信元: *Tutor* Tatsu Sato : Replying to "do you have any guid..."

There is no good instruction because there is no universally appropriate setting (i.e. case-by-case). During this lecture, we will explain a little bit, but you have to learn by yourself based on your own experience.
21:25:38	 発信元: *Tutor* Yusuke Matsuya : Replying to "do you have any guid..."

I think it depends on persons (users). I think the setting method differs depending on what kind of calculation you want to do.
21:27:02	 発信元: Anton Moroz : when we update t-wwg do we have to disable infl:{wwg.out} so that it does not interfere with the new weight window generation?
21:27:04	 発信元: *Tutor* Ohnishi Seiki : Replying to "do you have any guid..."

The mesh doesn’t necessarily have to cover the entire geometry. It should cover the "Important" region. If which part of geometry is important is not clear, wider ranges of meshes are used.
21:27:23	 発信元: 103 A.HUGUES : how do you know if phits re-generate wwg.out or update it ?
21:28:03	 発信元: 144 Margarita Tzivaki : Hat auf "The mesh doesn’t ..." mit 👍 reagiert
21:28:49	 発信元: *Tutor* Ohnishi Seiki : Replying to "when we update t-wwg..."

infl:{wwg.out} is needed to use previous WW set. The iterations make better weight window data progressively.
21:29:27	 発信元: Anton Moroz : Replying to "when we update t-wwg..."

I mean if we change wwg parameters like now by including elowthre
21:29:45	 発信元: *Tutor* Takuya Furuta : At the beginning of PHITS run the weight window input is stored in memory so that the newly created wwg.out file does not interfere. So no update of weight window is happening inside one PHITS run.
21:30:26	 発信元: 103 A.HUGUES : Reacted to "At the beginning of ..." with 🙏
21:31:54	 発信元: *Tutor* Ohnishi Seiki : Replying to "when we update t-wwg..."

If you start from scaratch, infl should be commented out. Actually, imperfect ww data may be better than nothing if the changes are small ( but it depends on probrems...).
21:32:33	 発信元: 144 Margarita Tzivaki : how would I go about setting up a non-uniform mesh? Or does this happen automatically during the repeat of the calculations?
21:32:33	 発信元: Anton Moroz : Reacted to "If you start from sc..." with 👍🏼
21:33:38	 発信元: Anton Moroz : Replying to "can a mesh be nonuni..."

I just checked that you can input your mesh manually with x-type=1
21:35:17	 発信元: *Tutor* Ohnishi Seiki : Replying to "can a mesh be nonuni..."

Yes, the meshes can be set in the same syntax as the tally mesh definition like x-type=, y-type, and so on. The mesh size cannot be cahged automatically.
21:35:30	 発信元: *Tutor* Takuya Furuta : Replying to "can a mesh be nonuni..."

Yes, that's right. With x-type=1, you can have non-equal size mesh. That can be used also in t-wwg. The mesh is written in weight window so the same mesh will be used in weight window and t-wwg.
21:35:48	 発信元: 144 Margarita Tzivaki : Hat auf "I just checked tha..." mit 👍 reagiert
21:36:02	 発信元: Anton Moroz : Reacted to "Yes, that's right. W..." with 👍🏼
21:36:05	 発信元: 144 Margarita Tzivaki : Hat auf "Yes, that's right...." mit 👍 reagiert
21:42:47	 発信元: 144 Margarita Tzivaki : I have had issues when using an iterative weight window generator in geometries with low intensity sources and high absorption shielding. I think the issue was that not enough particles were going all the way through to correctly calculate the ww. A couple of years ago someone in an MCNP course suggested to set the density to very low and then iteratively set up the density while running weight window generator runs to optimize. Is this still a valid strategy or has this become obsolete?
21:44:50	 発信元: 144 Margarita Tzivaki : Antworten an "I have had issues ..."

and follow-up to this: is there any external utilities or programs that allow us to use FW-CADIS with phits that you can recommend?
21:49:39	 発信元: *Tutor* Ohnishi Seiki : Replying to "I have had issues wh..."

Such simple method is still valid.  Of course, in PHITS you can adjust material densities easily by constants function like set:c1[0.1]. After all, empirical factors such as the density reduction width and the number of histories per iteration still remain. FW-CADIS eliminate such empirical factors by the deterministic transport calculation. If the geometry is simple, the simple empirical method do well, and for a complex geometry, somehow complex method like FW-CADIS may be better.
21:53:02	 発信元: *Tutor* Ohnishi Seiki : FW-CADIS tools for PHITS are not developed yet. Some researchers try to apply the results of the FW-CADIS code, ADVANTG, to PHITS. The geometry input of PHITS and MCNP is similar.
21:55:23	 発信元: 144 Margarita Tzivaki : Hat auf "FW-CADIS tools for..." mit 👍 reagiert
21:55:44	 発信元: 126 J.Omojola : If I am to calculate DPA of cladding materials use  in LWRs with operational life time of let say 20 years. How can I set this operational time for my calculation in PHITS?
21:56:07	 発信元: 144 Margarita Tzivaki : Antworten an "I have had issues ..."

thank you I will have to try the density set through a function. Are there any plans for FW-CADIS integration in phits?
21:56:35	 発信元: Anton Moroz : what is the nature of the Si ions? are they knocked off from the target? Does phits track those kinds of reactions?
21:56:37	 発信元: *Tutor* Tatsu Sato : Replying to "I have had issues wh..."

We are working with a private company for developing a tool like ADVANTG, but the outcome will be a commercial product.
21:57:06	 発信元: 103 A.HUGUES : Reacted to "FW-CADIS tools for P..." with 👍
21:57:14	 発信元: 144 Margarita Tzivaki : Hat auf "We are working wit..." mit 👍 reagiert
21:57:23	 発信元: 103 A.HUGUES : Reacted to "We are working with ..." with 👍
21:58:46	 発信元: *Tutor* Ohnishi Seiki : Replying to "If I am to calculate..."

The PHITS output is "DPA per source particles". So, you can multiply it by the total number of neutrons created during 20 years.
21:58:55	 発信元: 100 Herlan Setiawan : In the case of using t-cross output, when do we use current or flux, in the case of detecting a nucleus that is recoiled in all directions from the surface of materials, which one should be used?
21:59:32	 発信元: *Tutor* Tatsu Sato : Replying to "what is the nature o..."

Yes, but the range of knocked out Si ions are so short and they are mostly stopped in the target
21:59:49	 発信元: Anton Moroz : Reacted to "Yes, but the range o..." with 👍🏼
22:00:27	 発信元: 270 Shupti Sarker : Is it possible to apply the weight window technique in any type of PHITS calculation? For instance, if I aim to optimize the computation time in an absorbed dose calculation, would using a weight window be beneficial?
22:00:36	 発信元: 052 Dali Georgobiani : Thank you very much, Ogawa-san!!!
22:00:39	 発信元: 144 Margarita Tzivaki : Antworten an "If I am to calcula..."

is there a good way to account for changing neutron spectrum over time besides doing multiple phits runs?
22:01:39	 発信元: *Tutor* Tatsu Sato : If you simply want to know the number of particles incident to a tally region, you have to specify "output = current". If you want to know the fluence/flux at the region, you have to specify "output = flux".
22:01:41	 発信元: 056 David Ferreira : Could you share the link to donwload DChain files?
22:01:49	 発信元: *Tutor* Yusuke Matsuya : Replying to "In the case of using..."

In case of detecting the recoiled necleus in all direction, I think the [t-product] can be useful.
22:02:44	 発信元: 100 Herlan Setiawan : Replying to "In the case of using..."

ok thank you
22:02:46	 発信元: *Tutor* Hunter Ratliff : https://git.oecd-nea.org/phits/package/phits-lecture/-/tree/main/advanced/DCHAIN1
22:03:16	 発信元: *Tutor* Tatsu Sato : Replying to "Is it possible to ap..."

In principle yes, but practically not in most cases. If your [weight window] setting is inappropriate, the computational time becomes even longer than that without [weight window]
22:04:45	 発信元: *Tutor* Tatsu Sato : Replying to "If I am to calculate..."

Run PHITS separately and add tally results using sumtally function
22:05:32	 発信元: 056 David Ferreira : Reacted to "https://git.oecd-nea..." with 👍🏻
22:05:35	 発信元: 043 A.Bekouale : Thanks you Mr OGAWA
22:05:55	 発信元: 270 Shupti Sarker : Replying to "Is it possible to ap..."

I see. Thank you, sensei.
22:07:00	 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "Thanks you Mr OGAWA " with 👍
22:07:10	 発信元: 043 A.Bekouale : Reacted to Thanks you Mr OGAWA  with "👍"
22:07:11	 発信元: 043 A.Bekouale : Removed a 👍 reaction from "Thanks you Mr OGAWA "
22:10:55	 発信元: *Tutor* Hunter Ratliff : https://git.oecd-nea.org/phits/package/phits-lecture/-/tree/main/advanced/DCHAIN1
22:16:02	 発信元: 126 J.Omojola : Reacted to "The PHITS output is ..." with 👍
22:17:17	 発信元: 264 + Kebede Shogile : Are some nuclear libraries missing?
22:18:16	 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "Are some nuclear lib..."

As Hunter said, yes. But it is not an issue after all because nuclear reaction models work in place of nuclear data.
22:18:42	 発信元: 264 + Kebede Shogile : Replying to "Are some nuclear lib..."

Okay! Thank you!
22:18:51	 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "Okay! Thank you!" with 👍
22:19:16	 発信元: 072 E. Musacchio Gonzalez : Replying to "Are some nuclear lib..."

It doesn’t produce any files in my case
22:19:57	 発信元: 236 W.Asztalos : as a head's up for others, I needed to add a epsout=1 option to my t-dchain tally definition
22:20:00	 発信元: 264 + Kebede Shogile : Replying to "Are some nuclear lib..."

I think we need to install DCHAIN separetly?
22:20:12	 発信元: 144 Margarita Tzivaki : Hat auf "as a head's up for..." mit ❤️ reagiert
22:20:29	 発信元: *Tutor* Tatsu Sato : Replying to "Are some nuclear lib..."

No. If you install PHITS, DCHAIN is automatically installed. Please ask breakout room if you cannot get any output files
22:21:20	 発信元: 264 + Kebede Shogile : Replying to "Are some nuclear lib..."

Yeah! I got it!
22:21:22	 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "as a head's up for o..." with 👍
22:21:25	 発信元: 129 -K. Butalag : Replying to "Are some nuclear lib..."

It happened to me that I got permission denied, try to give permission to the dchain.sh
22:21:30	 発信元: 072 E. Musacchio Gonzalez : Replying to "Are some nuclear lib..."

I downloaded and copied again the file and now works, I think
22:21:41	 発信元: 129 -K. Butalag : Replying to "Are some nuclear lib..."

perfect
22:21:59	 発信元: *Tutor* Tatsu Sato : Reacted to "I downloaded and cop..." with 👍
22:22:09	 発信元: 043 A.Bekouale : I don't see the tdchain.out file in my phits folder.please..
22:23:05	 発信元: *Tutor* Tatsu Sato : Replying to "I don't see the tdch..."

Please ask breakout room if you cannot solve the issue. I guess that you run PHITS in the zipped file that you downloaded.
22:23:37	 発信元: *Tutor* Tatsu Sato : Replying to "I don't see the tdch..."

First, you have to extract them
22:26:29	 発信元: 043 A.Bekouale : Replying to "I don't see the tdch..." 

 I extracted them from the zip file, but so far I don't see them...
22:27:36	 発信元: *Tutor* Tatsu Sato : Replying to "I don't see the tdch..."

OK. We will take you to a breakout room
22:28:23	 発信元: 144 Margarita Tzivaki : Journal of Open Source Software: The PHITS Tools Python package for parsing, organizing, and analyzing results from the PHITS radiation transport and DCHAIN activation codes
22:28:55	 発信元: 096 Handy Nugraha : is there any tools to calculate reg volume in phits? like stochastically?
22:28:57	 発信元: 052 Dali Georgobiani : Replying to "Journal of Open Sour..."

Thank you, Margarita!
22:29:04	 発信元: *Tutor* Ohnishi Seiki : Reacted to "Journal of Open Sour..." with 👍
22:29:18	 発信元: 096 Handy Nugraha : Reacted to "Journal of Open Sour..." with 👍
22:29:54	 発信元: *Tutor* Tatsu Sato : Replying to "is there any tools t..."

PHITS has [t-volume] tally for the volume calculation purpose. Please see lecture/basic/lec03 for more detail.
22:31:10	 発信元: 103 A.HUGUES : Reacted to "Journal of Open Sour..." with 👍
22:31:36	 発信元: 144 Margarita Tzivaki : Hat auf "PHITS has [t-volum..." mit 👍 reagiert
22:35:33	 発信元: 060 Dewa Ngurah Yudhi Prasada : "amp" is number of the particles simulated/sec isn't it? so if I simulate neutron, we can just put the "amp" with the value without adding the "beam current" right?
22:36:14	 発信元: 096 Handy Nugraha : Replying to "is there any tools t..."

okay, thanks.
22:39:18	 発信元: *Tutor* Tatsu Sato : Replying to ""amp" is number of t..."

Yes. This is the sample file for proton accelerator
22:39:28	 発信元: 060 Dewa Ngurah Yudhi Prasada : Reacted to "Yes. This is the sam..." with 👍🏻
22:39:41	 発信元: *Tutor* Tatsu Sato : Replying to ""amp" is number of t..."

For neutron, we recommend directly specifying the neutron intensity per second
22:40:47	 発信元: 114 Gordeev Sub : Replying to ""amp" is number of t..."

And what about heavy ions? Do we need to consider charge? Or just provide number of ions per second?
22:41:03	 発信元: 060 Dewa Ngurah Yudhi Prasada : Reacted to "For neutron, we reco..." with ❤️
22:41:33	 発信元: *Tutor* Tatsu Sato : Replying to ""amp" is number of t..."

Yes. You have to consider charge. The important point is that you have to input "particle/sec" in amp parameter.
22:41:38	 発信元: 144 Margarita Tzivaki : So if I understood this correctly the source term (e.g. the energy) can be changed between different irradiation events in the irradiation Schedule but that needs to be done through the phits file?
22:41:59	 発信元: *Tutor* Hunter Ratliff : Replying to "So if I understood t..."

If you're changing the source energy, this should be done in the PHITS input
22:42:17	 発信元: *Tutor* Hunter Ratliff : Replying to "So if I understood t..."

If you're only changing the intensity, then you can choose to only update the DCHAIN input for that
22:42:29	 発信元: 264 + Kebede Shogile : Do we open tdchain.dout?
22:42:46	 発信元: 110 I. Casalbore : Which are the quantities evaluated by PHITS that are then used in DCHAIN?
22:42:56	 発信元: *Tutor* Hunter Ratliff : Replying to "Do we open tdchain.d..."

ah, no, that is a confusingly similarly named file. "tdchain.out" is the one you want
22:43:29	 発信元: 144 Margarita Tzivaki : Hat auf "If you're changing..." mit 👍 reagiert
22:44:17	 発信元: 110 I. Casalbore : Thank you!
22:44:43	 発信元: 096 Handy Nugraha : Reacted to "Which are the quanti..." with 👍
22:45:02	 発信元: 264 + Kebede Shogile : Replying to "Do we open tdchain.d..."

So, how can we get it after running?
22:46:07	 発信元: 264 + Kebede Shogile : Replying to "Do we open tdchain.d..."

Ooh I got it!
22:51:46	 発信元: 094 G. Franklin : is there a way to comment in the tdchain.out file?
22:52:36	 発信元: *Tutor* Tatsu Sato : Replying to "is there a way to co..."

You can use ! as a comment remark in DCHAIN input file
22:53:44	 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "is there a way to co..."

Also "*" works as a comment symbol.
22:53:52	 発信元: *Tutor* Tatsu Sato : Replying to "is there a way to co..."

You can also use $ in the same as PHITS input file
22:54:35	 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "is there a way to co..."

* works as a comment only when it is put at a beginning of a line. 
$ and ! are more robust in this sense.
22:55:23	 発信元: 060 Dewa Ngurah Yudhi Prasada : Can we determine which part of the specific cell is mostly activated using DCHAIN or do we assume all part of the specific cell activated uniformly?
22:56:47	 発信元: *Tutor* Tatsu Sato : Replying to "Can we determine whi..."

If you use mesh = reg, you cannot estimate the activation distribution inside a tally region. If you want to know more detailed spatial distribution, you have to use mesh = xyz.
22:57:04	 発信元: *Tutor* Tatsu Sato : Replying to "Can we determine whi..."

You will learn this difference later in this lecture
22:57:31	 発信元: 060 Dewa Ngurah Yudhi Prasada : Replying to "Can we determine whi..."

ahhhh i see, thank you
22:57:43	 発信元: 060 Dewa Ngurah Yudhi Prasada : Reacted to "You will learn this ..." with ❤️
22:57:46	 発信元: 060 Dewa Ngurah Yudhi Prasada : Reacted to "If you use mesh = re..." with ❤️
22:57:49	 発信元: 052 Dali Georgobiani : The beam energy number (3 MeV) on tdchain.act should be ignored, right?
22:58:25	 発信元: 144 Margarita Tzivaki : is there a way to set the energy groups in the output file of the photon spectrum?
22:58:29	 発信元: *Tutor* Tatsu Sato : Replying to "The beam energy numb..."

Yes! We should have removed the value...
22:59:07	 発信元: 052 Dali Georgobiani : Replying to "The beam energy numb..."

Thank you, Sato-san! (Also, 3 GeV, not 3 MeV 🙂)
23:00:44	 発信元: *Tutor* Tatsu Sato : Replying to "is there a way to se..."

As far as I know, there is no method to change the energy group because DCHAIN has the gamma-decay database based on the group table. I will check and tell you if I can find the way to change the group.
23:01:05	 発信元: 144 Margarita Tzivaki : Hat auf "As far as I know, ..." mit 👍 reagiert
23:01:18	 発信元: *Tutor* Hunter Ratliff : Replying to "is there a way to se..."

There are a few different presets for this though
23:01:58	 発信元: *Tutor* Hunter Ratliff : Replying to "is there a way to se..."

See the IGGRP parameter in the DCHAIN manual
23:02:36	 発信元: *Tutor* Hunter Ratliff : Replying to "is there a way to se..."

but yes, it is just selecting between predefined group structures, not inputting your own.
23:02:43	 発信元: *Tutor* Tatsu Sato : Replying to "is there a way to se..."

Thank you!
23:03:22	 発信元: *Tutor* Hunter Ratliff : Replying to "is there a way to se..."

DCHAIN can output PHITS RI [Source] sections, which we'll use shortly in the exercise, which of course is more precise than these binned energy groups
23:05:40	 発信元: 199 Renato Carvalho : Hi Prof Ratliff. When I tried to rum DCHAIN, i saw the message: DCHAIN-SP is started at  29/10/2025 11:04:17,02, input= "dchain.inp"
904
 *** error message from s.rdinpt ***
 irregular data key word.
 ===> [Title]  
What could be the error?
23:06:25	 発信元: *Tutor* Hunter Ratliff : Replying to "Hi Prof Ratliff. Whe..."

dchain.inp is a PHITS input file, not a DCHAIN input file. Please run dchain.inp with PHITS, then run the produced tdchain.out file with DCHAIN.
23:07:20	 発信元: 199 Renato Carvalho : Reacted to "dchain.inp is a PH..." with ❤️
23:07:30	 発信元: *Tutor* Hunter Ratliff : Replying to "Hi Prof Ratliff. Whe..."

Apologies if my tongue getting twisted between saying "[T-Dchain]" and "DCHAIN" many times earlier caused any confusion. 😅
23:07:33	 発信元: 199 Renato Carvalho : Right! Thank you veru much! Now it´s OK
23:07:55	 発信元: 199 Renato Carvalho : Reacted to "Apologies if my to..." with ❤️
23:08:22	 発信元: 144 Margarita Tzivaki : Hat auf "but yes, it is jus..." mit 👍 reagiert
23:08:56	 発信元: 144 Margarita Tzivaki : Antworten an "is there a way to ..."

oh that is super neat!
23:19:57	 発信元: Anton Moroz : how can we set up these phits tools for plotting?
23:21:23	 発信元: *Tutor* Tatsu Sato : Replying to "how can we set up th..."

Please see the URL below
https://meteor.nucl.kyushu-u.ac.jp/phitsforum/t/topic/3651
23:21:44	 発信元: Anton Moroz : Reacted to "Please see the URL b..." with 🙏🏼
23:21:47	 発信元: 144 Margarita Tzivaki : Antworten an "how can we set up ..."

there is also a pretty detailled explanation in the github repository where you get it from :)
23:21:55	 発信元: 208 S.Laca : I cant run the Fe_target.out with dchain. It says the system cant find the rute
23:22:00	 発信元: Anton Moroz : Reacted to "there is also a pret..." with 👍🏼
23:23:30	 発信元: *Tutor* Tatsu Sato : Replying to "I cant run the Fe_ta..."

A tutor will introduce you to a breakout room.
23:24:10	 発信元: 208 S.Laca : Respondiendo a "I cant run the Fe_..."

okei
23:24:57	 発信元: 144 Margarita Tzivaki : Antworten an "how can we set up ..."

GitHub - Lindt8/PHITS-Tools: Python package for parsing and post-processing output from the PHITS code
23:25:37	 発信元: 052 Dali Georgobiani : Reacted to "GitHub - Lindt8/PHIT..." with ❤️
23:26:59	 発信元: Anton Moroz : Reacted to "GitHub - Lindt8/PHIT..." with ❤️
23:27:00	 発信元: 144 Margarita Tzivaki : Hat auf "Please see the URL..." mit 👍 reagiert
23:27:45	 発信元: 060 Dewa Ngurah Yudhi Prasada : I missed a bit, which function that generate .pjt files?
23:28:01	 発信元: 060 Dewa Ngurah Yudhi Prasada : Replying to "I missed a bit, whic..."

.pht
23:28:46	 発信元: *Tutor* Tatsu Sato : Replying to "I missed a bit, whic..."

iphtout in PHITS / phitsout in DCHAIN. Please see page 47 in pptx
23:28:59	 発信元: *Tutor* Hunter Ratliff : Replying to "how can we set up th..."

There are a few different ways to do it, but what you're seeing on my end is described here:
https://github.com/Lindt8/PHITS-Tools?tab=readme-ov-file#automatic-processing-at-phits-runtime 
I have modified my phits.bat and dchain.bat files to automatically feed the input file also to PHITS Tools after PHITS/DCHAIN finishes to initiate post-processing automatically.
23:31:31	 発信元: 060 Dewa Ngurah Yudhi Prasada : Replying to "I missed a bit, whic..."

thaaankss Sato-san!
23:31:34	 発信元: 060 Dewa Ngurah Yudhi Prasada : Reacted to "iphtout in PHITS / p..." with ❤️
23:31:46	 発信元: *Tutor* Tatsu Sato : Reacted to "thaaankss Sato-san!" with 👍
23:34:54	 発信元: *Tutor* Hunter Ratliff : Replying to "how can we set up th..."

PHITS Tools also has a GUI if you prefer that which can be used for file processing/plotting.
23:35:13	 発信元: Anton Moroz : Reacted to "There are a few diff..." with 🙏🏼
23:45:46	 発信元: 103 A.HUGUES : can we use $OMP with DCHAIN ?
23:49:15	 発信元: *Tutor* Tatsuhiko OGAWA : Replying to "can we use $OMP with..."

From the next version, to be released in next March or April, OpenMP parallelization is going to be available.
23:49:49	 発信元: 103 A.HUGUES : Reacted to "From the next versio..." with 👍
23:49:58	 発信元: Anton Moroz : Reacted to "From the next versio..." with 👍🏼
23:50:10	 発信元: *Tutor* Hunter Ratliff : Reacted to "From the next versio..." with 🥳
23:50:44	 発信元: *Tutor* Tatsuhiko OGAWA : Reacted to "From the next versio..." with 😄
23:52:58	 発信元: *Tutor* Tatsu Sato : Replying to "can we use $OMP with..."

MPI too
23:56:36	 発信元: 060 Dewa Ngurah Yudhi Prasada : is there any ipltmode to generate xz plot instead :D?
00:00:28	 発信元: 060 Dewa Ngurah Yudhi Prasada : Reacted to "From the next versio..." with 😄
00:02:21	 発信元: 052 Dali Georgobiani : Thank you very much, Hunter!!!
00:02:53	 発信元: *Tutor* Tatsuhiko OGAWA : 1.Video/Chat files of lectures are available from here   http://phits.jaea.go.jp/tutorial/rec.html
2.Please answer a survey by Friday from here https://forms.gle/htPwzn7eSbFMKoAr5
00:03:31	 発信元: 052 Dali Georgobiani : Thank you, and see you tomorrow!
00:03:32	 発信元: *Tutor* Hunter Ratliff : Replying to "is there any ipltmod..."

Yes! In the DCHAIN manual, see IPLOTMODE parameter. You can output xy, xz, yz, and any combination of the three.
00:03:41	 発信元: *Tutor* Hunter Ratliff : Reacted to "Thank you very much,..." with 👍
00:03:46	 発信元: 114 Gordeev Main : Thank you, goodbye!
00:03:48	 発信元: 129 -K. Butalag : Thanks. See tomorrow. Bye!
00:03:49	 発信元: 103 A.HUGUES : Thank you !
See you tomorrow !!
00:03:50	 発信元: 100 Herlan Setiawan : thank you..
00:03:51	 発信元: 240 Youngran Kwon : Thank you See you tomorrow
00:03:51	 発信元: 160 M. Adithya P : thank you
00:03:52	 発信元: 124-John Paul Bustillo : Thank you. 🙂
00:04:01	 発信元: 214 S.Schmidt : Thank you
00:04:03	 発信元: 107 Iaroslav Bielov : Thank you! See you tomorrow!
00:04:09	 発信元: 211_Somchai_Tan : Thank you very much
00:04:11	 発信元: 157 MOHD HARIS RIDZUAN OOI BIN ABDULLAH : Thank you!
00:04:29	 発信元: 047-Chalit Muanglay-(Sub) : Thank you very much
00:04:29	 発信元: 251Kayo : Thanks
00:04:30	 発信元: 043 A.Bekouale : Thanks for you