@@ -281,14 +281,14 @@ struct Kstarqa {
281281 hPID.add (" Before/hNsigma_TPC_TOF_Ka_before" , " N #sigma Kaon TOF before" , kTH2F , {{50 , -5 .0f , 5 .0f }, {50 , -5 .0f , 5 .0f }});
282282 hPID.add (" Before/hNsigma_TPC_TOF_Pi_before" , " N #sigma Pion TOF before" , kTH2F , {{50 , -5 .0f , 5 .0f }, {50 , -5 .0f , 5 .0f }});
283283
284- hPID.add (" Before/h1PID_TPC_pos_kaon" , " Kaon PID distribution in data" , kTH1F , {{100 , -10 .0f , 10 .0f }});
285- hPID.add (" Before/h1PID_TPC_pos_pion" , " Pion PID distribution in data" , kTH1F , {{100 , -10 .0f , 10 .0f }});
286- hPID.add (" Before/h1PID_TOF_pos_kaon" , " Kaon PID distribution in data" , kTH1F , {{100 , -10 .0f , 10 .0f }});
287- hPID.add (" Before/h1PID_TOF_pos_pion" , " Pion PID distribution in data" , kTH1F , {{100 , -10 .0f , 10 .0f }});
288- hPID.add (" Before/h1PID_TPC_neg_kaon" , " Kaon PID distribution in data" , kTH1F , {{100 , -10 .0f , 10 .0f }});
289- hPID.add (" Before/h1PID_TPC_neg_pion" , " Pion PID distribution in data" , kTH1F , {{100 , -10 .0f , 10 .0f }});
290- hPID.add (" Before/h1PID_TOF_neg_kaon" , " Kaon PID distribution in data" , kTH1F , {{100 , -10 .0f , 10 .0f }});
291- hPID.add (" Before/h1PID_TOF_neg_pion" , " Pion PID distribution in data" , kTH1F , {{100 , -10 .0f , 10 .0f }});
284+ // hPID.add("Before/h1PID_TPC_pos_kaon", "Kaon PID distribution in data", kTH1F, {{100, -10.0f, 10.0f}});
285+ // hPID.add("Before/h1PID_TPC_pos_pion", "Pion PID distribution in data", kTH1F, {{100, -10.0f, 10.0f}});
286+ // hPID.add("Before/h1PID_TOF_pos_kaon", "Kaon PID distribution in data", kTH1F, {{100, -10.0f, 10.0f}});
287+ // hPID.add("Before/h1PID_TOF_pos_pion", "Pion PID distribution in data", kTH1F, {{100, -10.0f, 10.0f}});
288+ // hPID.add("Before/h1PID_TPC_neg_kaon", "Kaon PID distribution in data", kTH1F, {{100, -10.0f, 10.0f}});
289+ // hPID.add("Before/h1PID_TPC_neg_pion", "Pion PID distribution in data", kTH1F, {{100, -10.0f, 10.0f}});
290+ // hPID.add("Before/h1PID_TOF_neg_kaon", "Kaon PID distribution in data", kTH1F, {{100, -10.0f, 10.0f}});
291+ // hPID.add("Before/h1PID_TOF_neg_pion", "Pion PID distribution in data", kTH1F, {{100, -10.0f, 10.0f}});
292292
293293 hPID.add (" Before/hTPCnsigKa_mult_pt" , " TPC nsigma of Kaon brfore PID with pt and multiplicity" , kTH3F , {{100 , -10 .0f , 10 .0f }, multiplicityAxis, ptAxis});
294294 hPID.add (" Before/hTPCnsigPi_mult_pt" , " TPC nsigma of Pion brfore PID with pt and multiplicity" , kTH3F , {{100 , -10 .0f , 10 .0f }, multiplicityAxis, ptAxis});
@@ -1029,10 +1029,10 @@ struct Kstarqa {
10291029 continue ;
10301030
10311031 if (cQAplots) {
1032- hPID.fill (HIST (" Before/hNsigmaTPC_Ka_before" ), track1.pt (), track1.tpcNSigmaKa ());
1033- hPID.fill (HIST (" Before/hNsigmaTOF_Ka_before" ), track1.pt (), track1.tofNSigmaKa ());
1034- hPID.fill (HIST (" Before/hNsigmaTPC_Pi_before" ), track2.pt (), track2.tpcNSigmaPi ());
1035- hPID.fill (HIST (" Before/hNsigmaTOF_Pi_before" ), track2.pt (), track2.tofNSigmaPi ());
1032+ // hPID.fill(HIST("Before/hNsigmaTPC_Ka_before"), track1.pt(), track1.tpcNSigmaKa());
1033+ // hPID.fill(HIST("Before/hNsigmaTOF_Ka_before"), track1.pt(), track1.tofNSigmaKa());
1034+ // hPID.fill(HIST("Before/hNsigmaTPC_Pi_before"), track2.pt(), track2.tpcNSigmaPi());
1035+ // hPID.fill(HIST("Before/hNsigmaTOF_Pi_before"), track2.pt(), track2.tofNSigmaPi());
10361036 hPID.fill (HIST (" Before/hNsigma_TPC_TOF_Ka_before" ), track1.tpcNSigmaKa (), track1.tofNSigmaKa ());
10371037 hPID.fill (HIST (" Before/hNsigma_TPC_TOF_Pi_before" ), track2.tpcNSigmaPi (), track2.tofNSigmaPi ());
10381038
@@ -1042,20 +1042,20 @@ struct Kstarqa {
10421042 hPID.fill (HIST (" Before/hTOFnsigPi_mult_pt" ), track2.tofNSigmaPi (), multiplicity, track2.pt ());
10431043
10441044 hOthers.fill (HIST (" hCRFC_before" ), track1.tpcCrossedRowsOverFindableCls ());
1045- hOthers.fill (HIST (" dE_by_dx_TPC" ), track1.p (), track1.tpcSignal ());
1045+ // hOthers.fill(HIST("dE_by_dx_TPC"), track1.p(), track1.tpcSignal());
10461046 hOthers.fill (HIST (" hphi" ), track1.phi ());
10471047
1048- if (track1.sign () < 0 ) {
1049- hPID.fill (HIST (" Before/h1PID_TPC_neg_kaon" ), track1.tpcNSigmaKa ());
1050- hPID.fill (HIST (" Before/h1PID_TPC_neg_pion" ), track2.tpcNSigmaPi ());
1051- hPID.fill (HIST (" Before/h1PID_TOF_neg_kaon" ), track1.tofNSigmaKa ());
1052- hPID.fill (HIST (" Before/h1PID_TOF_neg_pion" ), track2.tofNSigmaPi ());
1053- } else {
1054- hPID.fill (HIST (" Before/h1PID_TPC_pos_kaon" ), track1.tpcNSigmaKa ());
1055- hPID.fill (HIST (" Before/h1PID_TPC_pos_pion" ), track2.tpcNSigmaPi ());
1056- hPID.fill (HIST (" Before/h1PID_TOF_pos_kaon" ), track1.tofNSigmaKa ());
1057- hPID.fill (HIST (" Before/h1PID_TOF_pos_pion" ), track2.tofNSigmaPi ());
1058- }
1048+ // if (track1.sign() < 0) {
1049+ // hPID.fill(HIST("Before/h1PID_TPC_neg_kaon"), track1.tpcNSigmaKa());
1050+ // hPID.fill(HIST("Before/h1PID_TPC_neg_pion"), track2.tpcNSigmaPi());
1051+ // hPID.fill(HIST("Before/h1PID_TOF_neg_kaon"), track1.tofNSigmaKa());
1052+ // hPID.fill(HIST("Before/h1PID_TOF_neg_pion"), track2.tofNSigmaPi());
1053+ // } else {
1054+ // hPID.fill(HIST("Before/h1PID_TPC_pos_kaon"), track1.tpcNSigmaKa());
1055+ // hPID.fill(HIST("Before/h1PID_TPC_pos_pion"), track2.tpcNSigmaPi());
1056+ // hPID.fill(HIST("Before/h1PID_TOF_pos_kaon"), track1.tofNSigmaKa());
1057+ // hPID.fill(HIST("Before/h1PID_TOF_pos_pion"), track2.tofNSigmaPi());
1058+ // }
10591059 }
10601060
10611061 if (cQAevents) {
@@ -1114,10 +1114,10 @@ struct Kstarqa {
11141114 hPID.fill (HIST (" After/hTOFnsigPi_mult_pt" ), track2.tofNSigmaPi (), multiplicity, track2.pt ());
11151115 hOthers.fill (HIST (" hEta_after" ), track1.eta ());
11161116 hOthers.fill (HIST (" hCRFC_after" ), track1.tpcCrossedRowsOverFindableCls ());
1117- hPID.fill (HIST (" After/hNsigmaKaonTPC_after" ), track1.pt (), track1.tpcNSigmaKa ());
1118- hPID.fill (HIST (" After/hNsigmaKaonTOF_after" ), track1.pt (), track1.tofNSigmaKa ());
1119- hPID.fill (HIST (" After/hNsigmaPionTPC_after" ), track2.pt (), track2.tpcNSigmaPi ());
1120- hPID.fill (HIST (" After/hNsigmaPionTOF_after" ), track2.pt (), track2.tofNSigmaPi ());
1117+ // hPID.fill(HIST("After/hNsigmaKaonTPC_after"), track1.pt(), track1.tpcNSigmaKa());
1118+ // hPID.fill(HIST("After/hNsigmaKaonTOF_after"), track1.pt(), track1.tofNSigmaKa());
1119+ // hPID.fill(HIST("After/hNsigmaPionTPC_after"), track2.pt(), track2.tpcNSigmaPi());
1120+ // hPID.fill(HIST("After/hNsigmaPionTOF_after"), track2.pt(), track2.tofNSigmaPi());
11211121 hPID.fill (HIST (" After/hNsigma_TPC_TOF_Ka_after" ), track1.tpcNSigmaKa (), track1.tofNSigmaKa ());
11221122 hPID.fill (HIST (" After/hNsigma_TPC_TOF_Pi_after" ), track2.tpcNSigmaPi (), track2.tofNSigmaPi ());
11231123 }
@@ -1407,20 +1407,20 @@ struct Kstarqa {
14071407 hPID.fill (HIST (" Before/hTOFnsigPi_mult_pt" ), track2.tofNSigmaPi (), multiplicity, track2.pt ());
14081408
14091409 hOthers.fill (HIST (" hCRFC_before" ), track1.tpcCrossedRowsOverFindableCls ());
1410- hOthers.fill (HIST (" dE_by_dx_TPC" ), track1.p (), track1.tpcSignal ());
1410+ // hOthers.fill(HIST("dE_by_dx_TPC"), track1.p(), track1.tpcSignal());
14111411 hOthers.fill (HIST (" hphi" ), track1.phi ());
14121412
1413- if (track1.sign () < 0 ) {
1414- hPID.fill (HIST (" Before/h1PID_TPC_neg_kaon" ), track1.tpcNSigmaKa ());
1415- hPID.fill (HIST (" Before/h1PID_TPC_neg_pion" ), track2.tpcNSigmaPi ());
1416- hPID.fill (HIST (" Before/h1PID_TOF_neg_kaon" ), track1.tofNSigmaKa ());
1417- hPID.fill (HIST (" Before/h1PID_TOF_neg_pion" ), track2.tofNSigmaPi ());
1418- } else {
1419- hPID.fill (HIST (" Before/h1PID_TPC_pos_kaon" ), track1.tpcNSigmaKa ());
1420- hPID.fill (HIST (" Before/h1PID_TPC_pos_pion" ), track2.tpcNSigmaPi ());
1421- hPID.fill (HIST (" Before/h1PID_TOF_pos_kaon" ), track1.tofNSigmaKa ());
1422- hPID.fill (HIST (" Before/h1PID_TOF_pos_pion" ), track2.tofNSigmaPi ());
1423- }
1413+ // if (track1.sign() < 0) {
1414+ // hPID.fill(HIST("Before/h1PID_TPC_neg_kaon"), track1.tpcNSigmaKa());
1415+ // hPID.fill(HIST("Before/h1PID_TPC_neg_pion"), track2.tpcNSigmaPi());
1416+ // hPID.fill(HIST("Before/h1PID_TOF_neg_kaon"), track1.tofNSigmaKa());
1417+ // hPID.fill(HIST("Before/h1PID_TOF_neg_pion"), track2.tofNSigmaPi());
1418+ // } else {
1419+ // hPID.fill(HIST("Before/h1PID_TPC_pos_kaon"), track1.tpcNSigmaKa());
1420+ // hPID.fill(HIST("Before/h1PID_TPC_pos_pion"), track2.tpcNSigmaPi());
1421+ // hPID.fill(HIST("Before/h1PID_TOF_pos_kaon"), track1.tofNSigmaKa());
1422+ // hPID.fill(HIST("Before/h1PID_TOF_pos_pion"), track2.tofNSigmaPi());
1423+ // }
14241424 }
14251425
14261426 if (cQAevents) {
@@ -1914,26 +1914,26 @@ struct Kstarqa {
19141914 int track2PDG = std::abs (mctrack2.pdgCode ());
19151915
19161916 if (cQAplots && (mctrack2.pdgCode () == PDG_t::kPiPlus )) { // pion
1917- hPID.fill (HIST (" Before/h1PID_TPC_pos_pion" ), track2.tpcNSigmaPi ());
1918- hPID.fill (HIST (" Before/h1PID_TOF_pos_pion" ), track2.tofNSigmaPi ());
1917+ // hPID.fill(HIST("Before/h1PID_TPC_pos_pion"), track2.tpcNSigmaPi());
1918+ // hPID.fill(HIST("Before/h1PID_TOF_pos_pion"), track2.tofNSigmaPi());
19191919 hPID.fill (HIST (" Before/hNsigmaTPC_Pi_before" ), track2.pt (), track2.tpcNSigmaPi ());
19201920 hPID.fill (HIST (" Before/hNsigmaTOF_Pi_before" ), track2.pt (), track2.tofNSigmaPi ());
19211921 }
19221922 if (cQAplots && (mctrack2.pdgCode () == PDG_t::kKPlus )) { // kaon
1923- hPID.fill (HIST (" Before/h1PID_TPC_pos_kaon" ), track2.tpcNSigmaKa ());
1924- hPID.fill (HIST (" Before/h1PID_TOF_pos_kaon" ), track2.tofNSigmaKa ());
1923+ // hPID.fill(HIST("Before/h1PID_TPC_pos_kaon"), track2.tpcNSigmaKa());
1924+ // hPID.fill(HIST("Before/h1PID_TOF_pos_kaon"), track2.tofNSigmaKa());
19251925 hPID.fill (HIST (" Before/hNsigmaTPC_Ka_before" ), track2.pt (), track2.tpcNSigmaKa ());
19261926 hPID.fill (HIST (" Before/hNsigmaTOF_Ka_before" ), track2.pt (), track2.tofNSigmaKa ());
19271927 }
19281928 if (cQAplots && (mctrack2.pdgCode () == PDG_t::kPiMinus )) { // negative track pion
1929- hPID.fill (HIST (" Before/h1PID_TPC_neg_pion" ), track2.tpcNSigmaPi ());
1930- hPID.fill (HIST (" Before/h1PID_TOF_neg_pion" ), track2.tofNSigmaPi ());
1929+ // hPID.fill(HIST("Before/h1PID_TPC_neg_pion"), track2.tpcNSigmaPi());
1930+ // hPID.fill(HIST("Before/h1PID_TOF_neg_pion"), track2.tofNSigmaPi());
19311931 hPID.fill (HIST (" Before/hNsigmaTPC_Pi_before" ), track2.pt (), track2.tpcNSigmaPi ());
19321932 hPID.fill (HIST (" Before/hNsigmaTOF_Pi_before" ), track2.pt (), track2.tofNSigmaPi ());
19331933 }
19341934 if (cQAplots && (mctrack2.pdgCode () == PDG_t::kKMinus )) { // negative track kaon
1935- hPID.fill (HIST (" Before/h1PID_TPC_neg_kaon" ), track2.tpcNSigmaKa ());
1936- hPID.fill (HIST (" Before/h1PID_TOF_neg_kaon" ), track2.tofNSigmaKa ());
1935+ // hPID.fill(HIST("Before/h1PID_TPC_neg_kaon"), track2.tpcNSigmaKa());
1936+ // hPID.fill(HIST("Before/h1PID_TOF_neg_kaon"), track2.tofNSigmaKa());
19371937 hPID.fill (HIST (" Before/hNsigmaTPC_Ka_before" ), track2.pt (), track2.tpcNSigmaKa ());
19381938 hPID.fill (HIST (" Before/hNsigmaTOF_Ka_before" ), track2.pt (), track2.tofNSigmaKa ());
19391939 }
@@ -2186,31 +2186,31 @@ struct Kstarqa {
21862186 int track1PDG = std::abs (mctrack1.pdgCode ());
21872187 int track2PDG = std::abs (mctrack2.pdgCode ());
21882188
2189- if (cQAplots && (mctrack2.pdgCode () == PDG_t::kPiPlus )) { // pion
2190- hPID.fill (HIST (" Before/h1PID_TPC_pos_pion" ), track2.tpcNSigmaPi ());
2191- hPID.fill (HIST (" Before/h1PID_TOF_pos_pion" ), track2.tofNSigmaPi ());
2192- hPID.fill (HIST (" Before/hNsigmaTPC_Pi_before" ), track2.pt (), track2.tpcNSigmaPi ());
2193- hPID.fill (HIST (" Before/hNsigmaTOF_Pi_before" ), track2.pt (), track2.tofNSigmaPi ());
2194- }
2195- if (cQAplots && (mctrack2.pdgCode () == PDG_t::kKPlus )) { // kaon
2196- hPID.fill (HIST (" Before/h1PID_TPC_pos_kaon" ), track2.tpcNSigmaKa ());
2197- hPID.fill (HIST (" Before/h1PID_TOF_pos_kaon" ), track2.tofNSigmaKa ());
2198- hPID.fill (HIST (" Before/hNsigmaTPC_Ka_before" ), track2.pt (), track2.tpcNSigmaKa ());
2199- hPID.fill (HIST (" Before/hNsigmaTOF_Ka_before" ), track2.pt (), track2.tofNSigmaKa ());
2200- }
2201- if (cQAplots && (mctrack2.pdgCode () == PDG_t::kPiMinus )) { // negative track pion
2202- hPID.fill (HIST (" Before/h1PID_TPC_neg_pion" ), track2.tpcNSigmaPi ());
2203- hPID.fill (HIST (" Before/h1PID_TOF_neg_pion" ), track2.tofNSigmaPi ());
2189+ if (cQAplots && (track2PDG == PDG_t::kPiPlus )) { // pion
2190+ // hPID.fill(HIST("Before/h1PID_TPC_pos_pion"), track2.tpcNSigmaPi());
2191+ // hPID.fill(HIST("Before/h1PID_TOF_pos_pion"), track2.tofNSigmaPi());
22042192 hPID.fill (HIST (" Before/hNsigmaTPC_Pi_before" ), track2.pt (), track2.tpcNSigmaPi ());
22052193 hPID.fill (HIST (" Before/hNsigmaTOF_Pi_before" ), track2.pt (), track2.tofNSigmaPi ());
22062194 }
2207- if (cQAplots && (mctrack2. pdgCode () == PDG_t::kKMinus )) { // negative track kaon
2208- hPID.fill (HIST (" Before/h1PID_TPC_neg_kaon " ), track2.tpcNSigmaKa ());
2209- hPID.fill (HIST (" Before/h1PID_TOF_neg_kaon " ), track2.tofNSigmaKa ());
2195+ if (cQAplots && (track2PDG == PDG_t::kKPlus )) { // kaon
2196+ // hPID.fill(HIST("Before/h1PID_TPC_pos_kaon "), track2.tpcNSigmaKa());
2197+ // hPID.fill(HIST("Before/h1PID_TOF_pos_kaon "), track2.tofNSigmaKa());
22102198 hPID.fill (HIST (" Before/hNsigmaTPC_Ka_before" ), track2.pt (), track2.tpcNSigmaKa ());
22112199 hPID.fill (HIST (" Before/hNsigmaTOF_Ka_before" ), track2.pt (), track2.tofNSigmaKa ());
22122200 }
2213- if (cQAplots && (std::abs (mctrack1.pdgCode ()) == PDG_t::kKPlus && std::abs (mctrack2.pdgCode ()) == PDG_t::kPiPlus )) {
2201+ // if (cQAplots && (mctrack2.pdgCode() == PDG_t::kPiMinus)) { // negative track pion
2202+ // // hPID.fill(HIST("Before/h1PID_TPC_neg_pion"), track2.tpcNSigmaPi());
2203+ // // hPID.fill(HIST("Before/h1PID_TOF_neg_pion"), track2.tofNSigmaPi());
2204+ // hPID.fill(HIST("Before/hNsigmaTPC_Pi_before"), track2.pt(), track2.tpcNSigmaPi());
2205+ // hPID.fill(HIST("Before/hNsigmaTOF_Pi_before"), track2.pt(), track2.tofNSigmaPi());
2206+ // }
2207+ // if (cQAplots && (mctrack2.pdgCode() == PDG_t::kKMinus)) { // negative track kaon
2208+ // // hPID.fill(HIST("Before/h1PID_TPC_neg_kaon"), track2.tpcNSigmaKa());
2209+ // // hPID.fill(HIST("Before/h1PID_TOF_neg_kaon"), track2.tofNSigmaKa());
2210+ // hPID.fill(HIST("Before/hNsigmaTPC_Ka_before"), track2.pt(), track2.tpcNSigmaKa());
2211+ // hPID.fill(HIST("Before/hNsigmaTOF_Ka_before"), track2.pt(), track2.tofNSigmaKa());
2212+ // }
2213+ if (cQAplots && track1PDG == PDG_t::kKPlus && track2PDG == PDG_t::kPiPlus ) {
22142214 hPID.fill (HIST (" Before/hNsigma_TPC_TOF_Ka_before" ), track1.tpcNSigmaKa (), track1.tofNSigmaKa ());
22152215 hPID.fill (HIST (" Before/hNsigma_TPC_TOF_Pi_before" ), track2.tpcNSigmaPi (), track2.tofNSigmaPi ());
22162216 hPID.fill (HIST (" Before/hTPCnsigKa_mult_pt" ), track1.tpcNSigmaKa (), multiplicity, track1.pt ());
@@ -2393,20 +2393,20 @@ struct Kstarqa {
23932393 hPID.fill (HIST (" Before/hTOFnsigPi_mult_pt" ), track2.tofNSigmaPi (), multiplicity, track2.pt ());
23942394
23952395 hOthers.fill (HIST (" hCRFC_before" ), track1.tpcCrossedRowsOverFindableCls ());
2396- hOthers.fill (HIST (" dE_by_dx_TPC" ), track1.p (), track1.tpcSignal ());
2396+ // hOthers.fill(HIST("dE_by_dx_TPC"), track1.p(), track1.tpcSignal());
23972397 hOthers.fill (HIST (" hphi" ), track1.phi ());
23982398
2399- if (track1.sign () < 0 ) {
2400- hPID.fill (HIST (" Before/h1PID_TPC_neg_kaon" ), track1.tpcNSigmaKa ());
2401- // hPID.fill(HIST("Before/h1PID_TPC_neg_pion"), track2.tpcNSigmaPi());
2402- hPID.fill (HIST (" Before/h1PID_TOF_neg_kaon" ), track1.tofNSigmaKa ());
2403- // hPID.fill(HIST("Before/h1PID_TOF_neg_pion"), track2.tofNSigmaPi());
2404- } else {
2405- hPID.fill (HIST (" Before/h1PID_TPC_pos_kaon" ), track1.tpcNSigmaKa ());
2406- // hPID.fill(HIST("Before/h1PID_TPC_pos_pion"), track2.tpcNSigmaPi());
2407- hPID.fill (HIST (" Before/h1PID_TOF_pos_kaon" ), track1.tofNSigmaKa ());
2408- // hPID.fill(HIST("Before/h1PID_TOF_pos_pion"), track2.tofNSigmaPi());
2409- }
2399+ // if (track1.sign() < 0) {
2400+ // hPID.fill(HIST("Before/h1PID_TPC_neg_kaon"), track1.tpcNSigmaKa());
2401+ // // hPID.fill(HIST("Before/h1PID_TPC_neg_pion"), track2.tpcNSigmaPi());
2402+ // hPID.fill(HIST("Before/h1PID_TOF_neg_kaon"), track1.tofNSigmaKa());
2403+ // // hPID.fill(HIST("Before/h1PID_TOF_neg_pion"), track2.tofNSigmaPi());
2404+ // } else {
2405+ // hPID.fill(HIST("Before/h1PID_TPC_pos_kaon"), track1.tpcNSigmaKa());
2406+ // // hPID.fill(HIST("Before/h1PID_TPC_pos_pion"), track2.tpcNSigmaPi());
2407+ // hPID.fill(HIST("Before/h1PID_TOF_pos_kaon"), track1.tofNSigmaKa());
2408+ // // hPID.fill(HIST("Before/h1PID_TOF_pos_pion"), track2.tofNSigmaPi());
2409+ // }
24102410 }
24112411
24122412 if (cQAevents) {
@@ -2489,8 +2489,10 @@ struct Kstarqa {
24892489 continue;
24902490 } */
24912491
2492- hOthers.fill (HIST (" hKstar_rap_pt" ), mother.Rapidity (), mother.Pt ());
2493- hOthers.fill (HIST (" hKstar_eta_pt" ), mother.Eta (), mother.Pt ());
2492+ if (cQAplots) {
2493+ hOthers.fill (HIST (" hKstar_rap_pt" ), mother.Rapidity (), mother.Pt ());
2494+ hOthers.fill (HIST (" hKstar_eta_pt" ), mother.Eta (), mother.Pt ());
2495+ }
24942496
24952497 isMix = false ;
24962498 fillInvMass (daughter1, daughter2, mother, multiplicity, isMix, track1, track2);
@@ -2799,25 +2801,25 @@ struct Kstarqa {
27992801 // hPID.fill(HIST("Before/hNsigmaTPC_Pi_before"), track2.pt(), track2.tpcNSigmaPi());
28002802 // hPID.fill(HIST("Before/hNsigmaTOF_Pi_before"), track2.pt(), track2.tofNSigmaPi());
28012803 // }
2802- if (cQAplots && (mctrack2.pdgCode () == PDG_t::kKPlus )) { // kaon
2803- hPID.fill (HIST (" Before/h1PID_TPC_pos_kaon" ), track2.tpcNSigmaKa ());
2804- hPID.fill (HIST (" Before/h1PID_TOF_pos_kaon" ), track2.tofNSigmaKa ());
2805- hPID.fill (HIST (" Before/hNsigmaTPC_Ka_before" ), track2.pt (), track2.tpcNSigmaKa ());
2806- hPID.fill (HIST (" Before/hNsigmaTOF_Ka_before" ), track2.pt (), track2.tofNSigmaKa ());
2807- }
2808- // if (cQAplots && (mctrack2.pdgCode() == - PDG_t::kPiMinus)) { // negative track pion
2804+ // if (cQAplots && (mctrack2.pdgCode() == PDG_t::kKPlus)) { // kaon
2805+ // // hPID.fill(HIST("Before/h1PID_TPC_pos_kaon"), track2.tpcNSigmaKa());
2806+ // // hPID.fill(HIST("Before/h1PID_TOF_pos_kaon"), track2.tofNSigmaKa());
2807+ // hPID.fill(HIST("Before/hNsigmaTPC_Ka_before"), track2.pt(), track2.tpcNSigmaKa());
2808+ // hPID.fill(HIST("Before/hNsigmaTOF_Ka_before"), track2.pt(), track2.tofNSigmaKa());
2809+ // }
2810+ // if (cQAplots && (mctrack2.pdgCode() == PDG_t::kPiMinus)) { // negative track pion
28092811 // hPID.fill(HIST("Before/h1PID_TPC_neg_pion"), track2.tpcNSigmaPi());
28102812 // hPID.fill(HIST("Before/h1PID_TOF_neg_pion"), track2.tofNSigmaPi());
28112813 // hPID.fill(HIST("Before/hNsigmaTPC_Pi_before"), track2.pt(), track2.tpcNSigmaPi());
28122814 // hPID.fill(HIST("Before/hNsigmaTOF_Pi_before"), track2.pt(), track2.tofNSigmaPi());
28132815 // }
2814- if (cQAplots && (mctrack2.pdgCode () == - PDG_t::kKMinus )) { // negative track kaon
2815- hPID.fill (HIST (" Before/h1PID_TPC_neg_kaon" ), track2.tpcNSigmaKa ());
2816- hPID.fill (HIST (" Before/h1PID_TOF_neg_kaon" ), track2.tofNSigmaKa ());
2817- hPID.fill (HIST (" Before/hNsigmaTPC_Ka_before" ), track2.pt (), track2.tpcNSigmaKa ());
2818- hPID.fill (HIST (" Before/hNsigmaTOF_Ka_before" ), track2.pt (), track2.tofNSigmaKa ());
2819- }
2820- if (cQAplots && (std::abs (mctrack1. pdgCode ()) == PDG_t::kKPlus && std::abs (mctrack2. pdgCode ()) == PDG_t::kKMinus )) {
2816+ // if (cQAplots && (mctrack2.pdgCode() == PDG_t::kKMinus)) { // negative track kaon
2817+ // // hPID.fill(HIST("Before/h1PID_TPC_neg_kaon"), track2.tpcNSigmaKa());
2818+ // // hPID.fill(HIST("Before/h1PID_TOF_neg_kaon"), track2.tofNSigmaKa());
2819+ // hPID.fill(HIST("Before/hNsigmaTPC_Ka_before"), track2.pt(), track2.tpcNSigmaKa());
2820+ // hPID.fill(HIST("Before/hNsigmaTOF_Ka_before"), track2.pt(), track2.tofNSigmaKa());
2821+ // }
2822+ if (cQAplots && (track1PDG == PDG_t::kKPlus && track2PDG == PDG_t::kKMinus )) {
28212823 hPID.fill (HIST (" Before/hNsigma_TPC_TOF_Ka_before" ), track1.tpcNSigmaKa (), track1.tofNSigmaKa ());
28222824 // hPID.fill(HIST("Before/hNsigma_TPC_TOF_Pi_before"), track2.tpcNSigmaPi(), track2.tofNSigmaPi());
28232825 }
0 commit comments