#include <DTt0DBValidation.h>

Inheritance diagram for DTt0DBValidation:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &setup) override

void	beginRun (const edm::Run &run, const edm::EventSetup &setup) override
	Operations. More...

	DTt0DBValidation (const edm::ParameterSet &pset)
	Constructor. More...

void	endJob () override

void	endRun (edm::Run const &, edm::EventSetup const &) override

	~DTt0DBValidation () override
	Destructor. More...

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Plot the t0 from the DB

Author: G. Mila - INFN Torino

Definition at line 31 of file DTt0DBValidation.h.

Constructor & Destructor Documentation

DTt0DBValidation::DTt0DBValidation ( const edm::ParameterSet & pset )

Constructor.

Definition at line 38 of file DTt0DBValidation.cc.

References dbe_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), Utilities::operator, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                            {
   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 ( )

override

Destructor.

Definition at line 61 of file DTt0DBValidation.cc.

61 {}

Member Function Documentation

void DTt0DBValidation::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

inlineoverride

Definition at line 44 of file DTt0DBValidation.h.

References bookHistos(), and makeMuonMisalignmentScenario::wheel.

44 {}

void DTt0DBValidation::beginRun	(	const edm::Run &	run,
		const edm::EventSetup &	setup
	)

override

Operations.

Definition at line 63 of file DTt0DBValidation.cc.

References bookHistos(), DTTimeUnits::counts, DTTopology::firstChannel(), edm::EventSetup::get(), DTTopology::lastChannel(), LogTrace, genParticles_cff::map, cscNeutronWriter_cfi::t0, tzero, DTT0::version(), and makeMuonMisalignmentScenario::wheel.

                                                                           {
   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 193 of file DTt0DBValidation.cc.

References DTSuperLayerId::chamberId(), dbe_, DTLayerId::layer(), LogTrace, DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, DTSuperLayerId::superlayer(), DTLayerId::superlayerId(), DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

Referenced by analyze().

                                                                             {
   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 225 of file DTt0DBValidation.cc.

References dbe_, and makeMuonMisalignmentScenario::wheel.

                                            {
   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 )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 186 of file DTt0DBValidation.cc.

References dbe_.

                               {
   // Write the histos on a file
   if (outputMEsInRootFile_)
     dbe_->save(outputFileName_);
 }

void DTt0DBValidation::endRun	(	edm::Run const &	run,
		edm::EventSetup const &	setup
	)

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 157 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());
       }
     }
   }
 }

Member Data Documentation

DQMStore* DTt0DBValidation::dbe_

private

Definition at line 50 of file DTt0DBValidation.h.

edm::ESHandle<DTGeometry> DTt0DBValidation::dtGeom_

private

Definition at line 64 of file DTt0DBValidation.h.

std::string DTt0DBValidation::labelDB_

private

Definition at line 55 of file DTt0DBValidation.h.

std::string DTt0DBValidation::labelDBRef_

private

Definition at line 54 of file DTt0DBValidation.h.

std::string DTt0DBValidation::metname_

private

Definition at line 52 of file DTt0DBValidation.h.

std::string DTt0DBValidation::outputFileName_

private

Definition at line 59 of file DTt0DBValidation.h.

bool DTt0DBValidation::outputMEsInRootFile_

private

Definition at line 58 of file DTt0DBValidation.h.

std::map<DTLayerId, MonitorElement *> DTt0DBValidation::t0DiffHistos_

private

Definition at line 75 of file DTt0DBValidation.h.

std::map<DTWireId, std::vector<float> > DTt0DBValidation::t0Map_

private

Definition at line 72 of file DTt0DBValidation.h.

std::map<DTWireId, std::vector<float> > DTt0DBValidation::t0RefMap_

private

Definition at line 71 of file DTt0DBValidation.h.

std::string DTt0DBValidation::t0TestName_

private

Definition at line 61 of file DTt0DBValidation.h.

const DTT0* DTt0DBValidation::tZeroMap_

private

Definition at line 67 of file DTt0DBValidation.h.

const DTT0* DTt0DBValidation::tZeroRefMap_

private

Definition at line 68 of file DTt0DBValidation.h.

std::map<int, MonitorElement *> DTt0DBValidation::wheelSummary_

private

Definition at line 78 of file DTt0DBValidation.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation