@@ -94,7 +94,6 @@ struct HfCorrelatorHfeHadrons {
9494
9595 {" hptElectron" , " hptElectron" , {HistType::kTH1F , {{100 , 0 , 100 }}}},
9696
97-
9897 {" hMixEventInclusiveEHCorrl" , " Sparse for mix event Delta phi and Delta eta Inclusive Electron with Hadron;p_{T}^{e} (GeV#it{/c});p_{T}^{h} (GeV#it{/c});#Delta#varphi;#Delta#eta;" , hCorrelSpec},
9998 {" hMixEventLSEHCorrel" , " Sparse for mix event Delta phi and Delta eta Like sign Electron pair with Hadron;p_{T}^{e} (GeV#it{/c});p_{T}^{h} (GeV#it{/c});#Delta#varphi;#Delta#eta;" , hCorrelSpec},
10099 {" hMixEventULSEHCorrel" , " Sparse for mix event Delta phi and Delta eta Unlike sign Electron pair with Hadron;p_{T}^{e} (GeV#it{/c});p_{T}^{h} (GeV#it{/c});#Delta#varphi;#Delta#eta;" , hCorrelSpec},
@@ -383,7 +382,6 @@ struct HfCorrelatorHfeHadrons {
383382 auto tracksTuple = std::make_tuple (electrons, mcParticles);
384383 Pair<McGenTableCollisions, aod::HfMcGenSelEl, aod::McParticles, BinningTypeMcGen> pairMcGen{corrBinningMcGen, 5 , -1 , mcCollision, tracksTuple, &cache};
385384
386-
387385 // loop over the rows of the new table
388386 double ptElectronMix = -999 ;
389387 double phiElectronMix = -999 ;
@@ -423,50 +421,48 @@ struct HfCorrelatorHfeHadrons {
423421 registry.fill (HIST (" hMixEventMcGenInclusiveEHCorrl" ), ptElectronMix, ptHadronMix, deltaPhiMix, deltaEtaMix);
424422 }
425423
426-
427- // loop over the rows of the new table
428- double ptElectronMix = -999 ;
429- double phiElectronMix = -999 ;
430- double etaElectronMix = -999 ;
431- double deltaPhiMix = -999 ;
432- double deltaEtaMix = -999 ;
433- double ptHadronMix = -999 ;
434- double etaHadronMix = -999 ;
435- double phiHadronMix = -999 ;
436- for (const auto & [c1, tracks1, c2, tracks2] : pairMcGen) {
437- int poolBin = corrBinningMcGen.getBin (std::make_tuple (c1.posZ (), c1.multMCFT0A ()));
438- for (const auto & [t1, t2] : combinations (CombinationsFullIndexPolicy (tracks1, tracks2))) {
439- ptHadronMix = t2.pt ();
440- ptElectronMix = t1.ptTrackMc ();
441- phiElectronMix = t1.phiTrackMc ();
442- phiHadronMix = t2.phi ();
443- etaElectronMix = t1.etaTrackMc ();
444- etaHadronMix = t2.eta ();
445- if (t2.eta () < etaTrackMin || t2.eta () > etaTrackMax) {
446- continue ;
447- }
448- if (t2.pt () < ptTrackMin) {
449- continue ;
450- }
451- if (ptCondition && (ptElectronMix < ptHadronMix)) {
452- continue ;
453- }
454-
455- deltaPhiMix = RecoDecay::constrainAngle (phiElectronMix - phiHadronMix, -o2::constants::math::PIHalf);
456- deltaEtaMix = etaElectronMix - etaHadronMix;
457- bool isNonHfeCorr = false ;
458- if (t1.isNonHfeMc ()) {
459- isNonHfeCorr = true ;
460- registry.fill (HIST (" hMixEventMcGenNonHfEHCorrl" ), ptElectronMix, ptHadronMix, deltaPhiMix, deltaEtaMix);
461- } else {
462-
463- registry.fill (HIST (" hMixEventMcGenInclusiveEHCorrl" ), ptElectronMix, ptHadronMix, deltaPhiMix, deltaEtaMix);
424+ // loop over the rows of the new table
425+ double ptElectronMix = -999 ;
426+ double phiElectronMix = -999 ;
427+ double etaElectronMix = -999 ;
428+ double deltaPhiMix = -999 ;
429+ double deltaEtaMix = -999 ;
430+ double ptHadronMix = -999 ;
431+ double etaHadronMix = -999 ;
432+ double phiHadronMix = -999 ;
433+ for (const auto & [c1, tracks1, c2, tracks2] : pairMcGen) {
434+ int poolBin = corrBinningMcGen.getBin (std::make_tuple (c1.posZ (), c1.multMCFT0A ()));
435+ for (const auto & [t1, t2] : combinations (CombinationsFullIndexPolicy (tracks1, tracks2))) {
436+ ptHadronMix = t2.pt ();
437+ ptElectronMix = t1.ptTrackMc ();
438+ phiElectronMix = t1.phiTrackMc ();
439+ phiHadronMix = t2.phi ();
440+ etaElectronMix = t1.etaTrackMc ();
441+ etaHadronMix = t2.eta ();
442+ if (t2.eta () < etaTrackMin || t2.eta () > etaTrackMax) {
443+ continue ;
444+ }
445+ if (t2.pt () < ptTrackMin) {
446+ continue ;
447+ }
448+ if (ptCondition && (ptElectronMix < ptHadronMix)) {
449+ continue ;
450+ }
451+
452+ deltaPhiMix = RecoDecay::constrainAngle (phiElectronMix - phiHadronMix, -o2::constants::math::PIHalf);
453+ deltaEtaMix = etaElectronMix - etaHadronMix;
454+ bool isNonHfeCorr = false ;
455+ if (t1.isNonHfeMc ()) {
456+ isNonHfeCorr = true ;
457+ registry.fill (HIST (" hMixEventMcGenNonHfEHCorrl" ), ptElectronMix, ptHadronMix, deltaPhiMix, deltaEtaMix);
458+ } else {
459+
460+ registry.fill (HIST (" hMixEventMcGenInclusiveEHCorrl" ), ptElectronMix, ptHadronMix, deltaPhiMix, deltaEtaMix);
461+ }
462+
463+ entryElectronHadronPairmcGen (deltaPhiMix, deltaEtaMix, ptElectronMix, ptHadronMix, poolBin, isNonHfeCorr);
464+ }
464465 }
465-
466-
467- entryElectronHadronPairmcGen (deltaPhiMix, deltaEtaMix, ptElectronMix, ptHadronMix, poolBin, isNonHfeCorr);
468- }
469- }
470466 }
471467 PROCESS_SWITCH (HfCorrelatorHfeHadrons, processMcGenMixedEvent, " Process Mixed Event MC Gen mode" , false );
472468};
0 commit comments