|
|
|
#include <DTt0DBValidation.h>
Public Member Functions | |
| void | analyze (const edm::Event &event, const edm::EventSetup &setup) |
| void | beginRun (const edm::Run &run, const edm::EventSetup &setup) |
| Operations. | |
| DTt0DBValidation (const edm::ParameterSet &pset) | |
| Constructor. | |
| void | endJob () |
| void | endRun (edm::Run const &, edm::EventSetup const &) |
| virtual | ~DTt0DBValidation () |
| Destructor. | |
Private Member Functions | |
| void | bookHistos (DTLayerId lId, int firstWire, int lastWire) |
| void | bookHistos (int wheel) |
Private Attributes | |
| DQMStore * | dbe_ |
| edm::ESHandle< DTGeometry > | dtGeom_ |
| std::string | labelDB_ |
| std::string | labelDBRef_ |
| std::string | metname_ |
| std::string | outputFileName_ |
| bool | outputMEsInRootFile_ |
| std::map< DTLayerId, MonitorElement * > | t0DiffHistos_ |
| std::map< DTWireId, std::vector< float > > | t0Map_ |
| std::map< DTWireId, std::vector< float > > | t0RefMap_ |
| std::string | t0TestName_ |
| const DTT0 * | tZeroMap_ |
| const DTT0 * | tZeroRefMap_ |
| std::map< int, MonitorElement * > | wheelSummary_ |
Plot the t0 from the DB
Definition at line 33 of file DTt0DBValidation.h.
| DTt0DBValidation::DTt0DBValidation | ( | const edm::ParameterSet & | pset | ) |
Constructor.
Definition at line 40 of file DTt0DBValidation.cc.
References dbe_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), cmsCodeRules::cppFunctionSkipper::operator, and DQMStore::setCurrentFolder().
{
metname_ = "InterChannelSynchDBValidation";
LogVerbatim(metname_) << "[DTt0DBValidation] Constructor called!";
// Get the DQM needed services
dbe_ = edm::Service<DQMStore>().operator->();
dbe_->setCurrentFolder("DT/DtCalib/InterChannelSynchDBValidation");
// Get dataBase label
labelDBRef_ = pset.getParameter<string>("labelDBRef");
labelDB_ = pset.getParameter<string>("labelDB");
t0TestName_ = "t0DifferenceInRange";
if( pset.exists("t0TestName") ) t0TestName_ = pset.getParameter<string>("t0TestName");
outputMEsInRootFile_ = false;
if( pset.exists("OutputFileName") ){
outputMEsInRootFile_ = true;
outputFileName_ = pset.getParameter<std::string>("OutputFileName");
}
}
| DTt0DBValidation::~DTt0DBValidation | ( | ) | [virtual] |
| void DTt0DBValidation::analyze | ( | const edm::Event & | event, |
| const edm::EventSetup & | setup | ||
| ) | [inline, virtual] |
| void DTt0DBValidation::beginRun | ( | const edm::Run & | run, |
| const edm::EventSetup & | setup | ||
| ) | [virtual] |
Operations.
Reimplemented from edm::EDAnalyzer.
Definition at line 67 of file DTt0DBValidation.cc.
References bookHistos(), DTTimeUnits::counts, DTTopology::firstChannel(), edm::EventSetup::get(), DTTopology::lastChannel(), LogTrace, and tzero.
{
metname_ = "InterChannelSynchDBValidation";
LogVerbatim(metname_) << "[DTt0DBValidation] Parameters initialization";
ESHandle<DTT0> t0_Ref;
setup.get<DTT0Rcd>().get(labelDBRef_, t0_Ref);
tZeroRefMap_ = &*t0_Ref;
LogVerbatim(metname_) << "[DTt0DBValidation] reference T0 version: " << t0_Ref->version();
ESHandle<DTT0> t0;
setup.get<DTT0Rcd>().get(labelDB_, t0);
tZeroMap_ = &*t0;
LogVerbatim(metname_) << "[DTt0DBValidation] T0 to validate version: " << t0->version();
//book&reset the summary histos
for(int wheel=-2; wheel<=2; wheel++){
bookHistos(wheel);
wheelSummary_[wheel]->Reset();
}
// Get the geometry
setup.get<MuonGeometryRecord>().get(dtGeom_);
// Loop over Ref DB entries
for(DTT0::const_iterator tzero = tZeroRefMap_->begin();
tzero != tZeroRefMap_->end(); tzero++) {
// t0s and rms are TDC counts
// @@@ NEW DTT0 FORMAT
// DTWireId wireId((*tzero).first.wheelId,
// (*tzero).first.stationId,
// (*tzero).first.sectorId,
// (*tzero).first.slId,
// (*tzero).first.layerId,
// (*tzero).first.cellId);
int channelId = tzero->channelId;
if ( channelId == 0 ) continue;
DTWireId wireId(channelId);
// @@@ NEW DTT0 END
float t0mean;
float t0rms;
tZeroRefMap_->get( wireId, t0mean, t0rms, DTTimeUnits::counts );
LogTrace(metname_)<< "Ref Wire: " << wireId <<endl
<< " T0 mean (TDC counts): " << t0mean
<< " T0_rms (TDC counts): " << t0rms;
t0RefMap_[wireId].push_back(t0mean);
t0RefMap_[wireId].push_back(t0rms);
}
// Loop over Ref DB entries
for(DTT0::const_iterator tzero = tZeroMap_->begin();
tzero != tZeroMap_->end(); tzero++) {
// t0s and rms are TDC counts
// @@@ NEW DTT0 FORMAT
// DTWireId wireId((*tzero).first.wheelId,
// (*tzero).first.stationId,
// (*tzero).first.sectorId,
// (*tzero).first.slId,
// (*tzero).first.layerId,
// (*tzero).first.cellId);
int channelId = tzero->channelId;
if ( channelId == 0 ) continue;
DTWireId wireId(channelId);
// @@@ NEW DTT0 END
float t0mean;
float t0rms;
tZeroMap_->get( wireId, t0mean, t0rms, DTTimeUnits::counts );
LogTrace(metname_)<< "Wire: " << wireId <<endl
<< " T0 mean (TDC counts): " << t0mean
<< " T0_rms (TDC counts): " << t0rms;
t0Map_[wireId].push_back(t0mean);
t0Map_[wireId].push_back(t0rms);
}
double difference = 0;
for(map<DTWireId, vector<float> >::const_iterator theMap = t0RefMap_.begin();
theMap != t0RefMap_.end();
theMap++) {
if(t0Map_.find((*theMap).first) != t0Map_.end()) {
// Compute the difference
difference = t0Map_[(*theMap).first][0]-(*theMap).second[0];
//book histo
DTLayerId layerId = (*theMap).first.layerId();
if(t0DiffHistos_.find(layerId) == t0DiffHistos_.end()) {
const DTTopology& dtTopo = dtGeom_->layer(layerId)->specificTopology();
const int firstWire = dtTopo.firstChannel();
const int lastWire = dtTopo.lastChannel();
bookHistos(layerId, firstWire, lastWire);
}
LogTrace(metname_)<< "Filling the histo for wire: "<<(*theMap).first
<<" difference: "<<difference;
t0DiffHistos_[layerId]->Fill((*theMap).first.wire(),difference);
}
} // Loop over the t0 map reference
}
| void DTt0DBValidation::bookHistos | ( | DTLayerId | lId, |
| int | firstWire, | ||
| int | lastWire | ||
| ) | [private] |
Definition at line 211 of file DTt0DBValidation.cc.
References DQMStore::book1D(), DTSuperLayerId::chamberId(), dbe_, DTLayerId::layer(), LogTrace, DTChamberId::sector(), DQMStore::setCurrentFolder(), relativeConstraints::station, DTChamberId::station(), DTSuperLayerId::superlayer(), DTLayerId::superlayerId(), and DTChamberId::wheel().
{
LogTrace(metname_)<< " Booking histos for L: " << lId;
// Compose the chamber name
stringstream wheel; wheel << lId.superlayerId().chamberId().wheel();
stringstream station; station << lId.superlayerId().chamberId().station();
stringstream sector; sector << lId.superlayerId().chamberId().sector();
stringstream superLayer; superLayer << lId.superlayerId().superlayer();
stringstream layer; layer << lId.layer();
string lHistoName =
"_W" + wheel.str() +
"_St" + station.str() +
"_Sec" + sector.str() +
"_SL" + superLayer.str()+
"_L" + layer.str();
dbe_->setCurrentFolder("DT/DtCalib/InterChannelSynchDBValidation/Wheel" + wheel.str() +
"/Station" + station.str() +
"/Sector" + sector.str() +
"/SuperLayer" +superLayer.str());
// Create the monitor elements
MonitorElement * hDifference;
hDifference = dbe_->book1D("T0Difference"+lHistoName, "difference between the two t0 values",lastWire-firstWire+1, firstWire-0.5, lastWire+0.5);
t0DiffHistos_[lId] = hDifference;
}
| void DTt0DBValidation::bookHistos | ( | int | wheel | ) | [private] |
Definition at line 241 of file DTt0DBValidation.cc.
References DQMStore::book2D(), dbe_, MonitorElement::setBinLabel(), and DQMStore::setCurrentFolder().
{
dbe_->setCurrentFolder("DT/DtCalib/InterChannelSynchDBValidation");
stringstream wh; wh << wheel;
wheelSummary_[wheel]= dbe_->book2D("SummaryWrongT0_W"+wh.str(), "W"+wh.str()+": summary of wrong t0 differences",44,1,45,14,1,15);
wheelSummary_[wheel]->setBinLabel(1,"M1L1",1);
wheelSummary_[wheel]->setBinLabel(2,"M1L2",1);
wheelSummary_[wheel]->setBinLabel(3,"M1L3",1);
wheelSummary_[wheel]->setBinLabel(4,"M1L4",1);
wheelSummary_[wheel]->setBinLabel(5,"M1L5",1);
wheelSummary_[wheel]->setBinLabel(6,"M1L6",1);
wheelSummary_[wheel]->setBinLabel(7,"M1L7",1);
wheelSummary_[wheel]->setBinLabel(8,"M1L8",1);
wheelSummary_[wheel]->setBinLabel(9,"M1L9",1);
wheelSummary_[wheel]->setBinLabel(10,"M1L10",1);
wheelSummary_[wheel]->setBinLabel(11,"M1L11",1);
wheelSummary_[wheel]->setBinLabel(12,"M1L12",1);
wheelSummary_[wheel]->setBinLabel(13,"M2L1",1);
wheelSummary_[wheel]->setBinLabel(14,"M2L2",1);
wheelSummary_[wheel]->setBinLabel(15,"M2L3",1);
wheelSummary_[wheel]->setBinLabel(16,"M2L4",1);
wheelSummary_[wheel]->setBinLabel(17,"M2L5",1);
wheelSummary_[wheel]->setBinLabel(18,"M2L6",1);
wheelSummary_[wheel]->setBinLabel(19,"M2L7",1);
wheelSummary_[wheel]->setBinLabel(20,"M2L8",1);
wheelSummary_[wheel]->setBinLabel(21,"M2L9",1);
wheelSummary_[wheel]->setBinLabel(22,"M2L10",1);
wheelSummary_[wheel]->setBinLabel(23,"M2L11",1);
wheelSummary_[wheel]->setBinLabel(24,"M2L12",1);
wheelSummary_[wheel]->setBinLabel(25,"M3L1",1);
wheelSummary_[wheel]->setBinLabel(26,"M3L2",1);
wheelSummary_[wheel]->setBinLabel(27,"M3L3",1);
wheelSummary_[wheel]->setBinLabel(28,"M3L4",1);
wheelSummary_[wheel]->setBinLabel(29,"M3L5",1);
wheelSummary_[wheel]->setBinLabel(30,"M3L6",1);
wheelSummary_[wheel]->setBinLabel(31,"M3L7",1);
wheelSummary_[wheel]->setBinLabel(32,"M3L8",1);
wheelSummary_[wheel]->setBinLabel(33,"M3L9",1);
wheelSummary_[wheel]->setBinLabel(34,"M3L10",1);
wheelSummary_[wheel]->setBinLabel(35,"M3L11",1);
wheelSummary_[wheel]->setBinLabel(36,"M3L12",1);
wheelSummary_[wheel]->setBinLabel(37,"M4L1",1);
wheelSummary_[wheel]->setBinLabel(38,"M4L2",1);
wheelSummary_[wheel]->setBinLabel(39,"M4L3",1);
wheelSummary_[wheel]->setBinLabel(40,"M4L4",1);
wheelSummary_[wheel]->setBinLabel(41,"M4L5",1);
wheelSummary_[wheel]->setBinLabel(42,"M4L6",1);
wheelSummary_[wheel]->setBinLabel(43,"M4L7",1);
wheelSummary_[wheel]->setBinLabel(44,"M4L8",1);
}
| void DTt0DBValidation::endJob | ( | void | ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 205 of file DTt0DBValidation.cc.
References dbe_, and DQMStore::save().
{
// Write the histos on a file
if(outputMEsInRootFile_) dbe_->save(outputFileName_);
}
| void DTt0DBValidation::endRun | ( | edm::Run const & | run, |
| edm::EventSetup const & | setup | ||
| ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 170 of file DTt0DBValidation.cc.
References QReport::getBadChannels(), QReport::getMessage(), QReport::getQTresult(), and QReport::getStatus().
{
// Check the histos
string testCriterionName = t0TestName_;
for(map<DTLayerId, MonitorElement*>::const_iterator hDiff = t0DiffHistos_.begin();
hDiff != t0DiffHistos_.end();
hDiff++) {
const QReport * theDiffQReport = (*hDiff).second->getQReport(testCriterionName);
if(theDiffQReport) {
int xBin = ((*hDiff).first.station()-1)*12+(*hDiff).first.layer()+4*((*hDiff).first.superlayer()-1);
if( (*hDiff).first.station()==4 && (*hDiff).first.superlayer()==3 )
xBin = ((*hDiff).first.station()-1)*12+(*hDiff).first.layer()+4*((*hDiff).first.superlayer()-2);
int qReportStatus = theDiffQReport->getStatus()/100;
wheelSummary_[(*hDiff).first.wheel()]->setBinContent(xBin,(*hDiff).first.sector(),qReportStatus);
LogVerbatim(metname_) << "-------- layer: " << (*hDiff).first << " " << theDiffQReport->getMessage()
<< " ------- " << theDiffQReport->getStatus()
<< " ------- " << setprecision(3) << theDiffQReport->getQTresult();
vector<dqm::me_util::Channel> badChannels = theDiffQReport->getBadChannels();
for (vector<dqm::me_util::Channel>::iterator channel = badChannels.begin();
channel != badChannels.end(); channel++) {
LogVerbatim(metname_) << "layer: " << (*hDiff).first << " Bad channel: "
<< (*channel).getBin() << " Contents : "
<< (*channel).getContents();
//wheelSummary_[(*hDiff).first.wheel()]->Fill(xBin,(*hDiff).first.sector());
}
}
}
}
DQMStore* DTt0DBValidation::dbe_ [private] |
Definition at line 52 of file DTt0DBValidation.h.
edm::ESHandle<DTGeometry> DTt0DBValidation::dtGeom_ [private] |
Definition at line 66 of file DTt0DBValidation.h.
std::string DTt0DBValidation::labelDB_ [private] |
Definition at line 57 of file DTt0DBValidation.h.
std::string DTt0DBValidation::labelDBRef_ [private] |
Definition at line 56 of file DTt0DBValidation.h.
std::string DTt0DBValidation::metname_ [private] |
Definition at line 54 of file DTt0DBValidation.h.
std::string DTt0DBValidation::outputFileName_ [private] |
Definition at line 61 of file DTt0DBValidation.h.
bool DTt0DBValidation::outputMEsInRootFile_ [private] |
Definition at line 60 of file DTt0DBValidation.h.
std::map<DTLayerId, MonitorElement* > DTt0DBValidation::t0DiffHistos_ [private] |
Definition at line 77 of file DTt0DBValidation.h.
std::map<DTWireId, std::vector<float> > DTt0DBValidation::t0Map_ [private] |
Definition at line 74 of file DTt0DBValidation.h.
std::map<DTWireId, std::vector<float> > DTt0DBValidation::t0RefMap_ [private] |
Definition at line 73 of file DTt0DBValidation.h.
std::string DTt0DBValidation::t0TestName_ [private] |
Definition at line 63 of file DTt0DBValidation.h.
const DTT0* DTt0DBValidation::tZeroMap_ [private] |
Definition at line 69 of file DTt0DBValidation.h.
const DTT0* DTt0DBValidation::tZeroRefMap_ [private] |
Definition at line 70 of file DTt0DBValidation.h.
std::map<int, MonitorElement* > DTt0DBValidation::wheelSummary_ [private] |
Definition at line 80 of file DTt0DBValidation.h.
1.7.3