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.
64 if (! eventInfo.isValid() ){
70 if (! hlt_electron_cont.isValid() ) {
75 if (! hlt_muon_cont.isValid() ) {
81 if (! hlt_topoclusters_cont.isValid() ) {
87 if (! hlt_tracks_cont.isValid() ) {
93 if (! hlt_vertex_cont.isValid() ) {
98 if (! offline_vertex_cont.isValid() ) {
104 if (! offline_met_cont.isValid() ) {
110 if (! l1_roi_cont.isValid() ) {
116 if (! l1_jFexMet_cont.isValid() ) {
120 if (! l1_jFexSumEt_cont.isValid() ) {
124 if (! l1_gFexJwojScalar_cont.isValid() ) {
128 if (! l1_gFexJwojMETComponents_cont.isValid() ) {
132 if (! l1_gFexJwojMHTComponents_cont.isValid() ) {
136 if (! l1_gFexJwojMSTComponents_cont.isValid() ) {
140 if (! l1_gFexNCMETScalar_cont.isValid() ) {
144 if (! l1_gFexNCMETComponents_cont.isValid() ) {
148 if (! l1_gFexRhoMETScalar_cont.isValid() ) {
152 if (! l1_gFexRhoMETComponents_cont.isValid() ) {
158 if (! hlt_cell_met_cont.isValid() ) {
163 if (! hlt_mht_met_cont.isValid() ) {
168 if (! hlt_tc_met_cont.isValid() ) {
173 if (! hlt_tc_em_met_cont.isValid() ) {
178 if (! hlt_tcpufit_met_cont.isValid() ) {
183 if (! hlt_tcpufit_sig30_met_cont.isValid() ) {
188 if (! hlt_trkmht_met_cont.isValid() ) {
193 if (! hlt_pfsum_met_cont.isValid() ) {
198 if (! hlt_pfsum_cssk_met_cont.isValid() ) {
203 if (! hlt_pfsum_vssk_met_cont.isValid() ) {
208 if (! hlt_pfopufit_met_cont.isValid() ) {
213 if (! hlt_pfopufit_sig30_met_cont.isValid() ) {
218 if (! hlt_cvfpufit_met_cont.isValid() ) {
223 if (! hlt_mhtpufit_pf_met_cont.isValid() ) {
228 if (! hlt_mhtpufit_em_met_cont.isValid() ) {
233 if (! hlt_met_nn_cont.isValid() ) {
278 double epsilon = 1.189;
286 act_IPBC = eventInfo->actualInteractionsPerCrossing();
290 bool passedPrimaryEl =
false;
293 passedPrimaryEl =
true;
297 if ( hlt_electron_cont.isValid() && passedPrimaryEl ){
298 hlt_el_mult = hlt_electron_cont->size();
300 if( hlt_electron_cont->size() > 0 ) {
301 for (
auto Electron: *hlt_electron_cont) {
309 bool passedPrimaryMu =
false;
312 passedPrimaryMu =
true;
316 if( hlt_muon_cont.isValid() &&passedPrimaryMu ){
317 hlt_mu_mult = hlt_muon_cont->size();
319 if ( hlt_muon_cont->size() > 0 ){
320 for(
auto Muon : *hlt_muon_cont){
328 if(hlt_topoclusters_cont.isValid() && hlt_topoclusters_cont->size() > 0){
329 hlt_topoclusters_mult = hlt_topoclusters_cont->size();
330 for(
auto topoclusters : *hlt_topoclusters_cont){
332 if(hlt_topoclusters_pt > 0){
337 if(hlt_topoclusters_mult > 0){
343 if( hlt_tracks_cont.isValid() && hlt_tracks_cont->size() > 0){
344 hlt_tracks_mult = hlt_tracks_cont->size();
345 float scalarSumPt = 0.0;
346 float scalarSumPx = 0.0;
347 float scalarSumPy = 0.0;
348 for(
auto tracks : *hlt_tracks_cont){
350 hlt_tracks_pt = i_track_pt;
351 scalarSumPt += i_track_pt;
352 scalarSumPx += (tracks->p4().Px());
353 scalarSumPy += (tracks->p4().Py());
355 if(hlt_tracks_pt > hlt_tracks_leading_pt){
356 hlt_tracks_leading_pt = i_track_pt;
359 if(hlt_tracks_pt > 0){
363 if(hlt_tracks_pt > 3){
364 hlt_tracks_eta = (tracks->eta());
365 hlt_tracks_phi = (tracks->phi());
366 fill(
tool, hlt_tracks_eta, hlt_tracks_phi);
370 hlt_tracks_vec_sumPt = std::sqrt(scalarSumPx*scalarSumPx + scalarSumPy*scalarSumPy)/
Gaudi::Units::GeV;
371 hlt_tracks_sca_sumPt = scalarSumPt;
373 fill(
tool, hlt_tracks_mult, hlt_tracks_leading_pt, hlt_tracks_vec_sumPt, hlt_tracks_sca_sumPt);
377 if(hlt_vertex_cont.isValid() && hlt_vertex_cont->size() > 0){
378 hlt_vertex_mult = hlt_vertex_cont->size();
379 hlt_vertex_mult_mu = hlt_vertex_cont->size();
381 for(
auto vertex : *hlt_vertex_cont){
389 hlt_vertex_z = hlt_vertex->
z();
392 if(offline_vertex_cont.isValid() && offline_vertex_cont->size() > 0){
394 for(
auto vertex : *offline_vertex_cont){
402 hlt_vertex_z_diff = hlt_vertex_z - offline_vertex->
z();
408 hlt_vertex_mult = -1;
409 hlt_vertex_mult_mu = 1;
413 fill(
tool, act_IPBC, hlt_vertex_mult_mu);
418 if ( offline_met_cont.isValid() && offline_met_cont->size() > 0 ) {
419 finalTrkMET = ((*offline_met_cont)[
"FinalTrk"]);
420 muonsMET = ((*offline_met_cont)[
"Muons"]);
426 offline_Et = std::sqrt(offline_Ex*offline_Ex + offline_Ey*offline_Ey);
427 offline_Et_eff = std::sqrt(offline_Ex*offline_Ex + offline_Ey*offline_Ey);
428 fill(
tool,offline_Ex,offline_Ey,offline_Et,offline_sumEt);
435 offline_NoMu_Et = std::sqrt(offline_NoMu_Ex*offline_NoMu_Ex + offline_NoMu_Ey*offline_NoMu_Ey);
436 offline_NoMu_Et_eff = std::sqrt(offline_NoMu_Ex*offline_NoMu_Ex + offline_NoMu_Ey*offline_NoMu_Ey);
437 fill(
tool,offline_NoMu_Ex,offline_NoMu_Ey,offline_NoMu_Et,offline_NoMu_sumEt);
445 if (
alg ==
"roi" && l1_roi_cont.isValid()) {
446 l1_met_cont = l1_roi_cont;
450 if ((l1_met_cont->
energyX())>-9e12 && (l1_met_cont->
energyX())<9e12 && (l1_met_cont->
energyY())>-9e12 && (l1_met_cont->
energyY())<9e12) {
453 float L1_met_Et = std::sqrt(L1_met_Ex*L1_met_Ex + L1_met_Ey*L1_met_Ey);
455 float L1_met_Ex_log =
signed_log(L1_met_Ex, epsilon);
456 float L1_met_Ey_log =
signed_log(L1_met_Ey, epsilon);
457 float L1_met_Et_log =
signed_log(L1_met_Et, epsilon);
458 float L1_met_sumEt_log =
signed_log(L1_met_sumEt, epsilon);
459 TVector3
v(L1_met_Ex, L1_met_Ey, 0.0);
460 float L1_met_phi =
v.Phi();
471 fill(
tool, L1_Ex, L1_Ey, L1_Et, L1_sumEt,
472 L1_Ex_log, L1_Ey_log, L1_Et_log, L1_sumEt_log, L1_phi);
478 if (l1_jFexMet_cont.isValid() && l1_jFexMet_cont->size() > 0) {
481 for (
const auto l1_jmet: *l1_jFexMet_cont) {
485 float L1_met_Et = std::sqrt(L1_met_Ex*L1_met_Ex + L1_met_Ey*L1_met_Ey);
486 float L1_met_Ex_log =
signed_log(L1_met_Ex, epsilon);
487 float L1_met_Ey_log =
signed_log(L1_met_Ey, epsilon);
488 float L1_met_Et_log =
signed_log(L1_met_Et, epsilon);
489 TVector3
v(L1_met_Ex, L1_met_Ey, 0.0);
490 float L1_met_phi =
v.Phi();
498 fill(
tool, L1_Ex, L1_Ey, L1_Et, L1_Ex_log, L1_Ey_log, L1_Et_log, L1_phi);
500 if (l1_jFexSumEt_cont.isValid() && l1_jFexSumEt_cont->size() > 0) {
501 float L1_met_sumEt = 0;
502 for (
const auto l1_jsumEt: *l1_jFexSumEt_cont) {
505 float L1_met_sumEt_log =
signed_log(L1_met_sumEt, epsilon);
516 if (l1_gFexJwojScalar_cont.isValid() && l1_gFexJwojScalar_cont->size() > 0) {
517 l1_gmet = l1_gFexJwojScalar_cont->at(0);
520 float L1_met_Et_log =
signed_log(L1_met_Et, epsilon);
521 float L1_met_sumEt_log =
signed_log(L1_met_sumEt, epsilon);
526 fill(
tool, L1_Et, L1_Et_log, L1_sumEt, L1_sumEt_log);
529 if (l1_gFexJwojMETComponents_cont.isValid() && l1_gFexJwojMETComponents_cont->size() > 0) {
530 l1_gmet = l1_gFexJwojMETComponents_cont->at(0);
532 float L1_met_Ex_log =
signed_log(L1_met_Ex, epsilon);
534 float L1_met_Ey_log =
signed_log(L1_met_Ey, epsilon);
535 TVector3
v(L1_met_Ex, L1_met_Ey, 0.0);
536 float L1_met_phi =
v.Phi();
542 fill(
tool, L1_Ex, L1_Ey, L1_Ex_log, L1_Ey_log, L1_phi);
545 if (l1_gFexJwojMHTComponents_cont.isValid() && l1_gFexJwojMHTComponents_cont->size() > 0) {
546 l1_gmet = l1_gFexJwojMHTComponents_cont->at(0);
548 float L1_met_HT_Ex_log =
signed_log(L1_met_HT_Ex, epsilon);
550 float L1_met_HT_Ey_log =
signed_log(L1_met_HT_Ey, epsilon);
551 TVector3
v(L1_met_HT_Ex, L1_met_HT_Ey, 0.0);
552 float L1_met_HT_phi =
v.Phi();
558 fill(
tool, L1_HT_Ex, L1_HT_Ey, L1_HT_Ex_log, L1_HT_Ey_log, L1_HT_phi);
561 if (l1_gFexJwojMSTComponents_cont.isValid() && l1_gFexJwojMSTComponents_cont->size() > 0) {
562 l1_gmet = l1_gFexJwojMSTComponents_cont->at(0);
564 float L1_met_ST_Ex_log =
signed_log(L1_met_ST_Ex, epsilon);
566 float L1_met_ST_Ey_log =
signed_log(L1_met_ST_Ey, epsilon);
567 TVector3
v(L1_met_ST_Ex, L1_met_ST_Ey, 0.0);
568 float L1_met_ST_phi =
v.Phi();
574 fill(
tool, L1_ST_Ex, L1_ST_Ey, L1_ST_Ex_log, L1_ST_Ey_log, L1_ST_phi);
577 if (l1_gFexNCMETScalar_cont.isValid() && l1_gFexNCMETScalar_cont->size() > 0) {
578 l1_gmet = l1_gFexNCMETScalar_cont->at(0);
581 float L1_met_Et_log =
signed_log(L1_met_Et, epsilon);
582 float L1_met_sumEt_log =
signed_log(L1_met_sumEt, epsilon);
587 fill(
tool, L1_Et, L1_Et_log, L1_sumEt, L1_sumEt_log);
590 if (l1_gFexNCMETComponents_cont.isValid() && l1_gFexNCMETComponents_cont->size() > 0) {
591 l1_gmet = l1_gFexNCMETComponents_cont->at(0);
593 float L1_met_Ex_log =
signed_log(L1_met_Ex, epsilon);
595 float L1_met_Ey_log =
signed_log(L1_met_Ey, epsilon);
596 TVector3
v(L1_met_Ex, L1_met_Ey, 0.0);
597 float L1_met_phi =
v.Phi();
603 fill(
tool, L1_Ex, L1_Ey, L1_Ex_log, L1_Ey_log, L1_phi);
606 if (l1_gFexRhoMETScalar_cont.isValid() && l1_gFexRhoMETScalar_cont->size() > 0) {
607 l1_gmet = l1_gFexRhoMETScalar_cont->at(0);
610 float L1_met_Et_log =
signed_log(L1_met_Et, epsilon);
611 float L1_met_sumEt_log =
signed_log(L1_met_sumEt, epsilon);
616 fill(
tool, L1_Et, L1_Et_log, L1_sumEt, L1_sumEt_log);
619 if (l1_gFexRhoMETComponents_cont.isValid() && l1_gFexRhoMETComponents_cont->size() > 0) {
620 l1_gmet = l1_gFexRhoMETComponents_cont->at(0);
622 float L1_met_Ex_log =
signed_log(L1_met_Ex, epsilon);
624 float L1_met_Ey_log =
signed_log(L1_met_Ey, epsilon);
625 TVector3
v(L1_met_Ex, L1_met_Ey, 0.0);
626 float L1_met_phi =
v.Phi();
632 fill(
tool, L1_Ex, L1_Ey, L1_Ex_log, L1_Ey_log, L1_phi);
641 if ( hlt_cell_met_cont.isValid() && hlt_cell_met_cont->size() > 0 ) {
642 hlt_met = hlt_cell_met_cont->at(0);
643 for (
int j=0; j<nStatus; ++j) {
644 unsigned mask = (1
u<<j);
654 for (
int i=0;
i<nComponent; ++
i) {
657 component_Et = sqrt(ex*ex+ey*ey);
662 for (
int i=0;
i<nComponent; ++
i) {
663 for (
int j=0; j<nStatus; ++j) {
664 unsigned mask = (1
u<<j);
666 component_status_weight = 1.;
668 component_status_weight = 0;
672 fill(
tool,mon_comp,mon_bit,component_status_weight);
678 float L1_roiMet_Et = 0;
679 if ( l1_roi_cont.isValid() ) {
680 if ((l1_roi_cont->energyX())>-9e12 && (l1_roi_cont->energyX())<9e12 && (l1_roi_cont->energyY())>-9e12 && (l1_roi_cont->energyY())<9e12) {
683 L1_roiMet_Et = std::sqrt(Ex*Ex + Ey*Ey);
686 float L1_jFexMet_Et = 0;
687 if (l1_jFexMet_cont.isValid() && l1_jFexMet_cont->size() > 0) {
690 for (
const auto l1_jmet: *l1_jFexMet_cont) {
694 L1_jFexMet_Et = std::sqrt(L1_met_Ex*L1_met_Ex + L1_met_Ey*L1_met_Ey);
696 float L1_gFexJWOJMet_Et = 0;
698 float L1_PreSelMet_Et = L1_roiMet_Et;
699 if (
m_L1MetAlg == 1) L1_PreSelMet_Et = L1_jFexMet_Et;
700 if (
m_L1MetAlg == 2) L1_PreSelMet_Et = L1_gFexJWOJMet_Et;
704 if (
alg ==
"cell" && hlt_cell_met_cont.isValid() && hlt_cell_met_cont->size() > 0) {
705 hlt_met = hlt_cell_met_cont->at(0);
706 }
else if (
alg ==
"mht" && hlt_mht_met_cont.isValid() && hlt_mht_met_cont->size() > 0) {
707 hlt_met = hlt_mht_met_cont->at(0);
708 }
else if (
alg ==
"tc" && hlt_tc_met_cont.isValid() && hlt_tc_met_cont->size() > 0) {
709 hlt_met = hlt_tc_met_cont->at(0);
710 }
else if (
alg ==
"tc_em" && hlt_tc_em_met_cont.isValid() && hlt_tc_em_met_cont->size() > 0) {
711 hlt_met = hlt_tc_em_met_cont->at(0);
712 }
else if (
alg ==
"tcpufit" && hlt_tcpufit_met_cont.isValid() && hlt_tcpufit_met_cont->size() > 0) {
713 hlt_met = hlt_tcpufit_met_cont->at(0);
714 }
else if (
alg ==
"tcpufit_sig30" && hlt_tcpufit_sig30_met_cont.isValid() && hlt_tcpufit_sig30_met_cont->size() > 0) {
715 hlt_met = hlt_tcpufit_sig30_met_cont->at(0);
716 }
else if (
alg ==
"trkmht" && hlt_trkmht_met_cont.isValid() && hlt_trkmht_met_cont->size() > 0) {
717 hlt_met = hlt_trkmht_met_cont->at(0);
718 }
else if (
alg ==
"pfsum" && hlt_pfsum_met_cont.isValid() && hlt_pfsum_met_cont->size() > 0) {
719 hlt_met = hlt_pfsum_met_cont->at(0);
720 }
else if (
alg ==
"pfsum_cssk" && hlt_pfsum_cssk_met_cont.isValid() && hlt_pfsum_cssk_met_cont->size() > 0) {
721 hlt_met = hlt_pfsum_cssk_met_cont->at(0);
722 }
else if (
alg ==
"pfsum_vssk" && hlt_pfsum_vssk_met_cont.isValid() && hlt_pfsum_vssk_met_cont->size() > 0) {
723 hlt_met = hlt_pfsum_vssk_met_cont->at(0);
724 }
else if (
alg ==
"pfopufit" && hlt_pfopufit_met_cont.isValid() && hlt_pfopufit_met_cont->size() > 0) {
725 hlt_met = hlt_pfopufit_met_cont->at(0);
726 }
else if (
alg ==
"pfopufit_sig30" && hlt_pfopufit_sig30_met_cont.isValid() && hlt_pfopufit_sig30_met_cont->size() > 0) {
727 hlt_met = hlt_pfopufit_sig30_met_cont->at(0);
728 }
else if (
alg ==
"cvfpufit" && hlt_cvfpufit_met_cont.isValid() && hlt_cvfpufit_met_cont->size() > 0) {
729 hlt_met = hlt_cvfpufit_met_cont->at(0);
730 }
else if (
alg ==
"mhtpufit_pf" && hlt_mhtpufit_pf_met_cont.isValid() && hlt_mhtpufit_pf_met_cont->size() > 0) {
731 hlt_met = hlt_mhtpufit_pf_met_cont->at(0);
732 }
else if (
alg ==
"mhtpufit_em" && hlt_mhtpufit_em_met_cont.isValid() && hlt_mhtpufit_em_met_cont->size() > 0) {
733 hlt_met = hlt_mhtpufit_em_met_cont->at(0);
734 }
else if (
alg ==
"met_nn" && hlt_met_nn_cont.isValid() && hlt_met_nn_cont->size() > 0) {
735 hlt_met = hlt_met_nn_cont->at(0);
744 float hlt_Et = std::sqrt(hlt_Ex*hlt_Ex + hlt_Ey*hlt_Ey);
746 float hlt_Ex_log =
signed_log(hlt_Ex, epsilon);
747 float hlt_Ey_log =
signed_log(hlt_Ey, epsilon);
748 float hlt_Et_log =
signed_log(hlt_Et, epsilon);
749 float hlt_sumEt_log =
signed_log(hlt_sumEt, epsilon);
750 TVector3
v(hlt_Ex, hlt_Ey, hlt_Ez);
751 float hlt_eta =
v.Eta();
752 float hlt_phi =
v.Phi();
764 fill(
tool,met_Ex,met_Ey,met_Et,met_sumEt,
765 met_Ex_log,met_Ey_log,met_Et_log,met_sumEt_log,
773 if (
alg ==
"cell" && hlt_cell_met_cont.isValid() && hlt_cell_met_cont->size() > 0) {
774 hlt_met = hlt_cell_met_cont->at(0);
775 }
else if (
alg ==
"mht" && hlt_mht_met_cont.isValid() && hlt_mht_met_cont->size() > 0) {
776 hlt_met = hlt_mht_met_cont->at(0);
777 }
else if (
alg ==
"tc" && hlt_tc_met_cont.isValid() && hlt_tc_met_cont->size() > 0) {
778 hlt_met = hlt_tc_met_cont->at(0);
779 }
else if (
alg ==
"tc_em" && hlt_tc_em_met_cont.isValid() && hlt_tc_em_met_cont->size() > 0) {
780 hlt_met = hlt_tc_em_met_cont->at(0);
781 }
else if (
alg ==
"tcpufit" && hlt_tcpufit_met_cont.isValid() && hlt_tcpufit_met_cont->size() > 0) {
782 hlt_met = hlt_tcpufit_met_cont->at(0);
783 }
else if (
alg ==
"tcpufit_sig30" && hlt_tcpufit_sig30_met_cont.isValid() && hlt_tcpufit_sig30_met_cont->size() > 0) {
784 hlt_met = hlt_tcpufit_sig30_met_cont->at(0);
785 }
else if (
alg ==
"trkmht" && hlt_trkmht_met_cont.isValid() && hlt_trkmht_met_cont->size() > 0) {
786 hlt_met = hlt_trkmht_met_cont->at(0);
787 }
else if (
alg ==
"pfsum" && hlt_pfsum_met_cont.isValid() && hlt_pfsum_met_cont->size() > 0) {
788 hlt_met = hlt_pfsum_met_cont->at(0);
789 }
else if (
alg ==
"pfsum_cssk" && hlt_pfsum_cssk_met_cont.isValid() && hlt_pfsum_cssk_met_cont->size() > 0) {
790 hlt_met = hlt_pfsum_cssk_met_cont->at(0);
791 }
else if (
alg ==
"pfsum_vssk" && hlt_pfsum_vssk_met_cont.isValid() && hlt_pfsum_vssk_met_cont->size() > 0) {
792 hlt_met = hlt_pfsum_vssk_met_cont->at(0);
793 }
else if (
alg ==
"pfopufit" && hlt_pfopufit_met_cont.isValid() && hlt_pfopufit_met_cont->size() > 0) {
794 hlt_met = hlt_pfopufit_met_cont->at(0);
795 }
else if (
alg ==
"pfopufit_sig30" && hlt_pfopufit_sig30_met_cont.isValid() && hlt_pfopufit_sig30_met_cont->size() > 0) {
796 hlt_met = hlt_pfopufit_sig30_met_cont->at(0);
797 }
else if (
alg ==
"cvfpufit" && hlt_cvfpufit_met_cont.isValid() && hlt_cvfpufit_met_cont->size() > 0) {
798 hlt_met = hlt_cvfpufit_met_cont->at(0);
799 }
else if (
alg ==
"mhtpufit_pf" && hlt_mhtpufit_pf_met_cont.isValid() && hlt_mhtpufit_pf_met_cont->size() > 0) {
800 hlt_met = hlt_mhtpufit_pf_met_cont->at(0);
801 }
else if (
alg ==
"mhtpufit_em" && hlt_mhtpufit_em_met_cont.isValid() && hlt_mhtpufit_em_met_cont->size() > 0) {
802 hlt_met = hlt_mhtpufit_em_met_cont->at(0);
803 }
else if (
alg ==
"met_nn" && hlt_met_nn_cont.isValid() && hlt_met_nn_cont->size() > 0) {
804 hlt_met = hlt_met_nn_cont->at(0);
812 float hlt_Et = std::sqrt(hlt_Ex*hlt_Ex + hlt_Ey*hlt_Ey);
813 if (L1_PreSelMet_Et >
m_L1MetCut && !std::isnan(hlt_Et)) {
821 if (
alg ==
"cell" && hlt_cell_met_cont.isValid() && hlt_cell_met_cont->size() > 0) {
822 hlt_met = hlt_cell_met_cont->at(0);
823 }
else if (
alg ==
"mht" && hlt_mht_met_cont.isValid() && hlt_mht_met_cont->size() > 0) {
824 hlt_met = hlt_mht_met_cont->at(0);
825 }
else if (
alg ==
"tc" && hlt_tc_met_cont.isValid() && hlt_tc_met_cont->size() > 0) {
826 hlt_met = hlt_tc_met_cont->at(0);
827 }
else if (
alg ==
"tc_em" && hlt_tc_em_met_cont.isValid() && hlt_tc_em_met_cont->size() > 0) {
828 hlt_met = hlt_tc_em_met_cont->at(0);
829 }
else if (
alg ==
"tcpufit" && hlt_tcpufit_met_cont.isValid() && hlt_tcpufit_met_cont->size() > 0) {
830 hlt_met = hlt_tcpufit_met_cont->at(0);
831 }
else if (
alg ==
"trkmht" && hlt_trkmht_met_cont.isValid() && hlt_trkmht_met_cont->size() > 0) {
832 hlt_met = hlt_trkmht_met_cont->at(0);
833 }
else if (
alg ==
"pfsum" && hlt_pfsum_met_cont.isValid() && hlt_pfsum_met_cont->size() > 0) {
834 hlt_met = hlt_pfsum_met_cont->at(0);
835 }
else if (
alg ==
"pfsum_cssk" && hlt_pfsum_cssk_met_cont.isValid() && hlt_pfsum_cssk_met_cont->size() > 0) {
836 hlt_met = hlt_pfsum_cssk_met_cont->at(0);
837 }
else if (
alg ==
"pfsum_vssk" && hlt_pfsum_vssk_met_cont.isValid() && hlt_pfsum_vssk_met_cont->size() > 0) {
838 hlt_met = hlt_pfsum_vssk_met_cont->at(0);
839 }
else if (
alg ==
"pfopufit" && hlt_pfopufit_met_cont.isValid() && hlt_pfopufit_met_cont->size() > 0) {
840 hlt_met = hlt_pfopufit_met_cont->at(0);
841 }
else if (
alg ==
"cvfpufit" && hlt_cvfpufit_met_cont.isValid() && hlt_cvfpufit_met_cont->size() > 0) {
842 hlt_met = hlt_cvfpufit_met_cont->at(0);
843 }
else if (
alg ==
"mhtpufit_pf" && hlt_mhtpufit_pf_met_cont.isValid() && hlt_mhtpufit_pf_met_cont->size() > 0) {
844 hlt_met = hlt_mhtpufit_pf_met_cont->at(0);
845 }
else if (
alg ==
"mhtpufit_em" && hlt_mhtpufit_em_met_cont.isValid() && hlt_mhtpufit_em_met_cont->size() > 0) {
846 hlt_met = hlt_mhtpufit_em_met_cont->at(0);
855 float hlt_Et = std::sqrt(hlt_Ex*hlt_Ex + hlt_Ey*hlt_Ey);
857 float hlt_Ex_log =
signed_log(hlt_Ex, epsilon);
858 float hlt_Ey_log =
signed_log(hlt_Ey, epsilon);
859 float hlt_Et_log =
signed_log(hlt_Et, epsilon);
860 float hlt_sumEt_log =
signed_log(hlt_sumEt, epsilon);
861 TVector3
v(hlt_Ex, hlt_Ey, hlt_Ez);
862 float hlt_eta =
v.Eta();
863 float hlt_phi =
v.Phi();
879 fill(
tool,met_Ex,met_Ey,met_Et,met_sumEt,
880 met_Ex_log,met_Ey_log,met_Et_log,met_sumEt_log,
888 if (
alg ==
"cell" && hlt_cell_met_cont.isValid() && hlt_cell_met_cont->size() > 0) {
889 hlt_met = hlt_cell_met_cont->at(0);
890 }
else if (
alg ==
"mht" && hlt_mht_met_cont.isValid() && hlt_mht_met_cont->size() > 0) {
891 hlt_met = hlt_mht_met_cont->at(0);
892 }
else if (
alg ==
"tc" && hlt_tc_met_cont.isValid() && hlt_tc_met_cont->size() > 0) {
893 hlt_met = hlt_tc_met_cont->at(0);
894 }
else if (
alg ==
"tc_em" && hlt_tc_em_met_cont.isValid() && hlt_tc_em_met_cont->size() > 0) {
895 hlt_met = hlt_tc_em_met_cont->at(0);
896 }
else if (
alg ==
"tcpufit" && hlt_tcpufit_met_cont.isValid() && hlt_tcpufit_met_cont->size() > 0) {
897 hlt_met = hlt_tcpufit_met_cont->at(0);
898 }
else if (
alg ==
"trkmht" && hlt_trkmht_met_cont.isValid() && hlt_trkmht_met_cont->size() > 0) {
899 hlt_met = hlt_trkmht_met_cont->at(0);
900 }
else if (
alg ==
"pfsum" && hlt_pfsum_met_cont.isValid() && hlt_pfsum_met_cont->size() > 0) {
901 hlt_met = hlt_pfsum_met_cont->at(0);
902 }
else if (
alg ==
"pfsum_cssk" && hlt_pfsum_cssk_met_cont.isValid() && hlt_pfsum_cssk_met_cont->size() > 0) {
903 hlt_met = hlt_pfsum_cssk_met_cont->at(0);
904 }
else if (
alg ==
"pfsum_vssk" && hlt_pfsum_vssk_met_cont.isValid() && hlt_pfsum_vssk_met_cont->size() > 0) {
905 hlt_met = hlt_pfsum_vssk_met_cont->at(0);
906 }
else if (
alg ==
"pfopufit" && hlt_pfopufit_met_cont.isValid() && hlt_pfopufit_met_cont->size() > 0) {
907 hlt_met = hlt_pfopufit_met_cont->at(0);
908 }
else if (
alg ==
"cvfpufit" && hlt_cvfpufit_met_cont.isValid() && hlt_cvfpufit_met_cont->size() > 0) {
909 hlt_met = hlt_cvfpufit_met_cont->at(0);
910 }
else if (
alg ==
"mhtpufit_pf" && hlt_mhtpufit_pf_met_cont.isValid() && hlt_mhtpufit_pf_met_cont->size() > 0) {
911 hlt_met = hlt_mhtpufit_pf_met_cont->at(0);
912 }
else if (
alg ==
"mhtpufit_em" && hlt_mhtpufit_em_met_cont.isValid() && hlt_mhtpufit_em_met_cont->size() > 0) {
913 hlt_met = hlt_mhtpufit_em_met_cont->at(0);
922 float hlt_Et = std::sqrt(hlt_Ex*hlt_Ex + hlt_Ey*hlt_Ey);
924 float hlt_Ex_log =
signed_log(hlt_Ex, epsilon);
925 float hlt_Ey_log =
signed_log(hlt_Ey, epsilon);
926 float hlt_Et_log =
signed_log(hlt_Et, epsilon);
927 float hlt_sumEt_log =
signed_log(hlt_sumEt, epsilon);
928 TVector3
v(hlt_Ex, hlt_Ey, hlt_Ez);
929 float hlt_eta =
v.Eta();
930 float hlt_phi =
v.Phi();
944 fill(
tool,met_Ex,met_Ey,met_Et,met_sumEt,
945 met_Ex_log,met_Ey_log,met_Et_log,met_sumEt_log,
961 fill(
tool,met_Ex,met_Ey,met_Et,met_sumEt,
962 met_Ex_log,met_Ey_log,met_Et_log,met_sumEt_log,
970 if (
alg ==
"cell" && hlt_cell_met_cont.isValid() && hlt_cell_met_cont->size() > 0) {
971 hlt_met = hlt_cell_met_cont->at(0);
972 }
else if (
alg ==
"mht" && hlt_mht_met_cont.isValid() && hlt_mht_met_cont->size() > 0) {
973 hlt_met = hlt_mht_met_cont->at(0);
974 }
else if (
alg ==
"tc" && hlt_tc_met_cont.isValid() && hlt_tc_met_cont->size() > 0) {
975 hlt_met = hlt_tc_met_cont->at(0);
976 }
else if (
alg ==
"tc_em" && hlt_tc_em_met_cont.isValid() && hlt_tc_em_met_cont->size() > 0) {
977 hlt_met = hlt_tc_em_met_cont->at(0);
978 }
else if (
alg ==
"tcpufit" && hlt_tcpufit_met_cont.isValid() && hlt_tcpufit_met_cont->size() > 0) {
979 hlt_met = hlt_tcpufit_met_cont->at(0);
980 }
else if (
alg ==
"tcpufit_sig30" && hlt_tcpufit_sig30_met_cont.isValid() && hlt_tcpufit_sig30_met_cont->size() > 0) {
981 hlt_met = hlt_tcpufit_sig30_met_cont->at(0);
982 }
else if (
alg ==
"trkmht" && hlt_trkmht_met_cont.isValid() && hlt_trkmht_met_cont->size() > 0) {
983 hlt_met = hlt_trkmht_met_cont->at(0);
984 }
else if (
alg ==
"pfsum" && hlt_pfsum_met_cont.isValid() && hlt_pfsum_met_cont->size() > 0) {
985 hlt_met = hlt_pfsum_met_cont->at(0);
986 }
else if (
alg ==
"pfsum_cssk" && hlt_pfsum_cssk_met_cont.isValid() && hlt_pfsum_cssk_met_cont->size() > 0) {
987 hlt_met = hlt_pfsum_cssk_met_cont->at(0);
988 }
else if (
alg ==
"pfsum_vssk" && hlt_pfsum_vssk_met_cont.isValid() && hlt_pfsum_vssk_met_cont->size() > 0) {
989 hlt_met = hlt_pfsum_vssk_met_cont->at(0);
990 }
else if (
alg ==
"pfopufit" && hlt_pfopufit_met_cont.isValid() && hlt_pfopufit_met_cont->size() > 0) {
991 hlt_met = hlt_pfopufit_met_cont->at(0);
992 }
else if (
alg ==
"pfopufit_sig30" && hlt_pfopufit_sig30_met_cont.isValid() && hlt_pfopufit_sig30_met_cont->size() > 0) {
993 hlt_met = hlt_pfopufit_sig30_met_cont->at(0);
994 }
else if (
alg ==
"cvfpufit" && hlt_cvfpufit_met_cont.isValid() && hlt_cvfpufit_met_cont->size() > 0) {
995 hlt_met = hlt_cvfpufit_met_cont->at(0);
996 }
else if (
alg ==
"mhtpufit_pf" && hlt_mhtpufit_pf_met_cont.isValid() && hlt_mhtpufit_pf_met_cont->size() > 0) {
997 hlt_met = hlt_mhtpufit_pf_met_cont->at(0);
998 }
else if (
alg ==
"mhtpufit_em" && hlt_mhtpufit_em_met_cont.isValid() && hlt_mhtpufit_em_met_cont->size() > 0) {
999 hlt_met = hlt_mhtpufit_em_met_cont->at(0);
1000 }
else if (
alg ==
"met_nn" && hlt_met_nn_cont.isValid() && hlt_met_nn_cont->size() > 0) {
1001 hlt_met = hlt_met_nn_cont->at(0);
1009 float hlt_Et = std::sqrt(hlt_Ex*hlt_Ex + hlt_Ey*hlt_Ey);
1011 if (!std::isnan(hlt_Et)) {
1016 fill(
tool,met_Ex,met_Ey,met_Et,met_sumEt);
1018 if (L1_roiMet_Et > 50.) {
1029 if (
alg ==
"cell" && hlt_cell_met_cont.isValid() && hlt_cell_met_cont->size() > 0) {
1030 hlt_met = hlt_cell_met_cont->at(0);
1031 }
else if (
alg ==
"mht" && hlt_mht_met_cont.isValid() && hlt_mht_met_cont->size() > 0) {
1032 hlt_met = hlt_mht_met_cont->at(0);
1033 }
else if (
alg ==
"tc" && hlt_tc_met_cont.isValid() && hlt_tc_met_cont->size() > 0) {
1034 hlt_met = hlt_tc_met_cont->at(0);
1035 }
else if (
alg ==
"tc_em" && hlt_tc_em_met_cont.isValid() && hlt_tc_em_met_cont->size() > 0) {
1036 hlt_met = hlt_tc_em_met_cont->at(0);
1037 }
else if (
alg ==
"trkmht_pf" && hlt_trkmht_met_cont.isValid() && hlt_trkmht_met_cont->size() > 0) {
1038 hlt_met = hlt_trkmht_met_cont->at(0);
1039 }
else if (
alg ==
"pfsum" && hlt_pfsum_met_cont.isValid() && hlt_pfsum_met_cont->size() > 0) {
1040 hlt_met = hlt_pfsum_met_cont->at(0);
1041 }
else if (
alg ==
"pfsum_cssk" && hlt_pfsum_cssk_met_cont.isValid() && hlt_pfsum_cssk_met_cont->size() > 0) {
1042 hlt_met = hlt_pfsum_cssk_met_cont->at(0);
1043 }
else if (
alg ==
"pfsum_vssk" && hlt_pfsum_vssk_met_cont.isValid() && hlt_pfsum_vssk_met_cont->size() > 0) {
1044 hlt_met = hlt_pfsum_vssk_met_cont->at(0);
1045 }
else if (
alg ==
"pfopufit" && hlt_pfopufit_met_cont.isValid() && hlt_pfopufit_met_cont->size() > 0) {
1046 hlt_met = hlt_pfopufit_met_cont->at(0);
1047 }
else if (
alg ==
"cvfpufit" && hlt_cvfpufit_met_cont.isValid() && hlt_cvfpufit_met_cont->size() > 0) {
1048 hlt_met = hlt_cvfpufit_met_cont->at(0);
1049 }
else if (
alg ==
"mhtpufit_pf" && hlt_mhtpufit_pf_met_cont.isValid() && hlt_mhtpufit_pf_met_cont->size() > 0) {
1050 hlt_met = hlt_mhtpufit_pf_met_cont->at(0);
1051 }
else if (
alg ==
"mhtpufit_em" && hlt_mhtpufit_em_met_cont.isValid() && hlt_mhtpufit_em_met_cont->size() > 0) {
1052 hlt_met = hlt_mhtpufit_em_met_cont->at(0);
1057 if(hlt_met && hlt_tcpufit_met_cont.isValid() && hlt_tcpufit_met_cont->size() > 0){
1058 hlt_tcpufit_met = hlt_tcpufit_met_cont->at(0);
1061 float hlt_Et = std::sqrt(hlt_Ex*hlt_Ex + hlt_Ey*hlt_Ey);
1066 float hlt_tcpufit_Et = std::sqrt(hlt_tcpufit_Ex*hlt_tcpufit_Ex + hlt_tcpufit_Ey*hlt_tcpufit_Ey);
1068 fill(
tool, met_Et, tcpufit_met_Et);
1075 fill(
tool, pass_chain,offline_NoMu_Et_eff);
1079 fill(
tool, pass_chain,offline_NoMu_Et_eff);
1083 fill(
tool, pass_chain,offline_NoMu_Et_eff);
1087 fill(
tool, pass_chain,offline_NoMu_Et_eff);
1090 return StatusCode::SUCCESS;
1.8.18