GRIFFINCollaboration / detectorSimulations_v10

Geant4 version 10 of the simulation code for the GRIFFIN array and it's suite of ancillary detection systems.
MIT License
9 stars 14 forks source link

Setting energy threshold for DESCANT detectors #88

Closed akanksha-sing closed 5 months ago

akanksha-sing commented 5 months ago

Hello,

I am plotting the neutron time-of-flight for 11Li decay, and I see a very long tail in the plot resulting from low-energy neutrons, which may be below the DESCANT detection threshold. In order to cut those and get a more realistic spectrum, I am trying to implement a threshold on the energy deposited in DESCANT detectors. I would like to set a threshold on the deposited energy and if the energy deposited in one of the DESCANT detectors is less than this, then it should not be counted. Is there a way I can do this? Any help is greatly appreciated.

VinzenzBildstein commented 5 months ago

If you use NTuple2EventTree (and I think this also works for NTuple), you can set a threshold for each detector using Descant.colour.number.Threshold.keV Descant.colour.number.ThresholdWidth.keV where colour is the colour of the detector, and number its number. This applies an error function centred at the threshold with the given width to the energies.

VinzenzBildstein commented 5 months ago

@akanksha-sing has this solved your issue? If so, please close it.

akanksha-sing commented 5 months ago

@VinzenzBildstein Yes, it works for me. Thanks a lot!