@@ -1222,34 +1222,33 @@ struct femtoUniversePairTaskTrackCascadeExtended {
12221222 int pdgCode = static_cast <int >(part.pidCut ());
12231223 if ((confCascType1 == 0 && pdgCode != kOmegaMinus ) || (confCascType1 == 2 && pdgCode != kOmegaPlusBar ) || (confCascType1 == 1 && pdgCode != kXiMinus ) || (confCascType1 == 3 && pdgCode != kXiPlusBar ))
12241224 continue ;
1225-
12261225 cascQAHistos.fillQA <false , false >(part);
1226+ }
12271227
1228- for (const auto & part : groupPartsOne) {
1229- int pdgCode = static_cast <int >(part.pidCut ());
1230- if (pdgCode != trackparticleconfigs.confTrkPDGCodePartOne )
1231- continue ;
1232- const auto & pdgTrackParticle = pdgMC->GetParticle (pdgCode);
1233- if (!pdgTrackParticle) {
1234- continue ;
1235- }
1236-
1237- if (pdgTrackParticle->Charge () > 0 ) {
1238- trackHistoPartOnePos.fillQA <false , false >(part);
1239- } else if (pdgTrackParticle->Charge () < 0 ) {
1240- trackHistoPartOneNeg.fillQA <false , false >(part);
1241- }
1228+ for (const auto & part : groupPartsOne) {
1229+ int pdgCode = static_cast <int >(part.pidCut ());
1230+ if (pdgCode != trackparticleconfigs.confTrkPDGCodePartOne )
1231+ continue ;
1232+ const auto & pdgTrackParticle = pdgMC->GetParticle (pdgCode);
1233+ if (!pdgTrackParticle) {
1234+ continue ;
12421235 }
12431236
1244- for (const auto & [p1, p2] : combinations (CombinationsFullIndexPolicy (groupPartsOne, groupPartsTwo))) {
1245- if (static_cast <int >(p1.pidCut ()) != trackparticleconfigs.confTrkPDGCodePartOne )
1246- continue ;
1247- int pdgCodeCasc = static_cast <int >(p2.pidCut ());
1248- if ((confCascType1 == 0 && pdgCodeCasc != kOmegaMinus ) || (confCascType1 == 2 && pdgCodeCasc != kOmegaPlusBar ) || (confCascType1 == 1 && pdgCodeCasc != kXiMinus ) || (confCascType1 == 3 && pdgCodeCasc != kXiPlusBar ))
1249- continue ;
1250- sameEventCont.setPair <false >(p1, p2, multCol, confUse3D, 1 .0f );
1237+ if (pdgTrackParticle->Charge () > 0 ) {
1238+ trackHistoPartOnePos.fillQA <false , false >(part);
1239+ } else if (pdgTrackParticle->Charge () < 0 ) {
1240+ trackHistoPartOneNeg.fillQA <false , false >(part);
12511241 }
12521242 }
1243+
1244+ for (const auto & [p1, p2] : combinations (CombinationsFullIndexPolicy (groupPartsOne, groupPartsTwo))) {
1245+ if (static_cast <int >(p1.pidCut ()) != trackparticleconfigs.confTrkPDGCodePartOne )
1246+ continue ;
1247+ int pdgCodeCasc = static_cast <int >(p2.pidCut ());
1248+ if ((confCascType1 == 0 && pdgCodeCasc != kOmegaMinus ) || (confCascType1 == 2 && pdgCodeCasc != kOmegaPlusBar ) || (confCascType1 == 1 && pdgCodeCasc != kXiMinus ) || (confCascType1 == 3 && pdgCodeCasc != kXiPlusBar ))
1249+ continue ;
1250+ sameEventCont.setPair <false >(p1, p2, multCol, confUse3D, 1 .0f );
1251+ }
12531252 }
12541253 PROCESS_SWITCH (femtoUniversePairTaskTrackCascadeExtended, processSameEventMCgen, " Enable processing same event MC truth for track - cascade" , false );
12551254
0 commit comments