GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Particle Identification » Charged Hits in Neutral List
icon5.gif  Charged Hits in Neutral List [message #17139] Fri, 08 August 2014 13:54 Go to next message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *ikp.kfa-juelich.de
Dear PID experts,
I simulated 10 charged kaon in box generator with the momentum (1,4) GeV/c and checked the photon multiplicity in this K+ box like this:
---event loop---
theAnalysis->FillList(neu,"Neutral");
hgeve->Fill(neu.GetLength()); 

index.php?t=getfile&id=8011&private=0
As shown in the attached .png, in the 10 evt K+, we got avg. 5.7 neutrals and 2 of those events got surprising 16 neutrals.

To understand this in a deeper level, Tobias helped to check the EMC cluster in MCTrack:
fMCMatch->CreateArtificialStage("MCTrack");
Int_t nCluster = fClusterArray->GetEntriesFast();
Int_t nTracks = fMCTrackArray->GetEntriesFast();
//loop to test Cluster
	for (Int_t iCluster=0; iCluster<nCluster; iCluster++)
	{
		PndEmcCluster* theCluster = (PndEmcCluster*) fClusterArray->At(iCluster);
		std::vector<Int_t> mcIndexes = theCluster->GetMcList();
		for (Int_t iMcIndex = 0; iMcIndex < mcIndexes.size(); iMcIndex++){
			std::cout << "Tracks in Cluster " << iCluster << " TrackId: " << mcIndexes.at(iMcIndex); // << std::endl;
			PndMCTrack* myTrack = (PndMCTrack*)fMCTrackArray->At(mcIndexes.at(iMcIndex));
			std::cout << " PdgCode: " << myTrack->GetPdgCode() << std::endl;
		}
	}


Then, we got this info:
Tracks in Cluster 1 TrackId: 0 PdgCode: 321
Tracks in Cluster 2 TrackId: 0 PdgCode: 321
Tracks in Cluster 3 TrackId: 0 PdgCode: 321
Tracks in Cluster 4 TrackId: 0 PdgCode: 321
Tracks in Cluster 5 TrackId: 0 PdgCode: 321
...


It looks some charged particles fired the EMC crystals but they are treated as "neutral candidates".

With best regards,
Lu
  • Attachment: Canvas_1.png
    (Size: 18.24KB, Downloaded 746 times)
Re: Charged Hits in Neutral List [message #17140 is a reply to message #17139] Fri, 08 August 2014 15:25 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *20-87-r.retail.telecomitalia.it
This means they are not correlated to the charged track. Try to check the GetEmcQuality, which is the distance?2 to the closest track. You should put a cut on minimum distance there.
Re: Charged Hits in Neutral List [message #17143 is a reply to message #17140] Fri, 08 August 2014 16:22 Go to previous messageGo to next message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *ikp.kfa-juelich.de
Dear Sefano,

I checked the EmcQuality in the output of pid, but I don't understand why they're all -1.
Do I need to enable any functions at somewhere? The code I'm using is a copy of macro/run/*_complete.C.
index.php?t=getfile&id=8012&private=0

And one more interesting thing we found within the eventdisplay, is that several photons (light pink) from the charged kaon (light yellow) fly backward to hit the emc:
index.php?t=getfile&id=8013&private=0
Best,
Lu
  • Attachment: Canvas_1.png
    (Size: 7.43KB, Downloaded 498 times)
  • Attachment: photons.png
    (Size: 72.15KB, Downloaded 617 times)
Re: Charged Hits in Neutral List [message #17146 is a reply to message #17143] Fri, 08 August 2014 18:56 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *250-95-r.retail.telecomitalia.it
Those numbers should be different from -1. Could you please check with higher statistics?
Re: Charged Hits in Neutral List [message #17147 is a reply to message #17146] Mon, 11 August 2014 09:14 Go to previous messageGo to next message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *ikp.kfa-juelich.de
Dear Stefano,

I tried 10k evt and got this following output:
index.php?t=getfile&id=8014&private=0
Should I cut at some value?

Best,
Lu
  • Attachment: Canvas_1.png
    (Size: 6.29KB, Downloaded 525 times)

[Updated on: Mon, 11 August 2014 09:14]

Report message to a moderator

Re: Charged Hits in Neutral List [message #17163 is a reply to message #17147] Tue, 12 August 2014 12:59 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *11-87-r.retail.telecomitalia.it
This should be the distance squared of the cluster position (in cm^2) from the projected track. This means that EmcQuality>2500 means that you consider neutral candidates only when the closest track is at at least 50 cm.
You could try with several selections, to see that happens.
Re: Charged Hits in Neutral List [message #17180 is a reply to message #17163] Thu, 14 August 2014 08:33 Go to previous messageGo to next message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *ikp.kfa-juelich.de
Many thanks. I'll have a try.
One more question...What's the meaning when it equals to -1 ?
Re: Charged Hits in Neutral List [message #17181 is a reply to message #17180] Thu, 14 August 2014 08:48 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *180-80-r.retail.telecomitalia.it
-1 should mean that there is no charged track nearby.
Re: Charged Hits in Neutral List [message #17204 is a reply to message #17181] Thu, 21 August 2014 16:21 Go to previous messageGo to next message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *ikp.kfa-juelich.de
Hi Stefano,
I just looked into the PidCorrelator code and found some reasons for -1:

* no MC or no MCindex found when ideal pid is used
* cluster energy below emc barrel threshold 20 MeV
* hit not come from emc or fsc
* no propagation at module when geane is used

** too low momenta (< 0.1) or too high momenta (> 15) // ?? why 0.1 GeV/c
** track flag <= 0
** no track info

And the "nearby" criteria of emcquality < 2500 is distinguishing the neutral candidate and the track-correlated hit. Then, the track-correlated hit will be constructed as a charged candidate.
Do I understand it right?

Best regards,
Lu
Re: Charged Hits in Neutral List [message #17208 is a reply to message #17204] Thu, 21 August 2014 22:06 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *211-62-r.retail.telecomitalia.it
These are criteria for the charged tracks, but you are interest in the neutral candidates.
Re: Charged Hits in Neutral List [message #17211 is a reply to message #17208] Thu, 21 August 2014 23:51 Go to previous messageGo to next message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *adsl.alicedsl.de
What I've checked is these criteria also affect on the neutral candidates. In the function of ConstructNeutralCandi, the emc quality is obtained from the map fClusterQ, which is stored and updated in the function of GetEmcInfo. GetEmcInfo is only called by ConstructChargedCandi.

This means that once the fClusterQ[i] is assigned as -1 in the step of constructing charged candidate, the emc quality of cluster/bump i for neutral candidate will be also -1.

Could you please have a look at the code, especially the operations regarding to fClusterQ? If it doesn't perform as expectation, there would be something "strange".

Thanks in advance.
Lu
Re: Charged Hits in Neutral List [message #17230 is a reply to message #17211] Mon, 01 September 2014 22:32 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *ip71.fastwebnet.it
This code was written as placeholder but never checked carefully. We need somebody taking care of neutral reconstruction.
Re: Charged Hits in Neutral List [message #17232 is a reply to message #17230] Tue, 02 September 2014 11:34 Go to previous messageGo to next message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *tagung.uni-giessen.de
I'd like to be a volunteer and have a discussion for more details in the coming CM.
Best regards,
Lu
Re: Charged Hits in Neutral List [message #17276 is a reply to message #17232] Wed, 10 September 2014 17:58 Go to previous messageGo to next message
Shyam Kumar is currently offline  Shyam Kumar
Messages: 78
Registered: September 2012
Location: Mumbai, Maharashtra India
continuous participant
From: 103.21.126*
Hi Lu,
I am also doing the same thing for the cut parameter optimization for Ftof, It will be useful in every case whichever algorithm I use. I have read in detail about it. I can give some idea what i understand stefano can tell whether it is correct.

Actually when a charged particle passes through medium and detector it deviates its actual trajectory (multiple scattering) so only from the actual hit it is very difficult to get the actual path. so kalman filter is used for this. It has three steps extrapolation, filtering, and smoothing. In the extrapolation we predict the theortical value and one is the measured hit on detector then in filtering we do the weighted mean of this to get kalman fitted value and then smoothing is used to get the true point.

So In my case Tof quality (as TofCut in all.par file) is the distance square of the measured position (actual hit) and extrapolated position(theortical predicted value) if it will be small then we will get the better result.

The all parameter is in all.par file https://subversion.gsi.de/trac/fairroot/browser/pandaroot/release/dec13/ macro/params/all.par (412-428). EmC quality should be Emc12Cut in all.par in line 418.

Thank you
Shyam

[Updated on: Wed, 10 September 2014 18:05]

Report message to a moderator

Re: Charged Hits in Neutral List [message #17277 is a reply to message #17276] Thu, 11 September 2014 09:23 Go to previous message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *lnf.infn.it
Hi Shyam,
good to know the correlator with tof, but the situation is quite different with the case of emc. For a more detailed info on the emc issue, you can have a look at my slides which has been presented in the collaboration meeting: https://indico.gsi.de/conferenceDisplay.py?confId=2943
The problem we have now is not only parameter optimization, but deeper in the level of algorithm itself.

With best regards,
Lu
Previous Topic: doubt in tutorial
Next Topic: Faster PndPidCorrelator with smart MVD Info
Goto Forum:
  


Current Time: Thu Mar 28 19:35:59 CET 2024

Total time taken to generate the page: 0.01047 seconds