DTDigiForNoiseTask Class Reference

#include <DTDigiForNoiseTask.h>

Inheritance diagram for DTDigiForNoiseTask:

Public Member Functions
	DTDigiForNoiseTask (const edm::ParameterSet &ps)
	Constructor. More...
virtual	~DTDigiForNoiseTask ()
	Destructor. More...
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c)
	Analyze. More...
void	beginJob ()
	BeginJob. More...
void	beginLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &context)
	To reset the MEs. More...
void	beginRun (const edm::Run &, const edm::EventSetup &)
	BeginRun. More...
void	bookHistos (const DTLayerId &dtSL)
	Book the ME. More...
void	endJob ()
	Endjob. More...
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 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)

Private Attributes
DQMStore *	dbe
bool	debug
std::map< DTLayerId, MonitorElement * >	digiHistos
edm::EDGetTokenT < DTDigiCollection >	dtDigisToken_
edm::ESHandle< DTGeometry >	muonGeom
int	nevents
edm::ParameterSet	parameters

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)

Detailed Description

Definition at line 40 of file DTDigiForNoiseTask.h.

Constructor & Destructor Documentation

DTDigiForNoiseTask::DTDigiForNoiseTask

(

const edm::ParameterSet &

ps

)

Constructor.

Definition at line 34 of file DTDigiForNoiseTask.cc.

References gather_cfg::cout, debug, edm::ParameterSet::getUntrackedParameter(), cppFunctionSkipper::operator, Parameters::parameters, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                               {

  debug = ps.getUntrackedParameter<bool>("debug", false);
  dtDigisToken_ = consumes<DTDigiCollection>(
      edm::InputTag(ps.getUntrackedParameter<std::string>("diDigisLabel", "dtunpacker")));

  if(debug)
    cout<<"[DTDigiForNoiseTask]: Constructor"<<endl;

  parameters = ps;

  dbe = edm::Service<DQMStore>().operator->();

}

DTDigiForNoiseTask::~DTDigiForNoiseTask ( )

virtual

Destructor.

Definition at line 50 of file DTDigiForNoiseTask.cc.

References gather_cfg::cout, debug, and nevents.

                                       {

  if(debug)
    cout << "DTDigiForNoiseTask: analyzed " << nevents << " events" << endl;

}

Member Function Documentation

void DTDigiForNoiseTask::analyze ( const edm::Event &  e, const edm::EventSetup &  c  )

protectedvirtual

Analyze.

Implements edm::EDAnalyzer.

Definition at line 144 of file DTDigiForNoiseTask.cc.

References bookHistos(), debug, DTTopology::firstChannel(), edm::Event::getByToken(), DTTopology::lastChannel(), and nevents.

                                                                         {

  nevents++;
  //  cout << "events:  " << nevents << endl;
  if (nevents%1000 == 0 && debug) {}

  edm::Handle<DTDigiCollection> dtdigis;
  e.getByToken(dtDigisToken_, dtdigis);

  std::map< int,int > DigiPerWirePerEvent;

  // Loop over all the chambers
  vector<DTChamber*>::const_iterator ch_it = muonGeom->chambers().begin();
  vector<DTChamber*>::const_iterator ch_end = muonGeom->chambers().end();
  // Loop over the SLs
  for (; ch_it != ch_end; ++ch_it) {
    //    DTChamberId ch = (*ch_it)->id();
    vector<const DTSuperLayer*>::const_iterator sl_it = (*ch_it)->superLayers().begin();
    vector<const DTSuperLayer*>::const_iterator sl_end = (*ch_it)->superLayers().end();
    // Loop over the SLs
    for(; sl_it != sl_end; ++sl_it) {
      vector<const DTLayer*>::const_iterator l_it = (*sl_it)->layers().begin();
      vector<const DTLayer*>::const_iterator l_end = (*sl_it)->layers().end();
      // Loop over the Ls
      for(; l_it != l_end; ++l_it) {
    DTLayerId layerId = (*l_it)->id();

    DTDigiCollection::Range layerDigi= dtdigis->get(layerId);
    if(layerDigi.first != layerDigi.second){

      const DTTopology& dtTopo = muonGeom->layer(layerId)->specificTopology();
      const int firstWire = dtTopo.firstChannel();
      const int lastWire = dtTopo.lastChannel();

      if (digiHistos.find(layerId) == digiHistos.end())
        bookHistos(layerId);

      if (digiHistos.find(layerId) != digiHistos.end()){
        for (int wire=firstWire; wire-lastWire <= 0; wire++) {
          DigiPerWirePerEvent[wire]= 0;
        }

        for (DTDigiCollection::const_iterator digi = layerDigi.first;
         digi!=layerDigi.second;
         ++digi){
          DigiPerWirePerEvent[(*digi).wire()]+=1;
        }

        for (int wire=firstWire; wire-lastWire<=0; wire++) {
          digiHistos.find(layerId)->second->Fill(wire,DigiPerWirePerEvent[wire]);
        }
      }
    }

      } //Loop Ls
    } //Loop SLs
  } //Loop over chambers

}

void DTDigiForNoiseTask::beginJob ( void  )

protectedvirtual

BeginJob.

Reimplemented from edm::EDAnalyzer.

Definition at line 68 of file DTDigiForNoiseTask.cc.

References gather_cfg::cout, debug, and nevents.

                                 {

  if(debug)
    cout<<"[DTDigiForNoiseTask]: BeginJob"<<endl;

  nevents = 0;

}

