Message: Re: How to add another sensitive detector type to the command-based scoring system? Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Note Re: How to add another sensitive detector type to the command-based scoring system? 

Forum: Biasing and Scoring
Re: Question How to add another sensitive detector type to the command-based scoring system? (Tim Meehan)
Date: 14 Sep, 2017
From: Tim Meehan <Tim Meehan>

I figured out how to add the command, but the problem I have now is that when I dump the the quantity to file, everything is in the (0, 0) location. I tried follow the form of the existing primitive scorer classes as closely as I could.

To summarize what I did, I started with a copy of the existing source tree in my home directory, and added the new primitive scorer classes (a "regular" PS class, and a "3D" PS class) to  “~/geant4.10.03.p02_working/source/digits_hits/scorer”. I also needed to modify the “sources.cmake” file there as well.

Next, I needed to make the interactive system aware of the new scorer. I had to modify “G4ScoreQuantityMessengerQCmd.cc” (which is the implementation of the function G4ScoreQuantityMessenger::QuantityCommands — I simply added my new command to the end of the function), and “G4ScoreQuantityMessenger.cc”.

After that, I was able to compile everything without error, ran “make” and then “make install” and then loaded up my problem. I was able to get the quantity loaded with:
/score/quantity/remDose rem_dose

I ran some particles through it, and then dumped it to file with:
/score/dumpQuantityToFile dose_slab rem_dose dosemap_rem.csv

However, the whole file has zeros in it, except for the first line:
# mesh name: dose_slab
# primitive scorer name: rem_dose
# iX, iY, iZ, total(value) [NoUnit], total(val^2), entry
0,0,0,35421788.46339133,6509346278844.741,10013
0,1,0,0,0,0
0,2,0,0,0,0
0,3,0,0,0,0
0,4,0,0,0,0
0,5,0,0,0,0

Any ideas out there?

 Add Message Add Message
to: "Re: How to add another sensitive detector type to the command-based scoring system?"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews