Home » PANDA » PandaRoot » Tracking » Tracking: Kalman Task with STT,(electron hypo)
|
|
Re: Tracking: Kalman Task with STT,(electron hypo) [message #11371 is a reply to message #11370] |
Tue, 21 December 2010 12:11 |
Anonymous Poster
|
|
From: *adsl.alicedsl.de
|
|
Hi,
this is a very interesting point. I know from my studies that the GEANE track follower which is used when you perform the Kalman fit performs quite nicely on electrons, meaning that the performance should be better if you make the right particle hypothesis. So this looks like some kind of bug in our code. Maybe Lia could confirm whether the PID information is passed correctly yo GEANE in all places where this is needed in GeaneTrackRep?
Could you tell us what the exact macros are you are using to get to these results?
Cheers, Christian
|
|
|
|
|
|
|
|
Re: Tracking: Kalman Task with STT,(electron hypo) [message #11385 is a reply to message #11368] |
Wed, 22 December 2010 09:58 |
Ronald Kunne
Messages: 32 Registered: October 2009
|
continuous participant |
From: *in2p3.fr
|
|
I don't know how relevant this is, as the version I use dates from october, but I did a comparison between Kalman reconstruction using the muon hypothesis for all particles versus Kalman using the correct particle hypothesis. Only the reconstruction of electrons supposing they are electrons gives strange results for the resolution.
- Pandaroot trunk version 9969
- Geant3
- pgun generator with vertex (0,0,0)
- theta = 5-140°, phi = 0-360° momentum = 0.2-5 GeV/c
- one primary track per event, only primary tracks are taken into account in the analyses
- 10,000 events per particle
- MVD+TPC+GEM tracking detectors
- macros:
run_sim_tpccombi_pgun.C
run_sim_tpccombi_evtgen.C
run_digi_tpccombi.C
run_reco_tpccombi.C
Attached pdf file with five plots:
p1: sigma(p)/p = (p_MC-p_reco)/p_MC as a function of p (in GeV/c) for four bins in theta. This is mu reconstucted as mu and sigma(p) is the width if the gaussian fitted to the residues. It has to be noted that dp/p ought to be proportional to p, at least for momenta where multiple scattering is not important, and this is not the case.
p2: sigma(p)/p (e- reconstructed as mu-) Note that for the first angle the points scatter a bit, as there are not many events.
p3: sigma(p)/p (e- reconstructed as e-) These are the reconstructions that come out wrong
p4: sigma(p)/p (pi- reconstructed as mu-)
p5: sigma(p)/p (pi- reconstructed as pi-) The behaviour is comparable in both reconstructions.
-
Attachment: kalman.pdf
(Size: 105.99KB, Downloaded 292 times)
|
|
|
Re: Tracking: Kalman Task with STT,(electron hypo) [message #11390 is a reply to message #11368] |
Wed, 22 December 2010 18:22 |
Lia Lavezzi
Messages: 291 Registered: May 2007 Location: Torino
|
first-grade participant |
From: *pv.infn.it
|
|
Hi Mohammad and others,
here I attach the results of some tests made with the PndTools/MVD/macro macros (except for the selectEvents.cpp, do I have to compile it?).
The first two figures show the momentum distributions for a simulated electron with momentum in the range 1. - 1.5 GeV/c. It is reconstructed with kalman with electron and muon mass hypothesis respectively:
The electron reconstructed as electron distribution has a more gaussian shape, the electron reconstructed as a muon is more peaked, but the tail is more evident.
My first question is: did you use the kalman flag check to throw away badly reconstructed tracks?
The other two figures are the pull distributions: (p_MC - p_RECO)/sigma_from_kalman in the two cases again:
here you can see that the electron reconstructed as electron is almost a gaussian with mean = 0 (from the fit I get 0.087) and sigma = 1 (from the fit again I get 0.85). Sorry, I did the test with just 1000 events, so the statistics is quite low... The electron reconstructed as muon has a very bad shape: here it should be more evident that the right choice is the electron mass.
This happens even thought the bremsstrahlung is not yet treated in a complete way.
Now, from Ronald post we can also see that pions behave the same way when reconstructed with muon and pion hypothesis, but this is ok, since they are not distinguishable from dedx, so they are expected to behave like this.
What is missing is just a test with kaons and protons...
Bye bye and Buon Natale to everyone!
Lia.
|
|
|
|
|
|
|
Re: Tracking: Kalman Task with STT,(electron hypo) [message #11397 is a reply to message #11396] |
Thu, 23 December 2010 13:30 |
StefanoSpataro
Messages: 2736 Registered: June 2005 Location: Torino
|
first-grade participant |
From: *0-87-r.retail.telecomitalia.it
|
|
I still do not understand; let's forget about the tails, for the moment:
Muon: p_rec - p_mc narrower, p_rec - p_mv / sigma broader
EWlectrons: p_rec - p_mc broader, p_rec - p_mc / sigma narrower
->
Muons: sigma smaller
Electrons: sigma bigger
Muons have a better residulas, but small errors (why?)
Electrons have larger residuals, and larger errors.
I would suppose that errors for muons should be larger, and also the residuals, considering that this is the wrong hypothesis. Electrons should fit better the points, with a reduced residual and properly calculated errors. But this does not appear.
I continue to think that there is something wrong somewhere for electrons.
Let's remember that track parameters are used not only for kinematical fit, but also for pid and for correlation, where we do not need covariance matrix.
|
|
|
|
Re: Tracking: Kalman Task with STT,(electron hypo) [message #11404 is a reply to message #11400] |
Thu, 23 December 2010 21:54 |
Lia Lavezzi
Messages: 291 Registered: May 2007 Location: Torino
|
first-grade participant |
From: *47-151.net24.it
|
|
Hi Soeren,
sorry for the late reply. I tested the macros without the added Urban, i.e. setting the value of gcalpha = 1 in Geane.C.
In the reconstruction of the electron as a muon nothing changes, while the distribution and the pull of the electron reconstructed as an electron changes:
The momentum distribution is much more similar to the one reconstructed as a muon, but the pull distribution is not good.
It seems that the reason for the electron behaviour is indeed the Urban model correction, but the shape of the momentum distribution and the pull distribution are better with the Urban correction, though the peak is larger.
Answering also to Stefano:
Quote: | I still do not understand; let's forget about the tails, for the moment:
Muon: p_rec - p_mc narrower, p_rec - p_mv / sigma broader
EWlectrons: p_rec - p_mc broader, p_rec - p_mc / sigma narrower
->
Muons: sigma smaller
Electrons: sigma bigger
Muons have a better residulas, but small errors (why?)
Electrons have larger residuals, and larger errors.
I would suppose that errors for muons should be larger, and also the residuals, considering that this is the wrong hypothesis. Electrons should fit better the points, with a reduced residual and properly calculated errors. But this does not appear.
I continue to think that there is something wrong somewhere for electrons.
|
the errors are smaller for the muons because the bremsstrahlung contribution is not taken into account, while it is for the electron.
The point is that here we are considering the kalman fitted result, so a different error means also a different result when you compare the extrapolated error to the measured hit one. Then if you have different errors, you get different results after the filtering step and then different momentum distributions. The momentum itself does not come only from the extrapolation of the mean value.
To understand completely the GEANE contribution, I think we should consider the pull distribution just after an extrapolation, without the kalman procedure, maybe just extrapolating the tracks on a plane, taking as starting momentum and position the MC ones.
Regards,
Lia.
|
|
|
|
|
Re: Tracking: Kalman Task with STT,(electron hypo) [message #11409 is a reply to message #11368] |
Fri, 24 December 2010 14:14 |
Alberto Rotondi
Messages: 9 Registered: June 2007
|
occasional visitor |
From: 151.57.156*
|
|
Dear All,
perhaps I can help the discussion by recalling the modifications that we made in GEANE concerning the Energy loss. They are summarized in the enclosed report paper.pdf.
The original GEANE deals correctly the average energy loss for every particle. Electrons and positrons are tracked also taking correctly into account the bremsstrahlung process (see the enclosed report, table 2).
Therefore, no new implementations are necessary for the total energy loss. I recall that, when one studies a detector response, usually only the ionization energy loss is taken into account, whereas, while tracking, it is essential to consider the total energy loss, that is ionization plus bremsstrahlung, which is the main source of dE/dx for high energy electrons. This is considered with electrons, but not for muons (for energies below 10 GeV). Hence, to track an electron as a muon underestimates heavily the current particle energy. Note that the curves (pmc – pGeane)/pmc when electron is tracked as muon, have all the wrong mean and the peak is shifted.
Concerning the errors, we found in the original GEANE only the energy loss straggling by ionization in the Gaussian approximation (eq . 4 of the report). This approximation is absolutely insufficient for a correct error propagation!
Therefore (see the report), we implemented the energy loss by collision straggling with the Urban model and added the bremsstrahlung straggling to electron and positrons. All the formulae implemented by us can be found in the routine ERLAND of GEANE or in the enclosed report.
The bremsstrahlung straggling deals correctly the sigma, but, since the shape is pathological and not Gaussian (see fig.5 of the report), the Kalman algorithm gives biased results. This is well known, and the remedy is to parametrize the bremsstrahlung straggling with a sum of Gaussian, to perform multiple tracking and to weight the results with the area of the Gaussians. The method is called Gaussian Sum Filter (GSF). At present this method is missing in our tracking. However, also with the actual implementations, the pulls of electrons are not so bad.
In the result reported in this discussion I see that the pulls for electrons and the other particles are acceptable when Urban fluctuations by ionization and the present version of bremsstrahlung straggling are on.
When an electron is tracked as a muon, the energy loss is wrong and the tracking error is underestimated. In the Kalman filter, the result is shifted more toward the tracking values than in the case of electrons. The errors used in Kalman is wrong and too small in any case. The results show that the pulls are bad (and this is expected, since the errors are wrong because the bremsstrahlung straggling is missing), and that the absolute differences are sharper. The mean is shifted (and I understand why) and the width is smaller (this is not completely clear to me).
I checked with Lia the code during these days and we did not find bugs.
Merry Christmas and Happy new year to all
Alberto
-
Attachment: paper.pdf
(Size: 140.26KB, Downloaded 354 times)
[Updated on: Fri, 24 December 2010 14:26] Report message to a moderator
|
|
|
Goto Forum:
Current Time: Sat Nov 09 11:39:17 CET 2024
Total time taken to generate the page: 0.01010 seconds
|