Inheritance diagram for DTTPAnalyzer:

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

	DTTPAnalyzer (const edm::ParameterSet &)

void	endJob () override

	~DTTPAnalyzer () override

Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () 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

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

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
std::string	getHistoName (const DTLayerId &)

Private Attributes
const edm::EDGetTokenT< DTDigiCollection >	digiToken_

edm::ESHandle< DTGeometry >	dtGeom_

const edm::ESGetToken< DTGeometry, MuonGeometryRecord >	dtGeomToken_

std::map< DTWireId, int >	nDigisPerWire_

TFile *	rootFile_

const bool	subtractT0_

std::map< DTWireId, double >	sumW2PerWire_

std::map< DTWireId, double >	sumWPerWire_

std::unique_ptr< DTTTrigBaseSync >	tTrigSync_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

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

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Author: A. Vilela Pereira

Definition at line 23 of file DTTPAnalyzer.cc.

Constructor & Destructor Documentation

◆ DTTPAnalyzer()

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

Definition at line 65 of file DTTPAnalyzer.cc.

References edm::EDConsumerBase::consumesCollector(), get, muonDTDigis_cfi::pset, rootFile_, CSCSkim_cfi::rootFileName, AlCaHLTBitMon_QueryRunRegistry::string, subtractT0_, and tTrigSync_.

     : subtractT0_(pset.getParameter<bool>("subtractT0")),
       digiToken_(consumes<DTDigiCollection>(pset.getParameter<edm::InputTag>("digiLabel"))),
       dtGeomToken_(esConsumes()) {
   std::string rootFileName = pset.getUntrackedParameter<std::string>("rootFileName");
   rootFile_ = new TFile(rootFileName.c_str(), "RECREATE");
   rootFile_->cd();
 
   if (subtractT0_)
     tTrigSync_ = DTTTrigSyncFactory::get()->create(pset.getParameter<std::string>("tTrigMode"),
                                                    pset.getParameter<edm::ParameterSet>("tTrigModeConfig"),
                                                    consumesCollector());
 }

◆ ~DTTPAnalyzer()

DTTPAnalyzer::~DTTPAnalyzer ( )

override

Definition at line 79 of file DTTPAnalyzer.cc.

References rootFile_.

79 { rootFile_->Close(); }

DTTPAnalyzer::rootFile_

TFile * rootFile_

Definition: DTTPAnalyzer.cc:37

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 81 of file DTTPAnalyzer.cc.

References digiToken_, dtGeom_, dtGeomToken_, nano_mu_digi_cff::layer, nDigisPerWire_, hltrates_dqm_sourceclient-live_cfg::offset, singleTopDQM_cfi::setup, subtractT0_, sumW2PerWire_, sumWPerWire_, FrontierCondition_GT_autoExpress_cfi::t0, and tTrigSync_.

                                                                             {
   // Get the digis from the event
   const edm::Handle<DTDigiCollection>& digis = event.getHandle(digiToken_);
 
   if (subtractT0_) {
     tTrigSync_->setES(setup);
   }
   // Get the DT Geometry
   dtGeom_ = setup.getHandle(dtGeomToken_);
 
   // Iterate through all digi collections ordered by LayerId
   DTDigiCollection::DigiRangeIterator dtLayerIt;
   for (dtLayerIt = digis->begin(); dtLayerIt != digis->end(); ++dtLayerIt) {
     // Get the iterators over the digis associated with this LayerId
     const DTDigiCollection::Range& digiRange = (*dtLayerIt).second;
 
     // Get the layerId
     const DTLayerId layerId = (*dtLayerIt).first;  //FIXME: check to be in the right sector
 
     // Loop over all digis in the given layer
     for (DTDigiCollection::const_iterator digi = digiRange.first; digi != digiRange.second; ++digi) {
       const DTWireId wireId(layerId, (*digi).wire());
 
       double t0 = (*digi).countsTDC();
 
       //FIXME: Reject digis not coming from TP
 
       if (subtractT0_) {
         const DTLayer* layer = nullptr;  //fake
         const GlobalPoint glPt;          //fake
         double offset = tTrigSync_->offset(layer, wireId, glPt);
         t0 -= offset;
       }
 
       if (nDigisPerWire_.find(wireId) == nDigisPerWire_.end()) {
         nDigisPerWire_[wireId] = 0;
         sumWPerWire_[wireId] = 0.;
         sumW2PerWire_[wireId] = 0.;
       }
 
       ++nDigisPerWire_[wireId];
       sumWPerWire_[wireId] += t0;
       sumW2PerWire_[wireId] += t0 * t0;
     }
   }
 }

◆ endJob()

void DTTPAnalyzer::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 128 of file DTTPAnalyzer.cc.

References DTTopology::channels(), dtGeom_, DTTopology::firstChannel(), getHistoName(), HltBtagPostValidation_cff::histoName, crabWrapper::key, DTGeometry::layer(), DTWireId::layerId(), SiStripPI::mean, PFElectronDQMAnalyzer_cfi::nBin, nDigisPerWire_, SiStripPI::rms, rootFile_, DTLayer::specificTopology(), mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, sumW2PerWire_, sumWPerWire_, and DTWireId::wire().

