#include <DTt0DBValidation.h>

Inheritance diagram for DTt0DBValidation:

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

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

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

void	endJob ()

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

virtual	~DTt0DBValidation ()
	Destructor. More...

Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

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

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

Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Plot the t0 from the DB

Date:: 2011/08/06 13:37:52

Revision:: 1.5

Author: G. Mila - INFN Torino

Definition at line 33 of file DTt0DBValidation.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 40 of file DTt0DBValidation.cc.

References dbe_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), cppFunctionSkipper::operator, and DQMStore::setCurrentFolder().

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

virtual

Destructor.

Definition at line 64 of file DTt0DBValidation.cc.

64 {}

Member Function Documentation

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

inlinevirtual

Implements edm::EDAnalyzer.

Definition at line 46 of file DTt0DBValidation.h.

46 {}

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

virtual

Operations.

Reimplemented from edm::EDAnalyzer.

Definition at line 67 of file DTt0DBValidation.cc.

References bookHistos(), DTTimeUnits::counts, DTTopology::firstChannel(), edm::EventSetup::get(), DTTopology::lastChannel(), LogTrace, Association::map, and tzero.

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

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

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

References DQMStore::book2D(), dbe_, MonitorElement::setBinLabel(), and DQMStore::setCurrentFolder().

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

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 205 of file DTt0DBValidation.cc.

References dbe_, and DQMStore::save().

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

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

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 170 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 52 of file DTt0DBValidation.h.

edm::ESHandle<DTGeometry> DTt0DBValidation::dtGeom_

private

Definition at line 66 of file DTt0DBValidation.h.

std::string DTt0DBValidation::labelDB_

private

Definition at line 57 of file DTt0DBValidation.h.

std::string DTt0DBValidation::labelDBRef_

private

Definition at line 56 of file DTt0DBValidation.h.

std::string DTt0DBValidation::metname_

private

Definition at line 54 of file DTt0DBValidation.h.

std::string DTt0DBValidation::outputFileName_

private

Definition at line 61 of file DTt0DBValidation.h.

bool DTt0DBValidation::outputMEsInRootFile_

private

Definition at line 60 of file DTt0DBValidation.h.

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

private

Definition at line 77 of file DTt0DBValidation.h.

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

private

Definition at line 74 of file DTt0DBValidation.h.

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

private

Definition at line 73 of file DTt0DBValidation.h.

std::string DTt0DBValidation::t0TestName_

private

Definition at line 63 of file DTt0DBValidation.h.

const DTT0* DTt0DBValidation::tZeroMap_

private

Definition at line 69 of file DTt0DBValidation.h.

const DTT0* DTt0DBValidation::tZeroRefMap_

private

Definition at line 70 of file DTt0DBValidation.h.

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

private

Definition at line 80 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