Hello this is my G4 application about a 45GeV positron beam hitting two circular targets
https://we.tl/t-JFgCYJ2LwB. My targets have diameter d=1.5 cm and thickness 3.00 mm. The first one is centered in (0,0,0) the second one in (0,0,2.3).
In the application there is the Scoring mesh function defined in run1.mac to save data about released energy in the two targets in two dump file.
This is the code for the target 1
/score/create/cylinderMesh boxMesh_1
/score/mesh/cylinderSize 0.75 0.15 cm
/score/mesh/translate/xyz 0.0 0.0 0.0 cm
/score/mesh/nBin 30 30 30
/score/quantity/energyDeposit eDep
/score/quantity/nOfStep nOfStepGamma
/score/filter/particle gammaFilter gamma
/score/quantity/nOfStep nOfStepEMinus
/score/filter/particle eMinusFilter e-
/score/filter/particle muMinusFilter mu-
/score/quantity/nOfStep nOfStepEPlus
/score/filter/particle ePlusFilter e+
/score/filter/particle muPlusFilter mu+
/score/close
/score/list
and this is fo the target 2
/score/create/cylinderMesh boxMesh_2
/score/mesh/cylinderSize 0.75 0.15 cm
/score/mesh/translate/xyz 0. 0. 2.3 cm
/score/mesh/nBin 30 30 30
/score/quantity/energyDeposit eDep
/score/quantity/nOfStep nOfStepGamma
/score/filter/particle gammaFilter gamma
/score/quantity/nOfStep nOfStepEMinus
/score/filter/particle eMinusFilter e-
/score/filter/particle muMinusFilter mu-
/score/quantity/nOfStep nOfStepEPlus
/score/filter/particle ePlusFilter e+
/score/filter/particle muPlusFilter mu+
/score/close
/score/list
You see, I setted
/score/mesh/cylinderSize 0.75 0.15 cm
/score/mesh/translate/xyz 0.0 0.0 0.0 cm
and
/score/mesh/cylinderSize 0.75 0.15 cm
/score/mesh/translate/xyz 0. 0. 2.3 cm
because I want the scoring mesh in my two targets
Here the dump file of the two targets
boxMesh_1.txt (1.8 MB)
boxMesh_2.txt (2.2 MB)
but it looks like that results are wrong.
Because
-
the z axis should be in the range [-1.5mm;1.5mm] for the first target and [21.5mm;24.5mm] for the second target, instead in the dump files iZ values are in the range [0mm;30mm]
-
The phi axis should be in the range [0;2pi] (i.e. the angle of a full circle) instead it is in the range [0;30]
-
The R axis should be in the range [-0.75cm;0.75cm] (because the diameter is 1.5 cm instead it is in the range [0;30]
-
Mostly of realeased energy should be at center of the two targets, instead this is for example the plot of
# iZ, iPHI, iR, total(value) [MeV],
for the first target
I got it by using this root macro
So, I can’t understand the results…they don’t look like right…
Please, can someone check how to fix the results?
Thank you