void DTDigiForNoiseTask::beginLuminosityBlock ( edm::LuminosityBlock const &  lumiSeg, edm::EventSetup const &  context  )

protectedvirtual

To reset the MEs.

Reimplemented from edm::EDAnalyzer.

Definition at line 87 of file DTDigiForNoiseTask.cc.

References gather_cfg::cout, debug, timingPdfMaker::histo, edm::LuminosityBlockBase::id(), edm::LuminosityBlockID::luminosityBlock(), and Parameters::parameters.

                                                                                                       {

  if(debug)
    cout<<"[DTDigiForNoiseTask]: Begin of LS transition"<<endl;

  if(lumiSeg.id().luminosityBlock()%parameters.getUntrackedParameter<int>("ResetCycle", 3) == 0) {
    for(map< DTLayerId, MonitorElement* > ::const_iterator histo = digiHistos.begin();
    histo != digiHistos.end();
    histo++) {
      (*histo).second->Reset();
    }
  }

}

void DTDigiForNoiseTask::beginRun ( const edm::Run &  run, const edm::EventSetup &  context  )

protectedvirtual

BeginRun.

Reimplemented from edm::EDAnalyzer.

Definition at line 77 of file DTDigiForNoiseTask.cc.

References gather_cfg::cout, debug, and edm::EventSetup::get().

                                                                                 {

  if(debug)
    cout<<"[DTDigiForNoiseTask]: BeginRun"<<endl;

  // Get the geometry
  context.get<MuonGeometryRecord>().get(muonGeom);

}

void DTDigiForNoiseTask::bookHistos ( const DTLayerId &  dtSL )

protected

Book the ME.

Definition at line 103 of file DTDigiForNoiseTask.cc.

References DTSuperLayerId::chamberId(), gather_cfg::cout, debug, DTTopology::firstChannel(), DTTopology::lastChannel(), DTLayerId::layer(), DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, DTSuperLayerId::superlayer(), DTLayerId::superlayerId(), and DTChamberId::wheel().

                                                        {

  if (debug) cout<<"[DTDigiForNoiseTask]: booking"<<endl;

  const  DTSuperLayerId dtSLId = lId.superlayerId();
  const  DTChamberId dtChId = dtSLId.chamberId();
  stringstream layer; layer << lId.layer();
  stringstream superLayer; superLayer << dtSLId.superlayer();
  stringstream wheel; wheel << dtChId.wheel();
  stringstream station; station << dtChId.station();
  stringstream sector; sector << dtChId.sector();

  dbe->setCurrentFolder("DT/DTDigiForNoiseTask/Wheel" + wheel.str() +
            "/Station" + station.str() +
            "/Sector" + sector.str() + "/DigiPerEvent");

  if (debug){
    cout<<"[DTDigiForNoiseTask]: folder "<< "DT/DTDigiTask/Wheel" + wheel.str() +
      "/Station" + station.str() +
      "/Sector" + sector.str() + "/DigiPerEvent"<<endl;
  }

  string histoName =
    "DigiPerEvent_W" + wheel.str()
    + "_St" + station.str()
    + "_Sec" + sector.str()
    + "_SL" + superLayer.str()
    + "_L" + layer.str();

  if (debug) cout<<"[DTDigiTask]: histoName "<<histoName<<endl;

  const DTTopology& dtTopo = muonGeom->layer(lId)->specificTopology();
  const int firstWire = dtTopo.firstChannel();
  const int lastWire = dtTopo.lastChannel();
  int nWires = lastWire-firstWire+1;

  digiHistos[lId] = dbe->book2D(histoName,histoName,nWires,firstWire,lastWire,10,-0.5,9.5);

}

void DTDigiForNoiseTask::endJob ( void  )

protectedvirtual

Endjob.

Reimplemented from edm::EDAnalyzer.

Definition at line 58 of file DTDigiForNoiseTask.cc.

References gather_cfg::cout, and debug.

                               {

  if(debug)
    cout<<"[DTDigiForNoiseTask] endjob called!"<<endl;

  dbe->rmdir("DT/DTDigiForNoiseTask");

}

Member Data Documentation

DQMStore* DTDigiForNoiseTask::dbe

private

Definition at line 75 of file DTDigiForNoiseTask.h.

bool DTDigiForNoiseTask::debug

private

Definition at line 72 of file DTDigiForNoiseTask.h.

Referenced by pkg.AbstractPkg::generate(), pkg.AbstractPkg::get_kwds(), and pkg.AbstractPkg::write().

std::map< DTLayerId, MonitorElement* > DTDigiForNoiseTask::digiHistos

private

Definition at line 82 of file DTDigiForNoiseTask.h.

edm::EDGetTokenT<DTDigiCollection> DTDigiForNoiseTask::dtDigisToken_

private

Definition at line 80 of file DTDigiForNoiseTask.h.

edm::ESHandle<DTGeometry> DTDigiForNoiseTask::muonGeom

private

Definition at line 79 of file DTDigiForNoiseTask.h.

int DTDigiForNoiseTask::nevents

private

Definition at line 73 of file DTDigiForNoiseTask.h.

edm::ParameterSet DTDigiForNoiseTask::parameters

private

Definition at line 77 of file DTDigiForNoiseTask.h.

Referenced by Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputTags(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::recursePSetProperties().