|
|
|
Definition at line 58 of file DQMHcalIsoTrackHLT.cc.
| DQMHcalIsoTrackHLT::DQMHcalIsoTrackHLT | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 135 of file DQMHcalIsoTrackHLT.cc.
References folderName_, edm::ParameterSet::getParameter(), hltAODEventTag_, hltProcess_, hltRAWEventTag_, l1filterLabel_, l2collectionLabel_, l2filterLabel_, l3collectionLabel_, l3filterLabel_, outRootFileName_, saveToRootFile_, and useHLTDebug_.
{
folderName_ = iConfig.getParameter<std::string>("folderName");
outRootFileName_=iConfig.getParameter<std::string>("outputRootFileName");
useHLTDebug_=iConfig.getParameter<bool>("useHLTDebug");
hltRAWEventTag_=iConfig.getParameter<std::string>("hltRAWTriggerEventLabel");
hltAODEventTag_=iConfig.getParameter<std::string>("hltAODTriggerEventLabel");
l2collectionLabel_=iConfig.getParameter<std::string>("l2collectionLabel");
l3collectionLabel_=iConfig.getParameter<std::string>("l3collectionLabel");
l3filterLabel_=iConfig.getParameter<std::string>("hltL3filterLabel");
l1filterLabel_=iConfig.getParameter<std::string>("hltL1filterLabel");
l2filterLabel_=iConfig.getParameter<std::string>("hltL2filterLabel");
hltProcess_=iConfig.getParameter<std::string>("hltProcessName");
saveToRootFile_=iConfig.getParameter<bool>("SaveToRootFile");
}
| DQMHcalIsoTrackHLT::~DQMHcalIsoTrackHLT | ( | ) |
Definition at line 157 of file DQMHcalIsoTrackHLT.cc.
{}
| void DQMHcalIsoTrackHLT::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [private, virtual] |
Implements edm::EDAnalyzer.
Definition at line 160 of file DQMHcalIsoTrackHLT.cc.
References eta(), trigger::TriggerObject::eta(), MonitorElement::Fill(), first, edm::Event::getByLabel(), getDist(), hisopT, hisopTvsEta, hL1pT, hL2eta, hL2L3acc, hL2phi, hL2pT, hL2TowerOccupancy, hL3candL2rat, hL3colEta, hL3colP, hL3eta, hL3etaAOD, hL3etaPureAOD, hL3L2rat, hL3L2trackMatch, hL3pAOD, hL3phi, hL3pPureAOD, hL3Pt, hltAODEventTag_, hltProcess_, hltRAWEventTag_, i, edm::HandleBase::isValid(), j, gen::k, l1filterLabel_, patRefSel_triggerSelection_cff::l1Tag, l2filterLabel_, l3collectionLabel_, l3filterLabel_, AlCaHLTBitMon_ParallelJobs::p, phi, trigger::TriggerObject::pt(), edm::second(), findQualityFiles::size, trigger::TriggerL1CenJet, trigger::TriggerL1ForJet, trigger::TriggerL1TauJet, trigger::TriggerTrack, and useHLTDebug_.
{
edm::Handle<trigger::TriggerEventWithRefs> triggerObj;
edm::InputTag toLab=edm::InputTag(hltRAWEventTag_,"",hltProcess_);
iEvent.getByLabel(toLab,triggerObj);
if(!triggerObj.isValid())
{
edm::LogWarning("DQMHcalIsoTrack") << "RAW-type HLT results not found, skipping event";
return;
}
std::vector<l1extra::L1JetParticleRef> l1CenJets;
std::vector<l1extra::L1JetParticleRef> l1ForJets;
std::vector<l1extra::L1JetParticleRef> l1TauJets;
edm::InputTag l1Tag = edm::InputTag(l1filterLabel_, "",hltProcess_);
trigger::size_type l1filterIndex=triggerObj->filterIndex(l1Tag);
if (l1filterIndex<triggerObj->size())
{
triggerObj->getObjects(l1filterIndex, trigger::TriggerL1CenJet, l1CenJets);
triggerObj->getObjects(l1filterIndex, trigger::TriggerL1ForJet, l1ForJets);
triggerObj->getObjects(l1filterIndex, trigger::TriggerL1TauJet, l1TauJets);
}
std::vector<reco::IsolatedPixelTrackCandidateRef> l2tracks;
edm::InputTag l2Tag = edm::InputTag(l2filterLabel_,"",hltProcess_);
trigger::size_type l2filterIndex=triggerObj->filterIndex(l2Tag);
if (l2filterIndex<triggerObj->size()) triggerObj->getObjects(l2filterIndex, trigger::TriggerTrack, l2tracks);
std::vector<reco::IsolatedPixelTrackCandidateRef> l3tracks;
edm::InputTag l3Tag = edm::InputTag(l3filterLabel_, "",hltProcess_);
trigger::size_type l3filterIndex=triggerObj->filterIndex(l3Tag);
if (l3filterIndex<triggerObj->size()) triggerObj->getObjects(l3filterIndex, trigger::TriggerTrack, l3tracks);
if (l1CenJets.size()>0||l1ForJets.size()>0||l1CenJets.size()>0) hL2L3acc->Fill(1+0.0001,1);
edm::Handle<reco::IsolatedPixelTrackCandidateCollection> l3col;
edm::InputTag l3colTag=edm::InputTag(l3collectionLabel_,"",hltProcess_);
if (l2tracks.size()>0)
{
hL2L3acc->Fill(2+0.0001,1);
if (useHLTDebug_)
{
iEvent.getByLabel(l3collectionLabel_,l3col);
for (reco::IsolatedPixelTrackCandidateCollection::const_iterator l3it=l3col->begin(); l3it!=l3col->end(); ++l3it)
{
hL3colP->Fill(l3it->track()->pt()*cosh(l3it->track()->eta()),1);
hL3colEta->Fill(l3it->track()->eta(),1);
}
}
}
if (l3tracks.size()>0) hL2L3acc->Fill(3+0.0001,1);
for (unsigned int i=0; i<l1CenJets.size(); i++)
{
hL1pT->Fill(l1CenJets[i]->pt(),1);
}
for (unsigned int i=0; i<l1ForJets.size(); i++)
{
hL1pT->Fill(l1ForJets[i]->pt(),1);
}
for (unsigned int i=0; i<l1TauJets.size(); i++)
{
hL1pT->Fill(l1TauJets[i]->pt(),1);
}
for (unsigned int i=0; i<l2tracks.size(); i++)
{
hL2eta->Fill(l2tracks[i]->track()->eta(),1);
hL2phi->Fill(l2tracks[i]->track()->phi(),1);
hL2pT->Fill(l2tracks[i]->track()->pt(),1);
hL2TowerOccupancy->Fill((l2tracks[i]->towerIndex()).first,(l2tracks[i]->towerIndex()).second,1);
hisopT->Fill(l2tracks[i]->maxPtPxl(),1);
hisopTvsEta->Fill(l2tracks[i]->track()->eta(),l2tracks[i]->maxPtPxl(),1);
if (useHLTDebug_)
{
reco::IsolatedPixelTrackCandidateCollection::const_iterator selTrIt;
double drmin=0.3;
for (reco::IsolatedPixelTrackCandidateCollection::const_iterator l3it=l3col->begin(); l3it!=l3col->end(); ++l3it)
{
double drl2l3=getDist(l3it->track()->eta(),l3it->track()->phi(),l2tracks[i]->track()->eta(),l2tracks[i]->track()->phi());
if (drl2l3<drmin)
{
drmin=drl2l3;
selTrIt=l3it;
}
}
if (drmin!=0.3) hL3candL2rat->Fill(selTrIt->track()->p()/l2tracks[i]->track()->p(),1);
}
}
for (unsigned int i=0; i<l3tracks.size(); i++)
{
hL3Pt->Fill(l3tracks[i]->track()->pt(),1);
hL3eta->Fill(l3tracks[i]->track()->eta(),1);
hL3phi->Fill(l3tracks[i]->track()->phi(),1);
int seltr=-1;
double drmin=100;
for (unsigned int j=0; j<l2tracks.size(); j++)
{
double drl2l3=getDist(l3tracks[i]->track()->eta(),l3tracks[i]->track()->phi(),l2tracks[i]->track()->eta(),l2tracks[i]->track()->phi());
if (drl2l3<drmin)
{
drmin=drl2l3;
seltr=j;
}
}
if (seltr!=-1)
{
hL3L2trackMatch->Fill(drmin,1);
hL3L2rat->Fill(l3tracks[i]->track()->p()/l2tracks[seltr]->track()->p(),1);
}
}
edm::Handle<trigger::TriggerEvent> trevt;
edm::InputTag taodLab=edm::InputTag(hltAODEventTag_,"",hltProcess_);
iEvent.getByLabel(taodLab,trevt);
const trigger::TriggerObjectCollection& TOCol(trevt->getObjects());
trigger::Keys KEYS;
const trigger::size_type nFilt(trevt->sizeFilters());
int nFired=0;
edm::InputTag hltFilterTag_=edm::InputTag(l3filterLabel_,"",hltProcess_);
for (trigger::size_type iFilt=0; iFilt!=nFilt; iFilt++)
{
trigger::Keys KEYS1=trevt->filterKeys(iFilt);
if (KEYS1.size()>0) nFired++;
if (trevt->filterTag(iFilt)==hltFilterTag_) KEYS=trevt->filterKeys(iFilt);
}
trigger::size_type nReg=KEYS.size();
for (trigger::size_type k=0; k<nReg; k++)
{
const trigger::TriggerObject& TObj(TOCol[KEYS[k]]);
hL3etaAOD->Fill(TObj.eta(),1);
hL3pAOD->Fill(TObj.pt()*cosh(TObj.eta()),1);
}
if (nFired==2&&nReg>0)
{
for (trigger::size_type iReg=0; iReg<nReg; iReg++)
{
const trigger::TriggerObject& TObj(TOCol[KEYS[iReg]]);
hL3etaPureAOD->Fill(TObj.eta(),1);
hL3pPureAOD->Fill(TObj.pt()*cosh(TObj.eta()),1);
}
}
}
| void DQMHcalIsoTrackHLT::beginJob | ( | void | ) | [private, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 313 of file DQMHcalIsoTrackHLT.cc.
References DQMStore::book1D(), DQMStore::book2D(), dbe_, folderName_, MonitorElement::getTH2F(), hisopT, hisopTvsEta, hL1pT, hL2eta, hL2L3acc, hL2phi, hL2pT, hL2TowerOccupancy, hL3candL2rat, hL3colEta, hL3colP, hL3eta, hL3etaAOD, hL3etaPureAOD, hL3L2rat, hL3L2trackMatch, hL3pAOD, hL3phi, hL3pPureAOD, hL3Pt, cmsCodeRules::cppFunctionSkipper::operator, MonitorElement::setAxisTitle(), and DQMStore::setCurrentFolder().
{
dbe_ = edm::Service<DQMStore>().operator->();
dbe_->setCurrentFolder(folderName_);
hL2TowerOccupancy=dbe_->book2D("hL2TowerOccupancy","L2 tower occupancy",48,-25,25,73,0,73);
hL2TowerOccupancy->setAxisTitle("ieta",1);
hL2TowerOccupancy->setAxisTitle("iphi",2);
hL2TowerOccupancy->getTH2F()->SetOption("colz");
hL2TowerOccupancy->getTH2F()->SetStats(kFALSE);
hL2L3acc=dbe_->book1D("hL2L3acc","number of L1, L2 and L3 accepts",3,1,4);
hL2L3acc->setAxisTitle("level",1);
hL3L2trackMatch=dbe_->book1D("hL3L2trackMatch","R from L3 object to L2 object ",1000,0,1);
hL3L2trackMatch->setAxisTitle("R(eta,phi)",1);
hL3colP=dbe_->book1D("hL3colP","P of L3 candidates",1000,0,100);
hL3colP->setAxisTitle("P (GeV)",1);
hL3colEta=dbe_->book1D("hL3colEta","eta of L3 candidates",100,-3,3);
hL3colEta->setAxisTitle("eta",1);
hL3candL2rat=dbe_->book1D("hL3candL2rat","ratio of L3 candidate to accepted L2",1000,0,10);
hL3candL2rat->setAxisTitle("P_L3/P_L2",1);
hL3L2rat=dbe_->book1D("hL3L2rat","ratio of L3 to L2 measurement",1000,0,10);
hL3L2rat->setAxisTitle("pT_L3/pT_L2",1);
hL3Pt=dbe_->book1D("hl3Pt","pT of L3 objects",1000,0,100);
hL3Pt->setAxisTitle("pT(GeV)",1);
hL3eta=dbe_->book1D("hl3eta","eta of L3 objects",50,-2.5,2.5);
hL3eta->setAxisTitle("eta",1);
hL3phi=dbe_->book1D("hl3phi","phi of L3 objects",70,-3.5,3.5);
hL3phi->setAxisTitle("phi(rad)",1);
hL3etaAOD=dbe_->book1D("hL3etaAOD","eta of L3 objects (AOD)",50,-2.5,2.5);
hL3etaAOD->setAxisTitle("eta",1);
hL3etaPureAOD=dbe_->book1D("hL3etaPureAOD","eta of L3 objects (AOD, pure)",50,-2.5,2.5);
hL3etaPureAOD->setAxisTitle("eta",1);
hL3pAOD=dbe_->book1D("hl3pAOD","p of L3 objects (AOD)",1000,0,100);
hL3pAOD->setAxisTitle("p(GeV)",1);
hL3pPureAOD=dbe_->book1D("hl3pPureAOD","p of L3 objects (AOD, pure)",1000,0,100);
hL3pPureAOD->setAxisTitle("p(GeV)",1);
hL1pT=dbe_->book1D("hl1pT","pT of L1 objects",1000,0,1000);
hL1pT->setAxisTitle("pT(GeV)",1);
hL2pT=dbe_->book1D("hl2pT","pT of L2 objects",1000,0,1000);
hL2pT->setAxisTitle("pT(GeV)",1);
hL2eta=dbe_->book1D("hl2eta","eta of L2 objects",50,-2.5,2.5);
hL2eta->setAxisTitle("eta",1);
hL2phi=dbe_->book1D("hl2phi","phi of L2 objects",70,-3.5,3.5);
hL2phi->setAxisTitle("phi(rad)",1);
hisopT=dbe_->book1D("hisopT","isolation pT",100,0,5.5);
hisopT->setAxisTitle("iso pT (GeV)",1);
hisopTvsEta=dbe_->book2D("hisopTvsEta","isolation pT vs Eta",8,-2,2,100,0,5.5);
hisopTvsEta->setAxisTitle("eta",1);
hisopTvsEta->setAxisTitle("iso pT (GeV)",2);
}
| void DQMHcalIsoTrackHLT::endJob | ( | void | ) | [private, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 383 of file DQMHcalIsoTrackHLT.cc.
References dbe_, outRootFileName_, DQMStore::save(), and saveToRootFile_.
{
if(dbe_&&saveToRootFile_)
{
dbe_->save(outRootFileName_);
}
}
| double DQMHcalIsoTrackHLT::getDist | ( | double | eta1, |
| double | phi1, | ||
| double | eta2, | ||
| double | phi2 | ||
| ) |
Definition at line 127 of file DQMHcalIsoTrackHLT.cc.
References funct::pow(), and mathSSE::sqrt().
Referenced by analyze().
DQMStore* DQMHcalIsoTrackHLT::dbe_ [private] |
Definition at line 68 of file DQMHcalIsoTrackHLT.cc.
Referenced by beginJob(), and endJob().
int DQMHcalIsoTrackHLT::evtBuf [private] |
Definition at line 66 of file DQMHcalIsoTrackHLT.cc.
std::string DQMHcalIsoTrackHLT::folderName_ [private] |
Definition at line 74 of file DQMHcalIsoTrackHLT.cc.
Referenced by beginJob(), and DQMHcalIsoTrackHLT().
MonitorElement* DQMHcalIsoTrackHLT::hisopT [private] |
Definition at line 121 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and beginJob().
MonitorElement* DQMHcalIsoTrackHLT::hisopTvsEta [private] |
Definition at line 122 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and beginJob().
MonitorElement* DQMHcalIsoTrackHLT::hL1pT [private] |
Definition at line 117 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and beginJob().
MonitorElement* DQMHcalIsoTrackHLT::hL2eta [private] |
Definition at line 118 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and beginJob().
MonitorElement* DQMHcalIsoTrackHLT::hL2L3acc [private] |
Definition at line 95 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and beginJob().
MonitorElement* DQMHcalIsoTrackHLT::hL2phi [private] |
Definition at line 119 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and beginJob().
MonitorElement* DQMHcalIsoTrackHLT::hL2pT [private] |
Definition at line 120 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and beginJob().
Definition at line 94 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and beginJob().
MonitorElement* DQMHcalIsoTrackHLT::hL3candL2rat [private] |
Definition at line 108 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and beginJob().
MonitorElement* DQMHcalIsoTrackHLT::hL3colEta [private] |
Definition at line 104 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and beginJob().
MonitorElement* DQMHcalIsoTrackHLT::hL3colP [private] |
Definition at line 103 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and beginJob().
MonitorElement* DQMHcalIsoTrackHLT::hL3eta [private] |
Definition at line 106 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and beginJob().
MonitorElement* DQMHcalIsoTrackHLT::hL3etaAOD [private] |
Definition at line 111 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and beginJob().
MonitorElement* DQMHcalIsoTrackHLT::hL3etaPureAOD [private] |
Definition at line 113 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and beginJob().
MonitorElement* DQMHcalIsoTrackHLT::hL3L2rat [private] |
Definition at line 96 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and beginJob().
Definition at line 123 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and beginJob().
MonitorElement* DQMHcalIsoTrackHLT::hL3pAOD [private] |
Definition at line 112 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and beginJob().
MonitorElement* DQMHcalIsoTrackHLT::hL3phi [private] |
Definition at line 107 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and beginJob().
MonitorElement* DQMHcalIsoTrackHLT::hL3pPureAOD [private] |
Definition at line 114 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and beginJob().
MonitorElement* DQMHcalIsoTrackHLT::hL3Pt [private] |
Definition at line 100 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and beginJob().
MonitorElement* DQMHcalIsoTrackHLT::hL3pVsEta [private] |
Definition at line 101 of file DQMHcalIsoTrackHLT.cc.
std::string DQMHcalIsoTrackHLT::hltAODEventTag_ [private] |
Definition at line 78 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and DQMHcalIsoTrackHLT().
std::string DQMHcalIsoTrackHLT::hltProcess_ [private] |
Definition at line 86 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and DQMHcalIsoTrackHLT().
std::string DQMHcalIsoTrackHLT::hltRAWEventTag_ [private] |
Definition at line 77 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and DQMHcalIsoTrackHLT().
std::string DQMHcalIsoTrackHLT::l1filterLabel_ [private] |
Definition at line 84 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and DQMHcalIsoTrackHLT().
std::string DQMHcalIsoTrackHLT::l2collectionLabel_ [private] |
Definition at line 80 of file DQMHcalIsoTrackHLT.cc.
Referenced by DQMHcalIsoTrackHLT().
std::string DQMHcalIsoTrackHLT::l2filterLabel_ [private] |
Definition at line 85 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and DQMHcalIsoTrackHLT().
std::string DQMHcalIsoTrackHLT::l3collectionLabel_ [private] |
Definition at line 81 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and DQMHcalIsoTrackHLT().
std::string DQMHcalIsoTrackHLT::l3filterLabel_ [private] |
Definition at line 83 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and DQMHcalIsoTrackHLT().
std::string DQMHcalIsoTrackHLT::outRootFileName_ [private] |
Definition at line 75 of file DQMHcalIsoTrackHLT.cc.
Referenced by DQMHcalIsoTrackHLT(), and endJob().
bool DQMHcalIsoTrackHLT::saveToRootFile_ [private] |
Definition at line 90 of file DQMHcalIsoTrackHLT.cc.
Referenced by DQMHcalIsoTrackHLT(), and endJob().
bool DQMHcalIsoTrackHLT::useHLTDebug_ [private] |
Definition at line 88 of file DQMHcalIsoTrackHLT.cc.
Referenced by analyze(), and DQMHcalIsoTrackHLT().
1.7.3