Home » PANDA » PandaRoot » EMC » EMC resolution
Re: EMC resolution [message #15671 is a reply to message #15567] |
Wed, 06 November 2013 12:35 |
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
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);
Dima
-
Attachment: e_vs_theta.png
(Size: 38.03KB, Downloaded 739 times)
-
Attachment: energy_1GeV.png
(Size: 17.17KB, Downloaded 695 times)
-
Attachment: mpi0.png
(Size: 25.53KB, Downloaded 680 times)
|
|
|
Goto Forum:
Current Time: Tue Nov 05 08:26:35 CET 2024
Total time taken to generate the page: 0.00807 seconds
|