@@ -2182,67 +2182,67 @@ struct HfTaskFlow {
21822182 } // end of trigger loop
21832183 } // end of fillCorrelationsFt0aFt0c
21842184
2185- template <typename TTarget, typename TTracksTrig, typename TTracksAssoc>
2186- void fillCorrelationsMonteCarlo (TTarget target, CorrelationContainer::CFStep step,
2187- TTracksTrig const & tracks1, TTracksAssoc const & tracks2,
2188- float multiplicity, float posZ, bool sameEvent)
2189- {
2190- auto triggerWeight = 1 ;
2191- auto associatedWeight = 1 ;
2192- auto loopCounter = 0 ; // To avoid filling associated tracks QA many times, I fill it only for the first trigger track of the collision
2193- int sampleIndex = gRandom ->Uniform (0 , configTask.nSamples );
2194-
2195- for (auto const & track1 : tracks1) {
2196-
2197- if (track1.eta () < configTask.etaMcParticlesTriggerMin || track1.eta () > configTask.etaMcParticlesTriggerMax ) {
2198- continue ;
2199- }
2200- if (track1.pt () < configTask.ptMcParticlesTriggerMin || track1.pt () > configTask.ptMcParticlesTriggerMax ) {
2201- continue ;
2202- }
2203- if (step >= CorrelationContainer::kCFStepTrackedOnlyPrim && !track1.isPhysicalPrimary ()) {
2204- continue ;
2205- }
2206-
2207- target->getTriggerHist ()->Fill (step, track1.pt (), multiplicity, posZ, triggerWeight);
2208- if (configTask.doEtaDependentFlow ) {
2209- registry.fill (HIST (" Trig_hist" ), sampleIndex, posZ, track1.eta (), triggerWeight);
2210- } else {
2211- registry.fill (HIST (" Trig_hist" ), sampleIndex, posZ, track1.pt (), triggerWeight);
2212- }
2185+ // template <typename TTarget, typename TTracksTrig, typename TTracksAssoc>
2186+ // void fillCorrelationsMonteCarlo(TTarget target, CorrelationContainer::CFStep step,
2187+ // TTracksTrig const& tracks1, TTracksAssoc const& tracks2,
2188+ // float multiplicity, float posZ, bool sameEvent)
2189+ // {
2190+ // auto triggerWeight = 1;
2191+ // auto associatedWeight = 1;
2192+ // auto loopCounter = 0; // To avoid filling associated tracks QA many times, I fill it only for the first trigger track of the collision
2193+ // int sampleIndex = gRandom->Uniform(0, configTask.nSamples);
22132194
2214- for (auto const & track2 : tracks2 ) {
2195+ // for (auto const& track1 : tracks1 ) {
22152196
2216- if (track1.globalIndex () == track2.globalIndex ()) {
2217- continue ;
2218- }
2219- if (track2.eta () < configTask.etaMcParticlesAssocMin || track2.eta () > configTask.etaMcParticlesAssocMax ) {
2220- continue ;
2221- }
2222- if (track2.pt () < configTask.ptMcParticlesAssocMin || track2.pt () > configTask.ptMcParticlesAssocMax ) {
2223- continue ;
2224- }
2225- if (step >= CorrelationContainer::kCFStepTrackedOnlyPrim && !track2.isPhysicalPrimary ()) {
2226- continue ;
2227- }
2228-
2229- float deltaPhi = RecoDecay::constrainAngle (track2.phi () - track1.phi (), -PIHalf);
2230- float deltaEta = track2.eta () - track1.eta ();
2197+ // if (track1.eta() < configTask.etaMcParticlesTriggerMin || track1.eta() > configTask.etaMcParticlesTriggerMax) {
2198+ // continue;
2199+ // }
2200+ // if (track1.pt() < configTask.ptMcParticlesTriggerMin || track1.pt() > configTask.ptMcParticlesTriggerMax) {
2201+ // continue;
2202+ // }
2203+ // if (step >= CorrelationContainer::kCFStepTrackedOnlyPrim && !track1.isPhysicalPrimary()) {
2204+ // continue;
2205+ // }
22312206
2232- if (!configTask.doEtaDependentFlow && !configTask.doVariationContainers ) {
2233- target->getPairHist ()->Fill (step, deltaEta, track2.pt (), track1.pt (), multiplicity, deltaPhi, posZ,
2234- triggerWeight * associatedWeight);
2235- } else if (configTask.doEtaDependentFlow ) {
2236- target->getPairHist ()->Fill (step, sampleIndex, posZ, track2.eta (), track1.eta (), deltaPhi, deltaEta,
2237- triggerWeight * associatedWeight);
2238- } else {
2239- target->getPairHist ()->Fill (step, sampleIndex, posZ, track1.pt (), multiplicity, deltaPhi, deltaEta,
2240- triggerWeight * associatedWeight);
2241- }
2242- } // end of loop over track2
2207+ // target->getTriggerHist()->Fill(step, track1.pt(), multiplicity, posZ, triggerWeight);
2208+ // if (configTask.doEtaDependentFlow) {
2209+ // registry.fill(HIST("Trig_hist"), sampleIndex, posZ, track1.eta(), triggerWeight);
2210+ // } else {
2211+ // registry.fill(HIST("Trig_hist"), sampleIndex, posZ, track1.pt(), triggerWeight);
2212+ // }
22432213
2244- } // end of loop over track1
2245- }
2214+ // for (auto const& track2 : tracks2) {
2215+
2216+ // if (track1.globalIndex() == track2.globalIndex()) {
2217+ // continue;
2218+ // }
2219+ // if (track2.eta() < configTask.etaMcParticlesAssocMin || track2.eta() > configTask.etaMcParticlesAssocMax) {
2220+ // continue;
2221+ // }
2222+ // if (track2.pt() < configTask.ptMcParticlesAssocMin || track2.pt() > configTask.ptMcParticlesAssocMax) {
2223+ // continue;
2224+ // }
2225+ // if (step >= CorrelationContainer::kCFStepTrackedOnlyPrim && !track2.isPhysicalPrimary()) {
2226+ // continue;
2227+ // }
2228+
2229+ // float deltaPhi = RecoDecay::constrainAngle(track2.phi() - track1.phi(), -PIHalf);
2230+ // float deltaEta = track2.eta() - track1.eta();
2231+
2232+ // if (!configTask.doEtaDependentFlow && !configTask.doVariationContainers) {
2233+ // target->getPairHist()->Fill(step, deltaEta, track2.pt(), track1.pt(), multiplicity, deltaPhi, posZ,
2234+ // triggerWeight * associatedWeight);
2235+ // } else if (configTask.doEtaDependentFlow) {
2236+ // target->getPairHist()->Fill(step, sampleIndex, posZ, track2.eta(), track1.eta(), deltaPhi, deltaEta,
2237+ // triggerWeight * associatedWeight);
2238+ // } else {
2239+ // target->getPairHist()->Fill(step, sampleIndex, posZ, track1.pt(), multiplicity, deltaPhi, deltaEta,
2240+ // triggerWeight * associatedWeight);
2241+ // }
2242+ // } // end of loop over track2
2243+
2244+ // } // end of loop over track1
2245+ // }
22462246
22472247 // ===============================================================================================================================================================================
22482248 // mixCollisions for RECONSTRUCTED events
0 commit comments