@@ -876,28 +876,6 @@ struct NchCumulantsId {
876876 fillIdentificationQA<mode, kPr , NoId>(hist, track); // Look at Proton
877877 }
878878
879- // template <int Mode, int pidMode, typename H, typename T>
880- // void fillGenTrackQA(H& histReg, const T& mcTrack)
881- // {
882- // histReg.fill(HIST(HistRegDire2[Mode]) + HIST(PidDire[pidMode]) + HIST("h12_p"), mcTrack.p());
883- // histReg.fill(HIST(HistRegDire2[Mode]) + HIST(PidDire[pidMode]) + HIST("h13_pt"), mcTrack.pt());
884- // histReg.fill(HIST(HistRegDire2[Mode]) + HIST(PidDire[pidMode]) + HIST("h14_eta"), mcTrack.eta());
885- // histReg.fill(HIST(HistRegDire2[Mode]) + HIST(PidDire[pidMode]) + HIST("h15_phi"), mcTrack.phi());
886- // histReg.fill(HIST(HistRegDire2[Mode]) + HIST(PidDire[pidMode]) + HIST("h16_rapidity"), mcTrack.y());
887- // histReg.fill(HIST(HistRegDire2[Mode]) + HIST(PidDire[pidMode]) + HIST("h20_pt_eta"), mcTrack.pt(), mcTrack.eta());
888- // }
889- //
890- // template <int Mode, int pidMode, typename H, typename T>
891- // void fillRecoTrackQA(H& histReg, const T& track)
892- // {
893- // histReg.fill(HIST(HistRegDire2[Mode]) + HIST(PidDire[pidMode]) + HIST("h12_p"), track.p());
894- // histReg.fill(HIST(HistRegDire2[Mode]) + HIST(PidDire[pidMode]) + HIST("h13_pt"), track.pt());
895- // histReg.fill(HIST(HistRegDire2[Mode]) + HIST(PidDire[pidMode]) + HIST("h14_eta"), track.eta());
896- // histReg.fill(HIST(HistRegDire2[Mode]) + HIST(PidDire[pidMode]) + HIST("h15_phi"), track.phi());
897- // histReg.fill(HIST(HistRegDire2[Mode]) + HIST(PidDire[pidMode]) + HIST("h16_rapidity"), track.y());
898- // histReg.fill(HIST(HistRegDire2[Mode]) + HIST(PidDire[pidMode]) + HIST("h20_pt_eta"), track.pt(), track.eta());
899- // }
900-
901879 template <HistRegEnum2 DIR, PidEnum P, ChargeEnum C, typename TrackType>
902880 void fillGenTrackQA (HistogramRegistry& r, const TrackType& t)
903881 {
@@ -1230,14 +1208,6 @@ struct NchCumulantsId {
12301208 trackIsProton = true ;
12311209 BITSET (trackIdTag, ID_BIT_PR);
12321210 }
1233- // if (selElectron(track, idMethodEl)) {
1234- // trackIsElectron = true;
1235- // BITSET(trackIdTag, ID_BIT_EL);
1236- // }
1237- // if (selDeuteron(track, idMethodDe)) {
1238- // trackIsDeuteron = true;
1239- // BITSET(trackIdTag, ID_BIT_DE);
1240- // }
12411211
12421212 if constexpr (analysisType == doPurityProcessing) {
12431213 if (trackIsPion) {
@@ -1350,7 +1320,6 @@ struct NchCumulantsId {
13501320 int pdg = mcTrack.pdgCode ();
13511321
13521322 if (pdg == kPiPlus || pdg == kKPlus || pdg == kProton || pdg == kPositron || pdg == kMuonPlus || pdg == kDeuteron ) {
1353- // fillGenTrackQA<genAnalysisDir, kPos>(genAnalysis, mcTrack);
13541323 nP++;
13551324
13561325 genAnalysis.fill (HIST (" genAnalysis/Charge/Pos/h12_p" ), mcTrack.p ());
@@ -1360,7 +1329,6 @@ struct NchCumulantsId {
13601329 genAnalysis.fill (HIST (" genAnalysis/Charge/Pos/h16_rapidity" ), mcTrack.y ());
13611330 genAnalysis.fill (HIST (" genAnalysis/Charge/Pos/h20_pt_eta" ), mcTrack.pt (), mcTrack.eta ());
13621331 } else if (pdg == kPiMinus || pdg == kKMinus || pdg == kProtonBar || pdg == kElectron || pdg == kMuonMinus || pdg == -kDeuteron ) {
1363- // fillGenTrackQA<genAnalysisDir, kNeg>(genAnalysis, mcTrack);
13641332 nM++;
13651333
13661334 genAnalysis.fill (HIST (" genAnalysis/Charge/Neg/h12_p" ), mcTrack.p ());
@@ -1372,33 +1340,21 @@ struct NchCumulantsId {
13721340 }
13731341 // ----- Pions -----
13741342 if (pdg == kPiPlus ) {
1375- // fillGenTrackQA<genAnalysisDir, kPi>(genAnalysis, mcTrack);
13761343 nPi++;
13771344 fillGenTrackQA<genAnalysisDir, kPi , kPos >(genAnalysis, mcTrack);
13781345 } else if (pdg == kPiMinus ) {
1379- // fillGenTrackQA<genAnalysisDir, kPi>(genAnalysis, mcTrack);
13801346 nAPi++;
13811347 fillGenTrackQA<genAnalysisDir, kPi , kNeg >(genAnalysis, mcTrack);
1382- }
1383-
1384- // ----- Kaons -----
1385- else if (pdg == kKPlus ) {
1386- // fillGenTrackQA<genAnalysisDir, kKa>(genAnalysis, mcTrack);
1348+ } else if (pdg == kKPlus ) {
13871349 nKa++;
13881350 fillGenTrackQA<genAnalysisDir, kKa , kPos >(genAnalysis, mcTrack);
13891351 } else if (pdg == kKMinus ) {
1390- // fillGenTrackQA<genAnalysisDir, kKa>(genAnalysis, mcTrack);
13911352 nAKa++;
13921353 fillGenTrackQA<genAnalysisDir, kKa , kNeg >(genAnalysis, mcTrack);
1393- }
1394-
1395- // ----- Protons -----
1396- else if (pdg == kProton ) {
1397- // fillGenTrackQA<genAnalysisDir, kPr>(genAnalysis, mcTrack);
1354+ } else if (pdg == kProton ) {
13981355 nPr++;
13991356 fillGenTrackQA<genAnalysisDir, kPr , kPos >(genAnalysis, mcTrack);
14001357 } else if (pdg == kProtonBar ) {
1401- // fillGenTrackQA<genAnalysisDir, kPr>(genAnalysis, mcTrack);
14021358 nAPr++;
14031359 fillGenTrackQA<genAnalysisDir, kPr , kNeg >(genAnalysis, mcTrack);
14041360 }
@@ -1463,7 +1419,6 @@ struct NchCumulantsId {
14631419 int pdg = mcTrack.pdgCode ();
14641420
14651421 if (pdg == kPiPlus || pdg == kKPlus || pdg == kProton || pdg == kPositron || pdg == kMuonPlus || pdg == kDeuteron ) {
1466- // fillGenTrackQA<genAnalysisDir, kPos>(genAnalysis, mcTrack);
14671422 nPGen++;
14681423
14691424 genAnalysis.fill (HIST (" genAnalysis/Charge/Pos/h12_p" ), mcTrack.p ());
@@ -1473,7 +1428,6 @@ struct NchCumulantsId {
14731428 genAnalysis.fill (HIST (" genAnalysis/Charge/Pos/h16_rapidity" ), mcTrack.y ());
14741429 genAnalysis.fill (HIST (" genAnalysis/Charge/Pos/h20_pt_eta" ), mcTrack.pt (), mcTrack.eta ());
14751430 } else if (pdg == kPiMinus || pdg == kKMinus || pdg == kProtonBar || pdg == kElectron || pdg == kMuonMinus || pdg == -kDeuteron ) {
1476- // fillGenTrackQA<genAnalysisDir, kNeg>(genAnalysis, mcTrack);
14771431 nMGen++;
14781432
14791433 genAnalysis.fill (HIST (" genAnalysis/Charge/Neg/h12_p" ), mcTrack.p ());
@@ -1483,35 +1437,23 @@ struct NchCumulantsId {
14831437 genAnalysis.fill (HIST (" genAnalysis/Charge/Neg/h16_rapidity" ), mcTrack.y ());
14841438 genAnalysis.fill (HIST (" genAnalysis/Charge/Neg/h20_pt_eta" ), mcTrack.pt (), mcTrack.eta ());
14851439 }
1486- // ----- Pions -----
1440+
14871441 if (pdg == kPiPlus ) {
1488- // fillGenTrackQA<genAnalysisDir, kPi>(genAnalysis, mcTrack);
14891442 nPiGen++;
14901443 fillGenTrackQA<genAnalysisDir, kPi , kPos >(genAnalysis, mcTrack);
14911444 } else if (pdg == kPiMinus ) {
1492- // fillGenTrackQA<genAnalysisDir, kPi>(genAnalysis, mcTrack);
14931445 nAPiGen++;
14941446 fillGenTrackQA<genAnalysisDir, kPi , kNeg >(genAnalysis, mcTrack);
1495- }
1496-
1497- // ----- Kaons -----
1498- else if (pdg == kKPlus ) {
1499- // fillGenTrackQA<genAnalysisDir, kKa>(genAnalysis, mcTrack);
1447+ } else if (pdg == kKPlus ) {
15001448 nKaGen++;
15011449 fillGenTrackQA<genAnalysisDir, kKa , kPos >(genAnalysis, mcTrack);
15021450 } else if (pdg == kKMinus ) {
1503- // fillGenTrackQA<genAnalysisDir, kKa>(genAnalysis, mcTrack);
15041451 nAKaGen++;
15051452 fillGenTrackQA<genAnalysisDir, kKa , kNeg >(genAnalysis, mcTrack);
1506- }
1507-
1508- // ----- Protons -----
1509- else if (pdg == kProton ) {
1510- // fillGenTrackQA<genAnalysisDir, kPr>(genAnalysis, mcTrack);
1453+ } else if (pdg == kProton ) {
15111454 nPrGen++;
15121455 fillGenTrackQA<genAnalysisDir, kPr , kPos >(genAnalysis, mcTrack);
15131456 } else if (pdg == kProtonBar ) {
1514- // fillGenTrackQA<genAnalysisDir, kPr>(genAnalysis, mcTrack);
15151457 nAPrGen++;
15161458 fillGenTrackQA<genAnalysisDir, kPr , kNeg >(genAnalysis, mcTrack);
15171459 }
0 commit comments