GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » EMC » EMC resolution
Re: EMC resolution [message #15671 is a reply to message #15567] Wed, 06 November 2013 12:35 Go to previous messageGo to previous message
Dima Melnychuk is currently offline  Dima Melnychuk
Messages: 213
Registered: April 2004
Location: National Centre for Nucle...
first-grade participant
From: *

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/

You use correction like

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

And then obtain calibrated energy like

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

Reconstructed energy of 1 GeV photons

Reconstructed pi0 invariant mass
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);


  • Attachment: e_vs_theta.png
    (Size: 38.03KB, Downloaded 650 times)
  • Attachment: energy_1GeV.png
    (Size: 17.17KB, Downloaded 612 times)
  • Attachment: mpi0.png
    (Size: 25.53KB, Downloaded 608 times)
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: [FIXED] About negative track id -2
Next Topic: An idea for EMC Error matix
Goto Forum:

Current Time: Tue Jun 25 12:05:14 CEST 2024

Total time taken to generate the page: 0.01724 seconds