@@ -3442,26 +3442,26 @@ struct AnalysisDileptonTrackTrack {
34423442 int indexLepton2 = dilepton.index1Id ();
34433443 auto lepton1 = tracks.rawIteratorAt (dilepton.index0Id ());
34443444 auto lepton2 = tracks.rawIteratorAt (dilepton.index1Id ());
3445- auto lepton1MC = lepton1.reducedMCTrack ();
3446- auto lepton2MC = lepton2.reducedMCTrack ();
3445+ auto lepton1MC = lepton1.mcParticle ();
3446+ auto lepton2MC = lepton2.mcParticle ();
34473447 // Check that the dilepton has zero charge
34483448 if (dilepton.sign () != 0 ) {
34493449 continue ;
34503450 }
34513451 VarManager::FillTrack<fgDileptonFillMap>(dilepton, fValuesQuadruplet );
34523452
3453- bool isAmbiguousLepton = (dilepton.filterMap_raw () & (static_cast <uint32_t >(1 ) << 28 )) ||
3454- (dilepton.filterMap_raw () & (static_cast <uint32_t >(1 ) << 29 )) ||
3455- (dilepton.filterMap_raw () & (static_cast <uint32_t >(1 ) << 30 )) ||
3456- (dilepton.filterMap_raw () & (static_cast <uint32_t >(1 ) << 31 ));
3457- // if (isAmbi && isAmbiguousLepton)
3458- // continue; // skip ambiguous dileptons
3459- if constexpr ((TTrackFillMap & VarManager::ObjTypes::AmbiTrack) > 0 ) {
3460- if (isAmbiguousLepton) {
3461- // LOGP(info, "Dilepton {} is ambiguous, skipping", dilepton.index0Id());
3462- continue ; // skip ambiguous dileptons
3463- }
3464- }
3453+ // bool isAmbiguousLepton = (dilepton.filterMap_raw() & (static_cast<uint32_t>(1) << 28)) ||
3454+ // (dilepton.filterMap_raw() & (static_cast<uint32_t>(1) << 29)) ||
3455+ // (dilepton.filterMap_raw() & (static_cast<uint32_t>(1) << 30)) ||
3456+ // (dilepton.filterMap_raw() & (static_cast<uint32_t>(1) << 31));
3457+ // // if (isAmbi && isAmbiguousLepton)
3458+ // // continue; // skip ambiguous dileptons
3459+ // if constexpr ((TTrackFillMap & VarManager::ObjTypes::AmbiTrack) > 0) {
3460+ // if (isAmbiguousLepton) {
3461+ // // LOGP(info, "Dilepton {} is ambiguous, skipping", dilepton.index0Id());
3462+ // continue; // skip ambiguous dileptons
3463+ // }
3464+ // }
34653465
34663466 // apply the dilepton cut
34673467 if (!fDileptonCut .IsSelected (fValuesQuadruplet ))
@@ -3479,8 +3479,8 @@ struct AnalysisDileptonTrackTrack {
34793479 }
34803480
34813481 // get the track from this association
3482- auto track1 = a1.template reducedtrack_as <TTracks>( );
3483- auto track2 = a2.template reducedtrack_as <TTracks>( );
3482+ auto track1 = tracks. rawIteratorAt ( a1.trackId () );
3483+ auto track2 = tracks. rawIteratorAt ( a2.trackId () );
34843484 if (track1.sign () < 0 && track2.sign () > 0 ) {
34853485 std::swap (track1, track2);
34863486 }
@@ -3511,8 +3511,8 @@ struct AnalysisDileptonTrackTrack {
35113511 // VarManager::FillBdtScore(fOutputMlQuad, fValuesQuadruplet);
35123512 // }
35133513
3514- auto track1MC = track1.reducedMCTrack ();
3515- auto track2MC = track2.reducedMCTrack ();
3514+ auto track1MC = track1.mcParticle ();
3515+ auto track2MC = track2.mcParticle ();
35163516 mcDecision = 0 ;
35173517 isig = 0 ;
35183518 for (auto sig = fRecMCSignals .begin (); sig != fRecMCSignals .end (); sig++, isig++) {
0 commit comments