Message: Energy loss seems to be overestimated. Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Energy loss seems to be overestimated. 

Forum: Electromagnetic Processes
Date: 22 Jan, 2015
From: Genie Jhang <Genie Jhang>

Hello.

I’m currently writing a simulation package for an experiment.
Prior to using it, we’re checking if the energy loss of proton and alpha is the same as the other simulation.
We’re comparing GEANT4 energy loss data with TRIM/LISE++ data.

Attached links are proton and alpha energy loss plot with TRIM/LISE++ result.

   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2015/01/22/15.07-13310-0atGMsC.png
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2015/01/22/15.07-54093-KDTlfHd.png

As you can see, in relatively low energy, energy loss curves doesn’t match each other.
We couldn’t find what the problem is.

Currently, we’re writing our code on top of FairROOT package.
So, GEANT4 is steered over G4VMC interface.
Cuts used for running simulation is following.

If you have any clue or comments or recommendation, please don’t hesitate to mail us.
Anything could be helpful

Thank you in advance.

Best regards,
Genie

==================================================================

g4config.in:
/mcPhysics/rangeCuts 0.001 mm

/cuts/setLowEdge 10 eV
/cuts/setHighEdge 1 MeV

g4Config.C:
 TG4RunConfiguration* runConfiguration = new TG4RunConfiguration("geomRoot",
"QGSP_BERT_HP_EMY", "stepLimiter+specialCuts+specialControls”);
TGeant4 -> SetMaxNStep(1.E7);

SetCuts.C
gMC->SetProcess("PAIR",1); /** pair production **/
gMC->SetProcess("COMP",1); /** Compton scattering **/
gMC->SetProcess("PHOT",1); /** photo electric effect **/
gMC->SetProcess("PFIS",1); /** photofission **/
gMC->SetProcess("DRAY",1); /** delta-ray **/
gMC->SetProcess("ANNI",1); /** annihilation **/
gMC->SetProcess("BREM",1); /** bremsstrahlung **/
gMC->SetProcess("HADR",1); /** hadronic process **/
gMC->SetProcess("MUNU",1); /** muon nuclear interaction **/
gMC->SetProcess("DCAY",1); /** decay **/
gMC->SetProcess("LOSS",1); /** energy loss **/
gMC->SetProcess("MULS",1); /** multiple scattering **/

Double_t cutC = 1.0E-5;
//Double_t cutB = 1.E4;
Double_t cutB = -1;
Double_t cut1 = 1.0E-8; // [GeV]
Double_t cutp = 0.01; // [GeV]
Double_t tofmax = 1.E10; // seconds

gMC->SetCut("CUTGAM",cutC); /** gammas (GeV)*/
gMC->SetCut("CUTELE",cut1); /** electrons (GeV)*/
gMC->SetCut("CUTNEU",cutC); /** neutral hadrons (GeV)*/
gMC->SetCut("CUTHAD",cutC); /** charged hadrons (GeV)*/
gMC->SetCut("CUTMUO",cutC); /** muons (GeV)*/
gMC->SetCut("BCUTE",cutB); /** electron bremsstrahlung (GeV)*/
gMC->SetCut("BCUTM",cutB); /** muon and hadron bremsstrahlung(GeV)*/
gMC->SetCut("DCUTE",cutB); /** delta-rays by electrons (GeV)*/
gMC->SetCut("DCUTM",cutB); /** delta-rays by muons (GeV)*/
gMC->SetCut("PPCUTM",cutp); /** direct pair production by muons (GeV)*/
gMC->SetCut("TOFMAX",tofmax); /**time of flight cut in seconds*/

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

1 Idea: Re: Energy loss seems to be overestimated.   (Vladimir Ivanchenko - 27 Jan, 2015)
(_ Question: Re: Energy loss seems to be overestimated.   (Genie Jhang - 27 Jan, 2015)
(_ Idea: Re: Energy loss seems to be overestimated.   (Vladimir Ivanchenko - 03 Mar, 2015)
(_ Feedback: Re: Energy loss seems to be overestimated.   (Genie Jhang - 21 Mar, 2015)
 Add Message Add Message
to: "Energy loss seems to be overestimated."

 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 ]