43 using namespace trigger;
48 LogDebug(
"JetMETHLTOfflineSource") <<
"constructor....";
52 LogDebug(
"JetMETHLTOfflineSource") <<
"unabel to get DQMStore service?";
101 std::vector<edm::ParameterSet>
paths = iConfig.
getParameter<std::vector<edm::ParameterSet> >(
"pathPairs");
102 for(std::vector<edm::ParameterSet>::iterator pathconf = paths.begin() ; pathconf != paths.end(); pathconf++) {
105 pathconf->getParameter<
std::string>(
"denompathname")
128 edm::LogInfo(
"FourVectorHLTOffline") <<
"TriggerResults not found, "
142 edm::LogInfo(
"CaloMETHLTOfflineSource") <<
"TriggerResults::HLT not found, "
143 "automatically select events";
153 edm::LogInfo(
"FourVectorHLTOffline") <<
"TriggerEvent not found, "
161 if(!ValidJetColl_)
return;
166 if(!ValidPFJetColl_)
return;
171 if(!ValidMETColl_)
return;
174 if(!ValidPFMETColl_)
return;
181 for(
int i=0;
i<2;
i++){
203 CaloJetCollection::const_iterator calojet_ =
calojet.begin();
204 for(; calojet_ !=
calojet.end(); ++calojet_){
244 PFJetCollection::const_iterator pfjet_ =
pfjet.begin();
245 for(; pfjet_ !=
pfjet.end(); ++pfjet_){
249 if(scale*pfjet_->pt()>
PFJetPt[0]){
260 PFJetPt[0] = scale*pfjet_->pt();
261 PFJetPx[0] = scale*pfjet_->px();
262 PFJetPy[0] = scale*pfjet_->py();
265 PFJetNHEF[0] = pfjet_->neutralHadronEnergyFraction();
266 PFJetCHEF[0] = pfjet_->chargedHadronEnergyFraction();
267 PFJetNEMF[0] = pfjet_->neutralEmEnergyFraction();
268 PFJetCEMF[0] = pfjet_->chargedEmEnergyFraction();
270 else if(scale*pfjet_->pt()<
PFJetPt[0] && scale*pfjet_->pt()>
PFJetPt[1] ){
271 PFJetPt[1] = scale*pfjet_->pt();
272 PFJetPx[1] = scale*pfjet_->px();
273 PFJetPy[1] = scale*pfjet_->py();
276 PFJetNHEF[1] = pfjet_->neutralHadronEnergyFraction();
277 PFJetCHEF[1] = pfjet_->chargedHadronEnergyFraction();
278 PFJetNEMF[1] = pfjet_->neutralEmEnergyFraction();
279 PFJetCEMF[1] = pfjet_->chargedEmEnergyFraction();
285 for(
int i = 0;
i<2;
i++){
289 for(
int i = 0;
i<2;
i++){
312 for (
unsigned int j=0;
j!=nPath; ++
j) {
317 if(
verbose_)
cout<<
"fillMEforMonTriggerSummary: Muon Match"<<endl;
328 for (
unsigned int j=0;
j!=nPath; ++
j) {
333 if(
verbose_)
cout<<
"fillMEforMonTriggerSummary: MinBias Match"<<endl;
343 bool trigFirst=
false;
346 if(!trigFirst)
continue;
359 for(PathInfoCollection::iterator
w = v+1;
w!=
hltPathsAll_.end(); ++
w ){
360 bool trigSec =
false;
363 if(trigSec && trigFirst){
368 if(!trigSec && trigFirst){
378 iEvent.
getByLabel (
"offlinePrimaryVertices",Vtx);
380 for (VertexCollection::const_iterator itv=Vtx->begin(); itv!=Vtx->end(); itv++){
402 v->getMEhisto_TriggerSummary()->Fill(0.);
405 bool l1found =
false;
406 if(l1Index < triggerObj_->sizeFilters() ) l1found =
true;
407 if(!l1found)
v->getMEhisto_TriggerSummary()->Fill(1.);
408 if(!l1found && !(
triggerResults_->accept(index)))
v->getMEhisto_TriggerSummary()->Fill(2.);
409 if(!l1found && (
triggerResults_->accept(index)))
v->getMEhisto_TriggerSummary()->Fill(3.);
410 if(l1found)
v->getMEhisto_TriggerSummary()->Fill(4.);
411 if(l1found && (
triggerResults_->accept(index)))
v->getMEhisto_TriggerSummary()->Fill(5.);
412 if(l1found && !(
triggerResults_->accept(index)))
v->getMEhisto_TriggerSummary()->Fill(6.);
416 CaloJetCollection::const_iterator
jet =
calojet.begin();
417 v->getMEhisto_JetPt()->Fill(jet->pt());
418 v->getMEhisto_EtavsPt()->Fill(jet->eta(),jet->pt());
419 v->getMEhisto_PhivsPt()->Fill(jet->phi(),jet->pt());
424 v->getMEhisto_JetSize()->Fill(
calojet.size());
426 CaloJetCollection::const_iterator
jet =
calojet.begin();
427 CaloJetCollection::const_iterator jet2=
calojet.begin(); jet2++;
430 CaloJetCollection::const_iterator jet3 = jet2++;
433 v->getMEhisto_Pt12()->Fill((jet->pt()+jet2->pt())/2.);
434 v->getMEhisto_Eta12()->Fill((jet->eta()+jet2->eta())/2.);
435 v->getMEhisto_Phi12()->Fill(
deltaPhi(jet->phi(),jet2->phi()));
436 v->getMEhisto_Pt3()->Fill(jet3pt);
437 v->getMEhisto_Pt12Pt3()->Fill((jet->pt()+jet2->pt())/2., jet3pt);
438 v->getMEhisto_Pt12Phi12()->Fill((jet->pt()+jet2->pt())/2.,
deltaPhi(jet->phi(),jet2->phi()));
442 if(((
v->getTriggerType().compare(
"MET_Trigger") == 0)|| (
v->getTriggerType().compare(
"TET_Trigger") == 0)) &&
calometColl_.
isValid() ){
445 v->getMEhisto_JetPt()->Fill(met.
pt());
465 std::vector<double>jetPtVec;
466 std::vector<double>jetPhiVec;
467 std::vector<double>jetEtaVec;
468 std::vector<double>jetPxVec;
469 std::vector<double>jetPyVec;
470 std::vector<double>hltPtVec;
471 std::vector<double>hltPhiVec;
472 std::vector<double>hltEtaVec;
473 std::vector<double>hltPxVec;
474 std::vector<double>hltPyVec;
485 bool hltTrigBool =
false;
486 bool diJetFire =
false;
490 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index "<< l1Index <<
" of that name "<<
l1Tag;
497 && v->getTriggerType().compare(
"SingleJet_Trigger") == 0)
498 v->getMEhisto_N_L1()->Fill(kl1.size());
500 trigger::Keys::const_iterator ki = kl1.begin();
501 for(; ki != kl1.end(); ++ki){
502 double l1TrigEta = -100;
503 double l1TrigPhi = -100;
506 l1TrigEta = toc[*ki].eta();
507 l1TrigPhi = toc[*ki].phi();
508 if(v->getTriggerType().compare(
"SingleJet_Trigger") == 0){
509 v->getMEhisto_Pt_L1()->Fill(toc[*ki].pt());
510 if (
isBarrel(toc[*ki].
eta())) v->getMEhisto_PtBarrel_L1()->Fill(toc[*ki].pt());
511 if (
isEndCap(toc[*ki].
eta())) v->getMEhisto_PtEndcap_L1()->Fill(toc[*ki].pt());
512 if (
isForward(toc[*ki].
eta())) v->getMEhisto_PtForward_L1()->Fill(toc[*ki].pt());
513 v->getMEhisto_Eta_L1()->Fill(toc[*ki].
eta());
514 v->getMEhisto_Phi_L1()->Fill(toc[*ki].
phi());
515 v->getMEhisto_EtaPhi_L1()->Fill(toc[*ki].
eta(),toc[*ki].
phi());
520 v->getMEhisto_Pt_L1()->Fill(toc[*ki].pt());
521 v->getMEhisto_Phi_L1()->Fill(toc[*ki].
phi());
526 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
532 && v->getTriggerType().compare(
"SingleJet_Trigger") == 0)
533 v->getMEhisto_N_HLT()->Fill(khlt.size());
535 trigger::Keys::const_iterator kj = khlt.begin();
537 for(;kj != khlt.end(); ++kj){
539 double hltTrigEta = -100;
540 double hltTrigPhi = -100;
541 hltTrigEta = toc[*kj].eta();
542 hltTrigPhi = toc[*kj].phi();
543 if((
deltaR(hltTrigEta, hltTrigPhi, l1TrigEta, l1TrigPhi)) < 0.4
544 && (v->getTriggerType().compare(
"DiJet_Trigger") == 0))
550 for(;kj != khlt.end(); ++kj){
551 double hltTrigEta = -100.;
552 double hltTrigPhi = -100.;
556 v->getMEhisto_Pt_HLT()->Fill(toc[*kj].pt());
557 v->getMEhisto_Phi_HLT()->Fill(toc[*kj].
phi());
558 v->getMEhisto_PtCorrelation_L1HLT()->Fill(toc[*ki].pt(),toc[*kj].pt());
559 v->getMEhisto_PhiCorrelation_L1HLT()->Fill(toc[*ki].
phi(),toc[*kj].
phi());
560 v->getMEhisto_PtResolution_L1HLT()->Fill((toc[*ki].pt()-toc[*kj].pt())/(toc[*ki].pt()));
561 v->getMEhisto_PhiResolution_L1HLT()->Fill((toc[*ki].
phi()-toc[*kj].
phi())/(toc[*ki].
phi()));
565 hltTrigEta = toc[*kj].eta();
566 hltTrigPhi = toc[*kj].phi();
567 if((
deltaR(hltTrigEta, hltTrigPhi, l1TrigEta, l1TrigPhi)) < 0.4){
568 if(v->getTriggerType().compare(
"SingleJet_Trigger") == 0){
569 v->getMEhisto_PtCorrelation_L1HLT()->Fill(toc[*ki].pt(),toc[*kj].pt());
570 v->getMEhisto_EtaCorrelation_L1HLT()->Fill(toc[*ki].
eta(),toc[*kj].
eta());
571 v->getMEhisto_PhiCorrelation_L1HLT()->Fill(toc[*ki].
phi(),toc[*kj].
phi());
572 v->getMEhisto_PtResolution_L1HLT()->Fill((toc[*ki].pt()-toc[*kj].pt())/(toc[*ki].pt()));
573 v->getMEhisto_EtaResolution_L1HLT()->Fill((toc[*ki].
eta()-toc[*kj].
eta())/(toc[*ki].
eta()));
574 v->getMEhisto_PhiResolution_L1HLT()->Fill((toc[*ki].
phi()-toc[*kj].
phi())/(toc[*ki].
phi()));
577 if(((
deltaR(hltTrigEta, hltTrigPhi, l1TrigEta, l1TrigPhi) < 0.4 )
578 || ((v->getTriggerType().compare(
"DiJet_Trigger") == 0) && hltTrigBool)) && !diJetFire){
579 if(v->getTriggerType().compare(
"SingleJet_Trigger") == 0){
580 v->getMEhisto_Pt_HLT()->Fill(toc[*kj].pt());
581 if (
isBarrel(toc[*kj].
eta())) v->getMEhisto_PtBarrel_HLT()->Fill(toc[*kj].pt());
582 if (
isEndCap(toc[*kj].
eta())) v->getMEhisto_PtEndcap_HLT()->Fill(toc[*kj].pt());
583 if (
isForward(toc[*kj].
eta())) v->getMEhisto_PtForward_HLT()->Fill(toc[*kj].pt());
584 v->getMEhisto_Eta_HLT()->Fill(toc[*kj].
eta());
585 v->getMEhisto_Phi_HLT()->Fill(toc[*kj].
phi());
586 v->getMEhisto_EtaPhi_HLT()->Fill(toc[*kj].
eta(),toc[*kj].
phi());
592 && (v->getPath().compare(
"PFJet") == 0)){
595 for(
int iCalo=0; iCalo<2; iCalo++){
598 if(v->getTriggerType().compare(
"SingleJet_Trigger") == 0){
599 v->getMEhisto_Pt()->Fill(
CaloJetPt[iCalo]);
608 v->getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*kj].pt(),
CaloJetPt[iCalo]);
609 v->getMEhisto_EtaCorrelation_HLTRecObj()->Fill(toc[*kj].
eta(),
CaloJetEta[iCalo]);
610 v->getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*kj].
phi(),
CaloJetPhi[iCalo]);
612 v->getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*kj].pt()-
CaloJetPt[iCalo])/(toc[*kj].pt()));
613 v->getMEhisto_EtaResolution_HLTRecObj()->Fill((toc[*kj].
eta()-
CaloJetEta[iCalo])/(toc[*kj].
eta()));
614 v->getMEhisto_PhiResolution_HLTRecObj()->Fill((toc[*kj].
phi()-
CaloJetPhi[iCalo])/(toc[*kj].
phi()));
618 if((v->getTriggerType().compare(
"DiJet_Trigger") == 0)){
625 hltPhiVec.push_back(toc[*kj].
phi());
626 hltPtVec.push_back(toc[*kj].pt());
627 hltEtaVec.push_back(toc[*kj].
eta());
628 hltPxVec.push_back(toc[*kj].px());
629 hltPyVec.push_back(toc[*kj].py());
638 && (v->getPath().compare(
"PFJet") != 0)){
641 for(
int iPF=0; iPF<2; iPF++){
644 if(v->getTriggerType().compare(
"SingleJet_Trigger") == 0){
645 v->getMEhisto_Pt()->Fill(
PFJetPt[iPF]);
650 v->getMEhisto_Eta()->Fill(
PFJetEta[iPF]);
651 v->getMEhisto_Phi()->Fill(
PFJetPhi[iPF]);
654 v->getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*kj].pt(),
PFJetPt[iPF]);
655 v->getMEhisto_EtaCorrelation_HLTRecObj()->Fill(toc[*kj].
eta(),
PFJetEta[iPF]);
656 v->getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*kj].
phi(),
PFJetPhi[iPF]);
658 v->getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*kj].pt()-
PFJetPt[iPF])/(toc[*kj].pt()));
659 v->getMEhisto_EtaResolution_HLTRecObj()->Fill((toc[*kj].
eta()-
PFJetEta[iPF])/(toc[*kj].
eta()));
660 v->getMEhisto_PhiResolution_HLTRecObj()->Fill((toc[*kj].
phi()-
PFJetPhi[iPF])/(toc[*kj].
phi()));
664 if((v->getTriggerType().compare(
"DiJet_Trigger") == 0)){
666 jetPtVec.push_back(
PFJetPt[iPF]);
668 jetPxVec.push_back(
PFJetPx[iPF]);
669 jetPyVec.push_back(
PFJetPy[iPF]);
671 hltPhiVec.push_back(toc[*kj].
phi());
672 hltPtVec.push_back(toc[*kj].pt());
673 hltEtaVec.push_back(toc[*kj].
eta());
674 hltPxVec.push_back(toc[*kj].px());
675 hltPyVec.push_back(toc[*kj].py());
687 && (v->getPath().find(
"HLT_PFMET")==std::string::npos)){
691 v->getMEhisto_Pt()->Fill(met.
et());
692 v->getMEhisto_Phi()->Fill(met.
phi());
694 v->getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*kj].et(),met.
et());
695 v->getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*kj].
phi(),met.
phi());
696 v->getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*kj].et()-met.
et())/(toc[*kj].et()));
697 v->getMEhisto_PhiResolution_HLTRecObj()->Fill((toc[*kj].
phi()-met.
phi())/(toc[*kj].
phi()));
703 && (v->getPath().find(
"HLT_PFMET")!=std::string::npos)){
705 const PFMET pfmet = pfmetcol->front();
707 v->getMEhisto_Pt()->Fill(pfmet.
et());
708 v->getMEhisto_Phi()->Fill(pfmet.
phi());
710 v->getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*kj].et(),pfmet.
et());
711 v->getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*kj].
phi(),pfmet.
phi());
712 v->getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*kj].et()-pfmet.
et())/(toc[*kj].et()));
713 v->getMEhisto_PhiResolution_HLTRecObj()->Fill((toc[*kj].
phi()-pfmet.
phi())/(toc[*kj].
phi()));
716 if((v->getTriggerType().compare(
"DiJet_Trigger") == 0)) diJetFire =
true;
720 v->getMEhisto_N()->Fill(jetsize);
723 if((v->getTriggerType().compare(
"DiJet_Trigger") == 0) && jetPtVec.size() >1){
724 double AveJetPt = (jetPtVec[0] + jetPtVec[1])/2;
725 double AveJetEta = (jetEtaVec[0] + jetEtaVec[1])/2;
726 double JetDelPhi =
deltaPhi(jetPhiVec[0],jetPhiVec[1]);
727 double AveHLTPt = (hltPtVec[0] + hltPtVec[1])/2;
728 double AveHLTEta = (hltEtaVec[0] + hltEtaVec[1])/2;
729 double HLTDelPhi =
deltaPhi(hltPhiVec[0],hltPhiVec[1]);
730 v->getMEhisto_AveragePt_RecObj()->Fill(AveJetPt);
731 v->getMEhisto_AverageEta_RecObj()->Fill(AveJetEta);
732 v->getMEhisto_DeltaPhi_RecObj()->Fill(JetDelPhi);
734 v->getMEhisto_AveragePt_HLTObj()->Fill(AveHLTPt);
735 v->getMEhisto_AverageEta_HLTObj()->Fill(AveHLTEta);
736 v->getMEhisto_DeltaPhi_HLTObj()->Fill(HLTDelPhi);
808 for (
unsigned int j=0;
j!=nPath; ++
j) {
813 if(
verbose_)
cout<<
"fillMEforMonAllTriggerwrtMuonTrigger: Muon Match"<<endl;
828 std::vector<double>jetPtVec;
829 std::vector<double>jetPhiVec;
830 std::vector<double>jetEtaVec;
831 std::vector<double>jetPxVec;
832 std::vector<double>jetPyVec;
833 std::vector<double>hltPtVec;
834 std::vector<double>hltPhiVec;
835 std::vector<double>hltEtaVec;
836 std::vector<double>hltPxVec;
837 std::vector<double>hltPyVec;
840 bool fillL1HLT =
false;
848 bool hltTrigBool =
false;
849 bool diJetFire =
false;
853 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index "<< l1Index <<
" of that name "<<
l1Tag;
859 trigger::Keys::const_iterator ki = kl1.begin();
860 for(; ki != kl1.end(); ++ki){
861 double l1TrigEta = -100;
862 double l1TrigPhi = -100;
866 l1TrigEta = toc[*ki].eta();
867 l1TrigPhi = toc[*ki].phi();
868 v->getMEhisto_Pt_L1()->Fill(toc[*ki].pt());
869 if (
isBarrel(toc[*ki].
eta())) v->getMEhisto_PtBarrel_L1()->Fill(toc[*ki].pt());
870 if (
isEndCap(toc[*ki].
eta())) v->getMEhisto_PtEndcap_L1()->Fill(toc[*ki].pt());
871 if (
isForward(toc[*ki].
eta())) v->getMEhisto_PtForward_L1()->Fill(toc[*ki].pt());
872 v->getMEhisto_Eta_L1()->Fill(toc[*ki].
eta());
873 v->getMEhisto_Phi_L1()->Fill(toc[*ki].
phi());
874 v->getMEhisto_EtaPhi_L1()->Fill(toc[*ki].
eta(),toc[*ki].
phi());
877 v->getMEhisto_Pt_L1()->Fill(toc[*ki].pt());
878 v->getMEhisto_Phi_L1()->Fill(toc[*ki].
phi());
883 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
888 v->getMEhisto_N_HLT()->Fill(khlt.size());
889 trigger::Keys::const_iterator kj = khlt.begin();
891 for(;kj != khlt.end(); ++kj){
893 double hltTrigEta = -100;
894 double hltTrigPhi = -100;
895 hltTrigEta = toc[*kj].eta();
896 hltTrigPhi = toc[*kj].phi();
897 if((
deltaR(hltTrigEta, hltTrigPhi, l1TrigEta, l1TrigPhi)) < 0.4
898 && (v->getTriggerType().compare(
"DiJet_Trigger") == 0))
904 for(;kj != khlt.end(); ++kj){
905 double hltTrigEta = -100.;
906 double hltTrigPhi = -100.;
910 v->getMEhisto_Pt_HLT()->Fill(toc[*kj].pt());
911 v->getMEhisto_Phi_HLT()->Fill(toc[*kj].
phi());
912 v->getMEhisto_PtCorrelation_L1HLT()->Fill(toc[*ki].pt(),toc[*kj].pt());
913 v->getMEhisto_PhiCorrelation_L1HLT()->Fill(toc[*ki].
phi(),toc[*kj].
phi());
914 v->getMEhisto_PtResolution_L1HLT()->Fill((toc[*ki].pt()-toc[*kj].pt())/(toc[*ki].pt()));
915 v->getMEhisto_PhiResolution_L1HLT()->Fill((toc[*ki].
phi()-toc[*kj].
phi())/(toc[*ki].
phi()));
919 hltTrigEta = toc[*kj].eta();
920 hltTrigPhi = toc[*kj].phi();
921 if((
deltaR(hltTrigEta, hltTrigPhi, l1TrigEta, l1TrigPhi)) < 0.4){
922 v->getMEhisto_PtCorrelation_L1HLT()->Fill(toc[*ki].pt(),toc[*kj].pt());
923 v->getMEhisto_EtaCorrelation_L1HLT()->Fill(toc[*ki].
eta(),toc[*kj].
eta());
924 v->getMEhisto_PhiCorrelation_L1HLT()->Fill(toc[*ki].
phi(),toc[*kj].
phi());
925 v->getMEhisto_PtResolution_L1HLT()->Fill((toc[*ki].pt()-toc[*kj].pt())/(toc[*ki].pt()));
926 v->getMEhisto_EtaResolution_L1HLT()->Fill((toc[*ki].
eta()-toc[*kj].
eta())/(toc[*ki].
eta()));
927 v->getMEhisto_PhiResolution_L1HLT()->Fill((toc[*ki].
phi()-toc[*kj].
phi())/(toc[*ki].
phi()));
929 if(((
deltaR(hltTrigEta, hltTrigPhi, l1TrigEta, l1TrigPhi) < 0.4 )
930 || ((v->getTriggerType().compare(
"DiJet_Trigger") == 0) && hltTrigBool)) && !diJetFire){
931 v->getMEhisto_Pt_HLT()->Fill(toc[*kj].pt());
932 if (
isBarrel(toc[*kj].
eta())) v->getMEhisto_PtBarrel_HLT()->Fill(toc[*kj].pt());
933 if (
isEndCap(toc[*kj].
eta())) v->getMEhisto_PtEndcap_HLT()->Fill(toc[*kj].pt());
934 if (
isForward(toc[*kj].
eta())) v->getMEhisto_PtForward_HLT()->Fill(toc[*kj].pt());
935 v->getMEhisto_Eta_HLT()->Fill(toc[*kj].
eta());
936 v->getMEhisto_Phi_HLT()->Fill(toc[*kj].
phi());
937 v->getMEhisto_EtaPhi_HLT()->Fill(toc[*kj].
eta(),toc[*kj].
phi());
942 && (v->getPath().compare(
"PFJet") == 0)){
945 for(
int iCalo=0; iCalo<2; iCalo++){
948 if(v->getTriggerType().compare(
"SingleJet_Trigger") == 0){
949 v->getMEhisto_Pt()->Fill(
CaloJetPt[iCalo]);
958 v->getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*kj].pt(),
CaloJetPt[iCalo]);
959 v->getMEhisto_EtaCorrelation_HLTRecObj()->Fill(toc[*kj].
eta(),
CaloJetEta[iCalo]);
960 v->getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*kj].
phi(),
CaloJetPhi[iCalo]);
962 v->getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*kj].pt()-
CaloJetPt[iCalo])/(toc[*kj].pt()));
963 v->getMEhisto_EtaResolution_HLTRecObj()->Fill((toc[*kj].
eta()-
CaloJetEta[iCalo])/(toc[*kj].
eta()));
964 v->getMEhisto_PhiResolution_HLTRecObj()->Fill((toc[*kj].
phi()-
CaloJetPhi[iCalo])/(toc[*kj].
phi()));
968 if((v->getTriggerType().compare(
"DiJet_Trigger") == 0)){
975 hltPhiVec.push_back(toc[*kj].
phi());
976 hltPtVec.push_back(toc[*kj].pt());
977 hltEtaVec.push_back(toc[*kj].
eta());
978 hltPxVec.push_back(toc[*kj].px());
979 hltPyVec.push_back(toc[*kj].py());
988 && (v->getPath().compare(
"PFJet") != 0)){
991 for(
int iPF=0; iPF<2; iPF++){
994 if(v->getTriggerType().compare(
"SingleJet_Trigger") == 0){
995 v->getMEhisto_Pt()->Fill(
PFJetPt[iPF]);
1000 v->getMEhisto_Eta()->Fill(
PFJetEta[iPF]);
1001 v->getMEhisto_Phi()->Fill(
PFJetPhi[iPF]);
1004 v->getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*kj].pt(),
PFJetPt[iPF]);
1005 v->getMEhisto_EtaCorrelation_HLTRecObj()->Fill(toc[*kj].
eta(),
PFJetEta[iPF]);
1006 v->getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*kj].
phi(),
PFJetPhi[iPF]);
1008 v->getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*kj].pt()-
PFJetPt[iPF])/(toc[*kj].pt()));
1009 v->getMEhisto_EtaResolution_HLTRecObj()->Fill((toc[*kj].
eta()-
PFJetEta[iPF])/(toc[*kj].
eta()));
1010 v->getMEhisto_PhiResolution_HLTRecObj()->Fill((toc[*kj].
phi()-
PFJetPhi[iPF])/(toc[*kj].
phi()));
1014 if((v->getTriggerType().compare(
"DiJet_Trigger") == 0)){
1015 jetPhiVec.push_back(
PFJetPhi[iPF]);
1016 jetPtVec.push_back(
PFJetPt[iPF]);
1017 jetEtaVec.push_back(
PFJetEta[iPF]);
1018 jetPxVec.push_back(
PFJetPx[iPF]);
1019 jetPyVec.push_back(
PFJetPy[iPF]);
1021 hltPhiVec.push_back(toc[*kj].
phi());
1022 hltPtVec.push_back(toc[*kj].pt());
1023 hltEtaVec.push_back(toc[*kj].
eta());
1024 hltPxVec.push_back(toc[*kj].px());
1025 hltPyVec.push_back(toc[*kj].py());
1037 && (v->getPath().find(
"HLT_PFMET")==std::string::npos)){
1041 v->getMEhisto_Pt()->Fill(met.
et());
1042 v->getMEhisto_Phi()->Fill(met.
phi());
1044 v->getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*kj].et(),met.
et());
1045 v->getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*kj].
phi(),met.
phi());
1046 v->getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*kj].et()-met.
et())/(toc[*kj].et()));
1047 v->getMEhisto_PhiResolution_HLTRecObj()->Fill((toc[*kj].
phi()-met.
phi())/(toc[*kj].
phi()));
1053 && (v->getPath().find(
"HLT_PFMET")!=std::string::npos)){
1055 const PFMET pfmet = pfmetcol->front();
1057 v->getMEhisto_Pt()->Fill(pfmet.
et());
1058 v->getMEhisto_Phi()->Fill(pfmet.
phi());
1060 v->getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*kj].et(),pfmet.
et());
1061 v->getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*kj].
phi(),pfmet.
phi());
1062 v->getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*kj].et()-pfmet.
et())/(toc[*kj].et()));
1063 v->getMEhisto_PhiResolution_HLTRecObj()->Fill((toc[*kj].
phi()-pfmet.
phi())/(toc[*kj].
phi()));
1066 if((v->getTriggerType().compare(
"DiJet_Trigger") == 0)) diJetFire =
true;
1070 v->getMEhisto_N()->Fill(jetsize);
1073 if((v->getTriggerType().compare(
"DiJet_Trigger") == 0) && jetPtVec.size() >1){
1074 double AveJetPt = (jetPtVec[0] + jetPtVec[1])/2;
1075 double AveJetEta = (jetEtaVec[0] + jetEtaVec[1])/2;
1076 double JetDelPhi =
deltaPhi(jetPhiVec[0],jetPhiVec[1]);
1077 double AveHLTPt = (hltPtVec[0] + hltPtVec[1])/2;
1078 double AveHLTEta = (hltEtaVec[0] + hltEtaVec[1])/2;
1079 double HLTDelPhi =
deltaPhi(hltPhiVec[0],hltPhiVec[1]);
1080 v->getMEhisto_AveragePt_RecObj()->Fill(AveJetPt);
1081 v->getMEhisto_AverageEta_RecObj()->Fill(AveJetEta);
1082 v->getMEhisto_DeltaPhi_RecObj()->Fill(JetDelPhi);
1084 v->getMEhisto_AveragePt_HLTObj()->Fill(AveHLTPt);
1085 v->getMEhisto_AverageEta_HLTObj()->Fill(AveHLTEta);
1086 v->getMEhisto_DeltaPhi_HLTObj()->Fill(HLTDelPhi);
1090 if(v->getPath().find(
"L1") != std::string::npos && !fillL1HLT){
1092 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index "<< l1Index <<
" of that name "<<
l1Tag;
1097 for( trigger::Keys::const_iterator ki = kl1.begin(); ki != kl1.end(); ++ki){
1098 double l1TrigEta = toc[*ki].eta();
1099 double l1TrigPhi = toc[*ki].phi();
1104 if(
deltaR(l1TrigEta, l1TrigPhi, jetEta, jetPhi) < 0.4){
1106 v->getMEhisto_Pt()->Fill(
jet->pt());
1107 if (
isBarrel(
jet->eta())) v->getMEhisto_PtBarrel()->Fill(
jet->pt());
1108 if (
isEndCap(
jet->eta())) v->getMEhisto_PtEndcap()->Fill(
jet->pt());
1109 if (
isForward(
jet->eta())) v->getMEhisto_PtForward()->Fill(
jet->pt());
1111 v->getMEhisto_Eta()->Fill(
jet->eta());
1112 v->getMEhisto_Phi()->Fill(
jet->phi());
1113 v->getMEhisto_EtaPhi()->Fill(
jet->eta(),
jet->phi());
1115 v->getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*ki].pt(),
jet->pt());
1116 v->getMEhisto_EtaCorrelation_HLTRecObj()->Fill(toc[*ki].
eta(),
jet->eta());
1117 v->getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*ki].
phi(),
jet->phi());
1119 v->getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*ki].pt()-
jet->pt())/(toc[*ki].pt()));
1120 v->getMEhisto_EtaResolution_HLTRecObj()->Fill((toc[*ki].
eta()-
jet->eta())/(toc[*ki].
eta()));
1121 v->getMEhisto_PhiResolution_HLTRecObj()->Fill((toc[*ki].
phi()-
jet->phi())/(toc[*ki].
phi()));
1126 v->getMEhisto_N()->Fill(jetsize);
1132 v->getMEhisto_Pt()->Fill(met.
pt());
1133 v->getMEhisto_Phi()->Fill(met.
phi());
1135 v->getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*ki].pt(),met.
pt());
1136 v->getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*ki].
phi(),met.
phi());
1137 v->getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*ki].pt()-met.
pt())/(toc[*ki].pt()));
1138 v->getMEhisto_PhiResolution_HLTRecObj()->Fill((toc[*ki].
phi()-met.
phi())/(toc[*ki].
phi()));
1156 bool denompassed =
false;
1157 bool numpassed =
false;
1163 denompassed =
false;
1172 if(denompassed==
false)
continue;
1183 double trigLowpTcut = 0;
1184 double trigMedpTcut = 0;
1185 double trigHighpTcut = 0;
1186 double trigLowpTcutFwd = 0;
1187 double trigMedpTcutFwd = 0;
1188 double trigHighpTcutFwd = 0;
1191 double trigLowpTPFcut = 0;
1192 double trigMedpTPFcut = 0;
1193 double trigHighpTPFcut = 0;
1194 double trigLowpTPFcutFwd = 0;
1195 double trigMedpTPFcutFwd = 0;
1196 double trigHighpTPFcutFwd = 0;
1199 size_t jetstrfound =
v->getPath().find(
"Jet");
1201 string tpath =
v->getPath();
1205 if(jetstrfound != string::npos){
1207 for(
int trig =
int(jetstrfound)+3; trig < int(jetstrfound)+7; trig++){
1208 if(!isdigit(tpath[trig]))
break;
1209 jetTrigVal+=tpath[trig];
1211 char *cjetTrigVal = (
char*)jetTrigVal.c_str();
1212 jetVal=atof(cjetTrigVal);
1217 trigMedpTcut = jetVal + 5.;
1218 trigHighpTcut = jetVal + 10.;
1220 trigLowpTcutFwd = jetVal + 9.;
1221 trigMedpTcutFwd = jetVal + 15.;
1222 trigHighpTcutFwd = jetVal + 21.;
1226 trigMedpTcut = jetVal + 2.;
1227 trigHighpTcut = jetVal + 60.;
1229 trigLowpTcutFwd = jetVal + 22.;
1230 trigMedpTcutFwd = jetVal + 25.;
1231 trigHighpTcutFwd = jetVal + 110.;
1233 trigLowpTcut = jetVal;
1239 trigMedpTPFcut = jetVal + 20.;
1240 trigHighpTPFcut = jetVal + 40.;
1242 trigLowpTPFcutFwd = jetVal + 60.;
1243 trigMedpTPFcutFwd = jetVal + 80.;
1244 trigHighpTPFcutFwd = jetVal + 100.;
1248 trigMedpTPFcut = jetVal + 40.;
1249 trigHighpTPFcut = jetVal + 140.;
1251 trigLowpTPFcutFwd = jetVal + 110.;
1252 trigMedpTPFcutFwd = jetVal + 130.;
1253 trigHighpTPFcutFwd = jetVal + 190.;
1255 trigLowpTPFcut = jetVal;
1261 if(
verbose_)
std::cout <<
"fillMEforEffAllTrigger: CaloJet -------------------" << std::endl;
1266 bool jetIDbool =
false;
1273 if((
v->getTriggerType().compare(
"SingleJet_Trigger") == 0) &&
calojet.size()){
1275 if(
verbose_)
cout<<
"Passed CaloJet ID -------------------" << endl;
1278 v->getMEhisto_DenominatorPt()->Fill(leadjpt);
1279 if (
isBarrel(leadjeta))
v->getMEhisto_DenominatorPtBarrel()->Fill(leadjpt);
1280 if (
isEndCap(leadjeta))
v->getMEhisto_DenominatorPtEndcap()->Fill(leadjpt);
1281 if (
isForward(leadjeta))
v->getMEhisto_DenominatorPtForward()->Fill(leadjpt);
1282 v->getMEhisto_DenominatorEta()->Fill(leadjeta);
1283 v->getMEhisto_DenominatorPhi()->Fill(leadjphi);
1284 v->getMEhisto_DenominatorEtaPhi()->Fill(leadjeta,leadjphi);
1286 v->getMEhisto_DenominatorEtaBarrel()->Fill(leadjeta);
1287 v->getMEhisto_DenominatorPhiBarrel()->Fill(leadjphi);
1290 v->getMEhisto_DenominatorEtaEndcap()->Fill(leadjeta);
1291 v->getMEhisto_DenominatorPhiEndcap()->Fill(leadjphi);
1294 v->getMEhisto_DenominatorEtaForward()->Fill(leadjeta);
1295 v->getMEhisto_DenominatorPhiForward()->Fill(leadjphi);
1297 if((leadjpt > trigLowpTcut && !
isForward(leadjeta)) || (leadjpt > trigLowpTcutFwd &&
isForward(leadjeta))){
1298 v->getMEhisto_DenominatorEta_LowpTcut()->Fill(leadjeta);
1299 v->getMEhisto_DenominatorPhi_LowpTcut()->Fill(leadjphi);
1300 v->getMEhisto_DenominatorEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
1302 if((leadjpt > trigMedpTcut && !
isForward(leadjeta)) || (leadjpt > trigMedpTcutFwd &&
isForward(leadjeta))){
1303 v->getMEhisto_DenominatorEta_MedpTcut()->Fill(leadjeta);
1304 v->getMEhisto_DenominatorPhi_MedpTcut()->Fill(leadjphi);
1305 v->getMEhisto_DenominatorEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
1307 if((leadjpt > trigHighpTcut && !
isForward(leadjeta)) || (leadjpt > trigHighpTcutFwd &&
isForward(leadjeta))){
1308 v->getMEhisto_DenominatorEta_HighpTcut()->Fill(leadjeta);
1309 v->getMEhisto_DenominatorPhi_HighpTcut()->Fill(leadjphi);
1310 v->getMEhisto_DenominatorEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
1316 double dRmin = 99999.;
1317 double dPhimin = 9999.;
1318 if(
v->getPath().find(
"L1") != std::string::npos){
1320 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
1324 for(trigger::Keys::const_iterator ki = kl1.begin();ki != kl1.end(); ++ki){
1325 double dR =
deltaR(toc[*ki].
eta(), toc[*ki].
phi(), leadjeta, leadjphi);
1334 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
1338 trigger::Keys::const_iterator kj = khlt.begin();
1339 for(;kj != khlt.end(); ++kj){
1341 leadjeta, leadjphi);
1351 v->getMEhisto_DeltaPhi()->Fill(dPhimin);
1352 v->getMEhisto_DeltaR()->Fill(dRmin);
1355 if(dRmin < 0.1 || (
v->getPath().find(
"L1") != std::string::npos && dRmin < 0.4)){
1356 v->getMEhisto_NumeratorPt()->Fill(leadjpt);
1357 if (
isBarrel(leadjeta))
v->getMEhisto_NumeratorPtBarrel()->Fill(leadjpt);
1358 if (
isEndCap(leadjeta))
v->getMEhisto_NumeratorPtEndcap()->Fill(leadjpt);
1359 if (
isForward(leadjeta))
v->getMEhisto_NumeratorPtForward()->Fill(leadjpt);
1360 v->getMEhisto_NumeratorEta()->Fill(leadjeta);
1361 v->getMEhisto_NumeratorPhi()->Fill(leadjphi);
1362 v->getMEhisto_NumeratorEtaPhi()->Fill(leadjeta,leadjphi);
1364 v->getMEhisto_NumeratorEtaBarrel()->Fill(leadjeta);
1365 v->getMEhisto_NumeratorPhiBarrel()->Fill(leadjphi);
1368 v->getMEhisto_NumeratorEtaEndcap()->Fill(leadjeta);
1369 v->getMEhisto_NumeratorPhiEndcap()->Fill(leadjphi);
1372 v->getMEhisto_NumeratorEtaForward()->Fill(leadjeta);
1373 v->getMEhisto_NumeratorPhiForward()->Fill(leadjphi);
1375 if((leadjpt > trigLowpTcut && !
isForward(leadjeta)) || (leadjpt > trigLowpTcutFwd &&
isForward(leadjeta))){
1376 v->getMEhisto_NumeratorEta_LowpTcut()->Fill(leadjeta);
1377 v->getMEhisto_NumeratorPhi_LowpTcut()->Fill(leadjphi);
1378 v->getMEhisto_NumeratorEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
1380 if((leadjpt > trigMedpTcut && !
isForward(leadjeta)) || (leadjpt > trigMedpTcutFwd &&
isForward(leadjeta))){
1381 v->getMEhisto_NumeratorEta_MedpTcut()->Fill(leadjeta);
1382 v->getMEhisto_NumeratorPhi_MedpTcut()->Fill(leadjphi);
1383 v->getMEhisto_NumeratorEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
1385 if((leadjpt > trigHighpTcut && !
isForward(leadjeta)) || (leadjpt > trigHighpTcutFwd &&
isForward(leadjeta))){
1386 v->getMEhisto_NumeratorEta_HighpTcut()->Fill(leadjeta);
1387 v->getMEhisto_NumeratorPhi_HighpTcut()->Fill(leadjphi);
1388 v->getMEhisto_NumeratorEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
1395 if(jetIDbool ==
true && (
v->getTriggerType().compare(
"DiJet_Trigger") == 0) &&
calojet.size()>1){
1400 if(numpassed==
true){
1409 if(
verbose_)
std::cout <<
"fillMEforEffAllTrigger: PFJet -------------------" << std::endl;
1414 bool jetIDbool =
false;
1433 if((
v->getTriggerType().compare(
"SingleJet_Trigger") == 0) &&
pfjet.size()){
1436 _pfMHT =
sqrt(pfMHTx*pfMHTx + pfMHTy*pfMHTy);
1437 v->getMEhisto_DenominatorPFMHT()->Fill(
_pfMHT);
1443 if(
verbose_)
cout<<
"Passed PFJet ID -------------------" << endl;
1445 v->getMEhisto_DenominatorPFPt()->Fill(leadjpt);
1446 if (
isBarrel(leadjeta))
v->getMEhisto_DenominatorPFPtBarrel()->Fill(leadjpt);
1447 if (
isEndCap(leadjeta))
v->getMEhisto_DenominatorPFPtEndcap()->Fill(leadjpt);
1448 if (
isForward(leadjeta))
v->getMEhisto_DenominatorPFPtForward()->Fill(leadjpt);
1449 v->getMEhisto_DenominatorPFEta()->Fill(leadjeta);
1450 v->getMEhisto_DenominatorPFPhi()->Fill(leadjphi);
1451 v->getMEhisto_DenominatorPFEtaPhi()->Fill(leadjeta,leadjphi);
1453 v->getMEhisto_DenominatorPFEtaBarrel()->Fill(leadjeta);
1454 v->getMEhisto_DenominatorPFPhiBarrel()->Fill(leadjphi);
1457 v->getMEhisto_DenominatorPFEtaEndcap()->Fill(leadjeta);
1458 v->getMEhisto_DenominatorPFPhiEndcap()->Fill(leadjphi);
1461 v->getMEhisto_DenominatorPFEtaForward()->Fill(leadjeta);
1462 v->getMEhisto_DenominatorPFPhiForward()->Fill(leadjphi);
1464 if((leadjpt > trigLowpTPFcut && !
isForward(leadjeta)) || (leadjpt > trigLowpTPFcutFwd &&
isForward(leadjeta))){
1465 v->getMEhisto_DenominatorPFEta_LowpTcut()->Fill(leadjeta);
1466 v->getMEhisto_DenominatorPFPhi_LowpTcut()->Fill(leadjphi);
1467 v->getMEhisto_DenominatorPFEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
1469 if((leadjpt > trigMedpTPFcut && !
isForward(leadjeta)) || (leadjpt > trigMedpTPFcutFwd &&
isForward(leadjeta))){
1470 v->getMEhisto_DenominatorPFEta_MedpTcut()->Fill(leadjeta);
1471 v->getMEhisto_DenominatorPFPhi_MedpTcut()->Fill(leadjphi);
1472 v->getMEhisto_DenominatorPFEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
1474 if((leadjpt > trigHighpTPFcut && !
isForward(leadjeta)) || (leadjpt > trigHighpTPFcutFwd &&
isForward(leadjeta))){
1475 v->getMEhisto_DenominatorPFEta_HighpTcut()->Fill(leadjeta);
1476 v->getMEhisto_DenominatorPFPhi_HighpTcut()->Fill(leadjphi);
1477 v->getMEhisto_DenominatorPFEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
1482 double dRmin = 99999.;
1483 double dPhimin = 9999.;
1484 if(
v->getPath().find(
"L1") != std::string::npos){
1486 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
1490 for(trigger::Keys::const_iterator ki = kl1.begin();ki != kl1.end(); ++ki){
1491 double dR =
deltaR(toc[*ki].
eta(), toc[*ki].
phi(), leadjeta, leadjphi);
1500 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
1504 for(trigger::Keys::const_iterator kj = khlt.begin();kj != khlt.end(); ++kj){
1505 double dR =
deltaR(toc[*kj].
eta(), toc[*kj].
phi(), leadjeta, leadjphi);
1514 v->getMEhisto_PFDeltaPhi()->Fill(dPhimin);
1515 v->getMEhisto_PFDeltaR()->Fill(dRmin);
1518 if(dRmin < 0.1 || (
v->getPath().find(
"L1") != std::string::npos && dRmin < 0.4)){
1519 v->getMEhisto_NumeratorPFPt()->Fill(leadjpt);
1520 if (
isBarrel(leadjeta))
v->getMEhisto_NumeratorPFPtBarrel()->Fill(leadjpt);
1521 if (
isEndCap(leadjeta))
v->getMEhisto_NumeratorPFPtEndcap()->Fill(leadjpt);
1522 if (
isForward(leadjeta))
v->getMEhisto_NumeratorPFPtForward()->Fill(leadjpt);
1523 v->getMEhisto_NumeratorPFEta()->Fill(leadjeta);
1524 v->getMEhisto_NumeratorPFPhi()->Fill(leadjphi);
1525 v->getMEhisto_NumeratorPFEtaPhi()->Fill(leadjeta,leadjphi);
1527 v->getMEhisto_NumeratorPFEtaBarrel()->Fill(leadjeta);
1528 v->getMEhisto_NumeratorPFPhiBarrel()->Fill(leadjphi);
1531 v->getMEhisto_NumeratorPFEtaEndcap()->Fill(leadjeta);
1532 v->getMEhisto_NumeratorPFPhiEndcap()->Fill(leadjphi);
1535 v->getMEhisto_NumeratorPFEtaForward()->Fill(leadjeta);
1536 v->getMEhisto_NumeratorPFPhiForward()->Fill(leadjphi);
1538 if((leadjpt > trigLowpTPFcut && !
isForward(leadjeta))
1539 || (leadjpt > trigLowpTPFcutFwd &&
isForward(leadjeta))){
1540 v->getMEhisto_NumeratorPFEta_LowpTcut()->Fill(leadjeta);
1541 v->getMEhisto_NumeratorPFPhi_LowpTcut()->Fill(leadjphi);
1542 v->getMEhisto_NumeratorPFEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
1544 if((leadjpt > trigMedpTPFcut && !
isForward(leadjeta))
1545 || (leadjpt > trigMedpTPFcutFwd &&
isForward(leadjeta))){
1546 v->getMEhisto_NumeratorPFEta_MedpTcut()->Fill(leadjeta);
1547 v->getMEhisto_NumeratorPFPhi_MedpTcut()->Fill(leadjphi);
1548 v->getMEhisto_NumeratorPFEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
1550 if((leadjpt > trigHighpTPFcut && !
isForward(leadjeta))
1551 || (leadjpt > trigHighpTPFcutFwd &&
isForward(leadjeta))){
1552 v->getMEhisto_NumeratorPFEta_HighpTcut()->Fill(leadjeta);
1553 v->getMEhisto_NumeratorPFPhi_HighpTcut()->Fill(leadjphi);
1554 v->getMEhisto_NumeratorPFEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
1560 if(jetIDbool ==
true && (
v->getTriggerType().compare(
"DiJet_Trigger") == 0) &&
pfjet.size()>1){
1580 if(
verbose_)
std::cout <<
"fillMEforEffAllTrigger: CaloMET -------------------" << std::endl;
1583 && (
v->getPath().find(
"HLT_PFMET")==std::string::npos)){
1586 v->getMEhisto_DenominatorPt()->Fill(met.
et());
1587 v->getMEhisto_DenominatorPhi()->Fill(met.
phi());
1589 v->getMEhisto_NumeratorPt()->Fill(met.
et());
1590 v->getMEhisto_NumeratorPhi()->Fill(met.
phi());
1592 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
1595 double dPhimin = 9999.;
1597 for(trigger::Keys::const_iterator kj = khlt.begin();kj != khlt.end(); ++kj){
1603 v->getMEhisto_DeltaPhi()->Fill(dPhimin);
1609 if(
verbose_)
std::cout <<
"fillMEforEffAllTrigger: PFMET -------------------" << std::endl;
1612 && (
v->getPath().find(
"HLT_PFMET")!=std::string::npos)){
1614 const PFMET met = pfmetcol->front();
1615 v->getMEhisto_DenominatorPt()->Fill(met.
et());
1616 v->getMEhisto_DenominatorPhi()->Fill(met.
phi());
1618 v->getMEhisto_NumeratorPt()->Fill(met.
et());
1619 v->getMEhisto_NumeratorPhi()->Fill(met.
phi());
1621 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
1624 double dPhimin = 9999.;
1626 for(trigger::Keys::const_iterator kj = khlt.begin();kj != khlt.end(); ++kj){
1632 v->getMEhisto_DeltaPhi()->Fill(dPhimin);
1655 bool muTrig =
false;
1656 bool numpassed =
false;
1660 for (
unsigned int j=0;
j!=nPath; ++
j) {
1665 if(
verbose_)
cout<<
"fillMEforEffWrtMuTrigger: Muon Match"<<endl;
1690 double trigLowpTcut = 0;
1691 double trigMedpTcut = 0;
1692 double trigHighpTcut = 0;
1693 double trigLowpTcutFwd = 0;
1694 double trigMedpTcutFwd = 0;
1695 double trigHighpTcutFwd = 0;
1698 double trigLowpTPFcut = 0;
1699 double trigMedpTPFcut = 0;
1700 double trigHighpTPFcut = 0;
1701 double trigLowpTPFcutFwd = 0;
1702 double trigMedpTPFcutFwd = 0;
1703 double trigHighpTPFcutFwd = 0;
1706 size_t jetstrfound = v->getPath().find(
"Jet");
1708 string tpath = v->getPath();
1712 if(jetstrfound != string::npos){
1714 for(
int trig =
int(jetstrfound)+3; trig < int(jetstrfound)+7; trig++){
1715 if(!isdigit(tpath[trig]))
break;
1716 jetTrigVal+=tpath[trig];
1718 char *cjetTrigVal = (
char*)jetTrigVal.c_str();
1719 jetVal=atof(cjetTrigVal);
1725 trigMedpTcut = jetVal + 5.;
1726 trigHighpTcut = jetVal + 10.;
1728 trigLowpTcutFwd = jetVal + 9.;
1729 trigMedpTcutFwd = jetVal + 15.;
1730 trigHighpTcutFwd = jetVal + 21.;
1734 trigMedpTcut = jetVal + 2.;
1735 trigHighpTcut = jetVal + 60.;
1737 trigLowpTcutFwd = jetVal + 22.;
1738 trigMedpTcutFwd = jetVal + 25.;
1739 trigHighpTcutFwd = jetVal + 110.;
1741 trigLowpTcut = jetVal;
1747 trigMedpTPFcut = jetVal + 20.;
1748 trigHighpTPFcut = jetVal + 40.;
1750 trigLowpTPFcutFwd = jetVal + 60.;
1751 trigMedpTPFcutFwd = jetVal + 80.;
1752 trigHighpTPFcutFwd = jetVal + 100.;
1756 trigMedpTPFcut = jetVal + 40.;
1757 trigHighpTPFcut = jetVal + 140.;
1759 trigLowpTPFcutFwd = jetVal + 110.;
1760 trigMedpTPFcutFwd = jetVal + 130.;
1761 trigHighpTPFcutFwd = jetVal + 190.;
1763 trigLowpTPFcut = jetVal;
1769 if(
verbose_)
std::cout <<
"fillMEforEffWrtMuTrigger: CaloJet -------------------" << std::endl;
1775 bool jetIDbool =
false;
1782 if((v->getTriggerType().compare(
"SingleJet_Trigger") == 0) &&
calojet.size()){
1784 if(
verbose_)
cout<<
"passed CaloJet ID -------------------" << endl;
1788 v->getMEhisto_DenominatorPt()->Fill(leadjpt);
1789 if (
isBarrel(leadjeta)) v->getMEhisto_DenominatorPtBarrel()->Fill(leadjpt);
1790 if (
isEndCap(leadjeta)) v->getMEhisto_DenominatorPtEndcap()->Fill(leadjpt);
1791 if (
isForward(leadjeta)) v->getMEhisto_DenominatorPtForward()->Fill(leadjpt);
1792 v->getMEhisto_DenominatorEta()->Fill(leadjeta);
1793 v->getMEhisto_DenominatorPhi()->Fill(leadjphi);
1794 v->getMEhisto_DenominatorEtaPhi()->Fill(leadjeta,leadjphi);
1796 v->getMEhisto_DenominatorEtaBarrel()->Fill(leadjeta);
1797 v->getMEhisto_DenominatorPhiBarrel()->Fill(leadjphi);
1800 v->getMEhisto_DenominatorEtaEndcap()->Fill(leadjeta);
1801 v->getMEhisto_DenominatorPhiEndcap()->Fill(leadjphi);
1804 v->getMEhisto_DenominatorEtaForward()->Fill(leadjeta);
1805 v->getMEhisto_DenominatorPhiForward()->Fill(leadjphi);
1807 if((leadjpt > trigLowpTcut && !
isForward(leadjeta)) || (leadjpt > trigLowpTcutFwd &&
isForward(leadjeta))){
1808 v->getMEhisto_DenominatorEta_LowpTcut()->Fill(leadjeta);
1809 v->getMEhisto_DenominatorPhi_LowpTcut()->Fill(leadjphi);
1810 v->getMEhisto_DenominatorEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
1812 if((leadjpt > trigMedpTcut && !
isForward(leadjeta)) || (leadjpt > trigMedpTcutFwd &&
isForward(leadjeta))){
1813 v->getMEhisto_DenominatorEta_MedpTcut()->Fill(leadjeta);
1814 v->getMEhisto_DenominatorPhi_MedpTcut()->Fill(leadjphi);
1815 v->getMEhisto_DenominatorEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
1817 if((leadjpt > trigHighpTcut && !
isForward(leadjeta)) || (leadjpt > trigHighpTcutFwd &&
isForward(leadjeta))){
1818 v->getMEhisto_DenominatorEta_HighpTcut()->Fill(leadjeta);
1819 v->getMEhisto_DenominatorPhi_HighpTcut()->Fill(leadjphi);
1820 v->getMEhisto_DenominatorEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
1826 double dRmin = 99999.;
1827 double dPhimin = 9999.;
1828 if(v->getPath().find(
"L1") != std::string::npos){
1830 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
1834 for(trigger::Keys::const_iterator ki = kl1.begin();ki != kl1.end(); ++ki){
1835 double dR =
deltaR(toc[*ki].
eta(), toc[*ki].
phi(), leadjeta, leadjphi);
1844 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
1848 trigger::Keys::const_iterator kj = khlt.begin();
1849 for(;kj != khlt.end(); ++kj){
1851 leadjeta, leadjphi);
1861 v->getMEhisto_DeltaPhi()->Fill(dPhimin);
1862 v->getMEhisto_DeltaR()->Fill(dRmin);
1865 if(dRmin < 0.1 || (v->getPath().find(
"L1") != std::string::npos && dRmin < 0.4)){
1866 v->getMEhisto_NumeratorPt()->Fill(leadjpt);
1867 if (
isBarrel(leadjeta)) v->getMEhisto_NumeratorPtBarrel()->Fill(leadjpt);
1868 if (
isEndCap(leadjeta)) v->getMEhisto_NumeratorPtEndcap()->Fill(leadjpt);
1869 if (
isForward(leadjeta)) v->getMEhisto_NumeratorPtForward()->Fill(leadjpt);
1870 v->getMEhisto_NumeratorEta()->Fill(leadjeta);
1871 v->getMEhisto_NumeratorPhi()->Fill(leadjphi);
1872 v->getMEhisto_NumeratorEtaPhi()->Fill(leadjeta,leadjphi);
1874 v->getMEhisto_NumeratorEtaBarrel()->Fill(leadjeta);
1875 v->getMEhisto_NumeratorPhiBarrel()->Fill(leadjphi);
1878 v->getMEhisto_NumeratorEtaEndcap()->Fill(leadjeta);
1879 v->getMEhisto_NumeratorPhiEndcap()->Fill(leadjphi);
1882 v->getMEhisto_NumeratorEtaForward()->Fill(leadjeta);
1883 v->getMEhisto_NumeratorPhiForward()->Fill(leadjphi);
1885 if((leadjpt > trigLowpTcut && !
isForward(leadjeta)) || (leadjpt > trigLowpTcutFwd &&
isForward(leadjeta))){
1886 v->getMEhisto_NumeratorEta_LowpTcut()->Fill(leadjeta);
1887 v->getMEhisto_NumeratorPhi_LowpTcut()->Fill(leadjphi);
1888 v->getMEhisto_NumeratorEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
1890 if((leadjpt > trigMedpTcut && !
isForward(leadjeta)) || (leadjpt > trigMedpTcutFwd &&
isForward(leadjeta))){
1891 v->getMEhisto_NumeratorEta_MedpTcut()->Fill(leadjeta);
1892 v->getMEhisto_NumeratorPhi_MedpTcut()->Fill(leadjphi);
1893 v->getMEhisto_NumeratorEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
1895 if((leadjpt > trigHighpTcut && !
isForward(leadjeta)) || (leadjpt > trigHighpTcutFwd &&
isForward(leadjeta))){
1896 v->getMEhisto_NumeratorEta_HighpTcut()->Fill(leadjeta);
1897 v->getMEhisto_NumeratorPhi_HighpTcut()->Fill(leadjphi);
1898 v->getMEhisto_NumeratorEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
1905 if(jetIDbool ==
true && (v->getTriggerType().compare(
"DiJet_Trigger") == 0) &&
calojet.size()>1){
1910 if(numpassed==
true){
1919 if(
verbose_)
std::cout <<
"fillMEforEffWrtMuTrigger: PFJet -------------------" << std::endl;
1925 bool jetIDbool =
false;
1944 if((v->getTriggerType().compare(
"SingleJet_Trigger") == 0) &&
pfjet.size()){
1947 _pfMHT =
sqrt(pfMHTx*pfMHTx + pfMHTy*pfMHTy);
1948 v->getMEhisto_DenominatorPFMHT()->Fill(
_pfMHT);
1954 if(
verbose_)
cout<<
"passed PFJet ID -------------------" << endl;
1956 v->getMEhisto_DenominatorPFPt()->Fill(leadjpt);
1957 if (
isBarrel(leadjeta)) v->getMEhisto_DenominatorPFPtBarrel()->Fill(leadjpt);
1958 if (
isEndCap(leadjeta)) v->getMEhisto_DenominatorPFPtEndcap()->Fill(leadjpt);
1959 if (
isForward(leadjeta)) v->getMEhisto_DenominatorPFPtForward()->Fill(leadjpt);
1960 v->getMEhisto_DenominatorPFEta()->Fill(leadjeta);
1961 v->getMEhisto_DenominatorPFPhi()->Fill(leadjphi);
1962 v->getMEhisto_DenominatorPFEtaPhi()->Fill(leadjeta,leadjphi);
1964 v->getMEhisto_DenominatorPFEtaBarrel()->Fill(leadjeta);
1965 v->getMEhisto_DenominatorPFPhiBarrel()->Fill(leadjphi);
1968 v->getMEhisto_DenominatorPFEtaEndcap()->Fill(leadjeta);
1969 v->getMEhisto_DenominatorPFPhiEndcap()->Fill(leadjphi);
1972 v->getMEhisto_DenominatorPFEtaForward()->Fill(leadjeta);
1973 v->getMEhisto_DenominatorPFPhiForward()->Fill(leadjphi);
1975 if((leadjpt > trigLowpTPFcut && !
isForward(leadjeta)) || (leadjpt > trigLowpTPFcutFwd &&
isForward(leadjeta))){
1976 v->getMEhisto_DenominatorPFEta_LowpTcut()->Fill(leadjeta);
1977 v->getMEhisto_DenominatorPFPhi_LowpTcut()->Fill(leadjphi);
1978 v->getMEhisto_DenominatorPFEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
1980 if((leadjpt > trigMedpTPFcut && !
isForward(leadjeta)) || (leadjpt > trigMedpTPFcutFwd &&
isForward(leadjeta))){
1981 v->getMEhisto_DenominatorPFEta_MedpTcut()->Fill(leadjeta);
1982 v->getMEhisto_DenominatorPFPhi_MedpTcut()->Fill(leadjphi);
1983 v->getMEhisto_DenominatorPFEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
1985 if((leadjpt > trigHighpTPFcut && !
isForward(leadjeta)) || (leadjpt > trigHighpTPFcutFwd &&
isForward(leadjeta))){
1986 v->getMEhisto_DenominatorPFEta_HighpTcut()->Fill(leadjeta);
1987 v->getMEhisto_DenominatorPFPhi_HighpTcut()->Fill(leadjphi);
1988 v->getMEhisto_DenominatorPFEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
1993 double dRmin = 99999.;
1994 double dPhimin = 9999.;
1995 if(v->getPath().find(
"L1") != std::string::npos){
1997 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
2001 for(trigger::Keys::const_iterator ki = kl1.begin();ki != kl1.end(); ++ki){
2002 double dR =
deltaR(toc[*ki].
eta(), toc[*ki].
phi(), leadjeta, leadjphi);
2011 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
2015 for(trigger::Keys::const_iterator kj = khlt.begin();kj != khlt.end(); ++kj){
2016 double dR =
deltaR(toc[*kj].
eta(), toc[*kj].
phi(), leadjeta, leadjphi);
2025 v->getMEhisto_PFDeltaPhi()->Fill(dPhimin);
2026 v->getMEhisto_PFDeltaR()->Fill(dRmin);
2029 if(dRmin < 0.1 || (v->getPath().find(
"L1") != std::string::npos && dRmin < 0.4)){
2030 v->getMEhisto_NumeratorPFPt()->Fill(leadjpt);
2031 if (
isBarrel(leadjeta)) v->getMEhisto_NumeratorPFPtBarrel()->Fill(leadjpt);
2032 if (
isEndCap(leadjeta)) v->getMEhisto_NumeratorPFPtEndcap()->Fill(leadjpt);
2033 if (
isForward(leadjeta)) v->getMEhisto_NumeratorPFPtForward()->Fill(leadjpt);
2034 v->getMEhisto_NumeratorPFEta()->Fill(leadjeta);
2035 v->getMEhisto_NumeratorPFPhi()->Fill(leadjphi);
2036 v->getMEhisto_NumeratorPFEtaPhi()->Fill(leadjeta,leadjphi);
2038 v->getMEhisto_NumeratorPFEtaBarrel()->Fill(leadjeta);
2039 v->getMEhisto_NumeratorPFPhiBarrel()->Fill(leadjphi);
2042 v->getMEhisto_NumeratorPFEtaEndcap()->Fill(leadjeta);
2043 v->getMEhisto_NumeratorPFPhiEndcap()->Fill(leadjphi);
2046 v->getMEhisto_NumeratorPFEtaForward()->Fill(leadjeta);
2047 v->getMEhisto_NumeratorPFPhiForward()->Fill(leadjphi);
2049 if((leadjpt > trigLowpTPFcut && !
isForward(leadjeta))
2050 || (leadjpt > trigLowpTPFcutFwd &&
isForward(leadjeta))){
2051 v->getMEhisto_NumeratorPFEta_LowpTcut()->Fill(leadjeta);
2052 v->getMEhisto_NumeratorPFPhi_LowpTcut()->Fill(leadjphi);
2053 v->getMEhisto_NumeratorPFEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
2055 if((leadjpt > trigMedpTPFcut && !
isForward(leadjeta))
2056 || (leadjpt > trigMedpTPFcutFwd &&
isForward(leadjeta))){
2057 v->getMEhisto_NumeratorPFEta_MedpTcut()->Fill(leadjeta);
2058 v->getMEhisto_NumeratorPFPhi_MedpTcut()->Fill(leadjphi);
2059 v->getMEhisto_NumeratorPFEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
2061 if((leadjpt > trigHighpTPFcut && !
isForward(leadjeta))
2062 || (leadjpt > trigHighpTPFcutFwd &&
isForward(leadjeta))){
2063 v->getMEhisto_NumeratorPFEta_HighpTcut()->Fill(leadjeta);
2064 v->getMEhisto_NumeratorPFPhi_HighpTcut()->Fill(leadjphi);
2065 v->getMEhisto_NumeratorPFEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
2071 if(jetIDbool ==
true && (v->getTriggerType().compare(
"DiJet_Trigger") == 0) &&
pfjet.size()>1){
2080 v->getMEhisto_DenominatorPFPt()->Fill((
PFJetPt[0] +
PFJetPt[1])/2.);
2083 v->getMEhisto_NumeratorPFPt()->Fill((
PFJetPt[0] +
PFJetPt[1])/2.);
2091 if(
verbose_)
std::cout <<
"fillMEforEffWrtMuTrigger: CaloMET -------------------" << std::endl;
2094 && (v->getPath().find(
"HLT_PFMET")==std::string::npos)){
2098 v->getMEhisto_DenominatorPt()->Fill(met.
et());
2099 v->getMEhisto_DenominatorPhi()->Fill(met.
phi());
2101 v->getMEhisto_NumeratorPt()->Fill(met.
et());
2102 v->getMEhisto_NumeratorPhi()->Fill(met.
phi());
2104 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
2107 double dPhimin = 9999.;
2109 for(trigger::Keys::const_iterator kj = khlt.begin();kj != khlt.end(); ++kj){
2115 v->getMEhisto_DeltaPhi()->Fill(dPhimin);
2121 if(
verbose_)
std::cout <<
"fillMEforEffWrtMuTrigger: PFMET -------------------" << std::endl;
2124 && (v->getPath().find(
"HLT_PFMET")!=std::string::npos)){
2127 const PFMET met = pfmetcol->front();
2128 v->getMEhisto_DenominatorPt()->Fill(met.
et());
2129 v->getMEhisto_DenominatorPhi()->Fill(met.
phi());
2131 v->getMEhisto_NumeratorPt()->Fill(met.
et());
2132 v->getMEhisto_NumeratorPhi()->Fill(met.
phi());
2134 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
2137 double dPhimin = 9999.;
2139 for(trigger::Keys::const_iterator kj = khlt.begin();kj != khlt.end(); ++kj){
2145 v->getMEhisto_DeltaPhi()->Fill(dPhimin);
2169 bool mbTrig =
false;
2170 bool numpassed =
false;
2174 for (
unsigned int j=0;
j!=nPath; ++
j) {
2179 if(
verbose_)
cout<<
"fillMEforEffWrtMBTrigger: MinBias Match"<<endl;
2204 double trigLowpTcut = 0;
2205 double trigMedpTcut = 0;
2206 double trigHighpTcut = 0;
2207 double trigLowpTcutFwd = 0;
2208 double trigMedpTcutFwd = 0;
2209 double trigHighpTcutFwd = 0;
2212 double trigLowpTPFcut = 0;
2213 double trigMedpTPFcut = 0;
2214 double trigHighpTPFcut = 0;
2215 double trigLowpTPFcutFwd = 0;
2216 double trigMedpTPFcutFwd = 0;
2217 double trigHighpTPFcutFwd = 0;
2220 size_t jetstrfound = v->getPath().find(
"Jet");
2222 string tpath = v->getPath();
2226 if(jetstrfound != string::npos){
2228 for(
int trig =
int(jetstrfound)+3; trig < int(jetstrfound)+7; trig++){
2229 if(!isdigit(tpath[trig]))
break;
2230 jetTrigVal+=tpath[trig];
2232 char *cjetTrigVal = (
char*)jetTrigVal.c_str();
2233 jetVal=atof(cjetTrigVal);
2239 trigMedpTcut = jetVal + 5.;
2240 trigHighpTcut = jetVal + 10.;
2242 trigLowpTcutFwd = jetVal + 9.;
2243 trigMedpTcutFwd = jetVal + 15.;
2244 trigHighpTcutFwd = jetVal + 21.;
2248 trigMedpTcut = jetVal + 2.;
2249 trigHighpTcut = jetVal + 60.;
2251 trigLowpTcutFwd = jetVal + 22.;
2252 trigMedpTcutFwd = jetVal + 25.;
2253 trigHighpTcutFwd = jetVal + 110.;
2255 trigLowpTcut = jetVal;
2262 trigMedpTPFcut = jetVal + 20.;
2263 trigHighpTPFcut = jetVal + 40.;
2265 trigLowpTPFcutFwd = jetVal + 60.;
2266 trigMedpTPFcutFwd = jetVal + 80.;
2267 trigHighpTPFcutFwd = jetVal + 100.;
2271 trigMedpTPFcut = jetVal + 40.;
2272 trigHighpTPFcut = jetVal + 140.;
2274 trigLowpTPFcutFwd = jetVal + 110.;
2275 trigMedpTPFcutFwd = jetVal + 130.;
2276 trigHighpTPFcutFwd = jetVal + 190.;
2278 trigLowpTPFcut = jetVal;
2285 if(
verbose_)
std::cout <<
"fillMEforEffWrtMuTrigger: CaloJet -------------------" << std::endl;
2291 bool jetIDbool =
false;
2298 if((v->getTriggerType().compare(
"SingleJet_Trigger") == 0) &&
calojet.size()){
2300 if(
verbose_)
cout<<
"passed CaloJet ID -------------------" << endl;
2304 v->getMEhisto_DenominatorPt()->Fill(leadjpt);
2305 if (
isBarrel(leadjeta)) v->getMEhisto_DenominatorPtBarrel()->Fill(leadjpt);
2306 if (
isEndCap(leadjeta)) v->getMEhisto_DenominatorPtEndcap()->Fill(leadjpt);
2307 if (
isForward(leadjeta)) v->getMEhisto_DenominatorPtForward()->Fill(leadjpt);
2308 v->getMEhisto_DenominatorEta()->Fill(leadjeta);
2309 v->getMEhisto_DenominatorPhi()->Fill(leadjphi);
2310 v->getMEhisto_DenominatorEtaPhi()->Fill(leadjeta,leadjphi);
2312 v->getMEhisto_DenominatorEtaBarrel()->Fill(leadjeta);
2313 v->getMEhisto_DenominatorPhiBarrel()->Fill(leadjphi);
2316 v->getMEhisto_DenominatorEtaEndcap()->Fill(leadjeta);
2317 v->getMEhisto_DenominatorPhiEndcap()->Fill(leadjphi);
2320 v->getMEhisto_DenominatorEtaForward()->Fill(leadjeta);
2321 v->getMEhisto_DenominatorPhiForward()->Fill(leadjphi);
2323 if((leadjpt > trigLowpTcut && !
isForward(leadjeta)) || (leadjpt > trigLowpTcutFwd &&
isForward(leadjeta))){
2324 v->getMEhisto_DenominatorEta_LowpTcut()->Fill(leadjeta);
2325 v->getMEhisto_DenominatorPhi_LowpTcut()->Fill(leadjphi);
2326 v->getMEhisto_DenominatorEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
2328 if((leadjpt > trigMedpTcut && !
isForward(leadjeta)) || (leadjpt > trigMedpTcutFwd &&
isForward(leadjeta))){
2329 v->getMEhisto_DenominatorEta_MedpTcut()->Fill(leadjeta);
2330 v->getMEhisto_DenominatorPhi_MedpTcut()->Fill(leadjphi);
2331 v->getMEhisto_DenominatorEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
2333 if((leadjpt > trigHighpTcut && !
isForward(leadjeta)) || (leadjpt > trigHighpTcutFwd &&
isForward(leadjeta))){
2334 v->getMEhisto_DenominatorEta_HighpTcut()->Fill(leadjeta);
2335 v->getMEhisto_DenominatorPhi_HighpTcut()->Fill(leadjphi);
2336 v->getMEhisto_DenominatorEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
2342 double dRmin = 99999.;
2343 double dPhimin = 9999.;
2344 if(v->getPath().find(
"L1") != std::string::npos){
2346 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
2350 for(trigger::Keys::const_iterator ki = kl1.begin();ki != kl1.end(); ++ki){
2351 double dR =
deltaR(toc[*ki].
eta(), toc[*ki].
phi(), leadjeta, leadjphi);
2360 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
2364 trigger::Keys::const_iterator kj = khlt.begin();
2365 for(;kj != khlt.end(); ++kj){
2367 leadjeta, leadjphi);
2377 v->getMEhisto_DeltaPhi()->Fill(dPhimin);
2378 v->getMEhisto_DeltaR()->Fill(dRmin);
2381 if(dRmin < 0.1 || (v->getPath().find(
"L1") != std::string::npos && dRmin < 0.4)){
2382 v->getMEhisto_NumeratorPt()->Fill(leadjpt);
2383 if (
isBarrel(leadjeta)) v->getMEhisto_NumeratorPtBarrel()->Fill(leadjpt);
2384 if (
isEndCap(leadjeta)) v->getMEhisto_NumeratorPtEndcap()->Fill(leadjpt);
2385 if (
isForward(leadjeta)) v->getMEhisto_NumeratorPtForward()->Fill(leadjpt);
2386 v->getMEhisto_NumeratorEta()->Fill(leadjeta);
2387 v->getMEhisto_NumeratorPhi()->Fill(leadjphi);
2388 v->getMEhisto_NumeratorEtaPhi()->Fill(leadjeta,leadjphi);
2390 v->getMEhisto_NumeratorEtaBarrel()->Fill(leadjeta);
2391 v->getMEhisto_NumeratorPhiBarrel()->Fill(leadjphi);
2394 v->getMEhisto_NumeratorEtaEndcap()->Fill(leadjeta);
2395 v->getMEhisto_NumeratorPhiEndcap()->Fill(leadjphi);
2398 v->getMEhisto_NumeratorEtaForward()->Fill(leadjeta);
2399 v->getMEhisto_NumeratorPhiForward()->Fill(leadjphi);
2401 if((leadjpt > trigLowpTcut && !
isForward(leadjeta)) || (leadjpt > trigLowpTcutFwd &&
isForward(leadjeta))){
2402 v->getMEhisto_NumeratorEta_LowpTcut()->Fill(leadjeta);
2403 v->getMEhisto_NumeratorPhi_LowpTcut()->Fill(leadjphi);
2404 v->getMEhisto_NumeratorEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
2406 if((leadjpt > trigMedpTcut && !
isForward(leadjeta)) || (leadjpt > trigMedpTcutFwd &&
isForward(leadjeta))){
2407 v->getMEhisto_NumeratorEta_MedpTcut()->Fill(leadjeta);
2408 v->getMEhisto_NumeratorPhi_MedpTcut()->Fill(leadjphi);
2409 v->getMEhisto_NumeratorEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
2411 if((leadjpt > trigHighpTcut && !
isForward(leadjeta)) || (leadjpt > trigHighpTcutFwd &&
isForward(leadjeta))){
2412 v->getMEhisto_NumeratorEta_HighpTcut()->Fill(leadjeta);
2413 v->getMEhisto_NumeratorPhi_HighpTcut()->Fill(leadjphi);
2414 v->getMEhisto_NumeratorEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
2421 if(jetIDbool ==
true && (v->getTriggerType().compare(
"DiJet_Trigger") == 0) &&
calojet.size()>1){
2426 if(numpassed==
true){
2435 if(
verbose_)
std::cout <<
"fillMEforEffWrtMuTrigger: PFJet -------------------" << std::endl;
2441 bool jetIDbool =
false;
2460 if((v->getTriggerType().compare(
"SingleJet_Trigger") == 0) &&
pfjet.size()){
2463 _pfMHT =
sqrt(pfMHTx*pfMHTx + pfMHTy*pfMHTy);
2464 v->getMEhisto_DenominatorPFMHT()->Fill(
_pfMHT);
2470 if(
verbose_)
cout<<
"passed PFJet ID -------------------" << endl;
2472 v->getMEhisto_DenominatorPFPt()->Fill(leadjpt);
2473 if (
isBarrel(leadjeta)) v->getMEhisto_DenominatorPFPtBarrel()->Fill(leadjpt);
2474 if (
isEndCap(leadjeta)) v->getMEhisto_DenominatorPFPtEndcap()->Fill(leadjpt);
2475 if (
isForward(leadjeta)) v->getMEhisto_DenominatorPFPtForward()->Fill(leadjpt);
2476 v->getMEhisto_DenominatorPFEta()->Fill(leadjeta);
2477 v->getMEhisto_DenominatorPFPhi()->Fill(leadjphi);
2478 v->getMEhisto_DenominatorPFEtaPhi()->Fill(leadjeta,leadjphi);
2480 v->getMEhisto_DenominatorPFEtaBarrel()->Fill(leadjeta);
2481 v->getMEhisto_DenominatorPFPhiBarrel()->Fill(leadjphi);
2484 v->getMEhisto_DenominatorPFEtaEndcap()->Fill(leadjeta);
2485 v->getMEhisto_DenominatorPFPhiEndcap()->Fill(leadjphi);
2488 v->getMEhisto_DenominatorPFEtaForward()->Fill(leadjeta);
2489 v->getMEhisto_DenominatorPFPhiForward()->Fill(leadjphi);
2491 if((leadjpt > trigLowpTPFcut && !
isForward(leadjeta)) || (leadjpt > trigLowpTPFcutFwd &&
isForward(leadjeta))){
2492 v->getMEhisto_DenominatorPFEta_LowpTcut()->Fill(leadjeta);
2493 v->getMEhisto_DenominatorPFPhi_LowpTcut()->Fill(leadjphi);
2494 v->getMEhisto_DenominatorPFEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
2496 if((leadjpt > trigMedpTPFcut && !
isForward(leadjeta)) || (leadjpt > trigMedpTPFcutFwd &&
isForward(leadjeta))){
2497 v->getMEhisto_DenominatorPFEta_MedpTcut()->Fill(leadjeta);
2498 v->getMEhisto_DenominatorPFPhi_MedpTcut()->Fill(leadjphi);
2499 v->getMEhisto_DenominatorPFEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
2501 if((leadjpt > trigHighpTPFcut && !
isForward(leadjeta)) || (leadjpt > trigHighpTPFcutFwd &&
isForward(leadjeta))){
2502 v->getMEhisto_DenominatorPFEta_HighpTcut()->Fill(leadjeta);
2503 v->getMEhisto_DenominatorPFPhi_HighpTcut()->Fill(leadjphi);
2504 v->getMEhisto_DenominatorPFEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
2509 double dRmin = 99999.;
2510 double dPhimin = 9999.;
2511 if(v->getPath().find(
"L1") != std::string::npos){
2513 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
2517 for(trigger::Keys::const_iterator ki = kl1.begin();ki != kl1.end(); ++ki){
2518 double dR =
deltaR(toc[*ki].
eta(), toc[*ki].
phi(), leadjeta, leadjphi);
2527 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
2531 for(trigger::Keys::const_iterator kj = khlt.begin();kj != khlt.end(); ++kj){
2532 double dR =
deltaR(toc[*kj].
eta(), toc[*kj].
phi(), leadjeta, leadjphi);
2541 v->getMEhisto_PFDeltaPhi()->Fill(dPhimin);
2542 v->getMEhisto_PFDeltaR()->Fill(dRmin);
2545 if(dRmin < 0.1 || (v->getPath().find(
"L1") != std::string::npos && dRmin < 0.4)){
2546 v->getMEhisto_NumeratorPFPt()->Fill(leadjpt);
2547 if (
isBarrel(leadjeta)) v->getMEhisto_NumeratorPFPtBarrel()->Fill(leadjpt);
2548 if (
isEndCap(leadjeta)) v->getMEhisto_NumeratorPFPtEndcap()->Fill(leadjpt);
2549 if (
isForward(leadjeta)) v->getMEhisto_NumeratorPFPtForward()->Fill(leadjpt);
2550 v->getMEhisto_NumeratorPFEta()->Fill(leadjeta);
2551 v->getMEhisto_NumeratorPFPhi()->Fill(leadjphi);
2552 v->getMEhisto_NumeratorPFEtaPhi()->Fill(leadjeta,leadjphi);
2554 v->getMEhisto_NumeratorPFEtaBarrel()->Fill(leadjeta);
2555 v->getMEhisto_NumeratorPFPhiBarrel()->Fill(leadjphi);
2558 v->getMEhisto_NumeratorPFEtaEndcap()->Fill(leadjeta);
2559 v->getMEhisto_NumeratorPFPhiEndcap()->Fill(leadjphi);
2562 v->getMEhisto_NumeratorPFEtaForward()->Fill(leadjeta);
2563 v->getMEhisto_NumeratorPFPhiForward()->Fill(leadjphi);
2565 if((leadjpt > trigLowpTPFcut && !
isForward(leadjeta))
2566 || (leadjpt > trigLowpTPFcutFwd &&
isForward(leadjeta))){
2567 v->getMEhisto_NumeratorPFEta_LowpTcut()->Fill(leadjeta);
2568 v->getMEhisto_NumeratorPFPhi_LowpTcut()->Fill(leadjphi);
2569 v->getMEhisto_NumeratorPFEtaPhi_LowpTcut()->Fill(leadjeta,leadjphi);
2571 if((leadjpt > trigMedpTPFcut && !
isForward(leadjeta))
2572 || (leadjpt > trigMedpTPFcutFwd &&
isForward(leadjeta))){
2573 v->getMEhisto_NumeratorPFEta_MedpTcut()->Fill(leadjeta);
2574 v->getMEhisto_NumeratorPFPhi_MedpTcut()->Fill(leadjphi);
2575 v->getMEhisto_NumeratorPFEtaPhi_MedpTcut()->Fill(leadjeta,leadjphi);
2577 if((leadjpt > trigHighpTPFcut && !
isForward(leadjeta))
2578 || (leadjpt > trigHighpTPFcutFwd &&
isForward(leadjeta))){
2579 v->getMEhisto_NumeratorPFEta_HighpTcut()->Fill(leadjeta);
2580 v->getMEhisto_NumeratorPFPhi_HighpTcut()->Fill(leadjphi);
2581 v->getMEhisto_NumeratorPFEtaPhi_HighpTcut()->Fill(leadjeta,leadjphi);
2587 if(jetIDbool ==
true && (v->getTriggerType().compare(
"DiJet_Trigger") == 0) &&
pfjet.size()>1){
2596 v->getMEhisto_DenominatorPFPt()->Fill((
PFJetPt[0] +
PFJetPt[1])/2.);
2599 v->getMEhisto_NumeratorPFPt()->Fill((
PFJetPt[0] +
PFJetPt[1])/2.);
2607 if(
verbose_)
std::cout <<
"fillMEforEffWrtMuTrigger: CaloMET -------------------" << std::endl;
2610 && (v->getPath().find(
"HLT_PFMET")==std::string::npos)){
2614 v->getMEhisto_DenominatorPt()->Fill(met.
et());
2615 v->getMEhisto_DenominatorPhi()->Fill(met.
phi());
2617 v->getMEhisto_NumeratorPt()->Fill(met.
et());
2618 v->getMEhisto_NumeratorPhi()->Fill(met.
phi());
2620 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
2623 double dPhimin = 9999.;
2625 for(trigger::Keys::const_iterator kj = khlt.begin();kj != khlt.end(); ++kj){
2631 v->getMEhisto_DeltaPhi()->Fill(dPhimin);
2637 if(
verbose_)
std::cout <<
"fillMEforEffWrtMuTrigger: PFMET -------------------" << std::endl;
2640 && (v->getPath().find(
"HLT_PFMET")!=std::string::npos)){
2643 const PFMET met = pfmetcol->front();
2644 v->getMEhisto_DenominatorPt()->Fill(met.
et());
2645 v->getMEhisto_DenominatorPhi()->Fill(met.
phi());
2647 v->getMEhisto_NumeratorPt()->Fill(met.
et());
2648 v->getMEhisto_NumeratorPhi()->Fill(met.
phi());
2650 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index hlt"<< hltIndex <<
" of that name ";
2653 double dPhimin = 9999.;
2655 for(trigger::Keys::const_iterator kj = khlt.begin();kj != khlt.end(); ++kj){
2661 v->getMEhisto_DeltaPhi()->Fill(dPhimin);
2704 LogDebug(
"HLTJetMETDQMSource") <<
"HLTConfigProvider failed to initialize.";
2720 for (
unsigned int i=0;
i!=
n; ++
i) {
2721 bool denomFound =
false;
2722 bool numFound =
false;
2723 bool mbFound =
false;
2724 bool muFound =
false;
2730 std::vector<std::string>::const_iterator controlPathname =
pathFilter_.begin();
2731 for(;controlPathname!=
pathFilter_.end(); ++controlPathname){
2732 if(pathname.find((*controlPathname)) != std::string::npos){
2737 if(checkPath==
false)
continue;
2740 std::vector<std::string>::const_iterator rejectPathname =
pathRejectKeyword_.begin();
2742 if(pathname.find((*rejectPathname)) != std::string::npos){
2747 if(checkPath==
false)
continue;
2755 unsigned int usedPrescale = 1;
2756 unsigned int objectType = 0;
2761 if(pathname.find(
"Jet") != std::string::npos
2762 && !(pathname.find(
"DoubleJet") != std::string::npos)
2763 && !(pathname.find(
"DiJet") != std::string::npos)
2764 && !(pathname.find(
"DiPFJet") != std::string::npos)
2765 && !(pathname.find(
"BTag") != std::string::npos)
2766 && !(pathname.find(
"Mu") != std::string::npos)
2767 && !(pathname.find(
"Fwd") != std::string::npos)){
2768 triggerType =
"SingleJet_Trigger";
2771 if(pathname.find(
"DiJet") != std::string::npos
2772 || pathname.find(
"DiPFJet") != std::string::npos
2773 || pathname.find(
"DoubleJet") != std::string::npos){
2774 triggerType =
"DiJet_Trigger";
2777 if(pathname.find(
"MET") != std::string::npos){
2778 triggerType =
"MET_Trigger";
2781 if(pathname.find(
"HT") != std::string::npos) {
2782 triggerType =
"TET_Trigger";
2788 && !(pathname.find(
"DiJet") != std::string::npos)
2789 && !(pathname.find(
"DiPFJet") != std::string::npos)
2790 && !(pathname.find(
"DoubleJet") != std::string::npos)){
2796 && ((pathname.find(
"DiJet") != std::string::npos)
2797 || (pathname.find(
"DiPFJet") != std::string::npos))){
2818 for (
unsigned int i=0;
i!=
n; ++
i){
2820 if(HLTname == pathname) numFound =
true;
2821 if(HLTname == dpathname) denomFound =
true;
2822 if(HLTname.find(
MBTrigPaths_[0]) != std::string::npos) mbFound =
true;
2823 if(HLTname.find(
MuonTrigPaths_[0]) != std::string::npos) muFound =
true;
2833 for(std::vector<std::string>::iterator numpathmodule = numpathmodules.begin(); numpathmodule!= numpathmodules.end(); ++numpathmodule ) {
2842 filtername = *numpathmodule;
2858 if(objectType != 0 && denomFound){
2860 for(std::vector<std::string>::iterator numpathmodule = numpathmodules.begin(); numpathmodule!= numpathmodules.end(); ++numpathmodule ) {
2869 Denomfiltername = *numpathmodule;
2873 if(objectType != 0 && numFound){
2875 cout<<
"=Pathname= "<<pathname
2876 <<
" | =Denompath= "<<dpathname
2877 <<
" | =Filtername= "<<filtername
2878 <<
" | =Denomfiltername= "<<Denomfiltername
2879 <<
" | =L1pathname= "<<l1pathname
2880 <<
" | =ObjectType= "<<objectType<<endl;
2881 if(!((pathname.find(
"HT") != std::string::npos) || (pathname.find(
"Quad") != std::string::npos))){
2884 filtername, Denomfiltername,
processname_, objectType, triggerType));
2886 filtername, Denomfiltername,
processname_, objectType, triggerType));
2888 filtername, Denomfiltername,
processname_, objectType, triggerType));
2890 filtername, Denomfiltername,
processname_, objectType, triggerType));
2893 filtername, Denomfiltername,
processname_, objectType, triggerType));
2897 if (
verbose_)
cout <<
"get names for efficicncy------------------"<< endl;
2902 unsigned int usedPrescale = 1;
2903 unsigned int objectType = 0;
2907 for (std::vector<std::pair<std::string, std::string> >::iterator custompathnamepair =
custompathnamepairs_.begin();
2910 std::string dpathname = custompathnamepair->second;
2911 bool numFound =
false;
2912 bool denomFound =
false;
2914 for (
unsigned int i=0;
i!=
n; ++
i) {
2916 if(HLTname.find(pathname)!=std::string::npos){numFound =
true; pathname = HLTname;}
2917 if(HLTname.find(dpathname)!=std::string::npos){denomFound =
true; dpathname = HLTname;}
2919 if(numFound && denomFound){
2920 if (pathname.find(
"Jet") != std::string::npos
2921 && !(pathname.find(
"DiJet") != std::string::npos)
2922 && !(pathname.find(
"DiPFJet") != std::string::npos)
2923 && !(pathname.find(
"DoubleJet") != std::string::npos)
2924 && !(pathname.find(
"BTag") != std::string::npos)
2925 && !(pathname.find(
"Mu") != std::string::npos)
2926 && !(pathname.find(
"Fwd") != std::string::npos)){
2927 triggerType =
"SingleJet_Trigger";
2930 if (pathname.find(
"DiJet") != std::string::npos
2931 || pathname.find(
"DiPFJet") != std::string::npos
2932 || pathname.find(
"DoubleJet") != std::string::npos ){
2933 triggerType =
"DiJet_Trigger";
2936 if (pathname.find(
"MET") != std::string::npos ){
2937 triggerType =
"MET_Trigger";
2940 if (pathname.find(
"TET") != std::string::npos ){
2941 triggerType =
"TET_Trigger";
2947 for(std::vector<std::string>::iterator numpathmodule = numpathmodules.begin(); numpathmodule!= numpathmodules.end(); ++numpathmodule ) {
2956 filtername = *numpathmodule;
2959 if(objectType != 0){
2961 for(std::vector<std::string>::iterator numpathmodule = numpathmodules.begin(); numpathmodule!= numpathmodules.end(); ++numpathmodule ) {
2970 Denomfiltername = *numpathmodule;
2974 <<
"==denompath=="<<dpathname
2975 <<
"==filtername=="<<filtername
2976 <<
"==denomfiltername=="<<Denomfiltername
2977 <<
"==l1pathname=="<<l1pathname
2978 <<
"==objectType=="<<objectType<<endl;
2980 filtername, Denomfiltername,
processname_, objectType, triggerType));
2986 if(
verbose_)
cout<<
"== end hltPathsEff_.push_back ======" << endl;
2996 double TrigMin_ = -0.5;
3001 rate_All = dbe->
book1D(histonm.c_str(),histot.c_str(),TrigBins_,TrigMin_,TrigMax_);
3003 histonm =
"JetMET_TriggerRate_Correlation";
3004 histot =
"JetMET TriggerRate Correlation Summary;y&&!x;x&&y";
3005 correlation_All = dbe->
book2D(histonm.c_str(),histot.c_str(),TrigBins_,TrigMin_,TrigMax_,TrigBins_,TrigMin_,TrigMax_);
3007 histonm =
"JetMET_NVertices";
3008 histot =
"No. of vertices";
3011 histonm =
"JetMET_PVZ";
3012 histot =
"Primary Vertex Z pos";
3013 PVZ = dbe->
book1D(histonm.c_str(),histot.c_str(),100,-50.,50.);
3016 histonm=
"JetMET_TriggerRate_WrtMuTrigger";
3017 histot=
"JetMET TriggerRate Summary Wrt Muon Trigger ";
3020 histonm=
"JetMET_TriggerRate_Correlation_WrtMuTrigger";
3021 histot=
"JetMET TriggerRate Correlation Summary Wrt Muon Trigger;y&&!x;x&&y";
3024 histonm=
"JetMET_TriggerRate_WrtMBTrigger";
3025 histot=
"JetMET TriggerRate Summary Wrt MB Trigger";
3028 histonm=
"JetMET_TriggerRate_Correlation_WrtMBTrigger";
3029 histot=
"JetMET TriggerRate Correlation Wrt MB Trigger;y&&!x;x&&y";
3038 labelnm =
v->getPath();
3040 for(
int ibin=1; ibin<nbins+1; ibin++){
3041 const char * binLabel =
rate_All->
getTH1()->GetXaxis()->GetBinLabel(ibin);
3043 if(binLabel_str.compare(labelnm)==0)
break;
3044 if(binLabel[0]==
'\0'){
3071 double PtMax_ = 1000.;
3075 double EtaMin_ = -5.;
3076 double EtaMax_ = 5.;
3079 double PhiMin_ = -3.5;
3080 double PhiMax_ = 3.5;
3083 double ResMin_ = -1.5;
3084 double ResMax_ = 1.5;
3101 if(
v->getObjectType() ==
trigger::TriggerJet &&
v->getTriggerType().compare(
"SingleJet_Trigger") == 0){
3103 histoname = labelname+
"_recObjN";
3104 title = labelname+
"_recObjN;Reco multiplicity()"+trigPath;
3108 histoname = labelname+
"_recObjPt";
3109 title = labelname+
"_recObjPt; Reco Pt[GeV/c]"+trigPath;
3113 histoname = labelname+
"_recObjPtBarrel";
3114 title = labelname+
"_recObjPtBarrel;Reco Pt[GeV/c]"+trigPath;
3115 MonitorElement * PtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3118 histoname = labelname+
"_recObjPtEndcap";
3119 title = labelname+
"_recObjPtEndcap;Reco Pt[GeV/c]"+trigPath;
3120 MonitorElement * PtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3123 histoname = labelname+
"_recObjPtForward";
3124 title = labelname+
"_recObjPtForward;Reco Pt[GeV/c]"+trigPath;
3125 MonitorElement * PtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3128 histoname = labelname+
"_recObjEta";
3129 title = labelname+
"_recObjEta;Reco #eta"+trigPath;
3133 histoname = labelname+
"_recObjPhi";
3134 title = labelname+
"_recObjPhi;Reco #Phi"+trigPath;
3138 histoname = labelname+
"_recObjEtaPhi";
3139 title = labelname+
"_recObjEtaPhi;Reco #eta;Reco #Phi"+trigPath;
3140 MonitorElement * EtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
3143 histoname = labelname+
"_l1ObjPt";
3144 title = labelname+
"_l1ObjPt;L1 Pt[GeV/c]"+trigPath;
3148 histoname = labelname+
"_l1ObjEta";
3149 title = labelname+
"_l1ObjEta;L1 #eta"+trigPath;
3150 MonitorElement * Eta_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
3153 histoname = labelname+
"_l1ObjPhi";
3154 title = labelname+
"_l1ObjPhi;L1 #Phi"+trigPath;
3155 MonitorElement * Phi_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3158 histoname = labelname+
"_l1ObjEtaPhi";
3159 title = labelname+
"_l1ObjEtaPhi;L1 #eta;L1 #Phi"+trigPath;
3160 MonitorElement * EtaPhi_L1 = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
3163 histoname = labelname+
"_l1ObjN";
3164 title = labelname+
"_l1ObjN;L1 multiplicity"+trigPath;
3168 histoname = labelname+
"_l1ObjPtBarrel";
3169 title = labelname+
"_l1ObjPtBarrel;L1 Pt[GeV/c]"+trigPath;
3170 MonitorElement * PtBarrel_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3173 histoname = labelname+
"_l1ObjPtEndcap";
3174 title = labelname+
"_l1ObjPtEndcap;L1 Pt[GeV/c]"+trigPath;
3175 MonitorElement * PtEndcap_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3178 histoname = labelname+
"_l1ObjPtForward";
3179 title = labelname+
"_l1ObjPtForward;L1 Pt[GeV/c]"+trigPath;
3180 MonitorElement * PtForward_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3183 histoname = labelname+
"_hltObjN";
3184 title = labelname+
"_hltObjN;HLT multiplicity"+trigPath;
3188 histoname = labelname+
"_hltObjPtBarrel";
3189 title = labelname+
"_hltObjPtBarrel;HLT Pt[GeV/c]"+trigPath;
3190 MonitorElement * PtBarrel_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3193 histoname = labelname+
"_hltObjPtEndcap";
3194 title = labelname+
"_hltObjPtEndcap;HLT Pt[GeV/c]"+trigPath;
3195 MonitorElement * PtEndcap_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3198 histoname = labelname+
"_hltObjPtForward";
3199 title = labelname+
"_hltObjPtForward;HLT Pt[GeV/c]"+trigPath;
3200 MonitorElement * PtForward_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3203 histoname = labelname+
"_hltObjPt";
3204 title = labelname+
"_hltObjPt;HLT Pt[GeV/c]"+trigPath;
3208 histoname = labelname+
"_hltObjEta";
3209 title = labelname+
"_hltObjEta;HLT #eta"+trigPath;
3210 MonitorElement * Eta_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
3213 histoname = labelname+
"_hltObjPhi";
3214 title = labelname+
"_hltObjPhi;HLT #Phi"+trigPath;
3215 MonitorElement * Phi_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3218 histoname = labelname+
"_hltObjEtaPhi";
3219 title = labelname+
"_hltObjEtaPhi;HLT #eta;HLT #Phi"+trigPath;
3220 MonitorElement * EtaPhi_HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
3223 histoname = labelname+
"_l1HLTPtResolution";
3224 title = labelname+
"_l1HLTPtResolution;(Pt(L1)-Pt(HLT))/Pt(L1)"+trigPath;
3225 MonitorElement * PtResolution_L1HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3226 PtResolution_L1HLT->
getTH1();
3228 histoname = labelname+
"_l1HLTEtaResolution";
3229 title = labelname+
"_l1HLTEtaResolution;(#eta(L1)-#eta(HLT))/#eta(L1)"+trigPath;
3230 MonitorElement * EtaResolution_L1HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3231 EtaResolution_L1HLT->
getTH1();
3233 histoname = labelname+
"_l1HLTPhiResolution";
3234 title = labelname+
"_l1HLTPhiResolution;(#Phi(L1)-#Phi(HLT))/#Phi(L1)"+trigPath;
3235 MonitorElement * PhiResolution_L1HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3236 PhiResolution_L1HLT->
getTH1();
3238 histoname = labelname+
"_l1HLTPtCorrelation";
3239 title = labelname+
"_l1HLTPtCorrelation;Pt(L1)[GeV/c];Pt(HLT)[GeV/c]"+trigPath;
3240 MonitorElement * PtCorrelation_L1HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_,Ptbins_,PtMin_,PtMax_);
3241 PtCorrelation_L1HLT->
getTH1();
3243 histoname = labelname+
"_l1HLTEtaCorrelation";
3244 title = labelname+
"_l1HLTEtaCorrelation;#eta(L1);#eta(HLT)"+trigPath;
3245 MonitorElement * EtaCorrelation_L1HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Etabins_,EtaMin_,EtaMax_);
3246 EtaCorrelation_L1HLT->
getTH1();
3248 histoname = labelname+
"_l1HLTPhiCorrelation";
3249 title = labelname+
"_l1HLTPhiCorrelation;#Phi(L1);#Phi(HLT)"+trigPath;
3250 MonitorElement * PhiCorrelation_L1HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_,Phibins_,PhiMin_,PhiMax_);
3251 PhiCorrelation_L1HLT->
getTH1();
3253 histoname = labelname+
"_hltRecObjPtResolution";
3254 title = labelname+
"_hltRecObjPtResolution;(Pt(HLT)-Pt(Reco))/Pt(HLT)"+trigPath;
3255 MonitorElement * PtResolution_HLTRecObj = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3256 PtResolution_HLTRecObj->
getTH1();
3258 histoname = labelname+
"_hltRecObjEtaResolution";
3259 title = labelname+
"_hltRecObjEtaResolution;(#eta(HLT)-#eta(Reco))/#eta(HLT)"+trigPath;
3260 MonitorElement * EtaResolution_HLTRecObj = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3261 EtaResolution_HLTRecObj->
getTH1();
3263 histoname = labelname+
"_hltRecObjPhiResolution";
3264 title = labelname+
"_hltRecObjPhiResolution;(#Phi(HLT)-#Phi(Reco))/#Phi(HLT)"+trigPath;
3265 MonitorElement * PhiResolution_HLTRecObj = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3266 PhiResolution_HLTRecObj->
getTH1();
3268 histoname = labelname+
"_hltRecObjPtCorrelation";
3269 title = labelname+
"_hltRecObjPtCorrelation;Pt(HLT)[GeV/c];Pt(Reco)[GeV/c]"+trigPath;
3270 MonitorElement * PtCorrelation_HLTRecObj = dbe->
book2D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_,Ptbins_,PtMin_,PtMax_);
3271 PtCorrelation_HLTRecObj->
getTH1();
3273 histoname = labelname+
"_hltRecObjEtaCorrelation";
3274 title = labelname+
"_hltRecObjEtaCorrelation;#eta(HLT);#eta(Reco)"+trigPath;
3275 MonitorElement * EtaCorrelation_HLTRecObj = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Etabins_,EtaMin_,EtaMax_);
3276 EtaCorrelation_HLTRecObj->
getTH1();
3278 histoname = labelname+
"_hltRecObjPhiCorrelation";
3279 title = labelname+
"_hltRecObjPhiCorrelation;#Phi(HLT);#Phi(Reco)"+trigPath;
3280 MonitorElement * PhiCorrelation_HLTRecObj = dbe->
book2D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_,Phibins_,PhiMin_,PhiMax_);
3281 PhiCorrelation_HLTRecObj->
getTH1();
3308 EtaResolution_L1HLT,
3309 PhiResolution_L1HLT,
3310 PtResolution_HLTRecObj,
3311 EtaResolution_HLTRecObj,
3312 PhiResolution_HLTRecObj,
3313 PtCorrelation_L1HLT,
3314 EtaCorrelation_L1HLT,
3315 PhiCorrelation_L1HLT,
3316 PtCorrelation_HLTRecObj,
3317 EtaCorrelation_HLTRecObj,
3318 PhiCorrelation_HLTRecObj,
3333 histoname = labelname+
"_RecObjAveragePt";
3334 title = labelname+
"_RecObjAveragePt;Reco Average Pt[GeV/c]"+trigPath;
3335 MonitorElement * jetAveragePt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3338 histoname = labelname+
"_RecObjAverageEta";
3339 title = labelname+
"_RecObjAverageEta;Reco Average #eta"+trigPath;
3340 MonitorElement * jetAverageEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
3343 histoname = labelname+
"_RecObjPhiDifference";
3344 title = labelname+
"_RecObjPhiDifference;Reco #Delta#Phi"+trigPath;
3345 MonitorElement * jetPhiDifference = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3346 jetPhiDifference->
getTH1();
3348 histoname = labelname+
"_hltObjAveragePt";
3349 title = labelname+
"_hltObjAveragePt;HLT Average Pt[GeV/c]"+trigPath;
3350 MonitorElement * hltAveragePt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3353 histoname = labelname+
"_hltObjAverageEta";
3354 title = labelname+
"_hltObjAverageEta;HLT Average #eta"+trigPath;
3355 MonitorElement * hltAverageEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
3358 histoname = labelname+
"_hltObjPhiDifference";
3359 title = labelname+
"_hltObjPhiDifference;Reco #Delta#Phi"+trigPath;
3360 MonitorElement * hltPhiDifference = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3361 hltPhiDifference->
getTH1();
3413 histoname = labelname+
"_recObjPt";
3414 title = labelname+
"_recObjPt;Reco Pt[GeV/c]"+trigPath;
3418 histoname = labelname+
"_recObjPhi";
3419 title = labelname+
"_recObjPhi;Reco #Phi"+trigPath;
3423 histoname = labelname+
"_l1ObjPt";
3424 title = labelname+
"_l1ObjPt;L1 Pt[GeV/c]"+trigPath;
3428 histoname = labelname+
"_l1ObjPhi";
3429 title = labelname+
"_l1ObjPhi;L1 #Phi"+trigPath;
3430 MonitorElement * Phi_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3433 histoname = labelname+
"_hltObjPt";
3434 title = labelname+
"_hltObjPt;HLT Pt[GeV/c]"+trigPath;
3438 histoname = labelname+
"_hltObjPhi";
3439 title = labelname+
"_hltObjPhi;HLT #Phi"+trigPath;
3440 MonitorElement * Phi_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3443 histoname = labelname+
"_l1HLTPtResolution";
3444 title = labelname+
"_l1HLTPtResolution;(Pt(L1)-Pt(HLT))/Pt(L1)"+trigPath;
3445 MonitorElement * PtResolution_L1HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3446 PtResolution_L1HLT->
getTH1();
3448 histoname = labelname+
"_l1HLTPhiResolution";
3449 title = labelname+
"_l1HLTPhiResolution;(#Phi(L1)-#Phi(HLT))/#Phi(L1)"+trigPath;
3450 MonitorElement * PhiResolution_L1HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3451 PhiResolution_L1HLT->
getTH1();
3453 histoname = labelname+
"_l1HLTPtCorrelation";
3454 title = labelname+
"_l1HLTPtCorrelation;Pt(L1)[GeV/c];Pt(HLT)[GeV/c]"+trigPath;
3455 MonitorElement * PtCorrelation_L1HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_,Ptbins_,PtMin_,PtMax_);
3456 PtCorrelation_L1HLT->
getTH1();
3458 histoname = labelname+
"_l1HLTPhiCorrelation";
3459 title = labelname+
"_l1HLTPhiCorrelation;#Phi(L1);#Phi(HLT)"+trigPath;
3460 MonitorElement * PhiCorrelation_L1HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_,Phibins_,PhiMin_,PhiMax_);
3461 PhiCorrelation_L1HLT->
getTH1();
3463 histoname = labelname+
"_hltRecObjPtResolution";
3464 title = labelname+
"_hltRecObjPtResolution;(Pt(HLT)-Pt(Reco))/Pt(HLT)"+trigPath;
3465 MonitorElement * PtResolution_HLTRecObj = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3466 PtResolution_HLTRecObj->
getTH1();
3468 histoname = labelname+
"_hltRecObjPhiResolution";
3469 title = labelname+
"_hltRecObjPhiResolution;(#Phi(HLT)-#Phi(Reco))/#Phi(HLT)"+trigPath;
3470 MonitorElement * PhiResolution_HLTRecObj = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3471 PhiResolution_HLTRecObj->
getTH1();
3473 histoname = labelname+
"_hltRecObjPtCorrelation";
3474 title = labelname+
"_hltRecObjPtCorrelation;Pt(HLT)[GeV/c];Pt(Reco)[GeV/c]"+trigPath;
3475 MonitorElement * PtCorrelation_HLTRecObj = dbe->
book2D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_,Ptbins_,PtMin_,PtMax_);
3476 PtCorrelation_HLTRecObj->
getTH1();
3478 histoname = labelname+
"_hltRecObjPhiCorrelation";
3479 title = labelname+
"_hltRecObjPhiCorrelation;#Phi(HLT);#Phi(Reco)"+trigPath;
3480 MonitorElement * PhiCorrelation_HLTRecObj = dbe->
book2D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_,Phibins_,PhiMin_,PhiMax_);
3481 PhiCorrelation_HLTRecObj->
getTH1();
3509 PhiResolution_L1HLT,
3510 PtResolution_HLTRecObj,
3512 PhiResolution_HLTRecObj,
3513 PtCorrelation_L1HLT,
3515 PhiCorrelation_L1HLT,
3516 PtCorrelation_HLTRecObj,
3518 PhiCorrelation_HLTRecObj,
3541 double PtMax_ = 1000.;
3544 double EtaMin_ = -5.;
3545 double EtaMax_ = 5.;
3548 double PhiMin_ = -3.5;
3549 double PhiMax_ = 3.5;
3552 double ResMin_ = -1.5;
3553 double ResMax_ = 1.5;
3571 if(
v->getObjectType() ==
trigger::TriggerJet &&
v->getTriggerType().compare(
"SingleJet_Trigger") == 0){
3573 histoname = labelname+
"_recObjN";
3574 title = labelname+
"_recObjN;Reco multiplicity()"+trigPath;
3578 histoname = labelname+
"_recObjPt";
3579 title = labelname+
"_recObjPt; Reco Pt[GeV/c]"+trigPath;
3583 histoname = labelname+
"_recObjPtBarrel";
3584 title = labelname+
"_recObjPtBarrel;Reco Pt[GeV/c]"+trigPath;
3585 MonitorElement * PtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3588 histoname = labelname+
"_recObjPtEndcap";
3589 title = labelname+
"_recObjPtEndcap;Reco Pt[GeV/c]"+trigPath;
3590 MonitorElement * PtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3593 histoname = labelname+
"_recObjPtForward";
3594 title = labelname+
"_recObjPtForward;Reco Pt[GeV/c]"+trigPath;
3595 MonitorElement * PtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3598 histoname = labelname+
"_recObjEta";
3599 title = labelname+
"_recObjEta;Reco #eta"+trigPath;
3603 histoname = labelname+
"_recObjPhi";
3604 title = labelname+
"_recObjPhi;Reco #Phi"+trigPath;
3608 histoname = labelname+
"_recObjEtaPhi";
3609 title = labelname+
"_recObjEtaPhi;Reco #eta;Reco #Phi"+trigPath;
3610 MonitorElement * EtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
3613 histoname = labelname+
"_l1ObjN";
3614 title = labelname+
"_l1ObjN;L1 multiplicity"+trigPath;
3618 histoname = labelname+
"_l1ObjPt";
3619 title = labelname+
"_l1ObjPt;L1 Pt[GeV/c]"+trigPath;
3623 histoname = labelname+
"_l1ObjPtBarrel";
3624 title = labelname+
"_l1ObjPtBarrel;L1 Pt[GeV/c]"+trigPath;
3625 MonitorElement * PtBarrel_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3628 histoname = labelname+
"_l1ObjPtEndcap";
3629 title = labelname+
"_l1ObjPtEndcap;L1 Pt[GeV/c]"+trigPath;
3630 MonitorElement * PtEndcap_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3633 histoname = labelname+
"_l1ObjPtForward";
3634 title = labelname+
"_l1ObjPtForward;L1 Pt[GeV/c]"+trigPath;
3635 MonitorElement * PtForward_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3638 histoname = labelname+
"_l1ObjEta";
3639 title = labelname+
"_l1ObjEta;L1 #eta"+trigPath;
3640 MonitorElement * Eta_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
3643 histoname = labelname+
"_l1ObjPhi";
3644 title = labelname+
"_l1ObjPhi;L1 #Phi"+trigPath;
3645 MonitorElement * Phi_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3648 histoname = labelname+
"_l1ObjEtaPhi";
3649 title = labelname+
"_l1ObjEtaPhi;L1 #eta;L1 #Phi"+trigPath;
3650 MonitorElement * EtaPhi_L1 = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
3653 histoname = labelname+
"_hltObjN";
3654 title = labelname+
"_hltObjN;HLT multiplicity"+trigPath;
3658 histoname = labelname+
"_hltObjPt";
3659 title = labelname+
"_hltObjPt;HLT Pt[GeV/c]"+trigPath;
3663 histoname = labelname+
"_hltObjPtBarrel";
3664 title = labelname+
"_hltObjPtBarrel;HLT Pt[GeV/c]"+trigPath;
3665 MonitorElement * PtBarrel_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3668 histoname = labelname+
"_hltObjPtEndcap";
3669 title = labelname+
"_hltObjPtEndcap;HLT Pt[GeV/c]"+trigPath;
3670 MonitorElement * PtEndcap_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3673 histoname = labelname+
"_hltObjPtForward";
3674 title = labelname+
"_hltObjPtForward;HLT Pt[GeV/c]"+trigPath;
3675 MonitorElement * PtForward_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3678 histoname = labelname+
"_hltObjEta";
3679 title = labelname+
"_hltObjEta;HLT #eta"+trigPath;
3680 MonitorElement * Eta_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
3683 histoname = labelname+
"_hltObjPhi";
3684 title = labelname+
"_hltObjPhi;HLT #Phi"+trigPath;
3685 MonitorElement * Phi_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3688 histoname = labelname+
"_hltObjEtaPhi";
3689 title = labelname+
"_hltObjEtaPhi;HLT #eta;HLT #Phi"+trigPath;
3690 MonitorElement * EtaPhi_HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
3693 histoname = labelname+
"_l1HLTPtResolution";
3694 title = labelname+
"_l1HLTPtResolution;(Pt(L1)-Pt(HLT))/Pt(L1)"+trigPath;
3695 MonitorElement * PtResolution_L1HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3696 PtResolution_L1HLT->
getTH1();
3698 histoname = labelname+
"_l1HLTEtaResolution";
3699 title = labelname+
"_l1HLTEtaResolution;(#eta(L1)-#eta(HLT))/#eta(L1)"+trigPath;
3700 MonitorElement * EtaResolution_L1HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3701 EtaResolution_L1HLT->
getTH1();
3703 histoname = labelname+
"_l1HLTPhiResolution";
3704 title = labelname+
"_l1HLTPhiResolution;(#Phi(L1)-#Phi(HLT))/#Phi(L1)"+trigPath;
3705 MonitorElement * PhiResolution_L1HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3706 PhiResolution_L1HLT->
getTH1();
3708 histoname = labelname+
"_l1HLTPtCorrelation";
3709 title = labelname+
"_l1HLTPtCorrelation;Pt(L1)[GeV/c];Pt(HLT)[GeV/c]"+trigPath;
3710 MonitorElement * PtCorrelation_L1HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_,Ptbins_,PtMin_,PtMax_);
3711 PtCorrelation_L1HLT->
getTH1();
3713 histoname = labelname+
"_l1HLTEtaCorrelation";
3714 title = labelname+
"_l1HLTEtaCorrelation;#eta(L1);#eta(HLT)"+trigPath;
3715 MonitorElement * EtaCorrelation_L1HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Etabins_,EtaMin_,EtaMax_);
3716 EtaCorrelation_L1HLT->
getTH1();
3718 histoname = labelname+
"_l1HLTPhiCorrelation";
3719 title = labelname+
"_l1HLTPhiCorrelation;#Phi(L1);#Phi(HLT)"+trigPath;
3720 MonitorElement * PhiCorrelation_L1HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_,Phibins_,PhiMin_,PhiMax_);
3721 PhiCorrelation_L1HLT->
getTH1();
3723 histoname = labelname+
"_hltRecObjPtResolution";
3724 title = labelname+
"_hltRecObjPtResolution;(Pt(HLT)-Pt(Reco))/Pt(HLT)"+trigPath;
3725 MonitorElement * PtResolution_HLTRecObj = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3726 PtResolution_HLTRecObj->
getTH1();
3728 histoname = labelname+
"_hltRecObjEtaResolution";
3729 title = labelname+
"_hltRecObjEtaResolution;(#eta(HLT)-#eta(Reco))/#eta(HLT)"+trigPath;
3730 MonitorElement * EtaResolution_HLTRecObj = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3731 EtaResolution_HLTRecObj->
getTH1();
3733 histoname = labelname+
"_hltRecObjPhiResolution";
3734 title = labelname+
"_hltRecObjPhiResolution;(#Phi(HLT)-#Phi(Reco))/#Phi(HLT)"+trigPath;
3735 MonitorElement * PhiResolution_HLTRecObj = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3736 PhiResolution_HLTRecObj->
getTH1();
3738 histoname = labelname+
"_hltRecObjPtCorrelation";
3739 title = labelname+
"_hltRecObjPtCorrelation;Pt(HLT)[GeV/c];Pt(Reco)[GeV/c]"+trigPath;
3740 MonitorElement * PtCorrelation_HLTRecObj = dbe->
book2D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_,Ptbins_,PtMin_,PtMax_);
3741 PtCorrelation_HLTRecObj->
getTH1();
3743 histoname = labelname+
"_hltRecObjEtaCorrelation";
3744 title = labelname+
"_hltRecObjEtaCorrelation;#eta(HLT);#eta(Reco)"+trigPath;
3745 MonitorElement * EtaCorrelation_HLTRecObj = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Etabins_,EtaMin_,EtaMax_);
3746 EtaCorrelation_HLTRecObj->
getTH1();
3748 histoname = labelname+
"_hltRecObjPhiCorrelation";
3749 title = labelname+
"_hltRecObjPhiCorrelation;#Phi(HLT);#Phi(Reco)"+trigPath;
3750 MonitorElement * PhiCorrelation_HLTRecObj = dbe->
book2D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_,Phibins_,PhiMin_,PhiMax_);
3751 PhiCorrelation_HLTRecObj->
getTH1();
3778 EtaResolution_L1HLT,
3779 PhiResolution_L1HLT,
3780 PtResolution_HLTRecObj,
3781 EtaResolution_HLTRecObj,
3782 PhiResolution_HLTRecObj,
3783 PtCorrelation_L1HLT,
3784 EtaCorrelation_L1HLT,
3785 PhiCorrelation_L1HLT,
3786 PtCorrelation_HLTRecObj,
3787 EtaCorrelation_HLTRecObj,
3788 PhiCorrelation_HLTRecObj,
3803 histoname = labelname+
"_RecObjAveragePt";
3804 title = labelname+
"_RecObjAveragePt;Reco Average Pt[GeV/c]"+trigPath;
3805 MonitorElement * jetAveragePt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3808 histoname = labelname+
"_RecObjAverageEta";
3809 title = labelname+
"_RecObjAverageEta;Reco Average #eta"+trigPath;
3810 MonitorElement * jetAverageEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
3813 histoname = labelname+
"_RecObjPhiDifference";
3814 title = labelname+
"_RecObjPhiDifference;Reco #Delta#Phi"+trigPath;
3815 MonitorElement * jetPhiDifference = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3816 jetPhiDifference->
getTH1();
3818 histoname = labelname+
"_hltObjAveragePt";
3819 title = labelname+
"_hltObjAveragePt;HLT Average Pt[GeV/c]"+trigPath;
3820 MonitorElement * hltAveragePt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
3823 histoname = labelname+
"_hltObjAverageEta";
3824 title = labelname+
"_hltObjAverageEta;HLT Average #eta"+trigPath;
3825 MonitorElement * hltAverageEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
3828 histoname = labelname+
"_hltObjPhiDifference";
3829 title = labelname+
"_hltObjPhiDifference;Reco #Delta#Phi"+trigPath;
3830 MonitorElement * hltPhiDifference = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3831 hltPhiDifference->
getTH1();
3883 histoname = labelname+
"_recObjPt";
3884 title = labelname+
"_recObjPt;Reco Pt[GeV/c]"+trigPath;
3888 histoname = labelname+
"_recObjPhi";
3889 title = labelname+
"_recObjPhi;Reco #Phi"+trigPath;
3893 histoname = labelname+
"_l1ObjPt";
3894 title = labelname+
"_l1ObjPt;L1 Pt[GeV/c]"+trigPath;
3898 histoname = labelname+
"_l1ObjPhi";
3899 title = labelname+
"_l1ObjPhi;L1 #Phi"+trigPath;
3900 MonitorElement * Phi_L1 = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3903 histoname = labelname+
"_hltObjPt";
3904 title = labelname+
"_hltObjPt;HLT Pt[GeV/c]"+trigPath;
3908 histoname = labelname+
"_hltObjPhi";
3909 title = labelname+
"_hltObjPhi;HLT #Phi"+trigPath;
3910 MonitorElement * Phi_HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
3913 histoname = labelname+
"_l1HLTPtResolution";
3914 title = labelname+
"_l1HLTPtResolution;(Pt(L1)-Pt(HLT))/Pt(L1)"+trigPath;
3915 MonitorElement * PtResolution_L1HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3916 PtResolution_L1HLT->
getTH1();
3918 histoname = labelname+
"_l1HLTPhiResolution";
3919 title = labelname+
"_l1HLTPhiResolution;(#Phi(L1)-#Phi(HLT))/#Phi(L1)"+trigPath;
3920 MonitorElement * PhiResolution_L1HLT = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3921 PhiResolution_L1HLT->
getTH1();
3923 histoname = labelname+
"_l1HLTPtCorrelation";
3924 title = labelname+
"_l1HLTPtCorrelation;Pt(L1)[GeV/c];Pt(HLT)[GeV/c]"+trigPath;
3925 MonitorElement * PtCorrelation_L1HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_,Ptbins_,PtMin_,PtMax_);
3926 PtCorrelation_L1HLT->
getTH1();
3928 histoname = labelname+
"_l1HLTPhiCorrelation";
3929 title = labelname+
"_l1HLTPhiCorrelation;#Phi(L1);#Phi(HLT)"+trigPath;
3930 MonitorElement * PhiCorrelation_L1HLT = dbe->
book2D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_,Phibins_,PhiMin_,PhiMax_);
3931 PhiCorrelation_L1HLT->
getTH1();
3933 histoname = labelname+
"_hltRecObjPtResolution";
3934 title = labelname+
"_hltRecObjPtResolution;(Pt(HLT)-Pt(Reco))/Pt(HLT)"+trigPath;
3935 MonitorElement * PtResolution_HLTRecObj = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3936 PtResolution_HLTRecObj->
getTH1();
3938 histoname = labelname+
"_hltRecObjPhiResolution";
3939 title = labelname+
"_hltRecObjPhiResolution;(#Phi(HLT)-#Phi(Reco))/#Phi(HLT)"+trigPath;
3940 MonitorElement * PhiResolution_HLTRecObj = dbe->
book1D(histoname.c_str(),title.c_str(),Resbins_,ResMin_,ResMax_);
3941 PhiResolution_HLTRecObj->
getTH1();
3943 histoname = labelname+
"_hltRecObjPtCorrelation";
3944 title = labelname+
"_hltRecObjPtCorrelation;Pt(HLT)[GeV/c];Pt(Reco)[GeV/c]"+trigPath;
3945 MonitorElement * PtCorrelation_HLTRecObj = dbe->
book2D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_,Ptbins_,PtMin_,PtMax_);
3946 PtCorrelation_HLTRecObj->
getTH1();
3948 histoname = labelname+
"_hltRecObjPhiCorrelation";
3949 title = labelname+
"_hltRecObjPhiCorrelation;#Phi(HLT);#Phi(Reco)"+trigPath;
3950 MonitorElement * PhiCorrelation_HLTRecObj = dbe->
book2D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_,Phibins_,PhiMin_,PhiMax_);
3951 PhiCorrelation_HLTRecObj->
getTH1();
3979 PhiResolution_L1HLT,
3980 PtResolution_HLTRecObj,
3982 PhiResolution_HLTRecObj,
3983 PtCorrelation_L1HLT,
3985 PhiCorrelation_L1HLT,
3986 PtCorrelation_HLTRecObj,
3988 PhiCorrelation_HLTRecObj,
4008 double PtMax_ = 1000.;
4011 double EtaMin_ = -5.;
4012 double EtaMax_ = 5.;
4015 double PhiMin_ = -3.5;
4016 double PhiMax_ = 3.5;
4025 std::string subdirName = dirName1 + trgPathName +
"_wrt_" + trgPathNameD;
4033 if((
v->getObjectType() ==
trigger::TriggerJet) && (
v->getTriggerType().compare(
"SingleJet_Trigger") == 0)){
4035 histoname = labelname+
"_NumeratorPt";
4036 title = labelname+
"NumeratorPt;Calo Pt[GeV/c]";
4037 MonitorElement * NumeratorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4040 histoname = labelname+
"_NumeratorPtBarrel";
4041 title = labelname+
"NumeratorPtBarrel;Calo Pt[GeV/c] ";
4042 MonitorElement * NumeratorPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4043 NumeratorPtBarrel->
getTH1();
4045 histoname = labelname+
"_NumeratorPtEndcap";
4046 title = labelname+
"NumeratorPtEndcap;Calo Pt[GeV/c]";
4047 MonitorElement * NumeratorPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4048 NumeratorPtEndcap->
getTH1();
4050 histoname = labelname+
"_NumeratorPtForward";
4051 title = labelname+
"NumeratorPtForward;Calo Pt[GeV/c]";
4052 MonitorElement * NumeratorPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4053 NumeratorPtForward->
getTH1();
4055 histoname = labelname+
"_NumeratorEta";
4056 title = labelname+
"NumeratorEta;Calo #eta ";
4057 MonitorElement * NumeratorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4060 histoname = labelname+
"_NumeratorPhi";
4061 title = labelname+
"NumeratorPhi;Calo #Phi";
4062 MonitorElement * NumeratorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4065 histoname = labelname+
"_NumeratorEtaPhi";
4066 title = labelname+
"NumeratorEtaPhi;Calo #eta;Calo #Phi";
4067 MonitorElement * NumeratorEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4068 NumeratorEtaPhi->
getTH1();
4070 histoname = labelname+
"_NumeratorEtaBarrel";
4071 title = labelname+
"NumeratorEtaBarrel;Calo #eta ";
4072 MonitorElement * NumeratorEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4073 NumeratorEtaBarrel->
getTH1();
4075 histoname = labelname+
"_NumeratorPhiBarrel";
4076 title = labelname+
"NumeratorPhiBarrel;Calo #Phi";
4077 MonitorElement * NumeratorPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4078 NumeratorPhiBarrel->
getTH1();
4080 histoname = labelname+
"_NumeratorEtaEndcap";
4081 title = labelname+
"NumeratorEtaEndcap;Calo #eta ";
4082 MonitorElement * NumeratorEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4083 NumeratorEtaEndcap->
getTH1();
4085 histoname = labelname+
"_NumeratorPhiEndcap";
4086 title = labelname+
"NumeratorPhiEndcap;Calo #Phi";
4087 MonitorElement * NumeratorPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4088 NumeratorPhiEndcap->
getTH1();
4090 histoname = labelname+
"_NumeratorEtaForward";
4091 title = labelname+
"NumeratorEtaForward;Calo #eta ";
4092 MonitorElement * NumeratorEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4093 NumeratorEtaForward->
getTH1();
4095 histoname = labelname+
"_NumeratorPhiForward";
4096 title = labelname+
"NumeratorPhiForward;Calo #Phi";
4097 MonitorElement * NumeratorPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4098 NumeratorPhiForward->
getTH1();
4100 histoname = labelname+
"_NumeratorEta_LowpTcut";
4101 title = labelname+
"NumeratorEta_LowpTcut;Calo #eta ";
4102 MonitorElement * NumeratorEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4103 NumeratorEta_LowpTcut->
getTH1();
4105 histoname = labelname+
"_NumeratorPhi_LowpTcut";
4106 title = labelname+
"NumeratorPhi_LowpTcut;Calo #Phi";
4107 MonitorElement * NumeratorPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4108 NumeratorPhi_LowpTcut->
getTH1();
4110 histoname = labelname+
"_NumeratorEtaPhi_LowpTcut";
4111 title = labelname+
"NumeratorEtaPhi_LowpTcut;Calo #eta;Calo #Phi";
4112 MonitorElement * NumeratorEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4113 NumeratorEtaPhi_LowpTcut->
getTH1();
4115 histoname = labelname+
"_NumeratorEta_MedpTcut";
4116 title = labelname+
"NumeratorEta_MedpTcut;Calo #eta ";
4117 MonitorElement * NumeratorEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4118 NumeratorEta_MedpTcut->
getTH1();
4120 histoname = labelname+
"_NumeratorPhi_MedpTcut";
4121 title = labelname+
"NumeratorPhi_MedpTcut;Calo #Phi";
4122 MonitorElement * NumeratorPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4123 NumeratorPhi_MedpTcut->
getTH1();
4125 histoname = labelname+
"_NumeratorEtaPhi_MedpTcut";
4126 title = labelname+
"NumeratorEtaPhi_MedpTcut;Calo #eta;Calo #Phi";
4127 MonitorElement * NumeratorEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4128 NumeratorEtaPhi_MedpTcut->
getTH1();
4130 histoname = labelname+
"_NumeratorEta_HighpTcut";
4131 title = labelname+
"NumeratorEta_HighpTcut;Calo #eta ";
4132 MonitorElement * NumeratorEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4133 NumeratorEta_HighpTcut->
getTH1();
4135 histoname = labelname+
"_NumeratorPhi_HighpTcut";
4136 title = labelname+
"NumeratorPhi_HighpTcut;Calo #Phi";
4137 MonitorElement * NumeratorPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4138 NumeratorPhi_HighpTcut->
getTH1();
4140 histoname = labelname+
"_NumeratorEtaPhi_HighpTcut";
4141 title = labelname+
"NumeratorEtaPhi_HighpTcut;Calo #eta;Calo #Phi";
4142 MonitorElement * NumeratorEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4143 NumeratorEtaPhi_HighpTcut->
getTH1();
4145 histoname = labelname+
"_DenominatorPt";
4146 title = labelname+
"DenominatorPt;Calo Pt[GeV/c]";
4147 MonitorElement * DenominatorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4150 histoname = labelname+
"_DenominatorPtBarrel";
4151 title = labelname+
"DenominatorPtBarrel;Calo Pt[GeV/c]";
4152 MonitorElement * DenominatorPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4153 DenominatorPtBarrel->
getTH1();
4155 histoname = labelname+
"_DenominatorPtEndcap";
4156 title = labelname+
"DenominatorPtEndcap;Calo Pt[GeV/c]";
4157 MonitorElement * DenominatorPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4158 DenominatorPtEndcap->
getTH1();
4160 histoname = labelname+
"_DenominatorPtForward";
4161 title = labelname+
"DenominatorPtForward;Calo Pt[GeV/c] ";
4162 MonitorElement * DenominatorPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4163 DenominatorPtForward->
getTH1();
4165 histoname = labelname+
"_DenominatorEta";
4166 title = labelname+
"DenominatorEta;Calo #eta ";
4167 MonitorElement * DenominatorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4168 DenominatorEta->
getTH1();
4170 histoname = labelname+
"_DenominatorPhi";
4171 title = labelname+
"DenominatorPhi;Calo #Phi";
4172 MonitorElement * DenominatorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4173 DenominatorPhi->
getTH1();
4175 histoname = labelname+
"_DenominatorEtaPhi";
4176 title = labelname+
"DenominatorEtaPhi;Calo #eta; Calo #Phi";
4177 MonitorElement * DenominatorEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4178 DenominatorEtaPhi->
getTH1();
4180 histoname = labelname+
"_DenominatorEtaBarrel";
4181 title = labelname+
"DenominatorEtaBarrel;Calo #eta ";
4182 MonitorElement * DenominatorEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4183 DenominatorEtaBarrel->
getTH1();
4185 histoname = labelname+
"_DenominatorPhiBarrel";
4186 title = labelname+
"DenominatorPhiBarrel;Calo #Phi";
4187 MonitorElement * DenominatorPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4188 DenominatorPhiBarrel->
getTH1();
4190 histoname = labelname+
"_DenominatorEtaEndcap";
4191 title = labelname+
"DenominatorEtaEndcap;Calo #eta ";
4192 MonitorElement * DenominatorEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4193 DenominatorEtaEndcap->
getTH1();
4195 histoname = labelname+
"_DenominatorPhiEndcap";
4196 title = labelname+
"DenominatorPhiEndcap;Calo #Phi";
4197 MonitorElement * DenominatorPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4198 DenominatorPhiEndcap->
getTH1();
4200 histoname = labelname+
"_DenominatorEtaForward";
4201 title = labelname+
"DenominatorEtaForward;Calo #eta ";
4202 MonitorElement * DenominatorEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4203 DenominatorEtaForward->
getTH1();
4205 histoname = labelname+
"_DenominatorPhiForward";
4206 title = labelname+
"DenominatorPhiForward;Calo #Phi";
4207 MonitorElement * DenominatorPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4208 DenominatorPhiForward->
getTH1();
4210 histoname = labelname+
"_DenominatorEta_LowpTcut";
4211 title = labelname+
"DenominatorEta_LowpTcut;Calo #eta ";
4212 MonitorElement * DenominatorEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4213 DenominatorEta_LowpTcut->
getTH1();
4215 histoname = labelname+
"_DenominatorPhi_LowpTcut";
4216 title = labelname+
"DenominatorPhi_LowpTcut;Calo #Phi";
4217 MonitorElement * DenominatorPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4218 DenominatorPhi_LowpTcut->
getTH1();
4220 histoname = labelname+
"_DenominatorEtaPhi_LowpTcut";
4221 title = labelname+
"DenominatorEtaPhi_LowpTcut;Calo #eta;Calo #Phi";
4222 MonitorElement * DenominatorEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4223 DenominatorEtaPhi_LowpTcut->
getTH1();
4225 histoname = labelname+
"_DenominatorEta_MedpTcut";
4226 title = labelname+
"DenominatorEta_MedpTcut;Calo #eta ";
4227 MonitorElement * DenominatorEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4228 DenominatorEta_MedpTcut->
getTH1();
4230 histoname = labelname+
"_DenominatorPhi_MedpTcut";
4231 title = labelname+
"DenominatorPhi_MedpTcut;Calo #Phi";
4232 MonitorElement * DenominatorPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4233 DenominatorPhi_MedpTcut->
getTH1();
4235 histoname = labelname+
"_DenominatorEtaPhi_MedpTcut";
4236 title = labelname+
"DenominatorEtaPhi_MedpTcut;Calo #eta;Calo #Phi";
4237 MonitorElement * DenominatorEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4238 DenominatorEtaPhi_MedpTcut->
getTH1();
4240 histoname = labelname+
"_DenominatorEta_HighpTcut";
4241 title = labelname+
"DenominatorEta_HighpTcut;Calo #eta ";
4242 MonitorElement * DenominatorEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4243 DenominatorEta_HighpTcut->
getTH1();
4245 histoname = labelname+
"_DenominatorPhi_HighpTcut";
4246 title = labelname+
"DenominatorPhi_HighpTcut;Calo #Phi";
4247 MonitorElement * DenominatorPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4248 DenominatorPhi_HighpTcut->
getTH1();
4250 histoname = labelname+
"_DenominatorEtaPhi_HighpTcut";
4251 title = labelname+
"DenominatorEtaPhi_HighpTcut;Calo #eta;Calo #Phi";
4252 MonitorElement * DenominatorEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4253 DenominatorEtaPhi_HighpTcut->
getTH1();
4255 histoname = labelname+
"_DeltaR";
4256 title = labelname+
"DeltaR;";
4260 histoname = labelname+
"_DeltaPhi";
4261 title = labelname+
"DeltaPhi;";
4265 histoname = labelname+
"_NumeratorPFMHT";
4266 title = labelname+
"NumeratorPFMHT;PFMHT[GeV/c]";
4267 MonitorElement * NumeratorPFMHT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4268 NumeratorPFMHT->
getTH1();
4270 histoname = labelname+
"_NumeratorPFPt";
4271 title = labelname+
"NumeratorPFPt;PF Pt[GeV/c]";
4272 MonitorElement * NumeratorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4275 histoname = labelname+
"_NumeratorPFPtBarrel";
4276 title = labelname+
"NumeratorPFPtBarrel;PF Pt[GeV/c] ";
4277 MonitorElement * NumeratorPFPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4278 NumeratorPFPtBarrel->
getTH1();
4280 histoname = labelname+
"_NumeratorPFPtEndcap";
4281 title = labelname+
"NumeratorPFPtEndcap;PF Pt[GeV/c]";
4282 MonitorElement * NumeratorPFPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4283 NumeratorPFPtEndcap->
getTH1();
4285 histoname = labelname+
"_NumeratorPFPtForward";
4286 title = labelname+
"NumeratorPFPtForward;PF Pt[GeV/c]";
4287 MonitorElement * NumeratorPFPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4288 NumeratorPFPtForward->
getTH1();
4290 histoname = labelname+
"_NumeratorPFEta";
4291 title = labelname+
"NumeratorPFEta;PF #eta ";
4292 MonitorElement * NumeratorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4293 NumeratorPFEta->
getTH1();
4295 histoname = labelname+
"_NumeratorPFPhi";
4296 title = labelname+
"NumeratorPFPhi;Calo #Phi";
4297 MonitorElement * NumeratorPFPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4298 NumeratorPFPhi->
getTH1();
4300 histoname = labelname+
"_NumeratorPFEtaPhi";
4301 title = labelname+
"NumeratorPFEtaPhi;PF #eta;Calo #Phi";
4302 MonitorElement * NumeratorPFEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4303 NumeratorPFEtaPhi->
getTH1();
4305 histoname = labelname+
"_NumeratorPFEtaBarrel";
4306 title = labelname+
"NumeratorPFEtaBarrel;PF #eta ";
4307 MonitorElement * NumeratorPFEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4308 NumeratorPFEtaBarrel->
getTH1();
4310 histoname = labelname+
"_NumeratorPFPhiBarrel";
4311 title = labelname+
"NumeratorPFPhiBarrel;PF #Phi";
4312 MonitorElement * NumeratorPFPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4313 NumeratorPFPhiBarrel->
getTH1();
4315 histoname = labelname+
"_NumeratorPFEtaEndcap";
4316 title = labelname+
"NumeratorPFEtaEndcap;Calo #eta ";
4317 MonitorElement * NumeratorPFEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4318 NumeratorPFEtaEndcap->
getTH1();
4320 histoname = labelname+
"_NumeratorPFPhiEndcap";
4321 title = labelname+
"NumeratorPFPhiEndcap;PF #Phi";
4322 MonitorElement * NumeratorPFPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4323 NumeratorPFPhiEndcap->
getTH1();
4325 histoname = labelname+
"_NumeratorPFEtaForward";
4326 title = labelname+
"NumeratorPFEtaForward;Calo #eta ";
4327 MonitorElement * NumeratorPFEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4328 NumeratorPFEtaForward->
getTH1();
4330 histoname = labelname+
"_NumeratorPFPhiForward";
4331 title = labelname+
"NumeratorPFPhiForward;PF #Phi";
4332 MonitorElement * NumeratorPFPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4333 NumeratorPFPhiForward->
getTH1();
4335 histoname = labelname+
"_NumeratorPFEta_LowpTcut";
4336 title = labelname+
"NumeratorPFEta_LowpTcut;PF #eta ";
4337 MonitorElement * NumeratorPFEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4338 NumeratorPFEta_LowpTcut->
getTH1();
4340 histoname = labelname+
"_NumeratorPFPhi_LowpTcut";
4341 title = labelname+
"NumeratorPFPhi_LowpTcut;PF #Phi";
4342 MonitorElement * NumeratorPFPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4343 NumeratorPFPhi_LowpTcut->
getTH1();
4345 histoname = labelname+
"_NumeratorPFEtaPhi_LowpTcut";
4346 title = labelname+
"NumeratorPFEtaPhi_LowpTcut;PF #eta;Calo #Phi";
4347 MonitorElement * NumeratorPFEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4348 NumeratorPFEtaPhi_LowpTcut->
getTH1();
4350 histoname = labelname+
"_NumeratorPFEta_MedpTcut";
4351 title = labelname+
"NumeratorPFEta_MedpTcut;PF #eta ";
4352 MonitorElement * NumeratorPFEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4353 NumeratorPFEta_MedpTcut->
getTH1();
4355 histoname = labelname+
"_NumeratorPFPhi_MedpTcut";
4356 title = labelname+
"NumeratorPFPhi_MedpTcut;PF #Phi";
4357 MonitorElement * NumeratorPFPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4358 NumeratorPFPhi_MedpTcut->
getTH1();
4360 histoname = labelname+
"_NumeratorPFEtaPhi_MedpTcut";
4361 title = labelname+
"NumeratorPFEtaPhi_MedpTcut;PF #eta;PF #Phi";
4362 MonitorElement * NumeratorPFEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4363 NumeratorPFEtaPhi_MedpTcut->
getTH1();
4365 histoname = labelname+
"_NumeratorPFEta_HighpTcut";
4366 title = labelname+
"NumeratorPFEta_HighpTcut;Calo #eta ";
4367 MonitorElement * NumeratorPFEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4368 NumeratorPFEta_HighpTcut->
getTH1();
4370 histoname = labelname+
"_NumeratorPFPhi_HighpTcut";
4371 title = labelname+
"NumeratorPFPhi_HighpTcut;PF #Phi";
4372 MonitorElement * NumeratorPFPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4373 NumeratorPFPhi_HighpTcut->
getTH1();
4375 histoname = labelname+
"_NumeratorPFEtaPhi_HighpTcut";
4376 title = labelname+
"NumeratorPFEtaPhi_HighpTcut;PF #eta;PF #Phi";
4377 MonitorElement * NumeratorPFEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4378 NumeratorPFEtaPhi_HighpTcut->
getTH1();
4380 histoname = labelname+
"_DenominatorPFMHT";
4381 title = labelname+
"DenominatorPFMHT;PF Pt[GeV/c]";
4382 MonitorElement * DenominatorPFMHT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4383 DenominatorPFMHT->
getTH1();
4385 histoname = labelname+
"_DenominatorPFPt";
4386 title = labelname+
"DenominatorPFPt;PF Pt[GeV/c]";
4387 MonitorElement * DenominatorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4388 DenominatorPFPt->
getTH1();
4390 histoname = labelname+
"_DenominatorPFPtBarrel";
4391 title = labelname+
"DenominatorPFPtBarrel;Calo Pt[GeV/c]";
4392 MonitorElement * DenominatorPFPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4393 DenominatorPFPtBarrel->
getTH1();
4395 histoname = labelname+
"_DenominatorPFPtEndcap";
4396 title = labelname+
"DenominatorPFPtEndcap;PF Pt[GeV/c]";
4397 MonitorElement * DenominatorPFPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4398 DenominatorPFPtEndcap->
getTH1();
4400 histoname = labelname+
"_DenominatorPFPtForward";
4401 title = labelname+
"DenominatorPFPtForward;PF Pt[GeV/c] ";
4402 MonitorElement * DenominatorPFPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4403 DenominatorPFPtForward->
getTH1();
4405 histoname = labelname+
"_DenominatorPFEta";
4406 title = labelname+
"DenominatorPFEta;PF #eta ";
4407 MonitorElement * DenominatorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4408 DenominatorPFEta->
getTH1();
4410 histoname = labelname+
"_DenominatorPFPhi";
4411 title = labelname+
"DenominatorPFPhi;PF #Phi";
4412 MonitorElement * DenominatorPFPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4413 DenominatorPFPhi->
getTH1();
4415 histoname = labelname+
"_DenominatorPFEtaPhi";
4416 title = labelname+
"DenominatorPFEtaPhi;PF #eta; Calo #Phi";
4417 MonitorElement * DenominatorPFEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4418 DenominatorPFEtaPhi->
getTH1();
4420 histoname = labelname+
"_DenominatorPFEtaBarrel";
4421 title = labelname+
"DenominatorPFEtaBarrel;Calo #eta ";
4422 MonitorElement * DenominatorPFEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4423 DenominatorPFEtaBarrel->
getTH1();
4425 histoname = labelname+
"_DenominatorPFPhiBarrel";
4426 title = labelname+
"DenominatorPFPhiBarrel;PF #Phi";
4427 MonitorElement * DenominatorPFPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4428 DenominatorPFPhiBarrel->
getTH1();
4430 histoname = labelname+
"_DenominatorPFEtaEndcap";
4431 title = labelname+
"DenominatorPFEtaEndcap;PF #eta ";
4432 MonitorElement * DenominatorPFEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4433 DenominatorPFEtaEndcap->
getTH1();
4435 histoname = labelname+
"_DenominatorPFPhiEndcap";
4436 title = labelname+
"DenominatorPFPhiEndcap;Calo #Phi";
4437 MonitorElement * DenominatorPFPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4438 DenominatorPFPhiEndcap->
getTH1();
4440 histoname = labelname+
"_DenominatorPFEtaForward";
4441 title = labelname+
"DenominatorPFEtaForward;PF #eta ";
4442 MonitorElement * DenominatorPFEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4443 DenominatorPFEtaForward->
getTH1();
4445 histoname = labelname+
"_DenominatorPFPhiForward";
4446 title = labelname+
"DenominatorPFPhiForward;PF #Phi";
4447 MonitorElement * DenominatorPFPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4448 DenominatorPFPhiForward->
getTH1();
4450 histoname = labelname+
"_DenominatorPFEta_LowpTcut";
4451 title = labelname+
"DenominatorPFEta_LowpTcut;PF #eta ";
4452 MonitorElement * DenominatorPFEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4453 DenominatorPFEta_LowpTcut->
getTH1();
4455 histoname = labelname+
"_DenominatorPFPhi_LowpTcut";
4456 title = labelname+
"DenominatorPFPhi_LowpTcut;PF #Phi";
4457 MonitorElement * DenominatorPFPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4458 DenominatorPFPhi_LowpTcut->
getTH1();
4460 histoname = labelname+
"_DenominatorPFEtaPhi_LowpTcut";
4461 title = labelname+
"DenominatorPFEtaPhi_LowpTcut;PF #eta;Calo #Phi";
4462 MonitorElement * DenominatorPFEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4463 DenominatorPFEtaPhi_LowpTcut->
getTH1();
4465 histoname = labelname+
"_DenominatorPFEta_MedpTcut";
4466 title = labelname+
"DenominatorPFEta_MedpTcut;PF #eta ";
4467 MonitorElement * DenominatorPFEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4468 DenominatorPFEta_MedpTcut->
getTH1();
4470 histoname = labelname+
"_DenominatorPFPhi_MedpTcut";
4471 title = labelname+
"DenominatorPFPhi_MedpTcut;PF #Phi";
4472 MonitorElement * DenominatorPFPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4473 DenominatorPFPhi_MedpTcut->
getTH1();
4475 histoname = labelname+
"_DenominatorPFEtaPhi_MedpTcut";
4476 title = labelname+
"DenominatorPFEtaPhi_MedpTcut;PF #eta;Calo #Phi";
4477 MonitorElement * DenominatorPFEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4478 DenominatorPFEtaPhi_MedpTcut->
getTH1();
4480 histoname = labelname+
"_DenominatorPFEta_HighpTcut";
4481 title = labelname+
"DenominatorPFEta_HighpTcut;PF #eta ";
4482 MonitorElement * DenominatorPFEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4483 DenominatorPFEta_HighpTcut->
getTH1();
4485 histoname = labelname+
"_DenominatorPFPhi_HighpTcut";
4486 title = labelname+
"DenominatorPFPhi_HighpTcut;PF #Phi";
4487 MonitorElement * DenominatorPFPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4488 DenominatorPFPhi_HighpTcut->
getTH1();
4490 histoname = labelname+
"_DenominatorPFEtaPhi_HighpTcut";
4491 title = labelname+
"DenominatorPFEtaPhi_HighpTcut;PF #eta;Calo #Phi";
4492 MonitorElement * DenominatorPFEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4493 DenominatorPFEtaPhi_HighpTcut->
getTH1();
4495 histoname = labelname+
"_PFDeltaR";
4496 title = labelname+
"PFDeltaR;";
4500 histoname = labelname+
"_PFDeltaPhi";
4501 title = labelname+
"PFDeltaPhi;";
4505 v->setEffHistos(NumeratorPt,
4517 NumeratorEtaForward,
4518 NumeratorPhiForward,
4519 NumeratorEta_LowpTcut,
4520 NumeratorPhi_LowpTcut,
4521 NumeratorEtaPhi_LowpTcut,
4522 NumeratorEta_MedpTcut,
4523 NumeratorPhi_MedpTcut,
4524 NumeratorEtaPhi_MedpTcut,
4525 NumeratorEta_HighpTcut,
4526 NumeratorPhi_HighpTcut,
4527 NumeratorEtaPhi_HighpTcut,
4530 DenominatorPtBarrel,
4531 DenominatorPtEndcap,
4532 DenominatorPtForward,
4537 DenominatorEtaBarrel,
4538 DenominatorPhiBarrel,
4539 DenominatorEtaEndcap,
4540 DenominatorPhiEndcap,
4541 DenominatorEtaForward,
4542 DenominatorPhiForward,
4543 DenominatorEta_LowpTcut,
4544 DenominatorPhi_LowpTcut,
4545 DenominatorEtaPhi_LowpTcut,
4546 DenominatorEta_MedpTcut,
4547 DenominatorPhi_MedpTcut,
4548 DenominatorEtaPhi_MedpTcut,
4549 DenominatorEta_HighpTcut,
4550 DenominatorPhi_HighpTcut,
4551 DenominatorEtaPhi_HighpTcut,
4557 NumeratorPFPtBarrel,
4558 NumeratorPFPtEndcap,
4559 NumeratorPFPtForward,
4563 NumeratorPFEtaBarrel,
4564 NumeratorPFPhiBarrel,
4565 NumeratorPFEtaEndcap,
4566 NumeratorPFPhiEndcap,
4567 NumeratorPFEtaForward,
4568 NumeratorPFPhiForward,
4569 NumeratorPFEta_LowpTcut,
4570 NumeratorPFPhi_LowpTcut,
4571 NumeratorPFEtaPhi_LowpTcut,
4572 NumeratorPFEta_MedpTcut,
4573 NumeratorPFPhi_MedpTcut,
4574 NumeratorPFEtaPhi_MedpTcut,
4575 NumeratorPFEta_HighpTcut,
4576 NumeratorPFPhi_HighpTcut,
4577 NumeratorPFEtaPhi_HighpTcut,
4580 DenominatorPFPtBarrel,
4581 DenominatorPFPtEndcap,
4582 DenominatorPFPtForward,
4585 DenominatorPFEtaPhi,
4586 DenominatorPFEtaBarrel,
4587 DenominatorPFPhiBarrel,
4588 DenominatorPFEtaEndcap,
4589 DenominatorPFPhiEndcap,
4590 DenominatorPFEtaForward,
4591 DenominatorPFPhiForward,
4592 DenominatorPFEta_LowpTcut,
4593 DenominatorPFPhi_LowpTcut,
4594 DenominatorPFEtaPhi_LowpTcut,
4595 DenominatorPFEta_MedpTcut,
4596 DenominatorPFPhi_MedpTcut,
4597 DenominatorPFEtaPhi_MedpTcut,
4598 DenominatorPFEta_HighpTcut,
4599 DenominatorPFPhi_HighpTcut,
4600 DenominatorPFEtaPhi_HighpTcut,
4607 if((
v->getObjectType() ==
trigger::TriggerJet) && (
v->getTriggerType().compare(
"DiJet_Trigger") == 0)){
4609 histoname = labelname+
"_NumeratorAvrgPt";
4610 title = labelname+
"NumeratorAvrgPt;Calo Pt[GeV/c]";
4611 MonitorElement * NumeratorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4614 histoname = labelname+
"_NumeratorAvrgEta";
4615 title = labelname+
"NumeratorAvrgEta;Calo #eta";
4616 MonitorElement * NumeratorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4619 histoname = labelname+
"_DenominatorAvrgPt";
4620 title = labelname+
"DenominatorAvrgPt;Calo Pt[GeV/c] ";
4621 MonitorElement * DenominatorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4624 histoname = labelname+
"_DenominatorAvrgEta";
4625 title = labelname+
"DenominatorAvrgEta;Calo #eta";
4626 MonitorElement * DenominatorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4627 DenominatorEta->
getTH1();
4629 histoname = labelname+
"_DeltaR";
4630 title = labelname+
"DeltaR;";
4634 histoname = labelname+
"_DeltaPhi";
4635 title = labelname+
"DeltaPhi;";
4640 histoname = labelname+
"_NumeratorAvrgPFPt";
4641 title = labelname+
"NumeratorAvrgPFPt;PF Pt[GeV/c]";
4642 MonitorElement * NumeratorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4645 histoname = labelname+
"_NumeratorAvrgPFEta";
4646 title = labelname+
"NumeratorAvrgPFEta;PF #eta";
4647 MonitorElement * NumeratorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4648 NumeratorPFEta->
getTH1();
4650 histoname = labelname+
"_DenominatorAvrgPFPt";
4651 title = labelname+
"DenominatorAvrgPFPt;PF Pt[GeV/c] ";
4652 MonitorElement * DenominatorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4653 DenominatorPFPt->
getTH1();
4655 histoname = labelname+
"_DenominatorAvrgPFEta";
4656 title = labelname+
"DenominatorAvrgPFEta;PF #eta";
4657 MonitorElement * DenominatorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4658 DenominatorPFEta->
getTH1();
4660 histoname = labelname+
"_PFDeltaR";
4661 title = labelname+
"PFDeltaR;";
4665 histoname = labelname+
"_PFDeltaPhi";
4666 title = labelname+
"PFDeltaPhi;";
4670 v->setEffHistos( dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4671 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4672 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4673 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4674 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
4685 histoname = labelname+
"_NumeratorPt";
4686 if(
v->getPath().find(
"HLT_PFMET")==std::string::npos)
4687 title = labelname+
"NumeratorPt; CaloMET[GeV/c]";
4689 title = labelname+
"NumeratorPt; PFMET[GeV/c]";
4690 MonitorElement * NumeratorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4693 histoname = labelname+
"_NumeratorPhi";
4694 title = labelname+
"NumeratorPhi; #Phi";
4695 MonitorElement * NumeratorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4698 histoname = labelname+
"_DenominatorPt";
4699 if(
v->getPath().find(
"HLT_PFMET")==std::string::npos)
4700 title = labelname+
"DenominatorPt; CaloMET[GeV/c]";
4702 title = labelname+
"DenominatorPt; PFMET[GeV/c]";
4703 MonitorElement * DenominatorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4706 histoname = labelname+
"_DenominatorPhi";
4707 title = labelname+
"DenominatorPhi; #Phi";
4708 MonitorElement * DenominatorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4709 DenominatorPhi->
getTH1();
4711 v->setEffHistos( NumeratorPt, dummy, dummy, dummy, dummy, NumeratorPhi, dummy, dummy, dummy, dummy,
4712 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4713 dummy, dummy, DenominatorPt, dummy, dummy, dummy, dummy, DenominatorPhi, dummy, dummy,
4714 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4715 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4716 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
4717 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, 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
4741 if((
v->getObjectType() ==
trigger::TriggerJet) && (
v->getTriggerType().compare(
"SingleJet_Trigger") == 0)){
4743 histoname = labelname+
"_NumeratorPt";
4744 title = labelname+
"NumeratorPt;Pt[GeV/c]";
4745 MonitorElement * NumeratorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4748 histoname = labelname+
"_NumeratorPtBarrel";
4749 title = labelname+
"NumeratorPtBarrel;Calo Pt[GeV/c]";
4750 MonitorElement * NumeratorPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4751 NumeratorPtBarrel->
getTH1();
4753 histoname = labelname+
"_NumeratorPtEndcap";
4754 title = labelname+
"NumeratorPtEndcap;Calo Pt[GeV/c]";
4755 MonitorElement * NumeratorPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4756 NumeratorPtEndcap->
getTH1();
4758 histoname = labelname+
"_NumeratorPtForward";
4759 title = labelname+
"NumeratorPtForward;Calo Pt[GeV/c]";
4760 MonitorElement * NumeratorPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4761 NumeratorPtForward->
getTH1();
4763 histoname = labelname+
"_NumeratorEta";
4764 title = labelname+
"NumeratorEta;Calo #eta ";
4765 MonitorElement * NumeratorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4768 histoname = labelname+
"_NumeratorPhi";
4769 title = labelname+
"NumeratorPhi;Calo #Phi";
4770 MonitorElement * NumeratorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4773 histoname = labelname+
"_NumeratorEtaPhi";
4774 title = labelname+
"NumeratorEtaPhi;Calo #eta;Calo #Phi";
4775 MonitorElement * NumeratorEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4776 NumeratorEtaPhi->
getTH1();
4778 histoname = labelname+
"_NumeratorEtaBarrel";
4779 title = labelname+
"NumeratorEtaBarrel;Calo #eta ";
4780 MonitorElement * NumeratorEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4781 NumeratorEtaBarrel->
getTH1();
4783 histoname = labelname+
"_NumeratorPhiBarrel";
4784 title = labelname+
"NumeratorPhiBarrel;Calo #Phi";
4785 MonitorElement * NumeratorPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4786 NumeratorPhiBarrel->
getTH1();
4788 histoname = labelname+
"_NumeratorEtaEndcap";
4789 title = labelname+
"NumeratorEtaEndcap;Calo #eta ";
4790 MonitorElement * NumeratorEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4791 NumeratorEtaEndcap->
getTH1();
4793 histoname = labelname+
"_NumeratorPhiEndcap";
4794 title = labelname+
"NumeratorPhiEndcap;Calo #Phi";
4795 MonitorElement * NumeratorPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4796 NumeratorPhiEndcap->
getTH1();
4798 histoname = labelname+
"_NumeratorEtaForward";
4799 title = labelname+
"NumeratorEtaForward;Calo #eta ";
4800 MonitorElement * NumeratorEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4801 NumeratorEtaForward->
getTH1();
4803 histoname = labelname+
"_NumeratorPhiForward";
4804 title = labelname+
"NumeratorPhiForward;Calo #Phi";
4805 MonitorElement * NumeratorPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4806 NumeratorPhiForward->
getTH1();
4808 histoname = labelname+
"_NumeratorEta_LowpTcut";
4809 title = labelname+
"NumeratorEta_LowpTcut;Calo #eta ";
4810 MonitorElement * NumeratorEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4811 NumeratorEta_LowpTcut->
getTH1();
4813 histoname = labelname+
"_NumeratorPhi_LowpTcut";
4814 title = labelname+
"NumeratorPhi_LowpTcut;Calo #Phi";
4815 MonitorElement * NumeratorPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4816 NumeratorPhi_LowpTcut->
getTH1();
4818 histoname = labelname+
"_NumeratorEtaPhi_LowpTcut";
4819 title = labelname+
"NumeratorEtaPhi_LowpTcut;Calo #eta;Calo #Phi";
4820 MonitorElement * NumeratorEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4821 NumeratorEtaPhi_LowpTcut->
getTH1();
4823 histoname = labelname+
"_NumeratorEta_MedpTcut";
4824 title = labelname+
"NumeratorEta_MedpTcut;Calo #eta ";
4825 MonitorElement * NumeratorEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4826 NumeratorEta_MedpTcut->
getTH1();
4828 histoname = labelname+
"_NumeratorPhi_MedpTcut";
4829 title = labelname+
"NumeratorPhi_MedpTcut;Calo #Phi";
4830 MonitorElement * NumeratorPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4831 NumeratorPhi_MedpTcut->
getTH1();
4833 histoname = labelname+
"_NumeratorEtaPhi_MedpTcut";
4834 title = labelname+
"NumeratorEtaPhi_MedpTcut;Calo #eta;Calo #Phi";
4835 MonitorElement * NumeratorEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4836 NumeratorEtaPhi_MedpTcut->
getTH1();
4838 histoname = labelname+
"_NumeratorEta_HighpTcut";
4839 title = labelname+
"NumeratorEta_HighpTcut;Calo #eta ";
4840 MonitorElement * NumeratorEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4841 NumeratorEta_HighpTcut->
getTH1();
4843 histoname = labelname+
"_NumeratorPhi_HighpTcut";
4844 title = labelname+
"NumeratorPhi_HighpTcut;Calo #Phi";
4845 MonitorElement * NumeratorPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4846 NumeratorPhi_HighpTcut->
getTH1();
4848 histoname = labelname+
"_NumeratorEtaPhi_HighpTcut";
4849 title = labelname+
"NumeratorEtaPhi_HighpTcut;Calo #eta;Calo #Phi";
4850 MonitorElement * NumeratorEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4851 NumeratorEtaPhi_HighpTcut->
getTH1();
4853 histoname = labelname+
"_DenominatorPt";
4854 title = labelname+
"DenominatorPt;Calo Pt[GeV/c]";
4855 MonitorElement * DenominatorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4858 histoname = labelname+
"_DenominatorPtBarrel";
4859 title = labelname+
"DenominatorPtBarrel;Calo Pt[GeV/c]";
4860 MonitorElement * DenominatorPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4861 DenominatorPtBarrel->
getTH1();
4863 histoname = labelname+
"_DenominatorPtEndcap";
4864 title = labelname+
"DenominatorPtEndcap;Calo Pt[GeV/c]";
4865 MonitorElement * DenominatorPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4866 DenominatorPtEndcap->
getTH1();
4868 histoname = labelname+
"_DenominatorPtForward";
4869 title = labelname+
"DenominatorPtForward;Calo Pt[GeV/c] ";
4870 MonitorElement * DenominatorPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4871 DenominatorPtForward->
getTH1();
4873 histoname = labelname+
"_DenominatorEta";
4874 title = labelname+
"DenominatorEta;Calo #eta";
4875 MonitorElement * DenominatorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4876 DenominatorEta->
getTH1();
4878 histoname = labelname+
"_DenominatorPhi";
4879 title = labelname+
"DenominatorPhi;Calo #Phi";
4880 MonitorElement * DenominatorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4881 DenominatorPhi->
getTH1();
4883 histoname = labelname+
"_DenominatorEtaPhi";
4884 title = labelname+
"DenominatorEtaPhi;Calo #eta (IC5);Calo #Phi ";
4885 MonitorElement * DenominatorEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4886 DenominatorEtaPhi->
getTH1();
4888 histoname = labelname+
"_DenominatorEtaBarrel";
4889 title = labelname+
"DenominatorEtaBarrel;Calo #eta ";
4890 MonitorElement * DenominatorEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4891 DenominatorEtaBarrel->
getTH1();
4893 histoname = labelname+
"_DenominatorPhiBarrel";
4894 title = labelname+
"DenominatorPhiBarrel;Calo #Phi";
4895 MonitorElement * DenominatorPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4896 DenominatorPhiBarrel->
getTH1();
4898 histoname = labelname+
"_DenominatorEtaEndcap";
4899 title = labelname+
"DenominatorEtaEndcap;Calo #eta ";
4900 MonitorElement * DenominatorEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4901 DenominatorEtaEndcap->
getTH1();
4903 histoname = labelname+
"_DenominatorPhiEndcap";
4904 title = labelname+
"DenominatorPhiEndcap;Calo #Phi";
4905 MonitorElement * DenominatorPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4906 DenominatorPhiEndcap->
getTH1();
4908 histoname = labelname+
"_DenominatorEtaForward";
4909 title = labelname+
"DenominatorEtaForward;Calo #eta ";
4910 MonitorElement * DenominatorEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4911 DenominatorEtaForward->
getTH1();
4913 histoname = labelname+
"_DenominatorPhiForward";
4914 title = labelname+
"DenominatorPhiForward;Calo #Phi";
4915 MonitorElement * DenominatorPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4916 DenominatorPhiForward->
getTH1();
4918 histoname = labelname+
"_DenominatorEta_LowpTcut";
4919 title = labelname+
"DenominatorEta_LowpTcut;Calo #eta ";
4920 MonitorElement * DenominatorEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4921 DenominatorEta_LowpTcut->
getTH1();
4923 histoname = labelname+
"_DenominatorPhi_LowpTcut";
4924 title = labelname+
"DenominatorPhi_LowpTcut;Calo #Phi";
4925 MonitorElement * DenominatorPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4926 DenominatorPhi_LowpTcut->
getTH1();
4928 histoname = labelname+
"_DenominatorEtaPhi_LowpTcut";
4929 title = labelname+
"DenominatorEtaPhi_LowpTcut;Calo #eta;Calo #Phi";
4930 MonitorElement * DenominatorEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4931 DenominatorEtaPhi_LowpTcut->
getTH1();
4933 histoname = labelname+
"_DenominatorEta_MedpTcut";
4934 title = labelname+
"DenominatorEta_MedpTcut;Calo #eta ";
4935 MonitorElement * DenominatorEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4936 DenominatorEta_MedpTcut->
getTH1();
4938 histoname = labelname+
"_DenominatorPhi_MedpTcut";
4939 title = labelname+
"DenominatorPhi_MedpTcut;Calo #Phi";
4940 MonitorElement * DenominatorPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4941 DenominatorPhi_MedpTcut->
getTH1();
4943 histoname = labelname+
"_DenominatorEtaPhi_MedpTcut";
4944 title = labelname+
"DenominatorEtaPhi_MedpTcut;Calo #eta;Calo #Phi";
4945 MonitorElement * DenominatorEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4946 DenominatorEtaPhi_MedpTcut->
getTH1();
4948 histoname = labelname+
"_DenominatorEta_HighpTcut";
4949 title = labelname+
"DenominatorEta_HighpTcut;Calo #eta ";
4950 MonitorElement * DenominatorEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
4951 DenominatorEta_HighpTcut->
getTH1();
4953 histoname = labelname+
"_DenominatorPhi_HighpTcut";
4954 title = labelname+
"DenominatorPhi_HighpTcut;Calo #Phi";
4955 MonitorElement * DenominatorPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
4956 DenominatorPhi_HighpTcut->
getTH1();
4958 histoname = labelname+
"_DenominatorEtaPhi_HighpTcut";
4959 title = labelname+
"DenominatorEtaPhi_HighpTcut;Calo #eta;Calo #Phi";
4960 MonitorElement * DenominatorEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
4961 DenominatorEtaPhi_HighpTcut->
getTH1();
4963 histoname = labelname+
"_DeltaR";
4964 title = labelname+
"DeltaR;";
4968 histoname = labelname+
"_DeltaPhi";
4969 title = labelname+
"DeltaPhi;";
4973 histoname = labelname+
"_NumeratorPFPt";
4974 title = labelname+
"NumeratorPFPt;PFPt[GeV/c]";
4975 MonitorElement * NumeratorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4978 histoname = labelname+
"_NumeratorPFMHT";
4979 title = labelname+
"NumeratorPFMHT;PFMHT[GeV/c]";
4980 MonitorElement * NumeratorPFMHT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4981 NumeratorPFMHT->
getTH1();
4983 histoname = labelname+
"_NumeratorPFPtBarrel";
4984 title = labelname+
"NumeratorPFPtBarrel;PF Pt[GeV/c]";
4985 MonitorElement * NumeratorPFPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4986 NumeratorPFPtBarrel->
getTH1();
4988 histoname = labelname+
"_NumeratorPFPtEndcap";
4989 title = labelname+
"NumeratorPFPtEndcap;PF Pt[GeV/c]";
4990 MonitorElement * NumeratorPFPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4991 NumeratorPFPtEndcap->
getTH1();
4993 histoname = labelname+
"_NumeratorPFPtForward";
4994 title = labelname+
"NumeratorPFPtForward;PF Pt[GeV/c]";
4995 MonitorElement * NumeratorPFPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
4996 NumeratorPFPtForward->
getTH1();
4998 histoname = labelname+
"_NumeratorPFEta";
4999 title = labelname+
"NumeratorPFEta;PF #eta ";
5000 MonitorElement * NumeratorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5001 NumeratorPFEta->
getTH1();
5003 histoname = labelname+
"_NumeratorPFPhi";
5004 title = labelname+
"NumeratorPFPhi;PF #Phi";
5005 MonitorElement * NumeratorPFPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5006 NumeratorPFPhi->
getTH1();
5008 histoname = labelname+
"_NumeratorPFEtaPhi";
5009 title = labelname+
"NumeratorPFEtaPhi;PF #eta;PF #Phi";
5010 MonitorElement * NumeratorPFEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5011 NumeratorPFEtaPhi->
getTH1();
5013 histoname = labelname+
"_NumeratorPFEtaBarrel";
5014 title = labelname+
"NumeratorPFEtaBarrel;PF #eta ";
5015 MonitorElement * NumeratorPFEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5016 NumeratorPFEtaBarrel->
getTH1();
5018 histoname = labelname+
"_NumeratorPFPhiBarrel";
5019 title = labelname+
"NumeratorPFPhiBarrel;PF #Phi";
5020 MonitorElement * NumeratorPFPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5021 NumeratorPFPhiBarrel->
getTH1();
5023 histoname = labelname+
"_NumeratorPFEtaEndcap";
5024 title = labelname+
"NumeratorPFEtaEndcap;PF #eta ";
5025 MonitorElement * NumeratorPFEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5026 NumeratorPFEtaEndcap->
getTH1();
5028 histoname = labelname+
"_NumeratorPFPhiEndcap";
5029 title = labelname+
"NumeratorPFPhiEndcap;PF #Phi";
5030 MonitorElement * NumeratorPFPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5031 NumeratorPFPhiEndcap->
getTH1();
5033 histoname = labelname+
"_NumeratorPFEtaForward";
5034 title = labelname+
"NumeratorPFEtaForward;PF #eta ";
5035 MonitorElement * NumeratorPFEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5036 NumeratorPFEtaForward->
getTH1();
5038 histoname = labelname+
"_NumeratorPFPhiForward";
5039 title = labelname+
"NumeratorPFPhiForward;PF #Phi";
5040 MonitorElement * NumeratorPFPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5041 NumeratorPFPhiForward->
getTH1();
5043 histoname = labelname+
"_NumeratorPFEta_LowpTcut";
5044 title = labelname+
"NumeratorPFEta_LowpTcut;PF #eta ";
5045 MonitorElement * NumeratorPFEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5046 NumeratorPFEta_LowpTcut->
getTH1();
5048 histoname = labelname+
"_NumeratorPFPhi_LowpTcut";
5049 title = labelname+
"NumeratorPFPhi_LowpTcut;PF #Phi";
5050 MonitorElement * NumeratorPFPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5051 NumeratorPFPhi_LowpTcut->
getTH1();
5053 histoname = labelname+
"_NumeratorPFEtaPhi_LowpTcut";
5054 title = labelname+
"NumeratorPFEtaPhi_LowpTcut;PF #eta;PF #Phi";
5055 MonitorElement * NumeratorPFEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5056 NumeratorPFEtaPhi_LowpTcut->
getTH1();
5058 histoname = labelname+
"_NumeratorPFEta_MedpTcut";
5059 title = labelname+
"NumeratorPFEta_MedpTcut;PF #eta ";
5060 MonitorElement * NumeratorPFEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5061 NumeratorPFEta_MedpTcut->
getTH1();
5063 histoname = labelname+
"_NumeratorPFPhi_MedpTcut";
5064 title = labelname+
"NumeratorPFPhi_MedpTcut;PF #Phi";
5065 MonitorElement * NumeratorPFPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5066 NumeratorPFPhi_MedpTcut->
getTH1();
5068 histoname = labelname+
"_NumeratorPFEtaPhi_MedpTcut";
5069 title = labelname+
"NumeratorPFEtaPhi_MedpTcut;PF #eta;PF #Phi";
5070 MonitorElement * NumeratorPFEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5071 NumeratorPFEtaPhi_MedpTcut->
getTH1();
5073 histoname = labelname+
"_NumeratorPFEta_HighpTcut";
5074 title = labelname+
"NumeratorPFEta_HighpTcut;PF #eta ";
5075 MonitorElement * NumeratorPFEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5076 NumeratorPFEta_HighpTcut->
getTH1();
5078 histoname = labelname+
"_NumeratorPFPhi_HighpTcut";
5079 title = labelname+
"NumeratorPFPhi_HighpTcut;PF #Phi";
5080 MonitorElement * NumeratorPFPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5081 NumeratorPFPhi_HighpTcut->
getTH1();
5083 histoname = labelname+
"_NumeratorPFEtaPhi_HighpTcut";
5084 title = labelname+
"NumeratorPFEtaPhi_HighpTcut;PF #eta;PF #Phi";
5085 MonitorElement * NumeratorPFEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5086 NumeratorPFEtaPhi_HighpTcut->
getTH1();
5088 histoname = labelname+
"_DenominatorPFMHT";
5089 title = labelname+
"DenominatorPFMHT;PF Pt[GeV/c]";
5090 MonitorElement * DenominatorPFMHT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5091 DenominatorPFMHT->
getTH1();
5093 histoname = labelname+
"_DenominatorPFPt";
5094 title = labelname+
"DenominatorPFPt;PF Pt[GeV/c]";
5095 MonitorElement * DenominatorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5096 DenominatorPFPt->
getTH1();
5098 histoname = labelname+
"_DenominatorPFPtBarrel";
5099 title = labelname+
"DenominatorPFPtBarrel;PF Pt[GeV/c]";
5100 MonitorElement * DenominatorPFPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5101 DenominatorPFPtBarrel->
getTH1();
5103 histoname = labelname+
"_DenominatorPFPtEndcap";
5104 title = labelname+
"DenominatorPFPtEndcap;PF Pt[GeV/c]";
5105 MonitorElement * DenominatorPFPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5106 DenominatorPFPtEndcap->
getTH1();
5108 histoname = labelname+
"_DenominatorPFPtForward";
5109 title = labelname+
"DenominatorPFPtForward;PF Pt[GeV/c] ";
5110 MonitorElement * DenominatorPFPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5111 DenominatorPFPtForward->
getTH1();
5113 histoname = labelname+
"_DenominatorPFEta";
5114 title = labelname+
"DenominatorPFEta;PF #eta";
5115 MonitorElement * DenominatorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5116 DenominatorPFEta->
getTH1();
5118 histoname = labelname+
"_DenominatorPFPhi";
5119 title = labelname+
"DenominatorPFPhi;PF #Phi";
5120 MonitorElement * DenominatorPFPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5121 DenominatorPFPhi->
getTH1();
5123 histoname = labelname+
"_DenominatorPFEtaPhi";
5124 title = labelname+
"DenominatorPFEtaPhi;PF #eta (IC5);PF #Phi ";
5125 MonitorElement * DenominatorPFEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5126 DenominatorPFEtaPhi->
getTH1();
5128 histoname = labelname+
"_DenominatorPFEtaBarrel";
5129 title = labelname+
"DenominatorPFEtaBarrel;PF #eta ";
5130 MonitorElement * DenominatorPFEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5131 DenominatorPFEtaBarrel->
getTH1();
5133 histoname = labelname+
"_DenominatorPFPhiBarrel";
5134 title = labelname+
"DenominatorPFPhiBarrel;PF #Phi";
5135 MonitorElement * DenominatorPFPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5136 DenominatorPFPhiBarrel->
getTH1();
5138 histoname = labelname+
"_DenominatorPFEtaEndcap";
5139 title = labelname+
"DenominatorPFEtaEndcap;PF #eta ";
5140 MonitorElement * DenominatorPFEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5141 DenominatorPFEtaEndcap->
getTH1();
5143 histoname = labelname+
"_DenominatorPFPhiEndcap";
5144 title = labelname+
"DenominatorPFPhiEndcap;PF #Phi";
5145 MonitorElement * DenominatorPFPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5146 DenominatorPFPhiEndcap->
getTH1();
5148 histoname = labelname+
"_DenominatorPFEtaForward";
5149 title = labelname+
"DenominatorPFEtaForward;PF #eta ";
5150 MonitorElement * DenominatorPFEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5151 DenominatorPFEtaForward->
getTH1();
5153 histoname = labelname+
"_DenominatorPFPhiForward";
5154 title = labelname+
"DenominatorPFPhiForward;PF #Phi";
5155 MonitorElement * DenominatorPFPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5156 DenominatorPFPhiForward->
getTH1();
5158 histoname = labelname+
"_DenominatorPFEta_LowpTcut";
5159 title = labelname+
"DenominatorPFEta_LowpTcut;PF #eta ";
5160 MonitorElement * DenominatorPFEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5161 DenominatorPFEta_LowpTcut->
getTH1();
5163 histoname = labelname+
"_DenominatorPFPhi_LowpTcut";
5164 title = labelname+
"DenominatorPFPhi_LowpTcut;PF #Phi";
5165 MonitorElement * DenominatorPFPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5166 DenominatorPFPhi_LowpTcut->
getTH1();
5168 histoname = labelname+
"_DenominatorPFEtaPhi_LowpTcut";
5169 title = labelname+
"DenominatorPFEtaPhi_LowpTcut;PF #eta;PF #Phi";
5170 MonitorElement * DenominatorPFEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5171 DenominatorPFEtaPhi_LowpTcut->
getTH1();
5173 histoname = labelname+
"_DenominatorPFEta_MedpTcut";
5174 title = labelname+
"DenominatorPFEta_MedpTcut;PF #eta ";
5175 MonitorElement * DenominatorPFEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5176 DenominatorPFEta_MedpTcut->
getTH1();
5178 histoname = labelname+
"_DenominatorPFPhi_MedpTcut";
5179 title = labelname+
"DenominatorPFPhi_MedpTcut;PF #Phi";
5180 MonitorElement * DenominatorPFPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5181 DenominatorPFPhi_MedpTcut->
getTH1();
5183 histoname = labelname+
"_DenominatorPFEtaPhi_MedpTcut";
5184 title = labelname+
"DenominatorPFEtaPhi_MedpTcut;PF #eta;PF #Phi";
5185 MonitorElement * DenominatorPFEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5186 DenominatorPFEtaPhi_MedpTcut->
getTH1();
5188 histoname = labelname+
"_DenominatorPFEta_HighpTcut";
5189 title = labelname+
"DenominatorPFEta_HighpTcut;PF #eta ";
5190 MonitorElement * DenominatorPFEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5191 DenominatorPFEta_HighpTcut->
getTH1();
5193 histoname = labelname+
"_DenominatorPFPhi_HighpTcut";
5194 title = labelname+
"DenominatorPFPhi_HighpTcut;PF #Phi";
5195 MonitorElement * DenominatorPFPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5196 DenominatorPFPhi_HighpTcut->
getTH1();
5198 histoname = labelname+
"_DenominatorPFEtaPhi_HighpTcut";
5199 title = labelname+
"DenominatorPFEtaPhi_HighpTcut;PF #eta;PF #Phi";
5200 MonitorElement * DenominatorPFEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5201 DenominatorPFEtaPhi_HighpTcut->
getTH1();
5203 histoname = labelname+
"_PFDeltaR";
5204 title = labelname+
"PFDeltaR;";
5208 histoname = labelname+
"_PFDeltaPhi";
5209 title = labelname+
"PFDeltaPhi;";
5213 v->setEffHistos(NumeratorPt,
5225 NumeratorEtaForward,
5226 NumeratorPhiForward,
5227 NumeratorEta_LowpTcut,
5228 NumeratorPhi_LowpTcut,
5229 NumeratorEtaPhi_LowpTcut,
5230 NumeratorEta_MedpTcut,
5231 NumeratorPhi_MedpTcut,
5232 NumeratorEtaPhi_MedpTcut,
5233 NumeratorEta_HighpTcut,
5234 NumeratorPhi_HighpTcut,
5235 NumeratorEtaPhi_HighpTcut,
5238 DenominatorPtBarrel,
5239 DenominatorPtEndcap,
5240 DenominatorPtForward,
5245 DenominatorEtaBarrel,
5246 DenominatorPhiBarrel,
5247 DenominatorEtaEndcap,
5248 DenominatorPhiEndcap,
5249 DenominatorEtaForward,
5250 DenominatorPhiForward,
5251 DenominatorEta_LowpTcut,
5252 DenominatorPhi_LowpTcut,
5253 DenominatorEtaPhi_LowpTcut,
5254 DenominatorEta_MedpTcut,
5255 DenominatorPhi_MedpTcut,
5256 DenominatorEtaPhi_MedpTcut,
5257 DenominatorEta_HighpTcut,
5258 DenominatorPhi_HighpTcut,
5259 DenominatorEtaPhi_HighpTcut,
5266 NumeratorPFPtBarrel,
5267 NumeratorPFPtEndcap,
5268 NumeratorPFPtForward,
5272 NumeratorPFEtaBarrel,
5273 NumeratorPFPhiBarrel,
5274 NumeratorPFEtaEndcap,
5275 NumeratorPFPhiEndcap,
5276 NumeratorPFEtaForward,
5277 NumeratorPFPhiForward,
5278 NumeratorPFEta_LowpTcut,
5279 NumeratorPFPhi_LowpTcut,
5280 NumeratorPFEtaPhi_LowpTcut,
5281 NumeratorPFEta_MedpTcut,
5282 NumeratorPFPhi_MedpTcut,
5283 NumeratorPFEtaPhi_MedpTcut,
5284 NumeratorPFEta_HighpTcut,
5285 NumeratorPFPhi_HighpTcut,
5286 NumeratorPFEtaPhi_HighpTcut,
5289 DenominatorPFPtBarrel,
5290 DenominatorPFPtEndcap,
5291 DenominatorPFPtForward,
5294 DenominatorPFEtaPhi,
5295 DenominatorPFEtaBarrel,
5296 DenominatorPFPhiBarrel,
5297 DenominatorPFEtaEndcap,
5298 DenominatorPFPhiEndcap,
5299 DenominatorPFEtaForward,
5300 DenominatorPFPhiForward,
5301 DenominatorPFEta_LowpTcut,
5302 DenominatorPFPhi_LowpTcut,
5303 DenominatorPFEtaPhi_LowpTcut,
5304 DenominatorPFEta_MedpTcut,
5305 DenominatorPFPhi_MedpTcut,
5306 DenominatorPFEtaPhi_MedpTcut,
5307 DenominatorPFEta_HighpTcut,
5308 DenominatorPFPhi_HighpTcut,
5309 DenominatorPFEtaPhi_HighpTcut,
5315 if((
v->getObjectType() ==
trigger::TriggerJet) && (
v->getTriggerType().compare(
"DiJet_Trigger") == 0)){
5317 histoname = labelname+
"_NumeratorAvrgPt";
5318 title = labelname+
"NumeratorAvrgPt;Calo Pt[GeV/c] ";
5319 MonitorElement * NumeratorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5322 histoname = labelname+
"_NumeratorAvrgEta";
5323 title = labelname+
"NumeratorAvrgEta;Calo #eta";
5324 MonitorElement * NumeratorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5327 histoname = labelname+
"_DenominatorAvrgPt";
5328 title = labelname+
"DenominatorAvrgPt;Calo Pt[GeV/c]";
5329 MonitorElement * DenominatorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5332 histoname = labelname+
"_DenominatorAvrgEta";
5333 title = labelname+
"DenominatorAvrgEta;Calo #eta ";
5334 MonitorElement * DenominatorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5335 DenominatorEta->
getTH1();
5337 histoname = labelname+
"_DeltaR";
5338 title = labelname+
"DeltaR;";
5342 histoname = labelname+
"_DeltaPhi";
5343 title = labelname+
"DeltaPhi;";
5347 histoname = labelname+
"_NumeratorAvrgPFPt";
5348 title = labelname+
"NumeratorAvrgPFPt;PF Pt[GeV/c] ";
5349 MonitorElement * NumeratorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5352 histoname = labelname+
"_NumeratorAvrgPFEta";
5353 title = labelname+
"NumeratorAvrgPFEta;PF #eta";
5354 MonitorElement * NumeratorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5355 NumeratorPFEta->
getTH1();
5357 histoname = labelname+
"_DenominatorAvrgPFPt";
5358 title = labelname+
"DenominatorAvrgPFPt;PF Pt[GeV/c]";
5359 MonitorElement * DenominatorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5360 DenominatorPFPt->
getTH1();
5362 histoname = labelname+
"_DenominatorAvrgPFEta";
5363 title = labelname+
"DenominatorAvrgPFEta;PF #eta ";
5364 MonitorElement * DenominatorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5365 DenominatorPFEta->
getTH1();
5367 histoname = labelname+
"_PFDeltaR";
5368 title = labelname+
"PFDeltaR;";
5372 histoname = labelname+
"_PFDeltaPhi";
5373 title = labelname+
"PFDeltaPhi;";
5377 v->setEffHistos( dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5378 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5379 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5380 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5381 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
5392 histoname = labelname+
"_NumeratorPt";
5393 if(
v->getPath().find(
"HLT_PFMET")==std::string::npos)
5394 title = labelname+
"NumeratorPt; CaloMET[GeV/c]";
5396 title = labelname+
"NumeratorPt; PFMET[GeV/c]";
5397 MonitorElement * NumeratorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5400 histoname = labelname+
"_NumeratorPhi";
5401 title = labelname+
"NumeratorPhi; #Phi";
5402 MonitorElement * NumeratorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5405 histoname = labelname+
"_DenominatorPt";
5406 if(
v->getPath().find(
"HLT_PFMET")==std::string::npos)
5407 title = labelname+
"DenominatorPt; CaloMET[GeV/c]";
5409 title = labelname+
"DenominatorPt; PFMET[GeV/c]";
5410 MonitorElement * DenominatorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5413 histoname = labelname+
"_DenominatorPhi";
5414 title = labelname+
"DenominatorPhi; #Phi";
5415 MonitorElement * DenominatorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5416 DenominatorPhi->
getTH1();
5418 v->setEffHistos( NumeratorPt, dummy, dummy, dummy, dummy, NumeratorPhi, dummy, dummy, dummy, dummy,
5419 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5420 dummy, dummy, DenominatorPt, dummy, dummy, dummy, dummy, DenominatorPhi, dummy, dummy,
5421 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5422 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5423 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
5424 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, 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
5449 if((
v->getObjectType() ==
trigger::TriggerJet) && (
v->getTriggerType().compare(
"SingleJet_Trigger") == 0)){
5451 histoname = labelname+
"_NumeratorPt";
5452 title = labelname+
"NumeratorPt;Calo Pt[GeV/c] ";
5453 MonitorElement * NumeratorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5456 histoname = labelname+
"_NumeratorPtBarrel";
5457 title = labelname+
"NumeratorPtBarrel;Calo Pt[GeV/c]";
5458 MonitorElement * NumeratorPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5459 NumeratorPtBarrel->
getTH1();
5461 histoname = labelname+
"_NumeratorPtEndcap";
5462 title = labelname+
"NumeratorPtEndcap; Calo Pt[GeV/c] ";
5463 MonitorElement * NumeratorPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5464 NumeratorPtEndcap->
getTH1();
5466 histoname = labelname+
"_NumeratorPtForward";
5467 title = labelname+
"NumeratorPtForward;Calo Pt[GeV/c]";
5468 MonitorElement * NumeratorPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5469 NumeratorPtForward->
getTH1();
5471 histoname = labelname+
"_NumeratorEta";
5472 title = labelname+
"NumeratorEta;Calo #eta ";
5473 MonitorElement * NumeratorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5476 histoname = labelname+
"_NumeratorPhi";
5477 title = labelname+
"NumeratorPhi;Calo #Phi";
5478 MonitorElement * NumeratorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5481 histoname = labelname+
"_NumeratorEtaPhi";
5482 title = labelname+
"NumeratorEtaPhi;Calo #eta;Calo #Phi ";
5483 MonitorElement * NumeratorEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5484 NumeratorEtaPhi->
getTH1();
5486 histoname = labelname+
"_NumeratorEtaBarrel";
5487 title = labelname+
"NumeratorEtaBarrel;Calo #eta ";
5488 MonitorElement * NumeratorEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5489 NumeratorEtaBarrel->
getTH1();
5491 histoname = labelname+
"_NumeratorPhiBarrel";
5492 title = labelname+
"NumeratorPhiBarrel;Calo #Phi";
5493 MonitorElement * NumeratorPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5494 NumeratorPhiBarrel->
getTH1();
5496 histoname = labelname+
"_NumeratorEtaEndcap";
5497 title = labelname+
"NumeratorEtaEndcap;Calo #eta ";
5498 MonitorElement * NumeratorEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5499 NumeratorEtaEndcap->
getTH1();
5501 histoname = labelname+
"_NumeratorPhiEndcap";
5502 title = labelname+
"NumeratorPhiEndcap;Calo #Phi";
5503 MonitorElement * NumeratorPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5504 NumeratorPhiEndcap->
getTH1();
5506 histoname = labelname+
"_NumeratorEtaForward";
5507 title = labelname+
"NumeratorEtaForward;Calo #eta ";
5508 MonitorElement * NumeratorEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5509 NumeratorEtaForward->
getTH1();
5511 histoname = labelname+
"_NumeratorPhiForward";
5512 title = labelname+
"NumeratorPhiForward;Calo #Phi";
5513 MonitorElement * NumeratorPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5514 NumeratorPhiForward->
getTH1();
5516 histoname = labelname+
"_NumeratorEta_LowpTcut";
5517 title = labelname+
"NumeratorEta_LowpTcut;Calo #eta ";
5518 MonitorElement * NumeratorEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5519 NumeratorEta_LowpTcut->
getTH1();
5521 histoname = labelname+
"_NumeratorPhi_LowpTcut";
5522 title = labelname+
"NumeratorPhi_LowpTcut;Calo #Phi";
5523 MonitorElement * NumeratorPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5524 NumeratorPhi_LowpTcut->
getTH1();
5526 histoname = labelname+
"_NumeratorEtaPhi_LowpTcut";
5527 title = labelname+
"NumeratorEtaPhi_LowpTcut;Calo #eta;Calo #Phi";
5528 MonitorElement * NumeratorEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5529 NumeratorEtaPhi_LowpTcut->
getTH1();
5531 histoname = labelname+
"_NumeratorEta_MedpTcut";
5532 title = labelname+
"NumeratorEta_MedpTcut;Calo #eta ";
5533 MonitorElement * NumeratorEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5534 NumeratorEta_MedpTcut->
getTH1();
5536 histoname = labelname+
"_NumeratorPhi_MedpTcut";
5537 title = labelname+
"NumeratorPhi_MedpTcut;Calo #Phi";
5538 MonitorElement * NumeratorPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5539 NumeratorPhi_MedpTcut->
getTH1();
5541 histoname = labelname+
"_NumeratorEtaPhi_MedpTcut";
5542 title = labelname+
"NumeratorEtaPhi_MedpTcut;Calo #eta;Calo #Phi";
5543 MonitorElement * NumeratorEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5544 NumeratorEtaPhi_MedpTcut->
getTH1();
5546 histoname = labelname+
"_NumeratorEta_HighpTcut";
5547 title = labelname+
"NumeratorEta_HighpTcut;Calo #eta ";
5548 MonitorElement * NumeratorEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5549 NumeratorEta_HighpTcut->
getTH1();
5551 histoname = labelname+
"_NumeratorPhi_HighpTcut";
5552 title = labelname+
"NumeratorPhi_HighpTcut;Calo #Phi";
5553 MonitorElement * NumeratorPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5554 NumeratorPhi_HighpTcut->
getTH1();
5556 histoname = labelname+
"_NumeratorEtaPhi_HighpTcut";
5557 title = labelname+
"NumeratorEtaPhi_HighpTcut;Calo #eta;Calo #Phi";
5558 MonitorElement * NumeratorEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5559 NumeratorEtaPhi_HighpTcut->
getTH1();
5561 histoname = labelname+
"_DenominatorPt";
5562 title = labelname+
"DenominatorPt;Calo Pt[GeV/c]";
5563 MonitorElement * DenominatorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5566 histoname = labelname+
"_DenominatorPtBarrel";
5567 title = labelname+
"DenominatorPtBarrel;Calo Pt[GeV/c]";
5568 MonitorElement * DenominatorPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5569 DenominatorPtBarrel->
getTH1();
5571 histoname = labelname+
"_DenominatorPtEndcap";
5572 title = labelname+
"DenominatorPtEndcap;Calo Pt[GeV/c]";
5573 MonitorElement * DenominatorPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5574 DenominatorPtEndcap->
getTH1();
5576 histoname = labelname+
"_DenominatorPtForward";
5577 title = labelname+
"DenominatorPtForward;Calo Pt[GeV/c]";
5578 MonitorElement * DenominatorPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5579 DenominatorPtForward->
getTH1();
5581 histoname = labelname+
"_DenominatorEta";
5582 title = labelname+
"DenominatorEta;Calo #eta ";
5583 MonitorElement * DenominatorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5584 DenominatorEta->
getTH1();
5586 histoname = labelname+
"_DenominatorPhi";
5587 title = labelname+
"DenominatorPhi;Calo #Phi";
5588 MonitorElement * DenominatorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5589 DenominatorPhi->
getTH1();
5591 histoname = labelname+
"_DenominatorEtaPhi";
5592 title = labelname+
"DenominatorEtaPhi;Calo #eta ;Calo #Phi ";
5593 MonitorElement * DenominatorEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5594 DenominatorEtaPhi->
getTH1();
5596 histoname = labelname+
"_DenominatorEtaBarrel";
5597 title = labelname+
"DenominatorEtaBarrel;Calo #eta ";
5598 MonitorElement * DenominatorEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5599 DenominatorEtaBarrel->
getTH1();
5601 histoname = labelname+
"_DenominatorPhiBarrel";
5602 title = labelname+
"DenominatorPhiBarrel;Calo #Phi";
5603 MonitorElement * DenominatorPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5604 DenominatorPhiBarrel->
getTH1();
5606 histoname = labelname+
"_DenominatorEtaEndcap";
5607 title = labelname+
"DenominatorEtaEndcap;Calo #eta ";
5608 MonitorElement * DenominatorEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5609 DenominatorEtaEndcap->
getTH1();
5611 histoname = labelname+
"_DenominatorPhiEndcap";
5612 title = labelname+
"DenominatorPhiEndcap;Calo #Phi";
5613 MonitorElement * DenominatorPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5614 DenominatorPhiEndcap->
getTH1();
5616 histoname = labelname+
"_DenominatorEtaForward";
5617 title = labelname+
"DenominatorEtaForward;Calo #eta ";
5618 MonitorElement * DenominatorEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5619 DenominatorEtaForward->
getTH1();
5621 histoname = labelname+
"_DenominatorPhiForward";
5622 title = labelname+
"DenominatorPhiForward;Calo #Phi";
5623 MonitorElement * DenominatorPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5624 DenominatorPhiForward->
getTH1();
5626 histoname = labelname+
"_DenominatorEta_LowpTcut";
5627 title = labelname+
"DenominatorEta_LowpTcut;Calo #eta ";
5628 MonitorElement * DenominatorEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5629 DenominatorEta_LowpTcut->
getTH1();
5631 histoname = labelname+
"_DenominatorPhi_LowpTcut";
5632 title = labelname+
"DenominatorPhi_LowpTcut;Calo #Phi";
5633 MonitorElement * DenominatorPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5634 DenominatorPhi_LowpTcut->
getTH1();
5636 histoname = labelname+
"_DenominatorEtaPhi_LowpTcut";
5637 title = labelname+
"DenominatorEtaPhi_LowpTcut;Calo #eta;Calo #Phi";
5638 MonitorElement * DenominatorEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5639 DenominatorEtaPhi_LowpTcut->
getTH1();
5641 histoname = labelname+
"_DenominatorEta_MedpTcut";
5642 title = labelname+
"DenominatorEta_MedpTcut;Calo #eta ";
5643 MonitorElement * DenominatorEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5644 DenominatorEta_MedpTcut->
getTH1();
5646 histoname = labelname+
"_DenominatorPhi_MedpTcut";
5647 title = labelname+
"DenominatorPhi_MedpTcut;Calo #Phi";
5648 MonitorElement * DenominatorPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5649 DenominatorPhi_MedpTcut->
getTH1();
5651 histoname = labelname+
"_DenominatorEtaPhi_MedpTcut";
5652 title = labelname+
"DenominatorEtaPhi_MedpTcut;Calo #eta;Calo #Phi";
5653 MonitorElement * DenominatorEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5654 DenominatorEtaPhi_MedpTcut->
getTH1();
5656 histoname = labelname+
"_DenominatorEta_HighpTcut";
5657 title = labelname+
"DenominatorEta_HighpTcut;Calo #eta ";
5658 MonitorElement * DenominatorEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5659 DenominatorEta_HighpTcut->
getTH1();
5661 histoname = labelname+
"_DenominatorPhi_HighpTcut";
5662 title = labelname+
"DenominatorPhi_HighpTcut;Calo #Phi";
5663 MonitorElement * DenominatorPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5664 DenominatorPhi_HighpTcut->
getTH1();
5666 histoname = labelname+
"_DenominatorEtaPhi_HighpTcut";
5667 title = labelname+
"DenominatorEtaPhi_HighpTcut;Calo #eta;Calo #Phi";
5668 MonitorElement * DenominatorEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5669 DenominatorEtaPhi_HighpTcut->
getTH1();
5671 histoname = labelname+
"_DeltaR";
5672 title = labelname+
"DeltaR;";
5676 histoname = labelname+
"_DeltaPhi";
5677 title = labelname+
"DeltaPhi;";
5681 histoname = labelname+
"_NumeratorPFPt";
5682 title = labelname+
"NumeratorPFPt;PF Pt[GeV/c] ";
5683 MonitorElement * NumeratorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5686 histoname = labelname+
"_NumeratorPFMHT";
5687 title = labelname+
"NumeratorPFPt;PFMHT[GeV/c] ";
5688 MonitorElement * NumeratorPFMHT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5689 NumeratorPFMHT->
getTH1();
5691 histoname = labelname+
"_NumeratorPFPtBarrel";
5692 title = labelname+
"NumeratorPFPtBarrel;PF Pt[GeV/c]";
5693 MonitorElement * NumeratorPFPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5694 NumeratorPFPtBarrel->
getTH1();
5696 histoname = labelname+
"_NumeratorPFPtEndcap";
5697 title = labelname+
"NumeratorPFPtEndcap; PF Pt[GeV/c] ";
5698 MonitorElement * NumeratorPFPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5699 NumeratorPFPtEndcap->
getTH1();
5701 histoname = labelname+
"_NumeratorPFPtForward";
5702 title = labelname+
"NumeratorPFPtForward;PF Pt[GeV/c]";
5703 MonitorElement * NumeratorPFPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5704 NumeratorPFPtForward->
getTH1();
5706 histoname = labelname+
"_NumeratorPFEta";
5707 title = labelname+
"NumeratorPFEta;PF #eta ";
5708 MonitorElement * NumeratorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5709 NumeratorPFEta->
getTH1();
5711 histoname = labelname+
"_NumeratorPFPhi";
5712 title = labelname+
"NumeratorPFPhi;PF #Phi";
5713 MonitorElement * NumeratorPFPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5714 NumeratorPFPhi->
getTH1();
5716 histoname = labelname+
"_NumeratorPFEtaPhi";
5717 title = labelname+
"NumeratorPFEtaPhi;PF #eta;PF #Phi ";
5718 MonitorElement * NumeratorPFEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5719 NumeratorPFEtaPhi->
getTH1();
5721 histoname = labelname+
"_NumeratorPFEtaBarrel";
5722 title = labelname+
"NumeratorPFEtaBarrel;PF #eta ";
5723 MonitorElement * NumeratorPFEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5724 NumeratorPFEtaBarrel->
getTH1();
5726 histoname = labelname+
"_NumeratorPFPhiBarrel";
5727 title = labelname+
"NumeratorPFPhiBarrel;PF #Phi";
5728 MonitorElement * NumeratorPFPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5729 NumeratorPFPhiBarrel->
getTH1();
5731 histoname = labelname+
"_NumeratorPFEtaEndcap";
5732 title = labelname+
"NumeratorPFEtaEndcap;PF #eta ";
5733 MonitorElement * NumeratorPFEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5734 NumeratorPFEtaEndcap->
getTH1();
5736 histoname = labelname+
"_NumeratorPFPhiEndcap";
5737 title = labelname+
"NumeratorPFPhiEndcap;PF #Phi";
5738 MonitorElement * NumeratorPFPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5739 NumeratorPFPhiEndcap->
getTH1();
5741 histoname = labelname+
"_NumeratorPFEtaForward";
5742 title = labelname+
"NumeratorPFEtaForward;PF #eta ";
5743 MonitorElement * NumeratorPFEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5744 NumeratorPFEtaForward->
getTH1();
5746 histoname = labelname+
"_NumeratorPFPhiForward";
5747 title = labelname+
"NumeratorPFPhiForward;PF #Phi";
5748 MonitorElement * NumeratorPFPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5749 NumeratorPFPhiForward->
getTH1();
5751 histoname = labelname+
"_NumeratorPFEta_LowpTcut";
5752 title = labelname+
"NumeratorPFEta_LowpTcut;PF #eta ";
5753 MonitorElement * NumeratorPFEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5754 NumeratorPFEta_LowpTcut->
getTH1();
5756 histoname = labelname+
"_NumeratorPFPhi_LowpTcut";
5757 title = labelname+
"NumeratorPFPhi_LowpTcut;PF #Phi";
5758 MonitorElement * NumeratorPFPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5759 NumeratorPFPhi_LowpTcut->
getTH1();
5761 histoname = labelname+
"_NumeratorPFEtaPhi_LowpTcut";
5762 title = labelname+
"NumeratorPFEtaPhi_LowpTcut;PF #eta;PF #Phi";
5763 MonitorElement * NumeratorPFEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5764 NumeratorPFEtaPhi_LowpTcut->
getTH1();
5766 histoname = labelname+
"_NumeratorPFEta_MedpTcut";
5767 title = labelname+
"NumeratorPFEta_MedpTcut;PF #eta ";
5768 MonitorElement * NumeratorPFEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5769 NumeratorPFEta_MedpTcut->
getTH1();
5771 histoname = labelname+
"_NumeratorPFPhi_MedpTcut";
5772 title = labelname+
"NumeratorPFPhi_MedpTcut;PF #Phi";
5773 MonitorElement * NumeratorPFPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5774 NumeratorPFPhi_MedpTcut->
getTH1();
5776 histoname = labelname+
"_NumeratorPFEtaPhi_MedpTcut";
5777 title = labelname+
"NumeratorPFEtaPhi_MedpTcut;PF #eta;PF #Phi";
5778 MonitorElement * NumeratorPFEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5779 NumeratorPFEtaPhi_MedpTcut->
getTH1();
5781 histoname = labelname+
"_NumeratorPFEta_HighpTcut";
5782 title = labelname+
"NumeratorPFEta_HighpTcut;PF #eta ";
5783 MonitorElement * NumeratorPFEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5784 NumeratorPFEta_HighpTcut->
getTH1();
5786 histoname = labelname+
"_NumeratorPFPhi_HighpTcut";
5787 title = labelname+
"NumeratorPFPhi_HighpTcut;PF #Phi";
5788 MonitorElement * NumeratorPFPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5789 NumeratorPFPhi_HighpTcut->
getTH1();
5791 histoname = labelname+
"_NumeratorPFEtaPhi_HighpTcut";
5792 title = labelname+
"NumeratorPFEtaPhi_HighpTcut;PF #eta;PF #Phi";
5793 MonitorElement * NumeratorPFEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5794 NumeratorPFEtaPhi_HighpTcut->
getTH1();
5796 histoname = labelname+
"_DenominatorPFMHT";
5797 title = labelname+
"DenominatorPFMHT;PFMHT[GeV/c]";
5798 MonitorElement * DenominatorPFMHT = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5799 DenominatorPFMHT->
getTH1();
5801 histoname = labelname+
"_DenominatorPFPt";
5802 title = labelname+
"DenominatorPFPt;PF Pt[GeV/c]";
5803 MonitorElement * DenominatorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5804 DenominatorPFPt->
getTH1();
5806 histoname = labelname+
"_DenominatorPFPtBarrel";
5807 title = labelname+
"DenominatorPFPtBarrel;PF Pt[GeV/c]";
5808 MonitorElement * DenominatorPFPtBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5809 DenominatorPFPtBarrel->
getTH1();
5811 histoname = labelname+
"_DenominatorPFPtEndcap";
5812 title = labelname+
"DenominatorPFPtEndcap;PF Pt[GeV/c]";
5813 MonitorElement * DenominatorPFPtEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5814 DenominatorPFPtEndcap->
getTH1();
5816 histoname = labelname+
"_DenominatorPFPtForward";
5817 title = labelname+
"DenominatorPFPtForward;PF Pt[GeV/c]";
5818 MonitorElement * DenominatorPFPtForward = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
5819 DenominatorPFPtForward->
getTH1();
5821 histoname = labelname+
"_DenominatorPFEta";
5822 title = labelname+
"DenominatorPFEta;PF #eta ";
5823 MonitorElement * DenominatorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5824 DenominatorPFEta->
getTH1();
5826 histoname = labelname+
"_DenominatorPFPhi";
5827 title = labelname+
"DenominatorPFPhi;PF #Phi";
5828 MonitorElement * DenominatorPFPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5829 DenominatorPFPhi->
getTH1();
5831 histoname = labelname+
"_DenominatorPFEtaPhi";
5832 title = labelname+
"DenominatorPFEtaPhi;PF #eta ;PF #Phi ";
5833 MonitorElement * DenominatorPFEtaPhi = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5834 DenominatorPFEtaPhi->
getTH1();
5836 histoname = labelname+
"_DenominatorPFEtaBarrel";
5837 title = labelname+
"DenominatorPFEtaBarrel;PF #eta ";
5838 MonitorElement * DenominatorPFEtaBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5839 DenominatorPFEtaBarrel->
getTH1();
5841 histoname = labelname+
"_DenominatorPFPhiBarrel";
5842 title = labelname+
"DenominatorPFPhiBarrel;PF #Phi";
5843 MonitorElement * DenominatorPFPhiBarrel = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5844 DenominatorPFPhiBarrel->
getTH1();
5846 histoname = labelname+
"_DenominatorPFEtaEndcap";
5847 title = labelname+
"DenominatorPFEtaEndcap;PF #eta ";
5848 MonitorElement * DenominatorPFEtaEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5849 DenominatorPFEtaEndcap->
getTH1();
5851 histoname = labelname+
"_DenominatorPFPhiEndcap";
5852 title = labelname+
"DenominatorPFPhiEndcap;PF #Phi";
5853 MonitorElement * DenominatorPFPhiEndcap = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5854 DenominatorPFPhiEndcap->
getTH1();
5856 histoname = labelname+
"_DenominatorPFEtaForward";
5857 title = labelname+
"DenominatorPFEtaForward;PF #eta ";
5858 MonitorElement * DenominatorPFEtaForward = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5859 DenominatorPFEtaForward->
getTH1();
5861 histoname = labelname+
"_DenominatorPFPhiForward";
5862 title = labelname+
"DenominatorPFPhiForward;PF #Phi";
5863 MonitorElement * DenominatorPFPhiForward = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5864 DenominatorPFPhiForward->
getTH1();
5866 histoname = labelname+
"_DenominatorPFEta_LowpTcut";
5867 title = labelname+
"DenominatorPFEta_LowpTcut;PF #eta ";
5868 MonitorElement * DenominatorPFEta_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5869 DenominatorPFEta_LowpTcut->
getTH1();
5871 histoname = labelname+
"_DenominatorPFPhi_LowpTcut";
5872 title = labelname+
"DenominatorPFPhi_LowpTcut;PF #Phi";
5873 MonitorElement * DenominatorPFPhi_LowpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5874 DenominatorPFPhi_LowpTcut->
getTH1();
5876 histoname = labelname+
"_DenominatorPFEtaPhi_LowpTcut";
5877 title = labelname+
"DenominatorPFEtaPhi_LowpTcut;PF #eta;PF #Phi";
5878 MonitorElement * DenominatorPFEtaPhi_LowpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5879 DenominatorPFEtaPhi_LowpTcut->
getTH1();
5881 histoname = labelname+
"_DenominatorPFEta_MedpTcut";
5882 title = labelname+
"DenominatorPFEta_MedpTcut;PF #eta ";
5883 MonitorElement * DenominatorPFEta_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5884 DenominatorPFEta_MedpTcut->
getTH1();
5886 histoname = labelname+
"_DenominatorPFPhi_MedpTcut";
5887 title = labelname+
"DenominatorPFPhi_MedpTcut;PF #Phi";
5888 MonitorElement * DenominatorPFPhi_MedpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5889 DenominatorPFPhi_MedpTcut->
getTH1();
5891 histoname = labelname+
"_DenominatorPFEtaPhi_MedpTcut";
5892 title = labelname+
"DenominatorPFEtaPhi_MedpTcut;PF #eta;PF #Phi";
5893 MonitorElement * DenominatorPFEtaPhi_MedpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5894 DenominatorPFEtaPhi_MedpTcut->
getTH1();
5896 histoname = labelname+
"_DenominatorPFEta_HighpTcut";
5897 title = labelname+
"DenominatorPFEta_HighpTcut;PF #eta ";
5898 MonitorElement * DenominatorPFEta_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
5899 DenominatorPFEta_HighpTcut->
getTH1();
5901 histoname = labelname+
"_DenominatorPFPhi_HighpTcut";
5902 title = labelname+
"DenominatorPFPhi_HighpTcut;PF #Phi";
5903 MonitorElement * DenominatorPFPhi_HighpTcut = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
5904 DenominatorPFPhi_HighpTcut->
getTH1();
5906 histoname = labelname+
"_DenominatorPFEtaPhi_HighpTcut";
5907 title = labelname+
"DenominatorPFEtaPhi_HighpTcut;PF #eta;PF #Phi";
5908 MonitorElement * DenominatorPFEtaPhi_HighpTcut = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Phibins_,PhiMin_,PhiMax_);
5909 DenominatorPFEtaPhi_HighpTcut->
getTH1();
5911 histoname = labelname+
"_PFDeltaR";
5912 title = labelname+
"PFDeltaR;";
5916 histoname = labelname+
"_PFDeltaPhi";
5917 title = labelname+
"PFDeltaPhi;";
5921 v->setEffHistos(NumeratorPt,
5933 NumeratorEtaForward,
5934 NumeratorPhiForward,
5935 NumeratorEta_LowpTcut,
5936 NumeratorPhi_LowpTcut,
5937 NumeratorEtaPhi_LowpTcut,
5938 NumeratorEta_MedpTcut,
5939 NumeratorPhi_MedpTcut,
5940 NumeratorEtaPhi_MedpTcut,
5941 NumeratorEta_HighpTcut,
5942 NumeratorPhi_HighpTcut,
5943 NumeratorEtaPhi_HighpTcut,
5946 DenominatorPtBarrel,
5947 DenominatorPtEndcap,
5948 DenominatorPtForward,
5953 DenominatorEtaBarrel,
5954 DenominatorPhiBarrel,
5955 DenominatorEtaEndcap,
5956 DenominatorPhiEndcap,
5957 DenominatorEtaForward,
5958 DenominatorPhiForward,
5959 DenominatorEta_LowpTcut,
5960 DenominatorPhi_LowpTcut,
5961 DenominatorEtaPhi_LowpTcut,
5962 DenominatorEta_MedpTcut,
5963 DenominatorPhi_MedpTcut,
5964 DenominatorEtaPhi_MedpTcut,
5965 DenominatorEta_HighpTcut,
5966 DenominatorPhi_HighpTcut,
5967 DenominatorEtaPhi_HighpTcut,
5973 NumeratorPFPtBarrel,
5974 NumeratorPFPtEndcap,
5975 NumeratorPFPtForward,
5979 NumeratorPFEtaBarrel,
5980 NumeratorPFPhiBarrel,
5981 NumeratorPFEtaEndcap,
5982 NumeratorPFPhiEndcap,
5983 NumeratorPFEtaForward,
5984 NumeratorPFPhiForward,
5985 NumeratorPFEta_LowpTcut,
5986 NumeratorPFPhi_LowpTcut,
5987 NumeratorPFEtaPhi_LowpTcut,
5988 NumeratorPFEta_MedpTcut,
5989 NumeratorPFPhi_MedpTcut,
5990 NumeratorPFEtaPhi_MedpTcut,
5991 NumeratorPFEta_HighpTcut,
5992 NumeratorPFPhi_HighpTcut,
5993 NumeratorPFEtaPhi_HighpTcut,
5996 DenominatorPFPtBarrel,
5997 DenominatorPFPtEndcap,
5998 DenominatorPFPtForward,
6001 DenominatorPFEtaPhi,
6002 DenominatorPFEtaBarrel,
6003 DenominatorPFPhiBarrel,
6004 DenominatorPFEtaEndcap,
6005 DenominatorPFPhiEndcap,
6006 DenominatorPFEtaForward,
6007 DenominatorPFPhiForward,
6008 DenominatorPFEta_LowpTcut,
6009 DenominatorPFPhi_LowpTcut,
6010 DenominatorPFEtaPhi_LowpTcut,
6011 DenominatorPFEta_MedpTcut,
6012 DenominatorPFPhi_MedpTcut,
6013 DenominatorPFEtaPhi_MedpTcut,
6014 DenominatorPFEta_HighpTcut,
6015 DenominatorPFPhi_HighpTcut,
6016 DenominatorPFEtaPhi_HighpTcut,
6023 if((
v->getObjectType() ==
trigger::TriggerJet) && (
v->getTriggerType().compare(
"DiJet_Trigger") == 0)){
6025 histoname = labelname+
"_NumeratorAvrgPt";
6026 title = labelname+
"NumeratorAvrgPt;Calo Pt[GeV/c] ";
6027 MonitorElement * NumeratorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
6030 histoname = labelname+
"_NumeratorAvrgEta";
6031 title = labelname+
"NumeratorAvrgEta;Calo #eta ";
6032 MonitorElement * NumeratorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
6035 histoname = labelname+
"_DenominatorAvrgPt";
6036 title = labelname+
"DenominatorAvrgPt;Calo Pt[GeV/c] ";
6037 MonitorElement * DenominatorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
6040 histoname = labelname+
"_DenominatorAvrgEta";
6041 title = labelname+
"DenominatorAvrgEta;Calo #eta ";
6042 MonitorElement * DenominatorEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
6043 DenominatorEta->
getTH1();
6045 histoname = labelname+
"_DeltaR";
6046 title = labelname+
"DeltaR;";
6050 histoname = labelname+
"_DeltaPhi";
6051 title = labelname+
"DeltaPhi;";
6055 histoname = labelname+
"_NumeratorAvrgPFPt";
6056 title = labelname+
"NumeratorAvrgPFPt;PF pT [GeV/c] ";
6057 MonitorElement * NumeratorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
6060 histoname = labelname+
"_NumeratorAvrgPFEta";
6061 title = labelname+
"NumeratorAvrgPFEta;PF #eta ";
6062 MonitorElement * NumeratorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
6063 NumeratorPFEta->
getTH1();
6065 histoname = labelname+
"_DenominatorAvrgPFPt";
6066 title = labelname+
"DenominatorAvrgPFPt;PF Pt[GeV/c] ";
6067 MonitorElement * DenominatorPFPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
6068 DenominatorPFPt->
getTH1();
6070 histoname = labelname+
"_DenominatorAvrgPFEta";
6071 title = labelname+
"DenominatorAvrgPFEta;PF #eta ";
6072 MonitorElement * DenominatorPFEta = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
6073 DenominatorPFEta->
getTH1();
6075 histoname = labelname+
"_PFDeltaR";
6076 title = labelname+
"PFDeltaR;";
6080 histoname = labelname+
"_PFDeltaPhi";
6081 title = labelname+
"PFDeltaPhi;";
6085 v->setEffHistos( dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6086 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6087 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6088 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6089 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
6100 histoname = labelname+
"_NumeratorPt";
6101 if(
v->getPath().find(
"HLT_PFMET")==std::string::npos)
6102 title = labelname+
"NumeratorPt;CaloMET[GeV/c]";
6104 title = labelname+
"NumeratorPt;PFMET[GeV/c]";
6105 MonitorElement * NumeratorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
6108 histoname = labelname+
"_NumeratorPhi";
6109 title = labelname+
"NumeratorPhi;#Phi";
6110 MonitorElement * NumeratorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
6113 histoname = labelname+
"_DenominatorPt";
6114 if(
v->getPath().find(
"HLT_PFMET")==std::string::npos)
6115 title = labelname+
"DenominatorPt;CaloMET[GeV/c]";
6117 title = labelname+
"DenominatorPt;PFMET[GeV/c]";
6118 MonitorElement * DenominatorPt = dbe->
book1D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_);
6121 histoname = labelname+
"_DenominatorPhi";
6122 title = labelname+
"DenominatorPhi;#Phi";
6123 MonitorElement * DenominatorPhi = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
6124 DenominatorPhi->
getTH1();
6126 v->setEffHistos( NumeratorPt, dummy, dummy, dummy, dummy, NumeratorPhi, dummy, dummy, dummy, dummy,
6127 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6128 dummy, dummy, DenominatorPt, dummy, dummy, dummy, dummy, DenominatorPhi, dummy, dummy,
6129 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6130 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6131 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy,
6132 dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy, 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
6153 double PtMax_ = 1000.;
6154 double EtaMin_ = -5.;
6155 double EtaMax_ = 5.;
6156 double PhiMin_ = -3.14159;
6157 double PhiMax_ = 3.14159;
6170 histoname = labelname+
"_TriggerSummary";
6171 title = labelname+
"Summary of trigger levels";
6174 std::vector<std::string> trigger;
6175 trigger.push_back(
"Nevt");
6176 trigger.push_back(
"L1 failed");
6177 trigger.push_back(
"L1 & HLT failed");
6178 trigger.push_back(
"L1 failed but not HLT");
6179 trigger.push_back(
"L1 passed");
6180 trigger.push_back(
"L1 & HLT passed");
6181 trigger.push_back(
"L1 passed but not HLT");
6183 for(
unsigned int i =0;
i < trigger.size();
i++)
6186 if((
v->getTriggerType().compare(
"SingleJet_Trigger") == 0)){
6187 histoname = labelname+
"_JetPt";
6188 title = labelname+
"Leading jet pT;Pt[GeV/c]";
6192 histoname = labelname+
"_JetEtaVsPt";
6193 title = labelname+
"Leading jet #eta vs pT;#eta;Pt[GeV/c]";
6194 MonitorElement * JetEtaVsPt = dbe->
book2D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_,Ptbins_,PtMin_,PtMax_);
6197 histoname = labelname+
"_JetPhiVsPt";
6198 title = labelname+
"Leading jet #Phi vs pT;#Phi;Pt[GeV/c]";
6199 MonitorElement * JetPhiVsPt = dbe->
book2D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_,Ptbins_,PtMin_,PtMax_);
6202 v->setDgnsHistos( TriggerSummary, dummy, JetPt, JetEtaVsPt, JetPhiVsPt, dummy, dummy, dummy, dummy, dummy, dummy);
6205 if((
v->getTriggerType().compare(
"DiJet_Trigger") == 0)){
6206 histoname = labelname+
"_JetSize";
6207 title = labelname+
"Jet Size;multiplicity";
6211 histoname = labelname+
"_AvergPt";
6212 title = labelname+
"Average Pt;Pt[GeV/c]";
6216 histoname = labelname+
"_AvergEta";
6217 title = labelname+
"Average Eta;#eta";
6218 MonitorElement * Eta12 = dbe->
book1D(histoname.c_str(),title.c_str(),Etabins_,EtaMin_,EtaMax_);
6221 histoname = labelname+
"_PhiDifference";
6222 title = labelname+
"#Delta#Phi;#Delta#Phi";
6223 MonitorElement * Phi12 = dbe->
book1D(histoname.c_str(),title.c_str(),Phibins_,PhiMin_,PhiMax_);
6226 histoname = labelname+
"_Pt3Jet";
6227 title = labelname+
"Pt of 3rd Jet;Pt[GeV/c]";
6231 histoname = labelname+
"_Pt12VsPt3Jet";
6232 title = labelname+
"Pt of 3rd Jet vs Average Pt of leading jets;Avergage Pt[GeV/c]; Pt of 3rd Jet [GeV/c]";
6233 MonitorElement * Pt12Pt3 = dbe->
book2D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_,Ptbins_,PtMin_,PtMax_);
6236 histoname = labelname+
"_Pt12VsPhi12";
6237 title = labelname+
"Average Pt of leading jets vs #Delta#Phi between leading jets;Avergage Pt[GeV/c]; #Delta#Phi";
6238 MonitorElement * Pt12Phi12 = dbe->
book2D(histoname.c_str(),title.c_str(),Ptbins_,PtMin_,PtMax_,Phibins_,PhiMin_,PhiMax_);
6241 v->setDgnsHistos( TriggerSummary, JetSize, dummy, dummy, dummy, Pt12, Eta12, Phi12, Pt3, Pt12Pt3, Pt12Phi12);
6244 if((
v->getTriggerType().compare(
"MET_Trigger") == 0)){
6245 histoname = labelname+
"_MET";
6246 title = labelname+
"MET;Pt[GeV/c]";
6250 v->setDgnsHistos(TriggerSummary, dummy, MET, dummy, dummy, dummy, dummy, dummy,dummy,dummy,dummy);
6253 if((
v->getTriggerType().compare(
"TET_Trigger") == 0)){
6254 histoname = labelname+
"_TET";
6255 title = labelname+
"TET;Pt[GeV/c]";
6259 v->setDgnsHistos(TriggerSummary, dummy, TET, dummy, dummy, dummy, dummy, dummy,dummy,dummy,dummy);
6287 for(std::vector<std::string>::iterator numpathmodule = numpathmodules.begin();
6288 numpathmodule!= numpathmodules.end(); ++numpathmodule ) {
6291 l1pathname = *numpathmodule;
6312 if (fabs(eta)<=1.3) output=
true;
6319 if (fabs(eta)<=3.0 && fabs(eta)>1.3) output=
true;
6326 if (fabs(eta)>3.0) output=
true;
6355 double binVal = -100;
6356 for(
int ibin=1; ibin<nbins+1; ibin++)
6358 const char * binLabel =
rate_All->
getTH1()->GetXaxis()->GetBinLabel(ibin);
6359 if(binLabel[0]==
'\0')
continue;
6362 if(trigName.compare(binLabel)!=0)
continue;
6364 if(trigName.compare(binLabel)==0){
6379 edm::LogInfo(
"JetMETHLTOfflineSource") <<
"no index "<< index <<
" of that name ";
6382 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
const std::string moduleType(const std::string &module) const
C++ class name of module.
virtual bool isEndCap(double eta)
edm::InputTag caloMETTag_
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)
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)
edm::InputTag caloJetsTag_
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_
~JetMETHLTOfflineSource()
virtual void fillMEforEffWrtMBTrigger(const edm::Event &iEvent, const edm::EventSetup &)
PathInfoCollection hltPathsEffWrtMB_
void setVerbose(unsigned level)
reco::CaloJetCollection calojet
std::vector< std::string > MuonTrigPaths_
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
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)
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_
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
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_
virtual void analyze(const edm::Event &, const edm::EventSetup &)
virtual void fillMEforEffAllTrigger(const edm::Event &iEvent, const edm::EventSetup &)
MonitorElement * rate_AllWrtMB
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_