DTTestPulsesTask Class Reference

#include <DTTestPulsesTask.h>

Inheritance diagram for DTTestPulsesTask:

Public Member Functions
	DTTestPulsesTask (const edm::ParameterSet &ps)
	Constructor. More...
virtual	~DTTestPulsesTask ()
	Destructor. More...
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndex	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) 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	beginRun (const edm::Run &, const edm::EventSetup &)
	BeginRun. More...
void	bookHistos (const DTLayerId &dtLayer, std::string folder, std::string histoTag)
	Book the ME. More...
Protected Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
CurrentProcessingContext const *	currentContext () const
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
edm::ESHandle< DTGeometry >	muonGeom
int	nevents
edm::ParameterSet	parameters
edm::ESHandle< DTRangeT0 >	t0RangeMap
std::pair< int, int >	t0sPeakRange
std::map< int, MonitorElement * >	testPulsesOccupancies
std::map< int, MonitorElement * >	testPulsesProfiles
std::map< int, MonitorElement * >	testPulsesTimeBoxes

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
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 37 of file DTTestPulsesTask.h.

Constructor & Destructor Documentation

DTTestPulsesTask::DTTestPulsesTask

(

const edm::ParameterSet &

ps

)

Constructor.

Definition at line 38 of file DTTestPulsesTask.cc.

References gather_cfg::cout, cppFunctionSkipper::operator, and Parameters::parameters.

                                                           {


  cout<<"[DTTestPulseTask]: Constructor"<<endl;

  parameters = ps;


  t0sPeakRange = make_pair( parameters.getUntrackedParameter<int>("t0sRangeLowerBound", -100), 
                parameters.getUntrackedParameter<int>("t0sRangeUpperBound", 100));

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

}

DTTestPulsesTask::~DTTestPulsesTask ( )

virtual

Destructor.

Definition at line 54 of file DTTestPulsesTask.cc.

References gather_cfg::cout, and nevents.

                                   {

  cout <<"[DTTestPulsesTask]: analyzed " << nevents << " events" << endl;
 
}

Member Function Documentation

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

protectedvirtual

Analyze.

Implements edm::EDAnalyzer.

Definition at line 186 of file DTTestPulsesTask.cc.

References bookHistos(), edm::Event::getByLabel(), nevents, and edm::second().

                                                                       {

  nevents++;
  
  edm::Handle<DTDigiCollection> dtdigis;
  e.getByLabel("dtunpacker", dtdigis);
  
  DTDigiCollection::DigiRangeIterator dtLayerId_It;
  for (dtLayerId_It=dtdigis->begin(); dtLayerId_It!=dtdigis->end(); ++dtLayerId_It){
    
    for (DTDigiCollection::const_iterator digiIt = ((*dtLayerId_It).second).first;
     digiIt!=((*dtLayerId_It).second).second; ++digiIt){
      
      // for clearness..
      int layerIndex = ((*dtLayerId_It).first).rawId();
      int chIndex = ((*dtLayerId_It).first).chamberId().rawId();


      if ((int)(*digiIt).countsTDC() > t0sPeakRange.first &&
      (int)(*digiIt).countsTDC() < t0sPeakRange.second ) {

    // Occupancies
    if (testPulsesOccupancies.find(layerIndex) != testPulsesOccupancies.end())
      testPulsesOccupancies.find(layerIndex)->second->Fill((*digiIt).wire());
    else {
      bookHistos( (*dtLayerId_It).first , string("TPOccupancy"), string("TestPulses") );
      testPulsesOccupancies.find(layerIndex)->second->Fill((*digiIt).wire());
    }

    // Profiles
    if (testPulsesProfiles.find(layerIndex) != testPulsesProfiles.end())
      testPulsesProfiles.find(layerIndex)->second->Fill((*digiIt).wire(),(*digiIt).countsTDC());
    else {
      bookHistos( (*dtLayerId_It).first , string("TPProfile"), string("TestPulses2D") );
      testPulsesProfiles.find(layerIndex)->second->Fill((*digiIt).wire(),(*digiIt).countsTDC());
    }
      }

      // Time Box
      if (testPulsesTimeBoxes.find(chIndex) != testPulsesTimeBoxes.end())
    testPulsesTimeBoxes.find(chIndex)->second->Fill((*digiIt).countsTDC());
      else {
    bookHistos( (*dtLayerId_It).first , string("TPTimeBox"), string("TestPulsesTB") );
    testPulsesTimeBoxes.find(chIndex)->second->Fill((*digiIt).countsTDC());
      }
    }
  }

}

void DTTestPulsesTask::beginJob ( void  )

protectedvirtual

BeginJob.

Reimplemented from edm::EDAnalyzer.

Definition at line 61 of file DTTestPulsesTask.cc.

References gather_cfg::cout, and nevents.

                               {

  cout<<"[DTTestPulsesTask]: BeginJob"<<endl;

  nevents = 0;

}

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

protectedvirtual

BeginRun.

Reimplemented from edm::EDAnalyzer.

Definition at line 69 of file DTTestPulsesTask.cc.

References edm::EventSetup::get().

                                                                               {

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

  // Get the pedestals tTrig (always get it, even if the TPRange is taken from conf)
  //context.get<DTRangeT0Rcd>().get(t0RangeMap);

}

void DTTestPulsesTask::bookHistos ( const DTLayerId &  dtLayer, std::string  folder, std::string  histoTag  )

protected

Book the ME.

Definition at line 79 of file DTTestPulsesTask.cc.

