GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » EMC » Cut in energy for EMC hit and digi
Cut in energy for EMC hit and digi [message #17665] Wed, 17 December 2014 21:55 Go to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *ip71.fastwebnet.it
Dear EMC experts,
I noticed running the standard simulation macros the following message coming from the Emc Hit producer:

-I- PndEmcHitProducer INITIALIZATION *********************
-I- PndEmcHitProducer: Using nonuniform lightoutput HitProducer has EnergyHitThreshold of 0.000001 GeV and Use_nonuniformity 1
-I- PndEmcHitProducer: Intialization successful


Indeed checking the emc params:

[PndEmcDigiPar]
EnergyHitThreshold:Double_t 0.000001


The energy threshold for each EMC hit is 1 KeV... Is it maybe a too small cut? I mean, I suppose the minimum energy for emc is at the level of MeV and not KeV.
Indeed later in the param file:

EnergyDigiThreshold:Double_t 2.0e-3


Then the minimum energy is 2 MeV. This low cut for energy in the hit could have effects on the data size (but I did not try to raise the cut and see the effect). Were there maybe particular reasons to use such low cut?

Speaking about minimum digi energy, I read also in the digi params :

Incoherent_elec_noise_width_GeV_APD:Double_t 1.5e-3     ### optimized based on the results of proto60 experiment
Incoherent_elec_noise_width_GeV_VPT:Double_t 1.5e-3


If the sigma of the elec noise is 1.5MeV, I would say that a reasonable energy cut could be 3 sigma -> 4.5 MeV and not 3 MeV, maybe with 3 sigma the cut is too close to the electronic noise... but I am not expert of emc electronics, these are just thought which came into my mind.
Re: Cut in energy for EMC hit and digi [message #17666 is a reply to message #17665] Thu, 18 December 2014 11:06 Go to previous messageGo to next message
Dima Melnychuk is currently offline  Dima Melnychuk
Messages: 213
Registered: April 2004
Location: National Centre for Nucle...
first-grade participant
From: *fuw.edu.pl
Hi Stefano,
I suppose that EnergyHitThreshold could be safely changed to the same value as EnergyDigiThreshold and indeed it can reduce data size but probably not so much.

The reason for two separate cuts is historical.
Initially where only PndEmcHitProducer was implemented the EnergyHitThreshold cut played the role of the threshold set for EMC electronics. But later when digitization classes were implemented this electronic threshold was moved to digi classes but separate parameter was introduced EnergyDigiThreshold.

So the low value of EnergyHitThreshold is effectively 0 and from PndEmcHit you can study the distribution of energy deposited in crystals even for low energy close or below the threshold and changing one parameter EnergyDigiThreshold you can see how energy resolution is affected by this threshold.
And the role of EnergyHitThreshold could be only the data size reduction. And if you want to reduce the data size by default you can increase it as high as EnergyDigiThreshold.

Concerning the value of the threshold it should be indeed at the level of 3 sigma of electronics noise for the whole calorimeter, which will give 1 per 1000 EmcDigi from electronics noise. But so far with prototype measurements where the number of crystals is lower the threshold was not set at 3 sigma level but was optimized by hand and digitization parameters in simulation are based more or less on prototype measurements. So I would leave in simulation it as it is now and at the end it will be modified according to values used in real Panda experiment.

Dima


Re: Cut in energy for EMC hit and digi [message #17667 is a reply to message #17666] Thu, 18 December 2014 12:02 Go to previous message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
I plotted together the hit energy (red) and the digi energy (in blue) in 10k events of psi(3686)->j/psipi+pi-:

index.php?t=getfile&id=8219&private=0

and now hit vs digi:

index.php?t=getfile&id=8220&private=0

Most probably a cut of 2 MeV could be enough to remove a lot of noise. What do you think?
Previous Topic: Correlation with charged tracks
Next Topic: MC Truth with EMC
Goto Forum:
  


Current Time: Fri Mar 29 13:06:56 CET 2024

Total time taken to generate the page: 0.00581 seconds