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 previous 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 833 times)
 
Read Message icon5.gif
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: doubt in tutorial
Next Topic: Faster PndPidCorrelator with smart MVD Info
Goto Forum:
  


Current Time: Sat Oct 12 05:08:44 CEST 2024

Total time taken to generate the page: 0.00839 seconds