Home » PANDA » PandaRoot » EMC » Mass calculation from vector<PndEmcDigi*> in EMC
Mass calculation from vector<PndEmcDigi*> in EMC [message #8592] |
Tue, 19 May 2009 13:51 |
donghee
Messages: 385 Registered: January 2009 Location: Germnay
|
first-grade participant |
From: *kph.uni-mainz.de
|
|
Dear all,
I have a question about Mass function in PndEmcCluster using the list of PndEmcDigi.
In this Function, a list of digitization in ONE cluster is called, and all stored digi iterate to build a momentum.
If I correctly understand now, a cluster has few digitization according certain space points, where is only covered region in one cluster. I cannot clearly understand why two or even more digitizations are needed to calculate momentum, how do they involve in one cluster?
Could you explain about that? I just want to know the functionality of this digi list.
Thank you for your help in advance
Donghee Kang.
Quote: | Double_t PndEmcCluster::Mass() const
00216 {
00217 Double_t mass;
00218 TVector3 clusterMomentum(0,0,0);
00219 vector<PndEmcDigi*>::const_iterator digi_iter;
00220 TVector3 digiDirection;
00221 Double_t digiEnergy;
00222 for (digi_iter=fDigiList.begin();digi_iter!=fDigiList.end();++digi_iter)
00223 {
00224 digiDirection=(*digi_iter)->where().Unit();
00225 digiEnergy=(*digi_iter)->GetEnergy();
00226 clusterMomentum=clusterMomentum+digiDirection*digiEnergy;
00227 }
00228
00229 Double_t clEnergy=energy();
00230
00231 mass=sqrt(clEnergy*clEnergy-clusterMomentum.Mag2());
00232
00233 return mass;
00234 }
|
|
|
|
|
Re: Mass calculation from vector<PndEmcDigi*> in EMC [message #8599 is a reply to message #8592] |
Wed, 20 May 2009 08:53 |
donghee
Messages: 385 Registered: January 2009 Location: Germnay
|
first-grade participant |
From: *dip.t-dialin.net
|
|
Dear Bertram,
I can clearly understand your explanation for the calculation of shower mass.
A question to what comes next is the mechanical construction of EMC. I think that a cluster is composed of few crystals, which connects with Front-End-Electronics. If it is correct, the concept of mass calculation can be easily understood.
Could you give me how is the EMC cluster built in current planning setup, in the level of crystals?
Second issue is the inconsistancy of mass function with comparing energy.
If I get the mass info from cluster, we need to correct the calculated mass as performed energy, as well.
Since deposit energy or energy loss can calibrate with GetEnergyCorrected() function, one has to apply same calibration to the mass calculation. But I think that is not considered.
I guess it is not possible to do with digi list.
After constructing the mass, the cluster mass will be corrected again.
If digi list can provide such process, than we need two different mass function with corrected energy and without correction.
Finally, I would like to show you a plot for electron, photon and proton event.
This plot shows a mass distribution of all EMC cluster.
You can see a electronic noise near 0 due to high energetic 15 GeV forward direction proton, second bump is exactly photon, and last bump is for electron. Because I have checked already last two bump corresponding photon and electron with comparing true MC info.
Unfortunately all electron is not located in the expected mass region.
I plot the mass with the unit in GeV scale, but I assume that the unit must be MeV, otherwise mass(e) = 0.511 MeV of electron cannot archive.
This is the reason why we need a correction for mass also.
Best wishes,
Donghee Kang
[Updated on: Wed, 20 May 2009 08:56] Report message to a moderator
|
|
|
|
Re: Mass calculation from vector<PndEmcDigi*> in EMC [message #8601 is a reply to message #8599] |
Wed, 20 May 2009 11:58 |
Bertram Kopf
Messages: 110 Registered: March 2006
|
continuous participant |
From: *ep1.ruhr-uni-bochum.de
|
|
Dear Donghee,
thank you very much, Dima, for you explanation. Dima is right, it doesn't make sense to correct the invariant shower mass.
donghee wrote on Wed, 20 May 2009 08:53 |
Finally, I would like to show you a plot for electron, photon and proton event.
This plot shows a mass distribution of all EMC cluster.
You can see a electronic noise near 0 due to high energetic 15 GeV forward direction proton, second bump is exactly photon, and last bump is for electron. Because I have checked already last two bump corresponding photon and electron with comparing true MC info.
Unfortunately all electron is not located in the expected mass region.
I plot the mass with the unit in GeV scale, but I assume that the unit must be MeV, otherwise mass(e) = 0.511 MeV of electron cannot archive.
|
As Dima wrote, the invariant shower mass is a property which can be used to discriminate clusters from real photons and cluster which are created by two photons from pi0 decay. I think it doesn't make sense to distinguish photons, electrons and protons with a cut one the shower mass. As I mentioned it in one of my last postings: The first step is to do a matching between clusters and charged tracks which allows to distinguish between neutral and charged particles. After that one can start with the PID of charged particles where also the shower shape informations of the cluster are helpful properties.
Best regards,
Bertram.
|
|
|
|
Re: Mass calculation from vector<PndEmcDigi*> in EMC [message #8605 is a reply to message #8603] |
Wed, 20 May 2009 13:53 |
Bertram Kopf
Messages: 110 Registered: March 2006
|
continuous participant |
From: *ep1.ruhr-uni-bochum.de
|
|
Hi Donghee,
donghee wrote on Wed, 20 May 2009 12:18 |
I think that can be also used to separate electron and single photon.
|
I think, it is in not possible (at least not easy) to seperate electrons from single photons by using only shower shape informations. The shower will be produced in a very similar way via electromagnetic processes and also the complete energy will be deposited for both particles. Therefore the shower shapes look very similar. Hadrons instead lose only a certain fraction of the energy by ionization processes and thus the shape of the shower looks different.
donghee wrote on Wed, 20 May 2009 12:18 |
Concerning geometry of EMC!
I have learn, current setup of EMC designed with missing region at 142 to 149 degree for theta angle in backward endcap.
I have heard actual design should be slightly different, this missing acceptance will be covered in the future, but it is not introduced just in the MC simulation.
How about the situation for the theta region 142 to 149?
Do you have plan to modify some acceptance for EMC in MC simulation?
|
This is a question which could answer the relevant EMC hardware group. As far as I know, the design of the backward endcap is still not fixed. It's also not definitely clear how much space can be provided for the endcap. Therefore I think is not unrealistic to simulate with a gap between 142 to 149 degrees. Am I right?
Best regards,
Bertram.
|
|
|
Goto Forum:
Current Time: Fri Nov 08 22:18:35 CET 2024
Total time taken to generate the page: 0.00712 seconds
|