diff --git a/reco/detectors/trd/pid/CbmTrdSetTracksPidLike.cxx b/reco/detectors/trd/pid/CbmTrdSetTracksPidLike.cxx index 12103694c3e369f8b84d15472764219dc32b619c..6100edd615b6caa343dc1c8950b80a99db842689 100644 --- a/reco/detectors/trd/pid/CbmTrdSetTracksPidLike.cxx +++ b/reco/detectors/trd/pid/CbmTrdSetTracksPidLike.cxx @@ -98,8 +98,8 @@ Bool_t CbmTrdSetTracksPidLike::ReadData() if (fMCinput) { /// mc pid method if (fMomDep) { - std::vector<TString> histnames {"MC_electron_p_eloss", "MC_pion_p_eloss", "MC_kaon_p_eloss", "MC_proton_p_eloss", - "MC_muon_p_eloss"}; + std::vector<TString> histnames{"MC_electron_p_eloss", "MC_pion_p_eloss", "MC_kaon_p_eloss", "MC_proton_p_eloss", + "MC_muon_p_eloss"}; inArr = new TObjArray(histnames.size()); for (size_t i = 0; i < histnames.size(); i++) { h[i] = histFile->Get<TH2F>(histnames[i]); @@ -124,8 +124,8 @@ Bool_t CbmTrdSetTracksPidLike::ReadData() } } if (!fMomDep) { - std::vector<TString> histnames {"MC_electron_eloss", "MC_pion_eloss", "MC_kaon_eloss", "MC_proton_eloss", - "MC_muon_eloss"}; + std::vector<TString> histnames{"MC_electron_eloss", "MC_pion_eloss", "MC_kaon_eloss", "MC_proton_eloss", + "MC_muon_eloss"}; inArr = new TObjArray(histnames.size()); for (size_t i = 0; i < histnames.size(); i++) { h[i] = histFile->Get<TH1F>(histnames[i]); @@ -146,8 +146,8 @@ Bool_t CbmTrdSetTracksPidLike::ReadData() } else { /// data driven method if (fMomDep) { - std::vector<TString> histnames {"ELE_electron_p_eloss", "PIO_pion_p_eloss", "ELE_kaon_p_eloss", - "ELE_proton_p_eloss", "ELE_muon_p_eloss"}; + std::vector<TString> histnames{"ELE_electron_p_eloss", "PIO_pion_p_eloss", "ELE_kaon_p_eloss", + "ELE_proton_p_eloss", "ELE_muon_p_eloss"}; inArr = new TObjArray(histnames.size()); for (size_t i = 0; i < histnames.size(); i++) { h[i] = histFile->Get<TH2F>(histnames[i]); @@ -173,8 +173,8 @@ Bool_t CbmTrdSetTracksPidLike::ReadData() } } if (!fMomDep) { - std::vector<TString> histnames {"ELE_electron_eloss", "PIO_pion_eloss", "ELE_kaon_eloss", "ELE_proton_eloss", - "ELE_muon_eloss"}; + std::vector<TString> histnames{"ELE_electron_eloss", "PIO_pion_eloss", "ELE_kaon_eloss", "ELE_proton_eloss", + "ELE_muon_eloss"}; inArr = new TObjArray(histnames.size()); for (size_t i = 0; i < histnames.size(); i++) { h[i] = histFile->Get<TH1F>(histnames[i]); @@ -203,7 +203,8 @@ Bool_t CbmTrdSetTracksPidLike::ReadData() LOG_IF(info, hist->GetEntries() < 1000) << "Input histogram is almost empty for" << name.Data(); // check particles - if (name.Contains("electron")) particle = CbmTrdSetTracksPidLike::kElectron; + if (name.Contains("electron")) + particle = CbmTrdSetTracksPidLike::kElectron; else if (name.Contains("pion")) particle = CbmTrdSetTracksPidLike::kPion; else if (name.Contains("kaon")) @@ -312,7 +313,8 @@ void CbmTrdSetTracksPidLike::Exec(Option_t*) } /// only trd tracks with mimimum 1 reconstructed point - if (pTrack->GetNofHits() < 1) continue; + if (pTrack->GetNofHits() < 1) + continue; else fNofTracks++; @@ -390,17 +392,25 @@ Double_t CbmTrdSetTracksPidLike::GetProbability(Int_t k, Double_t mom, Double_t // /// check for undefined dedx (dedx = -1. is sometimes passed in) - if (dedx < 0.) { return -999.; } + if (dedx < 0.) { + return -999.; + } /// useless protection - if (k < 0 || k > fgkNParts) { return -999.; } + if (k < 0 || k > fgkNParts) { + return -999.; + } /// histogram has TRD momentum at inner point vs. dedx sinal [keV] TH1* hist = (TH1*) fHistdEdx->At(k); - if (!hist) { return -999.; } + if (!hist) { + return -999.; + } /// check for entries/ non-empty histograms - if (hist->GetEntries() < 1000.) { return -999.; } + if (hist->GetEntries() < 1000.) { + return -999.; + } Int_t ndim = hist->GetDimension();