43 using namespace trigger;
48 LogDebug(
"JetMETHLTOfflineSource") <<
"constructor....";
52 LogDebug(
"JetMETHLTOfflineSource") <<
"unabel to get DQMStore service?";
107 std::vector<edm::ParameterSet>
paths = iConfig.
getParameter<std::vector<edm::ParameterSet> >(
"pathPairs");
108 for(std::vector<edm::ParameterSet>::iterator pathconf = paths.begin() ; pathconf != paths.end(); pathconf++) {
111 pathconf->getParameter<
std::string>(
"denompathname")
133 edm::LogInfo(
"FourVectorHLTOffline") <<
"TriggerResults not found, "
147 edm::LogInfo(
"CaloMETHLTOfflineSource") <<
"TriggerResults::HLT not found, "
148 "automatically select events";
157 edm::LogInfo(
"FourVectorHLTOffline") <<
"TriggerEvent not found, "
185 for(
int i=0;
i<2;
i++){
207 CaloJetCollection::const_iterator calojet_ =
calojet.begin();
208 for(; calojet_ !=
calojet.end(); ++calojet_){
248 PFJetCollection::const_iterator pfjet_ =
pfjet.begin();
249 for(; pfjet_ !=
pfjet.end(); ++pfjet_){
253 if(scale*pfjet_->pt()>
PFJetPt[0]){
264 PFJetPt[0] = scale*pfjet_->pt();
265 PFJetPx[0] = scale*pfjet_->px();
266 PFJetPy[0] = scale*pfjet_->py();
269 PFJetNHEF[0] = pfjet_->neutralHadronEnergyFraction();
270 PFJetCHEF[0] = pfjet_->chargedHadronEnergyFraction();
271 PFJetNEMF[0] = pfjet_->neutralEmEnergyFraction();
272 PFJetCEMF[0] = pfjet_->chargedEmEnergyFraction();
274 else if(scale*pfjet_->pt()<
PFJetPt[0] && scale*pfjet_->pt()>
PFJetPt[1] ){
275 PFJetPt[1] = scale*pfjet_->pt();
276 PFJetPx[1] = scale*pfjet_->px();
277 PFJetPy[1] = scale*pfjet_->py();
280 PFJetNHEF[1] = pfjet_->neutralHadronEnergyFraction();
281 PFJetCHEF[1] = pfjet_->chargedHadronEnergyFraction();
282 PFJetNEMF[1] = pfjet_->neutralEmEnergyFraction();
283 PFJetCEMF[1] = pfjet_->chargedEmEnergyFraction();
289 for(
int i = 0;
i<2;
i++){
293 for(
int i = 0;
i<2;
i++){
316 for (
unsigned int j=0;
j!=nPath; ++
j) {
321 if(
verbose_)
cout<<
"fillMEforMonTriggerSummary: Muon Match"<<endl;
332 for (
unsigned int j=0;
j!=nPath; ++
j) {
337 if(
verbose_)
cout<<
"fillMEforMonTriggerSummary: MinBias Match"<<endl;
347 bool trigFirst=
false;
350 if(!trigFirst)
continue;
363 for(PathInfoCollection::iterator
w = v+1;
w!=
hltPathsAll_.end(); ++
w ){
364 bool trigSec =
false;
367 if(trigSec && trigFirst){
372 if(!trigSec && trigFirst){
384 for (VertexCollection::const_iterator itv=Vtx->begin(); itv!=Vtx->end(); itv++){
406 v->getMEhisto_TriggerSummary()->Fill(0.);
409 bool l1found =
false;
410 if(l1Index < triggerObj_->sizeFilters() ) l1found =
true;
411 if(!l1found)
v->getMEhisto_TriggerSummary()->Fill(1.);
412 if(!l1found && !(
triggerResults_->accept(index)))
v->getMEhisto_TriggerSummary()->Fill(2.);
413 if(!l1found && (
triggerResults_->accept(index)))
v->getMEhisto_TriggerSummary()->Fill(3.);
414 if(l1found)
v->getMEhisto_TriggerSummary()->Fill(4.);
415 if(l1found && (
triggerResults_->accept(index)))
v->getMEhisto_TriggerSummary()->Fill(5.);
416 if(l1found && !(
triggerResults_->accept(index)))
v->getMEhisto_TriggerSummary()->Fill(6.);
420 CaloJetCollection::const_iterator
jet =
calojet.begin();
421 v->getMEhisto_JetPt()->Fill(jet->pt());
422 v->getMEhisto_EtavsPt()->Fill(jet->eta(),jet->pt());
423 v->getMEhisto_PhivsPt()->Fill(jet->phi(),jet->pt());
428 v->getMEhisto_JetSize()->Fill(
calojet.size());
430 CaloJetCollection::const_iterator
jet =
calojet.begin();
431 CaloJetCollection::const_iterator jet2=
calojet.begin(); jet2++;
434 CaloJetCollection::const_iterator jet3 = jet2++;
437 v->getMEhisto_Pt12()->Fill((jet->pt()+jet2->pt())/2.);
438 v->getMEhisto_Eta12()->Fill((jet->eta()+jet2->eta())/2.);
439 v->getMEhisto_Phi12()->Fill(
deltaPhi(jet->phi(),jet2->phi()));
440 v->getMEhisto_Pt3()->Fill(jet3pt);
441 v->getMEhisto_Pt12Pt3()->Fill((jet->pt()+jet2->pt())/2., jet3pt);
442 v->getMEhisto_Pt12Phi12()->Fill((jet->pt()+jet2->pt())/2.,
deltaPhi(jet->phi(),jet2->phi()));
446 if(((
v->getTriggerType().compare(
"MET_Trigger") == 0)|| (
v->getTriggerType().compare(
"TET_Trigger") == 0)) &&
calometColl_.
isValid() ){
449 v->getMEhisto_JetPt()->Fill(met.
pt());
469 std::vector<double>jetPtVec;
470 std::vector<double>jetPhiVec;
471 std::vector<double>jetEtaVec;
472 std::vector<double>jetPxVec;
473 std::vector<double>jetPyVec;
474 std::vector<double>hltPtVec;
475 std::vector<double>hltPhiVec;
476 std::vector<double>hltEtaVec;
477 std::vector<double>hltPxVec;
478 std::vector<double>hltPyVec;
489 bool hltTrigBool =
false;
490 bool diJetFire =
false;
494 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index "<< l1Index <<
" of that name "<<
l1Tag;
501 && v->getTriggerType().compare(
"SingleJet_Trigger") == 0)
502 v->getMEhisto_N_L1()->Fill(kl1.size());
504 trigger::Keys::const_iterator ki = kl1.begin();
505 for(; ki != kl1.end(); ++ki){
506 double l1TrigEta = -100;
507 double l1TrigPhi = -100;
510 l1TrigEta = toc[*ki].eta();
511 l1TrigPhi = toc[*ki].phi();
512 if(v->getTriggerType().compare(
"SingleJet_Trigger") == 0){
513 v->getMEhisto_Pt_L1()->Fill(toc[*ki].
pt());
514 if (
isBarrel(toc[*ki].
eta())) v->getMEhisto_PtBarrel_L1()->Fill(toc[*ki].
pt());
515 if (
isEndCap(toc[*ki].
eta())) v->getMEhisto_PtEndcap_L1()->Fill(toc[*ki].
pt());
516 if (
isForward(toc[*ki].
eta())) v->getMEhisto_PtForward_L1()->Fill(toc[*ki].
pt());
517 v->getMEhisto_Eta_L1()->Fill(toc[*ki].
eta());
518 v->getMEhisto_Phi_L1()->Fill(toc[*ki].
phi());
519 v->getMEhisto_EtaPhi_L1()->Fill(toc[*ki].
eta(),toc[*ki].
phi());
524 v->getMEhisto_Pt_L1()->Fill(toc[*ki].
pt());
525 v->getMEhisto_Phi_L1()->Fill(toc[*ki].
phi());
530 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
536 && v->getTriggerType().compare(
"SingleJet_Trigger") == 0)
537 v->getMEhisto_N_HLT()->Fill(khlt.size());
539 trigger::Keys::const_iterator kj = khlt.begin();
541 for(;kj != khlt.end(); ++kj){
543 double hltTrigEta = -100;
544 double hltTrigPhi = -100;
545 hltTrigEta = toc[*kj].eta();
546 hltTrigPhi = toc[*kj].phi();
547 if((
deltaR(hltTrigEta, hltTrigPhi, l1TrigEta, l1TrigPhi)) < 0.4
548 && (v->getTriggerType().compare(
"DiJet_Trigger") == 0))
554 for(;kj != khlt.end(); ++kj){
555 double hltTrigEta = -100.;
556 double hltTrigPhi = -100.;
560 v->getMEhisto_Pt_HLT()->Fill(toc[*kj].
pt());
561 v->getMEhisto_Phi_HLT()->Fill(toc[*kj].
phi());
562 v->getMEhisto_PtCorrelation_L1HLT()->Fill(toc[*ki].
pt(),toc[*kj].
pt());
563 v->getMEhisto_PhiCorrelation_L1HLT()->Fill(toc[*ki].
phi(),toc[*kj].
phi());
564 v->getMEhisto_PtResolution_L1HLT()->Fill((toc[*ki].
pt()-toc[*kj].
pt())/(toc[*ki].
pt()));
565 v->getMEhisto_PhiResolution_L1HLT()->Fill((toc[*ki].
phi()-toc[*kj].
phi())/(toc[*ki].
phi()));
569 hltTrigEta = toc[*kj].eta();
570 hltTrigPhi = toc[*kj].phi();
571 if((
deltaR(hltTrigEta, hltTrigPhi, l1TrigEta, l1TrigPhi)) < 0.4){
572 if(v->getTriggerType().compare(
"SingleJet_Trigger") == 0){
573 v->getMEhisto_PtCorrelation_L1HLT()->Fill(toc[*ki].
pt(),toc[*kj].
pt());
574 v->getMEhisto_EtaCorrelation_L1HLT()->Fill(toc[*ki].
eta(),toc[*kj].
eta());
575 v->getMEhisto_PhiCorrelation_L1HLT()->Fill(toc[*ki].
phi(),toc[*kj].
phi());
576 v->getMEhisto_PtResolution_L1HLT()->Fill((toc[*ki].
pt()-toc[*kj].
pt())/(toc[*ki].
pt()));
577 v->getMEhisto_EtaResolution_L1HLT()->Fill((toc[*ki].
eta()-toc[*kj].
eta())/(toc[*ki].
eta()));
578 v->getMEhisto_PhiResolution_L1HLT()->Fill((toc[*ki].
phi()-toc[*kj].
phi())/(toc[*ki].
phi()));
581 if(((
deltaR(hltTrigEta, hltTrigPhi, l1TrigEta, l1TrigPhi) < 0.4 )
582 || ((v->getTriggerType().compare(
"DiJet_Trigger") == 0) && hltTrigBool)) && !diJetFire){
583 if(v->getTriggerType().compare(
"SingleJet_Trigger") == 0){
584 v->getMEhisto_Pt_HLT()->Fill(toc[*kj].
pt());
585 if (
isBarrel(toc[*kj].
eta())) v->getMEhisto_PtBarrel_HLT()->Fill(toc[*kj].
pt());
586 if (
isEndCap(toc[*kj].
eta())) v->getMEhisto_PtEndcap_HLT()->Fill(toc[*kj].
pt());
587 if (
isForward(toc[*kj].
eta())) v->getMEhisto_PtForward_HLT()->Fill(toc[*kj].
pt());
588 v->getMEhisto_Eta_HLT()->Fill(toc[*kj].
eta());
589 v->getMEhisto_Phi_HLT()->Fill(toc[*kj].
phi());
590 v->getMEhisto_EtaPhi_HLT()->Fill(toc[*kj].
eta(),toc[*kj].
phi());
596 && (v->getPath().compare(
"PFJet") == 0)){
599 for(
int iCalo=0; iCalo<2; iCalo++){
602 if(v->getTriggerType().compare(
"SingleJet_Trigger") == 0){
603 v->getMEhisto_Pt()->Fill(
CaloJetPt[iCalo]);
612 v->getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*kj].
pt(),
CaloJetPt[iCalo]);
613 v->getMEhisto_EtaCorrelation_HLTRecObj()->Fill(toc[*kj].
eta(),
CaloJetEta[iCalo]);
614 v->getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*kj].
phi(),
CaloJetPhi[iCalo]);
616 v->getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*kj].
pt()-
CaloJetPt[iCalo])/(toc[*kj].
pt()));
617 v->getMEhisto_EtaResolution_HLTRecObj()->Fill((toc[*kj].
eta()-
CaloJetEta[iCalo])/(toc[*kj].
eta()));
618 v->getMEhisto_PhiResolution_HLTRecObj()->Fill((toc[*kj].
phi()-
CaloJetPhi[iCalo])/(toc[*kj].
phi()));
622 if((v->getTriggerType().compare(
"DiJet_Trigger") == 0)){
629 hltPhiVec.push_back(toc[*kj].
phi());
630 hltPtVec.push_back(toc[*kj].
pt());
631 hltEtaVec.push_back(toc[*kj].
eta());
632 hltPxVec.push_back(toc[*kj].px());
633 hltPyVec.push_back(toc[*kj].py());
642 && (v->getPath().compare(
"PFJet") != 0)){
645 for(
int iPF=0; iPF<2; iPF++){
648 if(v->getTriggerType().compare(
"SingleJet_Trigger") == 0){
649 v->getMEhisto_Pt()->Fill(
PFJetPt[iPF]);
654 v->getMEhisto_Eta()->Fill(
PFJetEta[iPF]);
655 v->getMEhisto_Phi()->Fill(
PFJetPhi[iPF]);
658 v->getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*kj].
pt(),
PFJetPt[iPF]);
659 v->getMEhisto_EtaCorrelation_HLTRecObj()->Fill(toc[*kj].
eta(),
PFJetEta[iPF]);
660 v->getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*kj].
phi(),
PFJetPhi[iPF]);
662 v->getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*kj].
pt()-
PFJetPt[iPF])/(toc[*kj].
pt()));
663 v->getMEhisto_EtaResolution_HLTRecObj()->Fill((toc[*kj].
eta()-
PFJetEta[iPF])/(toc[*kj].
eta()));
664 v->getMEhisto_PhiResolution_HLTRecObj()->Fill((toc[*kj].
phi()-
PFJetPhi[iPF])/(toc[*kj].
phi()));
668 if((v->getTriggerType().compare(
"DiJet_Trigger") == 0)){
670 jetPtVec.push_back(
PFJetPt[iPF]);
672 jetPxVec.push_back(
PFJetPx[iPF]);
673 jetPyVec.push_back(
PFJetPy[iPF]);
675 hltPhiVec.push_back(toc[*kj].
phi());
676 hltPtVec.push_back(toc[*kj].
pt());
677 hltEtaVec.push_back(toc[*kj].
eta());
678 hltPxVec.push_back(toc[*kj].px());
679 hltPyVec.push_back(toc[*kj].py());
691 && (v->getPath().find(
"HLT_PFMET")==std::string::npos)){
695 v->getMEhisto_Pt()->Fill(met.
et());
696 v->getMEhisto_Phi()->Fill(met.
phi());
698 v->getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*kj].et(),met.
et());
699 v->getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*kj].
phi(),met.
phi());
700 v->getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*kj].et()-met.
et())/(toc[*kj].et()));
701 v->getMEhisto_PhiResolution_HLTRecObj()->Fill((toc[*kj].
phi()-met.
phi())/(toc[*kj].
phi()));
707 && (v->getPath().find(
"HLT_PFMET")!=std::string::npos)){
709 const PFMET pfmet = pfmetcol->front();
711 v->getMEhisto_Pt()->Fill(pfmet.
et());
712 v->getMEhisto_Phi()->Fill(pfmet.
phi());
714 v->getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*kj].et(),pfmet.
et());
715 v->getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*kj].
phi(),pfmet.
phi());
716 v->getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*kj].et()-pfmet.
et())/(toc[*kj].et()));
717 v->getMEhisto_PhiResolution_HLTRecObj()->Fill((toc[*kj].
phi()-pfmet.
phi())/(toc[*kj].
phi()));
720 if((v->getTriggerType().compare(
"DiJet_Trigger") == 0)) diJetFire =
true;
724 v->getMEhisto_N()->Fill(jetsize);
727 if((v->getTriggerType().compare(
"DiJet_Trigger") == 0) && jetPtVec.size() >1){
728 double AveJetPt = (jetPtVec[0] + jetPtVec[1])/2;
729 double AveJetEta = (jetEtaVec[0] + jetEtaVec[1])/2;
730 double JetDelPhi =
deltaPhi(jetPhiVec[0],jetPhiVec[1]);
731 double AveHLTPt = (hltPtVec[0] + hltPtVec[1])/2;
732 double AveHLTEta = (hltEtaVec[0] + hltEtaVec[1])/2;
733 double HLTDelPhi =
deltaPhi(hltPhiVec[0],hltPhiVec[1]);
734 v->getMEhisto_AveragePt_RecObj()->Fill(AveJetPt);
735 v->getMEhisto_AverageEta_RecObj()->Fill(AveJetEta);
736 v->getMEhisto_DeltaPhi_RecObj()->Fill(JetDelPhi);
738 v->getMEhisto_AveragePt_HLTObj()->Fill(AveHLTPt);
739 v->getMEhisto_AverageEta_HLTObj()->Fill(AveHLTEta);
740 v->getMEhisto_DeltaPhi_HLTObj()->Fill(HLTDelPhi);
812 for (
unsigned int j=0;
j!=nPath; ++
j) {
817 if(
verbose_)
cout<<
"fillMEforMonAllTriggerwrtMuonTrigger: Muon Match"<<endl;
832 std::vector<double>jetPtVec;
833 std::vector<double>jetPhiVec;
834 std::vector<double>jetEtaVec;
835 std::vector<double>jetPxVec;
836 std::vector<double>jetPyVec;
837 std::vector<double>hltPtVec;
838 std::vector<double>hltPhiVec;
839 std::vector<double>hltEtaVec;
840 std::vector<double>hltPxVec;
841 std::vector<double>hltPyVec;
844 bool fillL1HLT =
false;
852 bool hltTrigBool =
false;
853 bool diJetFire =
false;
857 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index "<< l1Index <<
" of that name "<<
l1Tag;
863 trigger::Keys::const_iterator ki = kl1.begin();
864 for(; ki != kl1.end(); ++ki){
865 double l1TrigEta = -100;
866 double l1TrigPhi = -100;
870 l1TrigEta = toc[*ki].eta();
871 l1TrigPhi = toc[*ki].phi();
872 v->getMEhisto_Pt_L1()->Fill(toc[*ki].
pt());
873 if (
isBarrel(toc[*ki].
eta())) v->getMEhisto_PtBarrel_L1()->Fill(toc[*ki].
pt());
874 if (
isEndCap(toc[*ki].
eta())) v->getMEhisto_PtEndcap_L1()->Fill(toc[*ki].
pt());
875 if (
isForward(toc[*ki].
eta())) v->getMEhisto_PtForward_L1()->Fill(toc[*ki].
pt());
876 v->getMEhisto_Eta_L1()->Fill(toc[*ki].
eta());
877 v->getMEhisto_Phi_L1()->Fill(toc[*ki].
phi());
878 v->getMEhisto_EtaPhi_L1()->Fill(toc[*ki].
eta(),toc[*ki].
phi());
881 v->getMEhisto_Pt_L1()->Fill(toc[*ki].
pt());
882 v->getMEhisto_Phi_L1()->Fill(toc[*ki].
phi());
887 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
892 v->getMEhisto_N_HLT()->Fill(khlt.size());
893 trigger::Keys::const_iterator kj = khlt.begin();
895 for(;kj != khlt.end(); ++kj){
897 double hltTrigEta = -100;
898 double hltTrigPhi = -100;
899 hltTrigEta = toc[*kj].eta();
900 hltTrigPhi = toc[*kj].phi();
901 if((
deltaR(hltTrigEta, hltTrigPhi, l1TrigEta, l1TrigPhi)) < 0.4
902 && (v->getTriggerType().compare(
"DiJet_Trigger") == 0))
908 for(;kj != khlt.end(); ++kj){
909 double hltTrigEta = -100.;
910 double hltTrigPhi = -100.;
914 v->getMEhisto_Pt_HLT()->Fill(toc[*kj].
pt());
915 v->getMEhisto_Phi_HLT()->Fill(toc[*kj].
phi());
916 v->getMEhisto_PtCorrelation_L1HLT()->Fill(toc[*ki].
pt(),toc[*kj].
pt());
917 v->getMEhisto_PhiCorrelation_L1HLT()->Fill(toc[*ki].
phi(),toc[*kj].
phi());
918 v->getMEhisto_PtResolution_L1HLT()->Fill((toc[*ki].
pt()-toc[*kj].
pt())/(toc[*ki].
pt()));
919 v->getMEhisto_PhiResolution_L1HLT()->Fill((toc[*ki].
phi()-toc[*kj].
phi())/(toc[*ki].
phi()));
923 hltTrigEta = toc[*kj].eta();
924 hltTrigPhi = toc[*kj].phi();
925 if((
deltaR(hltTrigEta, hltTrigPhi, l1TrigEta, l1TrigPhi)) < 0.4){
926 v->getMEhisto_PtCorrelation_L1HLT()->Fill(toc[*ki].
pt(),toc[*kj].
pt());
927 v->getMEhisto_EtaCorrelation_L1HLT()->Fill(toc[*ki].
eta(),toc[*kj].
eta());
928 v->getMEhisto_PhiCorrelation_L1HLT()->Fill(toc[*ki].
phi(),toc[*kj].
phi());
929 v->getMEhisto_PtResolution_L1HLT()->Fill((toc[*ki].
pt()-toc[*kj].
pt())/(toc[*ki].
pt()));
930 v->getMEhisto_EtaResolution_L1HLT()->Fill((toc[*ki].
eta()-toc[*kj].
eta())/(toc[*ki].
eta()));
931 v->getMEhisto_PhiResolution_L1HLT()->Fill((toc[*ki].
phi()-toc[*kj].
phi())/(toc[*ki].
phi()));
933 if(((
deltaR(hltTrigEta, hltTrigPhi, l1TrigEta, l1TrigPhi) < 0.4 )
934 || ((v->getTriggerType().compare(
"DiJet_Trigger") == 0) && hltTrigBool)) && !diJetFire){
935 v->getMEhisto_Pt_HLT()->Fill(toc[*kj].
pt());
936 if (
isBarrel(toc[*kj].
eta())) v->getMEhisto_PtBarrel_HLT()->Fill(toc[*kj].
pt());
937 if (
isEndCap(toc[*kj].
eta())) v->getMEhisto_PtEndcap_HLT()->Fill(toc[*kj].
pt());
938 if (
isForward(toc[*kj].
eta())) v->getMEhisto_PtForward_HLT()->Fill(toc[*kj].
pt());
939 v->getMEhisto_Eta_HLT()->Fill(toc[*kj].
eta());
940 v->getMEhisto_Phi_HLT()->Fill(toc[*kj].
phi());
941 v->getMEhisto_EtaPhi_HLT()->Fill(toc[*kj].
eta(),toc[*kj].
phi());
946 && (v->getPath().compare(
"PFJet") == 0)){
949 for(
int iCalo=0; iCalo<2; iCalo++){
952 if(v->getTriggerType().compare(
"SingleJet_Trigger") == 0){
953 v->getMEhisto_Pt()->Fill(
CaloJetPt[iCalo]);
962 v->getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*kj].
pt(),
CaloJetPt[iCalo]);
963 v->getMEhisto_EtaCorrelation_HLTRecObj()->Fill(toc[*kj].
eta(),
CaloJetEta[iCalo]);
964 v->getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*kj].
phi(),
CaloJetPhi[iCalo]);
966 v->getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*kj].
pt()-
CaloJetPt[iCalo])/(toc[*kj].
pt()));
967 v->getMEhisto_EtaResolution_HLTRecObj()->Fill((toc[*kj].
eta()-
CaloJetEta[iCalo])/(toc[*kj].
eta()));
968 v->getMEhisto_PhiResolution_HLTRecObj()->Fill((toc[*kj].
phi()-
CaloJetPhi[iCalo])/(toc[*kj].
phi()));
972 if((v->getTriggerType().compare(
"DiJet_Trigger") == 0)){
979 hltPhiVec.push_back(toc[*kj].
phi());
980 hltPtVec.push_back(toc[*kj].
pt());
981 hltEtaVec.push_back(toc[*kj].
eta());
982 hltPxVec.push_back(toc[*kj].px());
983 hltPyVec.push_back(toc[*kj].py());
992 && (v->getPath().compare(
"PFJet") != 0)){
995 for(
int iPF=0; iPF<2; iPF++){
998 if(v->getTriggerType().compare(
"SingleJet_Trigger") == 0){
999 v->getMEhisto_Pt()->Fill(
PFJetPt[iPF]);
1004 v->getMEhisto_Eta()->Fill(
PFJetEta[iPF]);
1005 v->getMEhisto_Phi()->Fill(
PFJetPhi[iPF]);
1008 v->getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*kj].
pt(),
PFJetPt[iPF]);
1009 v->getMEhisto_EtaCorrelation_HLTRecObj()->Fill(toc[*kj].
eta(),
PFJetEta[iPF]);
1010 v->getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*kj].
phi(),
PFJetPhi[iPF]);
1012 v->getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*kj].
pt()-
PFJetPt[iPF])/(toc[*kj].
pt()));
1013 v->getMEhisto_EtaResolution_HLTRecObj()->Fill((toc[*kj].
eta()-
PFJetEta[iPF])/(toc[*kj].
eta()));
1014 v->getMEhisto_PhiResolution_HLTRecObj()->Fill((toc[*kj].
phi()-
PFJetPhi[iPF])/(toc[*kj].
phi()));
1018 if((v->getTriggerType().compare(
"DiJet_Trigger") == 0)){
1019 jetPhiVec.push_back(
PFJetPhi[iPF]);
1020 jetPtVec.push_back(
PFJetPt[iPF]);
1021 jetEtaVec.push_back(
PFJetEta[iPF]);
1022 jetPxVec.push_back(
PFJetPx[iPF]);
1023 jetPyVec.push_back(
PFJetPy[iPF]);
1025 hltPhiVec.push_back(toc[*kj].
phi());
1026 hltPtVec.push_back(toc[*kj].
pt());
1027 hltEtaVec.push_back(toc[*kj].
eta());
1028 hltPxVec.push_back(toc[*kj].px());
1029 hltPyVec.push_back(toc[*kj].py());
1041 && (v->getPath().find(
"HLT_PFMET")==std::string::npos)){
1045 v->getMEhisto_Pt()->Fill(met.
et());
1046 v->getMEhisto_Phi()->Fill(met.
phi());
1048 v->getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*kj].et(),met.
et());
1049 v->getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*kj].
phi(),met.
phi());
1050 v->getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*kj].et()-met.
et())/(toc[*kj].et()));
1051 v->getMEhisto_PhiResolution_HLTRecObj()->Fill((toc[*kj].
phi()-met.
phi())/(toc[*kj].
phi()));
1057 && (v->getPath().find(
"HLT_PFMET")!=std::string::npos)){
1059 const PFMET pfmet = pfmetcol->front();
1061 v->getMEhisto_Pt()->Fill(pfmet.
et());
1062 v->getMEhisto_Phi()->Fill(pfmet.
phi());
1064 v->getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*kj].et(),pfmet.
et());
1065 v->getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*kj].
phi(),pfmet.
phi());
1066 v->getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*kj].et()-pfmet.
et())/(toc[*kj].et()));
1067 v->getMEhisto_PhiResolution_HLTRecObj()->Fill((toc[*kj].
phi()-pfmet.
phi())/(toc[*kj].
phi()));
1070 if((v->getTriggerType().compare(
"DiJet_Trigger") == 0)) diJetFire =
true;
1074 v->getMEhisto_N()->Fill(jetsize);
1077 if((v->getTriggerType().compare(
"DiJet_Trigger") == 0) && jetPtVec.size() >1){
1078 double AveJetPt = (jetPtVec[0] + jetPtVec[1])/2;
1079 double AveJetEta = (jetEtaVec[0] + jetEtaVec[1])/2;
1080 double JetDelPhi =
deltaPhi(jetPhiVec[0],jetPhiVec[1]);
1081 double AveHLTPt = (hltPtVec[0] + hltPtVec[1])/2;
1082 double AveHLTEta = (hltEtaVec[0] + hltEtaVec[1])/2;
1083 double HLTDelPhi =
deltaPhi(hltPhiVec[0],hltPhiVec[1]);
1084 v->getMEhisto_AveragePt_RecObj()->Fill(AveJetPt);
1085 v->getMEhisto_AverageEta_RecObj()->Fill(AveJetEta);
1086 v->getMEhisto_DeltaPhi_RecObj()->Fill(JetDelPhi);
1088 v->getMEhisto_AveragePt_HLTObj()->Fill(AveHLTPt);
1089 v->getMEhisto_AverageEta_HLTObj()->Fill(AveHLTEta);
1090 v->getMEhisto_DeltaPhi_HLTObj()->Fill(HLTDelPhi);
1094 if(v->getPath().find(
"L1") != std::string::npos && !fillL1HLT){
1096 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index "<< l1Index <<
" of that name "<<
l1Tag;
1101 for( trigger::Keys::const_iterator ki = kl1.begin(); ki != kl1.end(); ++ki){
1102 double l1TrigEta = toc[*ki].eta();
1103 double l1TrigPhi = toc[*ki].phi();
1108 if(
deltaR(l1TrigEta, l1TrigPhi, jetEta, jetPhi) < 0.4){
1110 v->getMEhisto_Pt()->Fill(
jet->pt());
1111 if (
isBarrel(
jet->eta())) v->getMEhisto_PtBarrel()->Fill(
jet->pt());
1112 if (
isEndCap(
jet->eta())) v->getMEhisto_PtEndcap()->Fill(
jet->pt());
1113 if (
isForward(
jet->eta())) v->getMEhisto_PtForward()->Fill(
jet->pt());
1115 v->getMEhisto_Eta()->Fill(
jet->eta());
1116 v->getMEhisto_Phi()->Fill(
jet->phi());
1117 v->getMEhisto_EtaPhi()->Fill(
jet->eta(),
jet->phi());
1119 v->getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*ki].
pt(),
jet->pt());
1120 v->getMEhisto_EtaCorrelation_HLTRecObj()->Fill(toc[*ki].
eta(),
jet->eta());
1121 v->getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*ki].
phi(),
jet->phi());
1123 v->getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*ki].
pt()-
jet->pt())/(toc[*ki].
pt()));
1124 v->getMEhisto_EtaResolution_HLTRecObj()->Fill((toc[*ki].
eta()-
jet->eta())/(toc[*ki].
eta()));
1125 v->getMEhisto_PhiResolution_HLTRecObj()->Fill((toc[*ki].
phi()-
jet->phi())/(toc[*ki].
phi()));
1130 v->getMEhisto_N()->Fill(jetsize);
1136 v->getMEhisto_Pt()->Fill(met.
pt());
1137 v->getMEhisto_Phi()->Fill(met.
phi());
1139 v->getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*ki].
pt(),met.
pt());
1140 v->getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*ki].
phi(),met.
phi());
1141 v->getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*ki].
pt()-met.
pt())/(toc[*ki].
pt()));
1142 v->getMEhisto_PhiResolution_HLTRecObj()->Fill((toc[*ki].
phi()-met.
phi())/(toc[*ki].
phi()));
1160 bool denompassed =
false;
1161 bool numpassed =
false;
1167 denompassed =
false;
1176 if(denompassed==
false)
continue;
1187 double trigLowpTcut = 0;
1188 double trigMedpTcut = 0;
1189 double trigHighpTcut = 0;
1190 double trigLowpTcutFwd = 0;
1191 double trigMedpTcutFwd = 0;
1192 double trigHighpTcutFwd = 0;
1195 double trigLowpTPFcut = 0;
1196 double trigMedpTPFcut = 0;
1197 double trigHighpTPFcut = 0;
1198 double trigLowpTPFcutFwd = 0;
1199 double trigMedpTPFcutFwd = 0;
1200 double trigHighpTPFcutFwd = 0;
1203 size_t jetstrfound =
v->getPath().find(
"Jet");
1205 string tpath =
v->getPath();
1209 if(jetstrfound != string::npos){
1211 for(
int trig =
int(jetstrfound)+3; trig < int(jetstrfound)+7; trig++){
1212 if(!isdigit(tpath[trig]))
break;
1213 jetTrigVal+=tpath[trig];
1215 char *cjetTrigVal = (
char*)jetTrigVal.c_str();
1216 jetVal=atof(cjetTrigVal);
1221 trigMedpTcut = jetVal + 5.;
1222 trigHighpTcut = jetVal + 10.;
1224 trigLowpTcutFwd = jetVal + 9.;
1225 trigMedpTcutFwd = jetVal + 15.;
1226 trigHighpTcutFwd = jetVal + 21.;
1230 trigMedpTcut = jetVal + 2.;
1231 trigHighpTcut = jetVal + 60.;
1233 trigLowpTcutFwd = jetVal + 22.;
1234 trigMedpTcutFwd = jetVal + 25.;
1235 trigHighpTcutFwd = jetVal + 110.;
1237 trigLowpTcut = jetVal;
1243 trigMedpTPFcut = jetVal + 20.;
1244 trigHighpTPFcut = jetVal + 40.;
1246 trigLowpTPFcutFwd = jetVal + 60.;
1247 trigMedpTPFcutFwd = jetVal + 80.;
1248 trigHighpTPFcutFwd = jetVal + 100.;
1252 trigMedpTPFcut = jetVal + 40.;
1253 trigHighpTPFcut = jetVal + 140.;
1255 trigLowpTPFcutFwd = jetVal + 110.;
1256 trigMedpTPFcutFwd = jetVal + 130.;
1257 trigHighpTPFcutFwd = jetVal + 190.;
1259 trigLowpTPFcut = jetVal;
1265 if(
verbose_)
std::cout <<
"fillMEforEffAllTrigger: CaloJet -------------------" << std::endl;
1270 bool jetIDbool =
false;
1277 if((
v->getTriggerType().compare(
"SingleJet_Trigger") == 0) &&
calojet.size()){
1279 if(
verbose_)
cout<<
"Passed CaloJet ID -------------------" << endl;
1282 v->getMEhisto_DenominatorPt()->Fill(leadjpt);
1283 if (
isBarrel(leadjeta))
v->getMEhisto_DenominatorPtBarrel()->Fill(leadjpt);
1284 if (
isEndCap(leadjeta))
v->getMEhisto_DenominatorPtEndcap()->Fill(leadjpt);
1285 if (
isForward(leadjeta))
v->getMEhisto_DenominatorPtForward()->Fill(leadjpt);
1286 v->getMEhisto_DenominatorEta()->Fill(leadjeta);
1287 v->getMEhisto_DenominatorPhi()->Fill(leadjphi);
1288 v->getMEhisto_DenominatorEtaPhi()->Fill(leadjeta,leadjphi);
1290 v->getMEhisto_DenominatorEtaBarrel()->Fill(leadjeta);
1291 v->getMEhisto_DenominatorPhiBarrel()->Fill(leadjphi);
1294 v->getMEhisto_DenominatorEtaEndcap()->Fill(leadjeta);
1295 v->getMEhisto_DenominatorPhiEndcap()->Fill(leadjphi);
1298 v->getMEhisto_DenominatorEtaForward()->Fill(leadjeta);
1299 v->getMEhisto_DenominatorPhiForward()->Fill(leadjphi);
1301 if((leadjpt > trigLowpTcut && !
isForward(leadjeta)) || (leadjpt > trigLowpTcutFwd &&
isForward(leadjeta))){
1302 v->getMEhisto_DenominatorEta_LowpTcut()->Fill(leadjeta);
1303 v->getMEhisto_DenominatorPhi_LowpTcut()->Fill(leadjphi);
1304 v->getMEhisto_DenominatorEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
1306 if((leadjpt > trigMedpTcut && !
isForward(leadjeta)) || (leadjpt > trigMedpTcutFwd &&
isForward(leadjeta))){
1307 v->getMEhisto_DenominatorEta_MedpTcut()->Fill(leadjeta);
1308 v->getMEhisto_DenominatorPhi_MedpTcut()->Fill(leadjphi);
1309 v->getMEhisto_DenominatorEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
1311 if((leadjpt > trigHighpTcut && !
isForward(leadjeta)) || (leadjpt > trigHighpTcutFwd &&
isForward(leadjeta))){
1312 v->getMEhisto_DenominatorEta_HighpTcut()->Fill(leadjeta);
1313 v->getMEhisto_DenominatorPhi_HighpTcut()->Fill(leadjphi);
1314 v->getMEhisto_DenominatorEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
1320 double dRmin = 99999.;
1321 double dPhimin = 9999.;
1322 if(
v->getPath().find(
"L1") != std::string::npos){
1324 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
1328 for(trigger::Keys::const_iterator ki = kl1.begin();ki != kl1.end(); ++ki){
1329 double dR =
deltaR(toc[*ki].
eta(), toc[*ki].
phi(), leadjeta, leadjphi);
1338 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
1342 trigger::Keys::const_iterator kj = khlt.begin();
1343 for(;kj != khlt.end(); ++kj){
1345 leadjeta, leadjphi);
1355 v->getMEhisto_DeltaPhi()->Fill(dPhimin);
1356 v->getMEhisto_DeltaR()->Fill(dRmin);
1359 if(dRmin < 0.1 || (
v->getPath().find(
"L1") != std::string::npos && dRmin < 0.4)){
1360 v->getMEhisto_NumeratorPt()->Fill(leadjpt);
1361 if (
isBarrel(leadjeta))
v->getMEhisto_NumeratorPtBarrel()->Fill(leadjpt);
1362 if (
isEndCap(leadjeta))
v->getMEhisto_NumeratorPtEndcap()->Fill(leadjpt);
1363 if (
isForward(leadjeta))
v->getMEhisto_NumeratorPtForward()->Fill(leadjpt);
1364 v->getMEhisto_NumeratorEta()->Fill(leadjeta);
1365 v->getMEhisto_NumeratorPhi()->Fill(leadjphi);
1366 v->getMEhisto_NumeratorEtaPhi()->Fill(leadjeta,leadjphi);
1368 v->getMEhisto_NumeratorEtaBarrel()->Fill(leadjeta);
1369 v->getMEhisto_NumeratorPhiBarrel()->Fill(leadjphi);
1372 v->getMEhisto_NumeratorEtaEndcap()->Fill(leadjeta);
1373 v->getMEhisto_NumeratorPhiEndcap()->Fill(leadjphi);
1376 v->getMEhisto_NumeratorEtaForward()->Fill(leadjeta);
1377 v->getMEhisto_NumeratorPhiForward()->Fill(leadjphi);
1379 if((leadjpt > trigLowpTcut && !
isForward(leadjeta)) || (leadjpt > trigLowpTcutFwd &&
isForward(leadjeta))){
1380 v->getMEhisto_NumeratorEta_LowpTcut()->Fill(leadjeta);
1381 v->getMEhisto_NumeratorPhi_LowpTcut()->Fill(leadjphi);
1382 v->getMEhisto_NumeratorEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
1384 if((leadjpt > trigMedpTcut && !
isForward(leadjeta)) || (leadjpt > trigMedpTcutFwd &&
isForward(leadjeta))){
1385 v->getMEhisto_NumeratorEta_MedpTcut()->Fill(leadjeta);
1386 v->getMEhisto_NumeratorPhi_MedpTcut()->Fill(leadjphi);
1387 v->getMEhisto_NumeratorEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
1389 if((leadjpt > trigHighpTcut && !
isForward(leadjeta)) || (leadjpt > trigHighpTcutFwd &&
isForward(leadjeta))){
1390 v->getMEhisto_NumeratorEta_HighpTcut()->Fill(leadjeta);
1391 v->getMEhisto_NumeratorPhi_HighpTcut()->Fill(leadjphi);
1392 v->getMEhisto_NumeratorEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
1399 if(jetIDbool ==
true && (
v->getTriggerType().compare(
"DiJet_Trigger") == 0) &&
calojet.size()>1){
1404 if(numpassed==
true){
1413 if(
verbose_)
std::cout <<
"fillMEforEffAllTrigger: PFJet -------------------" << std::endl;
1418 bool jetIDbool =
false;
1437 if((
v->getTriggerType().compare(
"SingleJet_Trigger") == 0) &&
pfjet.size()){
1440 _pfMHT =
sqrt(pfMHTx*pfMHTx + pfMHTy*pfMHTy);
1441 v->getMEhisto_DenominatorPFMHT()->Fill(
_pfMHT);
1447 if(
verbose_)
cout<<
"Passed PFJet ID -------------------" << endl;
1449 v->getMEhisto_DenominatorPFPt()->Fill(leadjpt);
1450 if (
isBarrel(leadjeta))
v->getMEhisto_DenominatorPFPtBarrel()->Fill(leadjpt);
1451 if (
isEndCap(leadjeta))
v->getMEhisto_DenominatorPFPtEndcap()->Fill(leadjpt);
1452 if (
isForward(leadjeta))
v->getMEhisto_DenominatorPFPtForward()->Fill(leadjpt);
1453 v->getMEhisto_DenominatorPFEta()->Fill(leadjeta);
1454 v->getMEhisto_DenominatorPFPhi()->Fill(leadjphi);
1455 v->getMEhisto_DenominatorPFEtaPhi()->Fill(leadjeta,leadjphi);
1457 v->getMEhisto_DenominatorPFEtaBarrel()->Fill(leadjeta);
1458 v->getMEhisto_DenominatorPFPhiBarrel()->Fill(leadjphi);
1461 v->getMEhisto_DenominatorPFEtaEndcap()->Fill(leadjeta);
1462 v->getMEhisto_DenominatorPFPhiEndcap()->Fill(leadjphi);
1465 v->getMEhisto_DenominatorPFEtaForward()->Fill(leadjeta);
1466 v->getMEhisto_DenominatorPFPhiForward()->Fill(leadjphi);
1468 if((leadjpt > trigLowpTPFcut && !
isForward(leadjeta)) || (leadjpt > trigLowpTPFcutFwd &&
isForward(leadjeta))){
1469 v->getMEhisto_DenominatorPFEta_LowpTcut()->Fill(leadjeta);
1470 v->getMEhisto_DenominatorPFPhi_LowpTcut()->Fill(leadjphi);
1471 v->getMEhisto_DenominatorPFEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
1473 if((leadjpt > trigMedpTPFcut && !
isForward(leadjeta)) || (leadjpt > trigMedpTPFcutFwd &&
isForward(leadjeta))){
1474 v->getMEhisto_DenominatorPFEta_MedpTcut()->Fill(leadjeta);
1475 v->getMEhisto_DenominatorPFPhi_MedpTcut()->Fill(leadjphi);
1476 v->getMEhisto_DenominatorPFEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
1478 if((leadjpt > trigHighpTPFcut && !
isForward(leadjeta)) || (leadjpt > trigHighpTPFcutFwd &&
isForward(leadjeta))){
1479 v->getMEhisto_DenominatorPFEta_HighpTcut()->Fill(leadjeta);
1480 v->getMEhisto_DenominatorPFPhi_HighpTcut()->Fill(leadjphi);
1481 v->getMEhisto_DenominatorPFEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
1486 double dRmin = 99999.;
1487 double dPhimin = 9999.;
1488 if(
v->getPath().find(
"L1") != std::string::npos){
1490 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
1494 for(trigger::Keys::const_iterator ki = kl1.begin();ki != kl1.end(); ++ki){
1495 double dR =
deltaR(toc[*ki].
eta(), toc[*ki].
phi(), leadjeta, leadjphi);
1504 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
1508 for(trigger::Keys::const_iterator kj = khlt.begin();kj != khlt.end(); ++kj){
1509 double dR =
deltaR(toc[*kj].
eta(), toc[*kj].
phi(), leadjeta, leadjphi);
1518 v->getMEhisto_PFDeltaPhi()->Fill(dPhimin);
1519 v->getMEhisto_PFDeltaR()->Fill(dRmin);
1522 if(dRmin < 0.1 || (
v->getPath().find(
"L1") != std::string::npos && dRmin < 0.4)){
1523 v->getMEhisto_NumeratorPFPt()->Fill(leadjpt);
1524 if (
isBarrel(leadjeta))
v->getMEhisto_NumeratorPFPtBarrel()->Fill(leadjpt);
1525 if (
isEndCap(leadjeta))
v->getMEhisto_NumeratorPFPtEndcap()->Fill(leadjpt);
1526 if (
isForward(leadjeta))
v->getMEhisto_NumeratorPFPtForward()->Fill(leadjpt);
1527 v->getMEhisto_NumeratorPFEta()->Fill(leadjeta);
1528 v->getMEhisto_NumeratorPFPhi()->Fill(leadjphi);
1529 v->getMEhisto_NumeratorPFEtaPhi()->Fill(leadjeta,leadjphi);
1531 v->getMEhisto_NumeratorPFEtaBarrel()->Fill(leadjeta);
1532 v->getMEhisto_NumeratorPFPhiBarrel()->Fill(leadjphi);
1535 v->getMEhisto_NumeratorPFEtaEndcap()->Fill(leadjeta);
1536 v->getMEhisto_NumeratorPFPhiEndcap()->Fill(leadjphi);
1539 v->getMEhisto_NumeratorPFEtaForward()->Fill(leadjeta);
1540 v->getMEhisto_NumeratorPFPhiForward()->Fill(leadjphi);
1542 if((leadjpt > trigLowpTPFcut && !
isForward(leadjeta))
1543 || (leadjpt > trigLowpTPFcutFwd &&
isForward(leadjeta))){
1544 v->getMEhisto_NumeratorPFEta_LowpTcut()->Fill(leadjeta);
1545 v->getMEhisto_NumeratorPFPhi_LowpTcut()->Fill(leadjphi);
1546 v->getMEhisto_NumeratorPFEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
1548 if((leadjpt > trigMedpTPFcut && !
isForward(leadjeta))
1549 || (leadjpt > trigMedpTPFcutFwd &&
isForward(leadjeta))){
1550 v->getMEhisto_NumeratorPFEta_MedpTcut()->Fill(leadjeta);
1551 v->getMEhisto_NumeratorPFPhi_MedpTcut()->Fill(leadjphi);
1552 v->getMEhisto_NumeratorPFEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
1554 if((leadjpt > trigHighpTPFcut && !
isForward(leadjeta))
1555 || (leadjpt > trigHighpTPFcutFwd &&
isForward(leadjeta))){
1556 v->getMEhisto_NumeratorPFEta_HighpTcut()->Fill(leadjeta);
1557 v->getMEhisto_NumeratorPFPhi_HighpTcut()->Fill(leadjphi);
1558 v->getMEhisto_NumeratorPFEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
1564 if(jetIDbool ==
true && (
v->getTriggerType().compare(
"DiJet_Trigger") == 0) &&
pfjet.size()>1){
1584 if(
verbose_)
std::cout <<
"fillMEforEffAllTrigger: CaloMET -------------------" << std::endl;
1587 && (
v->getPath().find(
"HLT_PFMET")==std::string::npos)){
1590 v->getMEhisto_DenominatorPt()->Fill(met.
et());
1591 v->getMEhisto_DenominatorPhi()->Fill(met.
phi());
1593 v->getMEhisto_NumeratorPt()->Fill(met.
et());
1594 v->getMEhisto_NumeratorPhi()->Fill(met.
phi());
1596 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
1599 double dPhimin = 9999.;
1601 for(trigger::Keys::const_iterator kj = khlt.begin();kj != khlt.end(); ++kj){
1607 v->getMEhisto_DeltaPhi()->Fill(dPhimin);
1613 if(
verbose_)
std::cout <<
"fillMEforEffAllTrigger: PFMET -------------------" << std::endl;
1616 && (
v->getPath().find(
"HLT_PFMET")!=std::string::npos)){
1618 const PFMET met = pfmetcol->front();
1619 v->getMEhisto_DenominatorPt()->Fill(met.
et());
1620 v->getMEhisto_DenominatorPhi()->Fill(met.
phi());
1622 v->getMEhisto_NumeratorPt()->Fill(met.
et());
1623 v->getMEhisto_NumeratorPhi()->Fill(met.
phi());
1625 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
1628 double dPhimin = 9999.;
1630 for(trigger::Keys::const_iterator kj = khlt.begin();kj != khlt.end(); ++kj){
1636 v->getMEhisto_DeltaPhi()->Fill(dPhimin);
1659 bool muTrig =
false;
1660 bool numpassed =
false;
1664 for (
unsigned int j=0;
j!=nPath; ++
j) {
1669 if(
verbose_)
cout<<
"fillMEforEffWrtMuTrigger: Muon Match"<<endl;
1694 double trigLowpTcut = 0;
1695 double trigMedpTcut = 0;
1696 double trigHighpTcut = 0;
1697 double trigLowpTcutFwd = 0;
1698 double trigMedpTcutFwd = 0;
1699 double trigHighpTcutFwd = 0;
1702 double trigLowpTPFcut = 0;
1703 double trigMedpTPFcut = 0;
1704 double trigHighpTPFcut = 0;
1705 double trigLowpTPFcutFwd = 0;
1706 double trigMedpTPFcutFwd = 0;
1707 double trigHighpTPFcutFwd = 0;
1710 size_t jetstrfound = v->getPath().find(
"Jet");
1712 string tpath = v->getPath();
1716 if(jetstrfound != string::npos){
1718 for(
int trig =
int(jetstrfound)+3; trig < int(jetstrfound)+7; trig++){
1719 if(!isdigit(tpath[trig]))
break;
1720 jetTrigVal+=tpath[trig];
1722 char *cjetTrigVal = (
char*)jetTrigVal.c_str();
1723 jetVal=atof(cjetTrigVal);
1729 trigMedpTcut = jetVal + 5.;
1730 trigHighpTcut = jetVal + 10.;
1732 trigLowpTcutFwd = jetVal + 9.;
1733 trigMedpTcutFwd = jetVal + 15.;
1734 trigHighpTcutFwd = jetVal + 21.;
1738 trigMedpTcut = jetVal + 2.;
1739 trigHighpTcut = jetVal + 60.;
1741 trigLowpTcutFwd = jetVal + 22.;
1742 trigMedpTcutFwd = jetVal + 25.;
1743 trigHighpTcutFwd = jetVal + 110.;
1745 trigLowpTcut = jetVal;
1751 trigMedpTPFcut = jetVal + 20.;
1752 trigHighpTPFcut = jetVal + 40.;
1754 trigLowpTPFcutFwd = jetVal + 60.;
1755 trigMedpTPFcutFwd = jetVal + 80.;
1756 trigHighpTPFcutFwd = jetVal + 100.;
1760 trigMedpTPFcut = jetVal + 40.;
1761 trigHighpTPFcut = jetVal + 140.;
1763 trigLowpTPFcutFwd = jetVal + 110.;
1764 trigMedpTPFcutFwd = jetVal + 130.;
1765 trigHighpTPFcutFwd = jetVal + 190.;
1767 trigLowpTPFcut = jetVal;
1773 if(
verbose_)
std::cout <<
"fillMEforEffWrtMuTrigger: CaloJet -------------------" << std::endl;
1779 bool jetIDbool =
false;
1786 if((v->getTriggerType().compare(
"SingleJet_Trigger") == 0) &&
calojet.size()){
1788 if(
verbose_)
cout<<
"passed CaloJet ID -------------------" << endl;
1792 v->getMEhisto_DenominatorPt()->Fill(leadjpt);
1793 if (
isBarrel(leadjeta)) v->getMEhisto_DenominatorPtBarrel()->Fill(leadjpt);
1794 if (
isEndCap(leadjeta)) v->getMEhisto_DenominatorPtEndcap()->Fill(leadjpt);
1795 if (
isForward(leadjeta)) v->getMEhisto_DenominatorPtForward()->Fill(leadjpt);
1796 v->getMEhisto_DenominatorEta()->Fill(leadjeta);
1797 v->getMEhisto_DenominatorPhi()->Fill(leadjphi);
1798 v->getMEhisto_DenominatorEtaPhi()->Fill(leadjeta,leadjphi);
1800 v->getMEhisto_DenominatorEtaBarrel()->Fill(leadjeta);
1801 v->getMEhisto_DenominatorPhiBarrel()->Fill(leadjphi);
1804 v->getMEhisto_DenominatorEtaEndcap()->Fill(leadjeta);
1805 v->getMEhisto_DenominatorPhiEndcap()->Fill(leadjphi);
1808 v->getMEhisto_DenominatorEtaForward()->Fill(leadjeta);
1809 v->getMEhisto_DenominatorPhiForward()->Fill(leadjphi);
1811 if((leadjpt > trigLowpTcut && !
isForward(leadjeta)) || (leadjpt > trigLowpTcutFwd &&
isForward(leadjeta))){
1812 v->getMEhisto_DenominatorEta_LowpTcut()->Fill(leadjeta);
1813 v->getMEhisto_DenominatorPhi_LowpTcut()->Fill(leadjphi);
1814 v->getMEhisto_DenominatorEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
1816 if((leadjpt > trigMedpTcut && !
isForward(leadjeta)) || (leadjpt > trigMedpTcutFwd &&
isForward(leadjeta))){
1817 v->getMEhisto_DenominatorEta_MedpTcut()->Fill(leadjeta);
1818 v->getMEhisto_DenominatorPhi_MedpTcut()->Fill(leadjphi);
1819 v->getMEhisto_DenominatorEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
1821 if((leadjpt > trigHighpTcut && !
isForward(leadjeta)) || (leadjpt > trigHighpTcutFwd &&
isForward(leadjeta))){
1822 v->getMEhisto_DenominatorEta_HighpTcut()->Fill(leadjeta);
1823 v->getMEhisto_DenominatorPhi_HighpTcut()->Fill(leadjphi);
1824 v->getMEhisto_DenominatorEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
1830 double dRmin = 99999.;
1831 double dPhimin = 9999.;
1832 if(v->getPath().find(
"L1") != std::string::npos){
1834 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
1838 for(trigger::Keys::const_iterator ki = kl1.begin();ki != kl1.end(); ++ki){
1839 double dR =
deltaR(toc[*ki].
eta(), toc[*ki].
phi(), leadjeta, leadjphi);
1848 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
1852 trigger::Keys::const_iterator kj = khlt.begin();
1853 for(;kj != khlt.end(); ++kj){
1855 leadjeta, leadjphi);
1865 v->getMEhisto_DeltaPhi()->Fill(dPhimin);
1866 v->getMEhisto_DeltaR()->Fill(dRmin);
1869 if(dRmin < 0.1 || (v->getPath().find(
"L1") != std::string::npos && dRmin < 0.4)){
1870 v->getMEhisto_NumeratorPt()->Fill(leadjpt);
1871 if (
isBarrel(leadjeta)) v->getMEhisto_NumeratorPtBarrel()->Fill(leadjpt);
1872 if (
isEndCap(leadjeta)) v->getMEhisto_NumeratorPtEndcap()->Fill(leadjpt);
1873 if (
isForward(leadjeta)) v->getMEhisto_NumeratorPtForward()->Fill(leadjpt);
1874 v->getMEhisto_NumeratorEta()->Fill(leadjeta);
1875 v->getMEhisto_NumeratorPhi()->Fill(leadjphi);
1876 v->getMEhisto_NumeratorEtaPhi()->Fill(leadjeta,leadjphi);
1878 v->getMEhisto_NumeratorEtaBarrel()->Fill(leadjeta);
1879 v->getMEhisto_NumeratorPhiBarrel()->Fill(leadjphi);
1882 v->getMEhisto_NumeratorEtaEndcap()->Fill(leadjeta);
1883 v->getMEhisto_NumeratorPhiEndcap()->Fill(leadjphi);
1886 v->getMEhisto_NumeratorEtaForward()->Fill(leadjeta);
1887 v->getMEhisto_NumeratorPhiForward()->Fill(leadjphi);
1889 if((leadjpt > trigLowpTcut && !
isForward(leadjeta)) || (leadjpt > trigLowpTcutFwd &&
isForward(leadjeta))){
1890 v->getMEhisto_NumeratorEta_LowpTcut()->Fill(leadjeta);
1891 v->getMEhisto_NumeratorPhi_LowpTcut()->Fill(leadjphi);
1892 v->getMEhisto_NumeratorEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
1894 if((leadjpt > trigMedpTcut && !
isForward(leadjeta)) || (leadjpt > trigMedpTcutFwd &&
isForward(leadjeta))){
1895 v->getMEhisto_NumeratorEta_MedpTcut()->Fill(leadjeta);
1896 v->getMEhisto_NumeratorPhi_MedpTcut()->Fill(leadjphi);
1897 v->getMEhisto_NumeratorEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
1899 if((leadjpt > trigHighpTcut && !
isForward(leadjeta)) || (leadjpt > trigHighpTcutFwd &&
isForward(leadjeta))){
1900 v->getMEhisto_NumeratorEta_HighpTcut()->Fill(leadjeta);
1901 v->getMEhisto_NumeratorPhi_HighpTcut()->Fill(leadjphi);
1902 v->getMEhisto_NumeratorEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
1909 if(jetIDbool ==
true && (v->getTriggerType().compare(
"DiJet_Trigger") == 0) &&
calojet.size()>1){
1914 if(numpassed==
true){
1923 if(
verbose_)
std::cout <<
"fillMEforEffWrtMuTrigger: PFJet -------------------" << std::endl;
1929 bool jetIDbool =
false;
1948 if((v->getTriggerType().compare(
"SingleJet_Trigger") == 0) &&
pfjet.size()){
1951 _pfMHT =
sqrt(pfMHTx*pfMHTx + pfMHTy*pfMHTy);
1952 v->getMEhisto_DenominatorPFMHT()->Fill(
_pfMHT);
1958 if(
verbose_)
cout<<
"passed PFJet ID -------------------" << endl;
1960 v->getMEhisto_DenominatorPFPt()->Fill(leadjpt);
1961 if (
isBarrel(leadjeta)) v->getMEhisto_DenominatorPFPtBarrel()->Fill(leadjpt);
1962 if (
isEndCap(leadjeta)) v->getMEhisto_DenominatorPFPtEndcap()->Fill(leadjpt);
1963 if (
isForward(leadjeta)) v->getMEhisto_DenominatorPFPtForward()->Fill(leadjpt);
1964 v->getMEhisto_DenominatorPFEta()->Fill(leadjeta);
1965 v->getMEhisto_DenominatorPFPhi()->Fill(leadjphi);
1966 v->getMEhisto_DenominatorPFEtaPhi()->Fill(leadjeta,leadjphi);
1968 v->getMEhisto_DenominatorPFEtaBarrel()->Fill(leadjeta);
1969 v->getMEhisto_DenominatorPFPhiBarrel()->Fill(leadjphi);
1972 v->getMEhisto_DenominatorPFEtaEndcap()->Fill(leadjeta);
1973 v->getMEhisto_DenominatorPFPhiEndcap()->Fill(leadjphi);
1976 v->getMEhisto_DenominatorPFEtaForward()->Fill(leadjeta);
1977 v->getMEhisto_DenominatorPFPhiForward()->Fill(leadjphi);
1979 if((leadjpt > trigLowpTPFcut && !
isForward(leadjeta)) || (leadjpt > trigLowpTPFcutFwd &&
isForward(leadjeta))){
1980 v->getMEhisto_DenominatorPFEta_LowpTcut()->Fill(leadjeta);
1981 v->getMEhisto_DenominatorPFPhi_LowpTcut()->Fill(leadjphi);
1982 v->getMEhisto_DenominatorPFEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
1984 if((leadjpt > trigMedpTPFcut && !
isForward(leadjeta)) || (leadjpt > trigMedpTPFcutFwd &&
isForward(leadjeta))){
1985 v->getMEhisto_DenominatorPFEta_MedpTcut()->Fill(leadjeta);
1986 v->getMEhisto_DenominatorPFPhi_MedpTcut()->Fill(leadjphi);
1987 v->getMEhisto_DenominatorPFEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
1989 if((leadjpt > trigHighpTPFcut && !
isForward(leadjeta)) || (leadjpt > trigHighpTPFcutFwd &&
isForward(leadjeta))){
1990 v->getMEhisto_DenominatorPFEta_HighpTcut()->Fill(leadjeta);
1991 v->getMEhisto_DenominatorPFPhi_HighpTcut()->Fill(leadjphi);
1992 v->getMEhisto_DenominatorPFEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
1997 double dRmin = 99999.;
1998 double dPhimin = 9999.;
1999 if(v->getPath().find(
"L1") != std::string::npos){
2001 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
2005 for(trigger::Keys::const_iterator ki = kl1.begin();ki != kl1.end(); ++ki){
2006 double dR =
deltaR(toc[*ki].
eta(), toc[*ki].
phi(), leadjeta, leadjphi);
2015 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
2019 for(trigger::Keys::const_iterator kj = khlt.begin();kj != khlt.end(); ++kj){
2020 double dR =
deltaR(toc[*kj].
eta(), toc[*kj].
phi(), leadjeta, leadjphi);
2029 v->getMEhisto_PFDeltaPhi()->Fill(dPhimin);
2030 v->getMEhisto_PFDeltaR()->Fill(dRmin);
2033 if(dRmin < 0.1 || (v->getPath().find(
"L1") != std::string::npos && dRmin < 0.4)){
2034 v->getMEhisto_NumeratorPFPt()->Fill(leadjpt);
2035 if (
isBarrel(leadjeta)) v->getMEhisto_NumeratorPFPtBarrel()->Fill(leadjpt);
2036 if (
isEndCap(leadjeta)) v->getMEhisto_NumeratorPFPtEndcap()->Fill(leadjpt);
2037 if (
isForward(leadjeta)) v->getMEhisto_NumeratorPFPtForward()->Fill(leadjpt);
2038 v->getMEhisto_NumeratorPFEta()->Fill(leadjeta);
2039 v->getMEhisto_NumeratorPFPhi()->Fill(leadjphi);
2040 v->getMEhisto_NumeratorPFEtaPhi()->Fill(leadjeta,leadjphi);
2042 v->getMEhisto_NumeratorPFEtaBarrel()->Fill(leadjeta);
2043 v->getMEhisto_NumeratorPFPhiBarrel()->Fill(leadjphi);
2046 v->getMEhisto_NumeratorPFEtaEndcap()->Fill(leadjeta);
2047 v->getMEhisto_NumeratorPFPhiEndcap()->Fill(leadjphi);
2050 v->getMEhisto_NumeratorPFEtaForward()->Fill(leadjeta);
2051 v->getMEhisto_NumeratorPFPhiForward()->Fill(leadjphi);
2053 if((leadjpt > trigLowpTPFcut && !
isForward(leadjeta))
2054 || (leadjpt > trigLowpTPFcutFwd &&
isForward(leadjeta))){
2055 v->getMEhisto_NumeratorPFEta_LowpTcut()->Fill(leadjeta);
2056 v->getMEhisto_NumeratorPFPhi_LowpTcut()->Fill(leadjphi);
2057 v->getMEhisto_NumeratorPFEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
2059 if((leadjpt > trigMedpTPFcut && !
isForward(leadjeta))
2060 || (leadjpt > trigMedpTPFcutFwd &&
isForward(leadjeta))){
2061 v->getMEhisto_NumeratorPFEta_MedpTcut()->Fill(leadjeta);
2062 v->getMEhisto_NumeratorPFPhi_MedpTcut()->Fill(leadjphi);
2063 v->getMEhisto_NumeratorPFEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
2065 if((leadjpt > trigHighpTPFcut && !
isForward(leadjeta))
2066 || (leadjpt > trigHighpTPFcutFwd &&
isForward(leadjeta))){
2067 v->getMEhisto_NumeratorPFEta_HighpTcut()->Fill(leadjeta);
2068 v->getMEhisto_NumeratorPFPhi_HighpTcut()->Fill(leadjphi);
2069 v->getMEhisto_NumeratorPFEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
2075 if(jetIDbool ==
true && (v->getTriggerType().compare(
"DiJet_Trigger") == 0) &&
pfjet.size()>1){
2084 v->getMEhisto_DenominatorPFPt()->Fill((
PFJetPt[0] +
PFJetPt[1])/2.);
2087 v->getMEhisto_NumeratorPFPt()->Fill((
PFJetPt[0] +
PFJetPt[1])/2.);
2095 if(
verbose_)
std::cout <<
"fillMEforEffWrtMuTrigger: CaloMET -------------------" << std::endl;
2098 && (v->getPath().find(
"HLT_PFMET")==std::string::npos)){
2102 v->getMEhisto_DenominatorPt()->Fill(met.
et());
2103 v->getMEhisto_DenominatorPhi()->Fill(met.
phi());
2105 v->getMEhisto_NumeratorPt()->Fill(met.
et());
2106 v->getMEhisto_NumeratorPhi()->Fill(met.
phi());
2108 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
2111 double dPhimin = 9999.;
2113 for(trigger::Keys::const_iterator kj = khlt.begin();kj != khlt.end(); ++kj){
2119 v->getMEhisto_DeltaPhi()->Fill(dPhimin);
2125 if(
verbose_)
std::cout <<
"fillMEforEffWrtMuTrigger: PFMET -------------------" << std::endl;
2128 && (v->getPath().find(
"HLT_PFMET")!=std::string::npos)){
2131 const PFMET met = pfmetcol->front();
2132 v->getMEhisto_DenominatorPt()->Fill(met.
et());
2133 v->getMEhisto_DenominatorPhi()->Fill(met.
phi());
2135 v->getMEhisto_NumeratorPt()->Fill(met.
et());
2136 v->getMEhisto_NumeratorPhi()->Fill(met.
phi());
2138 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
2141 double dPhimin = 9999.;
2143 for(trigger::Keys::const_iterator kj = khlt.begin();kj != khlt.end(); ++kj){
2149 v->getMEhisto_DeltaPhi()->Fill(dPhimin);
2173 bool mbTrig =
false;
2174 bool numpassed =
false;
2178 for (
unsigned int j=0;
j!=nPath; ++
j) {
2183 if(
verbose_)
cout<<
"fillMEforEffWrtMBTrigger: MinBias Match"<<endl;
2208 double trigLowpTcut = 0;
2209 double trigMedpTcut = 0;
2210 double trigHighpTcut = 0;
2211 double trigLowpTcutFwd = 0;
2212 double trigMedpTcutFwd = 0;
2213 double trigHighpTcutFwd = 0;
2216 double trigLowpTPFcut = 0;
2217 double trigMedpTPFcut = 0;
2218 double trigHighpTPFcut = 0;
2219 double trigLowpTPFcutFwd = 0;
2220 double trigMedpTPFcutFwd = 0;
2221 double trigHighpTPFcutFwd = 0;
2224 size_t jetstrfound = v->getPath().find(
"Jet");
2226 string tpath = v->getPath();
2230 if(jetstrfound != string::npos){
2232 for(
int trig =
int(jetstrfound)+3; trig < int(jetstrfound)+7; trig++){
2233 if(!isdigit(tpath[trig]))
break;
2234 jetTrigVal+=tpath[trig];
2236 char *cjetTrigVal = (
char*)jetTrigVal.c_str();
2237 jetVal=atof(cjetTrigVal);
2243 trigMedpTcut = jetVal + 5.;
2244 trigHighpTcut = jetVal + 10.;
2246 trigLowpTcutFwd = jetVal + 9.;
2247 trigMedpTcutFwd = jetVal + 15.;
2248 trigHighpTcutFwd = jetVal + 21.;
2252 trigMedpTcut = jetVal + 2.;
2253 trigHighpTcut = jetVal + 60.;
2255 trigLowpTcutFwd = jetVal + 22.;
2256 trigMedpTcutFwd = jetVal + 25.;
2257 trigHighpTcutFwd = jetVal + 110.;
2259 trigLowpTcut = jetVal;
2266 trigMedpTPFcut = jetVal + 20.;
2267 trigHighpTPFcut = jetVal + 40.;
2269 trigLowpTPFcutFwd = jetVal + 60.;
2270 trigMedpTPFcutFwd = jetVal + 80.;
2271 trigHighpTPFcutFwd = jetVal + 100.;
2275 trigMedpTPFcut = jetVal + 40.;
2276 trigHighpTPFcut = jetVal + 140.;
2278 trigLowpTPFcutFwd = jetVal + 110.;
2279 trigMedpTPFcutFwd = jetVal + 130.;
2280 trigHighpTPFcutFwd = jetVal + 190.;
2282 trigLowpTPFcut = jetVal;
2289 if(
verbose_)
std::cout <<
"fillMEforEffWrtMuTrigger: CaloJet -------------------" << std::endl;
2295 bool jetIDbool =
false;
2302 if((v->getTriggerType().compare(
"SingleJet_Trigger") == 0) &&
calojet.size()){
2304 if(
verbose_)
cout<<
"passed CaloJet ID -------------------" << endl;
2308 v->getMEhisto_DenominatorPt()->Fill(leadjpt);
2309 if (
isBarrel(leadjeta)) v->getMEhisto_DenominatorPtBarrel()->Fill(leadjpt);
2310 if (
isEndCap(leadjeta)) v->getMEhisto_DenominatorPtEndcap()->Fill(leadjpt);
2311 if (
isForward(leadjeta)) v->getMEhisto_DenominatorPtForward()->Fill(leadjpt);
2312 v->getMEhisto_DenominatorEta()->Fill(leadjeta);
2313 v->getMEhisto_DenominatorPhi()->Fill(leadjphi);
2314 v->getMEhisto_DenominatorEtaPhi()->Fill(leadjeta,leadjphi);
2316 v->getMEhisto_DenominatorEtaBarrel()->Fill(leadjeta);
2317 v->getMEhisto_DenominatorPhiBarrel()->Fill(leadjphi);
2320 v->getMEhisto_DenominatorEtaEndcap()->Fill(leadjeta);
2321 v->getMEhisto_DenominatorPhiEndcap()->Fill(leadjphi);
2324 v->getMEhisto_DenominatorEtaForward()->Fill(leadjeta);
2325 v->getMEhisto_DenominatorPhiForward()->Fill(leadjphi);
2327 if((leadjpt > trigLowpTcut && !
isForward(leadjeta)) || (leadjpt > trigLowpTcutFwd &&
isForward(leadjeta))){
2328 v->getMEhisto_DenominatorEta_LowpTcut()->Fill(leadjeta);
2329 v->getMEhisto_DenominatorPhi_LowpTcut()->Fill(leadjphi);
2330 v->getMEhisto_DenominatorEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
2332 if((leadjpt > trigMedpTcut && !
isForward(leadjeta)) || (leadjpt > trigMedpTcutFwd &&
isForward(leadjeta))){
2333 v->getMEhisto_DenominatorEta_MedpTcut()->Fill(leadjeta);
2334 v->getMEhisto_DenominatorPhi_MedpTcut()->Fill(leadjphi);
2335 v->getMEhisto_DenominatorEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
2337 if((leadjpt > trigHighpTcut && !
isForward(leadjeta)) || (leadjpt > trigHighpTcutFwd &&
isForward(leadjeta))){
2338 v->getMEhisto_DenominatorEta_HighpTcut()->Fill(leadjeta);
2339 v->getMEhisto_DenominatorPhi_HighpTcut()->Fill(leadjphi);
2340 v->getMEhisto_DenominatorEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
2346 double dRmin = 99999.;
2347 double dPhimin = 9999.;
2348 if(v->getPath().find(
"L1") != std::string::npos){
2350 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
2354 for(trigger::Keys::const_iterator ki = kl1.begin();ki != kl1.end(); ++ki){
2355 double dR =
deltaR(toc[*ki].
eta(), toc[*ki].
phi(), leadjeta, leadjphi);
2364 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
2368 trigger::Keys::const_iterator kj = khlt.begin();
2369 for(;kj != khlt.end(); ++kj){
2371 leadjeta, leadjphi);
2381 v->getMEhisto_DeltaPhi()->Fill(dPhimin);
2382 v->getMEhisto_DeltaR()->Fill(dRmin);
2385 if(dRmin < 0.1 || (v->getPath().find(
"L1") != std::string::npos && dRmin < 0.4)){
2386 v->getMEhisto_NumeratorPt()->Fill(leadjpt);
2387 if (
isBarrel(leadjeta)) v->getMEhisto_NumeratorPtBarrel()->Fill(leadjpt);
2388 if (
isEndCap(leadjeta)) v->getMEhisto_NumeratorPtEndcap()->Fill(leadjpt);
2389 if (
isForward(leadjeta)) v->getMEhisto_NumeratorPtForward()->Fill(leadjpt);
2390 v->getMEhisto_NumeratorEta()->Fill(leadjeta);
2391 v->getMEhisto_NumeratorPhi()->Fill(leadjphi);
2392 v->getMEhisto_NumeratorEtaPhi()->Fill(leadjeta,leadjphi);
2394 v->getMEhisto_NumeratorEtaBarrel()->Fill(leadjeta);
2395 v->getMEhisto_NumeratorPhiBarrel()->Fill(leadjphi);
2398 v->getMEhisto_NumeratorEtaEndcap()->Fill(leadjeta);
2399 v->getMEhisto_NumeratorPhiEndcap()->Fill(leadjphi);
2402 v->getMEhisto_NumeratorEtaForward()->Fill(leadjeta);
2403 v->getMEhisto_NumeratorPhiForward()->Fill(leadjphi);
2405 if((leadjpt > trigLowpTcut && !
isForward(leadjeta)) || (leadjpt > trigLowpTcutFwd &&
isForward(leadjeta))){
2406 v->getMEhisto_NumeratorEta_LowpTcut()->Fill(leadjeta);
2407 v->getMEhisto_NumeratorPhi_LowpTcut()->Fill(leadjphi);
2408 v->getMEhisto_NumeratorEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
2410 if((leadjpt > trigMedpTcut && !
isForward(leadjeta)) || (leadjpt > trigMedpTcutFwd &&
isForward(leadjeta))){
2411 v->getMEhisto_NumeratorEta_MedpTcut()->Fill(leadjeta);
2412 v->getMEhisto_NumeratorPhi_MedpTcut()->Fill(leadjphi);
2413 v->getMEhisto_NumeratorEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
2415 if((leadjpt > trigHighpTcut && !
isForward(leadjeta)) || (leadjpt > trigHighpTcutFwd &&
isForward(leadjeta))){
2416 v->getMEhisto_NumeratorEta_HighpTcut()->Fill(leadjeta);
2417 v->getMEhisto_NumeratorPhi_HighpTcut()->Fill(leadjphi);
2418 v->getMEhisto_NumeratorEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
2425 if(jetIDbool ==
true && (v->getTriggerType().compare(
"DiJet_Trigger") == 0) &&
calojet.size()>1){
2430 if(numpassed==
true){
2439 if(
verbose_)
std::cout <<
"fillMEforEffWrtMuTrigger: PFJet -------------------" << std::endl;
2445 bool jetIDbool =
false;
2464 if((v->getTriggerType().compare(
"SingleJet_Trigger") == 0) &&
pfjet.size()){
2467 _pfMHT =
sqrt(pfMHTx*pfMHTx + pfMHTy*pfMHTy);
2468 v->getMEhisto_DenominatorPFMHT()->Fill(
_pfMHT);
2474 if(
verbose_)
cout<<
"passed PFJet ID -------------------" << endl;
2476 v->getMEhisto_DenominatorPFPt()->Fill(leadjpt);
2477 if (
isBarrel(leadjeta)) v->getMEhisto_DenominatorPFPtBarrel()->Fill(leadjpt);
2478 if (
isEndCap(leadjeta)) v->getMEhisto_DenominatorPFPtEndcap()->Fill(leadjpt);
2479 if (
isForward(leadjeta)) v->getMEhisto_DenominatorPFPtForward()->Fill(leadjpt);
2480 v->getMEhisto_DenominatorPFEta()->Fill(leadjeta);
2481 v->getMEhisto_DenominatorPFPhi()->Fill(leadjphi);
2482 v->getMEhisto_DenominatorPFEtaPhi()->Fill(leadjeta,leadjphi);
2484 v->getMEhisto_DenominatorPFEtaBarrel()->Fill(leadjeta);
2485 v->getMEhisto_DenominatorPFPhiBarrel()->Fill(leadjphi);
2488 v->getMEhisto_DenominatorPFEtaEndcap()->Fill(leadjeta);
2489 v->getMEhisto_DenominatorPFPhiEndcap()->Fill(leadjphi);
2492 v->getMEhisto_DenominatorPFEtaForward()->Fill(leadjeta);
2493 v->getMEhisto_DenominatorPFPhiForward()->Fill(leadjphi);
2495 if((leadjpt > trigLowpTPFcut && !
isForward(leadjeta)) || (leadjpt > trigLowpTPFcutFwd &&
isForward(leadjeta))){
2496 v->getMEhisto_DenominatorPFEta_LowpTcut()->Fill(leadjeta);
2497 v->getMEhisto_DenominatorPFPhi_LowpTcut()->Fill(leadjphi);
2498 v->getMEhisto_DenominatorPFEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
2500 if((leadjpt > trigMedpTPFcut && !
isForward(leadjeta)) || (leadjpt > trigMedpTPFcutFwd &&
isForward(leadjeta))){
2501 v->getMEhisto_DenominatorPFEta_MedpTcut()->Fill(leadjeta);
2502 v->getMEhisto_DenominatorPFPhi_MedpTcut()->Fill(leadjphi);
2503 v->getMEhisto_DenominatorPFEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
2505 if((leadjpt > trigHighpTPFcut && !
isForward(leadjeta)) || (leadjpt > trigHighpTPFcutFwd &&
isForward(leadjeta))){
2506 v->getMEhisto_DenominatorPFEta_HighpTcut()->Fill(leadjeta);
2507 v->getMEhisto_DenominatorPFPhi_HighpTcut()->Fill(leadjphi);
2508 v->getMEhisto_DenominatorPFEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
2513 double dRmin = 99999.;
2514 double dPhimin = 9999.;
2515 if(v->getPath().find(
"L1") != std::string::npos){
2517 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
2521 for(trigger::Keys::const_iterator ki = kl1.begin();ki != kl1.end(); ++ki){
2522 double dR =
deltaR(toc[*ki].
eta(), toc[*ki].
phi(), leadjeta, leadjphi);
2531 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
2535 for(trigger::Keys::const_iterator kj = khlt.begin();kj != khlt.end(); ++kj){
2536 double dR =
deltaR(toc[*kj].
eta(), toc[*kj].
phi(), leadjeta, leadjphi);
2545 v->getMEhisto_PFDeltaPhi()->Fill(dPhimin);
2546 v->getMEhisto_PFDeltaR()->Fill(dRmin);
2549 if(dRmin < 0.1 || (v->getPath().find(
"L1") != std::string::npos && dRmin < 0.4)){
2550 v->getMEhisto_NumeratorPFPt()->Fill(leadjpt);
2551 if (
isBarrel(leadjeta)) v->getMEhisto_NumeratorPFPtBarrel()->Fill(leadjpt);
2552 if (
isEndCap(leadjeta)) v->getMEhisto_NumeratorPFPtEndcap()->Fill(leadjpt);
2553 if (
isForward(leadjeta)) v->getMEhisto_NumeratorPFPtForward()->Fill(leadjpt);
2554 v->getMEhisto_NumeratorPFEta()->Fill(leadjeta);
2555 v->getMEhisto_NumeratorPFPhi()->Fill(leadjphi);
2556 v->getMEhisto_NumeratorPFEtaPhi()->Fill(leadjeta,leadjphi);
2558 v->getMEhisto_NumeratorPFEtaBarrel()->Fill(leadjeta);
2559 v->getMEhisto_NumeratorPFPhiBarrel()->Fill(leadjphi);
2562 v->getMEhisto_NumeratorPFEtaEndcap()->Fill(leadjeta);
2563 v->getMEhisto_NumeratorPFPhiEndcap()->Fill(leadjphi);
2566 v->getMEhisto_NumeratorPFEtaForward()->Fill(leadjeta);
2567 v->getMEhisto_NumeratorPFPhiForward()->Fill(leadjphi);
2569 if((leadjpt > trigLowpTPFcut && !
isForward(leadjeta))
2570 || (leadjpt > trigLowpTPFcutFwd &&
isForward(leadjeta))){
2571 v->getMEhisto_NumeratorPFEta_LowpTcut()->Fill(leadjeta);
2572 v->getMEhisto_NumeratorPFPhi_LowpTcut()->Fill(leadjphi);
2573 v->getMEhisto_NumeratorPFEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
2575 if((leadjpt > trigMedpTPFcut && !
isForward(leadjeta))
2576 || (leadjpt > trigMedpTPFcutFwd &&
isForward(leadjeta))){
2577 v->getMEhisto_NumeratorPFEta_MedpTcut()->Fill(leadjeta);
2578 v->getMEhisto_NumeratorPFPhi_MedpTcut()->Fill(leadjphi);
2579 v->getMEhisto_NumeratorPFEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
2581 if((leadjpt > trigHighpTPFcut && !
isForward(leadjeta))
2582 || (leadjpt > trigHighpTPFcutFwd &&
isForward(leadjeta))){
2583 v->getMEhisto_NumeratorPFEta_HighpTcut()->Fill(leadjeta);
2584 v->getMEhisto_NumeratorPFPhi_HighpTcut()->Fill(leadjphi);
2585 v->getMEhisto_NumeratorPFEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
2591 if(jetIDbool ==
true && (v->getTriggerType().compare(
"DiJet_Trigger") == 0) &&
pfjet.size()>1){
2600 v->getMEhisto_DenominatorPFPt()->Fill((
PFJetPt[0] +
PFJetPt[1])/2.);
2603 v->getMEhisto_NumeratorPFPt()->Fill((
PFJetPt[0] +
PFJetPt[1])/2.);
2611 if(
verbose_)
std::cout <<
"fillMEforEffWrtMuTrigger: CaloMET -------------------" << std::endl;
2614 && (v->getPath().find(
"HLT_PFMET")==std::string::npos)){
2618 v->getMEhisto_DenominatorPt()->Fill(met.
et());
2619 v->getMEhisto_DenominatorPhi()->Fill(met.
phi());
2621 v->getMEhisto_NumeratorPt()->Fill(met.
et());
2622 v->getMEhisto_NumeratorPhi()->Fill(met.
phi());
2624 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
2627 double dPhimin = 9999.;
2629 for(trigger::Keys::const_iterator kj = khlt.begin();kj != khlt.end(); ++kj){
2635 v->getMEhisto_DeltaPhi()->Fill(dPhimin);
2641 if(
verbose_)
std::cout <<
"fillMEforEffWrtMuTrigger: PFMET -------------------" << std::endl;
2644 && (v->getPath().find(
"HLT_PFMET")!=std::string::npos)){
2647 const PFMET met = pfmetcol->front();
2648 v->getMEhisto_DenominatorPt()->Fill(met.
et());
2649 v->getMEhisto_DenominatorPhi()->Fill(met.
phi());
2651 v->getMEhisto_NumeratorPt()->Fill(met.
et());
2652 v->getMEhisto_NumeratorPhi()->Fill(met.
phi());
2654 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
2657 double dPhimin = 9999.;
2659 for(trigger::Keys::const_iterator kj = khlt.begin();kj != khlt.end(); ++kj){
2665 v->getMEhisto_DeltaPhi()->Fill(dPhimin);
2708 LogDebug(
"HLTJetMETDQMSource") <<
"HLTConfigProvider failed to initialize.";
2724 for (
unsigned int i=0;
i!=
n; ++
i) {
2725 bool denomFound =
false;
2726 bool numFound =
false;
2727 bool mbFound =
false;
2728 bool muFound =
false;
2734 std::vector<std::string>::const_iterator controlPathname =
pathFilter_.begin();
2735 for(;controlPathname!=
pathFilter_.end(); ++controlPathname){
2736 if(pathname.find((*controlPathname)) != std::string::npos){
2741 if(checkPath==
false)
continue;
2744 std::vector<std::string>::const_iterator rejectPathname =
pathRejectKeyword_.begin();
2746 if(pathname.find((*rejectPathname)) != std::string::npos){
2751 if(checkPath==
false)
continue;
2759 unsigned int usedPrescale = 1;
2760 unsigned int objectType = 0;
2765 if(pathname.find(
"Jet") != std::string::npos
2766 && !(pathname.find(
"DoubleJet") != std::string::npos)
2767 && !(pathname.find(
"DiJet") != std::string::npos)
2768 && !(pathname.find(
"DiPFJet") != std::string::npos)
2769 && !(pathname.find(
"BTag") != std::string::npos)
2770 && !(pathname.find(
"Mu") != std::string::npos)
2771 && !(pathname.find(
"Fwd") != std::string::npos)){
2772 triggerType =
"SingleJet_Trigger";
2775 if(pathname.find(
"DiJet") != std::string::npos
2776 || pathname.find(
"DiPFJet") != std::string::npos
2777 || pathname.find(
"DoubleJet") != std::string::npos){
2778 triggerType =
"DiJet_Trigger";
2781 if(pathname.find(
"MET") != std::string::npos){
2782 triggerType =
"MET_Trigger";
2785 if(pathname.find(
"HT") != std::string::npos) {
2786 triggerType =
"TET_Trigger";
2792 && !(pathname.find(
"DiJet") != std::string::npos)
2793 && !(pathname.find(
"DiPFJet") != std::string::npos)
2794 && !(pathname.find(
"DoubleJet") != std::string::npos)){
2800 && ((pathname.find(
"DiJet") != std::string::npos)
2801 || (pathname.find(
"DiPFJet") != std::string::npos))){
2822 for (
unsigned int i=0;
i!=
n; ++
i){
2824 if(HLTname == pathname) numFound =
true;
2825 if(HLTname == dpathname) denomFound =
true;
2826 if(HLTname.find(
MBTrigPaths_[0]) != std::string::npos) mbFound =
true;
2827 if(HLTname.find(
MuonTrigPaths_[0]) != std::string::npos) muFound =
true;
2837 for(std::vector<std::string>::iterator numpathmodule = numpathmodules.begin(); numpathmodule!= numpathmodules.end(); ++numpathmodule ) {
2846 filtername = *numpathmodule;
2862 if(objectType != 0 && denomFound){
2864 for(std::vector<std::string>::iterator numpathmodule = numpathmodules.begin(); numpathmodule!= numpathmodules.end(); ++numpathmodule ) {
2873 Denomfiltername = *numpathmodule;
2877 if(objectType != 0 && numFound){
2879 cout<<
"=Pathname= "<<pathname
2880 <<
" | =Denompath= "<<dpathname
2881 <<
" | =Filtername= "<<filtername
2882 <<
" | =Denomfiltername= "<<Denomfiltername
2883 <<
" | =L1pathname= "<<l1pathname
2884 <<
" | =ObjectType= "<<objectType<<endl;
2885 if(!((pathname.find(
"HT") != std::string::npos) || (pathname.find(
"Quad") != std::string::npos))){
2888 filtername, Denomfiltername,
processname_, objectType, triggerType));
2890 filtername, Denomfiltername,
processname_, objectType, triggerType));
2892 filtername, Denomfiltername,
processname_, objectType, triggerType));
2894 filtername, Denomfiltername,
processname_, objectType, triggerType));
2897 filtername, Denomfiltername,
processname_, objectType, triggerType));
2901 if (
verbose_)
cout <<
"get names for efficicncy------------------"<< endl;
2906 unsigned int usedPrescale = 1;
2907 unsigned int objectType = 0;
2911 for (std::vector<std::pair<std::string, std::string> >::iterator custompathnamepair =
custompathnamepairs_.begin();
2914 std::string dpathname = custompathnamepair->second;
2915 bool numFound =
false;
2916 bool denomFound =
false;
2918 for (
unsigned int i=0;
i!=
n; ++
i) {
2920 if(HLTname.find(pathname)!=std::string::npos){numFound =
true; pathname = HLTname;}
2921 if(HLTname.find(dpathname)!=std::string::npos){denomFound =
true; dpathname = HLTname;}
2923 if(numFound && denomFound){
2924 if (pathname.find(
"Jet") != std::string::npos
2925 && !(pathname.find(
"DiJet") != std::string::npos)
2926 && !(pathname.find(
"DiPFJet") != std::string::npos)
2927 && !(pathname.find(
"DoubleJet") != std::string::npos)
2928 && !(pathname.find(
"BTag") != std::string::npos)
2929 && !(pathname.find(
"Mu") != std::string::npos)
2930 && !(pathname.find(
"Fwd") != std::string::npos)){
2931 triggerType =
"SingleJet_Trigger";
2934 if (pathname.find(
"DiJet") != std::string::npos
2935 || pathname.find(
"DiPFJet") != std::string::npos
2936 || pathname.find(
"DoubleJet") != std::string::npos ){
2937 triggerType =
"DiJet_Trigger";
2940 if (pathname.find(
"MET") != std::string::npos ){
2941 triggerType =
"MET_Trigger";
2944 if (pathname.find(
"TET") != std::string::npos ){
2945 triggerType =
"TET_Trigger";
2951 for(std::vector<std::string>::iterator numpathmodule = numpathmodules.begin(); numpathmodule!= numpathmodules.end(); ++numpathmodule ) {
2960 filtername = *numpathmodule;
2963 if(objectType != 0){
2965 for(std::vector<std::string>::iterator numpathmodule = numpathmodules.begin(); numpathmodule!= numpathmodules.end(); ++numpathmodule ) {
2974 Denomfiltername = *numpathmodule;
2978 <<
"==denompath=="<<dpathname
2979 <<
"==filtername=="<<filtername
2980 <<
"==denomfiltername=="<<Denomfiltername
2981 <<
"==l1pathname=="<<l1pathname
2982 <<
"==objectType=="<<objectType<<endl;
2984 filtername, Denomfiltername,
processname_, objectType, triggerType));
2990 if(
verbose_)
cout<<
"== end hltPathsEff_.push_back ======" << endl;
3000 double TrigMin_ = -0.5;
3005 rate_All = dbe->
book1D(histonm.c_str(),histot.c_str(),TrigBins_,TrigMin_,TrigMax_);
3007 histonm =
"JetMET_TriggerRate_Correlation";
3008 histot =
"JetMET TriggerRate Correlation Summary;y&&!x;x&&y";
3009 correlation_All = dbe->
book2D(histonm.c_str(),histot.c_str(),TrigBins_,TrigMin_,TrigMax_,TrigBins_,TrigMin_,TrigMax_);
3011 histonm =
"JetMET_NVertices";
3012 histot =
"No. of vertices";
3015 histonm =
"JetMET_PVZ";
3016 histot =
"Primary Vertex Z pos";
3017 PVZ = dbe->
book1D(histonm.c_str(),histot.c_str(),100,-50.,50.);
3020 histonm=
"JetMET_TriggerRate_WrtMuTrigger";
3021 histot=
"JetMET TriggerRate Summary Wrt Muon Trigger ";
3024 histonm=
"JetMET_TriggerRate_Correlation_WrtMuTrigger";
3025 histot=
"JetMET TriggerRate Correlation Summary Wrt Muon Trigger;y&&!x;x&&y";
3028 histonm=
"JetMET_TriggerRate_WrtMBTrigger";
3029 histot=
"JetMET TriggerRate Summary Wrt MB Trigger";
3032 histonm=
"JetMET_TriggerRate_Correlation_WrtMBTrigger";
3033 histot=
"JetMET TriggerRate Correlation Wrt MB Trigger;y&&!x;x&&y";
3042 labelnm =
v->getPath();
3044 for(
int ibin=1; ibin<nbins+1; ibin++){
3045 const char * binLabel =
rate_All->
getTH1()->GetXaxis()->GetBinLabel(ibin);
3047 if(binLabel_str.compare(labelnm)==0)
break;
3048 if(binLabel[0]==
'\0'){
3075 double PtMax_ = 1000.;
3079 double EtaMin_ = -5.;
3080 double EtaMax_ = 5.;
3083 double PhiMin_ = -3.5;
3084 double PhiMax_ = 3.5;
3087 double ResMin_ = -1.5;
3088 double ResMax_ = 1.5;
3105 if(
v->getObjectType() ==
trigger::TriggerJet &&
v->getTriggerType().compare(
"SingleJet_Trigger") == 0){
3107 histoname = labelname+
"_recObjN";
3108 title = labelname+
"_recObjN;Reco multiplicity()"+trigPath;
3112 histoname = labelname+
"_recObjPt";
3113 title = labelname+
"_recObjPt; Reco Pt[GeV/c]"+trigPath;
3117 histoname = labelname+
"_recObjPtBarrel";
3118 title = labelname+
"_recObjPtBarrel;Reco Pt[GeV/c]"+trigPath;
3119 MonitorElement * PtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3122 histoname = labelname+
"_recObjPtEndcap";
3123 title = labelname+
"_recObjPtEndcap;Reco Pt[GeV/c]"+trigPath;
3124 MonitorElement * PtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3127 histoname = labelname+
"_recObjPtForward";
3128 title = labelname+
"_recObjPtForward;Reco Pt[GeV/c]"+trigPath;
3129 MonitorElement * PtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3132 histoname = labelname+
"_recObjEta";
3133 title = labelname+
"_recObjEta;Reco #eta"+trigPath;
3137 histoname = labelname+
"_recObjPhi";
3138 title = labelname+
"_recObjPhi;Reco #Phi"+trigPath;
3142 histoname = labelname+
"_recObjEtaPhi";
3143 title = labelname+
"_recObjEtaPhi;Reco #eta;Reco #Phi"+trigPath;
3144 MonitorElement * EtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
3147 histoname = labelname+
"_l1ObjPt";
3148 title = labelname+
"_l1ObjPt;L1 Pt[GeV/c]"+trigPath;
3152 histoname = labelname+
"_l1ObjEta";
3153 title = labelname+
"_l1ObjEta;L1 #eta"+trigPath;
3154 MonitorElement * Eta_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
3157 histoname = labelname+
"_l1ObjPhi";
3158 title = labelname+
"_l1ObjPhi;L1 #Phi"+trigPath;
3159 MonitorElement * Phi_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3162 histoname = labelname+
"_l1ObjEtaPhi";
3163 title = labelname+
"_l1ObjEtaPhi;L1 #eta;L1 #Phi"+trigPath;
3164 MonitorElement * EtaPhi_L1 = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
3167 histoname = labelname+
"_l1ObjN";
3168 title = labelname+
"_l1ObjN;L1 multiplicity"+trigPath;
3172 histoname = labelname+
"_l1ObjPtBarrel";
3173 title = labelname+
"_l1ObjPtBarrel;L1 Pt[GeV/c]"+trigPath;
3174 MonitorElement * PtBarrel_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3177 histoname = labelname+
"_l1ObjPtEndcap";
3178 title = labelname+
"_l1ObjPtEndcap;L1 Pt[GeV/c]"+trigPath;
3179 MonitorElement * PtEndcap_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3182 histoname = labelname+
"_l1ObjPtForward";
3183 title = labelname+
"_l1ObjPtForward;L1 Pt[GeV/c]"+trigPath;
3184 MonitorElement * PtForward_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3187 histoname = labelname+
"_hltObjN";
3188 title = labelname+
"_hltObjN;HLT multiplicity"+trigPath;
3192 histoname = labelname+
"_hltObjPtBarrel";
3193 title = labelname+
"_hltObjPtBarrel;HLT Pt[GeV/c]"+trigPath;
3194 MonitorElement * PtBarrel_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3197 histoname = labelname+
"_hltObjPtEndcap";
3198 title = labelname+
"_hltObjPtEndcap;HLT Pt[GeV/c]"+trigPath;
3199 MonitorElement * PtEndcap_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3202 histoname = labelname+
"_hltObjPtForward";
3203 title = labelname+
"_hltObjPtForward;HLT Pt[GeV/c]"+trigPath;
3204 MonitorElement * PtForward_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3207 histoname = labelname+
"_hltObjPt";
3208 title = labelname+
"_hltObjPt;HLT Pt[GeV/c]"+trigPath;
3212 histoname = labelname+
"_hltObjEta";
3213 title = labelname+
"_hltObjEta;HLT #eta"+trigPath;
3214 MonitorElement * Eta_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
3217 histoname = labelname+
"_hltObjPhi";
3218 title = labelname+
"_hltObjPhi;HLT #Phi"+trigPath;
3219 MonitorElement * Phi_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3222 histoname = labelname+
"_hltObjEtaPhi";
3223 title = labelname+
"_hltObjEtaPhi;HLT #eta;HLT #Phi"+trigPath;
3224 MonitorElement * EtaPhi_HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
3227 histoname = labelname+
"_l1HLTPtResolution";
3228 title = labelname+
"_l1HLTPtResolution;(Pt(L1)-Pt(HLT))/Pt(L1)"+trigPath;
3229 MonitorElement * PtResolution_L1HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3230 PtResolution_L1HLT->
getTH1();
3232 histoname = labelname+
"_l1HLTEtaResolution";
3233 title = labelname+
"_l1HLTEtaResolution;(#eta(L1)-#eta(HLT))/#eta(L1)"+trigPath;
3234 MonitorElement * EtaResolution_L1HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3235 EtaResolution_L1HLT->
getTH1();
3237 histoname = labelname+
"_l1HLTPhiResolution";
3238 title = labelname+
"_l1HLTPhiResolution;(#Phi(L1)-#Phi(HLT))/#Phi(L1)"+trigPath;
3239 MonitorElement * PhiResolution_L1HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3240 PhiResolution_L1HLT->
getTH1();
3242 histoname = labelname+
"_l1HLTPtCorrelation";
3243 title = labelname+
"_l1HLTPtCorrelation;Pt(L1)[GeV/c];Pt(HLT)[GeV/c]"+trigPath;
3244 MonitorElement * PtCorrelation_L1HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_,Ptbins_,PtMin_,PtMax_);
3245 PtCorrelation_L1HLT->
getTH1();
3247 histoname = labelname+
"_l1HLTEtaCorrelation";
3248 title = labelname+
"_l1HLTEtaCorrelation;#eta(L1);#eta(HLT)"+trigPath;
3249 MonitorElement * EtaCorrelation_L1HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Etabins_,EtaMin_,EtaMax_);
3250 EtaCorrelation_L1HLT->
getTH1();
3252 histoname = labelname+
"_l1HLTPhiCorrelation";
3253 title = labelname+
"_l1HLTPhiCorrelation;#Phi(L1);#Phi(HLT)"+trigPath;
3254 MonitorElement * PhiCorrelation_L1HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_,Phibins_,PhiMin_,PhiMax_);
3255 PhiCorrelation_L1HLT->
getTH1();
3257 histoname = labelname+
"_hltRecObjPtResolution";
3258 title = labelname+
"_hltRecObjPtResolution;(Pt(HLT)-Pt(Reco))/Pt(HLT)"+trigPath;
3259 MonitorElement * PtResolution_HLTRecObj = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3260 PtResolution_HLTRecObj->
getTH1();
3262 histoname = labelname+
"_hltRecObjEtaResolution";
3263 title = labelname+
"_hltRecObjEtaResolution;(#eta(HLT)-#eta(Reco))/#eta(HLT)"+trigPath;
3264 MonitorElement * EtaResolution_HLTRecObj = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3265 EtaResolution_HLTRecObj->
getTH1();
3267 histoname = labelname+
"_hltRecObjPhiResolution";
3268 title = labelname+
"_hltRecObjPhiResolution;(#Phi(HLT)-#Phi(Reco))/#Phi(HLT)"+trigPath;
3269 MonitorElement * PhiResolution_HLTRecObj = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3270 PhiResolution_HLTRecObj->
getTH1();
3272 histoname = labelname+
"_hltRecObjPtCorrelation";
3273 title = labelname+
"_hltRecObjPtCorrelation;Pt(HLT)[GeV/c];Pt(Reco)[GeV/c]"+trigPath;
3274 MonitorElement * PtCorrelation_HLTRecObj = dbe->
book2D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_,Ptbins_,PtMin_,PtMax_);
3275 PtCorrelation_HLTRecObj->
getTH1();
3277 histoname = labelname+
"_hltRecObjEtaCorrelation";
3278 title = labelname+
"_hltRecObjEtaCorrelation;#eta(HLT);#eta(Reco)"+trigPath;
3279 MonitorElement * EtaCorrelation_HLTRecObj = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Etabins_,EtaMin_,EtaMax_);
3280 EtaCorrelation_HLTRecObj->
getTH1();
3282 histoname = labelname+
"_hltRecObjPhiCorrelation";
3283 title = labelname+
"_hltRecObjPhiCorrelation;#Phi(HLT);#Phi(Reco)"+trigPath;
3284 MonitorElement * PhiCorrelation_HLTRecObj = dbe->
book2D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_,Phibins_,PhiMin_,PhiMax_);
3285 PhiCorrelation_HLTRecObj->
getTH1();
3312 EtaResolution_L1HLT,
3313 PhiResolution_L1HLT,
3314 PtResolution_HLTRecObj,
3315 EtaResolution_HLTRecObj,
3316 PhiResolution_HLTRecObj,
3317 PtCorrelation_L1HLT,
3318 EtaCorrelation_L1HLT,
3319 PhiCorrelation_L1HLT,
3320 PtCorrelation_HLTRecObj,
3321 EtaCorrelation_HLTRecObj,
3322 PhiCorrelation_HLTRecObj,
3337 histoname = labelname+
"_RecObjAveragePt";
3338 title = labelname+
"_RecObjAveragePt;Reco Average Pt[GeV/c]"+trigPath;
3339 MonitorElement * jetAveragePt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3342 histoname = labelname+
"_RecObjAverageEta";
3343 title = labelname+
"_RecObjAverageEta;Reco Average #eta"+trigPath;
3344 MonitorElement * jetAverageEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
3347 histoname = labelname+
"_RecObjPhiDifference";
3348 title = labelname+
"_RecObjPhiDifference;Reco #Delta#Phi"+trigPath;
3349 MonitorElement * jetPhiDifference = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3350 jetPhiDifference->
getTH1();
3352 histoname = labelname+
"_hltObjAveragePt";
3353 title = labelname+
"_hltObjAveragePt;HLT Average Pt[GeV/c]"+trigPath;
3354 MonitorElement * hltAveragePt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3357 histoname = labelname+
"_hltObjAverageEta";
3358 title = labelname+
"_hltObjAverageEta;HLT Average #eta"+trigPath;
3359 MonitorElement * hltAverageEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
3362 histoname = labelname+
"_hltObjPhiDifference";
3363 title = labelname+
"_hltObjPhiDifference;Reco #Delta#Phi"+trigPath;
3364 MonitorElement * hltPhiDifference = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3365 hltPhiDifference->
getTH1();
3417 histoname = labelname+
"_recObjPt";
3418 title = labelname+
"_recObjPt;Reco Pt[GeV/c]"+trigPath;
3422 histoname = labelname+
"_recObjPhi";
3423 title = labelname+
"_recObjPhi;Reco #Phi"+trigPath;
3427 histoname = labelname+
"_l1ObjPt";
3428 title = labelname+
"_l1ObjPt;L1 Pt[GeV/c]"+trigPath;
3432 histoname = labelname+
"_l1ObjPhi";
3433 title = labelname+
"_l1ObjPhi;L1 #Phi"+trigPath;
3434 MonitorElement * Phi_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3437 histoname = labelname+
"_hltObjPt";
3438 title = labelname+
"_hltObjPt;HLT Pt[GeV/c]"+trigPath;
3442 histoname = labelname+
"_hltObjPhi";
3443 title = labelname+
"_hltObjPhi;HLT #Phi"+trigPath;
3444 MonitorElement * Phi_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3447 histoname = labelname+
"_l1HLTPtResolution";
3448 title = labelname+
"_l1HLTPtResolution;(Pt(L1)-Pt(HLT))/Pt(L1)"+trigPath;
3449 MonitorElement * PtResolution_L1HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3450 PtResolution_L1HLT->
getTH1();
3452 histoname = labelname+
"_l1HLTPhiResolution";
3453 title = labelname+
"_l1HLTPhiResolution;(#Phi(L1)-#Phi(HLT))/#Phi(L1)"+trigPath;
3454 MonitorElement * PhiResolution_L1HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3455 PhiResolution_L1HLT->
getTH1();
3457 histoname = labelname+
"_l1HLTPtCorrelation";
3458 title = labelname+
"_l1HLTPtCorrelation;Pt(L1)[GeV/c];Pt(HLT)[GeV/c]"+trigPath;
3459 MonitorElement * PtCorrelation_L1HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_,Ptbins_,PtMin_,PtMax_);
3460 PtCorrelation_L1HLT->
getTH1();
3462 histoname = labelname+
"_l1HLTPhiCorrelation";
3463 title = labelname+
"_l1HLTPhiCorrelation;#Phi(L1);#Phi(HLT)"+trigPath;
3464 MonitorElement * PhiCorrelation_L1HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_,Phibins_,PhiMin_,PhiMax_);
3465 PhiCorrelation_L1HLT->
getTH1();
3467 histoname = labelname+
"_hltRecObjPtResolution";
3468 title = labelname+
"_hltRecObjPtResolution;(Pt(HLT)-Pt(Reco))/Pt(HLT)"+trigPath;
3469 MonitorElement * PtResolution_HLTRecObj = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3470 PtResolution_HLTRecObj->
getTH1();
3472 histoname = labelname+
"_hltRecObjPhiResolution";
3473 title = labelname+
"_hltRecObjPhiResolution;(#Phi(HLT)-#Phi(Reco))/#Phi(HLT)"+trigPath;
3474 MonitorElement * PhiResolution_HLTRecObj = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3475 PhiResolution_HLTRecObj->
getTH1();
3477 histoname = labelname+
"_hltRecObjPtCorrelation";
3478 title = labelname+
"_hltRecObjPtCorrelation;Pt(HLT)[GeV/c];Pt(Reco)[GeV/c]"+trigPath;
3479 MonitorElement * PtCorrelation_HLTRecObj = dbe->
book2D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_,Ptbins_,PtMin_,PtMax_);
3480 PtCorrelation_HLTRecObj->
getTH1();
3482 histoname = labelname+
"_hltRecObjPhiCorrelation";
3483 title = labelname+
"_hltRecObjPhiCorrelation;#Phi(HLT);#Phi(Reco)"+trigPath;
3484 MonitorElement * PhiCorrelation_HLTRecObj = dbe->
book2D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_,Phibins_,PhiMin_,PhiMax_);
3485 PhiCorrelation_HLTRecObj->
getTH1();
3513 PhiResolution_L1HLT,
3514 PtResolution_HLTRecObj,
3516 PhiResolution_HLTRecObj,
3517 PtCorrelation_L1HLT,
3519 PhiCorrelation_L1HLT,
3520 PtCorrelation_HLTRecObj,
3522 PhiCorrelation_HLTRecObj,
3545 double PtMax_ = 1000.;
3548 double EtaMin_ = -5.;
3549 double EtaMax_ = 5.;
3552 double PhiMin_ = -3.5;
3553 double PhiMax_ = 3.5;
3556 double ResMin_ = -1.5;
3557 double ResMax_ = 1.5;
3575 if(
v->getObjectType() ==
trigger::TriggerJet &&
v->getTriggerType().compare(
"SingleJet_Trigger") == 0){
3577 histoname = labelname+
"_recObjN";
3578 title = labelname+
"_recObjN;Reco multiplicity()"+trigPath;
3582 histoname = labelname+
"_recObjPt";
3583 title = labelname+
"_recObjPt; Reco Pt[GeV/c]"+trigPath;
3587 histoname = labelname+
"_recObjPtBarrel";
3588 title = labelname+
"_recObjPtBarrel;Reco Pt[GeV/c]"+trigPath;
3589 MonitorElement * PtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3592 histoname = labelname+
"_recObjPtEndcap";
3593 title = labelname+
"_recObjPtEndcap;Reco Pt[GeV/c]"+trigPath;
3594 MonitorElement * PtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3597 histoname = labelname+
"_recObjPtForward";
3598 title = labelname+
"_recObjPtForward;Reco Pt[GeV/c]"+trigPath;
3599 MonitorElement * PtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3602 histoname = labelname+
"_recObjEta";
3603 title = labelname+
"_recObjEta;Reco #eta"+trigPath;
3607 histoname = labelname+
"_recObjPhi";
3608 title = labelname+
"_recObjPhi;Reco #Phi"+trigPath;
3612 histoname = labelname+
"_recObjEtaPhi";
3613 title = labelname+
"_recObjEtaPhi;Reco #eta;Reco #Phi"+trigPath;
3614 MonitorElement * EtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
3617 histoname = labelname+
"_l1ObjN";
3618 title = labelname+
"_l1ObjN;L1 multiplicity"+trigPath;
3622 histoname = labelname+
"_l1ObjPt";
3623 title = labelname+
"_l1ObjPt;L1 Pt[GeV/c]"+trigPath;
3627 histoname = labelname+
"_l1ObjPtBarrel";
3628 title = labelname+
"_l1ObjPtBarrel;L1 Pt[GeV/c]"+trigPath;
3629 MonitorElement * PtBarrel_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3632 histoname = labelname+
"_l1ObjPtEndcap";
3633 title = labelname+
"_l1ObjPtEndcap;L1 Pt[GeV/c]"+trigPath;
3634 MonitorElement * PtEndcap_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3637 histoname = labelname+
"_l1ObjPtForward";
3638 title = labelname+
"_l1ObjPtForward;L1 Pt[GeV/c]"+trigPath;
3639 MonitorElement * PtForward_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3642 histoname = labelname+
"_l1ObjEta";
3643 title = labelname+
"_l1ObjEta;L1 #eta"+trigPath;
3644 MonitorElement * Eta_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
3647 histoname = labelname+
"_l1ObjPhi";
3648 title = labelname+
"_l1ObjPhi;L1 #Phi"+trigPath;
3649 MonitorElement * Phi_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3652 histoname = labelname+
"_l1ObjEtaPhi";
3653 title = labelname+
"_l1ObjEtaPhi;L1 #eta;L1 #Phi"+trigPath;
3654 MonitorElement * EtaPhi_L1 = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
3657 histoname = labelname+
"_hltObjN";
3658 title = labelname+
"_hltObjN;HLT multiplicity"+trigPath;
3662 histoname = labelname+
"_hltObjPt";
3663 title = labelname+
"_hltObjPt;HLT Pt[GeV/c]"+trigPath;
3667 histoname = labelname+
"_hltObjPtBarrel";
3668 title = labelname+
"_hltObjPtBarrel;HLT Pt[GeV/c]"+trigPath;
3669 MonitorElement * PtBarrel_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3672 histoname = labelname+
"_hltObjPtEndcap";
3673 title = labelname+
"_hltObjPtEndcap;HLT Pt[GeV/c]"+trigPath;
3674 MonitorElement * PtEndcap_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3677 histoname = labelname+
"_hltObjPtForward";
3678 title = labelname+
"_hltObjPtForward;HLT Pt[GeV/c]"+trigPath;
3679 MonitorElement * PtForward_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3682 histoname = labelname+
"_hltObjEta";
3683 title = labelname+
"_hltObjEta;HLT #eta"+trigPath;
3684 MonitorElement * Eta_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
3687 histoname = labelname+
"_hltObjPhi";
3688 title = labelname+
"_hltObjPhi;HLT #Phi"+trigPath;
3689 MonitorElement * Phi_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3692 histoname = labelname+
"_hltObjEtaPhi";
3693 title = labelname+
"_hltObjEtaPhi;HLT #eta;HLT #Phi"+trigPath;
3694 MonitorElement * EtaPhi_HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
3697 histoname = labelname+
"_l1HLTPtResolution";
3698 title = labelname+
"_l1HLTPtResolution;(Pt(L1)-Pt(HLT))/Pt(L1)"+trigPath;
3699 MonitorElement * PtResolution_L1HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3700 PtResolution_L1HLT->
getTH1();
3702 histoname = labelname+
"_l1HLTEtaResolution";
3703 title = labelname+
"_l1HLTEtaResolution;(#eta(L1)-#eta(HLT))/#eta(L1)"+trigPath;
3704 MonitorElement * EtaResolution_L1HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3705 EtaResolution_L1HLT->
getTH1();
3707 histoname = labelname+
"_l1HLTPhiResolution";
3708 title = labelname+
"_l1HLTPhiResolution;(#Phi(L1)-#Phi(HLT))/#Phi(L1)"+trigPath;
3709 MonitorElement * PhiResolution_L1HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3710 PhiResolution_L1HLT->
getTH1();
3712 histoname = labelname+
"_l1HLTPtCorrelation";
3713 title = labelname+
"_l1HLTPtCorrelation;Pt(L1)[GeV/c];Pt(HLT)[GeV/c]"+trigPath;
3714 MonitorElement * PtCorrelation_L1HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_,Ptbins_,PtMin_,PtMax_);
3715 PtCorrelation_L1HLT->
getTH1();
3717 histoname = labelname+
"_l1HLTEtaCorrelation";
3718 title = labelname+
"_l1HLTEtaCorrelation;#eta(L1);#eta(HLT)"+trigPath;
3719 MonitorElement * EtaCorrelation_L1HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Etabins_,EtaMin_,EtaMax_);
3720 EtaCorrelation_L1HLT->
getTH1();
3722 histoname = labelname+
"_l1HLTPhiCorrelation";
3723 title = labelname+
"_l1HLTPhiCorrelation;#Phi(L1);#Phi(HLT)"+trigPath;
3724 MonitorElement * PhiCorrelation_L1HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_,Phibins_,PhiMin_,PhiMax_);
3725 PhiCorrelation_L1HLT->
getTH1();
3727 histoname = labelname+
"_hltRecObjPtResolution";
3728 title = labelname+
"_hltRecObjPtResolution;(Pt(HLT)-Pt(Reco))/Pt(HLT)"+trigPath;
3729 MonitorElement * PtResolution_HLTRecObj = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3730 PtResolution_HLTRecObj->
getTH1();
3732 histoname = labelname+
"_hltRecObjEtaResolution";
3733 title = labelname+
"_hltRecObjEtaResolution;(#eta(HLT)-#eta(Reco))/#eta(HLT)"+trigPath;
3734 MonitorElement * EtaResolution_HLTRecObj = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3735 EtaResolution_HLTRecObj->
getTH1();
3737 histoname = labelname+
"_hltRecObjPhiResolution";
3738 title = labelname+
"_hltRecObjPhiResolution;(#Phi(HLT)-#Phi(Reco))/#Phi(HLT)"+trigPath;
3739 MonitorElement * PhiResolution_HLTRecObj = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3740 PhiResolution_HLTRecObj->
getTH1();
3742 histoname = labelname+
"_hltRecObjPtCorrelation";
3743 title = labelname+
"_hltRecObjPtCorrelation;Pt(HLT)[GeV/c];Pt(Reco)[GeV/c]"+trigPath;
3744 MonitorElement * PtCorrelation_HLTRecObj = dbe->
book2D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_,Ptbins_,PtMin_,PtMax_);
3745 PtCorrelation_HLTRecObj->
getTH1();
3747 histoname = labelname+
"_hltRecObjEtaCorrelation";
3748 title = labelname+
"_hltRecObjEtaCorrelation;#eta(HLT);#eta(Reco)"+trigPath;
3749 MonitorElement * EtaCorrelation_HLTRecObj = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Etabins_,EtaMin_,EtaMax_);
3750 EtaCorrelation_HLTRecObj->
getTH1();
3752 histoname = labelname+
"_hltRecObjPhiCorrelation";
3753 title = labelname+
"_hltRecObjPhiCorrelation;#Phi(HLT);#Phi(Reco)"+trigPath;
3754 MonitorElement * PhiCorrelation_HLTRecObj = dbe->
book2D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_,Phibins_,PhiMin_,PhiMax_);
3755 PhiCorrelation_HLTRecObj->
getTH1();
3782 EtaResolution_L1HLT,
3783 PhiResolution_L1HLT,
3784 PtResolution_HLTRecObj,
3785 EtaResolution_HLTRecObj,
3786 PhiResolution_HLTRecObj,
3787 PtCorrelation_L1HLT,
3788 EtaCorrelation_L1HLT,
3789 PhiCorrelation_L1HLT,
3790 PtCorrelation_HLTRecObj,
3791 EtaCorrelation_HLTRecObj,
3792 PhiCorrelation_HLTRecObj,
3807 histoname = labelname+
"_RecObjAveragePt";
3808 title = labelname+
"_RecObjAveragePt;Reco Average Pt[GeV/c]"+trigPath;
3809 MonitorElement * jetAveragePt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3812 histoname = labelname+
"_RecObjAverageEta";
3813 title = labelname+
"_RecObjAverageEta;Reco Average #eta"+trigPath;
3814 MonitorElement * jetAverageEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
3817 histoname = labelname+
"_RecObjPhiDifference";
3818 title = labelname+
"_RecObjPhiDifference;Reco #Delta#Phi"+trigPath;
3819 MonitorElement * jetPhiDifference = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3820 jetPhiDifference->
getTH1();
3822 histoname = labelname+
"_hltObjAveragePt";
3823 title = labelname+
"_hltObjAveragePt;HLT Average Pt[GeV/c]"+trigPath;
3824 MonitorElement * hltAveragePt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3827 histoname = labelname+
"_hltObjAverageEta";
3828 title = labelname+
"_hltObjAverageEta;HLT Average #eta"+trigPath;
3829 MonitorElement * hltAverageEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
3832 histoname = labelname+
"_hltObjPhiDifference";
3833 title = labelname+
"_hltObjPhiDifference;Reco #Delta#Phi"+trigPath;
3834 MonitorElement * hltPhiDifference = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3835 hltPhiDifference->
getTH1();
3887 histoname = labelname+
"_recObjPt";
3888 title = labelname+
"_recObjPt;Reco Pt[GeV/c]"+trigPath;
3892 histoname = labelname+
"_recObjPhi";
3893 title = labelname+
"_recObjPhi;Reco #Phi"+trigPath;
3897 histoname = labelname+
"_l1ObjPt";
3898 title = labelname+
"_l1ObjPt;L1 Pt[GeV/c]"+trigPath;
3902 histoname = labelname+
"_l1ObjPhi";
3903 title = labelname+
"_l1ObjPhi;L1 #Phi"+trigPath;
3904 MonitorElement * Phi_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3907 histoname = labelname+
"_hltObjPt";
3908 title = labelname+
"_hltObjPt;HLT Pt[GeV/c]"+trigPath;
3912 histoname = labelname+
"_hltObjPhi";
3913 title = labelname+
"_hltObjPhi;HLT #Phi"+trigPath;
3914 MonitorElement * Phi_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3917 histoname = labelname+
"_l1HLTPtResolution";
3918 title = labelname+
"_l1HLTPtResolution;(Pt(L1)-Pt(HLT))/Pt(L1)"+trigPath;
3919 MonitorElement * PtResolution_L1HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3920 PtResolution_L1HLT->
getTH1();
3922 histoname = labelname+
"_l1HLTPhiResolution";
3923 title = labelname+
"_l1HLTPhiResolution;(#Phi(L1)-#Phi(HLT))/#Phi(L1)"+trigPath;
3924 MonitorElement * PhiResolution_L1HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3925 PhiResolution_L1HLT->
getTH1();
3927 histoname = labelname+
"_l1HLTPtCorrelation";
3928 title = labelname+
"_l1HLTPtCorrelation;Pt(L1)[GeV/c];Pt(HLT)[GeV/c]"+trigPath;
3929 MonitorElement * PtCorrelation_L1HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_,Ptbins_,PtMin_,PtMax_);
3930 PtCorrelation_L1HLT->
getTH1();
3932 histoname = labelname+
"_l1HLTPhiCorrelation";
3933 title = labelname+
"_l1HLTPhiCorrelation;#Phi(L1);#Phi(HLT)"+trigPath;
3934 MonitorElement * PhiCorrelation_L1HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_,Phibins_,PhiMin_,PhiMax_);
3935 PhiCorrelation_L1HLT->
getTH1();
3937 histoname = labelname+
"_hltRecObjPtResolution";
3938 title = labelname+
"_hltRecObjPtResolution;(Pt(HLT)-Pt(Reco))/Pt(HLT)"+trigPath;
3939 MonitorElement * PtResolution_HLTRecObj = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3940 PtResolution_HLTRecObj->
getTH1();
3942 histoname = labelname+
"_hltRecObjPhiResolution";
3943 title = labelname+
"_hltRecObjPhiResolution;(#Phi(HLT)-#Phi(Reco))/#Phi(HLT)"+trigPath;
3944 MonitorElement * PhiResolution_HLTRecObj = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3945 PhiResolution_HLTRecObj->
getTH1();
3947 histoname = labelname+
"_hltRecObjPtCorrelation";
3948 title = labelname+
"_hltRecObjPtCorrelation;Pt(HLT)[GeV/c];Pt(Reco)[GeV/c]"+trigPath;
3949 MonitorElement * PtCorrelation_HLTRecObj = dbe->
book2D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_,Ptbins_,PtMin_,PtMax_);
3950 PtCorrelation_HLTRecObj->
getTH1();
3952 histoname = labelname+
"_hltRecObjPhiCorrelation";
3953 title = labelname+
"_hltRecObjPhiCorrelation;#Phi(HLT);#Phi(Reco)"+trigPath;
3954 MonitorElement * PhiCorrelation_HLTRecObj = dbe->
book2D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_,Phibins_,PhiMin_,PhiMax_);
3955 PhiCorrelation_HLTRecObj->
getTH1();
3983 PhiResolution_L1HLT,
3984 PtResolution_HLTRecObj,
3986 PhiResolution_HLTRecObj,
3987 PtCorrelation_L1HLT,
3989 PhiCorrelation_L1HLT,
3990 PtCorrelation_HLTRecObj,
3992 PhiCorrelation_HLTRecObj,
4012 double PtMax_ = 1000.;
4015 double EtaMin_ = -5.;
4016 double EtaMax_ = 5.;
4019 double PhiMin_ = -3.5;
4020 double PhiMax_ = 3.5;
4029 std::string subdirName = dirName1 + trgPathName +
"_wrt_" + trgPathNameD;
4037 if((
v->getObjectType() ==
trigger::TriggerJet) && (
v->getTriggerType().compare(
"SingleJet_Trigger") == 0)){
4039 histoname = labelname+
"_NumeratorPt";
4040 title = labelname+
"NumeratorPt;Calo Pt[GeV/c]";
4041 MonitorElement * NumeratorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4044 histoname = labelname+
"_NumeratorPtBarrel";
4045 title = labelname+
"NumeratorPtBarrel;Calo Pt[GeV/c] ";
4046 MonitorElement * NumeratorPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4047 NumeratorPtBarrel->
getTH1();
4049 histoname = labelname+
"_NumeratorPtEndcap";
4050 title = labelname+
"NumeratorPtEndcap;Calo Pt[GeV/c]";
4051 MonitorElement * NumeratorPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4052 NumeratorPtEndcap->
getTH1();
4054 histoname = labelname+
"_NumeratorPtForward";
4055 title = labelname+
"NumeratorPtForward;Calo Pt[GeV/c]";
4056 MonitorElement * NumeratorPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4057 NumeratorPtForward->
getTH1();
4059 histoname = labelname+
"_NumeratorEta";
4060 title = labelname+
"NumeratorEta;Calo #eta ";
4061 MonitorElement * NumeratorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4064 histoname = labelname+
"_NumeratorPhi";
4065 title = labelname+
"NumeratorPhi;Calo #Phi";
4066 MonitorElement * NumeratorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4069 histoname = labelname+
"_NumeratorEtaPhi";
4070 title = labelname+
"NumeratorEtaPhi;Calo #eta;Calo #Phi";
4071 MonitorElement * NumeratorEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4072 NumeratorEtaPhi->
getTH1();
4074 histoname = labelname+
"_NumeratorEtaBarrel";
4075 title = labelname+
"NumeratorEtaBarrel;Calo #eta ";
4076 MonitorElement * NumeratorEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4077 NumeratorEtaBarrel->
getTH1();
4079 histoname = labelname+
"_NumeratorPhiBarrel";
4080 title = labelname+
"NumeratorPhiBarrel;Calo #Phi";
4081 MonitorElement * NumeratorPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4082 NumeratorPhiBarrel->
getTH1();
4084 histoname = labelname+
"_NumeratorEtaEndcap";
4085 title = labelname+
"NumeratorEtaEndcap;Calo #eta ";
4086 MonitorElement * NumeratorEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4087 NumeratorEtaEndcap->
getTH1();
4089 histoname = labelname+
"_NumeratorPhiEndcap";
4090 title = labelname+
"NumeratorPhiEndcap;Calo #Phi";
4091 MonitorElement * NumeratorPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4092 NumeratorPhiEndcap->
getTH1();
4094 histoname = labelname+
"_NumeratorEtaForward";
4095 title = labelname+
"NumeratorEtaForward;Calo #eta ";
4096 MonitorElement * NumeratorEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4097 NumeratorEtaForward->
getTH1();
4099 histoname = labelname+
"_NumeratorPhiForward";
4100 title = labelname+
"NumeratorPhiForward;Calo #Phi";
4101 MonitorElement * NumeratorPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4102 NumeratorPhiForward->
getTH1();
4104 histoname = labelname+
"_NumeratorEta_LowpTcut";
4105 title = labelname+
"NumeratorEta_LowpTcut;Calo #eta ";
4106 MonitorElement * NumeratorEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4107 NumeratorEta_LowpTcut->
getTH1();
4109 histoname = labelname+
"_NumeratorPhi_LowpTcut";
4110 title = labelname+
"NumeratorPhi_LowpTcut;Calo #Phi";
4111 MonitorElement * NumeratorPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4112 NumeratorPhi_LowpTcut->
getTH1();
4114 histoname = labelname+
"_NumeratorEtaPhi_LowpTcut";
4115 title = labelname+
"NumeratorEtaPhi_LowpTcut;Calo #eta;Calo #Phi";
4116 MonitorElement * NumeratorEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4117 NumeratorEtaPhi_LowpTcut->
getTH1();
4119 histoname = labelname+
"_NumeratorEta_MedpTcut";
4120 title = labelname+
"NumeratorEta_MedpTcut;Calo #eta ";
4121 MonitorElement * NumeratorEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4122 NumeratorEta_MedpTcut->
getTH1();
4124 histoname = labelname+
"_NumeratorPhi_MedpTcut";
4125 title = labelname+
"NumeratorPhi_MedpTcut;Calo #Phi";
4126 MonitorElement * NumeratorPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4127 NumeratorPhi_MedpTcut->
getTH1();
4129 histoname = labelname+
"_NumeratorEtaPhi_MedpTcut";
4130 title = labelname+
"NumeratorEtaPhi_MedpTcut;Calo #eta;Calo #Phi";
4131 MonitorElement * NumeratorEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4132 NumeratorEtaPhi_MedpTcut->
getTH1();
4134 histoname = labelname+
"_NumeratorEta_HighpTcut";
4135 title = labelname+
"NumeratorEta_HighpTcut;Calo #eta ";
4136 MonitorElement * NumeratorEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4137 NumeratorEta_HighpTcut->
getTH1();
4139 histoname = labelname+
"_NumeratorPhi_HighpTcut";
4140 title = labelname+
"NumeratorPhi_HighpTcut;Calo #Phi";
4141 MonitorElement * NumeratorPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4142 NumeratorPhi_HighpTcut->
getTH1();
4144 histoname = labelname+
"_NumeratorEtaPhi_HighpTcut";
4145 title = labelname+
"NumeratorEtaPhi_HighpTcut;Calo #eta;Calo #Phi";
4146 MonitorElement * NumeratorEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4147 NumeratorEtaPhi_HighpTcut->
getTH1();
4149 histoname = labelname+
"_DenominatorPt";
4150 title = labelname+
"DenominatorPt;Calo Pt[GeV/c]";
4151 MonitorElement * DenominatorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4154 histoname = labelname+
"_DenominatorPtBarrel";
4155 title = labelname+
"DenominatorPtBarrel;Calo Pt[GeV/c]";
4156 MonitorElement * DenominatorPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4157 DenominatorPtBarrel->
getTH1();
4159 histoname = labelname+
"_DenominatorPtEndcap";
4160 title = labelname+
"DenominatorPtEndcap;Calo Pt[GeV/c]";
4161 MonitorElement * DenominatorPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4162 DenominatorPtEndcap->
getTH1();
4164 histoname = labelname+
"_DenominatorPtForward";
4165 title = labelname+
"DenominatorPtForward;Calo Pt[GeV/c] ";
4166 MonitorElement * DenominatorPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4167 DenominatorPtForward->
getTH1();
4169 histoname = labelname+
"_DenominatorEta";
4170 title = labelname+
"DenominatorEta;Calo #eta ";
4171 MonitorElement * DenominatorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4172 DenominatorEta->
getTH1();
4174 histoname = labelname+
"_DenominatorPhi";
4175 title = labelname+
"DenominatorPhi;Calo #Phi";
4176 MonitorElement * DenominatorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4177 DenominatorPhi->
getTH1();
4179 histoname = labelname+
"_DenominatorEtaPhi";
4180 title = labelname+
"DenominatorEtaPhi;Calo #eta; Calo #Phi";
4181 MonitorElement * DenominatorEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4182 DenominatorEtaPhi->
getTH1();
4184 histoname = labelname+
"_DenominatorEtaBarrel";
4185 title = labelname+
"DenominatorEtaBarrel;Calo #eta ";
4186 MonitorElement * DenominatorEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4187 DenominatorEtaBarrel->
getTH1();
4189 histoname = labelname+
"_DenominatorPhiBarrel";
4190 title = labelname+
"DenominatorPhiBarrel;Calo #Phi";
4191 MonitorElement * DenominatorPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4192 DenominatorPhiBarrel->
getTH1();
4194 histoname = labelname+
"_DenominatorEtaEndcap";
4195 title = labelname+
"DenominatorEtaEndcap;Calo #eta ";
4196 MonitorElement * DenominatorEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4197 DenominatorEtaEndcap->
getTH1();
4199 histoname = labelname+
"_DenominatorPhiEndcap";
4200 title = labelname+
"DenominatorPhiEndcap;Calo #Phi";
4201 MonitorElement * DenominatorPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4202 DenominatorPhiEndcap->
getTH1();
4204 histoname = labelname+
"_DenominatorEtaForward";
4205 title = labelname+
"DenominatorEtaForward;Calo #eta ";
4206 MonitorElement * DenominatorEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4207 DenominatorEtaForward->
getTH1();
4209 histoname = labelname+
"_DenominatorPhiForward";
4210 title = labelname+
"DenominatorPhiForward;Calo #Phi";
4211 MonitorElement * DenominatorPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4212 DenominatorPhiForward->
getTH1();
4214 histoname = labelname+
"_DenominatorEta_LowpTcut";
4215 title = labelname+
"DenominatorEta_LowpTcut;Calo #eta ";
4216 MonitorElement * DenominatorEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4217 DenominatorEta_LowpTcut->
getTH1();
4219 histoname = labelname+
"_DenominatorPhi_LowpTcut";
4220 title = labelname+
"DenominatorPhi_LowpTcut;Calo #Phi";
4221 MonitorElement * DenominatorPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4222 DenominatorPhi_LowpTcut->
getTH1();
4224 histoname = labelname+
"_DenominatorEtaPhi_LowpTcut";
4225 title = labelname+
"DenominatorEtaPhi_LowpTcut;Calo #eta;Calo #Phi";
4226 MonitorElement * DenominatorEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4227 DenominatorEtaPhi_LowpTcut->
getTH1();
4229 histoname = labelname+
"_DenominatorEta_MedpTcut";
4230 title = labelname+
"DenominatorEta_MedpTcut;Calo #eta ";
4231 MonitorElement * DenominatorEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4232 DenominatorEta_MedpTcut->
getTH1();
4234 histoname = labelname+
"_DenominatorPhi_MedpTcut";
4235 title = labelname+
"DenominatorPhi_MedpTcut;Calo #Phi";
4236 MonitorElement * DenominatorPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4237 DenominatorPhi_MedpTcut->
getTH1();
4239 histoname = labelname+
"_DenominatorEtaPhi_MedpTcut";
4240 title = labelname+
"DenominatorEtaPhi_MedpTcut;Calo #eta;Calo #Phi";
4241 MonitorElement * DenominatorEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4242 DenominatorEtaPhi_MedpTcut->
getTH1();
4244 histoname = labelname+
"_DenominatorEta_HighpTcut";
4245 title = labelname+
"DenominatorEta_HighpTcut;Calo #eta ";
4246 MonitorElement * DenominatorEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4247 DenominatorEta_HighpTcut->
getTH1();
4249 histoname = labelname+
"_DenominatorPhi_HighpTcut";
4250 title = labelname+
"DenominatorPhi_HighpTcut;Calo #Phi";
4251 MonitorElement * DenominatorPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4252 DenominatorPhi_HighpTcut->
getTH1();
4254 histoname = labelname+
"_DenominatorEtaPhi_HighpTcut";
4255 title = labelname+
"DenominatorEtaPhi_HighpTcut;Calo #eta;Calo #Phi";
4256 MonitorElement * DenominatorEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4257 DenominatorEtaPhi_HighpTcut->
getTH1();
4259 histoname = labelname+
"_DeltaR";
4260 title = labelname+
"DeltaR;";
4264 histoname = labelname+
"_DeltaPhi";
4265 title = labelname+
"DeltaPhi;";
4269 histoname = labelname+
"_NumeratorPFMHT";
4270 title = labelname+
"NumeratorPFMHT;PFMHT[GeV/c]";
4271 MonitorElement * NumeratorPFMHT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4272 NumeratorPFMHT->
getTH1();
4274 histoname = labelname+
"_NumeratorPFPt";
4275 title = labelname+
"NumeratorPFPt;PF Pt[GeV/c]";
4276 MonitorElement * NumeratorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4279 histoname = labelname+
"_NumeratorPFPtBarrel";
4280 title = labelname+
"NumeratorPFPtBarrel;PF Pt[GeV/c] ";
4281 MonitorElement * NumeratorPFPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4282 NumeratorPFPtBarrel->
getTH1();
4284 histoname = labelname+
"_NumeratorPFPtEndcap";
4285 title = labelname+
"NumeratorPFPtEndcap;PF Pt[GeV/c]";
4286 MonitorElement * NumeratorPFPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4287 NumeratorPFPtEndcap->
getTH1();
4289 histoname = labelname+
"_NumeratorPFPtForward";
4290 title = labelname+
"NumeratorPFPtForward;PF Pt[GeV/c]";
4291 MonitorElement * NumeratorPFPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4292 NumeratorPFPtForward->
getTH1();
4294 histoname = labelname+
"_NumeratorPFEta";
4295 title = labelname+
"NumeratorPFEta;PF #eta ";
4296 MonitorElement * NumeratorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4297 NumeratorPFEta->
getTH1();
4299 histoname = labelname+
"_NumeratorPFPhi";
4300 title = labelname+
"NumeratorPFPhi;Calo #Phi";
4301 MonitorElement * NumeratorPFPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4302 NumeratorPFPhi->
getTH1();
4304 histoname = labelname+
"_NumeratorPFEtaPhi";
4305 title = labelname+
"NumeratorPFEtaPhi;PF #eta;Calo #Phi";
4306 MonitorElement * NumeratorPFEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4307 NumeratorPFEtaPhi->
getTH1();
4309 histoname = labelname+
"_NumeratorPFEtaBarrel";
4310 title = labelname+
"NumeratorPFEtaBarrel;PF #eta ";
4311 MonitorElement * NumeratorPFEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4312 NumeratorPFEtaBarrel->
getTH1();
4314 histoname = labelname+
"_NumeratorPFPhiBarrel";
4315 title = labelname+
"NumeratorPFPhiBarrel;PF #Phi";
4316 MonitorElement * NumeratorPFPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4317 NumeratorPFPhiBarrel->
getTH1();
4319 histoname = labelname+
"_NumeratorPFEtaEndcap";
4320 title = labelname+
"NumeratorPFEtaEndcap;Calo #eta ";
4321 MonitorElement * NumeratorPFEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4322 NumeratorPFEtaEndcap->
getTH1();
4324 histoname = labelname+
"_NumeratorPFPhiEndcap";
4325 title = labelname+
"NumeratorPFPhiEndcap;PF #Phi";
4326 MonitorElement * NumeratorPFPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4327 NumeratorPFPhiEndcap->
getTH1();
4329 histoname = labelname+
"_NumeratorPFEtaForward";
4330 title = labelname+
"NumeratorPFEtaForward;Calo #eta ";
4331 MonitorElement * NumeratorPFEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4332 NumeratorPFEtaForward->
getTH1();
4334 histoname = labelname+
"_NumeratorPFPhiForward";
4335 title = labelname+
"NumeratorPFPhiForward;PF #Phi";
4336 MonitorElement * NumeratorPFPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4337 NumeratorPFPhiForward->
getTH1();
4339 histoname = labelname+
"_NumeratorPFEta_LowpTcut";
4340 title = labelname+
"NumeratorPFEta_LowpTcut;PF #eta ";
4341 MonitorElement * NumeratorPFEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4342 NumeratorPFEta_LowpTcut->
getTH1();
4344 histoname = labelname+
"_NumeratorPFPhi_LowpTcut";
4345 title = labelname+
"NumeratorPFPhi_LowpTcut;PF #Phi";
4346 MonitorElement * NumeratorPFPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4347 NumeratorPFPhi_LowpTcut->
getTH1();
4349 histoname = labelname+
"_NumeratorPFEtaPhi_LowpTcut";
4350 title = labelname+
"NumeratorPFEtaPhi_LowpTcut;PF #eta;Calo #Phi";
4351 MonitorElement * NumeratorPFEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4352 NumeratorPFEtaPhi_LowpTcut->
getTH1();
4354 histoname = labelname+
"_NumeratorPFEta_MedpTcut";
4355 title = labelname+
"NumeratorPFEta_MedpTcut;PF #eta ";
4356 MonitorElement * NumeratorPFEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4357 NumeratorPFEta_MedpTcut->
getTH1();
4359 histoname = labelname+
"_NumeratorPFPhi_MedpTcut";
4360 title = labelname+
"NumeratorPFPhi_MedpTcut;PF #Phi";
4361 MonitorElement * NumeratorPFPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4362 NumeratorPFPhi_MedpTcut->
getTH1();
4364 histoname = labelname+
"_NumeratorPFEtaPhi_MedpTcut";
4365 title = labelname+
"NumeratorPFEtaPhi_MedpTcut;PF #eta;PF #Phi";
4366 MonitorElement * NumeratorPFEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4367 NumeratorPFEtaPhi_MedpTcut->
getTH1();
4369 histoname = labelname+
"_NumeratorPFEta_HighpTcut";
4370 title = labelname+
"NumeratorPFEta_HighpTcut;Calo #eta ";
4371 MonitorElement * NumeratorPFEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4372 NumeratorPFEta_HighpTcut->
getTH1();
4374 histoname = labelname+
"_NumeratorPFPhi_HighpTcut";
4375 title = labelname+
"NumeratorPFPhi_HighpTcut;PF #Phi";
4376 MonitorElement * NumeratorPFPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4377 NumeratorPFPhi_HighpTcut->
getTH1();
4379 histoname = labelname+
"_NumeratorPFEtaPhi_HighpTcut";
4380 title = labelname+
"NumeratorPFEtaPhi_HighpTcut;PF #eta;PF #Phi";
4381 MonitorElement * NumeratorPFEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4382 NumeratorPFEtaPhi_HighpTcut->
getTH1();
4384 histoname = labelname+
"_DenominatorPFMHT";
4385 title = labelname+
"DenominatorPFMHT;PF Pt[GeV/c]";
4386 MonitorElement * DenominatorPFMHT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4387 DenominatorPFMHT->
getTH1();
4389 histoname = labelname+
"_DenominatorPFPt";
4390 title = labelname+
"DenominatorPFPt;PF Pt[GeV/c]";
4391 MonitorElement * DenominatorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4392 DenominatorPFPt->
getTH1();
4394 histoname = labelname+
"_DenominatorPFPtBarrel";
4395 title = labelname+
"DenominatorPFPtBarrel;Calo Pt[GeV/c]";
4396 MonitorElement * DenominatorPFPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4397 DenominatorPFPtBarrel->
getTH1();
4399 histoname = labelname+
"_DenominatorPFPtEndcap";
4400 title = labelname+
"DenominatorPFPtEndcap;PF Pt[GeV/c]";
4401 MonitorElement * DenominatorPFPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4402 DenominatorPFPtEndcap->
getTH1();
4404 histoname = labelname+
"_DenominatorPFPtForward";
4405 title = labelname+
"DenominatorPFPtForward;PF Pt[GeV/c] ";
4406 MonitorElement * DenominatorPFPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4407 DenominatorPFPtForward->
getTH1();
4409 histoname = labelname+
"_DenominatorPFEta";
4410 title = labelname+
"DenominatorPFEta;PF #eta ";
4411 MonitorElement * DenominatorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4412 DenominatorPFEta->
getTH1();
4414 histoname = labelname+
"_DenominatorPFPhi";
4415 title = labelname+
"DenominatorPFPhi;PF #Phi";
4416 MonitorElement * DenominatorPFPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4417 DenominatorPFPhi->
getTH1();
4419 histoname = labelname+
"_DenominatorPFEtaPhi";
4420 title = labelname+
"DenominatorPFEtaPhi;PF #eta; Calo #Phi";
4421 MonitorElement * DenominatorPFEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4422 DenominatorPFEtaPhi->
getTH1();
4424 histoname = labelname+
"_DenominatorPFEtaBarrel";
4425 title = labelname+
"DenominatorPFEtaBarrel;Calo #eta ";
4426 MonitorElement * DenominatorPFEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4427 DenominatorPFEtaBarrel->
getTH1();
4429 histoname = labelname+
"_DenominatorPFPhiBarrel";
4430 title = labelname+
"DenominatorPFPhiBarrel;PF #Phi";
4431 MonitorElement * DenominatorPFPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4432 DenominatorPFPhiBarrel->
getTH1();
4434 histoname = labelname+
"_DenominatorPFEtaEndcap";
4435 title = labelname+
"DenominatorPFEtaEndcap;PF #eta ";
4436 MonitorElement * DenominatorPFEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4437 DenominatorPFEtaEndcap->
getTH1();
4439 histoname = labelname+
"_DenominatorPFPhiEndcap";
4440 title = labelname+
"DenominatorPFPhiEndcap;Calo #Phi";
4441 MonitorElement * DenominatorPFPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4442 DenominatorPFPhiEndcap->
getTH1();
4444 histoname = labelname+
"_DenominatorPFEtaForward";
4445 title = labelname+
"DenominatorPFEtaForward;PF #eta ";
4446 MonitorElement * DenominatorPFEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4447 DenominatorPFEtaForward->
getTH1();
4449 histoname = labelname+
"_DenominatorPFPhiForward";
4450 title = labelname+
"DenominatorPFPhiForward;PF #Phi";
4451 MonitorElement * DenominatorPFPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4452 DenominatorPFPhiForward->
getTH1();
4454 histoname = labelname+
"_DenominatorPFEta_LowpTcut";
4455 title = labelname+
"DenominatorPFEta_LowpTcut;PF #eta ";
4456 MonitorElement * DenominatorPFEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4457 DenominatorPFEta_LowpTcut->
getTH1();
4459 histoname = labelname+
"_DenominatorPFPhi_LowpTcut";
4460 title = labelname+
"DenominatorPFPhi_LowpTcut;PF #Phi";
4461 MonitorElement * DenominatorPFPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4462 DenominatorPFPhi_LowpTcut->
getTH1();
4464 histoname = labelname+
"_DenominatorPFEtaPhi_LowpTcut";
4465 title = labelname+
"DenominatorPFEtaPhi_LowpTcut;PF #eta;Calo #Phi";
4466 MonitorElement * DenominatorPFEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4467 DenominatorPFEtaPhi_LowpTcut->
getTH1();
4469 histoname = labelname+
"_DenominatorPFEta_MedpTcut";
4470 title = labelname+
"DenominatorPFEta_MedpTcut;PF #eta ";
4471 MonitorElement * DenominatorPFEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4472 DenominatorPFEta_MedpTcut->
getTH1();
4474 histoname = labelname+
"_DenominatorPFPhi_MedpTcut";
4475 title = labelname+
"DenominatorPFPhi_MedpTcut;PF #Phi";
4476 MonitorElement * DenominatorPFPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4477 DenominatorPFPhi_MedpTcut->
getTH1();
4479 histoname = labelname+
"_DenominatorPFEtaPhi_MedpTcut";
4480 title = labelname+
"DenominatorPFEtaPhi_MedpTcut;PF #eta;Calo #Phi";
4481 MonitorElement * DenominatorPFEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4482 DenominatorPFEtaPhi_MedpTcut->
getTH1();
4484 histoname = labelname+
"_DenominatorPFEta_HighpTcut";
4485 title = labelname+
"DenominatorPFEta_HighpTcut;PF #eta ";
4486 MonitorElement * DenominatorPFEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4487 DenominatorPFEta_HighpTcut->
getTH1();
4489 histoname = labelname+
"_DenominatorPFPhi_HighpTcut";
4490 title = labelname+
"DenominatorPFPhi_HighpTcut;PF #Phi";
4491 MonitorElement * DenominatorPFPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4492 DenominatorPFPhi_HighpTcut->
getTH1();
4494 histoname = labelname+
"_DenominatorPFEtaPhi_HighpTcut";
4495 title = labelname+
"DenominatorPFEtaPhi_HighpTcut;PF #eta;Calo #Phi";
4496 MonitorElement * DenominatorPFEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4497 DenominatorPFEtaPhi_HighpTcut->
getTH1();
4499 histoname = labelname+
"_PFDeltaR";
4500 title = labelname+
"PFDeltaR;";
4504 histoname = labelname+
"_PFDeltaPhi";
4505 title = labelname+
"PFDeltaPhi;";
4509 v->setEffHistos(NumeratorPt,
4521 NumeratorEtaForward,
4522 NumeratorPhiForward,
4523 NumeratorEta_LowpTcut,
4524 NumeratorPhi_LowpTcut,
4525 NumeratorEtaPhi_LowpTcut,
4526 NumeratorEta_MedpTcut,
4527 NumeratorPhi_MedpTcut,
4528 NumeratorEtaPhi_MedpTcut,
4529 NumeratorEta_HighpTcut,
4530 NumeratorPhi_HighpTcut,
4531 NumeratorEtaPhi_HighpTcut,
4534 DenominatorPtBarrel,
4535 DenominatorPtEndcap,
4536 DenominatorPtForward,
4541 DenominatorEtaBarrel,
4542 DenominatorPhiBarrel,
4543 DenominatorEtaEndcap,
4544 DenominatorPhiEndcap,
4545 DenominatorEtaForward,
4546 DenominatorPhiForward,
4547 DenominatorEta_LowpTcut,
4548 DenominatorPhi_LowpTcut,
4549 DenominatorEtaPhi_LowpTcut,
4550 DenominatorEta_MedpTcut,
4551 DenominatorPhi_MedpTcut,
4552 DenominatorEtaPhi_MedpTcut,
4553 DenominatorEta_HighpTcut,
4554 DenominatorPhi_HighpTcut,
4555 DenominatorEtaPhi_HighpTcut,
4561 NumeratorPFPtBarrel,
4562 NumeratorPFPtEndcap,
4563 NumeratorPFPtForward,
4567 NumeratorPFEtaBarrel,
4568 NumeratorPFPhiBarrel,
4569 NumeratorPFEtaEndcap,
4570 NumeratorPFPhiEndcap,
4571 NumeratorPFEtaForward,
4572 NumeratorPFPhiForward,
4573 NumeratorPFEta_LowpTcut,
4574 NumeratorPFPhi_LowpTcut,
4575 NumeratorPFEtaPhi_LowpTcut,
4576 NumeratorPFEta_MedpTcut,
4577 NumeratorPFPhi_MedpTcut,
4578 NumeratorPFEtaPhi_MedpTcut,
4579 NumeratorPFEta_HighpTcut,
4580 NumeratorPFPhi_HighpTcut,
4581 NumeratorPFEtaPhi_HighpTcut,
4584 DenominatorPFPtBarrel,
4585 DenominatorPFPtEndcap,
4586 DenominatorPFPtForward,
4589 DenominatorPFEtaPhi,
4590 DenominatorPFEtaBarrel,
4591 DenominatorPFPhiBarrel,
4592 DenominatorPFEtaEndcap,
4593 DenominatorPFPhiEndcap,
4594 DenominatorPFEtaForward,
4595 DenominatorPFPhiForward,
4596 DenominatorPFEta_LowpTcut,
4597 DenominatorPFPhi_LowpTcut,
4598 DenominatorPFEtaPhi_LowpTcut,
4599 DenominatorPFEta_MedpTcut,
4600 DenominatorPFPhi_MedpTcut,
4601 DenominatorPFEtaPhi_MedpTcut,
4602 DenominatorPFEta_HighpTcut,
4603 DenominatorPFPhi_HighpTcut,
4604 DenominatorPFEtaPhi_HighpTcut,
4611 if((
v->getObjectType() ==
trigger::TriggerJet) && (
v->getTriggerType().compare(
"DiJet_Trigger") == 0)){
4613 histoname = labelname+
"_NumeratorAvrgPt";
4614 title = labelname+
"NumeratorAvrgPt;Calo Pt[GeV/c]";
4615 MonitorElement * NumeratorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4618 histoname = labelname+
"_NumeratorAvrgEta";
4619 title = labelname+
"NumeratorAvrgEta;Calo #eta";
4620 MonitorElement * NumeratorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4623 histoname = labelname+
"_DenominatorAvrgPt";
4624 title = labelname+
"DenominatorAvrgPt;Calo Pt[GeV/c] ";
4625 MonitorElement * DenominatorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4628 histoname = labelname+
"_DenominatorAvrgEta";
4629 title = labelname+
"DenominatorAvrgEta;Calo #eta";
4630 MonitorElement * DenominatorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4631 DenominatorEta->
getTH1();
4633 histoname = labelname+
"_DeltaR";
4634 title = labelname+
"DeltaR;";
4638 histoname = labelname+
"_DeltaPhi";
4639 title = labelname+
"DeltaPhi;";
4644 histoname = labelname+
"_NumeratorAvrgPFPt";
4645 title = labelname+
"NumeratorAvrgPFPt;PF Pt[GeV/c]";
4646 MonitorElement * NumeratorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4649 histoname = labelname+
"_NumeratorAvrgPFEta";
4650 title = labelname+
"NumeratorAvrgPFEta;PF #eta";
4651 MonitorElement * NumeratorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4652 NumeratorPFEta->
getTH1();
4654 histoname = labelname+
"_DenominatorAvrgPFPt";
4655 title = labelname+
"DenominatorAvrgPFPt;PF Pt[GeV/c] ";
4656 MonitorElement * DenominatorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4657 DenominatorPFPt->
getTH1();
4659 histoname = labelname+
"_DenominatorAvrgPFEta";
4660 title = labelname+
"DenominatorAvrgPFEta;PF #eta";
4661 MonitorElement * DenominatorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4662 DenominatorPFEta->
getTH1();
4664 histoname = labelname+
"_PFDeltaR";
4665 title = labelname+
"PFDeltaR;";
4669 histoname = labelname+
"_PFDeltaPhi";
4670 title = labelname+
"PFDeltaPhi;";
4674 v->setEffHistos( dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4675 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4676 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4677 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4678 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4679 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4680 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4681 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4682 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4683 dummy, dummy, dummy, dummy
4689 histoname = labelname+
"_NumeratorPt";
4690 if(
v->getPath().find(
"HLT_PFMET")==std::string::npos)
4691 title = labelname+
"NumeratorPt; CaloMET[GeV/c]";
4693 title = labelname+
"NumeratorPt; PFMET[GeV/c]";
4694 MonitorElement * NumeratorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4697 histoname = labelname+
"_NumeratorPhi";
4698 title = labelname+
"NumeratorPhi; #Phi";
4699 MonitorElement * NumeratorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4702 histoname = labelname+
"_DenominatorPt";
4703 if(
v->getPath().find(
"HLT_PFMET")==std::string::npos)
4704 title = labelname+
"DenominatorPt; CaloMET[GeV/c]";
4706 title = labelname+
"DenominatorPt; PFMET[GeV/c]";
4707 MonitorElement * DenominatorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4710 histoname = labelname+
"_DenominatorPhi";
4711 title = labelname+
"DenominatorPhi; #Phi";
4712 MonitorElement * DenominatorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4713 DenominatorPhi->
getTH1();
4715 v->setEffHistos( NumeratorPt, dummy, dummy, dummy, dummy, NumeratorPhi, dummy, dummy, dummy, dummy,
4716 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4717 dummy, dummy, DenominatorPt, dummy, dummy, dummy, dummy, DenominatorPhi, dummy, dummy,
4718 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4719 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4720 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4721 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4722 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4723 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4724 dummy, dummy, dummy, dummy
4745 if((
v->getObjectType() ==
trigger::TriggerJet) && (
v->getTriggerType().compare(
"SingleJet_Trigger") == 0)){
4747 histoname = labelname+
"_NumeratorPt";
4748 title = labelname+
"NumeratorPt;Pt[GeV/c]";
4749 MonitorElement * NumeratorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4752 histoname = labelname+
"_NumeratorPtBarrel";
4753 title = labelname+
"NumeratorPtBarrel;Calo Pt[GeV/c]";
4754 MonitorElement * NumeratorPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4755 NumeratorPtBarrel->
getTH1();
4757 histoname = labelname+
"_NumeratorPtEndcap";
4758 title = labelname+
"NumeratorPtEndcap;Calo Pt[GeV/c]";
4759 MonitorElement * NumeratorPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4760 NumeratorPtEndcap->
getTH1();
4762 histoname = labelname+
"_NumeratorPtForward";
4763 title = labelname+
"NumeratorPtForward;Calo Pt[GeV/c]";
4764 MonitorElement * NumeratorPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4765 NumeratorPtForward->
getTH1();
4767 histoname = labelname+
"_NumeratorEta";
4768 title = labelname+
"NumeratorEta;Calo #eta ";
4769 MonitorElement * NumeratorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4772 histoname = labelname+
"_NumeratorPhi";
4773 title = labelname+
"NumeratorPhi;Calo #Phi";
4774 MonitorElement * NumeratorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4777 histoname = labelname+
"_NumeratorEtaPhi";
4778 title = labelname+
"NumeratorEtaPhi;Calo #eta;Calo #Phi";
4779 MonitorElement * NumeratorEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4780 NumeratorEtaPhi->
getTH1();
4782 histoname = labelname+
"_NumeratorEtaBarrel";
4783 title = labelname+
"NumeratorEtaBarrel;Calo #eta ";
4784 MonitorElement * NumeratorEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4785 NumeratorEtaBarrel->
getTH1();
4787 histoname = labelname+
"_NumeratorPhiBarrel";
4788 title = labelname+
"NumeratorPhiBarrel;Calo #Phi";
4789 MonitorElement * NumeratorPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4790 NumeratorPhiBarrel->
getTH1();
4792 histoname = labelname+
"_NumeratorEtaEndcap";
4793 title = labelname+
"NumeratorEtaEndcap;Calo #eta ";
4794 MonitorElement * NumeratorEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4795 NumeratorEtaEndcap->
getTH1();
4797 histoname = labelname+
"_NumeratorPhiEndcap";
4798 title = labelname+
"NumeratorPhiEndcap;Calo #Phi";
4799 MonitorElement * NumeratorPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4800 NumeratorPhiEndcap->
getTH1();
4802 histoname = labelname+
"_NumeratorEtaForward";
4803 title = labelname+
"NumeratorEtaForward;Calo #eta ";
4804 MonitorElement * NumeratorEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4805 NumeratorEtaForward->
getTH1();
4807 histoname = labelname+
"_NumeratorPhiForward";
4808 title = labelname+
"NumeratorPhiForward;Calo #Phi";
4809 MonitorElement * NumeratorPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4810 NumeratorPhiForward->
getTH1();
4812 histoname = labelname+
"_NumeratorEta_LowpTcut";
4813 title = labelname+
"NumeratorEta_LowpTcut;Calo #eta ";
4814 MonitorElement * NumeratorEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4815 NumeratorEta_LowpTcut->
getTH1();
4817 histoname = labelname+
"_NumeratorPhi_LowpTcut";
4818 title = labelname+
"NumeratorPhi_LowpTcut;Calo #Phi";
4819 MonitorElement * NumeratorPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4820 NumeratorPhi_LowpTcut->
getTH1();
4822 histoname = labelname+
"_NumeratorEtaPhi_LowpTcut";
4823 title = labelname+
"NumeratorEtaPhi_LowpTcut;Calo #eta;Calo #Phi";
4824 MonitorElement * NumeratorEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4825 NumeratorEtaPhi_LowpTcut->
getTH1();
4827 histoname = labelname+
"_NumeratorEta_MedpTcut";
4828 title = labelname+
"NumeratorEta_MedpTcut;Calo #eta ";
4829 MonitorElement * NumeratorEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4830 NumeratorEta_MedpTcut->
getTH1();
4832 histoname = labelname+
"_NumeratorPhi_MedpTcut";
4833 title = labelname+
"NumeratorPhi_MedpTcut;Calo #Phi";
4834 MonitorElement * NumeratorPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4835 NumeratorPhi_MedpTcut->
getTH1();
4837 histoname = labelname+
"_NumeratorEtaPhi_MedpTcut";
4838 title = labelname+
"NumeratorEtaPhi_MedpTcut;Calo #eta;Calo #Phi";
4839 MonitorElement * NumeratorEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4840 NumeratorEtaPhi_MedpTcut->
getTH1();
4842 histoname = labelname+
"_NumeratorEta_HighpTcut";
4843 title = labelname+
"NumeratorEta_HighpTcut;Calo #eta ";
4844 MonitorElement * NumeratorEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4845 NumeratorEta_HighpTcut->
getTH1();
4847 histoname = labelname+
"_NumeratorPhi_HighpTcut";
4848 title = labelname+
"NumeratorPhi_HighpTcut;Calo #Phi";
4849 MonitorElement * NumeratorPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4850 NumeratorPhi_HighpTcut->
getTH1();
4852 histoname = labelname+
"_NumeratorEtaPhi_HighpTcut";
4853 title = labelname+
"NumeratorEtaPhi_HighpTcut;Calo #eta;Calo #Phi";
4854 MonitorElement * NumeratorEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4855 NumeratorEtaPhi_HighpTcut->
getTH1();
4857 histoname = labelname+
"_DenominatorPt";
4858 title = labelname+
"DenominatorPt;Calo Pt[GeV/c]";
4859 MonitorElement * DenominatorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4862 histoname = labelname+
"_DenominatorPtBarrel";
4863 title = labelname+
"DenominatorPtBarrel;Calo Pt[GeV/c]";
4864 MonitorElement * DenominatorPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4865 DenominatorPtBarrel->
getTH1();
4867 histoname = labelname+
"_DenominatorPtEndcap";
4868 title = labelname+
"DenominatorPtEndcap;Calo Pt[GeV/c]";
4869 MonitorElement * DenominatorPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4870 DenominatorPtEndcap->
getTH1();
4872 histoname = labelname+
"_DenominatorPtForward";
4873 title = labelname+
"DenominatorPtForward;Calo Pt[GeV/c] ";
4874 MonitorElement * DenominatorPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4875 DenominatorPtForward->
getTH1();
4877 histoname = labelname+
"_DenominatorEta";
4878 title = labelname+
"DenominatorEta;Calo #eta";
4879 MonitorElement * DenominatorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4880 DenominatorEta->
getTH1();
4882 histoname = labelname+
"_DenominatorPhi";
4883 title = labelname+
"DenominatorPhi;Calo #Phi";
4884 MonitorElement * DenominatorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4885 DenominatorPhi->
getTH1();
4887 histoname = labelname+
"_DenominatorEtaPhi";
4888 title = labelname+
"DenominatorEtaPhi;Calo #eta (IC5);Calo #Phi ";
4889 MonitorElement * DenominatorEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4890 DenominatorEtaPhi->
getTH1();
4892 histoname = labelname+
"_DenominatorEtaBarrel";
4893 title = labelname+
"DenominatorEtaBarrel;Calo #eta ";
4894 MonitorElement * DenominatorEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4895 DenominatorEtaBarrel->
getTH1();
4897 histoname = labelname+
"_DenominatorPhiBarrel";
4898 title = labelname+
"DenominatorPhiBarrel;Calo #Phi";
4899 MonitorElement * DenominatorPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4900 DenominatorPhiBarrel->
getTH1();
4902 histoname = labelname+
"_DenominatorEtaEndcap";
4903 title = labelname+
"DenominatorEtaEndcap;Calo #eta ";
4904 MonitorElement * DenominatorEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4905 DenominatorEtaEndcap->
getTH1();
4907 histoname = labelname+
"_DenominatorPhiEndcap";
4908 title = labelname+
"DenominatorPhiEndcap;Calo #Phi";
4909 MonitorElement * DenominatorPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4910 DenominatorPhiEndcap->
getTH1();
4912 histoname = labelname+
"_DenominatorEtaForward";
4913 title = labelname+
"DenominatorEtaForward;Calo #eta ";
4914 MonitorElement * DenominatorEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4915 DenominatorEtaForward->
getTH1();
4917 histoname = labelname+
"_DenominatorPhiForward";
4918 title = labelname+
"DenominatorPhiForward;Calo #Phi";
4919 MonitorElement * DenominatorPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4920 DenominatorPhiForward->
getTH1();
4922 histoname = labelname+
"_DenominatorEta_LowpTcut";
4923 title = labelname+
"DenominatorEta_LowpTcut;Calo #eta ";
4924 MonitorElement * DenominatorEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4925 DenominatorEta_LowpTcut->
getTH1();
4927 histoname = labelname+
"_DenominatorPhi_LowpTcut";
4928 title = labelname+
"DenominatorPhi_LowpTcut;Calo #Phi";
4929 MonitorElement * DenominatorPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4930 DenominatorPhi_LowpTcut->
getTH1();
4932 histoname = labelname+
"_DenominatorEtaPhi_LowpTcut";
4933 title = labelname+
"DenominatorEtaPhi_LowpTcut;Calo #eta;Calo #Phi";
4934 MonitorElement * DenominatorEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4935 DenominatorEtaPhi_LowpTcut->
getTH1();
4937 histoname = labelname+
"_DenominatorEta_MedpTcut";
4938 title = labelname+
"DenominatorEta_MedpTcut;Calo #eta ";
4939 MonitorElement * DenominatorEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4940 DenominatorEta_MedpTcut->
getTH1();
4942 histoname = labelname+
"_DenominatorPhi_MedpTcut";
4943 title = labelname+
"DenominatorPhi_MedpTcut;Calo #Phi";
4944 MonitorElement * DenominatorPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4945 DenominatorPhi_MedpTcut->
getTH1();
4947 histoname = labelname+
"_DenominatorEtaPhi_MedpTcut";
4948 title = labelname+
"DenominatorEtaPhi_MedpTcut;Calo #eta;Calo #Phi";
4949 MonitorElement * DenominatorEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4950 DenominatorEtaPhi_MedpTcut->
getTH1();
4952 histoname = labelname+
"_DenominatorEta_HighpTcut";
4953 title = labelname+
"DenominatorEta_HighpTcut;Calo #eta ";
4954 MonitorElement * DenominatorEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4955 DenominatorEta_HighpTcut->
getTH1();
4957 histoname = labelname+
"_DenominatorPhi_HighpTcut";
4958 title = labelname+
"DenominatorPhi_HighpTcut;Calo #Phi";
4959 MonitorElement * DenominatorPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4960 DenominatorPhi_HighpTcut->
getTH1();
4962 histoname = labelname+
"_DenominatorEtaPhi_HighpTcut";
4963 title = labelname+
"DenominatorEtaPhi_HighpTcut;Calo #eta;Calo #Phi";
4964 MonitorElement * DenominatorEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4965 DenominatorEtaPhi_HighpTcut->
getTH1();
4967 histoname = labelname+
"_DeltaR";
4968 title = labelname+
"DeltaR;";
4972 histoname = labelname+
"_DeltaPhi";
4973 title = labelname+
"DeltaPhi;";
4977 histoname = labelname+
"_NumeratorPFPt";
4978 title = labelname+
"NumeratorPFPt;PFPt[GeV/c]";
4979 MonitorElement * NumeratorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4982 histoname = labelname+
"_NumeratorPFMHT";
4983 title = labelname+
"NumeratorPFMHT;PFMHT[GeV/c]";
4984 MonitorElement * NumeratorPFMHT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4985 NumeratorPFMHT->
getTH1();
4987 histoname = labelname+
"_NumeratorPFPtBarrel";
4988 title = labelname+
"NumeratorPFPtBarrel;PF Pt[GeV/c]";
4989 MonitorElement * NumeratorPFPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4990 NumeratorPFPtBarrel->
getTH1();
4992 histoname = labelname+
"_NumeratorPFPtEndcap";
4993 title = labelname+
"NumeratorPFPtEndcap;PF Pt[GeV/c]";
4994 MonitorElement * NumeratorPFPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4995 NumeratorPFPtEndcap->
getTH1();
4997 histoname = labelname+
"_NumeratorPFPtForward";
4998 title = labelname+
"NumeratorPFPtForward;PF Pt[GeV/c]";
4999 MonitorElement * NumeratorPFPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5000 NumeratorPFPtForward->
getTH1();
5002 histoname = labelname+
"_NumeratorPFEta";
5003 title = labelname+
"NumeratorPFEta;PF #eta ";
5004 MonitorElement * NumeratorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5005 NumeratorPFEta->
getTH1();
5007 histoname = labelname+
"_NumeratorPFPhi";
5008 title = labelname+
"NumeratorPFPhi;PF #Phi";
5009 MonitorElement * NumeratorPFPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5010 NumeratorPFPhi->
getTH1();
5012 histoname = labelname+
"_NumeratorPFEtaPhi";
5013 title = labelname+
"NumeratorPFEtaPhi;PF #eta;PF #Phi";
5014 MonitorElement * NumeratorPFEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5015 NumeratorPFEtaPhi->
getTH1();
5017 histoname = labelname+
"_NumeratorPFEtaBarrel";
5018 title = labelname+
"NumeratorPFEtaBarrel;PF #eta ";
5019 MonitorElement * NumeratorPFEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5020 NumeratorPFEtaBarrel->
getTH1();
5022 histoname = labelname+
"_NumeratorPFPhiBarrel";
5023 title = labelname+
"NumeratorPFPhiBarrel;PF #Phi";
5024 MonitorElement * NumeratorPFPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5025 NumeratorPFPhiBarrel->
getTH1();
5027 histoname = labelname+
"_NumeratorPFEtaEndcap";
5028 title = labelname+
"NumeratorPFEtaEndcap;PF #eta ";
5029 MonitorElement * NumeratorPFEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5030 NumeratorPFEtaEndcap->
getTH1();
5032 histoname = labelname+
"_NumeratorPFPhiEndcap";
5033 title = labelname+
"NumeratorPFPhiEndcap;PF #Phi";
5034 MonitorElement * NumeratorPFPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5035 NumeratorPFPhiEndcap->
getTH1();
5037 histoname = labelname+
"_NumeratorPFEtaForward";
5038 title = labelname+
"NumeratorPFEtaForward;PF #eta ";
5039 MonitorElement * NumeratorPFEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5040 NumeratorPFEtaForward->
getTH1();
5042 histoname = labelname+
"_NumeratorPFPhiForward";
5043 title = labelname+
"NumeratorPFPhiForward;PF #Phi";
5044 MonitorElement * NumeratorPFPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5045 NumeratorPFPhiForward->
getTH1();
5047 histoname = labelname+
"_NumeratorPFEta_LowpTcut";
5048 title = labelname+
"NumeratorPFEta_LowpTcut;PF #eta ";
5049 MonitorElement * NumeratorPFEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5050 NumeratorPFEta_LowpTcut->
getTH1();
5052 histoname = labelname+
"_NumeratorPFPhi_LowpTcut";
5053 title = labelname+
"NumeratorPFPhi_LowpTcut;PF #Phi";
5054 MonitorElement * NumeratorPFPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5055 NumeratorPFPhi_LowpTcut->
getTH1();
5057 histoname = labelname+
"_NumeratorPFEtaPhi_LowpTcut";
5058 title = labelname+
"NumeratorPFEtaPhi_LowpTcut;PF #eta;PF #Phi";
5059 MonitorElement * NumeratorPFEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5060 NumeratorPFEtaPhi_LowpTcut->
getTH1();
5062 histoname = labelname+
"_NumeratorPFEta_MedpTcut";
5063 title = labelname+
"NumeratorPFEta_MedpTcut;PF #eta ";
5064 MonitorElement * NumeratorPFEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5065 NumeratorPFEta_MedpTcut->
getTH1();
5067 histoname = labelname+
"_NumeratorPFPhi_MedpTcut";
5068 title = labelname+
"NumeratorPFPhi_MedpTcut;PF #Phi";
5069 MonitorElement * NumeratorPFPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5070 NumeratorPFPhi_MedpTcut->
getTH1();
5072 histoname = labelname+
"_NumeratorPFEtaPhi_MedpTcut";
5073 title = labelname+
"NumeratorPFEtaPhi_MedpTcut;PF #eta;PF #Phi";
5074 MonitorElement * NumeratorPFEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5075 NumeratorPFEtaPhi_MedpTcut->
getTH1();
5077 histoname = labelname+
"_NumeratorPFEta_HighpTcut";
5078 title = labelname+
"NumeratorPFEta_HighpTcut;PF #eta ";
5079 MonitorElement * NumeratorPFEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5080 NumeratorPFEta_HighpTcut->
getTH1();
5082 histoname = labelname+
"_NumeratorPFPhi_HighpTcut";
5083 title = labelname+
"NumeratorPFPhi_HighpTcut;PF #Phi";
5084 MonitorElement * NumeratorPFPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5085 NumeratorPFPhi_HighpTcut->
getTH1();
5087 histoname = labelname+
"_NumeratorPFEtaPhi_HighpTcut";
5088 title = labelname+
"NumeratorPFEtaPhi_HighpTcut;PF #eta;PF #Phi";
5089 MonitorElement * NumeratorPFEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5090 NumeratorPFEtaPhi_HighpTcut->
getTH1();
5092 histoname = labelname+
"_DenominatorPFMHT";
5093 title = labelname+
"DenominatorPFMHT;PF Pt[GeV/c]";
5094 MonitorElement * DenominatorPFMHT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5095 DenominatorPFMHT->
getTH1();
5097 histoname = labelname+
"_DenominatorPFPt";
5098 title = labelname+
"DenominatorPFPt;PF Pt[GeV/c]";
5099 MonitorElement * DenominatorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5100 DenominatorPFPt->
getTH1();
5102 histoname = labelname+
"_DenominatorPFPtBarrel";
5103 title = labelname+
"DenominatorPFPtBarrel;PF Pt[GeV/c]";
5104 MonitorElement * DenominatorPFPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5105 DenominatorPFPtBarrel->
getTH1();
5107 histoname = labelname+
"_DenominatorPFPtEndcap";
5108 title = labelname+
"DenominatorPFPtEndcap;PF Pt[GeV/c]";
5109 MonitorElement * DenominatorPFPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5110 DenominatorPFPtEndcap->
getTH1();
5112 histoname = labelname+
"_DenominatorPFPtForward";
5113 title = labelname+
"DenominatorPFPtForward;PF Pt[GeV/c] ";
5114 MonitorElement * DenominatorPFPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5115 DenominatorPFPtForward->
getTH1();
5117 histoname = labelname+
"_DenominatorPFEta";
5118 title = labelname+
"DenominatorPFEta;PF #eta";
5119 MonitorElement * DenominatorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5120 DenominatorPFEta->
getTH1();
5122 histoname = labelname+
"_DenominatorPFPhi";
5123 title = labelname+
"DenominatorPFPhi;PF #Phi";
5124 MonitorElement * DenominatorPFPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5125 DenominatorPFPhi->
getTH1();
5127 histoname = labelname+
"_DenominatorPFEtaPhi";
5128 title = labelname+
"DenominatorPFEtaPhi;PF #eta (IC5);PF #Phi ";
5129 MonitorElement * DenominatorPFEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5130 DenominatorPFEtaPhi->
getTH1();
5132 histoname = labelname+
"_DenominatorPFEtaBarrel";
5133 title = labelname+
"DenominatorPFEtaBarrel;PF #eta ";
5134 MonitorElement * DenominatorPFEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5135 DenominatorPFEtaBarrel->
getTH1();
5137 histoname = labelname+
"_DenominatorPFPhiBarrel";
5138 title = labelname+
"DenominatorPFPhiBarrel;PF #Phi";
5139 MonitorElement * DenominatorPFPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5140 DenominatorPFPhiBarrel->
getTH1();
5142 histoname = labelname+
"_DenominatorPFEtaEndcap";
5143 title = labelname+
"DenominatorPFEtaEndcap;PF #eta ";
5144 MonitorElement * DenominatorPFEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5145 DenominatorPFEtaEndcap->
getTH1();
5147 histoname = labelname+
"_DenominatorPFPhiEndcap";
5148 title = labelname+
"DenominatorPFPhiEndcap;PF #Phi";
5149 MonitorElement * DenominatorPFPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5150 DenominatorPFPhiEndcap->
getTH1();
5152 histoname = labelname+
"_DenominatorPFEtaForward";
5153 title = labelname+
"DenominatorPFEtaForward;PF #eta ";
5154 MonitorElement * DenominatorPFEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5155 DenominatorPFEtaForward->
getTH1();
5157 histoname = labelname+
"_DenominatorPFPhiForward";
5158 title = labelname+
"DenominatorPFPhiForward;PF #Phi";
5159 MonitorElement * DenominatorPFPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5160 DenominatorPFPhiForward->
getTH1();
5162 histoname = labelname+
"_DenominatorPFEta_LowpTcut";
5163 title = labelname+
"DenominatorPFEta_LowpTcut;PF #eta ";
5164 MonitorElement * DenominatorPFEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5165 DenominatorPFEta_LowpTcut->
getTH1();
5167 histoname = labelname+
"_DenominatorPFPhi_LowpTcut";
5168 title = labelname+
"DenominatorPFPhi_LowpTcut;PF #Phi";
5169 MonitorElement * DenominatorPFPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5170 DenominatorPFPhi_LowpTcut->
getTH1();
5172 histoname = labelname+
"_DenominatorPFEtaPhi_LowpTcut";
5173 title = labelname+
"DenominatorPFEtaPhi_LowpTcut;PF #eta;PF #Phi";
5174 MonitorElement * DenominatorPFEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5175 DenominatorPFEtaPhi_LowpTcut->
getTH1();
5177 histoname = labelname+
"_DenominatorPFEta_MedpTcut";
5178 title = labelname+
"DenominatorPFEta_MedpTcut;PF #eta ";
5179 MonitorElement * DenominatorPFEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5180 DenominatorPFEta_MedpTcut->
getTH1();
5182 histoname = labelname+
"_DenominatorPFPhi_MedpTcut";
5183 title = labelname+
"DenominatorPFPhi_MedpTcut;PF #Phi";
5184 MonitorElement * DenominatorPFPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5185 DenominatorPFPhi_MedpTcut->
getTH1();
5187 histoname = labelname+
"_DenominatorPFEtaPhi_MedpTcut";
5188 title = labelname+
"DenominatorPFEtaPhi_MedpTcut;PF #eta;PF #Phi";
5189 MonitorElement * DenominatorPFEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5190 DenominatorPFEtaPhi_MedpTcut->
getTH1();
5192 histoname = labelname+
"_DenominatorPFEta_HighpTcut";
5193 title = labelname+
"DenominatorPFEta_HighpTcut;PF #eta ";
5194 MonitorElement * DenominatorPFEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5195 DenominatorPFEta_HighpTcut->
getTH1();
5197 histoname = labelname+
"_DenominatorPFPhi_HighpTcut";
5198 title = labelname+
"DenominatorPFPhi_HighpTcut;PF #Phi";
5199 MonitorElement * DenominatorPFPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5200 DenominatorPFPhi_HighpTcut->
getTH1();
5202 histoname = labelname+
"_DenominatorPFEtaPhi_HighpTcut";
5203 title = labelname+
"DenominatorPFEtaPhi_HighpTcut;PF #eta;PF #Phi";
5204 MonitorElement * DenominatorPFEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5205 DenominatorPFEtaPhi_HighpTcut->
getTH1();
5207 histoname = labelname+
"_PFDeltaR";
5208 title = labelname+
"PFDeltaR;";
5212 histoname = labelname+
"_PFDeltaPhi";
5213 title = labelname+
"PFDeltaPhi;";
5217 v->setEffHistos(NumeratorPt,
5229 NumeratorEtaForward,
5230 NumeratorPhiForward,
5231 NumeratorEta_LowpTcut,
5232 NumeratorPhi_LowpTcut,
5233 NumeratorEtaPhi_LowpTcut,
5234 NumeratorEta_MedpTcut,
5235 NumeratorPhi_MedpTcut,
5236 NumeratorEtaPhi_MedpTcut,
5237 NumeratorEta_HighpTcut,
5238 NumeratorPhi_HighpTcut,
5239 NumeratorEtaPhi_HighpTcut,
5242 DenominatorPtBarrel,
5243 DenominatorPtEndcap,
5244 DenominatorPtForward,
5249 DenominatorEtaBarrel,
5250 DenominatorPhiBarrel,
5251 DenominatorEtaEndcap,
5252 DenominatorPhiEndcap,
5253 DenominatorEtaForward,
5254 DenominatorPhiForward,
5255 DenominatorEta_LowpTcut,
5256 DenominatorPhi_LowpTcut,
5257 DenominatorEtaPhi_LowpTcut,
5258 DenominatorEta_MedpTcut,
5259 DenominatorPhi_MedpTcut,
5260 DenominatorEtaPhi_MedpTcut,
5261 DenominatorEta_HighpTcut,
5262 DenominatorPhi_HighpTcut,
5263 DenominatorEtaPhi_HighpTcut,
5270 NumeratorPFPtBarrel,
5271 NumeratorPFPtEndcap,
5272 NumeratorPFPtForward,
5276 NumeratorPFEtaBarrel,
5277 NumeratorPFPhiBarrel,
5278 NumeratorPFEtaEndcap,
5279 NumeratorPFPhiEndcap,
5280 NumeratorPFEtaForward,
5281 NumeratorPFPhiForward,
5282 NumeratorPFEta_LowpTcut,
5283 NumeratorPFPhi_LowpTcut,
5284 NumeratorPFEtaPhi_LowpTcut,
5285 NumeratorPFEta_MedpTcut,
5286 NumeratorPFPhi_MedpTcut,
5287 NumeratorPFEtaPhi_MedpTcut,
5288 NumeratorPFEta_HighpTcut,
5289 NumeratorPFPhi_HighpTcut,
5290 NumeratorPFEtaPhi_HighpTcut,
5293 DenominatorPFPtBarrel,
5294 DenominatorPFPtEndcap,
5295 DenominatorPFPtForward,
5298 DenominatorPFEtaPhi,
5299 DenominatorPFEtaBarrel,
5300 DenominatorPFPhiBarrel,
5301 DenominatorPFEtaEndcap,
5302 DenominatorPFPhiEndcap,
5303 DenominatorPFEtaForward,
5304 DenominatorPFPhiForward,
5305 DenominatorPFEta_LowpTcut,
5306 DenominatorPFPhi_LowpTcut,
5307 DenominatorPFEtaPhi_LowpTcut,
5308 DenominatorPFEta_MedpTcut,
5309 DenominatorPFPhi_MedpTcut,
5310 DenominatorPFEtaPhi_MedpTcut,
5311 DenominatorPFEta_HighpTcut,
5312 DenominatorPFPhi_HighpTcut,
5313 DenominatorPFEtaPhi_HighpTcut,
5319 if((
v->getObjectType() ==
trigger::TriggerJet) && (
v->getTriggerType().compare(
"DiJet_Trigger") == 0)){
5321 histoname = labelname+
"_NumeratorAvrgPt";
5322 title = labelname+
"NumeratorAvrgPt;Calo Pt[GeV/c] ";
5323 MonitorElement * NumeratorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5326 histoname = labelname+
"_NumeratorAvrgEta";
5327 title = labelname+
"NumeratorAvrgEta;Calo #eta";
5328 MonitorElement * NumeratorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5331 histoname = labelname+
"_DenominatorAvrgPt";
5332 title = labelname+
"DenominatorAvrgPt;Calo Pt[GeV/c]";
5333 MonitorElement * DenominatorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5336 histoname = labelname+
"_DenominatorAvrgEta";
5337 title = labelname+
"DenominatorAvrgEta;Calo #eta ";
5338 MonitorElement * DenominatorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5339 DenominatorEta->
getTH1();
5341 histoname = labelname+
"_DeltaR";
5342 title = labelname+
"DeltaR;";
5346 histoname = labelname+
"_DeltaPhi";
5347 title = labelname+
"DeltaPhi;";
5351 histoname = labelname+
"_NumeratorAvrgPFPt";
5352 title = labelname+
"NumeratorAvrgPFPt;PF Pt[GeV/c] ";
5353 MonitorElement * NumeratorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5356 histoname = labelname+
"_NumeratorAvrgPFEta";
5357 title = labelname+
"NumeratorAvrgPFEta;PF #eta";
5358 MonitorElement * NumeratorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5359 NumeratorPFEta->
getTH1();
5361 histoname = labelname+
"_DenominatorAvrgPFPt";
5362 title = labelname+
"DenominatorAvrgPFPt;PF Pt[GeV/c]";
5363 MonitorElement * DenominatorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5364 DenominatorPFPt->
getTH1();
5366 histoname = labelname+
"_DenominatorAvrgPFEta";
5367 title = labelname+
"DenominatorAvrgPFEta;PF #eta ";
5368 MonitorElement * DenominatorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5369 DenominatorPFEta->
getTH1();
5371 histoname = labelname+
"_PFDeltaR";
5372 title = labelname+
"PFDeltaR;";
5376 histoname = labelname+
"_PFDeltaPhi";
5377 title = labelname+
"PFDeltaPhi;";
5381 v->setEffHistos( dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5382 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5383 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5384 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5385 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5386 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5387 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5388 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5389 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5390 dummy, dummy, dummy, dummy
5396 histoname = labelname+
"_NumeratorPt";
5397 if(
v->getPath().find(
"HLT_PFMET")==std::string::npos)
5398 title = labelname+
"NumeratorPt; CaloMET[GeV/c]";
5400 title = labelname+
"NumeratorPt; PFMET[GeV/c]";
5401 MonitorElement * NumeratorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5404 histoname = labelname+
"_NumeratorPhi";
5405 title = labelname+
"NumeratorPhi; #Phi";
5406 MonitorElement * NumeratorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5409 histoname = labelname+
"_DenominatorPt";
5410 if(
v->getPath().find(
"HLT_PFMET")==std::string::npos)
5411 title = labelname+
"DenominatorPt; CaloMET[GeV/c]";
5413 title = labelname+
"DenominatorPt; PFMET[GeV/c]";
5414 MonitorElement * DenominatorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5417 histoname = labelname+
"_DenominatorPhi";
5418 title = labelname+
"DenominatorPhi; #Phi";
5419 MonitorElement * DenominatorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5420 DenominatorPhi->
getTH1();
5422 v->setEffHistos( NumeratorPt, dummy, dummy, dummy, dummy, NumeratorPhi, dummy, dummy, dummy, dummy,
5423 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5424 dummy, dummy, DenominatorPt, dummy, dummy, dummy, dummy, DenominatorPhi, dummy, dummy,
5425 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5426 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5427 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5428 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5429 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5430 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5431 dummy, dummy, dummy, dummy
5453 if((
v->getObjectType() ==
trigger::TriggerJet) && (
v->getTriggerType().compare(
"SingleJet_Trigger") == 0)){
5455 histoname = labelname+
"_NumeratorPt";
5456 title = labelname+
"NumeratorPt;Calo Pt[GeV/c] ";
5457 MonitorElement * NumeratorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5460 histoname = labelname+
"_NumeratorPtBarrel";
5461 title = labelname+
"NumeratorPtBarrel;Calo Pt[GeV/c]";
5462 MonitorElement * NumeratorPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5463 NumeratorPtBarrel->
getTH1();
5465 histoname = labelname+
"_NumeratorPtEndcap";
5466 title = labelname+
"NumeratorPtEndcap; Calo Pt[GeV/c] ";
5467 MonitorElement * NumeratorPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5468 NumeratorPtEndcap->
getTH1();
5470 histoname = labelname+
"_NumeratorPtForward";
5471 title = labelname+
"NumeratorPtForward;Calo Pt[GeV/c]";
5472 MonitorElement * NumeratorPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5473 NumeratorPtForward->
getTH1();
5475 histoname = labelname+
"_NumeratorEta";
5476 title = labelname+
"NumeratorEta;Calo #eta ";
5477 MonitorElement * NumeratorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5480 histoname = labelname+
"_NumeratorPhi";
5481 title = labelname+
"NumeratorPhi;Calo #Phi";
5482 MonitorElement * NumeratorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5485 histoname = labelname+
"_NumeratorEtaPhi";
5486 title = labelname+
"NumeratorEtaPhi;Calo #eta;Calo #Phi ";
5487 MonitorElement * NumeratorEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5488 NumeratorEtaPhi->
getTH1();
5490 histoname = labelname+
"_NumeratorEtaBarrel";
5491 title = labelname+
"NumeratorEtaBarrel;Calo #eta ";
5492 MonitorElement * NumeratorEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5493 NumeratorEtaBarrel->
getTH1();
5495 histoname = labelname+
"_NumeratorPhiBarrel";
5496 title = labelname+
"NumeratorPhiBarrel;Calo #Phi";
5497 MonitorElement * NumeratorPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5498 NumeratorPhiBarrel->
getTH1();
5500 histoname = labelname+
"_NumeratorEtaEndcap";
5501 title = labelname+
"NumeratorEtaEndcap;Calo #eta ";
5502 MonitorElement * NumeratorEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5503 NumeratorEtaEndcap->
getTH1();
5505 histoname = labelname+
"_NumeratorPhiEndcap";
5506 title = labelname+
"NumeratorPhiEndcap;Calo #Phi";
5507 MonitorElement * NumeratorPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5508 NumeratorPhiEndcap->
getTH1();
5510 histoname = labelname+
"_NumeratorEtaForward";
5511 title = labelname+
"NumeratorEtaForward;Calo #eta ";
5512 MonitorElement * NumeratorEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5513 NumeratorEtaForward->
getTH1();
5515 histoname = labelname+
"_NumeratorPhiForward";
5516 title = labelname+
"NumeratorPhiForward;Calo #Phi";
5517 MonitorElement * NumeratorPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5518 NumeratorPhiForward->
getTH1();
5520 histoname = labelname+
"_NumeratorEta_LowpTcut";
5521 title = labelname+
"NumeratorEta_LowpTcut;Calo #eta ";
5522 MonitorElement * NumeratorEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5523 NumeratorEta_LowpTcut->
getTH1();
5525 histoname = labelname+
"_NumeratorPhi_LowpTcut";
5526 title = labelname+
"NumeratorPhi_LowpTcut;Calo #Phi";
5527 MonitorElement * NumeratorPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5528 NumeratorPhi_LowpTcut->
getTH1();
5530 histoname = labelname+
"_NumeratorEtaPhi_LowpTcut";
5531 title = labelname+
"NumeratorEtaPhi_LowpTcut;Calo #eta;Calo #Phi";
5532 MonitorElement * NumeratorEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5533 NumeratorEtaPhi_LowpTcut->
getTH1();
5535 histoname = labelname+
"_NumeratorEta_MedpTcut";
5536 title = labelname+
"NumeratorEta_MedpTcut;Calo #eta ";
5537 MonitorElement * NumeratorEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5538 NumeratorEta_MedpTcut->
getTH1();
5540 histoname = labelname+
"_NumeratorPhi_MedpTcut";
5541 title = labelname+
"NumeratorPhi_MedpTcut;Calo #Phi";
5542 MonitorElement * NumeratorPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5543 NumeratorPhi_MedpTcut->
getTH1();
5545 histoname = labelname+
"_NumeratorEtaPhi_MedpTcut";
5546 title = labelname+
"NumeratorEtaPhi_MedpTcut;Calo #eta;Calo #Phi";
5547 MonitorElement * NumeratorEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5548 NumeratorEtaPhi_MedpTcut->
getTH1();
5550 histoname = labelname+
"_NumeratorEta_HighpTcut";
5551 title = labelname+
"NumeratorEta_HighpTcut;Calo #eta ";
5552 MonitorElement * NumeratorEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5553 NumeratorEta_HighpTcut->
getTH1();
5555 histoname = labelname+
"_NumeratorPhi_HighpTcut";
5556 title = labelname+
"NumeratorPhi_HighpTcut;Calo #Phi";
5557 MonitorElement * NumeratorPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5558 NumeratorPhi_HighpTcut->
getTH1();
5560 histoname = labelname+
"_NumeratorEtaPhi_HighpTcut";
5561 title = labelname+
"NumeratorEtaPhi_HighpTcut;Calo #eta;Calo #Phi";
5562 MonitorElement * NumeratorEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5563 NumeratorEtaPhi_HighpTcut->
getTH1();
5565 histoname = labelname+
"_DenominatorPt";
5566 title = labelname+
"DenominatorPt;Calo Pt[GeV/c]";
5567 MonitorElement * DenominatorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5570 histoname = labelname+
"_DenominatorPtBarrel";
5571 title = labelname+
"DenominatorPtBarrel;Calo Pt[GeV/c]";
5572 MonitorElement * DenominatorPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5573 DenominatorPtBarrel->
getTH1();
5575 histoname = labelname+
"_DenominatorPtEndcap";
5576 title = labelname+
"DenominatorPtEndcap;Calo Pt[GeV/c]";
5577 MonitorElement * DenominatorPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5578 DenominatorPtEndcap->
getTH1();
5580 histoname = labelname+
"_DenominatorPtForward";
5581 title = labelname+
"DenominatorPtForward;Calo Pt[GeV/c]";
5582 MonitorElement * DenominatorPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5583 DenominatorPtForward->
getTH1();
5585 histoname = labelname+
"_DenominatorEta";
5586 title = labelname+
"DenominatorEta;Calo #eta ";
5587 MonitorElement * DenominatorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5588 DenominatorEta->
getTH1();
5590 histoname = labelname+
"_DenominatorPhi";
5591 title = labelname+
"DenominatorPhi;Calo #Phi";
5592 MonitorElement * DenominatorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5593 DenominatorPhi->
getTH1();
5595 histoname = labelname+
"_DenominatorEtaPhi";
5596 title = labelname+
"DenominatorEtaPhi;Calo #eta ;Calo #Phi ";
5597 MonitorElement * DenominatorEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5598 DenominatorEtaPhi->
getTH1();
5600 histoname = labelname+
"_DenominatorEtaBarrel";
5601 title = labelname+
"DenominatorEtaBarrel;Calo #eta ";
5602 MonitorElement * DenominatorEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5603 DenominatorEtaBarrel->
getTH1();
5605 histoname = labelname+
"_DenominatorPhiBarrel";
5606 title = labelname+
"DenominatorPhiBarrel;Calo #Phi";
5607 MonitorElement * DenominatorPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5608 DenominatorPhiBarrel->
getTH1();
5610 histoname = labelname+
"_DenominatorEtaEndcap";
5611 title = labelname+
"DenominatorEtaEndcap;Calo #eta ";
5612 MonitorElement * DenominatorEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5613 DenominatorEtaEndcap->
getTH1();
5615 histoname = labelname+
"_DenominatorPhiEndcap";
5616 title = labelname+
"DenominatorPhiEndcap;Calo #Phi";
5617 MonitorElement * DenominatorPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5618 DenominatorPhiEndcap->
getTH1();
5620 histoname = labelname+
"_DenominatorEtaForward";
5621 title = labelname+
"DenominatorEtaForward;Calo #eta ";
5622 MonitorElement * DenominatorEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5623 DenominatorEtaForward->
getTH1();
5625 histoname = labelname+
"_DenominatorPhiForward";
5626 title = labelname+
"DenominatorPhiForward;Calo #Phi";
5627 MonitorElement * DenominatorPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5628 DenominatorPhiForward->
getTH1();
5630 histoname = labelname+
"_DenominatorEta_LowpTcut";
5631 title = labelname+
"DenominatorEta_LowpTcut;Calo #eta ";
5632 MonitorElement * DenominatorEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5633 DenominatorEta_LowpTcut->
getTH1();
5635 histoname = labelname+
"_DenominatorPhi_LowpTcut";
5636 title = labelname+
"DenominatorPhi_LowpTcut;Calo #Phi";
5637 MonitorElement * DenominatorPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5638 DenominatorPhi_LowpTcut->
getTH1();
5640 histoname = labelname+
"_DenominatorEtaPhi_LowpTcut";
5641 title = labelname+
"DenominatorEtaPhi_LowpTcut;Calo #eta;Calo #Phi";
5642 MonitorElement * DenominatorEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5643 DenominatorEtaPhi_LowpTcut->
getTH1();
5645 histoname = labelname+
"_DenominatorEta_MedpTcut";
5646 title = labelname+
"DenominatorEta_MedpTcut;Calo #eta ";
5647 MonitorElement * DenominatorEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5648 DenominatorEta_MedpTcut->
getTH1();
5650 histoname = labelname+
"_DenominatorPhi_MedpTcut";
5651 title = labelname+
"DenominatorPhi_MedpTcut;Calo #Phi";
5652 MonitorElement * DenominatorPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5653 DenominatorPhi_MedpTcut->
getTH1();
5655 histoname = labelname+
"_DenominatorEtaPhi_MedpTcut";
5656 title = labelname+
"DenominatorEtaPhi_MedpTcut;Calo #eta;Calo #Phi";
5657 MonitorElement * DenominatorEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5658 DenominatorEtaPhi_MedpTcut->
getTH1();
5660 histoname = labelname+
"_DenominatorEta_HighpTcut";
5661 title = labelname+
"DenominatorEta_HighpTcut;Calo #eta ";
5662 MonitorElement * DenominatorEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5663 DenominatorEta_HighpTcut->
getTH1();
5665 histoname = labelname+
"_DenominatorPhi_HighpTcut";
5666 title = labelname+
"DenominatorPhi_HighpTcut;Calo #Phi";
5667 MonitorElement * DenominatorPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5668 DenominatorPhi_HighpTcut->
getTH1();
5670 histoname = labelname+
"_DenominatorEtaPhi_HighpTcut";
5671 title = labelname+
"DenominatorEtaPhi_HighpTcut;Calo #eta;Calo #Phi";
5672 MonitorElement * DenominatorEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5673 DenominatorEtaPhi_HighpTcut->
getTH1();
5675 histoname = labelname+
"_DeltaR";
5676 title = labelname+
"DeltaR;";
5680 histoname = labelname+
"_DeltaPhi";
5681 title = labelname+
"DeltaPhi;";
5685 histoname = labelname+
"_NumeratorPFPt";
5686 title = labelname+
"NumeratorPFPt;PF Pt[GeV/c] ";
5687 MonitorElement * NumeratorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5690 histoname = labelname+
"_NumeratorPFMHT";
5691 title = labelname+
"NumeratorPFPt;PFMHT[GeV/c] ";
5692 MonitorElement * NumeratorPFMHT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5693 NumeratorPFMHT->
getTH1();
5695 histoname = labelname+
"_NumeratorPFPtBarrel";
5696 title = labelname+
"NumeratorPFPtBarrel;PF Pt[GeV/c]";
5697 MonitorElement * NumeratorPFPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5698 NumeratorPFPtBarrel->
getTH1();
5700 histoname = labelname+
"_NumeratorPFPtEndcap";
5701 title = labelname+
"NumeratorPFPtEndcap; PF Pt[GeV/c] ";
5702 MonitorElement * NumeratorPFPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5703 NumeratorPFPtEndcap->
getTH1();
5705 histoname = labelname+
"_NumeratorPFPtForward";
5706 title = labelname+
"NumeratorPFPtForward;PF Pt[GeV/c]";
5707 MonitorElement * NumeratorPFPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5708 NumeratorPFPtForward->
getTH1();
5710 histoname = labelname+
"_NumeratorPFEta";
5711 title = labelname+
"NumeratorPFEta;PF #eta ";
5712 MonitorElement * NumeratorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5713 NumeratorPFEta->
getTH1();
5715 histoname = labelname+
"_NumeratorPFPhi";
5716 title = labelname+
"NumeratorPFPhi;PF #Phi";
5717 MonitorElement * NumeratorPFPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5718 NumeratorPFPhi->
getTH1();
5720 histoname = labelname+
"_NumeratorPFEtaPhi";
5721 title = labelname+
"NumeratorPFEtaPhi;PF #eta;PF #Phi ";
5722 MonitorElement * NumeratorPFEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5723 NumeratorPFEtaPhi->
getTH1();
5725 histoname = labelname+
"_NumeratorPFEtaBarrel";
5726 title = labelname+
"NumeratorPFEtaBarrel;PF #eta ";
5727 MonitorElement * NumeratorPFEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5728 NumeratorPFEtaBarrel->
getTH1();
5730 histoname = labelname+
"_NumeratorPFPhiBarrel";
5731 title = labelname+
"NumeratorPFPhiBarrel;PF #Phi";
5732 MonitorElement * NumeratorPFPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5733 NumeratorPFPhiBarrel->
getTH1();
5735 histoname = labelname+
"_NumeratorPFEtaEndcap";
5736 title = labelname+
"NumeratorPFEtaEndcap;PF #eta ";
5737 MonitorElement * NumeratorPFEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5738 NumeratorPFEtaEndcap->
getTH1();
5740 histoname = labelname+
"_NumeratorPFPhiEndcap";
5741 title = labelname+
"NumeratorPFPhiEndcap;PF #Phi";
5742 MonitorElement * NumeratorPFPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5743 NumeratorPFPhiEndcap->
getTH1();
5745 histoname = labelname+
"_NumeratorPFEtaForward";
5746 title = labelname+
"NumeratorPFEtaForward;PF #eta ";
5747 MonitorElement * NumeratorPFEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5748 NumeratorPFEtaForward->
getTH1();
5750 histoname = labelname+
"_NumeratorPFPhiForward";
5751 title = labelname+
"NumeratorPFPhiForward;PF #Phi";
5752 MonitorElement * NumeratorPFPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5753 NumeratorPFPhiForward->
getTH1();
5755 histoname = labelname+
"_NumeratorPFEta_LowpTcut";
5756 title = labelname+
"NumeratorPFEta_LowpTcut;PF #eta ";
5757 MonitorElement * NumeratorPFEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5758 NumeratorPFEta_LowpTcut->
getTH1();
5760 histoname = labelname+
"_NumeratorPFPhi_LowpTcut";
5761 title = labelname+
"NumeratorPFPhi_LowpTcut;PF #Phi";
5762 MonitorElement * NumeratorPFPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5763 NumeratorPFPhi_LowpTcut->
getTH1();
5765 histoname = labelname+
"_NumeratorPFEtaPhi_LowpTcut";
5766 title = labelname+
"NumeratorPFEtaPhi_LowpTcut;PF #eta;PF #Phi";
5767 MonitorElement * NumeratorPFEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5768 NumeratorPFEtaPhi_LowpTcut->
getTH1();
5770 histoname = labelname+
"_NumeratorPFEta_MedpTcut";
5771 title = labelname+
"NumeratorPFEta_MedpTcut;PF #eta ";
5772 MonitorElement * NumeratorPFEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5773 NumeratorPFEta_MedpTcut->
getTH1();
5775 histoname = labelname+
"_NumeratorPFPhi_MedpTcut";
5776 title = labelname+
"NumeratorPFPhi_MedpTcut;PF #Phi";
5777 MonitorElement * NumeratorPFPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5778 NumeratorPFPhi_MedpTcut->
getTH1();
5780 histoname = labelname+
"_NumeratorPFEtaPhi_MedpTcut";
5781 title = labelname+
"NumeratorPFEtaPhi_MedpTcut;PF #eta;PF #Phi";
5782 MonitorElement * NumeratorPFEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5783 NumeratorPFEtaPhi_MedpTcut->
getTH1();
5785 histoname = labelname+
"_NumeratorPFEta_HighpTcut";
5786 title = labelname+
"NumeratorPFEta_HighpTcut;PF #eta ";
5787 MonitorElement * NumeratorPFEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5788 NumeratorPFEta_HighpTcut->
getTH1();
5790 histoname = labelname+
"_NumeratorPFPhi_HighpTcut";
5791 title = labelname+
"NumeratorPFPhi_HighpTcut;PF #Phi";
5792 MonitorElement * NumeratorPFPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5793 NumeratorPFPhi_HighpTcut->
getTH1();
5795 histoname = labelname+
"_NumeratorPFEtaPhi_HighpTcut";
5796 title = labelname+
"NumeratorPFEtaPhi_HighpTcut;PF #eta;PF #Phi";
5797 MonitorElement * NumeratorPFEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5798 NumeratorPFEtaPhi_HighpTcut->
getTH1();
5800 histoname = labelname+
"_DenominatorPFMHT";
5801 title = labelname+
"DenominatorPFMHT;PFMHT[GeV/c]";
5802 MonitorElement * DenominatorPFMHT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5803 DenominatorPFMHT->
getTH1();
5805 histoname = labelname+
"_DenominatorPFPt";
5806 title = labelname+
"DenominatorPFPt;PF Pt[GeV/c]";
5807 MonitorElement * DenominatorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5808 DenominatorPFPt->
getTH1();
5810 histoname = labelname+
"_DenominatorPFPtBarrel";
5811 title = labelname+
"DenominatorPFPtBarrel;PF Pt[GeV/c]";
5812 MonitorElement * DenominatorPFPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5813 DenominatorPFPtBarrel->
getTH1();
5815 histoname = labelname+
"_DenominatorPFPtEndcap";
5816 title = labelname+
"DenominatorPFPtEndcap;PF Pt[GeV/c]";
5817 MonitorElement * DenominatorPFPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5818 DenominatorPFPtEndcap->
getTH1();
5820 histoname = labelname+
"_DenominatorPFPtForward";
5821 title = labelname+
"DenominatorPFPtForward;PF Pt[GeV/c]";
5822 MonitorElement * DenominatorPFPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5823 DenominatorPFPtForward->
getTH1();
5825 histoname = labelname+
"_DenominatorPFEta";
5826 title = labelname+
"DenominatorPFEta;PF #eta ";
5827 MonitorElement * DenominatorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5828 DenominatorPFEta->
getTH1();
5830 histoname = labelname+
"_DenominatorPFPhi";
5831 title = labelname+
"DenominatorPFPhi;PF #Phi";
5832 MonitorElement * DenominatorPFPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5833 DenominatorPFPhi->
getTH1();
5835 histoname = labelname+
"_DenominatorPFEtaPhi";
5836 title = labelname+
"DenominatorPFEtaPhi;PF #eta ;PF #Phi ";
5837 MonitorElement * DenominatorPFEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5838 DenominatorPFEtaPhi->
getTH1();
5840 histoname = labelname+
"_DenominatorPFEtaBarrel";
5841 title = labelname+
"DenominatorPFEtaBarrel;PF #eta ";
5842 MonitorElement * DenominatorPFEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5843 DenominatorPFEtaBarrel->
getTH1();
5845 histoname = labelname+
"_DenominatorPFPhiBarrel";
5846 title = labelname+
"DenominatorPFPhiBarrel;PF #Phi";
5847 MonitorElement * DenominatorPFPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5848 DenominatorPFPhiBarrel->
getTH1();
5850 histoname = labelname+
"_DenominatorPFEtaEndcap";
5851 title = labelname+
"DenominatorPFEtaEndcap;PF #eta ";
5852 MonitorElement * DenominatorPFEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5853 DenominatorPFEtaEndcap->
getTH1();
5855 histoname = labelname+
"_DenominatorPFPhiEndcap";
5856 title = labelname+
"DenominatorPFPhiEndcap;PF #Phi";
5857 MonitorElement * DenominatorPFPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5858 DenominatorPFPhiEndcap->
getTH1();
5860 histoname = labelname+
"_DenominatorPFEtaForward";
5861 title = labelname+
"DenominatorPFEtaForward;PF #eta ";
5862 MonitorElement * DenominatorPFEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5863 DenominatorPFEtaForward->
getTH1();
5865 histoname = labelname+
"_DenominatorPFPhiForward";
5866 title = labelname+
"DenominatorPFPhiForward;PF #Phi";
5867 MonitorElement * DenominatorPFPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5868 DenominatorPFPhiForward->
getTH1();
5870 histoname = labelname+
"_DenominatorPFEta_LowpTcut";
5871 title = labelname+
"DenominatorPFEta_LowpTcut;PF #eta ";
5872 MonitorElement * DenominatorPFEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5873 DenominatorPFEta_LowpTcut->
getTH1();
5875 histoname = labelname+
"_DenominatorPFPhi_LowpTcut";
5876 title = labelname+
"DenominatorPFPhi_LowpTcut;PF #Phi";
5877 MonitorElement * DenominatorPFPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5878 DenominatorPFPhi_LowpTcut->
getTH1();
5880 histoname = labelname+
"_DenominatorPFEtaPhi_LowpTcut";
5881 title = labelname+
"DenominatorPFEtaPhi_LowpTcut;PF #eta;PF #Phi";
5882 MonitorElement * DenominatorPFEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5883 DenominatorPFEtaPhi_LowpTcut->
getTH1();
5885 histoname = labelname+
"_DenominatorPFEta_MedpTcut";
5886 title = labelname+
"DenominatorPFEta_MedpTcut;PF #eta ";
5887 MonitorElement * DenominatorPFEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5888 DenominatorPFEta_MedpTcut->
getTH1();
5890 histoname = labelname+
"_DenominatorPFPhi_MedpTcut";
5891 title = labelname+
"DenominatorPFPhi_MedpTcut;PF #Phi";
5892 MonitorElement * DenominatorPFPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5893 DenominatorPFPhi_MedpTcut->
getTH1();
5895 histoname = labelname+
"_DenominatorPFEtaPhi_MedpTcut";
5896 title = labelname+
"DenominatorPFEtaPhi_MedpTcut;PF #eta;PF #Phi";
5897 MonitorElement * DenominatorPFEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5898 DenominatorPFEtaPhi_MedpTcut->
getTH1();
5900 histoname = labelname+
"_DenominatorPFEta_HighpTcut";
5901 title = labelname+
"DenominatorPFEta_HighpTcut;PF #eta ";
5902 MonitorElement * DenominatorPFEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5903 DenominatorPFEta_HighpTcut->
getTH1();
5905 histoname = labelname+
"_DenominatorPFPhi_HighpTcut";
5906 title = labelname+
"DenominatorPFPhi_HighpTcut;PF #Phi";
5907 MonitorElement * DenominatorPFPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5908 DenominatorPFPhi_HighpTcut->
getTH1();
5910 histoname = labelname+
"_DenominatorPFEtaPhi_HighpTcut";
5911 title = labelname+
"DenominatorPFEtaPhi_HighpTcut;PF #eta;PF #Phi";
5912 MonitorElement * DenominatorPFEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5913 DenominatorPFEtaPhi_HighpTcut->
getTH1();
5915 histoname = labelname+
"_PFDeltaR";
5916 title = labelname+
"PFDeltaR;";
5920 histoname = labelname+
"_PFDeltaPhi";
5921 title = labelname+
"PFDeltaPhi;";
5925 v->setEffHistos(NumeratorPt,
5937 NumeratorEtaForward,
5938 NumeratorPhiForward,
5939 NumeratorEta_LowpTcut,
5940 NumeratorPhi_LowpTcut,
5941 NumeratorEtaPhi_LowpTcut,
5942 NumeratorEta_MedpTcut,
5943 NumeratorPhi_MedpTcut,
5944 NumeratorEtaPhi_MedpTcut,
5945 NumeratorEta_HighpTcut,
5946 NumeratorPhi_HighpTcut,
5947 NumeratorEtaPhi_HighpTcut,
5950 DenominatorPtBarrel,
5951 DenominatorPtEndcap,
5952 DenominatorPtForward,
5957 DenominatorEtaBarrel,
5958 DenominatorPhiBarrel,
5959 DenominatorEtaEndcap,
5960 DenominatorPhiEndcap,
5961 DenominatorEtaForward,
5962 DenominatorPhiForward,
5963 DenominatorEta_LowpTcut,
5964 DenominatorPhi_LowpTcut,
5965 DenominatorEtaPhi_LowpTcut,
5966 DenominatorEta_MedpTcut,
5967 DenominatorPhi_MedpTcut,
5968 DenominatorEtaPhi_MedpTcut,
5969 DenominatorEta_HighpTcut,
5970 DenominatorPhi_HighpTcut,
5971 DenominatorEtaPhi_HighpTcut,
5977 NumeratorPFPtBarrel,
5978 NumeratorPFPtEndcap,
5979 NumeratorPFPtForward,
5983 NumeratorPFEtaBarrel,
5984 NumeratorPFPhiBarrel,
5985 NumeratorPFEtaEndcap,
5986 NumeratorPFPhiEndcap,
5987 NumeratorPFEtaForward,
5988 NumeratorPFPhiForward,
5989 NumeratorPFEta_LowpTcut,
5990 NumeratorPFPhi_LowpTcut,
5991 NumeratorPFEtaPhi_LowpTcut,
5992 NumeratorPFEta_MedpTcut,
5993 NumeratorPFPhi_MedpTcut,
5994 NumeratorPFEtaPhi_MedpTcut,
5995 NumeratorPFEta_HighpTcut,
5996 NumeratorPFPhi_HighpTcut,
5997 NumeratorPFEtaPhi_HighpTcut,
6000 DenominatorPFPtBarrel,
6001 DenominatorPFPtEndcap,
6002 DenominatorPFPtForward,
6005 DenominatorPFEtaPhi,
6006 DenominatorPFEtaBarrel,
6007 DenominatorPFPhiBarrel,
6008 DenominatorPFEtaEndcap,
6009 DenominatorPFPhiEndcap,
6010 DenominatorPFEtaForward,
6011 DenominatorPFPhiForward,
6012 DenominatorPFEta_LowpTcut,
6013 DenominatorPFPhi_LowpTcut,
6014 DenominatorPFEtaPhi_LowpTcut,
6015 DenominatorPFEta_MedpTcut,
6016 DenominatorPFPhi_MedpTcut,
6017 DenominatorPFEtaPhi_MedpTcut,
6018 DenominatorPFEta_HighpTcut,
6019 DenominatorPFPhi_HighpTcut,
6020 DenominatorPFEtaPhi_HighpTcut,
6027 if((
v->getObjectType() ==
trigger::TriggerJet) && (
v->getTriggerType().compare(
"DiJet_Trigger") == 0)){
6029 histoname = labelname+
"_NumeratorAvrgPt";
6030 title = labelname+
"NumeratorAvrgPt;Calo Pt[GeV/c] ";
6031 MonitorElement * NumeratorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
6034 histoname = labelname+
"_NumeratorAvrgEta";
6035 title = labelname+
"NumeratorAvrgEta;Calo #eta ";
6036 MonitorElement * NumeratorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
6039 histoname = labelname+
"_DenominatorAvrgPt";
6040 title = labelname+
"DenominatorAvrgPt;Calo Pt[GeV/c] ";
6041 MonitorElement * DenominatorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
6044 histoname = labelname+
"_DenominatorAvrgEta";
6045 title = labelname+
"DenominatorAvrgEta;Calo #eta ";
6046 MonitorElement * DenominatorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
6047 DenominatorEta->
getTH1();
6049 histoname = labelname+
"_DeltaR";
6050 title = labelname+
"DeltaR;";
6054 histoname = labelname+
"_DeltaPhi";
6055 title = labelname+
"DeltaPhi;";
6059 histoname = labelname+
"_NumeratorAvrgPFPt";
6060 title = labelname+
"NumeratorAvrgPFPt;PF pT [GeV/c] ";
6061 MonitorElement * NumeratorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
6064 histoname = labelname+
"_NumeratorAvrgPFEta";
6065 title = labelname+
"NumeratorAvrgPFEta;PF #eta ";
6066 MonitorElement * NumeratorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
6067 NumeratorPFEta->
getTH1();
6069 histoname = labelname+
"_DenominatorAvrgPFPt";
6070 title = labelname+
"DenominatorAvrgPFPt;PF Pt[GeV/c] ";
6071 MonitorElement * DenominatorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
6072 DenominatorPFPt->
getTH1();
6074 histoname = labelname+
"_DenominatorAvrgPFEta";
6075 title = labelname+
"DenominatorAvrgPFEta;PF #eta ";
6076 MonitorElement * DenominatorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
6077 DenominatorPFEta->
getTH1();
6079 histoname = labelname+
"_PFDeltaR";
6080 title = labelname+
"PFDeltaR;";
6084 histoname = labelname+
"_PFDeltaPhi";
6085 title = labelname+
"PFDeltaPhi;";
6089 v->setEffHistos( dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6090 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6091 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6092 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6093 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6094 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6095 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6096 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6097 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6098 dummy, dummy, dummy, dummy
6104 histoname = labelname+
"_NumeratorPt";
6105 if(
v->getPath().find(
"HLT_PFMET")==std::string::npos)
6106 title = labelname+
"NumeratorPt;CaloMET[GeV/c]";
6108 title = labelname+
"NumeratorPt;PFMET[GeV/c]";
6109 MonitorElement * NumeratorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
6112 histoname = labelname+
"_NumeratorPhi";
6113 title = labelname+
"NumeratorPhi;#Phi";
6114 MonitorElement * NumeratorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
6117 histoname = labelname+
"_DenominatorPt";
6118 if(
v->getPath().find(
"HLT_PFMET")==std::string::npos)
6119 title = labelname+
"DenominatorPt;CaloMET[GeV/c]";
6121 title = labelname+
"DenominatorPt;PFMET[GeV/c]";
6122 MonitorElement * DenominatorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
6125 histoname = labelname+
"_DenominatorPhi";
6126 title = labelname+
"DenominatorPhi;#Phi";
6127 MonitorElement * DenominatorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
6128 DenominatorPhi->
getTH1();
6130 v->setEffHistos( NumeratorPt, dummy, dummy, dummy, dummy, NumeratorPhi, dummy, dummy, dummy, dummy,
6131 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6132 dummy, dummy, DenominatorPt, dummy, dummy, dummy, dummy, DenominatorPhi, dummy, dummy,
6133 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6134 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6135 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6136 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6137 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6138 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6139 dummy, dummy, dummy, dummy
6157 double PtMax_ = 1000.;
6158 double EtaMin_ = -5.;
6159 double EtaMax_ = 5.;
6160 double PhiMin_ = -3.14159;
6161 double PhiMax_ = 3.14159;
6174 histoname = labelname+
"_TriggerSummary";
6175 title = labelname+
"Summary of trigger levels";
6178 std::vector<std::string> trigger;
6179 trigger.push_back(
"Nevt");
6180 trigger.push_back(
"L1 failed");
6181 trigger.push_back(
"L1 & HLT failed");
6182 trigger.push_back(
"L1 failed but not HLT");
6183 trigger.push_back(
"L1 passed");
6184 trigger.push_back(
"L1 & HLT passed");
6185 trigger.push_back(
"L1 passed but not HLT");
6187 for(
unsigned int i =0;
i < trigger.size();
i++)
6190 if((
v->getTriggerType().compare(
"SingleJet_Trigger") == 0)){
6191 histoname = labelname+
"_JetPt";
6192 title = labelname+
"Leading jet pT;Pt[GeV/c]";
6196 histoname = labelname+
"_JetEtaVsPt";
6197 title = labelname+
"Leading jet #eta vs pT;#eta;Pt[GeV/c]";
6198 MonitorElement * JetEtaVsPt = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Ptbins_,PtMin_,PtMax_);
6201 histoname = labelname+
"_JetPhiVsPt";
6202 title = labelname+
"Leading jet #Phi vs pT;#Phi;Pt[GeV/c]";
6203 MonitorElement * JetPhiVsPt = dbe->
book2D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_,Ptbins_,PtMin_,PtMax_);
6206 v->setDgnsHistos( TriggerSummary, dummy, JetPt, JetEtaVsPt, JetPhiVsPt, dummy, dummy, dummy, dummy, dummy, dummy);
6209 if((
v->getTriggerType().compare(
"DiJet_Trigger") == 0)){
6210 histoname = labelname+
"_JetSize";
6211 title = labelname+
"Jet Size;multiplicity";
6215 histoname = labelname+
"_AvergPt";
6216 title = labelname+
"Average Pt;Pt[GeV/c]";
6220 histoname = labelname+
"_AvergEta";
6221 title = labelname+
"Average Eta;#eta";
6222 MonitorElement * Eta12 = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
6225 histoname = labelname+
"_PhiDifference";
6226 title = labelname+
"#Delta#Phi;#Delta#Phi";
6227 MonitorElement * Phi12 = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
6230 histoname = labelname+
"_Pt3Jet";
6231 title = labelname+
"Pt of 3rd Jet;Pt[GeV/c]";
6235 histoname = labelname+
"_Pt12VsPt3Jet";
6236 title = labelname+
"Pt of 3rd Jet vs Average Pt of leading jets;Avergage Pt[GeV/c]; Pt of 3rd Jet [GeV/c]";
6237 MonitorElement * Pt12Pt3 = dbe->
book2D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_,Ptbins_,PtMin_,PtMax_);
6240 histoname = labelname+
"_Pt12VsPhi12";
6241 title = labelname+
"Average Pt of leading jets vs #Delta#Phi between leading jets;Avergage Pt[GeV/c]; #Delta#Phi";
6242 MonitorElement * Pt12Phi12 = dbe->
book2D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_,Phibins_,PhiMin_,PhiMax_);
6245 v->setDgnsHistos( TriggerSummary, JetSize, dummy, dummy, dummy, Pt12, Eta12, Phi12, Pt3, Pt12Pt3, Pt12Phi12);
6248 if((
v->getTriggerType().compare(
"MET_Trigger") == 0)){
6249 histoname = labelname+
"_MET";
6250 title = labelname+
"MET;Pt[GeV/c]";
6254 v->setDgnsHistos(TriggerSummary, dummy, MET, dummy, dummy, dummy, dummy, dummy,dummy,dummy,dummy);
6257 if((
v->getTriggerType().compare(
"TET_Trigger") == 0)){
6258 histoname = labelname+
"_TET";
6259 title = labelname+
"TET;Pt[GeV/c]";
6263 v->setDgnsHistos(TriggerSummary, dummy, TET, dummy, dummy, dummy, dummy, dummy,dummy,dummy,dummy);
6291 for(std::vector<std::string>::iterator numpathmodule = numpathmodules.begin();
6292 numpathmodule!= numpathmodules.end(); ++numpathmodule ) {
6295 l1pathname = *numpathmodule;
6316 if (fabs(eta)<=1.3) output=
true;
6323 if (fabs(eta)<=3.0 && fabs(eta)>1.3) output=
true;
6330 if (fabs(eta)>3.0) output=
true;
6359 double binVal = -100;
6360 for(
int ibin=1; ibin<nbins+1; ibin++)
6362 const char * binLabel =
rate_All->
getTH1()->GetXaxis()->GetBinLabel(ibin);
6363 if(binLabel[0]==
'\0')
continue;
6366 if(trigName.compare(binLabel)!=0)
continue;
6368 if(trigName.compare(binLabel)==0){
6383 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index "<< index <<
" of that name ";
6386 if (k.size()) output=
true;
unsigned int size() const
number of trigger paths in trigger table
PathInfoCollection hltPathsAllWrtMu_
T getParameter(std::string const &) const
EventNumber_t event() const
T getUntrackedParameter(std::string const &, T const &) const
virtual double et() const GCC11_FINAL
transverse energy
virtual edm::TriggerNames const & triggerNames(edm::TriggerResults const &triggerResults) const
edm::EDGetTokenT< edm::TriggerResults > triggerResultsFUToken
edm::EDGetTokenT< reco::VertexCollection > vertexToken
edm::EDGetTokenT< reco::PFJetCollection > pfJetsToken
tuple met
____________________________________________________________________________||
const std::string moduleType(const std::string &module) const
C++ class name of module.
virtual bool isEndCap(double eta)
const std::string & triggerName(unsigned int triggerIndex) const
JetMETHLTOfflineSource(const edm::ParameterSet &)
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
void rmdir(const std::string &fullpath)
edm::TriggerNames triggerNames_
edm::Handle< reco::PFJetCollection > pfjetColl_
void beginLuminosityBlock(const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)
bool getByToken(EDGetToken token, Handle< PROD > &result) const
virtual double TriggerPosition(std::string trigName)
virtual double correction(const LorentzVector &fJet) const =0
get correction using Jet information only
void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
edm::Handle< edm::TriggerResults > triggerResults_
edm::LuminosityBlockNumber_t luminosityBlock() const
edm::InputTag triggerResultsLabel_
reco::PFJetCollection pfjet
Strings::size_type size() const
edm::Handle< reco::PFMETCollection > pfmetColl_
PathInfoCollection hltPathsAllTriggerSummary_
virtual bool validPathHLT(std::string path)
edm::InputTag triggerSummaryLabel_
static const std::string removeVersion(const std::string &trigger)
MonitorElement * bookFloat(const char *name)
Book float.
virtual void fillMEforMonAllTrigger(const edm::Event &iEvent, const edm::EventSetup &)
virtual bool isBarrel(double eta)
MonitorElement * NVertices
virtual float phi() const GCC11_FINAL
momentum azimuthal angle
std::string PFJetCorService_
std::vector< std::string > pathFilter_
virtual void fillMEforTriggerNTfired()
unsigned int triggerIndex(std::string const &name) const
std::vector< std::string > MBTrigPaths_
double dPhi(double phi1, double phi2)
virtual bool isTriggerObjectFound(std::string objectName)
ConsumesCollector consumesCollector()
Use a ConsumesCollector to gather consumes information from helper functions.
virtual bool isForward(double eta)
MonitorElement * correlation_AllWrtMu
std::vector< std::string > pathRejectKeyword_
virtual void fillMEforMonAllTriggerwrtMuonTrigger(const edm::Event &iEvent, const edm::EventSetup &)
std::string CaloJetCorService_
Abs< T >::type abs(const T &t)
~JetMETHLTOfflineSource()
virtual void fillMEforEffWrtMBTrigger(const edm::Event &iEvent, const edm::EventSetup &)
PathInfoCollection hltPathsEffWrtMB_
void setVerbose(unsigned level)
edm::EDGetTokenT< reco::CaloMETCollection > caloMetToken
reco::CaloJetCollection calojet
std::vector< std::string > MuonTrigPaths_
virtual bool isHLTPathAccepted(std::string pathName)
edm::Handle< reco::CaloJetCollection > calojetColl_
const std::vector< std::string > & moduleLabels(unsigned int trigger) const
label(s) of module(s) on a trigger path
double deltaR(double eta1, double eta2, double phi1, double phi2)
edm::EDGetTokenT< trigger::TriggerEvent > triggerSummaryFUToken
std::vector< TriggerObject > TriggerObjectCollection
collection of trigger physics objects (e.g., all isolated muons)
virtual void fillMEforMonTriggerSummary(const edm::Event &iEvent, const edm::EventSetup &)
std::vector< std::pair< std::string, std::string > > custompathnamepairs_
edm::EDGetTokenT< reco::PFMETCollection > pfMetToken
std::vector< size_type > Keys
PathInfoCollection hltPathsAll_
bool init(const edm::Run &iRun, const edm::EventSetup &iSetup, const std::string &processName, bool &changed)
d'tor
static const JetCorrector * getJetCorrector(const std::string &fName, const edm::EventSetup &fSetup)
retrieve corrector from the event setup. troughs exception if something is missing ...
const std::string getL1ConditionModuleName(const std::string &pathname)
edm::Handle< trigger::TriggerEvent > triggerObj_
void endLuminosityBlock(const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)
DQM Client Diagnostic.
MonitorElement * correlation_AllWrtMB
T const * product() const
edm::EDGetTokenT< reco::CaloJetCollection > caloJetsToken
PathInfoCollection hltPathsEffWrtMu_
void endRun(const edm::Run &run, const edm::EventSetup &c)
void beginRun(const edm::Run &run, const edm::EventSetup &c)
HLTConfigProvider hltConfig_
MonitorElement * correlation_All
virtual void fillMEforEffWrtMuTrigger(const edm::Event &iEvent, const edm::EventSetup &)
edm::Handle< reco::CaloMETCollection > calometColl_
edm::EDGetTokenT< edm::TriggerResults > triggerResultsToken
edm::EDGetTokenT< trigger::TriggerEvent > triggerSummaryToken
virtual void analyze(const edm::Event &, const edm::EventSetup &)
virtual void fillMEforEffAllTrigger(const edm::Event &iEvent, const edm::EventSetup &)
MonitorElement * rate_AllWrtMB
volatile std::atomic< bool > shutdown_flag false
MonitorElement * book2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
Book 2D histogram.
virtual float pt() const GCC11_FINAL
transverse momentum
MonitorElement * rate_AllWrtMu
reco::helper::JetIDHelper * jetID
void setCurrentFolder(const std::string &fullpath)
void calculate(const edm::Event &event, const reco::CaloJet &jet, const int iDbg=0)
MonitorElement * rate_All
PathInfoCollection hltPathsEff_