Message: Re: Question about flux scorer Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Question about flux scorer 

Forum: Hits, Digitization and Pileup
Re: None Question about flux scorer (Haori Yang)
Date: 21 Apr, 2011
From: Tsukasa Aso <Tsukasa Aso>

Hi,

  First of all, we found trouble for a combination of cellflux scorer
with cylinder shape. If you apply R segmentation, the value is not correct.
This is caused by wrong calculation of cubic volume inside scorer. 
The bug will be fixed in next patch. 

  Concerning to getting value with unit, your code is good.
Quantity inside scorer is maintained in default unit of Geant4.
(The GetUnit() and GetUnitValue() are prepared for the purpose of 
command-based scorers.) 

 Basically, you will specify the unit as,

  flux = theRun->GetTotalFlux(i)/(1./cm2);   // same as percm2 (/cm2)
  flux = theRun->GetTotalFlux(i)/(1./mm2);   // same as permm2 (/mm2).

 In case of output file in command-based scorer, each value
is calculated by "value/GetUnitValue()".
This should give same result with your code. 

 But, according to your post, do you mean the value in command-based scoring 
is 100 times larger than that of your expectation?

> 
> for the command-based method, the output file looks like this:
> 
> ****************************
> 
> # mesh name: boxMesh_1
> # primitive scorer name: cFlux
> # iZ, iPHI, iR, value [percm2]
> 
> 0,0,0,2426.526585771724 ****************************
> 
> I had to multiply the 'value' by 0.01 to get the expected flux in unit
> of percm2
> 

Tsukasa

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: Question about flux scorer   (Haori Yang - 21 Apr, 2011)
2 None: Re: Question about flux scorer   (Vito Di Benedetto - 30 Apr, 2011)
 Add Message Add Message
to: "Re: Question about flux scorer"

 Subscribe Subscribe

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


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]