@@ -113,6 +113,12 @@ struct JetTaggerHFQA {
113113 ConfigurableAxis binSigmaLxyz{" binSigmaLxyz" , {100 , 0 ., 0.1 }, " " };
114114
115115 int numberOfJetFlavourSpecies = 6 ;
116+ int firstTaggerForTrackCounting = 0 ;
117+ int secondTaggerForTrackCounting = 1 ;
118+ int thirdTaggerForTrackCounting = 2 ;
119+ float kUnsetJetAreaFraction = -98.0 ;
120+ float kUnsetConstituentPtMin = -98.0 ;
121+ float kUnsetConstituentPtMax = 9998.0 ;
116122 std::vector<int > eventSelectionBits;
117123 int trackSelection = -1 ;
118124
@@ -446,14 +452,14 @@ struct JetTaggerHFQA {
446452 template <typename T, typename U>
447453 bool isAcceptedJet (U const & jet)
448454 {
449- if (jetAreaFractionMin > - 98.0 ) {
455+ if (jetAreaFractionMin > kUnsetJetAreaFraction ) {
450456 if (jet.area () < jetAreaFractionMin * o2::constants::math::PI * (jet.r () / 100.0 ) * (jet.r () / 100.0 )) {
451457 return false ;
452458 }
453459 }
454460 bool checkConstituentPt = true ;
455- bool checkConstituentMinPt = (leadingConstituentPtMin > - 98.0 );
456- bool checkConstituentMaxPt = (leadingConstituentPtMax < 9998.0 );
461+ bool checkConstituentMinPt = (leadingConstituentPtMin > kUnsetConstituentPtMin );
462+ bool checkConstituentMaxPt = (leadingConstituentPtMax < kUnsetConstituentPtMax );
457463 if (!checkConstituentMinPt && !checkConstituentMaxPt) {
458464 checkConstituentPt = false ;
459465 }
@@ -596,23 +602,23 @@ struct JetTaggerHFQA {
596602 if (fillIPxyz)
597603 std::sort (vecSignImpXYZSig.begin (), vecSignImpXYZSig.end (), sortImp);
598604
599- if (vecSignImpXYSig.size () > 0 ) { // N1
605+ if (vecSignImpXYSig.size () > firstTaggerForTrackCounting ) { // N1
600606 if (fillIPxy)
601607 registry.fill (HIST (" h2_jet_pt_sign_impact_parameter_xy_significance_N1" ), jet.pt (), vecSignImpXYSig[0 ][0 ]);
602608 if (fillIPz)
603609 registry.fill (HIST (" h2_jet_pt_sign_impact_parameter_z_significance_N1" ), jet.pt (), vecSignImpZSig[0 ][0 ]);
604610 if (fillIPxyz)
605611 registry.fill (HIST (" h2_jet_pt_sign_impact_parameter_xyz_significance_N1" ), jet.pt (), vecSignImpXYZSig[0 ][0 ]);
606612 }
607- if (vecSignImpXYSig.size () > 1 ) { // N2
613+ if (vecSignImpXYSig.size () > secondTaggerForTrackCounting ) { // N2
608614 if (fillIPxy)
609615 registry.fill (HIST (" h2_jet_pt_sign_impact_parameter_xy_significance_N2" ), jet.pt (), vecSignImpXYSig[1 ][0 ]);
610616 if (fillIPz)
611617 registry.fill (HIST (" h2_jet_pt_sign_impact_parameter_z_significance_N2" ), jet.pt (), vecSignImpZSig[1 ][0 ]);
612618 if (fillIPxyz)
613619 registry.fill (HIST (" h2_jet_pt_sign_impact_parameter_xyz_significance_N2" ), jet.pt (), vecSignImpXYZSig[1 ][0 ]);
614620 }
615- if (vecSignImpXYSig.size () > 2 ) { // N3
621+ if (vecSignImpXYSig.size () > thirdTaggerForTrackCounting ) { // N3
616622 if (fillIPxy)
617623 registry.fill (HIST (" h2_jet_pt_sign_impact_parameter_xy_significance_N3" ), jet.pt (), vecSignImpXYSig[2 ][0 ]);
618624 if (fillIPz)
@@ -750,23 +756,23 @@ struct JetTaggerHFQA {
750756 sort (vecSignImpZSig[jetflavour].begin (), vecSignImpZSig[jetflavour].end (), std::greater<float >());
751757 sort (vecSignImpXYZSig[jetflavour].begin (), vecSignImpXYZSig[jetflavour].end (), std::greater<float >());
752758
753- if (vecImpXY[jetflavour].size () > 0 ) { // N1
759+ if (vecImpXY[jetflavour].size () > firstTaggerForTrackCounting ) { // N1
754760 if (fillIPxy)
755761 registry.fill (HIST (" h3_jet_pt_sign_impact_parameter_xy_significance_flavour_N1" ), mcdjet.pt (), vecSignImpXYSig[jetflavour][0 ], jetflavour, eventWeight);
756762 if (fillIPz)
757763 registry.fill (HIST (" h3_jet_pt_sign_impact_parameter_z_significance_flavour_N1" ), mcdjet.pt (), vecSignImpZSig[jetflavour][0 ], jetflavour, eventWeight);
758764 if (fillIPxyz)
759765 registry.fill (HIST (" h3_jet_pt_sign_impact_parameter_xyz_significance_flavour_N1" ), mcdjet.pt (), vecSignImpXYZSig[jetflavour][0 ], jetflavour, eventWeight);
760766 }
761- if (vecImpXY[jetflavour].size () > 1 ) { // N2
767+ if (vecImpXY[jetflavour].size () > secondTaggerForTrackCounting ) { // N2
762768 if (fillIPxy)
763769 registry.fill (HIST (" h3_jet_pt_sign_impact_parameter_xy_significance_flavour_N2" ), mcdjet.pt (), vecSignImpXYSig[jetflavour][1 ], jetflavour, eventWeight);
764770 if (fillIPz)
765771 registry.fill (HIST (" h3_jet_pt_sign_impact_parameter_z_significance_flavour_N2" ), mcdjet.pt (), vecSignImpZSig[jetflavour][1 ], jetflavour, eventWeight);
766772 if (fillIPxyz)
767773 registry.fill (HIST (" h3_jet_pt_sign_impact_parameter_xyz_significance_flavour_N2" ), mcdjet.pt (), vecSignImpXYZSig[jetflavour][1 ], jetflavour, eventWeight);
768774 }
769- if (vecImpXY[jetflavour].size () > 2 ) { // N3
775+ if (vecImpXY[jetflavour].size () > thirdTaggerForTrackCounting ) { // N3
770776 if (fillIPxy)
771777 registry.fill (HIST (" h3_jet_pt_sign_impact_parameter_xy_significance_flavour_N3" ), mcdjet.pt (), vecSignImpXYSig[jetflavour][2 ], jetflavour, eventWeight);
772778 if (fillIPz)
@@ -1142,10 +1148,11 @@ struct JetTaggerHFQA {
11421148 }
11431149 int jetflavour = mcdjet.origin ();
11441150 float secondaryPt = 0 ;
1145- float totalJetPt = 0 ;
1151+ float totalJetPt = 0 ;
11461152 for (auto const & track : mcdjet.template tracks_as <JetTagTracksMCD>()) {
11471153 float varImpXY = track.dcaXY () * jettaggingutilities::cmTomum;
1148- if (!track.has_mcParticle ()) continue ;
1154+ if (!track.has_mcParticle ())
1155+ continue ;
11491156 auto mcParticle = track.mcParticle ();
11501157 totalJetPt += track.pt ();
11511158 if (mcParticle.isPhysicalPrimary ()) {
@@ -1633,5 +1640,5 @@ using JetTaggerhfQACharged = JetTaggerHFQA<JetTaggerQAChargedDataJets, aod::Char
16331640WorkflowSpec defineDataProcessing (ConfigContext const & cfgc)
16341641{
16351642 return WorkflowSpec{
1636- adaptAnalysisTask<JetTaggerhfQACharged>(cfgc, TaskName{" jet-taggerhf-qa-charged" })}; // o2-linter: disable=name/o2-task
1643+ adaptAnalysisTask<JetTaggerhfQACharged>(cfgc, TaskName{" jet-taggerhf-qa-charged" })}; // o2-linter: disable=name/o2-task (wrong hyphenation)
16371644}
0 commit comments