@@ -47,6 +47,7 @@ using namespace o2::framework::expressions;
4747
4848struct JetSpectraCharged {
4949
50+ using McParticleCollision = soa::Join<aod::JetMcCollisions, aod::BkgChargedMcRhos>;
5051 using ChargedMCDMatchedJets = soa::Join<aod::ChargedMCDetectorLevelJets, aod::ChargedMCDetectorLevelJetConstituents, aod::ChargedMCDetectorLevelJetsMatchedToChargedMCParticleLevelJets>;
5152 using ChargedMCPMatchedJets = soa::Join<aod::ChargedMCParticleLevelJets, aod::ChargedMCParticleLevelJetConstituents, aod::ChargedMCParticleLevelJetsMatchedToChargedMCDetectorLevelJets>;
5253 using ChargedMCDMatchedJetsWeighted = soa::Join<aod::ChargedMCDetectorLevelJets, aod::ChargedMCDetectorLevelJetConstituents, aod::ChargedMCDetectorLevelJetsMatchedToChargedMCParticleLevelJets, aod::ChargedMCDetectorLevelJetEventWeights>;
@@ -65,6 +66,7 @@ struct JetSpectraCharged {
6566 Configurable<float > trackPtMax{" trackPtMax" , 100.0 , " maximum pT acceptance for tracks" };
6667 Configurable<std::string> trackSelections{" trackSelections" , " globalTracks" , " set track selections" };
6768 Configurable<float > pTHatMaxMCD{" pTHatMaxMCD" , 999.0 , " maximum fraction of hard scattering for jet acceptance in detector MC" };
69+ Configurable<float > pTHatMaxMCP{" pTHatMaxMCP" , 999.0 , " maximum fraction of hard scattering for jet acceptance in particle MC" };
6870 Configurable<float > pTHatExponent{" pTHatExponent" , 6.0 , " exponent of the event weight for the calculation of pTHat" };
6971 Configurable<double > jetPtMax{" jetPtMax" , 200 ., " set jet pT bin max" };
7072 Configurable<float > jetEtaMin{" jetEtaMin" , -0.7 , " minimum jet pseudorapidity" };
@@ -76,8 +78,9 @@ struct JetSpectraCharged {
7678 Configurable<int > trackOccupancyInTimeRangeMax{" trackOccupancyInTimeRangeMax" , 999999 , " maximum track occupancy of tracks in neighbouring collisions in a given time range; only applied to reconstructed collisions (data and mcd jets), not mc collisions (mcp jets)" };
7779 Configurable<int > trackOccupancyInTimeRangeMin{" trackOccupancyInTimeRangeMin" , -999999 , " minimum track occupancy of tracks in neighbouring collisions in a given time range; only applied to reconstructed collisions (data and mcd jets), not mc collisions (mcp jets)" };
7880 Configurable<bool > checkGeoMatched{" checkGeoMatched" , true , " 0: turn off geometry matching, 1: do geometry matching " };
79- Configurable<bool > checkPtMatched{" checkPtMatched" , true , " 0: turn off pT matching, 1: do pT matching" };
80- Configurable<bool > checkGeoPtMatched{" checkGeoPtMatched" , true , " 0: turn off geometry and pT matching, 1: do geometry and pT matching" };
81+ Configurable<bool > checkPtMatched{" checkPtMatched" , false , " 0: turn off pT matching, 1: do pT matching" };
82+ Configurable<bool > checkGeoPtMatched{" checkGeoPtMatched" , false , " 0: turn off geometry and pT matching, 1: do geometry and pT matching" };
83+ Configurable<bool > doMCPJetSubtraction{" doMCPJetSubtraction" , false , " 0: turn off subtraction for MCP jets, 1: do subtraction for MCP jets" };
8184
8285 int eventSelection = -1 ;
8386 int trackSelection = -1 ;
@@ -124,6 +127,28 @@ struct JetSpectraCharged {
124127 }
125128 }
126129
130+ if (doprocessSpectraMCP || doprocessSpectraMCPWeighted) {
131+ registry.add (" h_mc_zvertex" , " position of collision ;#it{Z} (cm)" , {HistType::kTH1F , {{300 , -15.0 , 15.0 }}});
132+ registry.add (" h_mcp_jet_pt" , " jet pT;#it{p}_{T,jet} (GeV/#it{c}); counts" , {HistType::kTH1F , {jetPtAxis}});
133+ registry.add (" h_mcp_jet_eta" , " jet eta;#eta; counts" , {HistType::kTH1F , {jetEtaAxis}});
134+ registry.add (" h_mcp_jet_phi" , " jet phi;#phi; counts" , {HistType::kTH1F , {phiAxis}});
135+ registry.add (" h2_mcp_jet_pt_jet_area" , " jet #it{p}_{T,jet} vs. Area_{jet}; #it{p}_{T,jet} (GeV/#it{c}); Area_{jet}" , {HistType::kTH2F , {jetPtAxis, {150 , 0 ., 1.5 }}});
136+ registry.add (" h2_mcp_jet_pt_jet_ntracks" , " jet #it{p}_{T,jet} vs. N_{jet tracks}; #it{p}_{T,jet} (GeV/#it{c}); N_{jet, tracks}" , {HistType::kTH2F , {jetPtAxis, {200 , -0.5 , 199.5 }}});
137+ registry.add (" h2_mcp_jet_pt_track_pt" , " jet #it{p}_{T,jet} vs. #it{p}_{T,track}; #it{p}_{T,jet} (GeV/#it{c}); #it{p}_{T,track} (GeV/#it{c})" , {HistType::kTH2F , {jetPtAxisRhoAreaSub, trackPtAxis}});
138+ registry.add (" hsparse_mcp_jet_pt_eta_phi" , " mcp jet_pt_eta_phi" , HistType::kTHnSparseF , {jetPtAxis, jetEtaAxis, phiAxis});
139+ if (doMCPJetSubtraction) {
140+ registry.add (" h_mcp_rho" , " mc collision rho;#rho (GeV/#it{c}); counts" , {HistType::kTH1F , {{500 , 0.0 , 500.0 }}});
141+ registry.add (" h_mcp_jet_corr_pt" , " mcp jet corr pT;#it{p}_{T,jet} (GeV/#it{c}); counts" , {HistType::kTH1F , {jetPtAxisRhoAreaSub}});
142+ registry.add (" h2_mcp_jet_corr_pt_jet_area" , " mcp jet #it{p}_{T,jet} vs. Area_{jet}; #it{p}_{T,jet} (GeV/#it{c}); Area_{jet}" , {HistType::kTH2F , {jetPtAxisRhoAreaSub, {150 , 0 ., 1.5 }}});
143+ registry.add (" h2_mcp_jet_corr_pt_jet_ntracks" , " mcp jet #it{p}_{T,jet} vs. N_{jet tracks}; #it{p}_{T,jet} (GeV/#it{c}); N_{jet, tracks}" , {HistType::kTH2F , {jetPtAxisRhoAreaSub, {200 , -0.5 , 199.5 }}});
144+ registry.add (" hsparse_mcp_jet_corr_pt_eta_phi" , " mcp jet_corr_pt_eta_phi" , HistType::kTHnSparseF , {jetPtAxisRhoAreaSub, jetEtaAxis, phiAxis});
145+ }
146+ if (doprocessSpectraMCPWeighted) {
147+ registry.add (" h_mcp_jet_ptcut_part" , " jet #hat{p};#hat{p} (GeV/#it{c});entries" , {HistType::kTH1F , {{1000 , 0 , 1000 }}});
148+ registry.add (" h_mcp_jet_phat_weighted" , " jet #hat{p};#hat{p} (GeV/#it{c});entries" , {HistType::kTH1F , {{1000 , 0 , 1000 }}});
149+ }
150+ }
151+
127152 if (doprocessEvtWiseConstSubJetsData || doprocessEvtWiseConstSubJetsMCD) {
128153 registry.add (" h2_centrality_jet_pt_eventwiseconstituentsubtracted" , " centrality vs. jet pT;centrality;#it{p}_{T,jet} (GeV/#it{c});entries" , {HistType::kTH2F , {centralityAxis, jetPtAxis}});
129154 registry.add (" jet_observables_eventwiseconstituentsubtracted" , " jet_observables_eventwiseconstituentsubtracted" , HistType::kTHnSparseF , {jetPtAxis, jetEtaAxis, phiAxis});
@@ -160,14 +185,19 @@ struct JetSpectraCharged {
160185 }
161186
162187 if (doprocessJetsMatchedSubtracted) {
163- registry.add (" h2_jet_pt_reco_corr_jet_pt_gen_matchedgeo" , " corr pT reco vs. pT gen;#it{p}_{T,jet}^{reco} (GeV/#it{c});#it{p}_{T,jet}^{gen} (GeV/#it{c})" , {HistType::kTH2F , {jetPtAxis, jetPtAxis}});
164- registry.add (" h2_jet_pt_gen_jet_pt_diff_corr_matchedgeo" , " jet gen pT vs. corr delta pT / jet gen pt;#it{p}_{T,jet}^{gen} (GeV/#it{c}); (#it{p}_{T,jet}^{gen} (GeV/#it{c}) - #it{p}_{T,jet}^{reco} (GeV/#it{c})) / #it{p}_{T,jet}^{gen} (GeV/#it{c})" , {HistType::kTH2F , {jetPtAxis, {1000 , -5.0 , 2.0 }}});
165- registry.add (" h2_jet_pt_reco_jet_pt_diff_corr_matchedgeo" , " jet reco pT vs. corr delta pT / jet reco pt;#it{p}_{T,jet}^{reco} (GeV/#it{c}); (#it{p}_{T,jet}^{reco} (GeV/#it{c}) - #it{p}_{T,jet}^{gen} (GeV/#it{c})) / #it{p}_{T,jet}^{reco} (GeV/#it{c})" , {HistType::kTH2F , {jetPtAxis, {1000 , -5.0 , 2.0 }}});
188+ registry.add (" h_mc_collisions_matched" , " mc collisions status;event status;entries" , {HistType::kTH1F , {{4 , 0.0 , 4.0 }}});
189+ registry.add (" h_mcd_events_matched" , " mcd event status;event status;entries" , {HistType::kTH1F , {{4 , 0.0 , 4.0 }}});
190+ registry.add (" h_mc_rho_matched" , " mc collision rho;#rho (GeV/#it{c}); counts" , {HistType::kTH1F , {{500 , -100.0 , 500.0 }}});
191+ registry.add (" h2_reco_jet_corr_pt_gen_jet_corr_pt_matchedgeo" , " corr pT reco vs. corr cpT gen;#it{p}_{T,jet}^{reco} (GeV/#it{c});#it{p}_{T,jet}^{gen} (GeV/#it{c})" , {HistType::kTH2F , {jetPtAxisRhoAreaSub, jetPtAxisRhoAreaSub}});
192+ registry.add (" h2_gen_jet_corr_pt_jet_pt_diff_corr_matchedgeo" , " jet gen corr pT vs. corr delta pT / jet gen corr pt;#it{p}_{T,jet}^{gen} (GeV/#it{c}); (#it{p}_{T,jet}^{gen} (GeV/#it{c}) - #it{p}_{T,jet}^{reco} (GeV/#it{c})) / #it{p}_{T,jet}^{gen} (GeV/#it{c})" , {HistType::kTH2F , {jetPtAxisRhoAreaSub, {1000 , -5.0 , 5.0 }}});
193+ registry.add (" h2_reco_jet_corr_pt_jet_pt_diff_corr_matchedgeo" , " jet reco corr pT vs. corr delta pT / jet reco corr pt;#it{p}_{T,jet}^{reco} (GeV/#it{c}); (#it{p}_{T,jet}^{reco} (GeV/#it{c}) - #it{p}_{T,jet}^{gen} (GeV/#it{c})) / #it{p}_{T,jet}^{reco} (GeV/#it{c})" , {HistType::kTH2F , {jetPtAxisRhoAreaSub, {1000 , -5.0 , 5.0 }}});
194+ registry.add (" h2_gen_jet_corr_pt_jet_corr_pt_ratio_matchedgeo" , " jet gen corr pT vs. jet reco corr pT / jet gen corr pt;#it{p}_{T,jet}^{gen} (GeV/#it{c}); #it{p}_{T,jet}^{reco} (GeV/#it{c}) / #it{p}_{T,jet}^{gen} (GeV/#it{c})" , {HistType::kTH2F , {jetPtAxisRhoAreaSub, {1000 , -5.0 , 5.0 }}});
166195 }
167196 }
168197
169198 Filter trackCuts = (aod::jtrack::pt >= trackPtMin && aod::jtrack::pt < trackPtMax && aod::jtrack::eta > trackEtaMin && aod::jtrack::eta < trackEtaMax);
170199 Filter eventCuts = (nabs(aod::jcollision::posZ) < vertexZCut && aod::jcollision::centrality >= centralityMin && aod::jcollision::centrality < centralityMax);
200+ Preslice<ChargedMCDMatchedJets> mcdjetsPerJCollision = o2::aod::jet::collisionId;
171201
172202 template <typename TTracks, typename TJets>
173203 bool isAcceptedJet (TJets const & jet)
@@ -229,6 +259,33 @@ struct JetSpectraCharged {
229259 }
230260 }
231261
262+ template <typename TJets>
263+ void fillMCPHistograms (TJets const & jet, float rho = 0.0 , float weight = 1.0 )
264+ {
265+ if (jet.r () == round (selectedJetsRadius * 100 .0f )) {
266+ // fill mcp jet histograms
267+ registry.fill (HIST (" h_mcp_jet_pt" ), jet.pt (), weight);
268+ registry.fill (HIST (" h_mcp_jet_eta" ), jet.eta (), weight);
269+ registry.fill (HIST (" h_mcp_jet_phi" ), jet.phi (), weight);
270+ registry.fill (HIST (" hsparse_mcp_jet_pt_eta_phi" ), jet.pt (), jet.eta (), jet.phi (), weight);
271+ registry.fill (HIST (" h2_mcp_jet_pt_jet_area" ), jet.pt (), jet.area (), weight);
272+ registry.fill (HIST (" h2_mcp_jet_pt_jet_ntracks" ), jet.pt (), jet.tracksIds ().size (), weight);
273+ if (doMCPJetSubtraction) {
274+ double jetcorrpt = jet.pt () - (rho * jet.area ());
275+ registry.fill (HIST (" h_mcp_jet_corr_pt" ), jetcorrpt, weight);
276+ registry.fill (HIST (" hsparse_mcp_jet_corr_pt_eta_phi" ), jet.pt (), jet.eta (), jet.phi (), weight);
277+ if (jetcorrpt > 0 ) {
278+ registry.fill (HIST (" h2_mcp_jet_corr_pt_jet_area" ), jetcorrpt, jet.area (), weight);
279+ registry.fill (HIST (" h2_mcp_jet_corr_pt_jet_ntracks" ), jetcorrpt, jet.tracksIds ().size (), weight);
280+ }
281+ }
282+ }
283+
284+ for (const auto & constituent : jet.template tracks_as <aod::JetParticles>()) {
285+ registry.fill (HIST (" h2_mcp_jet_pt_track_pt" ), jet.pt (), constituent.pt (), weight);
286+ }
287+ }
288+
232289 template <typename TJets>
233290 void fillEventWiseConstituentSubtractedHistograms (TJets const & jet, float centrality, float weight = 1.0 )
234291 {
@@ -315,7 +372,7 @@ struct JetSpectraCharged {
315372 }
316373
317374 template <typename TBase, typename TTag>
318- void fillGeoMatchedCorrHistograms (TBase const & jetBase, float rho, float weight = 1.0 )
375+ void fillGeoMatchedCorrHistograms (TBase const & jetBase, float rho, float mcrho = 0.0 , float weight = 1.0 )
319376 {
320377 float pTHat = 10 . / (std::pow (weight, 1.0 / pTHatExponent));
321378 if (jetBase.pt () > pTHatMaxMCD * pTHat) {
@@ -327,11 +384,13 @@ struct JetSpectraCharged {
327384 continue ;
328385 }
329386 if (jetBase.r () == round (selectedJetsRadius * 100 .0f )) {
387+ double corrTagjetpt = jetTag.pt () - (mcrho * jetTag.area ());
330388 double corrBasejetpt = jetBase.pt () - (rho * jetBase.area ());
331- double dcorrpt = jetTag.pt () - corrBasejetpt;
332- registry.fill (HIST (" h2_jet_pt_reco_corr_jet_pt_gen_matchedgeo" ), corrBasejetpt, jetTag.pt (), weight);
333- registry.fill (HIST (" h2_jet_pt_gen_jet_pt_diff_corr_matchedgeo" ), jetTag.pt (), dcorrpt / jetTag.pt (), weight);
334- registry.fill (HIST (" h2_jet_pt_reco_jet_pt_diff_corr_matchedgeo" ), corrBasejetpt, dcorrpt / corrBasejetpt, weight);
389+ double dcorrpt = corrTagjetpt - corrBasejetpt;
390+ registry.fill (HIST (" h2_reco_jet_corr_pt_gen_jet_corr_pt_matchedgeo" ), corrBasejetpt, corrTagjetpt, weight);
391+ registry.fill (HIST (" h2_gen_jet_corr_pt_jet_pt_diff_corr_matchedgeo" ), corrTagjetpt, dcorrpt / corrTagjetpt, weight);
392+ registry.fill (HIST (" h2_reco_jet_corr_pt_jet_pt_diff_corr_matchedgeo" ), corrBasejetpt, dcorrpt / corrBasejetpt, weight);
393+ registry.fill (HIST (" h2_gen_jet_corr_pt_jet_corr_pt_ratio_matchedgeo" ), corrTagjetpt, corrBasejetpt / corrTagjetpt, weight);
335394 }
336395 }
337396 }
@@ -462,6 +521,59 @@ struct JetSpectraCharged {
462521 }
463522 PROCESS_SWITCH (JetSpectraCharged, processSpectraMCDWeighted, " jet finder QA mcd with weighted events" , false );
464523
524+ void processSpectraMCP (McParticleCollision::iterator const & mccollision,
525+ soa::Join<aod::ChargedMCParticleLevelJets, aod::ChargedMCParticleLevelJetConstituents> const & jets,
526+ aod::JetParticles const &)
527+ {
528+ if (std::abs (mccollision.posZ ()) > vertexZCut) {
529+ return ;
530+ }
531+ registry.fill (HIST (" h_mc_zvertex" ), mccollision.posZ ());
532+ registry.fill (HIST (" h_mcp_rho" ), mccollision.rho ());
533+ for (auto const & jet : jets) {
534+ if (!jetfindingutilities::isInEtaAcceptance (jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
535+ continue ;
536+ }
537+ if (!isAcceptedJet<aod::JetParticles>(jet)) {
538+ continue ;
539+ }
540+ if (doMCPJetSubtraction) {
541+ fillMCPHistograms (jet, mccollision.rho ());
542+ } else {
543+ fillMCPHistograms (jet);
544+ }
545+ }
546+ }
547+ PROCESS_SWITCH (JetSpectraCharged, processSpectraMCP, " jet spectra for MC particle level" , false );
548+
549+ void processSpectraMCPWeighted (aod::JetMcCollisions::iterator const & mccollision,
550+ soa::Join<aod::ChargedMCParticleLevelJets, aod::ChargedMCParticleLevelJetConstituents, aod::ChargedMCParticleLevelJetEventWeights> const & jets,
551+ aod::JetParticles const &)
552+ {
553+ if (!(std::abs (mccollision.posZ ()) < vertexZCut)) {
554+ return ;
555+ }
556+ registry.fill (HIST (" h_mc_vertexZ" ), mccollision.posZ ());
557+ for (auto const & jet : jets) {
558+ if (!jetfindingutilities::isInEtaAcceptance (jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
559+ continue ;
560+ }
561+ if (!isAcceptedJet<aod::JetParticles>(jet)) {
562+ continue ;
563+ }
564+ float jetweight = jet.eventWeight ();
565+ double pTHat = 10 . / (std::pow (jetweight, 1.0 / pTHatExponent));
566+ for (int N = 1 ; N < 21 ; N++) {
567+ if (jet.pt () < N * 0.25 * pTHat && jet.r () == round (selectedJetsRadius * 100 .0f )) {
568+ registry.fill (HIST (" h_jet_ptcut_part" ), jet.pt (), N * 0.25 , jetweight);
569+ }
570+ }
571+ registry.fill (HIST (" h_mcp_jet_phat_weighted" ), pTHat, jetweight);
572+ fillMCPHistograms (jet, 0.0 , jetweight);
573+ }
574+ }
575+ PROCESS_SWITCH (JetSpectraCharged, processSpectraMCPWeighted, " jet spectra for MC particle level weighted" , false );
576+
465577 void processEvtWiseConstSubJetsData (soa::Filtered<aod::JetCollisions>::iterator const & collision,
466578 soa::Join<aod::ChargedEventWiseSubtractedJets, aod::ChargedEventWiseSubtractedJetConstituents> const & jets,
467579 aod::JetTracksSub const &)
@@ -553,33 +665,54 @@ struct JetSpectraCharged {
553665 }
554666 PROCESS_SWITCH (JetSpectraCharged, processJetsMatchedWeighted, " matched mcp and mcd jets with weighted events" , false );
555667
556- void processJetsMatchedSubtracted (soa::Filtered<soa::Join<aod::JetCollisions, aod::BkgChargedRhos>>::iterator const & collision,
668+ void processJetsMatchedSubtracted (McParticleCollision::iterator const & mccollision,
669+ soa::SmallGroups<soa::Join<aod::JetCollisionsMCD, aod::BkgChargedRhos>> const & collisions,
557670 ChargedMCDMatchedJets const & mcdjets,
558671 ChargedMCPMatchedJets const &,
559672 aod::JetTracks const &, aod::JetParticles const &)
560673 {
561- if (!jetderiveddatautilities::selectCollision (collision, eventSelection)) {
674+ registry.fill (HIST (" h_mc_collisions_matched" ), 0.5 );
675+ if (mccollision.size () < 1 ) {
562676 return ;
563677 }
564- if (collision.trackOccupancyInTimeRange () < trackOccupancyInTimeRangeMin || trackOccupancyInTimeRangeMax < collision.trackOccupancyInTimeRange ()) {
678+ registry.fill (HIST (" h_mc_collisions_matched" ), 1.5 );
679+ if (!(std::abs (mccollision.posZ ()) < vertexZCut)) {
565680 return ;
566681 }
567-
568- for (const auto & mcdjet : mcdjets) {
569- if (!jetfindingutilities::isInEtaAcceptance (mcdjet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
682+ registry.fill (HIST (" h_mc_collisions_matched" ), 2.5 );
683+ double mcrho = mccollision.rho ();
684+ registry.fill (HIST (" h_mc_rho_matched" ), mcrho);
685+ for (const auto & collision : collisions) {
686+ registry.fill (HIST (" h_mcd_events_matched" ), 0.5 );
687+ if (!jetderiveddatautilities::selectCollision (collision, eventSelection) || !(std::abs (collision.posZ ()) < vertexZCut)) {
570688 continue ;
571689 }
572- if (!isAcceptedJet<aod::JetTracks>(mcdjet)) {
573- continue ;
690+ registry.fill (HIST (" h_mcd_events_matched" ), 1.5 );
691+ if (collision.trackOccupancyInTimeRange () < trackOccupancyInTimeRangeMin || trackOccupancyInTimeRangeMax < collision.trackOccupancyInTimeRange ()) {
692+ return ;
693+ }
694+ registry.fill (HIST (" h_mcd_events_matched" ), 2.5 );
695+
696+ auto collmcdjets = mcdjets.sliceBy (mcdjetsPerJCollision, collision.globalIndex ());
697+ for (const auto & mcdjet : collmcdjets) {
698+ if (!jetfindingutilities::isInEtaAcceptance (mcdjet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
699+ continue ;
700+ }
701+ if (!isAcceptedJet<aod::JetTracks>(mcdjet)) {
702+ continue ;
703+ }
704+ if (doMCPJetSubtraction) {
705+ fillGeoMatchedCorrHistograms<ChargedMCDMatchedJets::iterator, ChargedMCPMatchedJets>(mcdjet, collision.rho (), mcrho);
706+ } else {
707+ fillGeoMatchedCorrHistograms<ChargedMCDMatchedJets::iterator, ChargedMCPMatchedJets>(mcdjet, collision.rho ());
708+ }
574709 }
575- // now only do subtraction for MCD jets, need to add subtraction for MCP jets
576- fillGeoMatchedCorrHistograms<ChargedMCDMatchedJets::iterator, ChargedMCPMatchedJets>(mcdjet, collision.rho ());
577710 }
578711 }
579712 PROCESS_SWITCH (JetSpectraCharged, processJetsMatchedSubtracted, " matched mcp and mcd jets after subtraction" , false );
580713};
581714
582715WorkflowSpec defineDataProcessing (ConfigContext const & cfgc)
583716{
584- return WorkflowSpec{adaptAnalysisTask<JetSpectraCharged>(cfgc, TaskName{ " jet-spectra-charged " } )};
717+ return WorkflowSpec{adaptAnalysisTask<JetSpectraCharged>(cfgc)};
585718}
0 commit comments