@@ -212,11 +212,11 @@ struct BJetTaggingML {
212212 registry.fill (HIST (" h2_SVMass_jetpT" ), myJet.pt (), massSV, eventweight);
213213
214214 if (doprocessMCJets) {
215- if (jetFlavor == 2 ) {
215+ if (jetFlavor == JetTaggingSpecies::beauty ) {
216216 registry.fill (HIST (" h2_LxyS_jetpT_bjet" ), myJet.pt (), candSV.decayLengthXY () / candSV.errorDecayLengthXY (), eventweight);
217217 registry.fill (HIST (" h2_Dispersion_jetpT_bjet" ), myJet.pt (), candSV.dispersion (), eventweight);
218218 registry.fill (HIST (" h2_SVMass_jetpT_bjet" ), myJet.pt (), massSV, eventweight);
219- } else if (jetFlavor == 1 ) {
219+ } else if (jetFlavor == JetTaggingSpecies::charm ) {
220220 registry.fill (HIST (" h2_LxyS_jetpT_cjet" ), myJet.pt (), candSV.decayLengthXY () / candSV.errorDecayLengthXY (), eventweight);
221221 registry.fill (HIST (" h2_Dispersion_jetpT_cjet" ), myJet.pt (), candSV.dispersion (), eventweight);
222222 registry.fill (HIST (" h2_SVMass_jetpT_cjet" ), myJet.pt (), massSV, eventweight);
@@ -250,10 +250,10 @@ struct BJetTaggingML {
250250 registry.fill (HIST (" h2_SIPs3D_jetpT" ), analysisJet.pt (), sign * std::abs (constituent.dcaXYZ ()) / constituent.sigmadcaXYZ (), eventweight);
251251
252252 if (doprocessMCJets) {
253- if (jetFlavor == 2 ) {
253+ if (jetFlavor == JetTaggingSpecies::beauty ) {
254254 registry.fill (HIST (" h2_SIPs2D_jetpT_bjet" ), analysisJet.pt (), sign * std::abs (constituent.dcaXY ()) / constituent.sigmadcaXY (), eventweight);
255255 registry.fill (HIST (" h2_SIPs3D_jetpT_bjet" ), analysisJet.pt (), sign * std::abs (constituent.dcaXYZ ()) / constituent.sigmadcaXYZ (), eventweight);
256- } else if (jetFlavor == 1 ) {
256+ } else if (jetFlavor == JetTaggingSpecies::charm ) {
257257 registry.fill (HIST (" h2_SIPs2D_jetpT_cjet" ), analysisJet.pt (), sign * std::abs (constituent.dcaXY ()) / constituent.sigmadcaXY (), eventweight);
258258 registry.fill (HIST (" h2_SIPs3D_jetpT_cjet" ), analysisJet.pt (), sign * std::abs (constituent.dcaXYZ ()) / constituent.sigmadcaXYZ (), eventweight);
259259 } else {
@@ -385,23 +385,23 @@ struct BJetTaggingML {
385385
386386 if (doDataDriven) {
387387 registry.fill (HIST (" hSparse_Incljets" ), analysisJet.pt (), analysisJet.scoreML (), useDb ? 0 : -1 * std::log (1 - analysisJet.scoreML ()), analysisJet.mass (), -1 * std::log (analysisJet.jetProb ()), svsParams[0 ].svMass , svsParams[0 ].svfE , eventWeight);
388- if (jetFlavor == 2 ) {
388+ if (jetFlavor == JetTaggingSpecies::beauty ) {
389389 registry.fill (HIST (" hSparse_bjets" ), analysisJet.pt (), analysisJet.scoreML (), useDb ? 0 : -1 * std::log (1 - analysisJet.scoreML ()), analysisJet.mass (), -1 * std::log (analysisJet.jetProb ()), svsParams[0 ].svMass , svsParams[0 ].svfE , eventWeight);
390- } else if (jetFlavor == 1 ) {
390+ } else if (jetFlavor == JetTaggingSpecies::charm ) {
391391 registry.fill (HIST (" hSparse_cjets" ), analysisJet.pt (), analysisJet.scoreML (), useDb ? 0 : -1 * std::log (1 - analysisJet.scoreML ()), analysisJet.mass (), -1 * std::log (analysisJet.jetProb ()), svsParams[0 ].svMass , svsParams[0 ].svfE , eventWeight);
392392 } else {
393393 registry.fill (HIST (" hSparse_lfjets" ), analysisJet.pt (), analysisJet.scoreML (), useDb ? 0 : -1 * std::log (1 - analysisJet.scoreML ()), analysisJet.mass (), -1 * std::log (analysisJet.jetProb ()), svsParams[0 ].svMass , svsParams[0 ].svfE , eventWeight);
394394 }
395395 }
396396
397- if (jetFlavor == 2 ) {
397+ if (jetFlavor == JetTaggingSpecies::beauty ) {
398398 registry.fill (HIST (" h2_score_jetpT_bjet" ), analysisJet.pt (), analysisJet.scoreML (), eventWeight);
399399 if (!useDb) {
400400 registry.fill (HIST (" h2_logscore_jetpT_bjet" ), analysisJet.pt (), -1 * std::log (1 - analysisJet.scoreML ()), eventWeight);
401401 }
402402 registry.fill (HIST (" h2_jetMass_jetpT_bjet" ), analysisJet.pt (), analysisJet.mass (), eventWeight);
403403 registry.fill (HIST (" h_jetpT_detector_bjet" ), analysisJet.pt (), eventWeight);
404- } else if (jetFlavor == 1 ) {
404+ } else if (jetFlavor == JetTaggingSpecies::charm ) {
405405 registry.fill (HIST (" h2_score_jetpT_cjet" ), analysisJet.pt (), analysisJet.scoreML (), eventWeight);
406406 if (!useDb) {
407407 registry.fill (HIST (" h2_logscore_jetpT_cjet" ), analysisJet.pt (), -1 * std::log (1 - analysisJet.scoreML ()), eventWeight);
@@ -422,9 +422,9 @@ struct BJetTaggingML {
422422 continue ;
423423 }
424424
425- if (jetFlavor == 2 ) {
425+ if (jetFlavor == JetTaggingSpecies::beauty ) {
426426 registry.fill (HIST (" h2_Response_DetjetpT_PartjetpT_bjet" ), analysisJet.pt (), mcpjet.pt (), eventWeight);
427- } else if (jetFlavor == 1 ) {
427+ } else if (jetFlavor == JetTaggingSpecies::charm ) {
428428 registry.fill (HIST (" h2_Response_DetjetpT_PartjetpT_cjet" ), analysisJet.pt (), mcpjet.pt (), eventWeight);
429429 } else {
430430 registry.fill (HIST (" h2_Response_DetjetpT_PartjetpT_lfjet" ), analysisJet.pt (), mcpjet.pt (), eventWeight);
@@ -457,9 +457,9 @@ struct BJetTaggingML {
457457
458458 registry.fill (HIST (" h_jetpT_particle_DetColl" ), mcpjet.pt (), eventWeight);
459459
460- if (jetFlavor == 2 ) {
460+ if (jetFlavor == JetTaggingSpecies::beauty ) {
461461 registry.fill (HIST (" h_jetpT_particle_DetColl_bjet" ), mcpjet.pt (), eventWeight);
462- } else if (jetFlavor == 1 ) {
462+ } else if (jetFlavor == JetTaggingSpecies::charm ) {
463463 registry.fill (HIST (" h_jetpT_particle_DetColl_cjet" ), mcpjet.pt (), eventWeight);
464464 } else {
465465 registry.fill (HIST (" h_jetpT_particle_DetColl_lfjet" ), mcpjet.pt (), eventWeight);
@@ -496,9 +496,9 @@ struct BJetTaggingML {
496496
497497 int8_t jetFlavor = mcpjet.origin ();
498498
499- if (jetFlavor == 2 ) {
499+ if (jetFlavor == JetTaggingSpecies::beauty ) {
500500 registry.fill (HIST (" h_jetpT_particle_bjet" ), mcpjet.pt (), eventWeight);
501- } else if (jetFlavor == 1 ) {
501+ } else if (jetFlavor == JetTaggingSpecies::charm ) {
502502 registry.fill (HIST (" h_jetpT_particle_cjet" ), mcpjet.pt (), eventWeight);
503503 } else {
504504 registry.fill (HIST (" h_jetpT_particle_lfjet" ), mcpjet.pt (), eventWeight);
0 commit comments