GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Particle Identification » PID combiner with different detector
PID combiner with different detector [message #15658] Mon, 04 November 2013 22:02 Go to next message
donghee is currently offline  donghee
Messages: 385
Registered: January 2009
Location: Germnay
first-grade participant
From: *dip0.t-ipconnect.de
Hi all,

I have a fundamental question about PID combiner.
In our analysis we are using PID combiner, which combines the probability values from different detectors.

Let assume a simple situation.
I want to identify electron from EMC+MUO+STT+DRC combination.
In some cases, I assume that the probability from MUO should be zero due to absorbing the electron already in EMCalorimeter.
In practice, MUO doesn't contribute electron PID.
If I multiply P(EMC)XP(MUO), then total probability should be zero because of P(MUO)=0 and will be set as 0.2 which is an eqaul probability of 5 hypothesis.

So effective way to identify the electron should be EMC+STT+DRC combination without MUO.
This means that one need to define best combination for 5 different particles.

Is there some study on this issue? or can we recommend simply EMC+STT+DRC+MUO+DISC+MVD combination for each particle hypothesis in practice at PID analysis?

Best regards,
Donghee




[Updated on: Mon, 04 November 2013 22:02]

Report message to a moderator

Re: PID combiner with different detector [message #15659 is a reply to message #15658] Mon, 04 November 2013 22:32 Go to previous messageGo to next message
Ronald Kunne is currently offline  Ronald Kunne
Messages: 32
Registered: October 2009
continuous participant
From: *fbx.proxad.net
> In practice, MUO doesn't contribute electron PID.

Why not ?
The fact that MUO doesn't give a signal *is* information.

But in PandaRoot by convention "no signal" results in P=0.2 for all five particles. As a result Pe(EMC) * Pe(MUO) will be equal to Pe(EMC).

Greetings,
Ronald Kunne
Re: PID combiner with different detector [message #15660 is a reply to message #15659] Mon, 04 November 2013 23:16 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*
In particular, if there is no MUO signal, the probability should be 0 to be muon and 0.25 to be pion kaon electron proton.
Re: PID combiner with different detector [message #15662 is a reply to message #15659] Tue, 05 November 2013 11:05 Go to previous messageGo to next message
donghee is currently offline  donghee
Messages: 385
Registered: January 2009
Location: Germnay
first-grade participant
From: *scc.kit.edu
Dear Ronald,

Now I am clear for the eqaul probability for absent PID info in certain detector.

If I use EMC and MUO, and an electron will identify with a single piece of detector as like
Pe(EMC) = 0.9 (90% probability at EMC)
Pe(MUO) = 0 ( no infomation at MUO)

For other particles with EMC,
Pmu(EMC) = 0.1
Ppi,k,p(EMC) = 0.5

And for MUO detector, the probabilities will be reset as 0.2 even for all other particles.
Pe(MUO) = 0.2
Pmu,pi,k,p(MUO) = 0.2


Then will calculate a global probability as like
Pe(EMC,MUO) = Pe(EMC) *Pe(MUO) = 0.9*0.2 = 0.18
Pmu(EMC,MUO) = Pmu(EMC)*Pmu(MUO) = 0.1*0.2 = 0.02
Ppi(EMC,MUO) = Ppi(EMC)*Ppi(MUO) = 0.5*0.2 = 0.1
Pka(EMC,MUO) = Pka(EMC)*Pka(MUO) = 0.5*0.2 = 0.1
Ppr(EMC,MUO) = Ppr(EMC)*Ppr(MUO) = 0.5*0.2 = 0.1
and so on.
After that will be normalized with
Pe(EMC,MUO) + Pmu(EMC,MUO) + Ppi(EMC,MUO) + Pka(EMC,MUO) + Ppr(EMC,MUO) = 0.5

So finally I can have normalized global PID probabilities
Pe(EMC,MUO) = 0.18/0.5 = 0.36
Pmu(EMC,MUO) = 0.02/0.5 = 0.04
Ppi(EMC,MUO) = 0.1/0.5 = 0.20
Pka(EMC,MUO) = 0.1/0.5 = 0.20
Ppr(EMC,MUO) = 0.1/0.5 = 0.20

This is a story of PID!
If I am wrong, correct me again.

Thanks,
Donghee


[Updated on: Tue, 05 November 2013 11:12]

Report message to a moderator

Re: PID combiner with different detector [message #15663 is a reply to message #15662] Tue, 05 November 2013 11:26 Go to previous messageGo to next message
donghee is currently offline  donghee
Messages: 385
Registered: January 2009
Location: Germnay
first-grade participant
From: *scc.kit.edu
Dear Ronald and stefano,

If I see a band at 0.2 in PID with usage of global probability and many detector types, that means
there are very poor information from all detector or
are most likely ghost tracks and low energetic electrons.

(see attached plot from Klaus Gotzen)

Best wishes,
Donghee
Re: PID combiner with different detector [message #15664 is a reply to message #15663] Tue, 05 November 2013 11:29 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*
This means that most probably some information was missing in the candidate, i.e. the track was not correlated to EMC DRC etc. Most probably they will be forward tracks, where no PID info is present.
Re: PID combiner with different detector [message #15665 is a reply to message #15664] Tue, 05 November 2013 11:31 Go to previous messageGo to next message
donghee is currently offline  donghee
Messages: 385
Registered: January 2009
Location: Germnay
first-grade participant
From: *scc.kit.edu
Thanks to all teachers.

close session!
Re: PID combiner with different detector [message #15666 is a reply to message #15662] Tue, 05 November 2013 11:39 Go to previous message
Ronald Kunne is currently offline  Ronald Kunne
Messages: 32
Registered: October 2009
continuous participant
From: *in2p3.fr
Your example is a bit off, as all the probas should add up to 1.

p(EMC) = 0.9 for electron, 0.025 for each of the others
p(MUO) = 0.2 for all particle.

Then we have:
p(EMC)*p(MUO) = 0.9 * 0.2 =0.18 for the electron
p(EMC)*p(MUO) = 0.025 * 0.2 =0.005 for the others
This adds up to 0.2, so the final result is

p(EMC)*p(MUO) = 0.18/0.2 = 0.9 for the electron,
0.005/0.2 = 0.025 for each of the others, as expected.

Quote:

If I see a band at 0.2 in PID with usage of global probability and many detector types, that means
there are very poor information from all detector or
are most likely ghost tracks and low energetic electrons.


Or particles falling outside the acceptance of the detector, or outside the momentum range 0.2 < p < 5 GeV/c for which the calculation was made.





[Updated on: Tue, 05 November 2013 11:41]

Report message to a moderator

Previous Topic: An option in the PID selector
Next Topic: doubt in tutorial
Goto Forum:
  


Current Time: Thu Mar 28 13:00:52 CET 2024

Total time taken to generate the page: 0.00911 seconds