GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » EMC » EMC resolution
EMC resolution [message #15567] Fri, 11 October 2013 14:05 Go to next message
binsong is currently offline  binsong
Messages: 10
Registered: December 2011
occasional visitor
From: *in2p3.fr
hi,
i find some problems with the EMC energy resolution. Maybe there are some bugs.

I do the simulations for electron and photon with box generator with the latest trunk at both barrel region (80deg--100deg) and forward region (12deg--20deg).

i get the reconstructed energy(E_rec) by using GetEmcRawEnergy(),

then i plot (E_mc-E_rec)/E_mc for each simulation.

but i find the resolution sigma/E much smaller then the result in TDR:

electron_barrel: 1.18%
electron_forward: 1.58%
photon_barrel: 1.07%
photon_forward: 1.39%

and in TDR, the resolution for photon at 1GeV is more than 2%.

you can find all the four plots in attachment.

does someone know this problem?

Binsong

Re: EMC resolution [message #15570 is a reply to message #15567] Mon, 14 October 2013 14:47 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,

Some additional comments from my side on EMC resolution.

I have more or less the same numbers as Binsong and I tried to look what are main constributions to these numbers.

By the way I used apr13 release of pandaroot for these studies and did simulation for barrel in 30-130 degree range.

So first of all I made a comparison of resolution which includes and does not include digitization.

The resolution without digitization can be obtained if PndEmcMakeDigi class is used for digitization and in emc.par
UseDigiEffectiveSmearing:Int_t 0

In this case the resolution is only due to energy leakage below threshold. With single crystal threshold 3 MeV the resolution is 1.08 %

The resolution in my case is obtained from Novosibirsk fit
(root/pandaroot/trunk/macro/emc/dedicated/fit_resolution.C)

When digitization is included two additional constributions to resolution are photostatistics and electronics noise. With current parameters in emc.par/all.par photostatistics for case of barrel for 1 GeV photon should give 0.43 % (94 p.e. per MeV and 1.7 excess noise factor). 1.5 MeV noise gives 0.15% contribution accordingly.
Adding these three contrubution quadratically
sqrt(1.08^2+0.15^2+0.43^2)=1.17 %

Simulation of 50 k photons with digitization included gives 1.20 % resolution which is in agreement with simple estimation.

If we want to compare the results of simulation with prototype measurements the non-uniformity of ligh yield should be included in simulation and it will give a sizeble effect.

But otherwise it seems that electronics noise and photon statistics are treated reasonably well and I can only imagine that the energy leakage could be underestimated.

Dima
Re: EMC resolution [message #15573 is a reply to message #15570] Mon, 14 October 2013 17:29 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: SURFnet-CIDR-90-145-invalid
If I remember well, after Hossein update, we switched off the non uniformity because there were some problems, but I do not remember which ones.
Re: EMC resolution [message #15616 is a reply to message #15573] Wed, 30 October 2013 10:50 Go to previous messageGo to next message
binsong is currently offline  binsong
Messages: 10
Registered: December 2011
occasional visitor
From: *in2p3.fr
Hi,

I try to re-put the nonuniformity to 1.

But there is no difference for the resolution.

it is always the same, 1.15%

Binsong
Re: EMC resolution [message #15618 is a reply to message #15616] Wed, 30 October 2013 14:57 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
I notived Dima has done some modifications to have the old digitization. Could you please update emc folder and try again?
Re: EMC resolution [message #15619 is a reply to message #15567] Wed, 30 October 2013 15:23 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
Some additional results on emc resolution from my side.

First of all after Jifeng put a latest version of digitization by default, (based on code for EMC feature extraction developed at KVI) the energy resolution is different.

I did study for barrel 30-130 degree range, 10 k events, 1 GeV photons.

For new digitization sigma=1.6%
index.php?t=getfile&id=7612&private=0

For previous digitization sigma=1.2%
index.php?t=getfile&id=7613&private=0

There are two options for non-uniformity, first is based on measurements and implemented by Christian Hammann, sigma=2.6%
index.php?t=getfile&id=7614&private=0

And second one (linear non-uniformity) is based on fit to reproduce prototype data and implemented by Hossein Moeini, sigma=2.5%, so both options are very close.
index.php?t=getfile&id=7615&private=0
You can use this option with
PndEmcHitProducer* emcHitProd = new PndEmcHitProducer();
TString nonuniformityFile=gSystem->Getenv("VMCWORKDIR");
nonuniformityFile+="/macro/params/EmcDigiNoniformityPars2.root";
emcHitProd->SetNonuniformityFile(nonuniformityFile.Data());


So non-uniformity definitely affect the result as well as newer digitization.

But you can see if you use non-uniform response the peak position is shifted and calibrartion/energy correction is necessary. And shift is in different direction for two implementations of non-uniformity. Hossein provided correction for his case but it was for older digitization and newer digitization shift a peak position itself a little bit. So it should be redone.

So I plan to redo energy correction for new digitization with both non-uniformity options and I hope rather soon.

2.6% resolution for 1 GeV photon is what Christian Hammann quoted before as a result in agreement with prototype data.

So when energy correction will be ready the question can be closed.

Dima


  • Attachment: energy1.png
    (Size: 14.68KB, Downloaded 568 times)
  • Attachment: energy2.png
    (Size: 15.18KB, Downloaded 551 times)
  • Attachment: energy3.png
    (Size: 15.80KB, Downloaded 486 times)
  • Attachment: energy4.png
    (Size: 16.12KB, Downloaded 606 times)
Re: EMC resolution [message #15670 is a reply to message #15619] Wed, 06 November 2013 11:13 Go to previous messageGo to next message
Ronald Kunne is currently offline  Ronald Kunne
Messages: 32
Registered: October 2009
continuous participant
From: *in2p3.fr
index.php?t=getfile&id=7637&private=0Some remarks on the EMC resolution problem.

1) the present EMC PID is two years old, release 14213.

2) I think what should be available in PandaRoot ideally is two versions for the cluster energy Eraw et Ecorrected, where Eraw is the cluster energy as close to the value as we will measure it in the experiment, and Ecorrected a software corrected value to be used for gamma tracks and Binsong's electron correction

3) At present there are three values Eraw (non-uniformity off), an Eraw(non-uniformity on) et an Ecalibrated (which is probably wrong, if it uses Eraw). Of these: Eraw(non-uniformity on) represents the best the data as will be measured in the experiment

4) I checked the resolutions using the tuples made recently by Gosia. Eraw and Ecalibrated give resolutions that practically identical and both are too low.

Description of the plot added.
Upper raw: data using Eraw
Lower raw: data using Ecalibrated

1st column: pMC versus (Eemc-pMC)/pMC in %
2st column: profile plot of pMC versus (Eemc-pMC)/pMC in %
3st column: resolution obtained by fitting Gaussians to slices of the 2D histos.


Greetings,
Ronald Kunne
Re: EMC resolution [message #15671 is a reply to message #15567] Wed, 06 November 2013 12:35 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,

I have just update emc energy correction parameters for the case of non-uniformity switched on with the latest version of digitization.

I also switched on the use of non-uniformity by default in all.par

The files with correction are in /macro/params/
emc_correction_hist_gamma_2.root
emc_correction_hist_gamma_3.root
emc_correction_hist_gamma_4.root
emc_correction_par_gamma_2.root
emc_correction_par_gamma_3.root
emc_correction_par_gamma_4.root

You use correction like

PndEmcAbsClusterCalibrator * calibrator1= PndEmcClusterCalibrator::MakeEmcClusterCalibrator(1, 3);

And then obtain calibrated energy like
cluster_energy_calibrated=calibrator1->Energy(cluster);

And here first number stands for the method applied (1 - correction from histogram, 2 - correction from parametrization)

Second number stands for version and here
1 - previous version of digitization, no non-uniformity (deprecated)
2 - current version of digitization, no non-uniformity
3 - current version of digitization, with non-uniformity
(should be used by default)
4 - current version of digitization, with non-uniformity calculated by Hossein (linear non-uniformity 1.5%, set from the file /macro/params/EmcDigiNoniformityPars2.root)

By the way the macro which produces the last non-uniformity file is /macro/emc/dedicated/fill_nonuniformity_param.C


Here some plots which demonstrate how applying non-uniformity change reconstructed emc energy. I do not know if it's by accident or by intention of Christian Hammann who provided non-uniformity parameters from measured data in case of barrel EMC correction is almost not needed, i.e. 1 GeV energy peak is centred around 1 GeV. But for endcaps correction is still needed anyway.

For energies in range 0-10 GeV energy versus polar angle
index.php?t=getfile&id=7638&private=0

Reconstructed energy of 1 GeV photons
index.php?t=getfile&id=7639&private=0

Reconstructed pi0 invariant mass
index.php?t=getfile&id=7640&private=0
In this last case the corrected energy gives even slightly worse results.

In PndPidCorrelator.cxx
fEmcCalibrator= PndEmcClusterCalibrator::MakeEmcClusterCalibrator(2, 1);

should be modified to

fEmcCalibrator= PndEmcClusterCalibrator::MakeEmcClusterCalibrator(2, 3);

Dima

  • Attachment: e_vs_theta.png
    (Size: 38.03KB, Downloaded 624 times)
  • Attachment: energy_1GeV.png
    (Size: 17.17KB, Downloaded 589 times)
  • Attachment: mpi0.png
    (Size: 25.53KB, Downloaded 578 times)
Re: EMC resolution [message #15679 is a reply to message #15671] Wed, 06 November 2013 17:11 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
I have just updated the PndPidCorrelator.
Can somebody check the output of the simulation now, to be sure that everything is fine and coherent?
Thanks.
Re: EMC resolution [message #15720 is a reply to message #15679] Fri, 15 November 2013 15:10 Go to previous messageGo to next message
binsong is currently offline  binsong
Messages: 10
Registered: December 2011
occasional visitor
From: 134.158.196*
hi,

I have checked with the new correlator for electron and photon at 1 GeV/c at Barrel region.

The calibrate energy in EMC seems a little strange, I plot the (E_MC-E_Cal)/E_MC, there is a gap in the peak. (see the figures)

Binsong
Re: EMC resolution [message #15722 is a reply to message #15720] Fri, 15 November 2013 15:24 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 Binsong,

Could you make the same plot not with calibrated but with raw energy?
It would help to understand at what point the gap in the energy distribution appear.

Dima
Re: EMC resolution [message #15723 is a reply to message #15722] Fri, 15 November 2013 15:30 Go to previous messageGo to next message
binsong is currently offline  binsong
Messages: 10
Registered: December 2011
occasional visitor
From: 134.158.196*
yes, this is the same plots with raw energy.

Binsong
Re: EMC resolution [message #15725 is a reply to message #15567] Fri, 15 November 2013 15:39 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
One more clarification.
For which polar angle range in barrel are these plots,
the whole barrel or more narrow angle range?

Dima
Re: EMC resolution [message #15726 is a reply to message #15725] Fri, 15 November 2013 15:49 Go to previous messageGo to next message
binsong is currently offline  binsong
Messages: 10
Registered: December 2011
occasional visitor
From: 134.158.196*
the theta is from 30deg to 130deg, and phi is from 0 to 360.
Re: EMC resolution [message #15728 is a reply to message #15720] Fri, 15 November 2013 16:23 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
Are you selecting all the clusters or only the ones coming from the primary particles?
Re: EMC resolution [message #15729 is a reply to message #15728] Fri, 15 November 2013 16:30 Go to previous messageGo to next message
binsong is currently offline  binsong
Messages: 10
Registered: December 2011
occasional visitor
From: 134.158.196*
only the one from primary track
Re: EMC resolution [message #15731 is a reply to message #15729] Fri, 15 November 2013 16:35 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
Could you please do a 2D plot as a function of theta?
Re: EMC resolution [message #15871 is a reply to message #15567] Thu, 13 February 2014 15:16 Go to previous messageGo to next message
binsong is currently offline  binsong
Messages: 10
Registered: December 2011
occasional visitor
From: *in2p3.fr
Dear all,

now the problem is clear. it is due to the nonuniformity.

To obtain a realistic resolution (validated by the PROTO60 results),
simulations must be done using following parameters in all.par:

Use_nonuniformity:1
Use-Shaped_noise 0

The first flag needs to be set to 1 to take into account the observed
non linearity of the light yield
in the barrel EMC crystals as a function of the interaction depth.
The effect is an increase of the sigma of a Nivosibirsk fit from
1.6 (Use_nonuniformity:0 ) to 2.5 % (with Use_nonuniformity:1) for
photons at 1 GeV with a flat distribution in theta from 30 to
130 deg.

With the second flag set to 0, a more precise description of the
electronic noise is obtained, while faster but less accurate results are
obtained with the flag set to 1. The effect is an increase of the
sigma of a Novosibirsk fit from
2.2 (Use-Shaped_noise 1) to 2.5 % (Use-Shaped_noise 0) for a photon
at 1 GeV

Before end of november 2013, the standard parameters were
"Use_nonuniformity 0" and "Use-Shaped_noise 1", leading to a too
optimistic resolution, by about a factor 2. Now, the standard setting is
"Use_nonuniformity 1" and "Use-Shaped_noise 0")

The raw energy includes both effects of the non-uniformity and of the
electronic noise. The cluster calibrated energy is deduced using
calibration
parameters adjusted as a function of theta and energy to reproduce the
pi0 mass.
By default, the calibration parameter file is the one corresponding to
Use_nonuniformity:1 and "Use-Shaped_noise 0".

The two attached figures are the results for photon with the new version of pandaroot.
index.php?t=getfile&id=7696&private=0

[Updated on: Thu, 13 February 2014 16:54]

Report message to a moderator

Re: EMC resolution [message #15872 is a reply to message #15567] Thu, 13 February 2014 15:17 Go to previous message
binsong is currently offline  binsong
Messages: 10
Registered: December 2011
occasional visitor
From: *in2p3.fr
here is another figure for the cal energy

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

[Updated on: Thu, 13 February 2014 16:55]

Report message to a moderator

Previous Topic: [FIXED] About negative track id -2
Next Topic: An idea for EMC Error matix
Goto Forum:
  


Current Time: Tue Mar 19 08:50:54 CET 2024

Total time taken to generate the page: 0.00957 seconds