GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Particle Identification » MisID vs Impurity
MisID vs Impurity [message #14221] Wed, 21 November 2012 15:48 Go to next message
Klaus Götzen is currently offline  Klaus Götzen
Messages: 293
Registered: June 2006
Location: GSI
first-grade participant
From: *gsi.de
Dear all,


since my mic died during the EVO meeting yesterday, I'll try to explain what I wanted to say concerning misID and impurity.

The issue I wanted to point out is, that the impurity is a quantity which is not independent of the fluxes, whereas the misID is. The misID is the fraction of false positive identified particles of a certain species. E.g. the pion-misID of a kaon selector is

misID(pi|K) = #selected pi / #total pi

Obviously this quantity is flux independent, since the flux would go into nominator and denominator. Please note, that there is not information about the selected number of kaons in this quantity, although it is a property of the kaon selector. It's just the probability for another particle species to be accepted by the kaon selector.

On the other hand the impurity is defined as

impurity(pi|K) = 1 - purity = 1 - #sel. K / (#sel. pi + #sel. K) = #sel. pi / (#sel K + #sel pi)

First of all, these two things are different quantities. The other issue is, that a change in relative fluxes would change the impurity (as well as the purity of course), since (with F being the relative pi/K flux factor change in the upper equation)

impurity = F·#sel pi / (#sel. K + F·#sel pi) != #sel. pi / (#sel K + #sel pi)


I think the quality measure of a particular selection algorithm should not depend on the current physical environment (like relative fluxes).

Of course I'm open for discussions in that respect.


Cheers,
Klaus


Re: MisID vs Impurity [message #14224 is a reply to message #14221] Thu, 22 November 2012 10:17 Go to previous messageGo to next message
Malgorzata Gumberidze is currently offline  Malgorzata Gumberidze
Messages: 98
Registered: June 2005
continuous participant
From: *gsi.de
Hi,

Unfortunatly i have missed full discussion after i have presented my slides ... In principle what i do i think is
correct, but i simple call it differently. Reading entry of
Klaus i realized, that what i call impurities, he is calling
mis-identification.

Just to be sure I recall what i do. I run simulation for of the
electrons and pions for example and then what i do:

X is some value of the cut on the PID for given particle to be
in this case electron.

electron efficiency:

ele_eff = PID_{e} > X / e_all

and than to study what i call impurity (but probably should be called mis-identification):

pi_imp = PID_{e} > X / pi_all

PID_{e} - i call here probablity of being an electron for given particle. In this case for pion.

greetings,
gosia
Re: MisID vs Impurity [message #14227 is a reply to message #14224] Thu, 22 November 2012 16:48 Go to previous messageGo to next message
Klaus Götzen is currently offline  Klaus Götzen
Messages: 293
Registered: June 2006
Location: GSI
first-grade participant
From: 140.181.9*
Dear Gosia,


thanks for pointing that out! Then your impurity is indeed the same I refer to as 'mis-ID'. And thank you for studying the PID performance in such detail, showing possible bugs in the mechanism!

Why sometimes the mis-ID levels get higher when adding more information I nevertheless don't understand.

I think something is going wrong with the computation and normalization of the PID values.


Best regards,
Klaus
Re: MisID vs Impurity [message #14228 is a reply to message #14227] Thu, 22 November 2012 17:21 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: 2.235.190*
Hi Klaus,
I believe that if you add more detectors, the meaning of the cut in the value we call "probability" has a different meaning, then cutting P>90% for EMC is different from P>90% in EMC+STT. If you compare the plots, once you add more detectors your efficiency with the same cut increases -> this means that you have more signal candidates in your selection -> your selection is looser, the misId increases.
In order to crosscheck really the numbers, one should select two algorithms cutting a value which provides the same efficiency. In such case one can compare purity and understand what is good and what is not good. The ROC curves could be helpful,maybe for the moment just integrated in separated momentum ranges.
Re: MisID vs Impurity [message #14231 is a reply to message #14224] Fri, 23 November 2012 17:29 Go to previous messageGo to next message
donghee is currently offline  donghee
Messages: 385
Registered: January 2009
Location: Germnay
first-grade participant
From: *kph.uni-mainz.de
Dear Gosia,

First of all, thank you for your great effort.
You showed to me the good guide line for PID study.

But as usual, I have still few question to better understand your definition of impurity.
pi_imp = PID_{e} > X / pi_all
I can understand about the nominator, that is reconstructed true pion, which is tested MC true PID matching.

I'm wondering about the PID_{e} > X.
Is it "selected true pion after doing MC PID match and requiring PID probability"
or
"any kind of tracks just passed given probability condition"?
I am not clear for this which one have to be applied to see whatever impurity or misID.

concerning single and more tracks
You have selected only one reconstructed particle, which has a closest momentum value to MC one.
We can have usually more than one track after reconstruction even from one event.
That means, all values of impurity shown in your categories are too ideal in some point. Could you tell me the number, how much % of such event can we have from tracking?

Thank you for your teaching.
Donghee



Re: MisID vs Impurity [message #14232 is a reply to message #14231] Fri, 23 November 2012 17:59 Go to previous messageGo to next message
donghee is currently offline  donghee
Messages: 385
Registered: January 2009
Location: Germnay
first-grade participant
From: *kph.uni-mainz.de
Dear all,

I have a idea about the definition of impurity and misID.
before getting the answer from Gosia

Gorsia showed us...
pi_imp = PID_{e} > X / pi_all
PID_{e} - it calls here probablity of being an electron for given particle. In this case for pion.

If we define PID_{e} > X as a
"selected true pion after doing MC PID match and requiring PID probability", then it refer to impurity.
or
"any kind of tracks passed required PID probability", then this quantity should be misID.

Above one is impurity as 1-purity and tell us how much % of wrong type particles are contributed in given PID selection.
Below quantity can have a meaning, how much % of a particle type can participate to other list of track candidate according given particle type.

Could we think two definition with this way?
Have a nice weekend,

Donghee





[Updated on: Fri, 23 November 2012 18:01]

Report message to a moderator

Re: MisID vs Impurity [message #14233 is a reply to message #14231] Sat, 24 November 2012 07:04 Go to previous messageGo to next message
Malgorzata Gumberidze is currently offline  Malgorzata Gumberidze
Messages: 98
Registered: June 2005
continuous participant
From: *pool.mediaWays.net
Hello


I'm wondering about the PID_{e} > X. 
Is it "selected true pion after doing MC PID match and requiring PID probability" 

Yes this is exactly what am I doing. I ask, that the particles is really pion by cross-checking MC id and than i apply to the primary PION cut on PID to be an electron.

You have selected only one reconstructed particle, which has a closest momentum value to MC one. 
We can have usually more than one track after reconstruction even from one event.
That means, all values of impurity shown in your categories are too ideal in some point. 
Could you tell me the number, how much % of such event can we have from tracking?


Of course we have more than one particle per event, but ideally you will try later in the experimental analysis also select one particle per event, assuming that you want to analyze one particle per event. If you do not do such selection, than in case of efficiency you will get more than 1.

For example in our case pbarp->e+e- in the old framework we were doing all combination of +- pair in the event, and than we were selecting to have only one
per event. The best one, looking to CHi2 from kinematical fit for example.

So I THINK that in this case it is fine what i do.

gosia
Re: MisID vs Impurity [message #14234 is a reply to message #14231] Sat, 24 November 2012 21:29 Go to previous messageGo to next message
Malgorzata Gumberidze is currently offline  Malgorzata Gumberidze
Messages: 98
Registered: June 2005
continuous participant
From: *fbx.proxad.net
hello

please find here plot which is showing multiplicity of reconstructed tracks for one true MC electron track :

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

and here are some numbers (number of primary electrons=972500):

mult ==0 : 83628
mult ==1 : 863344
mult ==2 : 24060
mult ==3 : 1116
mult ==4 : 239
mult ==5 : 64

and here is a plot for negative pions

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

and numbers (number of primary pions = 986500):

mult ==0 : 110975
mult ==1 : 854022
mult ==2 : 20824
mult ==3 : 550
mult ==4 : 90
mult ==5 : 29
mult ==6 : 5

3% of the events has more than 1 particles reconstructed per MC track in both cases for electrons and pions. So in this way i would say that selection of one good track per MC true one does not introduce big bias.

Re: MisID vs Impurity [message #14235 is a reply to message #14231] Sat, 24 November 2012 21:37 Go to previous messageGo to next message
Malgorzata Gumberidze is currently offline  Malgorzata Gumberidze
Messages: 98
Registered: June 2005
continuous participant
From: *fbx.proxad.net

few more explanations:

But as usual, I have still few question to better understand your definition of impurity.
pi_imp = PID_{e} > X / pi_all
I can understand about the nominator, that is reconstructed true pion, which is tested MC true PID matching. 

I'm wondering about the PID_{e} > X. 
Is it "selected true pion after doing MC PID match and requiring PID probability" 
or 
"any kind of tracks just passed given probability condition"? 
I am not clear for this which one have to be applied to see whatever impurity or misID. 


in both case i select MC true pions.
pi_all - are all TRUE MC pions without any condition on PID
PID_{e} - are TRUE MC pions with condition on PID of being electron

i hope that soon we will come to the common point with definitions Smile
Re: MisID vs Impurity [message #14236 is a reply to message #14221] Sun, 25 November 2012 10:07 Go to previous message
Ronald Kunne is currently offline  Ronald Kunne
Messages: 32
Registered: October 2009
continuous participant
From: *fbx.proxad.net
A physicist in Kuala Lumpur
said that her sample was rather impure:
"There's no denyin'
I mis-ID'd that pion."
Klaus Goetzen is right, for sure!
Laughing
Previous Topic: PID macro
Next Topic: SciTil implementation
Goto Forum:
  


Current Time: Tue Dec 03 04:16:01 CET 2024

Total time taken to generate the page: 0.00751 seconds