References gather_cfg::cout, DTLayerId, DTLayerId::layer(), Parameters::parameters, DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, DTSuperLayerId::superlayer(), DTLayerId::superlayerId(), and DTChamberId::wheel().

                                                                                          {


  stringstream wheel; wheel << dtLayer.wheel(); 
  stringstream station; station << dtLayer.station();   
  stringstream sector; sector << dtLayer.sector();  
  stringstream superLayer; superLayer << dtLayer.superlayer();  
  stringstream layer; layer << dtLayer.layer(); 

  cout<<"[DTTestPulseTask]: booking"<<endl;

  // TP Profiles  
  if ( folder == "TPProfile" ) {

    const int nWires = muonGeom->layer(DTLayerId(dtLayer.wheel(),
                         dtLayer.station(),
                         dtLayer.sector(),
                         dtLayer.superlayer(),
                         dtLayer.layer()))->specificTopology().channels();

    dbe->setCurrentFolder("DT/DTTestPulsesTask/Wheel" + wheel.str() +
              "/Station" + station.str() +
              "/Sector" + sector.str() + 
              "/SuperLayer" + superLayer.str() + 
              "/" +folder);
    
    string histoName = histoTag 
      + "_W" + wheel.str() 
      + "_St" + station.str() 
      + "_Sec" + sector.str() 
      + "_SL" + superLayer.str() 
      + "_L" + layer.str();

    // Setting the range 
    if ( parameters.getUntrackedParameter<bool>("readDB", false) ) {
      t0RangeMap->slRangeT0( dtLayer.superlayerId() , t0sPeakRange.first, t0sPeakRange.second);
    }
    
    
    cout<<"t0sRangeLowerBound "<<t0sPeakRange.first<<"; "
    <<"t0sRangeUpperBound "<<t0sPeakRange.second<<endl;
    
    
    testPulsesProfiles[int(DTLayerId(dtLayer.wheel(),
                   dtLayer.station(),
                   dtLayer.sector(),
                   dtLayer.superlayer(),
                   dtLayer.layer()).rawId())] =
      dbe->bookProfile(histoName,histoName,
               nWires, 0, nWires, // Xaxis: channels
               t0sPeakRange.first - t0sPeakRange.second, t0sPeakRange.first, t0sPeakRange.second); // Yaxis: times
  }

  // TP Occupancies
  else if ( folder == "TPOccupancy" ) {

    dbe->setCurrentFolder("DT/DTTestPulsesTask/Wheel" + wheel.str() +
              "/Station" + station.str() +
              "/Sector" + sector.str() + 
              "/SuperLayer" + superLayer.str() + 
              "/" +folder);
    
    string histoName = histoTag 
      + "_W" + wheel.str() 
      + "_St" + station.str() 
      + "_Sec" + sector.str() 
      + "_SL" + superLayer.str() 
      + "_L" + layer.str();
    
    const int nWires = muonGeom->layer(DTLayerId(dtLayer.wheel(),
                         dtLayer.station(),
                         dtLayer.sector(),
                         dtLayer.superlayer(),
                         dtLayer.layer()))->specificTopology().channels();
    
    testPulsesOccupancies[int(DTLayerId(dtLayer.wheel(),
                    dtLayer.station(),
                    dtLayer.sector(),
                    dtLayer.superlayer(),
                    dtLayer.layer()).rawId())] =
      dbe->book1D(histoName, histoName, nWires, 0, nWires); 
  }

  // Time Box per Chamber
  else if ( folder == "TPTimeBox" ) {

    dbe->setCurrentFolder("DT/DTTestPulsesTask/Wheel" + wheel.str() +
              "/Station" + station.str() +
              "/Sector" + sector.str() + 
              "/" +folder);
    
    string histoName = histoTag 
      + "_W" + wheel.str() 
      + "_St" + station.str() 
      + "_Sec" + sector.str(); 
    
    testPulsesTimeBoxes[int( DTLayerId(dtLayer.wheel(),
                       dtLayer.station(),
                       dtLayer.sector(),
                       dtLayer.superlayer(),
                       dtLayer.layer()).chamberId().rawId())] = 
      dbe->book1D(histoName, histoName, 10000, 0, 10000); // Overview of the TP (and noise) times
  }

}

Member Data Documentation

DQMStore* DTTestPulsesTask::dbe

private

Definition at line 66 of file DTTestPulsesTask.h.

edm::ESHandle<DTGeometry> DTTestPulsesTask::muonGeom

private

Definition at line 70 of file DTTestPulsesTask.h.

int DTTestPulsesTask::nevents

private

Definition at line 64 of file DTTestPulsesTask.h.

edm::ParameterSet DTTestPulsesTask::parameters

private

Definition at line 68 of file DTTestPulsesTask.h.

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

edm::ESHandle<DTRangeT0> DTTestPulsesTask::t0RangeMap

private

Definition at line 72 of file DTTestPulsesTask.h.

std::pair<int, int> DTTestPulsesTask::t0sPeakRange

private

Definition at line 74 of file DTTestPulsesTask.h.

std::map<int, MonitorElement*> DTTestPulsesTask::testPulsesOccupancies

private

Definition at line 78 of file DTTestPulsesTask.h.

std::map<int, MonitorElement*> DTTestPulsesTask::testPulsesProfiles

private

Definition at line 77 of file DTTestPulsesTask.h.

std::map<int, MonitorElement*> DTTestPulsesTask::testPulsesTimeBoxes

private

Definition at line 79 of file DTTestPulsesTask.h.