User will overwrite this function. Histogram booking is no longer done in C++. This function is called in execute once the filters are all passed.
63 if(
m_trigDecTool->ExperimentalAndExpertMethods().isHLTTruncated()){
65 return StatusCode::SUCCESS;
71 if (! eventInfo.isValid() ){
77 if (! hlt_electron_cont.isValid() ) {
82 if (! hlt_muon_cont.isValid() ) {
88 if (! hlt_topoclusters_cont.isValid() ) {
94 if (! hlt_tracks_cont.isValid() ) {
100 if (! hlt_vertex_cont.isValid() ) {
105 if (! offline_vertex_cont.isValid() ) {
111 if (! offline_met_cont.isValid() ) {
117 if (! l1_roi_cont.isValid() ) {
123 if (! l1_jFexMet_cont.isValid() ) {
127 if (! l1_jFexSumEt_cont.isValid() ) {
131 if (! l1_gFexJwojScalar_cont.isValid() ) {
135 if (! l1_gFexJwojMETComponents_cont.isValid() ) {
139 if (! l1_gFexJwojMHTComponents_cont.isValid() ) {
143 if (! l1_gFexJwojMSTComponents_cont.isValid() ) {
147 if (! l1_gFexNCMETScalar_cont.isValid() ) {
151 if (! l1_gFexNCMETComponents_cont.isValid() ) {
155 if (! l1_gFexRhoMETScalar_cont.isValid() ) {
159 if (! l1_gFexRhoMETComponents_cont.isValid() ) {
165 if (! hlt_cell_met_cont.isValid() ) {
170 if (! hlt_mht_met_cont.isValid() ) {
175 if (! hlt_tc_met_cont.isValid() ) {
180 if (! hlt_tc_em_met_cont.isValid() ) {
185 if (! hlt_tcpufit_met_cont.isValid() ) {
190 if (! hlt_tcpufit_sig30_met_cont.isValid() ) {
195 if (! hlt_trkmht_met_cont.isValid() ) {
200 if (! hlt_pfsum_met_cont.isValid() ) {
205 if (! hlt_pfsum_cssk_met_cont.isValid() ) {
210 if (! hlt_pfsum_vssk_met_cont.isValid() ) {
215 if (! hlt_pfopufit_met_cont.isValid() ) {
220 if (! hlt_pfopufit_sig30_met_cont.isValid() ) {
225 if (! hlt_cvfpufit_met_cont.isValid() ) {
230 if (! hlt_mhtpufit_pf_met_cont.isValid() ) {
235 if (! hlt_mhtpufit_em_met_cont.isValid() ) {
240 if (! hlt_met_nn_cont.isValid() ) {
285 double epsilon = 1.189;
293 act_IPBC = eventInfo->actualInteractionsPerCrossing();
297 bool passedPrimaryEl =
false;
300 passedPrimaryEl =
true;
304 if ( hlt_electron_cont.isValid() && passedPrimaryEl ){
305 hlt_el_mult = hlt_electron_cont->size();
307 if( hlt_electron_cont->size() > 0 ) {
308 for (
auto Electron: *hlt_electron_cont) {
316 bool passedPrimaryMu =
false;
319 passedPrimaryMu =
true;
323 if( hlt_muon_cont.isValid() &&passedPrimaryMu ){
324 hlt_mu_mult = hlt_muon_cont->size();
326 if ( hlt_muon_cont->size() > 0 ){
327 for(
auto Muon : *hlt_muon_cont){
335 if(hlt_topoclusters_cont.isValid() && hlt_topoclusters_cont->size() > 0){
336 hlt_topoclusters_mult = hlt_topoclusters_cont->size();
337 for(
auto topoclusters : *hlt_topoclusters_cont){
339 if(hlt_topoclusters_pt > 0){
344 if(hlt_topoclusters_mult > 0){
350 if( hlt_tracks_cont.isValid() && hlt_tracks_cont->size() > 0){
351 hlt_tracks_mult = hlt_tracks_cont->size();
352 float scalarSumPt = 0.0;
353 float scalarSumPx = 0.0;
354 float scalarSumPy = 0.0;
355 for(
auto tracks : *hlt_tracks_cont){
357 hlt_tracks_pt = i_track_pt;
358 scalarSumPt += i_track_pt;
359 scalarSumPx += (tracks->p4().Px());
360 scalarSumPy += (tracks->p4().Py());
362 if(hlt_tracks_pt > hlt_tracks_leading_pt){
363 hlt_tracks_leading_pt = i_track_pt;
366 if(hlt_tracks_pt > 0){
370 if(hlt_tracks_pt > 3){
371 hlt_tracks_eta = (tracks->eta());
372 hlt_tracks_phi = (tracks->phi());
373 fill(
tool, hlt_tracks_eta, hlt_tracks_phi);
377 hlt_tracks_vec_sumPt = std::sqrt(scalarSumPx*scalarSumPx + scalarSumPy*scalarSumPy)/
Gaudi::Units::GeV;
378 hlt_tracks_sca_sumPt = scalarSumPt;
380 fill(
tool, hlt_tracks_mult, hlt_tracks_leading_pt, hlt_tracks_vec_sumPt, hlt_tracks_sca_sumPt);
384 if(hlt_vertex_cont.isValid() && hlt_vertex_cont->size() > 0){
385 hlt_vertex_mult = hlt_vertex_cont->size();
386 hlt_vertex_mult_mu = hlt_vertex_cont->size();
388 for(
auto vertex : *hlt_vertex_cont){
396 hlt_vertex_z = hlt_vertex->
z();
399 if(offline_vertex_cont.isValid() && offline_vertex_cont->size() > 0){
401 for(
auto vertex : *offline_vertex_cont){
409 hlt_vertex_z_diff = hlt_vertex_z - offline_vertex->
z();
415 hlt_vertex_mult = -1;
416 hlt_vertex_mult_mu = 1;
420 fill(
tool, act_IPBC, hlt_vertex_mult_mu);
425 if ( offline_met_cont.isValid() && offline_met_cont->size() > 0 ) {
426 finalTrkMET = ((*offline_met_cont)[
"FinalTrk"]);
427 muonsMET = ((*offline_met_cont)[
"Muons"]);
433 offline_Et = std::sqrt(offline_Ex*offline_Ex + offline_Ey*offline_Ey);
434 offline_Et_eff = std::sqrt(offline_Ex*offline_Ex + offline_Ey*offline_Ey);
435 fill(
tool,offline_Ex,offline_Ey,offline_Et,offline_sumEt);
442 offline_NoMu_Et = std::sqrt(offline_NoMu_Ex*offline_NoMu_Ex + offline_NoMu_Ey*offline_NoMu_Ey);
443 offline_NoMu_Et_eff = std::sqrt(offline_NoMu_Ex*offline_NoMu_Ex + offline_NoMu_Ey*offline_NoMu_Ey);
444 fill(
tool,offline_NoMu_Ex,offline_NoMu_Ey,offline_NoMu_Et,offline_NoMu_sumEt);
452 if (
alg ==
"roi" && l1_roi_cont.isValid()) {
453 l1_met_cont = l1_roi_cont;
457 if ((l1_met_cont->
energyX())>-9e12 && (l1_met_cont->
energyX())<9e12 && (l1_met_cont->
energyY())>-9e12 && (l1_met_cont->
energyY())<9e12) {
460 float L1_met_Et = std::sqrt(L1_met_Ex*L1_met_Ex + L1_met_Ey*L1_met_Ey);
462 float L1_met_Ex_log =
signed_log(L1_met_Ex, epsilon);
463 float L1_met_Ey_log =
signed_log(L1_met_Ey, epsilon);
464 float L1_met_Et_log =
signed_log(L1_met_Et, epsilon);
465 float L1_met_sumEt_log =
signed_log(L1_met_sumEt, epsilon);
466 TVector3
v(L1_met_Ex, L1_met_Ey, 0.0);
467 float L1_met_phi =
v.Phi();
478 fill(
tool, L1_Ex, L1_Ey, L1_Et, L1_sumEt,
479 L1_Ex_log, L1_Ey_log, L1_Et_log, L1_sumEt_log, L1_phi);
485 if (l1_jFexMet_cont.isValid() && l1_jFexMet_cont->size() > 0) {
488 for (
const auto l1_jmet: *l1_jFexMet_cont) {
492 float L1_met_Et = std::sqrt(L1_met_Ex*L1_met_Ex + L1_met_Ey*L1_met_Ey);
493 float L1_met_Ex_log =
signed_log(L1_met_Ex, epsilon);
494 float L1_met_Ey_log =
signed_log(L1_met_Ey, epsilon);
495 float L1_met_Et_log =
signed_log(L1_met_Et, epsilon);
496 TVector3
v(L1_met_Ex, L1_met_Ey, 0.0);
497 float L1_met_phi =
v.Phi();
505 fill(
tool, L1_Ex, L1_Ey, L1_Et, L1_Ex_log, L1_Ey_log, L1_Et_log, L1_phi);
507 if (l1_jFexSumEt_cont.isValid() && l1_jFexSumEt_cont->size() > 0) {
508 float L1_met_sumEt = 0;
509 for (
const auto l1_jsumEt: *l1_jFexSumEt_cont) {
512 float L1_met_sumEt_log =
signed_log(L1_met_sumEt, epsilon);
523 if (l1_gFexJwojScalar_cont.isValid() && l1_gFexJwojScalar_cont->size() > 0) {
524 l1_gmet = l1_gFexJwojScalar_cont->at(0);
527 float L1_met_Et_log =
signed_log(L1_met_Et, epsilon);
528 float L1_met_sumEt_log =
signed_log(L1_met_sumEt, epsilon);
533 fill(
tool, L1_Et, L1_Et_log, L1_sumEt, L1_sumEt_log);
536 if (l1_gFexJwojMETComponents_cont.isValid() && l1_gFexJwojMETComponents_cont->size() > 0) {
537 l1_gmet = l1_gFexJwojMETComponents_cont->at(0);
539 float L1_met_Ex_log =
signed_log(L1_met_Ex, epsilon);
541 float L1_met_Ey_log =
signed_log(L1_met_Ey, epsilon);
542 TVector3
v(L1_met_Ex, L1_met_Ey, 0.0);
543 float L1_met_phi =
v.Phi();
549 fill(
tool, L1_Ex, L1_Ey, L1_Ex_log, L1_Ey_log, L1_phi);
552 if (l1_gFexJwojMHTComponents_cont.isValid() && l1_gFexJwojMHTComponents_cont->size() > 0) {
553 l1_gmet = l1_gFexJwojMHTComponents_cont->at(0);
555 float L1_met_HT_Ex_log =
signed_log(L1_met_HT_Ex, epsilon);
557 float L1_met_HT_Ey_log =
signed_log(L1_met_HT_Ey, epsilon);
558 TVector3
v(L1_met_HT_Ex, L1_met_HT_Ey, 0.0);
559 float L1_met_HT_phi =
v.Phi();
565 fill(
tool, L1_HT_Ex, L1_HT_Ey, L1_HT_Ex_log, L1_HT_Ey_log, L1_HT_phi);
568 if (l1_gFexJwojMSTComponents_cont.isValid() && l1_gFexJwojMSTComponents_cont->size() > 0) {
569 l1_gmet = l1_gFexJwojMSTComponents_cont->at(0);
571 float L1_met_ST_Ex_log =
signed_log(L1_met_ST_Ex, epsilon);
573 float L1_met_ST_Ey_log =
signed_log(L1_met_ST_Ey, epsilon);
574 TVector3
v(L1_met_ST_Ex, L1_met_ST_Ey, 0.0);
575 float L1_met_ST_phi =
v.Phi();
581 fill(
tool, L1_ST_Ex, L1_ST_Ey, L1_ST_Ex_log, L1_ST_Ey_log, L1_ST_phi);
584 if (l1_gFexNCMETScalar_cont.isValid() && l1_gFexNCMETScalar_cont->size() > 0) {
585 l1_gmet = l1_gFexNCMETScalar_cont->at(0);
588 float L1_met_Et_log =
signed_log(L1_met_Et, epsilon);
589 float L1_met_sumEt_log =
signed_log(L1_met_sumEt, epsilon);
594 fill(
tool, L1_Et, L1_Et_log, L1_sumEt, L1_sumEt_log);
597 if (l1_gFexNCMETComponents_cont.isValid() && l1_gFexNCMETComponents_cont->size() > 0) {
598 l1_gmet = l1_gFexNCMETComponents_cont->at(0);
600 float L1_met_Ex_log =
signed_log(L1_met_Ex, epsilon);
602 float L1_met_Ey_log =
signed_log(L1_met_Ey, epsilon);
603 TVector3
v(L1_met_Ex, L1_met_Ey, 0.0);
604 float L1_met_phi =
v.Phi();
610 fill(
tool, L1_Ex, L1_Ey, L1_Ex_log, L1_Ey_log, L1_phi);
613 if (l1_gFexRhoMETScalar_cont.isValid() && l1_gFexRhoMETScalar_cont->size() > 0) {
614 l1_gmet = l1_gFexRhoMETScalar_cont->at(0);
617 float L1_met_Et_log =
signed_log(L1_met_Et, epsilon);
618 float L1_met_sumEt_log =
signed_log(L1_met_sumEt, epsilon);
623 fill(
tool, L1_Et, L1_Et_log, L1_sumEt, L1_sumEt_log);
626 if (l1_gFexRhoMETComponents_cont.isValid() && l1_gFexRhoMETComponents_cont->size() > 0) {
627 l1_gmet = l1_gFexRhoMETComponents_cont->at(0);
629 float L1_met_Ex_log =
signed_log(L1_met_Ex, epsilon);
631 float L1_met_Ey_log =
signed_log(L1_met_Ey, epsilon);
632 TVector3
v(L1_met_Ex, L1_met_Ey, 0.0);
633 float L1_met_phi =
v.Phi();
639 fill(
tool, L1_Ex, L1_Ey, L1_Ex_log, L1_Ey_log, L1_phi);
648 if ( hlt_cell_met_cont.isValid() && hlt_cell_met_cont->size() > 0 ) {
649 hlt_met = hlt_cell_met_cont->at(0);
650 for (
int j=0; j<nStatus; ++j) {
651 unsigned mask = (1
u<<j);
661 for (
int i=0;
i<nComponent; ++
i) {
664 component_Et = sqrt(ex*ex+ey*ey);
669 for (
int i=0;
i<nComponent; ++
i) {
670 for (
int j=0; j<nStatus; ++j) {
671 unsigned mask = (1
u<<j);
673 component_status_weight = 1.;
675 component_status_weight = 0;
679 fill(
tool,mon_comp,mon_bit,component_status_weight);
685 float L1_roiMet_Et = 0;
686 if ( l1_roi_cont.isValid() ) {
687 if ((l1_roi_cont->energyX())>-9e12 && (l1_roi_cont->energyX())<9e12 && (l1_roi_cont->energyY())>-9e12 && (l1_roi_cont->energyY())<9e12) {
690 L1_roiMet_Et = std::sqrt(Ex*Ex + Ey*Ey);
693 float L1_jFexMet_Et = 0;
694 if (l1_jFexMet_cont.isValid() && l1_jFexMet_cont->size() > 0) {
697 for (
const auto l1_jmet: *l1_jFexMet_cont) {
701 L1_jFexMet_Et = std::sqrt(L1_met_Ex*L1_met_Ex + L1_met_Ey*L1_met_Ey);
703 float L1_gFexJWOJMet_Et = 0;
705 float L1_PreSelMet_Et = L1_roiMet_Et;
706 if (
m_L1MetAlg == 1) L1_PreSelMet_Et = L1_jFexMet_Et;
707 if (
m_L1MetAlg == 2) L1_PreSelMet_Et = L1_gFexJWOJMet_Et;
711 if (
alg ==
"cell" && hlt_cell_met_cont.isValid() && hlt_cell_met_cont->size() > 0) {
712 hlt_met = hlt_cell_met_cont->at(0);
713 }
else if (
alg ==
"mht" && hlt_mht_met_cont.isValid() && hlt_mht_met_cont->size() > 0) {
714 hlt_met = hlt_mht_met_cont->at(0);
715 }
else if (
alg ==
"tc" && hlt_tc_met_cont.isValid() && hlt_tc_met_cont->size() > 0) {
716 hlt_met = hlt_tc_met_cont->at(0);
717 }
else if (
alg ==
"tc_em" && hlt_tc_em_met_cont.isValid() && hlt_tc_em_met_cont->size() > 0) {
718 hlt_met = hlt_tc_em_met_cont->at(0);
719 }
else if (
alg ==
"tcpufit" && hlt_tcpufit_met_cont.isValid() && hlt_tcpufit_met_cont->size() > 0) {
720 hlt_met = hlt_tcpufit_met_cont->at(0);
721 }
else if (
alg ==
"tcpufit_sig30" && hlt_tcpufit_sig30_met_cont.isValid() && hlt_tcpufit_sig30_met_cont->size() > 0) {
722 hlt_met = hlt_tcpufit_sig30_met_cont->at(0);
723 }
else if (
alg ==
"trkmht" && hlt_trkmht_met_cont.isValid() && hlt_trkmht_met_cont->size() > 0) {
724 hlt_met = hlt_trkmht_met_cont->at(0);
725 }
else if (
alg ==
"pfsum" && hlt_pfsum_met_cont.isValid() && hlt_pfsum_met_cont->size() > 0) {
726 hlt_met = hlt_pfsum_met_cont->at(0);
727 }
else if (
alg ==
"pfsum_cssk" && hlt_pfsum_cssk_met_cont.isValid() && hlt_pfsum_cssk_met_cont->size() > 0) {
728 hlt_met = hlt_pfsum_cssk_met_cont->at(0);
729 }
else if (
alg ==
"pfsum_vssk" && hlt_pfsum_vssk_met_cont.isValid() && hlt_pfsum_vssk_met_cont->size() > 0) {
730 hlt_met = hlt_pfsum_vssk_met_cont->at(0);
731 }
else if (
alg ==
"pfopufit" && hlt_pfopufit_met_cont.isValid() && hlt_pfopufit_met_cont->size() > 0) {
732 hlt_met = hlt_pfopufit_met_cont->at(0);
733 }
else if (
alg ==
"pfopufit_sig30" && hlt_pfopufit_sig30_met_cont.isValid() && hlt_pfopufit_sig30_met_cont->size() > 0) {
734 hlt_met = hlt_pfopufit_sig30_met_cont->at(0);
735 }
else if (
alg ==
"cvfpufit" && hlt_cvfpufit_met_cont.isValid() && hlt_cvfpufit_met_cont->size() > 0) {
736 hlt_met = hlt_cvfpufit_met_cont->at(0);
737 }
else if (
alg ==
"mhtpufit_pf" && hlt_mhtpufit_pf_met_cont.isValid() && hlt_mhtpufit_pf_met_cont->size() > 0) {
738 hlt_met = hlt_mhtpufit_pf_met_cont->at(0);
739 }
else if (
alg ==
"mhtpufit_em" && hlt_mhtpufit_em_met_cont.isValid() && hlt_mhtpufit_em_met_cont->size() > 0) {
740 hlt_met = hlt_mhtpufit_em_met_cont->at(0);
741 }
else if (
alg ==
"met_nn" && hlt_met_nn_cont.isValid() && hlt_met_nn_cont->size() > 0) {
742 hlt_met = hlt_met_nn_cont->at(0);
751 float hlt_Et = std::sqrt(hlt_Ex*hlt_Ex + hlt_Ey*hlt_Ey);
753 float hlt_Ex_log =
signed_log(hlt_Ex, epsilon);
754 float hlt_Ey_log =
signed_log(hlt_Ey, epsilon);
755 float hlt_Et_log =
signed_log(hlt_Et, epsilon);
756 float hlt_sumEt_log =
signed_log(hlt_sumEt, epsilon);
757 TVector3
v(hlt_Ex, hlt_Ey, hlt_Ez);
758 float hlt_eta =
v.Eta();
759 float hlt_phi =
v.Phi();
771 fill(
tool,met_Ex,met_Ey,met_Et,met_sumEt,
772 met_Ex_log,met_Ey_log,met_Et_log,met_sumEt_log,
780 if (
alg ==
"cell" && hlt_cell_met_cont.isValid() && hlt_cell_met_cont->size() > 0) {
781 hlt_met = hlt_cell_met_cont->at(0);
782 }
else if (
alg ==
"mht" && hlt_mht_met_cont.isValid() && hlt_mht_met_cont->size() > 0) {
783 hlt_met = hlt_mht_met_cont->at(0);
784 }
else if (
alg ==
"tc" && hlt_tc_met_cont.isValid() && hlt_tc_met_cont->size() > 0) {
785 hlt_met = hlt_tc_met_cont->at(0);
786 }
else if (
alg ==
"tc_em" && hlt_tc_em_met_cont.isValid() && hlt_tc_em_met_cont->size() > 0) {
787 hlt_met = hlt_tc_em_met_cont->at(0);
788 }
else if (
alg ==
"tcpufit" && hlt_tcpufit_met_cont.isValid() && hlt_tcpufit_met_cont->size() > 0) {
789 hlt_met = hlt_tcpufit_met_cont->at(0);
790 }
else if (
alg ==
"tcpufit_sig30" && hlt_tcpufit_sig30_met_cont.isValid() && hlt_tcpufit_sig30_met_cont->size() > 0) {
791 hlt_met = hlt_tcpufit_sig30_met_cont->at(0);
792 }
else if (
alg ==
"trkmht" && hlt_trkmht_met_cont.isValid() && hlt_trkmht_met_cont->size() > 0) {
793 hlt_met = hlt_trkmht_met_cont->at(0);
794 }
else if (
alg ==
"pfsum" && hlt_pfsum_met_cont.isValid() && hlt_pfsum_met_cont->size() > 0) {
795 hlt_met = hlt_pfsum_met_cont->at(0);
796 }
else if (
alg ==
"pfsum_cssk" && hlt_pfsum_cssk_met_cont.isValid() && hlt_pfsum_cssk_met_cont->size() > 0) {
797 hlt_met = hlt_pfsum_cssk_met_cont->at(0);
798 }
else if (
alg ==
"pfsum_vssk" && hlt_pfsum_vssk_met_cont.isValid() && hlt_pfsum_vssk_met_cont->size() > 0) {
799 hlt_met = hlt_pfsum_vssk_met_cont->at(0);
800 }
else if (
alg ==
"pfopufit" && hlt_pfopufit_met_cont.isValid() && hlt_pfopufit_met_cont->size() > 0) {
801 hlt_met = hlt_pfopufit_met_cont->at(0);
802 }
else if (
alg ==
"pfopufit_sig30" && hlt_pfopufit_sig30_met_cont.isValid() && hlt_pfopufit_sig30_met_cont->size() > 0) {
803 hlt_met = hlt_pfopufit_sig30_met_cont->at(0);
804 }
else if (
alg ==
"cvfpufit" && hlt_cvfpufit_met_cont.isValid() && hlt_cvfpufit_met_cont->size() > 0) {
805 hlt_met = hlt_cvfpufit_met_cont->at(0);
806 }
else if (
alg ==
"mhtpufit_pf" && hlt_mhtpufit_pf_met_cont.isValid() && hlt_mhtpufit_pf_met_cont->size() > 0) {
807 hlt_met = hlt_mhtpufit_pf_met_cont->at(0);
808 }
else if (
alg ==
"mhtpufit_em" && hlt_mhtpufit_em_met_cont.isValid() && hlt_mhtpufit_em_met_cont->size() > 0) {
809 hlt_met = hlt_mhtpufit_em_met_cont->at(0);
810 }
else if (
alg ==
"met_nn" && hlt_met_nn_cont.isValid() && hlt_met_nn_cont->size() > 0) {
811 hlt_met = hlt_met_nn_cont->at(0);
819 float hlt_Et = std::sqrt(hlt_Ex*hlt_Ex + hlt_Ey*hlt_Ey);
820 if (L1_PreSelMet_Et >
m_L1MetCut && !std::isnan(hlt_Et)) {
828 if (
alg ==
"cell" && hlt_cell_met_cont.isValid() && hlt_cell_met_cont->size() > 0) {
829 hlt_met = hlt_cell_met_cont->at(0);
830 }
else if (
alg ==
"mht" && hlt_mht_met_cont.isValid() && hlt_mht_met_cont->size() > 0) {
831 hlt_met = hlt_mht_met_cont->at(0);
832 }
else if (
alg ==
"tc" && hlt_tc_met_cont.isValid() && hlt_tc_met_cont->size() > 0) {
833 hlt_met = hlt_tc_met_cont->at(0);
834 }
else if (
alg ==
"tc_em" && hlt_tc_em_met_cont.isValid() && hlt_tc_em_met_cont->size() > 0) {
835 hlt_met = hlt_tc_em_met_cont->at(0);
836 }
else if (
alg ==
"tcpufit" && hlt_tcpufit_met_cont.isValid() && hlt_tcpufit_met_cont->size() > 0) {
837 hlt_met = hlt_tcpufit_met_cont->at(0);
838 }
else if (
alg ==
"trkmht" && hlt_trkmht_met_cont.isValid() && hlt_trkmht_met_cont->size() > 0) {
839 hlt_met = hlt_trkmht_met_cont->at(0);
840 }
else if (
alg ==
"pfsum" && hlt_pfsum_met_cont.isValid() && hlt_pfsum_met_cont->size() > 0) {
841 hlt_met = hlt_pfsum_met_cont->at(0);
842 }
else if (
alg ==
"pfsum_cssk" && hlt_pfsum_cssk_met_cont.isValid() && hlt_pfsum_cssk_met_cont->size() > 0) {
843 hlt_met = hlt_pfsum_cssk_met_cont->at(0);
844 }
else if (
alg ==
"pfsum_vssk" && hlt_pfsum_vssk_met_cont.isValid() && hlt_pfsum_vssk_met_cont->size() > 0) {
845 hlt_met = hlt_pfsum_vssk_met_cont->at(0);
846 }
else if (
alg ==
"pfopufit" && hlt_pfopufit_met_cont.isValid() && hlt_pfopufit_met_cont->size() > 0) {
847 hlt_met = hlt_pfopufit_met_cont->at(0);
848 }
else if (
alg ==
"cvfpufit" && hlt_cvfpufit_met_cont.isValid() && hlt_cvfpufit_met_cont->size() > 0) {
849 hlt_met = hlt_cvfpufit_met_cont->at(0);
850 }
else if (
alg ==
"mhtpufit_pf" && hlt_mhtpufit_pf_met_cont.isValid() && hlt_mhtpufit_pf_met_cont->size() > 0) {
851 hlt_met = hlt_mhtpufit_pf_met_cont->at(0);
852 }
else if (
alg ==
"mhtpufit_em" && hlt_mhtpufit_em_met_cont.isValid() && hlt_mhtpufit_em_met_cont->size() > 0) {
853 hlt_met = hlt_mhtpufit_em_met_cont->at(0);
862 float hlt_Et = std::sqrt(hlt_Ex*hlt_Ex + hlt_Ey*hlt_Ey);
864 float hlt_Ex_log =
signed_log(hlt_Ex, epsilon);
865 float hlt_Ey_log =
signed_log(hlt_Ey, epsilon);
866 float hlt_Et_log =
signed_log(hlt_Et, epsilon);
867 float hlt_sumEt_log =
signed_log(hlt_sumEt, epsilon);
868 TVector3
v(hlt_Ex, hlt_Ey, hlt_Ez);
869 float hlt_eta =
v.Eta();
870 float hlt_phi =
v.Phi();
886 fill(
tool,met_Ex,met_Ey,met_Et,met_sumEt,
887 met_Ex_log,met_Ey_log,met_Et_log,met_sumEt_log,
895 if (
alg ==
"cell" && hlt_cell_met_cont.isValid() && hlt_cell_met_cont->size() > 0) {
896 hlt_met = hlt_cell_met_cont->at(0);
897 }
else if (
alg ==
"mht" && hlt_mht_met_cont.isValid() && hlt_mht_met_cont->size() > 0) {
898 hlt_met = hlt_mht_met_cont->at(0);
899 }
else if (
alg ==
"tc" && hlt_tc_met_cont.isValid() && hlt_tc_met_cont->size() > 0) {
900 hlt_met = hlt_tc_met_cont->at(0);
901 }
else if (
alg ==
"tc_em" && hlt_tc_em_met_cont.isValid() && hlt_tc_em_met_cont->size() > 0) {
902 hlt_met = hlt_tc_em_met_cont->at(0);
903 }
else if (
alg ==
"tcpufit" && hlt_tcpufit_met_cont.isValid() && hlt_tcpufit_met_cont->size() > 0) {
904 hlt_met = hlt_tcpufit_met_cont->at(0);
905 }
else if (
alg ==
"trkmht" && hlt_trkmht_met_cont.isValid() && hlt_trkmht_met_cont->size() > 0) {
906 hlt_met = hlt_trkmht_met_cont->at(0);
907 }
else if (
alg ==
"pfsum" && hlt_pfsum_met_cont.isValid() && hlt_pfsum_met_cont->size() > 0) {
908 hlt_met = hlt_pfsum_met_cont->at(0);
909 }
else if (
alg ==
"pfsum_cssk" && hlt_pfsum_cssk_met_cont.isValid() && hlt_pfsum_cssk_met_cont->size() > 0) {
910 hlt_met = hlt_pfsum_cssk_met_cont->at(0);
911 }
else if (
alg ==
"pfsum_vssk" && hlt_pfsum_vssk_met_cont.isValid() && hlt_pfsum_vssk_met_cont->size() > 0) {
912 hlt_met = hlt_pfsum_vssk_met_cont->at(0);
913 }
else if (
alg ==
"pfopufit" && hlt_pfopufit_met_cont.isValid() && hlt_pfopufit_met_cont->size() > 0) {
914 hlt_met = hlt_pfopufit_met_cont->at(0);
915 }
else if (
alg ==
"cvfpufit" && hlt_cvfpufit_met_cont.isValid() && hlt_cvfpufit_met_cont->size() > 0) {
916 hlt_met = hlt_cvfpufit_met_cont->at(0);
917 }
else if (
alg ==
"mhtpufit_pf" && hlt_mhtpufit_pf_met_cont.isValid() && hlt_mhtpufit_pf_met_cont->size() > 0) {
918 hlt_met = hlt_mhtpufit_pf_met_cont->at(0);
919 }
else if (
alg ==
"mhtpufit_em" && hlt_mhtpufit_em_met_cont.isValid() && hlt_mhtpufit_em_met_cont->size() > 0) {
920 hlt_met = hlt_mhtpufit_em_met_cont->at(0);
929 float hlt_Et = std::sqrt(hlt_Ex*hlt_Ex + hlt_Ey*hlt_Ey);
931 float hlt_Ex_log =
signed_log(hlt_Ex, epsilon);
932 float hlt_Ey_log =
signed_log(hlt_Ey, epsilon);
933 float hlt_Et_log =
signed_log(hlt_Et, epsilon);
934 float hlt_sumEt_log =
signed_log(hlt_sumEt, epsilon);
935 TVector3
v(hlt_Ex, hlt_Ey, hlt_Ez);
936 float hlt_eta =
v.Eta();
937 float hlt_phi =
v.Phi();
951 fill(
tool,met_Ex,met_Ey,met_Et,met_sumEt,
952 met_Ex_log,met_Ey_log,met_Et_log,met_sumEt_log,
968 fill(
tool,met_Ex,met_Ey,met_Et,met_sumEt,
969 met_Ex_log,met_Ey_log,met_Et_log,met_sumEt_log,
977 if (
alg ==
"cell" && hlt_cell_met_cont.isValid() && hlt_cell_met_cont->size() > 0) {
978 hlt_met = hlt_cell_met_cont->at(0);
979 }
else if (
alg ==
"mht" && hlt_mht_met_cont.isValid() && hlt_mht_met_cont->size() > 0) {
980 hlt_met = hlt_mht_met_cont->at(0);
981 }
else if (
alg ==
"tc" && hlt_tc_met_cont.isValid() && hlt_tc_met_cont->size() > 0) {
982 hlt_met = hlt_tc_met_cont->at(0);
983 }
else if (
alg ==
"tc_em" && hlt_tc_em_met_cont.isValid() && hlt_tc_em_met_cont->size() > 0) {
984 hlt_met = hlt_tc_em_met_cont->at(0);
985 }
else if (
alg ==
"tcpufit" && hlt_tcpufit_met_cont.isValid() && hlt_tcpufit_met_cont->size() > 0) {
986 hlt_met = hlt_tcpufit_met_cont->at(0);
987 }
else if (
alg ==
"tcpufit_sig30" && hlt_tcpufit_sig30_met_cont.isValid() && hlt_tcpufit_sig30_met_cont->size() > 0) {
988 hlt_met = hlt_tcpufit_sig30_met_cont->at(0);
989 }
else if (
alg ==
"trkmht" && hlt_trkmht_met_cont.isValid() && hlt_trkmht_met_cont->size() > 0) {
990 hlt_met = hlt_trkmht_met_cont->at(0);
991 }
else if (
alg ==
"pfsum" && hlt_pfsum_met_cont.isValid() && hlt_pfsum_met_cont->size() > 0) {
992 hlt_met = hlt_pfsum_met_cont->at(0);
993 }
else if (
alg ==
"pfsum_cssk" && hlt_pfsum_cssk_met_cont.isValid() && hlt_pfsum_cssk_met_cont->size() > 0) {
994 hlt_met = hlt_pfsum_cssk_met_cont->at(0);
995 }
else if (
alg ==
"pfsum_vssk" && hlt_pfsum_vssk_met_cont.isValid() && hlt_pfsum_vssk_met_cont->size() > 0) {
996 hlt_met = hlt_pfsum_vssk_met_cont->at(0);
997 }
else if (
alg ==
"pfopufit" && hlt_pfopufit_met_cont.isValid() && hlt_pfopufit_met_cont->size() > 0) {
998 hlt_met = hlt_pfopufit_met_cont->at(0);
999 }
else if (
alg ==
"pfopufit_sig30" && hlt_pfopufit_sig30_met_cont.isValid() && hlt_pfopufit_sig30_met_cont->size() > 0) {
1000 hlt_met = hlt_pfopufit_sig30_met_cont->at(0);
1001 }
else if (
alg ==
"cvfpufit" && hlt_cvfpufit_met_cont.isValid() && hlt_cvfpufit_met_cont->size() > 0) {
1002 hlt_met = hlt_cvfpufit_met_cont->at(0);
1003 }
else if (
alg ==
"mhtpufit_pf" && hlt_mhtpufit_pf_met_cont.isValid() && hlt_mhtpufit_pf_met_cont->size() > 0) {
1004 hlt_met = hlt_mhtpufit_pf_met_cont->at(0);
1005 }
else if (
alg ==
"mhtpufit_em" && hlt_mhtpufit_em_met_cont.isValid() && hlt_mhtpufit_em_met_cont->size() > 0) {
1006 hlt_met = hlt_mhtpufit_em_met_cont->at(0);
1007 }
else if (
alg ==
"met_nn" && hlt_met_nn_cont.isValid() && hlt_met_nn_cont->size() > 0) {
1008 hlt_met = hlt_met_nn_cont->at(0);
1016 float hlt_Et = std::sqrt(hlt_Ex*hlt_Ex + hlt_Ey*hlt_Ey);
1018 if (!std::isnan(hlt_Et)) {
1023 fill(
tool,met_Ex,met_Ey,met_Et,met_sumEt);
1025 if (L1_roiMet_Et > 50.) {
1036 if (
alg ==
"cell" && hlt_cell_met_cont.isValid() && hlt_cell_met_cont->size() > 0) {
1037 hlt_met = hlt_cell_met_cont->at(0);
1038 }
else if (
alg ==
"mht" && hlt_mht_met_cont.isValid() && hlt_mht_met_cont->size() > 0) {
1039 hlt_met = hlt_mht_met_cont->at(0);
1040 }
else if (
alg ==
"tc" && hlt_tc_met_cont.isValid() && hlt_tc_met_cont->size() > 0) {
1041 hlt_met = hlt_tc_met_cont->at(0);
1042 }
else if (
alg ==
"tc_em" && hlt_tc_em_met_cont.isValid() && hlt_tc_em_met_cont->size() > 0) {
1043 hlt_met = hlt_tc_em_met_cont->at(0);
1044 }
else if (
alg ==
"trkmht_pf" && hlt_trkmht_met_cont.isValid() && hlt_trkmht_met_cont->size() > 0) {
1045 hlt_met = hlt_trkmht_met_cont->at(0);
1046 }
else if (
alg ==
"pfsum" && hlt_pfsum_met_cont.isValid() && hlt_pfsum_met_cont->size() > 0) {
1047 hlt_met = hlt_pfsum_met_cont->at(0);
1048 }
else if (
alg ==
"pfsum_cssk" && hlt_pfsum_cssk_met_cont.isValid() && hlt_pfsum_cssk_met_cont->size() > 0) {
1049 hlt_met = hlt_pfsum_cssk_met_cont->at(0);
1050 }
else if (
alg ==
"pfsum_vssk" && hlt_pfsum_vssk_met_cont.isValid() && hlt_pfsum_vssk_met_cont->size() > 0) {
1051 hlt_met = hlt_pfsum_vssk_met_cont->at(0);
1052 }
else if (
alg ==
"pfopufit" && hlt_pfopufit_met_cont.isValid() && hlt_pfopufit_met_cont->size() > 0) {
1053 hlt_met = hlt_pfopufit_met_cont->at(0);
1054 }
else if (
alg ==
"cvfpufit" && hlt_cvfpufit_met_cont.isValid() && hlt_cvfpufit_met_cont->size() > 0) {
1055 hlt_met = hlt_cvfpufit_met_cont->at(0);
1056 }
else if (
alg ==
"mhtpufit_pf" && hlt_mhtpufit_pf_met_cont.isValid() && hlt_mhtpufit_pf_met_cont->size() > 0) {
1057 hlt_met = hlt_mhtpufit_pf_met_cont->at(0);
1058 }
else if (
alg ==
"mhtpufit_em" && hlt_mhtpufit_em_met_cont.isValid() && hlt_mhtpufit_em_met_cont->size() > 0) {
1059 hlt_met = hlt_mhtpufit_em_met_cont->at(0);
1064 if(hlt_met && hlt_tcpufit_met_cont.isValid() && hlt_tcpufit_met_cont->size() > 0){
1065 hlt_tcpufit_met = hlt_tcpufit_met_cont->at(0);
1068 float hlt_Et = std::sqrt(hlt_Ex*hlt_Ex + hlt_Ey*hlt_Ey);
1073 float hlt_tcpufit_Et = std::sqrt(hlt_tcpufit_Ex*hlt_tcpufit_Ex + hlt_tcpufit_Ey*hlt_tcpufit_Ey);
1075 fill(
tool, met_Et, tcpufit_met_Et);
1082 fill(
tool, pass_chain,offline_NoMu_Et_eff);
1086 fill(
tool, pass_chain,offline_NoMu_Et_eff);
1090 fill(
tool, pass_chain,offline_NoMu_Et_eff);
1094 fill(
tool, pass_chain,offline_NoMu_Et_eff);
1097 return StatusCode::SUCCESS;
1.8.18