|
|
|
#include <MuCorrMETAnalyzer.h>
DQM monitoring source for CaloMET
Definition at line 38 of file MuCorrMETAnalyzer.h.
| MuCorrMETAnalyzer::MuCorrMETAnalyzer | ( | const edm::ParameterSet & | pSet | ) |
Constructor.
Definition at line 30 of file MuCorrMETAnalyzer.cc.
References Parameters::parameters.
{
parameters = pSet;
}
| MuCorrMETAnalyzer::~MuCorrMETAnalyzer | ( | ) | [virtual] |
| void MuCorrMETAnalyzer::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup, | ||
| const edm::TriggerResults & | triggerResults | ||
| ) |
Get the analysis.
Definition at line 264 of file MuCorrMETAnalyzer.cc.
References edm::HLTGlobalStatus::accept(), FamosSequences_cff::caloJets, gather_cfg::cout, cmsPerfPublish::DirName, edm::Event::getByLabel(), i, edm::HandleBase::isValid(), bTagSequences_cff::jetID, LogDebug, LogTrace, metname, edm::Handle< T >::product(), edm::HLTGlobalStatus::size(), edm::TriggerNames::size(), AlCaHLTBitMon_QueryRunRegistry::string, edm::TriggerNames::triggerIndex(), and edm::Event::triggerNames().
{
if (_verbose) std::cout << "MuCorrMETAnalyzer analyze" << std::endl;
LogTrace(metname)<<"[MuCorrMETAnalyzer] Analyze MuCorrMET";
metME->Fill(4);
// ==========================================================
// Trigger information
//
_trig_JetMB=0;
_trig_HighPtJet=0;
_trig_LowPtJet=0;
_trig_HighMET=0;
// _trig_LowMET=0;
if (&triggerResults) {
//
//
// Check how many HLT triggers are in triggerResults
int ntrigs = triggerResults.size();
if (_verbose) std::cout << "ntrigs=" << ntrigs << std::endl;
//
//
// If index=ntrigs, this HLT trigger doesn't exist in the HLT table for this data.
const edm::TriggerNames & triggerNames = iEvent.triggerNames(triggerResults);
//
//
// count number of requested Jet or MB HLT paths which have fired
for (unsigned int i=0; i!=HLTPathsJetMBByName_.size(); i++) {
unsigned int triggerIndex = triggerNames.triggerIndex(HLTPathsJetMBByName_[i]);
if (triggerIndex<triggerResults.size()) {
if (triggerResults.accept(triggerIndex)) {
_trig_JetMB++;
}
}
}
// for empty input vectors (n==0), take all HLT triggers!
if (HLTPathsJetMBByName_.size()==0) _trig_JetMB=triggerResults.size()-1;
//
if (_verbose) std::cout << "triggerNames size" << " " << triggerNames.size() << std::endl;
if (_verbose) std::cout << _hlt_HighPtJet << " " << triggerNames.triggerIndex(_hlt_HighPtJet) << std::endl;
if (_verbose) std::cout << _hlt_LowPtJet << " " << triggerNames.triggerIndex(_hlt_LowPtJet) << std::endl;
if (_verbose) std::cout << _hlt_HighMET << " " << triggerNames.triggerIndex(_hlt_HighMET) << std::endl;
// if (_verbose) std::cout << _hlt_LowMET << " " << triggerNames.triggerIndex(_hlt_LowMET) << std::endl;
if (_verbose) std::cout << _hlt_Ele << " " << triggerNames.triggerIndex(_hlt_Ele) << std::endl;
if (_verbose) std::cout << _hlt_Muon << " " << triggerNames.triggerIndex(_hlt_Muon) << std::endl;
if (triggerNames.triggerIndex(_hlt_HighPtJet) != triggerNames.size() &&
triggerResults.accept(triggerNames.triggerIndex(_hlt_HighPtJet))) _trig_HighPtJet=1;
if (triggerNames.triggerIndex(_hlt_LowPtJet) != triggerNames.size() &&
triggerResults.accept(triggerNames.triggerIndex(_hlt_LowPtJet))) _trig_LowPtJet=1;
if (triggerNames.triggerIndex(_hlt_HighMET) != triggerNames.size() &&
triggerResults.accept(triggerNames.triggerIndex(_hlt_HighMET))) _trig_HighMET=1;
// if (triggerNames.triggerIndex(_hlt_LowMET) != triggerNames.size() &&
// triggerResults.accept(triggerNames.triggerIndex(_hlt_LowMET))) _trig_LowMET=1;
if (triggerNames.triggerIndex(_hlt_Ele) != triggerNames.size() &&
triggerResults.accept(triggerNames.triggerIndex(_hlt_Ele))) _trig_Ele=1;
if (triggerNames.triggerIndex(_hlt_Muon) != triggerNames.size() &&
triggerResults.accept(triggerNames.triggerIndex(_hlt_Muon))) _trig_Muon=1;
} else {
edm::LogInfo("MuCorrMetAnalyzer") << "TriggerResults::HLT not found, "
"automatically select events";
// TriggerResults object not found. Look at all events.
_trig_JetMB=1;
}
// ==========================================================
// MuCorrMET information
// **** Get the MET container
edm::Handle<reco::CaloMETCollection> muCorrmetcoll;
iEvent.getByLabel("corMetGlobalMuons", muCorrmetcoll);
if(!muCorrmetcoll.isValid()) return;
const CaloMETCollection *muCorrmetcol = muCorrmetcoll.product();
const CaloMET *muCorrmet;
muCorrmet = &(muCorrmetcol->front());
LogTrace(metname)<<"[MuCorrMETAnalyzer] Call to the MuCorrMET analyzer";
// ==========================================================
//
edm::Handle<HcalNoiseRBXCollection> HRBXCollection;
iEvent.getByLabel(HcalNoiseRBXCollectionTag,HRBXCollection);
if (!HRBXCollection.isValid()) {
LogDebug("") << "MuCorrMETAnalyzer: Could not find HcalNoiseRBX Collection" << std::endl;
if (_verbose) std::cout << "MuCorrMETAnalyzer: Could not find HcalNoiseRBX Collection" << std::endl;
}
edm::Handle<bool> HBHENoiseFilterResultHandle;
iEvent.getByLabel(HBHENoiseFilterResultTag, HBHENoiseFilterResultHandle);
bool HBHENoiseFilterResult = *HBHENoiseFilterResultHandle;
if (!HBHENoiseFilterResultHandle.isValid()) {
LogDebug("") << "MuCorrMETAnalyzer: Could not find HBHENoiseFilterResult" << std::endl;
if (_verbose) std::cout << "MuCorrMETAnalyzer: Could not find HBHENoiseFilterResult" << std::endl;
}
edm::Handle<reco::CaloJetCollection> caloJets;
iEvent.getByLabel(theJetCollectionLabel, caloJets);
if (!caloJets.isValid()) {
LogDebug("") << "MuCorrMETAnalyzer: Could not find jet product" << std::endl;
if (_verbose) std::cout << "MuCorrMETAnalyzer: Could not find jet product" << std::endl;
}
// ==========================================================
// MuCorrMET sanity check
// if (_source=="MuCorrMET") validateMET(*muCorrmet, tcCandidates);
// ==========================================================
// JetID
if (_verbose) std::cout << "JetID starts" << std::endl;
//
// --- Loose cuts, not Muon specific for now!
//
bool bJetID=true;
for (reco::CaloJetCollection::const_iterator cal = caloJets->begin();
cal!=caloJets->end(); ++cal){
jetID->calculate(iEvent, *cal);
if (_verbose) std::cout << jetID->n90Hits() << " "
<< jetID->restrictedEMF() << " "
<< cal->pt() << std::endl;
if (cal->pt()>10.){
//
// for all regions
if (jetID->n90Hits()<2) bJetID=false;
if (jetID->fHPD()>=0.98) bJetID=false;
//if (jetID->restrictedEMF()<0.01) bJetID=false;
//
// for non-forward
if (fabs(cal->eta())<2.55){
if (cal->emEnergyFraction()<=0.01) bJetID=false;
}
// for forward
else {
if (cal->emEnergyFraction()<=-0.9) bJetID=false;
if (cal->pt()>80.){
if (cal->emEnergyFraction()>= 1.0) bJetID=false;
}
} // forward vs non-forward
} // pt>10 GeV/c
} // calor-jets loop
//
// --- Tight cuts
//
bool bJetIDTight=true;
bJetIDTight=bJetID;
for (reco::CaloJetCollection::const_iterator cal = caloJets->begin();
cal!=caloJets->end(); ++cal){
jetID->calculate(iEvent, *cal);
if (cal->pt()>25.){
//
// for all regions
if (jetID->fHPD()>=0.95) bJetIDTight=false;
//
// for 1.0<|eta|<1.75
if (fabs(cal->eta())>=1.00 && fabs(cal->eta())<1.75){
if (cal->pt()>80. && cal->emEnergyFraction()>=1.) bJetIDTight=false;
}
//
// for 1.75<|eta|<2.55
else if (fabs(cal->eta())>=1.75 && fabs(cal->eta())<2.55){
if (cal->pt()>80. && cal->emEnergyFraction()>=1.) bJetIDTight=false;
}
//
// for 2.55<|eta|<3.25
else if (fabs(cal->eta())>=2.55 && fabs(cal->eta())<3.25){
if (cal->pt()< 50. && cal->emEnergyFraction()<=-0.3) bJetIDTight=false;
if (cal->pt()>=50. && cal->pt()< 80. && cal->emEnergyFraction()<=-0.2) bJetIDTight=false;
if (cal->pt()>=80. && cal->pt()<340. && cal->emEnergyFraction()<=-0.1) bJetIDTight=false;
if (cal->pt()>=340. && cal->emEnergyFraction()<=-0.1
&& cal->emEnergyFraction()>=0.95) bJetIDTight=false;
}
//
// for 3.25<|eta|
else if (fabs(cal->eta())>=3.25){
if (cal->pt()< 50. && cal->emEnergyFraction()<=-0.3
&& cal->emEnergyFraction()>=0.90) bJetIDTight=false;
if (cal->pt()>=50. && cal->pt()<130. && cal->emEnergyFraction()<=-0.2
&& cal->emEnergyFraction()>=0.80) bJetIDTight=false;
if (cal->pt()>=130. && cal->emEnergyFraction()<=-0.1
&& cal->emEnergyFraction()>=0.70) bJetIDTight=false;
}
} // pt>10 GeV/c
} // calor-jets loop
if (_verbose) std::cout << "JetID ends" << std::endl;
// ==========================================================
// HCAL Noise filter
bool bHcalNoiseFilter = HBHENoiseFilterResult;
// ==========================================================
// Reconstructed MET Information - fill MonitorElements
std::string DirName = "JetMET/MET/"+_source;
for (std::vector<std::string>::const_iterator ic = _FolderNames.begin();
ic != _FolderNames.end(); ic++){
if (*ic=="All") fillMESet(iEvent, DirName+"/"+*ic, *muCorrmet);
if (*ic=="Cleanup" && bHcalNoiseFilter && bJetID) fillMESet(iEvent, DirName+"/"+*ic, *muCorrmet);
if (_allSelection) {
if (*ic=="HcalNoiseFilter" && bHcalNoiseFilter ) fillMESet(iEvent, DirName+"/"+*ic, *muCorrmet);
if (*ic=="JetID" && bJetID) fillMESet(iEvent, DirName+"/"+*ic, *muCorrmet);
if (*ic=="JetIDTight" && bJetIDTight) fillMESet(iEvent, DirName+"/"+*ic, *muCorrmet);
}
}
}
| void MuCorrMETAnalyzer::beginJob | ( | DQMStore * | dbe | ) | [virtual] |
Inizialize parameters for histo binning.
Implements MuCorrMETAnalyzerBase.
Definition at line 39 of file MuCorrMETAnalyzer.cc.
References DQMStore::book1D(), cmsPerfPublish::DirName, bTagSequences_cff::jetID, LogTrace, metname, Parameters::parameters, MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.
{
evtCounter = 0;
metname = "muonMETAnalyzer";
// trigger information
HLTPathsJetMBByName_ = parameters.getParameter<std::vector<std::string > >("HLTPathsJetMB");
_hlt_HighPtJet = parameters.getParameter<std::string>("HLT_HighPtJet");
_hlt_LowPtJet = parameters.getParameter<std::string>("HLT_LowPtJet");
_hlt_HighMET = parameters.getParameter<std::string>("HLT_HighMET");
// _hlt_LowMET = parameters.getParameter<std::string>("HLT_LowMET");
_hlt_Ele = parameters.getParameter<std::string>("HLT_Ele");
_hlt_Muon = parameters.getParameter<std::string>("HLT_Muon");
// MuCorrMET information
theMuCorrMETCollectionLabel = parameters.getParameter<edm::InputTag>("MuCorrMETCollectionLabel");
_source = parameters.getParameter<std::string>("Source");
// Other data collections
HcalNoiseRBXCollectionTag = parameters.getParameter<edm::InputTag>("HcalNoiseRBXCollection");
theJetCollectionLabel = parameters.getParameter<edm::InputTag>("JetCollectionLabel");
HBHENoiseFilterResultTag = parameters.getParameter<edm::InputTag>("HBHENoiseFilterResultLabel");
// misc
_verbose = parameters.getParameter<int>("verbose");
_etThreshold = parameters.getParameter<double>("etThreshold"); // MET threshold
_allhist = parameters.getParameter<bool>("allHist"); // Full set of monitoring histograms
_allSelection= parameters.getParameter<bool>("allSelection"); // Plot with all sets of event selection
_highPtMuCorrJetThreshold = parameters.getParameter<double>("HighPtMuCorrJetThreshold"); // High Pt Jet threshold
_lowPtMuCorrJetThreshold = parameters.getParameter<double>("LowPtMuCorrJetThreshold"); // Low Pt Jet threshold
_highMuCorrMETThreshold = parameters.getParameter<double>("HighMuCorrMETThreshold"); // High MET threshold
// _lowMuCorrMETThreshold = parameters.getParameter<double>("LowMuCorrMETThreshold"); // Low MET threshold
//
jetID = new reco::helper::JetIDHelper(parameters.getParameter<ParameterSet>("JetIDParams"));
// DQStore stuff
LogTrace(metname)<<"[MuCorrMETAnalyzer] Parameters initialization";
std::string DirName = "JetMET/MET/"+_source;
dbe->setCurrentFolder(DirName);
metME = dbe->book1D("metReco", "metReco", 4, 1, 5);
metME->setBinLabel(4,"MuCorrMET",1);
_dbe = dbe;
_FolderNames.push_back("All");
_FolderNames.push_back("Cleanup");
_FolderNames.push_back("HcalNoiseFilter");
_FolderNames.push_back("JetID");
_FolderNames.push_back("JetIDTight");
for (std::vector<std::string>::const_iterator ic = _FolderNames.begin();
ic != _FolderNames.end(); ic++){
if (*ic=="All") bookMESet(DirName+"/"+*ic);
if (*ic=="Cleanup") bookMESet(DirName+"/"+*ic);
if (_allSelection){
if (*ic=="HcalNoiseFilter") bookMESet(DirName+"/"+*ic);
if (*ic=="JetID") bookMESet(DirName+"/"+*ic);
if (*ic=="JetIDTight") bookMESet(DirName+"/"+*ic);
}
}
}
| void MuCorrMETAnalyzer::beginRun | ( | const edm::Run & | iRun, |
| const edm::EventSetup & | iSetup | ||
| ) |
| void MuCorrMETAnalyzer::bookMESet | ( | std::string | DirName | ) |
Definition at line 113 of file MuCorrMETAnalyzer.cc.
{
bool bLumiSecPlot=false;
if (DirName.find("All")!=std::string::npos) bLumiSecPlot=true;
bookMonitorElement(DirName,bLumiSecPlot);
if (_hlt_HighPtJet.size()){
bookMonitorElement(DirName+"/"+"HighPtJet",false);
meTriggerName_HighPtJet = _dbe->bookString("triggerName_HighPtJet", _hlt_HighPtJet);
}
if (_hlt_LowPtJet.size()){
bookMonitorElement(DirName+"/"+"LowPtJet",false);
meTriggerName_LowPtJet = _dbe->bookString("triggerName_LowPtJet", _hlt_LowPtJet);
}
if (_hlt_HighMET.size()){
bookMonitorElement(DirName+"/"+"HighMET",false);
meTriggerName_HighMET = _dbe->bookString("triggerName_HighMET", _hlt_HighMET);
}
// if (_hlt_LowMET.size()){
// bookMonitorElement(DirName+"/"+"LowMET",false);
// meTriggerName_LowMET = _dbe->bookString("triggerName_LowMET", _hlt_LowMET);
// }
if (_hlt_Ele.size()){
bookMonitorElement(DirName+"/"+"Ele",false);
meTriggerName_Ele = _dbe->bookString("triggerName_Ele", _hlt_Ele);
}
if (_hlt_Muon.size()){
bookMonitorElement(DirName+"/"+"Muon",false);
meTriggerName_Muon = _dbe->bookString("triggerName_Muon", _hlt_Muon);
}
}
| void MuCorrMETAnalyzer::bookMonitorElement | ( | std::string | DirName, |
| bool | bLumiSecPlot = false |
||
| ) |
Definition at line 155 of file MuCorrMETAnalyzer.cc.
References gather_cfg::cout.
{
if (_verbose) std::cout << "booMonitorElement " << DirName << std::endl;
_dbe->setCurrentFolder(DirName);
meMuCorrMEx = _dbe->book1D("METTask_MuCorrMEx", "METTask_MuCorrMEx" , 200, -500, 500);
meMuCorrMEy = _dbe->book1D("METTask_MuCorrMEy", "METTask_MuCorrMEy" , 200, -500, 500);
meMuCorrMET = _dbe->book1D("METTask_MuCorrMET", "METTask_MuCorrMET" , 200, 0, 1000);
meMuCorrSumET = _dbe->book1D("METTask_MuCorrSumET", "METTask_MuCorrSumET" , 400, 0, 4000);
meMuCorrMETSig = _dbe->book1D("METTask_MuCorrMETSig", "METTask_MuCorrMETSig", 51, 0, 51);
meMuCorrMETPhi = _dbe->book1D("METTask_MuCorrMETPhi", "METTask_MuCorrMETPhi", 60, -3.2, 3.2);
meMuCorrMETIonFeedbck = _dbe->book1D("METTask_MuCorrMETIonFeedbck", "METTask_MuCorrMETIonFeedbck", 200, 0, 1000);
meMuCorrMETHPDNoise = _dbe->book1D("METTask_MuCorrMETHPDNoise", "METTask_MuCorrMETHPDNoise", 200, 0, 1000);
meMuCorrMETRBXNoise = _dbe->book1D("METTask_MuCorrMETRBXNoise", "METTask_MuCorrMETRBXNoise", 200, 0, 1000);
if (_allhist) {
if (bLumiSecPlot) {
meMuCorrMExLS = _dbe->book2D("METTask_MuCorrMEx_LS","METTask_MuCorrMEx_LS", 200, -200, 200, 50, 0, 500);
meMuCorrMEyLS = _dbe->book2D("METTask_MuCorrMEy_LS","METTask_MuCorrMEy_LS", 200, -200, 200, 50, 0, 500);
}
}
}
| void MuCorrMETAnalyzer::endJob | ( | void | ) |
Finish up a job.
Definition at line 106 of file MuCorrMETAnalyzer.cc.
References bTagSequences_cff::jetID.
{
delete jetID;
}
| void MuCorrMETAnalyzer::endRun | ( | const edm::Run & | iRun, |
| const edm::EventSetup & | iSetup, | ||
| DQMStore * | dbe | ||
| ) |
Finish up a run.
Definition at line 188 of file MuCorrMETAnalyzer.cc.
References cmsPerfPublish::DirName, TrackerOfflineValidation_Dqm_cff::dirName, MonitorElement::getRootObject(), MonitorElement::getTH1F(), i, and AlCaHLTBitMon_QueryRunRegistry::string.
{
//
//--- Check the time length of the Run from the lumi section plots
std::string dirName = "JetMET/MET/"+_source+"/";
_dbe->setCurrentFolder(dirName);
TH1F* tlumisec;
MonitorElement *meLumiSec = _dbe->get("aaa");
meLumiSec = _dbe->get("JetMET/lumisec");
int totlsec=0;
double totltime=0.;
if ( meLumiSec->getRootObject() ) {
tlumisec = meLumiSec->getTH1F();
for (int i=0; i<500; i++){
if (tlumisec->GetBinContent(i+1)) totlsec++;
}
totltime = double(totlsec*90); // one lumi sec ~ 90 (sec)
}
if (totltime==0.) totltime=1.;
//
//--- Make the integrated plots with rate (Hz)
for (std::vector<std::string>::const_iterator ic = _FolderNames.begin(); ic != _FolderNames.end(); ic++)
{
std::string DirName;
DirName = dirName+*ic;
makeRatePlot(DirName,totltime);
if (_hlt_HighPtJet.size()) makeRatePlot(DirName+"/"+_hlt_HighPtJet,totltime);
if (_hlt_LowPtJet.size()) makeRatePlot(DirName+"/"+_hlt_LowPtJet,totltime);
if (_hlt_HighMET.size()) makeRatePlot(DirName+"/"+_hlt_HighMET,totltime);
// if (_hlt_LowMET.size()) makeRatePlot(DirName+"/"+_hlt_LowMET,totltime);
if (_hlt_Ele.size()) makeRatePlot(DirName+"/"+_hlt_Ele,totltime);
if (_hlt_Muon.size()) makeRatePlot(DirName+"/"+_hlt_Muon,totltime);
}
}
| void MuCorrMETAnalyzer::fillMESet | ( | const edm::Event & | iEvent, |
| std::string | DirName, | ||
| const reco::CaloMET & | muCorrmet | ||
| ) |
Definition at line 499 of file MuCorrMETAnalyzer.cc.
{
_dbe->setCurrentFolder(DirName);
bool bLumiSecPlot=false;
if (DirName.find("All")) bLumiSecPlot=true;
if (_trig_JetMB) fillMonitorElement(iEvent,DirName,"",muCorrmet, bLumiSecPlot);
if (_hlt_HighPtJet.size() && _trig_HighPtJet) fillMonitorElement(iEvent,DirName,"HighPtJet",muCorrmet,false);
if (_hlt_LowPtJet.size() && _trig_LowPtJet) fillMonitorElement(iEvent,DirName,"LowPtJet",muCorrmet,false);
if (_hlt_HighMET.size() && _trig_HighMET) fillMonitorElement(iEvent,DirName,"HighMET",muCorrmet,false);
// if (_hlt_LowMET.size() && _trig_LowMET) fillMonitorElement(iEvent,DirName,"LowMET",muCorrmet,false);
if (_hlt_Ele.size() && _trig_Ele) fillMonitorElement(iEvent,DirName,"Ele",muCorrmet,false);
if (_hlt_Muon.size() && _trig_Muon) fillMonitorElement(iEvent,DirName,"Muon",muCorrmet,false);
}
| void MuCorrMETAnalyzer::fillMonitorElement | ( | const edm::Event & | iEvent, |
| std::string | DirName, | ||
| std::string | TriggerTypeName, | ||
| const reco::CaloMET & | muCorrmet, | ||
| bool | bLumiSecPlot | ||
| ) |
Definition at line 518 of file MuCorrMETAnalyzer.cc.
References gather_cfg::cout, edm::EventBase::luminosityBlock(), reco::MET::mEtSig(), reco::LeafCandidate::phi(), reco::LeafCandidate::pt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), and reco::MET::sumEt().
{
if (TriggerTypeName=="HighPtJet") {
if (!selectHighPtJetEvent(iEvent)) return;
}
else if (TriggerTypeName=="LowPtJet") {
if (!selectLowPtJetEvent(iEvent)) return;
}
else if (TriggerTypeName=="HighMET") {
if (muCorrmet.pt()<_highMuCorrMETThreshold) return;
}
// else if (TriggerTypeName=="LowMET") {
// if (muCorrmet.pt()<_lowMuCorrMETThreshold) return;
// }
else if (TriggerTypeName=="Ele") {
if (!selectWElectronEvent(iEvent)) return;
}
else if (TriggerTypeName=="Muon") {
if (!selectWMuonEvent(iEvent)) return;
}
// Reconstructed MET Information
double mucorrSumET = muCorrmet.sumEt();
double mucorrmetSig = muCorrmet.mEtSig();
// double mucorrEz = muCorrmet.e_longitudinal();
double mucorrmet = muCorrmet.pt();
double mucorrMEx = muCorrmet.px();
double mucorrMEy = muCorrmet.py();
double mucorrmetPhi = muCorrmet.phi();
//
int myLuminosityBlock;
// myLuminosityBlock = (evtCounter++)/1000;
myLuminosityBlock = iEvent.luminosityBlock();
//
if (TriggerTypeName!="") DirName = DirName +"/"+TriggerTypeName;
if (_verbose) std::cout << "_etThreshold = " << _etThreshold << std::endl;
if (mucorrmet>_etThreshold){
meMuCorrMEx = _dbe->get(DirName+"/"+"METTask_MuCorrMEx"); if (meMuCorrMEx && meMuCorrMEx->getRootObject()) meMuCorrMEx->Fill(mucorrMEx);
meMuCorrMEy = _dbe->get(DirName+"/"+"METTask_MuCorrMEy"); if (meMuCorrMEy && meMuCorrMEy->getRootObject()) meMuCorrMEy->Fill(mucorrMEy);
meMuCorrMET = _dbe->get(DirName+"/"+"METTask_MuCorrMET"); if (meMuCorrMET && meMuCorrMET->getRootObject()) meMuCorrMET->Fill(mucorrmet);
meMuCorrMETPhi = _dbe->get(DirName+"/"+"METTask_MuCorrMETPhi"); if (meMuCorrMETPhi && meMuCorrMETPhi->getRootObject()) meMuCorrMETPhi->Fill(mucorrmetPhi);
meMuCorrSumET = _dbe->get(DirName+"/"+"METTask_MuCorrSumET"); if (meMuCorrSumET && meMuCorrSumET->getRootObject()) meMuCorrSumET->Fill(mucorrSumET);
meMuCorrMETSig = _dbe->get(DirName+"/"+"METTask_MuCorrMETSig"); if (meMuCorrMETSig && meMuCorrMETSig->getRootObject()) meMuCorrMETSig->Fill(mucorrmetSig);
meMuCorrMETIonFeedbck = _dbe->get(DirName+"/"+"METTask_MuCorrMETIonFeedbck"); if (meMuCorrMETIonFeedbck && meMuCorrMETIonFeedbck->getRootObject()) meMuCorrMETIonFeedbck->Fill(mucorrmet);
meMuCorrMETHPDNoise = _dbe->get(DirName+"/"+"METTask_MuCorrMETHPDNoise"); if (meMuCorrMETHPDNoise && meMuCorrMETHPDNoise->getRootObject()) meMuCorrMETHPDNoise->Fill(mucorrmet);
meMuCorrMETRBXNoise = _dbe->get(DirName+"/"+"METTask_MuCorrMETRBXNoise"); if (meMuCorrMETRBXNoise && meMuCorrMETRBXNoise->getRootObject()) meMuCorrMETRBXNoise->Fill(mucorrmet);
if (_allhist){
if (bLumiSecPlot){
meMuCorrMExLS = _dbe->get(DirName+"/"+"METTask_MuCorrMExLS"); if (meMuCorrMExLS && meMuCorrMExLS->getRootObject()) meMuCorrMExLS->Fill(mucorrMEx,myLuminosityBlock);
meMuCorrMEyLS = _dbe->get(DirName+"/"+"METTask_MuCorrMEyLS"); if (meMuCorrMEyLS && meMuCorrMEyLS->getRootObject()) meMuCorrMEyLS->Fill(mucorrMEy,myLuminosityBlock);
}
} // _allhist
} // et threshold cut
}
| void MuCorrMETAnalyzer::makeRatePlot | ( | std::string | DirName, |
| double | totltime | ||
| ) |
Definition at line 235 of file MuCorrMETAnalyzer.cc.
References MonitorElement::getRootObject(), MonitorElement::getTH1F(), and i.
{
_dbe->setCurrentFolder(DirName);
MonitorElement *meMuCorrMET = _dbe->get(DirName+"/"+"METTask_MuCorrMET");
TH1F* tMuCorrMET;
TH1F* tMuCorrMETRate;
if ( meMuCorrMET )
if ( meMuCorrMET->getRootObject() ) {
tMuCorrMET = meMuCorrMET->getTH1F();
// Integral plot & convert number of events to rate (hz)
tMuCorrMETRate = (TH1F*) tMuCorrMET->Clone("METTask_MuCorrMETRate");
for (int i = tMuCorrMETRate->GetNbinsX()-1; i>=0; i--){
tMuCorrMETRate->SetBinContent(i+1,tMuCorrMETRate->GetBinContent(i+2)+tMuCorrMET->GetBinContent(i+1));
}
for (int i = 0; i<tMuCorrMETRate->GetNbinsX(); i++){
tMuCorrMETRate->SetBinContent(i+1,tMuCorrMETRate->GetBinContent(i+1)/double(totltime));
}
meMuCorrMETRate = _dbe->book1D("METTask_MuCorrMETRate",tMuCorrMETRate);
}
}
| bool MuCorrMETAnalyzer::selectHighPtJetEvent | ( | const edm::Event & | iEvent | ) |
Definition at line 584 of file MuCorrMETAnalyzer.cc.
References FamosSequences_cff::caloJets, gather_cfg::cout, edm::Event::getByLabel(), edm::HandleBase::isValid(), and LogDebug.
{
bool return_value=false;
edm::Handle<reco::CaloJetCollection> caloJets;
iEvent.getByLabel(theJetCollectionLabel, caloJets);
if (!caloJets.isValid()) {
LogDebug("") << "MuCorrMETAnalyzer: Could not find jet product" << std::endl;
if (_verbose) std::cout << "MuCorrMETAnalyzer: Could not find jet product" << std::endl;
}
for (reco::CaloJetCollection::const_iterator cal = caloJets->begin();
cal!=caloJets->end(); ++cal){
if (cal->pt()>_highPtMuCorrJetThreshold){
return_value=true;
}
}
return return_value;
}
| bool MuCorrMETAnalyzer::selectLowPtJetEvent | ( | const edm::Event & | iEvent | ) |
Definition at line 606 of file MuCorrMETAnalyzer.cc.
References FamosSequences_cff::caloJets, gather_cfg::cout, edm::Event::getByLabel(), edm::HandleBase::isValid(), and LogDebug.
{
bool return_value=false;
edm::Handle<reco::CaloJetCollection> caloJets;
iEvent.getByLabel(theJetCollectionLabel, caloJets);
if (!caloJets.isValid()) {
LogDebug("") << "MuCorrMETAnalyzer: Could not find jet product" << std::endl;
if (_verbose) std::cout << "MuCorrMETAnalyzer: Could not find jet product" << std::endl;
}
for (reco::CaloJetCollection::const_iterator cal = caloJets->begin();
cal!=caloJets->end(); ++cal){
if (cal->pt()>_lowPtMuCorrJetThreshold){
return_value=true;
}
}
return return_value;
}
| bool MuCorrMETAnalyzer::selectWElectronEvent | ( | const edm::Event & | iEvent | ) |
Definition at line 629 of file MuCorrMETAnalyzer.cc.
{
bool return_value=false;
/*
W-electron event selection comes here
*/
return return_value;
}
| bool MuCorrMETAnalyzer::selectWMuonEvent | ( | const edm::Event & | iEvent | ) |
Definition at line 642 of file MuCorrMETAnalyzer.cc.
{
bool return_value=false;
/*
W-muon event selection comes here
*/
return return_value;
}
| void MuCorrMETAnalyzer::setSource | ( | std::string | source | ) | [inline] |
Definition at line 77 of file MuCorrMETAnalyzer.h.
References _source, and LaserTracksInput_cfi::source.
bool MuCorrMETAnalyzer::_allhist [private] |
Definition at line 133 of file MuCorrMETAnalyzer.h.
bool MuCorrMETAnalyzer::_allSelection [private] |
Definition at line 134 of file MuCorrMETAnalyzer.h.
DQMStore* MuCorrMETAnalyzer::_dbe [private] |
Definition at line 140 of file MuCorrMETAnalyzer.h.
double MuCorrMETAnalyzer::_etThreshold [private] |
Definition at line 124 of file MuCorrMETAnalyzer.h.
std::vector<std::string> MuCorrMETAnalyzer::_FolderNames [private] |
Definition at line 137 of file MuCorrMETAnalyzer.h.
double MuCorrMETAnalyzer::_highMuCorrMETThreshold [private] |
Definition at line 120 of file MuCorrMETAnalyzer.h.
double MuCorrMETAnalyzer::_highPtMuCorrJetThreshold [private] |
Definition at line 118 of file MuCorrMETAnalyzer.h.
std::string MuCorrMETAnalyzer::_hlt_Ele [private] |
Definition at line 106 of file MuCorrMETAnalyzer.h.
std::string MuCorrMETAnalyzer::_hlt_HighMET [private] |
Definition at line 104 of file MuCorrMETAnalyzer.h.
std::string MuCorrMETAnalyzer::_hlt_HighPtJet [private] |
Definition at line 102 of file MuCorrMETAnalyzer.h.
std::string MuCorrMETAnalyzer::_hlt_LowPtJet [private] |
Definition at line 103 of file MuCorrMETAnalyzer.h.
std::string MuCorrMETAnalyzer::_hlt_Muon [private] |
Definition at line 107 of file MuCorrMETAnalyzer.h.
double MuCorrMETAnalyzer::_lowPtMuCorrJetThreshold [private] |
Definition at line 119 of file MuCorrMETAnalyzer.h.
std::string MuCorrMETAnalyzer::_source [private] |
Definition at line 91 of file MuCorrMETAnalyzer.h.
Referenced by setSource().
int MuCorrMETAnalyzer::_trig_Ele [private] |
Definition at line 114 of file MuCorrMETAnalyzer.h.
int MuCorrMETAnalyzer::_trig_HighMET [private] |
Definition at line 112 of file MuCorrMETAnalyzer.h.
int MuCorrMETAnalyzer::_trig_HighPtJet [private] |
Definition at line 110 of file MuCorrMETAnalyzer.h.
int MuCorrMETAnalyzer::_trig_JetMB [private] |
Definition at line 109 of file MuCorrMETAnalyzer.h.
int MuCorrMETAnalyzer::_trig_LowPtJet [private] |
Definition at line 111 of file MuCorrMETAnalyzer.h.
int MuCorrMETAnalyzer::_trig_Muon [private] |
Definition at line 115 of file MuCorrMETAnalyzer.h.
int MuCorrMETAnalyzer::_verbose [private] |
Definition at line 88 of file MuCorrMETAnalyzer.h.
Definition at line 81 of file MuCorrMETAnalyzer.h.
Definition at line 95 of file MuCorrMETAnalyzer.h.
Definition at line 94 of file MuCorrMETAnalyzer.h.
double MuCorrMETAnalyzer::hfCalibFactor_ [private] |
Definition at line 127 of file MuCorrMETAnalyzer.h.
std::vector<std::string > MuCorrMETAnalyzer::HLTPathsJetMBByName_ [private] |
Definition at line 100 of file MuCorrMETAnalyzer.h.
Definition at line 130 of file MuCorrMETAnalyzer.h.
MonitorElement* MuCorrMETAnalyzer::meMuCorrMET [private] |
Definition at line 155 of file MuCorrMETAnalyzer.h.
Definition at line 162 of file MuCorrMETAnalyzer.h.
Definition at line 161 of file MuCorrMETAnalyzer.h.
MonitorElement* MuCorrMETAnalyzer::meMuCorrMETPhi [private] |
Definition at line 156 of file MuCorrMETAnalyzer.h.
MonitorElement* MuCorrMETAnalyzer::meMuCorrMETRate [private] |
Definition at line 165 of file MuCorrMETAnalyzer.h.
Definition at line 163 of file MuCorrMETAnalyzer.h.
MonitorElement* MuCorrMETAnalyzer::meMuCorrMETSig [private] |
Definition at line 154 of file MuCorrMETAnalyzer.h.
MonitorElement* MuCorrMETAnalyzer::meMuCorrMEx [private] |
Definition at line 152 of file MuCorrMETAnalyzer.h.
MonitorElement* MuCorrMETAnalyzer::meMuCorrMExLS [private] |
Definition at line 158 of file MuCorrMETAnalyzer.h.
MonitorElement* MuCorrMETAnalyzer::meMuCorrMEy [private] |
Definition at line 153 of file MuCorrMETAnalyzer.h.
MonitorElement* MuCorrMETAnalyzer::meMuCorrMEyLS [private] |
Definition at line 159 of file MuCorrMETAnalyzer.h.
MonitorElement* MuCorrMETAnalyzer::meMuCorrSumET [private] |
Definition at line 157 of file MuCorrMETAnalyzer.h.
MonitorElement* MuCorrMETAnalyzer::metME [private] |
Definition at line 143 of file MuCorrMETAnalyzer.h.
std::string MuCorrMETAnalyzer::metname [private] |
Definition at line 90 of file MuCorrMETAnalyzer.h.
Definition at line 149 of file MuCorrMETAnalyzer.h.
Definition at line 147 of file MuCorrMETAnalyzer.h.
Definition at line 145 of file MuCorrMETAnalyzer.h.
Definition at line 146 of file MuCorrMETAnalyzer.h.
Definition at line 150 of file MuCorrMETAnalyzer.h.
Definition at line 86 of file MuCorrMETAnalyzer.h.
Definition at line 96 of file MuCorrMETAnalyzer.h.
Definition at line 93 of file MuCorrMETAnalyzer.h.
Definition at line 97 of file MuCorrMETAnalyzer.h.
1.7.3