|
|
|
#include <DTSegmentAnalysisTask.h>
DQM Analysis of 4D DT segments, it produces plots about:
Definition at line 36 of file DTSegmentAnalysisTask.h.
| DTSegmentAnalysisTask::DTSegmentAnalysisTask | ( | const edm::ParameterSet & | pset | ) |
Constructor.
Definition at line 44 of file DTSegmentAnalysisTask.cc.
References checkNoisyChannels, detailedAnalysis, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hltDQMMode, nhitsCut, nLSTimeBin, nTimeBins, cppFunctionSkipper::operator, phiSegmCut, slideTimeBins, theDbe, theRecHits4DLabel, and topHistoFolder.
: nevents(0) , nEventsInLS(0), hNevtPerLS(0) { edm::LogVerbatim ("DTDQM|DTMonitorModule|DTSegmentAnalysisTask") << "[DTSegmentAnalysisTask] Constructor called!"; // switch for detailed analysis detailedAnalysis = pset.getUntrackedParameter<bool>("detailedAnalysis",false); // the name of the 4D rec hits collection theRecHits4DLabel = pset.getParameter<string>("recHits4DLabel"); // Get the map of noisy channels checkNoisyChannels = pset.getUntrackedParameter<bool>("checkNoisyChannels",false); // # of bins in the time histos nTimeBins = pset.getUntrackedParameter<int>("nTimeBins",100); // # of LS per bin in the time histos nLSTimeBin = pset.getUntrackedParameter<int>("nLSTimeBin",2); // switch on/off sliding bins in time histos slideTimeBins = pset.getUntrackedParameter<bool>("slideTimeBins",true); phiSegmCut = pset.getUntrackedParameter<double>("phiSegmCut",30.); nhitsCut = pset.getUntrackedParameter<int>("nhitsCut",12); // Get the DQM needed services theDbe = edm::Service<DQMStore>().operator->(); // top folder for the histograms in DQMStore topHistoFolder = pset.getUntrackedParameter<string>("topHistoFolder","DT/02-Segments"); // hlt DQM mode hltDQMMode = pset.getUntrackedParameter<bool>("hltDQMMode",false); }
| DTSegmentAnalysisTask::~DTSegmentAnalysisTask | ( | ) | [virtual] |
Destructor.
Definition at line 74 of file DTSegmentAnalysisTask.cc.
{
edm::LogVerbatim ("DTDQM|DTMonitorModule|DTSegmentAnalysisTask") << "[DTSegmentAnalysisTask] Destructor called!";
}
| void DTSegmentAnalysisTask::analyze | ( | const edm::Event & | event, |
| const edm::EventSetup & | setup | ||
| ) | [virtual] |
Implements edm::EDAnalyzer.
Definition at line 134 of file DTSegmentAnalysisTask.cc.
References checkNoisyChannels, edm::EventID::event(), MonitorElement::Fill(), fillHistos(), edm::EventSetup::get(), edm::EventBase::id(), edm::HandleBase::isValid(), nEventMonitor, nevents, nEventsInLS, nhitsCut, phiSegmCut, Pi, findQualityFiles::size, DTRecSegment2D::specificRecHits(), crabStatusFromReport::statusMap, theRecHits4DLabel, and xdir.
{
nevents++;
nEventMonitor->Fill(nevents);
nEventsInLS++;
edm::LogVerbatim ("DTDQM|DTMonitorModule|DTSegmentAnalysisTask") << "[DTSegmentAnalysisTask] Analyze #Run: " << event.id().run()
<< " #Event: " << event.id().event();
if(!(event.id().event()%1000))
edm::LogVerbatim ("DTDQM|DTMonitorModule|DTSegmentAnalysisTask") << "[DTSegmentAnalysisTask] Analyze #Run: " << event.id().run()
<< " #Event: " << event.id().event();
ESHandle<DTStatusFlag> statusMap;
if(checkNoisyChannels) {
setup.get<DTStatusFlagRcd>().get(statusMap);
}
// -- 4D segment analysis -----------------------------------------------------
// Get the 4D segment collection from the event
edm::Handle<DTRecSegment4DCollection> all4DSegments;
event.getByLabel(theRecHits4DLabel, all4DSegments);
if(!all4DSegments.isValid()) return;
// Loop over all chambers containing a segment
DTRecSegment4DCollection::id_iterator chamberId;
for (chamberId = all4DSegments->id_begin();
chamberId != all4DSegments->id_end();
++chamberId){
// Get the range for the corresponding ChamerId
DTRecSegment4DCollection::range range = all4DSegments->get(*chamberId);
edm::LogVerbatim ("DTDQM|DTMonitorModule|DTSegmentAnalysisTask") << " Chamber: " << *chamberId << " has " << distance(range.first, range.second)
<< " 4D segments";
// Loop over the rechits of this ChamerId
for (DTRecSegment4DCollection::const_iterator segment4D = range.first;
segment4D!=range.second;
++segment4D){
//FOR NOISY CHANNELS////////////////////////////////
bool segmNoisy = false;
if(checkNoisyChannels) {
if((*segment4D).hasPhi()){
const DTChamberRecSegment2D* phiSeg = (*segment4D).phiSegment();
vector<DTRecHit1D> phiHits = phiSeg->specificRecHits();
map<DTSuperLayerId,vector<DTRecHit1D> > hitsBySLMap;
for(vector<DTRecHit1D>::const_iterator hit = phiHits.begin();
hit != phiHits.end(); ++hit) {
DTWireId wireId = (*hit).wireId();
// Check for noisy channels to skip them
bool isNoisy = false;
bool isFEMasked = false;
bool isTDCMasked = false;
bool isTrigMask = false;
bool isDead = false;
bool isNohv = false;
statusMap->cellStatus(wireId, isNoisy, isFEMasked, isTDCMasked, isTrigMask, isDead, isNohv);
if(isNoisy) {
edm::LogVerbatim ("DTDQM|DTMonitorModule|DTSegmentAnalysisTask") << "Wire: " << wireId << " is noisy, skipping!";
segmNoisy = true;
}
}
}
if((*segment4D).hasZed()) {
const DTSLRecSegment2D* zSeg = (*segment4D).zSegment(); // zSeg lives in the SL RF
// Check for noisy channels to skip them
vector<DTRecHit1D> zHits = zSeg->specificRecHits();
for(vector<DTRecHit1D>::const_iterator hit = zHits.begin();
hit != zHits.end(); ++hit) {
DTWireId wireId = (*hit).wireId();
bool isNoisy = false;
bool isFEMasked = false;
bool isTDCMasked = false;
bool isTrigMask = false;
bool isDead = false;
bool isNohv = false;
statusMap->cellStatus(wireId, isNoisy, isFEMasked, isTDCMasked, isTrigMask, isDead, isNohv);
if(isNoisy) {
edm::LogVerbatim ("DTDQM|DTMonitorModule|DTSegmentAnalysisTask") << "Wire: " << wireId << " is noisy, skipping!";
segmNoisy = true;
}
}
}
} // end of switch on noisy channels
if (segmNoisy) {
edm::LogVerbatim ("DTDQM|DTMonitorModule|DTSegmentAnalysisTask")<<"skipping the segment: it contains noisy cells";
continue;
}
//END FOR NOISY CHANNELS////////////////////////////////
int nHits=0;
if((*segment4D).hasPhi())
nHits = (((*segment4D).phiSegment())->specificRecHits()).size();
if((*segment4D).hasZed())
nHits = nHits + ((((*segment4D).zSegment())->specificRecHits()).size());
double anglePhiSegm(0.);
if( (*segment4D).hasPhi() ) {
double xdir = (*segment4D).phiSegment()->localDirection().x();
double zdir = (*segment4D).phiSegment()->localDirection().z();
anglePhiSegm = atan(xdir/zdir)*180./TMath::Pi();
}
if( fabs(anglePhiSegm) > phiSegmCut ) continue;
// If the segment is in Wh+-2/SecX/MB1, get the DT chambers just above and check if there is a segment
// to validate the segment present in MB1
if( fabs((*chamberId).wheel()) == 2 && (*chamberId).station() == 1 ) {
bool segmOk=false;
int mb(2);
while( mb < 4 ) {
DTChamberId checkMB((*chamberId).wheel(),mb,(*chamberId).sector());
DTRecSegment4DCollection::range ckrange = all4DSegments->get(checkMB);
for (DTRecSegment4DCollection::const_iterator cksegment4D = ckrange.first;
cksegment4D!=ckrange.second;
++cksegment4D){
int nHits=0;
if((*cksegment4D).hasPhi())
nHits = (((*cksegment4D).phiSegment())->specificRecHits()).size();
if((*cksegment4D).hasZed())
nHits = nHits + ((((*cksegment4D).zSegment())->specificRecHits()).size());
if( nHits >= nhitsCut ) segmOk=true;
}
mb++;
}
if( !segmOk ) continue;
}
fillHistos(*chamberId,
nHits,
(*segment4D).chi2()/(*segment4D).degreesOfFreedom());
}
}
// -----------------------------------------------------------------------------
}
| void DTSegmentAnalysisTask::beginLuminosityBlock | ( | edm::LuminosityBlock const & | lumiSeg, |
| edm::EventSetup const & | eSetup | ||
| ) | [virtual] |
Summary.
Reimplemented from edm::EDAnalyzer.
Definition at line 381 of file DTSegmentAnalysisTask.cc.
References nEventsInLS.
{
nEventsInLS = 0;
}
| void DTSegmentAnalysisTask::beginRun | ( | const edm::Run & | run, |
| const edm::EventSetup & | context | ||
| ) | [virtual] |
BeginRun.
Reimplemented from edm::EDAnalyzer.
Definition at line 79 of file DTSegmentAnalysisTask.cc.
References DQMStore::bookFloat(), bookHistos(), chambers, dtGeom, edm::EventSetup::get(), histoTimeEvol, hltDQMMode, hNevtPerLS, nEventMonitor, nLSTimeBin, nTimeBins, DQMStore::setCurrentFolder(), slideTimeBins, theDbe, and topHistoFolder.
{
// Get the DT Geometry
context.get<MuonGeometryRecord>().get(dtGeom);
if (!hltDQMMode) {
theDbe->setCurrentFolder("DT/EventInfo/Counters");
nEventMonitor = theDbe->bookFloat("nProcessedEventsSegment");
}
// loop over all the DT chambers & book the histos
vector<DTChamber*> chambers = dtGeom->chambers();
vector<DTChamber*>::const_iterator ch_it = chambers.begin();
vector<DTChamber*>::const_iterator ch_end = chambers.end();
for (; ch_it != ch_end; ++ch_it) {
bookHistos((*ch_it)->id());
}
// book sector time-evolution histos
int modeTimeHisto = 0;
if(!slideTimeBins) modeTimeHisto = 1;
for(int wheel = -2; wheel != 3; ++wheel) { // loop over wheels
for(int sector = 1; sector <= 12; ++sector) { // loop over sectors
stringstream wheelstr; wheelstr << wheel;
stringstream sectorstr; sectorstr << sector;
string sectorHistoName = "NSegmPerEvent_W" + wheelstr.str()
+ "_Sec" + sectorstr.str();
string sectorHistoTitle = "# segm. W" + wheelstr.str() + " Sect." + sectorstr.str();
theDbe->setCurrentFolder(topHistoFolder + "/Wheel" + wheelstr.str() +
"/Sector" + sectorstr.str());
histoTimeEvol[wheel][sector] = new DTTimeEvolutionHisto(&(*theDbe),sectorHistoName,sectorHistoTitle,
nTimeBins,nLSTimeBin,slideTimeBins,modeTimeHisto);
}
}
if(hltDQMMode) theDbe->setCurrentFolder(topHistoFolder);
else theDbe->setCurrentFolder("DT/EventInfo/");
hNevtPerLS = new DTTimeEvolutionHisto(&(*theDbe),"NevtPerLS","# evt.",nTimeBins,nLSTimeBin,slideTimeBins,2);
}
| void DTSegmentAnalysisTask::bookHistos | ( | DTChamberId | chamberId | ) | [private] |
Definition at line 285 of file DTSegmentAnalysisTask.cc.
References DQMStore::book1D(), DQMStore::book2D(), detailedAnalysis, DQMStore::get(), mergeVDriftHistosByStation::histos, histosPerCh, DTChamberId::sector(), DQMStore::setCurrentFolder(), relativeConstraints::station, DTChamberId::station(), summaryHistos, theDbe, topHistoFolder, and DTChamberId::wheel().
Referenced by beginRun().
{
edm::LogVerbatim ("DTDQM|DTMonitorModule|DTSegmentAnalysisTask") << " Booking histos for chamber: " << chamberId;
// Compose the chamber name
stringstream wheel; wheel << chamberId.wheel();
stringstream station; station << chamberId.station();
stringstream sector; sector << chamberId.sector();
string chamberHistoName =
"_W" + wheel.str() +
"_St" + station.str() +
"_Sec" + sector.str();
for(int wh=-2; wh<=2; wh++){
stringstream wheel; wheel << wh;
theDbe->setCurrentFolder(topHistoFolder + "/Wheel" + wheel.str());
string histoName = "numberOfSegments_W" + wheel.str();
if (theDbe->get(topHistoFolder + "/Wheel" + wheel.str() + "/" + histoName))
continue;
summaryHistos[wh] = theDbe->book2D(histoName.c_str(),histoName.c_str(),12,1,13,4,1,5);
summaryHistos[wh]->setAxisTitle("Sector",1);
summaryHistos[wh]->setBinLabel(1,"1",1);
summaryHistos[wh]->setBinLabel(2,"2",1);
summaryHistos[wh]->setBinLabel(3,"3",1);
summaryHistos[wh]->setBinLabel(4,"4",1);
summaryHistos[wh]->setBinLabel(5,"5",1);
summaryHistos[wh]->setBinLabel(6,"6",1);
summaryHistos[wh]->setBinLabel(7,"7",1);
summaryHistos[wh]->setBinLabel(8,"8",1);
summaryHistos[wh]->setBinLabel(9,"9",1);
summaryHistos[wh]->setBinLabel(10,"10",1);
summaryHistos[wh]->setBinLabel(11,"11",1);
summaryHistos[wh]->setBinLabel(12,"12",1);
summaryHistos[wh]->setBinLabel(1,"MB1",2);
summaryHistos[wh]->setBinLabel(2,"MB2",2);
summaryHistos[wh]->setBinLabel(3,"MB3",2);
summaryHistos[wh]->setBinLabel(4,"MB4",2);
}
theDbe->setCurrentFolder(topHistoFolder + "/Wheel" + wheel.str() +
"/Sector" + sector.str() +
"/Station" + station.str());
// Create the monitor elements
vector<MonitorElement *> histos;
histos.push_back(theDbe->book1D("h4DSegmNHits"+chamberHistoName,
"# of hits per segment",
16, 0.5, 16.5));
if(detailedAnalysis){
histos.push_back(theDbe->book1D("h4DChi2"+chamberHistoName,
"4D Segment reduced Chi2",
20, 0, 20));
}
histosPerCh[chamberId] = histos;
}
| void DTSegmentAnalysisTask::endJob | ( | void | ) | [virtual] |
Endjob.
Reimplemented from edm::EDAnalyzer.
Definition at line 123 of file DTSegmentAnalysisTask.cc.
References hNevtPerLS.
{
edm::LogVerbatim ("DTDQM|DTMonitorModule|DTSegmentAnalysisTask") <<"[DTSegmentAnalysisTask] endjob called!";
delete hNevtPerLS;
//theDbe->save("testMonitoring.root");
}
| void DTSegmentAnalysisTask::endLuminosityBlock | ( | edm::LuminosityBlock const & | lumiSeg, |
| edm::EventSetup const & | eSetup | ||
| ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 369 of file DTSegmentAnalysisTask.cc.
References histoTimeEvol, hNevtPerLS, edm::LuminosityBlockBase::luminosityBlock(), nEventsInLS, and DTTimeEvolutionHisto::updateTimeSlot().
{
hNevtPerLS->updateTimeSlot(lumiSeg.luminosityBlock(), nEventsInLS);
// book sector time-evolution histos
for(int wheel = -2; wheel != 3; ++wheel) {
for(int sector = 1; sector <= 12; ++sector) {
histoTimeEvol[wheel][sector]->updateTimeSlot(lumiSeg.luminosityBlock(), nEventsInLS);
}
}
}
| void DTSegmentAnalysisTask::fillHistos | ( | DTChamberId | chamberId, |
| int | nHits, | ||
| float | chi2 | ||
| ) | [private] |
Definition at line 347 of file DTSegmentAnalysisTask.cc.
References detailedAnalysis, mergeVDriftHistosByStation::histos, histosPerCh, histoTimeEvol, DTChamberId::sector(), DTChamberId::station(), summaryHistos, and DTChamberId::wheel().
Referenced by analyze().
{
int sector = chamberId.sector();
if(chamberId.sector()==13) {
sector = 4;
} else if(chamberId.sector()==14) {
sector = 10;
}
summaryHistos[chamberId.wheel()]->Fill(sector,chamberId.station());
histoTimeEvol[chamberId.wheel()][sector]->accumulateValueTimeSlot(1);
vector<MonitorElement *> histos = histosPerCh[chamberId];
histos[0]->Fill(nHits);
if(detailedAnalysis){
histos[1]->Fill(chi2);
}
}
bool DTSegmentAnalysisTask::checkNoisyChannels [private] |
Definition at line 78 of file DTSegmentAnalysisTask.h.
Referenced by analyze(), and DTSegmentAnalysisTask().
bool DTSegmentAnalysisTask::detailedAnalysis [private] |
Definition at line 69 of file DTSegmentAnalysisTask.h.
Referenced by bookHistos(), DTSegmentAnalysisTask(), and fillHistos().
edm::ESHandle<DTGeometry> DTSegmentAnalysisTask::dtGeom [private] |
Definition at line 72 of file DTSegmentAnalysisTask.h.
Referenced by beginRun().
std::map<DTChamberId, std::vector<MonitorElement*> > DTSegmentAnalysisTask::histosPerCh [private] |
Definition at line 90 of file DTSegmentAnalysisTask.h.
Referenced by bookHistos(), and fillHistos().
std::map<int, std::map<int, DTTimeEvolutionHisto*> > DTSegmentAnalysisTask::histoTimeEvol [private] |
Definition at line 92 of file DTSegmentAnalysisTask.h.
Referenced by beginRun(), endLuminosityBlock(), and fillHistos().
bool DTSegmentAnalysisTask::hltDQMMode [private] |
Definition at line 107 of file DTSegmentAnalysisTask.h.
Referenced by beginRun(), and DTSegmentAnalysisTask().
Definition at line 96 of file DTSegmentAnalysisTask.h.
Referenced by beginRun(), endJob(), and endLuminosityBlock().
Definition at line 113 of file DTSegmentAnalysisTask.h.
Referenced by analyze(), and beginRun().
int DTSegmentAnalysisTask::nevents [private] |
Definition at line 94 of file DTSegmentAnalysisTask.h.
Referenced by analyze().
int DTSegmentAnalysisTask::nEventsInLS [private] |
Definition at line 95 of file DTSegmentAnalysisTask.h.
Referenced by analyze(), beginLuminosityBlock(), and endLuminosityBlock().
int DTSegmentAnalysisTask::nhitsCut [private] |
Definition at line 111 of file DTSegmentAnalysisTask.h.
Referenced by analyze(), and DTSegmentAnalysisTask().
int DTSegmentAnalysisTask::nLSTimeBin [private] |
Definition at line 101 of file DTSegmentAnalysisTask.h.
Referenced by beginRun(), and DTSegmentAnalysisTask().
int DTSegmentAnalysisTask::nTimeBins [private] |
Definition at line 99 of file DTSegmentAnalysisTask.h.
Referenced by beginRun(), and DTSegmentAnalysisTask().
Definition at line 80 of file DTSegmentAnalysisTask.h.
double DTSegmentAnalysisTask::phiSegmCut [private] |
Definition at line 109 of file DTSegmentAnalysisTask.h.
Referenced by analyze(), and DTSegmentAnalysisTask().
bool DTSegmentAnalysisTask::slideTimeBins [private] |
Definition at line 103 of file DTSegmentAnalysisTask.h.
Referenced by beginRun(), and DTSegmentAnalysisTask().
std::map< int, MonitorElement* > DTSegmentAnalysisTask::summaryHistos [private] |
Definition at line 91 of file DTSegmentAnalysisTask.h.
Referenced by bookHistos(), and fillHistos().
DQMStore* DTSegmentAnalysisTask::theDbe [private] |
Definition at line 66 of file DTSegmentAnalysisTask.h.
Referenced by beginRun(), bookHistos(), and DTSegmentAnalysisTask().
std::string DTSegmentAnalysisTask::theRecHits4DLabel [private] |
Definition at line 75 of file DTSegmentAnalysisTask.h.
Referenced by analyze(), and DTSegmentAnalysisTask().
std::string DTSegmentAnalysisTask::topHistoFolder [private] |
Definition at line 105 of file DTSegmentAnalysisTask.h.
Referenced by beginRun(), bookHistos(), and DTSegmentAnalysisTask().
1.7.3