DTTPAnalyzer Class Reference

Inheritance diagram for DTTPAnalyzer:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	beginRun (const edm::Run &, const edm::EventSetup &) override
	DTTPAnalyzer (const edm::ParameterSet &)
void	endJob () override
	~DTTPAnalyzer () override
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 &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > 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 &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

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

Private Attributes
edm::InputTag	digiLabel_
edm::ESHandle< DTGeometry >	dtGeom_
std::map< DTWireId, int >	nDigisPerWire_
TFile *	rootFile_
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::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	wantsInputProcessBlocks ()
static bool	wantsProcessBlocks ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Author

A. Vilela Pereira

Definition at line 20 of file DTTPAnalyzer.cc.

Constructor & Destructor Documentation

DTTPAnalyzer()

DTTPAnalyzer::DTTPAnalyzer

(

const edm::ParameterSet &

pset

)

Definition at line 68 of file DTTPAnalyzer.cc.

    : subtractT0_(pset.getParameter<bool>("subtractT0")), digiLabel_(pset.getParameter<edm::InputTag>("digiLabel")) {
  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"));
}

References get, muonDTDigis_cfi::pset, rootFile_, CSCSkim_cfi::rootFileName, AlCaHLTBitMon_QueryRunRegistry::string, subtractT0_, and tTrigSync_.

~DTTPAnalyzer()

DTTPAnalyzer::~DTTPAnalyzer ( )

override

Definition at line 79 of file DTTPAnalyzer.cc.

{ rootFile_->Close(); }

References rootFile_.

Member Function Documentation

analyze()

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

overridevirtual

Implements edm::EDAnalyzer.

Definition at line 93 of file DTTPAnalyzer.cc.

                                                                            {
  // Get the digis from the event
  edm::Handle<DTDigiCollection> digis;
  event.getByLabel(digiLabel_, digis);
 
  // 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;
    }
  }
}

References digiLabel_, phase1PixelTopology::layer, nDigisPerWire_, hltrates_dqm_sourceclient-live_cfg::offset, subtractT0_, sumW2PerWire_, sumWPerWire_, FrontierCondition_GT_autoExpress_cfi::t0, and tTrigSync_.

beginRun()

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

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 81 of file DTTPAnalyzer.cc.

                                                                         {
  // Get the t0 map from the DB
  if (subtractT0_) {
    /*ESHandle<DTT0> t0;
     setup.get<DTT0Rcd>().get(t0);
     tZeroMap_ = &*t0;*/
    tTrigSync_->setES(setup);
  }
  // Get the DT Geometry
  setup.get<MuonGeometryRecord>().get(dtGeom_);
}

References dtGeom_, get, singleTopDQM_cfi::setup, subtractT0_, and tTrigSync_.

endJob()

void DTTPAnalyzer::endJob ( void  )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 135 of file DTTPAnalyzer.cc.

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

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

getHistoName()

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

private

Definition at line 181 of file DTTPAnalyzer.cc.

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

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

Referenced by endJob().

Member Data Documentation

digiLabel_

edm::InputTag DTTPAnalyzer::digiLabel_

private

Definition at line 34 of file DTTPAnalyzer.cc.

Referenced by analyze().

dtGeom_

edm::ESHandle<DTGeometry> DTTPAnalyzer::dtGeom_

private

Definition at line 38 of file DTTPAnalyzer.cc.

Referenced by beginRun(), and endJob().

nDigisPerWire_

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

private

Definition at line 42 of file DTTPAnalyzer.cc.

Referenced by analyze(), and endJob().

rootFile_

TFile* DTTPAnalyzer::rootFile_

private

Definition at line 36 of file DTTPAnalyzer.cc.

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

subtractT0_

bool DTTPAnalyzer::subtractT0_

private

Definition at line 33 of file DTTPAnalyzer.cc.

Referenced by analyze(), beginRun(), and DTTPAnalyzer().

sumW2PerWire_

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

private

Definition at line 44 of file DTTPAnalyzer.cc.

Referenced by analyze(), and endJob().

sumWPerWire_

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

private

Definition at line 43 of file DTTPAnalyzer.cc.

Referenced by analyze(), and endJob().

tTrigSync_

std::unique_ptr<DTTTrigBaseSync> DTTPAnalyzer::tTrigSync_

private

Definition at line 39 of file DTTPAnalyzer.cc.

Referenced by analyze(), beginRun(), and DTTPAnalyzer().