Skip to content
Snippets Groups Projects
Commit a837de70 authored by Sergey Gorbunov's avatar Sergey Gorbunov Committed by Volker Friese
Browse files

KF: cosmetic changes in CbmKFParticleFinderPID

parent d404c622
No related branches found
No related tags found
1 merge request!1066KF: fixes in KFParticleFinder according to the new TOF geometry
Pipeline #20883 passed
...@@ -79,7 +79,7 @@ CbmKFParticleFinderPID::CbmKFParticleFinderPID(const char* name, Int_t iVerbose) ...@@ -79,7 +79,7 @@ CbmKFParticleFinderPID::CbmKFParticleFinderPID(const char* name, Int_t iVerbose)
, fMuchMode(0) , fMuchMode(0)
, fUseSTSdEdX(kFALSE) , fUseSTSdEdX(kFALSE)
, fUseTRDdEdX(kFALSE) , fUseTRDdEdX(kFALSE)
, fTimeSliceMode(0) , fLegacyEventMode(0)
, fPID(0) , fPID(0)
{ {
//MuCh cuts //MuCh cuts
...@@ -104,9 +104,9 @@ InitStatus CbmKFParticleFinderPID::Init() ...@@ -104,9 +104,9 @@ InitStatus CbmKFParticleFinderPID::Init()
} }
//check the mode //check the mode
fTimeSliceMode = 0; fLegacyEventMode = 1;
if (ioman->CheckBranch("CbmEvent")) { if (ioman->CheckBranch("CbmEvent")) {
fTimeSliceMode = 1; fLegacyEventMode = 0;
LOG(info) << GetName() << ": Running in the timeslice mode."; LOG(info) << GetName() << ": Running in the timeslice mode.";
} }
else else
...@@ -122,16 +122,12 @@ InitStatus CbmKFParticleFinderPID::Init() ...@@ -122,16 +122,12 @@ InitStatus CbmKFParticleFinderPID::Init()
CbmMCDataManager* mcManager = 0; CbmMCDataManager* mcManager = 0;
if (fTimeSliceMode) mcManager = (CbmMCDataManager*) fManger->GetObject("MCDataManager"); if (!fLegacyEventMode) {
mcManager = (CbmMCDataManager*) fManger->GetObject("MCDataManager");
if (fTimeSliceMode) {
if (mcManager == 0) { if (mcManager == 0) {
Fatal("CbmKFParticleFinderPID::Init", "MC Data Manager is not found!"); Fatal("CbmKFParticleFinderPID::Init", "MC Data Manager is not found!");
return kERROR; return kERROR;
} }
}
if (fTimeSliceMode) {
fMCTracks = mcManager->InitBranch("MCTrack"); fMCTracks = mcManager->InitBranch("MCTrack");
if (fMCTracks == 0) { if (fMCTracks == 0) {
Fatal("CbmKFParticleFinderPID::Init", "MC track array not found!"); Fatal("CbmKFParticleFinderPID::Init", "MC track array not found!");
...@@ -255,7 +251,7 @@ void CbmKFParticleFinderPID::SetMCPID() ...@@ -255,7 +251,7 @@ void CbmKFParticleFinderPID::SetMCPID()
{ {
Int_t nTracks = fTrackArray->GetEntriesFast(); Int_t nTracks = fTrackArray->GetEntriesFast();
Int_t nMCTracks = 0; Int_t nMCTracks = 0;
if (!fTimeSliceMode) nMCTracks = fMCTrackArray->GetEntriesFast(); if (fLegacyEventMode) nMCTracks = fMCTrackArray->GetEntriesFast();
for (int iTr = 0; iTr < nTracks; iTr++) { for (int iTr = 0; iTr < nTracks; iTr++) {
fPID[iTr] = -2; fPID[iTr] = -2;
...@@ -273,7 +269,7 @@ void CbmKFParticleFinderPID::SetMCPID() ...@@ -273,7 +269,7 @@ void CbmKFParticleFinderPID::SetMCPID()
if (stsTrackMatch->GetLink(iLink).GetWeight() > bestWeight) { if (stsTrackMatch->GetLink(iLink).GetWeight() > bestWeight) {
bestWeight = stsTrackMatch->GetLink(iLink).GetWeight(); bestWeight = stsTrackMatch->GetLink(iLink).GetWeight();
mcTrackId = stsTrackMatch->GetLink(iLink).GetIndex(); mcTrackId = stsTrackMatch->GetLink(iLink).GetIndex();
if (fTimeSliceMode) { if (!fLegacyEventMode) {
iFile = stsTrackMatch->GetLink(iLink).GetFile(); iFile = stsTrackMatch->GetLink(iLink).GetFile();
iEvent = stsTrackMatch->GetLink(iLink).GetEntry(); iEvent = stsTrackMatch->GetLink(iLink).GetEntry();
} }
...@@ -281,7 +277,7 @@ void CbmKFParticleFinderPID::SetMCPID() ...@@ -281,7 +277,7 @@ void CbmKFParticleFinderPID::SetMCPID()
} }
if (bestWeight / totalWeight < 0.7) continue; if (bestWeight / totalWeight < 0.7) continue;
if ((!fTimeSliceMode) && (mcTrackId >= nMCTracks || mcTrackId < 0)) continue; if ((fLegacyEventMode) && (mcTrackId >= nMCTracks || mcTrackId < 0)) continue;
// if(mcTrackId >= nMCTracks || mcTrackId < 0) // if(mcTrackId >= nMCTracks || mcTrackId < 0)
// { // {
// LOG(info) << "Sts Matching is wrong! StsTrackId = " << mcTrackId << " N mc tracks = " << nMCTracks; // LOG(info) << "Sts Matching is wrong! StsTrackId = " << mcTrackId << " N mc tracks = " << nMCTracks;
...@@ -293,10 +289,10 @@ void CbmKFParticleFinderPID::SetMCPID() ...@@ -293,10 +289,10 @@ void CbmKFParticleFinderPID::SetMCPID()
CbmMCTrack* cbmMCTrack = 0; CbmMCTrack* cbmMCTrack = 0;
if (fTimeSliceMode) cbmMCTrack = dynamic_cast<CbmMCTrack*>(fMCTracks->Get(iFile, iEvent, mcTrackId)); if (!fLegacyEventMode) { cbmMCTrack = dynamic_cast<CbmMCTrack*>(fMCTracks->Get(iFile, iEvent, mcTrackId)); }
else else {
cbmMCTrack = (CbmMCTrack*) fMCTrackArray->At(mcTrackId); cbmMCTrack = (CbmMCTrack*) fMCTrackArray->At(mcTrackId);
}
if (!(TMath::Abs(cbmMCTrack->GetPdgCode()) == 11 || TMath::Abs(cbmMCTrack->GetPdgCode()) == 13 if (!(TMath::Abs(cbmMCTrack->GetPdgCode()) == 11 || TMath::Abs(cbmMCTrack->GetPdgCode()) == 13
|| TMath::Abs(cbmMCTrack->GetPdgCode()) == 211 || TMath::Abs(cbmMCTrack->GetPdgCode()) == 321 || TMath::Abs(cbmMCTrack->GetPdgCode()) == 211 || TMath::Abs(cbmMCTrack->GetPdgCode()) == 321
......
...@@ -125,7 +125,7 @@ private: ...@@ -125,7 +125,7 @@ private:
Bool_t fUseSTSdEdX; Bool_t fUseSTSdEdX;
Bool_t fUseTRDdEdX; Bool_t fUseTRDdEdX;
bool fTimeSliceMode; bool fLegacyEventMode;
//MuCh cuts //MuCh cuts
float fMuchCutsFloat[2]; float fMuchCutsFloat[2];
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment