GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » EMC » Energy correction for gammas
Energy correction for gammas [message #7578] Fri, 28 November 2008 18:22 Go 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 all,

I made an attempt to implement the energy correction in emc just porting correction from Panda Emc in Babar framework.

The results are not very promising.

I simply have taken parametrisation of energy correction from EmcCalibToo/EmcPhotonSimpleCalib.cc and implemented it as a method of PndEmcBump::GetEnergyCorrected().

The results are the following.

For the barrel emc, gamma with energy 1 GeV emitted in theta [25;135] degree, phi [0;360] degree.

The following plot is obtained with macro\emc\reco_analys4.C.
The top plot is without correction and bottom with correction.

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

So the correction is too large. It seems than that correction for 3 percent energy loss would be closer to generated energy.

For the shaslyk emc, gamma with energy 1 GeV emmited in theta [0;5] degree, phi [0;360] degree.

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

And correction is too small.

And I am really sorry but I do not have time at the moment to investigate the discrepancy (i.e. why this parametrisation of the energy correction works in Babar framework but not in pandaroot).

And one general remark how to implement energy correction in the proper way. Of course the EmcBump is not the proper place to apply energy correction. It should be applied at least after Cluster-Track matching to reject clusters produced by the charged particles.

But it seems than it cannot be applied in exactly the same way as in Babar framework, since as I can see there it is strongly coupled with PID. At least correction is applied to the object which is produced after identification of the Bump/Cluster as Neutral Candidate. And in pandaroot the PID procedure is foreseen to be different and I do not see now how to couple energy correction with it at the moment. So it's a point for discussion.

Best regards,
Dima
Re: Energy correction for gammas [message #7581 is a reply to message #7578] Sun, 30 November 2008 20:55 Go to previous messageGo to next message
Johan Messchendorp is currently offline  Johan Messchendorp
Messages: 693
Registered: April 2007
Location: University of Groningen
first-grade participant

From: *xs4all.nl
Hi,

I ran this weekend simulations with more than two million photon events (boxgenerator, G3&G4) with only the EMC as geometry using the Nov08 release of PandaRoot. I looked at the averaged deposited energy (note: cluster energy!) with respect to the incident photon energy for bins over the full range in polar angle (0-180 degrees, 62 bins) and for an incident energy interval from 0-15 GeV in 34 bins (binning inhomogeously chosen according to variations in the response). The results, e.g. mean cluster energy/incident photon energy, for G3 are depicted in the figure. In the attached root file, you can find the corresponding 2D histogram together with a simulation using G4 and simulations for the Slashlyk calorimeter. In principle, one could use this file to correct the data accordingly. If there is a need for other information, let me know.

Kind wishes,

Johan.

[Updated on: Wed, 03 December 2008 00:21]

Report message to a moderator

Re: Energy correction for gammas [message #7582 is a reply to message #7578] Sun, 30 November 2008 22:20 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: *mobile.playmobile.pl
Hi all,

I have an update to my previuos post.

The previous plots were shown for G3 simulation, but in case of G4 simulation the energy correction taken from Babar framework works better.

The barrel emc - 25-135 degree, 1 GeV gamma.

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

At least the mean value of corrected histogramm is at 1 GeV.

And shashlyk calorimeter, 0-5 degree, 1 GeV.

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

The match is not perfect but much closer than in G3 case.

Best regards,
Dima
Re: Energy correction for gammas [message #7585 is a reply to message #7582] Sun, 30 November 2008 23:07 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *dip.t-dialin.net
As I have already told to Dima,
I think that the Shashlick, which is made of different thin layers of lead+scintillator, suffers from the same problem of MVD in geant3, due to the straggling option and then delta electrons below geant cuts which are absorbed by the material.
Re: Energy correction for gammas [message #7586 is a reply to message #7585] Mon, 01 December 2008 00:36 Go to previous messageGo to next message
Johan Messchendorp is currently offline  Johan Messchendorp
Messages: 693
Registered: April 2007
Location: University of Groningen
first-grade participant

From: *xs4all.nl
Hi,

something we need to discuss and settle on on how to proceed here...(as we already planning to do for a long time) What about the upcoming collaboration meeting?

Johan.
Re: Energy correction for gammas [message #7589 is a reply to message #7578] Mon, 01 December 2008 13:39 Go to previous messageGo to next message
Bertram Kopf is currently offline  Bertram Kopf
Messages: 110
Registered: March 2006
continuous participant
From: *ep1.ruhr-uni-bochum.de
Hi Dima,

thanks for checking the emc correction in the PandaRoot software. The results - at least for the TS EMC with VM G4 - look reasonable to me. As far as I know Jan corrected the energies to the mean value and this seems to match more or less well. The correction for the shasklyk calorimeter cannot be reproduced perfectly but I think one can also use it as a work around for the PB gamma gamma studies.
This would require that Irina uses VM G4 as simulation engine together with the correction from the BaBar like software.

I think the following general things have to be done in the near future:

1. the VM G3 results (especially for the shashlyk calorimeter) look very strange and - in case that one would also make use of VM G3 - it is necessary to fix it.

2. To implement a propper calibration. A quick solution would be to calibrate it with mc truth informations. For a long term project one should get rid of the mc truth calibration and to replace it by a more realistic one, e.g. based on DPM events (pi0-, eta- calibration and maybe also calibration with electrons).

3. In case that different simulation engines should be supported it is necessary to calibrate the EMC seperately for all the supported engines (VM G3, VM G4, fluka, etc.).

Best regards,
Bertram.
Re: Energy correction for gammas [message #7597 is a reply to message #7578] Wed, 03 December 2008 00:36 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: 89.108.200*
No Message Body
Re: Energy correction for gammas [message #7598 is a reply to message #7578] Wed, 03 December 2008 12:38 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 all,

I have just made the small modification in the parameter file /macro/params/emc.par (svn revision 4091).

The threshold for shashlyk was there by default 3 MeV. I changed it to 8 MeV, to be consistent with what is written in Physics book and correspondingly is used in Babar framework.

Best regards,
Dima
Re: Energy correction for gammas [message #7599 is a reply to message #7598] Wed, 03 December 2008 13:20 Go to previous message
Aleksandra Biegun is currently offline  Aleksandra Biegun
Messages: 64
Registered: May 2007
Location: Groningen
continuous participant
From: *KVI.nl
Hi Dima,

could you please change the map version (MapperVersion:i 1)
in that file:

/macro/params/emc.par

from 1 to 2, since now the new version of the Forward End-Cap
is used, see: macro/emc/sim_emc.C

Emc->SetGeometryFileNameDouble("emc_module1245.dat","emc_module3new.root "); // if you want to use new geometry for FwEndCap

There was version 2 some time ago, but now,
version 1 appeared back.


Thanks a lot!

Ola.


Aleksandra Biegun
University of Groningen/KVI
Zernikelaan 25, 9747 AA Groningen
tel. +31 50 363 3630
fax. +31 50 363 4003
Previous Topic: error matrix calculation for EmcCluster
Next Topic: Storage of points, hits, digis, etc...
Goto Forum:
  


Current Time: Thu Mar 28 14:42:10 CET 2024

Total time taken to generate the page: 0.00972 seconds