Referenced by o2olib.O2ORunMgr::executeJob().

                           {
   rootFile_->cd();
   std::map<DTLayerId, TH1F*> meanHistoMap;
   std::map<DTLayerId, TH1F*> sigmaHistoMap;
   for (std::map<DTWireId, int>::const_iterator wireIdIt = nDigisPerWire_.begin(); wireIdIt != nDigisPerWire_.end();
        ++wireIdIt) {
     DTWireId wireId((*wireIdIt).first);
 
     int nDigis = nDigisPerWire_[wireId];
     double sumW = sumWPerWire_[wireId];
     double sumW2 = sumW2PerWire_[wireId];
 
     double mean = sumW / nDigis;
     double rms = sumW2 / nDigis - mean * mean;
     rms = sqrt(rms);
 
     DTLayerId layerId = wireId.layerId();
     if (meanHistoMap.find(layerId) == meanHistoMap.end()) {
       std::string histoName = getHistoName(layerId);
       const int firstChannel = dtGeom_->layer(layerId)->specificTopology().firstChannel();
       const int nWires = dtGeom_->layer(layerId)->specificTopology().channels();
       TH1F* meanHistoTP = new TH1F((histoName + "_tpMean").c_str(),
                                    "mean from test pulses by channel",
                                    nWires,
                                    firstChannel,
                                    (firstChannel + nWires));
       TH1F* sigmaHistoTP = new TH1F((histoName + "_tpSigma").c_str(),
                                     "sigma from test pulses by channel",
                                     nWires,
                                     firstChannel,
                                     (firstChannel + nWires));
       meanHistoMap[layerId] = meanHistoTP;
       sigmaHistoMap[layerId] = sigmaHistoTP;
     }
     // Fill the histograms
     int nBin = meanHistoMap[layerId]->GetXaxis()->FindFixBin(wireId.wire());
     meanHistoMap[layerId]->SetBinContent(nBin, mean);
     sigmaHistoMap[layerId]->SetBinContent(nBin, rms);
   }
 
   for (std::map<DTLayerId, TH1F*>::const_iterator key = meanHistoMap.begin(); key != meanHistoMap.end(); ++key) {
     meanHistoMap[(*key).first]->Write();
     sigmaHistoMap[(*key).first]->Write();
   }
 }

◆ getHistoName()

std::string DTTPAnalyzer::getHistoName ( const DTLayerId & lId )

private

Definition at line 174 of file DTTPAnalyzer.cc.

References HltBtagPostValidation_cff::histoName, DTLayerId::layer(), DTChamberId::sector(), DTChamberId::station(), AlCaHLTBitMon_QueryRunRegistry::string, DTSuperLayerId::superlayer(), and DTChamberId::wheel().

Referenced by endJob().

                                                          {
   std::string histoName;
   std::stringstream theStream;
   theStream << "Ch_" << lId.wheel() << "_" << lId.station() << "_" << lId.sector() << "_SL" << lId.superlayer() << "_L"
             << lId.layer();
   theStream >> histoName;
   return histoName;
 }

Member Data Documentation

◆ digiToken_

const edm::EDGetTokenT<DTDigiCollection> DTTPAnalyzer::digiToken_

private

Definition at line 35 of file DTTPAnalyzer.cc.

Referenced by analyze().

◆ dtGeom_

edm::ESHandle<DTGeometry> DTTPAnalyzer::dtGeom_

private

Definition at line 38 of file DTTPAnalyzer.cc.

Referenced by analyze(), and endJob().

◆ dtGeomToken_

const edm::ESGetToken<DTGeometry, MuonGeometryRecord> DTTPAnalyzer::dtGeomToken_

private

Definition at line 39 of file DTTPAnalyzer.cc.

Referenced by analyze().

◆ nDigisPerWire_

std::map<DTWireId, int> DTTPAnalyzer::nDigisPerWire_

private

Definition at line 43 of file DTTPAnalyzer.cc.

Referenced by analyze(), and endJob().

◆ rootFile_

TFile* DTTPAnalyzer::rootFile_

private

Definition at line 37 of file DTTPAnalyzer.cc.

Referenced by DTTPAnalyzer(), endJob(), and ~DTTPAnalyzer().

◆ subtractT0_

const bool DTTPAnalyzer::subtractT0_

private

Definition at line 34 of file DTTPAnalyzer.cc.

Referenced by analyze(), and DTTPAnalyzer().

◆ sumW2PerWire_

std::map<DTWireId, double> DTTPAnalyzer::sumW2PerWire_

private

Definition at line 45 of file DTTPAnalyzer.cc.

Referenced by analyze(), and endJob().

◆ sumWPerWire_

std::map<DTWireId, double> DTTPAnalyzer::sumWPerWire_

private

Definition at line 44 of file DTTPAnalyzer.cc.

Referenced by analyze(), and endJob().

◆ tTrigSync_

std::unique_ptr<DTTTrigBaseSync> DTTPAnalyzer::tTrigSync_

private

Definition at line 40 of file DTTPAnalyzer.cc.

Referenced by analyze(), and DTTPAnalyzer().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ DTTPAnalyzer()

◆ ~DTTPAnalyzer()

Member Function Documentation

◆ analyze()

◆ endJob()

◆ getHistoName()

Member Data Documentation

◆ digiToken_

◆ dtGeom_

◆ dtGeomToken_

◆ nDigisPerWire_

◆ rootFile_

◆ subtractT0_

◆ sumW2PerWire_

◆ sumWPerWire_

◆ tTrigSync_