[gate-users] Constant Energy resolution for solid-state detectors

[Eric Vandervoort]

Hello,
The most simple way might be to change the source code. If you are using
the simple blurring digitizer module, i.e.,

/gate/digitizer/Singles/blurring
/gate/digitizer/Singles/blurring/setResolution X
/gate/digitizer/Singles/blurring/setEnergyOfReference Y keV

You could change the following scrap of code in GateBlurring.cc:

m_coeff = m_resolution * sqrt(m_eref);
GatePulse* outputPulse = new GatePulse(*inputPulse);        
outputPulse->SetEnergy(G4RandGauss::shoot(inputPulse->GetEnergy(),m_coeff*sqrt(inputPulse->GetEnergy())/2.35));
outputPulseList.push_back(outputPulse);

Just change it to:

m_coeff = m_resolution;
GatePulse* outputPulse = new GatePulse(*inputPulse);        
outputPulse->SetEnergy(G4RandGauss::shoot(inputPulse->GetEnergy(),m_coeff*sqrt(inputPulse->GetEnergy())/2.35));
outputPulseList.push_back(outputPulse);

That will give you a constant FWHM for the energy resolution of X keV.
Cheers,
Eric

> Hi all,
>
> Since the electronic noise contribution is dominant over the uncertainty
> in
> the number of electron hole pairs generated from each interaction in
> solid-state detectors (for examples, HPGe, Si, CZT , etc.),  a constant
> energy resolution at all energy range deposited in the detector volume
> should be required to simulate blurring of the energy spectrum instead of
> using FWHM = FWHMo x sqrt(Eo/E) for scintillation detectors (LSO, BGO,
> NaI,
> etc.).
> If you know how to use constant blurring for all energy range, please let
> me
> know that.
>
> Thank you for your help in advance.
> Best,
> Sang-June Park
> _______________________________________________
> gate-users mailing list
> gate-users at lphe1pet1.epfl.ch
> http://lphe1pet1.epfl.ch/mailman/listinfo/gate-users
>