#include <DTTTrigCalibration.h>

Inheritance diagram for DTTTrigCalibration:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &eventSetup) override
	Fill the time boxes. More...

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

void	endJob () override
	Fit the time box rising edge and write the resulting ttrig to the DB. More...

	~DTTTrigCalibration () override
	Destructor. More...

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

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 const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

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
void	dumpTTrigMap (const DTTtrig *tTrig) const

std::string	getOccupancyName (const DTLayerId &slId) const

std::string	getTBoxName (const DTSuperLayerId &slId) const

void	plotTTrig (const DTTtrig *tTrig) const

Private Attributes
bool	checkNoisyChannels

bool	debug

std::string	digiLabel

bool	doSubtractT0

bool	findTMeanAndSigma

double	kFactor

int	maxDigiPerLayer

int	maxTDCCounts

TFile *	theFile

std::unique_ptr< DTTimeBoxFitter >	theFitter

std::map< DTSuperLayerId, TH1F * >	theHistoMap

std::map< DTLayerId, TH1F * >	theOccupancyMap

edm::ESGetToken< DTStatusFlag, DTStatusFlagRcd >	theStatusMapToken

std::unique_ptr< DTTTrigBaseSync >	theSync

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

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

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

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

Analyzer class which fills time box plots with SL granularity for t_trig computation, fits the rising edge and write results to DB. The time boxes are written to file.

Author: G. Cerminara - INFN Torino

Definition at line 33 of file DTTTrigCalibration.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 40 of file DTTTrigCalibration.cc.

References gather_cfg::cout, debug, dtT0WireCalibration_cfg::digiLabel, DeDxTools::esConsumes(), get, edm::ParameterSet::getUntrackedParameter(), HLT_FULL_cff::InputTag, dtT0Analyzer_cfg::rootFileName, and interactiveExample::theFile.

                                                                   {
   // Get the debug parameter for verbose output
   debug = pset.getUntrackedParameter<bool>("debug");
 
   // Get the label to retrieve digis from the event
   digiLabel = pset.getUntrackedParameter<string>("digiLabel");
   consumes<DTDigiCollection>(edm::InputTag(digiLabel));
 
   // Switch on/off the DB writing
   findTMeanAndSigma = pset.getUntrackedParameter<bool>("fitAndWrite", false);
 
   // The TDC time-window (ns)
   maxTDCCounts = 5000 * pset.getUntrackedParameter<int>("tdcRescale", 1);
   //The maximum number of digis per layer
   maxDigiPerLayer = pset.getUntrackedParameter<int>("maxDigiPerLayer", 10);
 
   // The root file which will contain the histos
   string rootFileName = pset.getUntrackedParameter<string>("rootFileName");
   theFile = new TFile(rootFileName.c_str(), "RECREATE");
   theFile->cd();
   theFitter = std::make_unique<DTTimeBoxFitter>();
   if (debug)
     theFitter->setVerbosity(1);
 
   double sigmaFit = pset.getUntrackedParameter<double>("sigmaTTrigFit", 10.);
   theFitter->setFitSigma(sigmaFit);
 
   doSubtractT0 = pset.getUntrackedParameter<bool>("doSubtractT0", "false");
   // Get the synchronizer
   if (doSubtractT0) {
     theSync = DTTTrigSyncFactory::get()->create(pset.getUntrackedParameter<string>("tTrigMode"),
                                                 pset.getUntrackedParameter<ParameterSet>("tTrigModeConfig"),
                                                 consumesCollector());
   }
 
   checkNoisyChannels = pset.getUntrackedParameter<bool>("checkNoisyChannels", "false");
 
   // the kfactor to be uploaded in the ttrig DB
   kFactor = pset.getUntrackedParameter<double>("kFactor", -0.7);
 
   if (debug)
     cout << "[DTTTrigCalibration]Constructor called!" << endl;
 
   if (checkNoisyChannels) {
     theStatusMapToken = esConsumes();
   }
 }

DTTTrigCalibration::~DTTTrigCalibration ( )

override

Destructor.

Definition at line 89 of file DTTTrigCalibration.cc.

References gather_cfg::cout, debug, and interactiveExample::theFile.

                                         {
   if (debug)
     cout << "[DTTTrigCalibration]Destructor called!" << endl;
 
   //   // Delete all histos
   //   for(map<DTSuperLayerId, TH1F*>::const_iterator slHisto = theHistoMap.begin();
   //       slHisto != theHistoMap.end();
   //       slHisto++) {
   //     delete (*slHisto).second;
   //   }
 
   theFile->Close();
 }

Member Function Documentation

void DTTTrigCalibration::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)

overridevirtual

Fill the time boxes.

Perform the real analysis.

Implements edm::one::EDAnalyzerBase.

Definition at line 104 of file DTTTrigCalibration.cc.

References DTLayer::chamber(), DTSuperLayerId::chamberId(), gather_cfg::cout, debug, dtT0WireCalibration_cfg::digiLabel, edm::EventSetup::getHandle(), phase1PixelTopology::layer, hltrates_dqm_sourceclient-live_cfg::offset, DTLayerId::superlayerId(), and interactiveExample::theFile.

                                                                                        {
   if (debug)
     cout << "[DTTTrigCalibration] #Event: " << event.id().event() << endl;
 
   // Get the digis from the event
   Handle<DTDigiCollection> digis;
   event.getByLabel(digiLabel, digis);
 
   ESHandle<DTStatusFlag> statusMap;
   if (checkNoisyChannels) {
     // Get the map of noisy channels
     statusMap = eventSetup.getHandle(theStatusMapToken);
   }
 
   if (doSubtractT0)
     theSync->setES(eventSetup);
 
   //The chambers too noisy in this event
   vector<DTChamberId> badChambers;
 
   // 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;
 
     const DTLayerId layerId = (*dtLayerIt).first;
     const DTSuperLayerId slId = layerId.superlayerId();
     const DTChamberId chId = slId.chamberId();
     bool badChamber = false;
 
     if (debug)
       cout << "----------- Layer " << layerId << " -------------" << endl;
 
     //Check if the layer is inside a noisy chamber
     for (vector<DTChamberId>::const_iterator chamber = badChambers.begin(); chamber != badChambers.end(); ++chamber) {
       if ((*chamber) == chId) {
         badChamber = true;
         break;
       }
     }
     if (badChamber)
       continue;
 
     //Check if the layer has too many digis
     if ((digiRange.second - digiRange.first) > maxDigiPerLayer) {
       if (debug)
         cout << "Layer " << layerId << "has too many digis (" << (digiRange.second - digiRange.first) << ")" << endl;
       badChambers.push_back(chId);
       continue;
     }
 
     // Get the histo from the map
     TH1F* hTBox = theHistoMap[slId];
     if (hTBox == nullptr) {
       // Book the histogram
       theFile->cd();
       hTBox =
           new TH1F(getTBoxName(slId).c_str(), "Time box (ns)", int(0.25 * 32.0 * maxTDCCounts / 25.0), 0, maxTDCCounts);
       if (debug)
         cout << "  New Time Box: " << hTBox->GetName() << endl;
       theHistoMap[slId] = hTBox;
     }
     TH1F* hO = theOccupancyMap[layerId];
     if (hO == nullptr) {
       // Book the histogram
       theFile->cd();
       hO = new TH1F(getOccupancyName(layerId).c_str(), "Occupancy", 100, 0, 100);
       if (debug)
         cout << "  New Time Box: " << hO->GetName() << endl;
       theOccupancyMap[layerId] = hO;
     }
 
     // Loop over all digis in the given range
     for (DTDigiCollection::const_iterator digi = digiRange.first; digi != digiRange.second; ++digi) {
       const DTWireId wireId(layerId, (*digi).wire());
 
       // Check for noisy channels and skip them
       if (checkNoisyChannels) {
         bool isNoisy = false;
         bool isFEMasked = false;
         bool isTDCMasked = false;
         bool isTrigMask = false;
         bool isDead = false;
         bool isNohv = false;
         statusMap->cellStatus(wireId, isNoisy, isFEMasked, isTDCMasked, isTrigMask, isDead, isNohv);
         if (isNoisy) {
           if (debug)
             cout << "Wire: " << wireId << " is noisy, skipping!" << endl;
           continue;
         }
       }
       theFile->cd();
       double offset = 0;
       if (doSubtractT0) {
         const DTLayer* layer = nullptr;  //fake
         const GlobalPoint glPt;          //fake
         offset = theSync->offset(layer, wireId, glPt);
       }
       hTBox->Fill((*digi).time() - offset);
       if (debug) {
         cout << "   Filling Time Box:   " << hTBox->GetName() << endl;
         cout << "           offset (ns): " << offset << endl;
         cout << "           time(ns):   " << (*digi).time() - offset << endl;
       }
       hO->Fill((*digi).wire());
     }
   }
 }

void DTTTrigCalibration::dumpTTrigMap ( const DTTtrig * tTrig ) const

private

Definition at line 280 of file DTTTrigCalibration.cc.

References DTTtrig::begin(), gather_cfg::cout, and DTTtrig::end().

                                                                 {
   static const double convToNs = 25. / 32.;
   for (DTTtrig::const_iterator ttrig = tTrig->begin(); ttrig != tTrig->end(); ++ttrig) {
     cout << "Wh: " << (*ttrig).first.wheelId << " St: " << (*ttrig).first.stationId
          << " Sc: " << (*ttrig).first.sectorId << " Sl: " << (*ttrig).first.slId
          << " TTrig mean (ns): " << (*ttrig).second.tTrig * convToNs
          << " TTrig sigma (ns): " << (*ttrig).second.tTrms * convToNs << endl;
   }
 }

void DTTTrigCalibration::endJob ( void )

overridevirtual

Fit the time box rising edge and write the resulting ttrig to the DB.

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 214 of file DTTTrigCalibration.cc.

References gather_cfg::cout, debug, DTTimeUnits::ns, DTTtrig::set(), interactiveExample::theFile, and DTCalibDBUtils::writeToDB().

Referenced by o2olib.O2ORunMgr::executeJob().

                                 {
   if (debug)
     cout << "[DTTTrigCalibration]Writing histos to file!" << endl;
 
   // Write all time boxes to file
   theFile->cd();
   for (map<DTSuperLayerId, TH1F*>::const_iterator slHisto = theHistoMap.begin(); slHisto != theHistoMap.end();
        ++slHisto) {
     (*slHisto).second->Write();
   }
   for (map<DTLayerId, TH1F*>::const_iterator slHisto = theOccupancyMap.begin(); slHisto != theOccupancyMap.end();
        ++slHisto) {
     (*slHisto).second->Write();
   }
 
   if (findTMeanAndSigma) {
     // Create the object to be written to DB
     DTTtrig* tTrig = new DTTtrig();
 
     // Loop over the map, fit the histos and write the resulting values to the DB
     for (map<DTSuperLayerId, TH1F*>::const_iterator slHisto = theHistoMap.begin(); slHisto != theHistoMap.end();
          ++slHisto) {
       pair<double, double> meanAndSigma = theFitter->fitTimeBox((*slHisto).second);
       tTrig->set((*slHisto).first, meanAndSigma.first, meanAndSigma.second, kFactor, DTTimeUnits::ns);
 
       if (debug) {
         cout << " SL: " << (*slHisto).first << " mean = " << meanAndSigma.first << " sigma = " << meanAndSigma.second
              << endl;
       }
     }
 
     // Print the ttrig map
     dumpTTrigMap(tTrig);
 
     // Plot the tTrig
     plotTTrig(tTrig);
 
     if (debug)
       cout << "[DTTTrigCalibration]Writing ttrig object to DB!" << endl;
 
     // FIXME: to be read from cfg?
     string tTrigRecord = "DTTtrigRcd";
 
     // Write the object to DB
     DTCalibDBUtils::writeToDB(tTrigRecord, tTrig);
   }
 }

string DTTTrigCalibration::getOccupancyName ( const DTLayerId & slId ) const

private

Definition at line 271 of file DTTTrigCalibration.cc.

References DTLayerId::layer(), DTChamberId::sector(), DTChamberId::station(), DTSuperLayerId::superlayer(), and DTChamberId::wheel().

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

string DTTTrigCalibration::getTBoxName ( const DTSuperLayerId & slId ) const

private

Definition at line 262 of file DTTTrigCalibration.cc.

References DTChamberId::sector(), DTChamberId::station(), DTSuperLayerId::superlayer(), and DTChamberId::wheel().

                                                                        {
   string histoName;
   stringstream theStream;
   theStream << "Ch_" << slId.wheel() << "_" << slId.station() << "_" << slId.sector() << "_SL" << slId.superlayer()
             << "_hTimeBox";
   theStream >> histoName;
   return histoName;
 }

void DTTTrigCalibration::plotTTrig ( const DTTtrig * tTrig ) const

private

Definition at line 290 of file DTTTrigCalibration.cc.

References DTTtrig::begin(), and DTTtrig::end().

                                                              {
   TH1F* tTrig_YB1_Se10 = new TH1F("tTrig_YB1_Se10", "tTrig YB1_Se10", 15, 1, 16);
   TH1F* tTrig_YB2_Se10 = new TH1F("tTrig_YB2_Se10", "tTrig YB2_Se10", 15, 1, 16);
   TH1F* tTrig_YB2_Se11 = new TH1F("tTrig_YB2_Se11", "tTrig YB2_Se11", 12, 1, 13);
 
   static const double convToNs = 25. / 32.;
   for (DTTtrig::const_iterator ttrig = tTrig->begin(); ttrig != tTrig->end(); ++ttrig) {
     // avoid to have wired numbers in the plot
     float tTrigValue = 0;
     float tTrmsValue = 0;
     if ((*ttrig).second.tTrig * convToNs > 0 && (*ttrig).second.tTrig * convToNs < 32000) {
       tTrigValue = (*ttrig).second.tTrig * convToNs;
       tTrmsValue = (*ttrig).second.tTrms * convToNs;
     }
 
     int binx;
     string binLabel;
     stringstream binLabelStream;
     if ((*ttrig).first.sectorId != 14) {
       binx = ((*ttrig).first.stationId - 1) * 3 + (*ttrig).first.slId;
       binLabelStream << "MB" << (*ttrig).first.stationId << "_SL" << (*ttrig).first.slId;
     } else {
       binx = 12 + (*ttrig).first.slId;
       binLabelStream << "MB14_SL" << (*ttrig).first.slId;
     }
     binLabelStream >> binLabel;
 
     if ((*ttrig).first.wheelId == 2) {
       if ((*ttrig).first.sectorId == 10 || (*ttrig).first.sectorId == 14) {
         tTrig_YB2_Se10->Fill(binx, tTrigValue);
         tTrig_YB2_Se10->SetBinError(binx, tTrmsValue);
         tTrig_YB2_Se10->GetXaxis()->SetBinLabel(binx, binLabel.c_str());
         tTrig_YB2_Se10->GetYaxis()->SetTitle("ns");
       } else {
         tTrig_YB2_Se11->Fill(binx, tTrigValue);
         tTrig_YB2_Se11->SetBinError(binx, tTrmsValue);
         tTrig_YB2_Se11->GetXaxis()->SetBinLabel(binx, binLabel.c_str());
         tTrig_YB2_Se11->GetYaxis()->SetTitle("ns");
       }
     } else {
       tTrig_YB1_Se10->Fill(binx, tTrigValue);
       tTrig_YB1_Se10->SetBinError(binx, tTrmsValue);
       tTrig_YB1_Se10->GetXaxis()->SetBinLabel(binx, binLabel.c_str());
       tTrig_YB1_Se10->GetYaxis()->SetTitle("ns");
     }
   }
 
   tTrig_YB1_Se10->Write();
   tTrig_YB2_Se10->Write();
   tTrig_YB2_Se11->Write();
 }

Member Data Documentation

bool DTTTrigCalibration::checkNoisyChannels

private

Definition at line 83 of file DTTTrigCalibration.h.

bool DTTTrigCalibration::debug

private

Definition at line 63 of file DTTTrigCalibration.h.

Referenced by util.rrapi.RRApi::dprint(), rrapi.RRApi::dprint(), pkg.AbstractPkg::generate(), rrapi.RRApi::get(), util.rrapi.RRApi::get(), pkg.AbstractPkg::get_kwds(), runTauIdMVA.TauIDEmbedder::loadMVA_WPs_run2_2017(), runTauIdMVA.TauIDEmbedder::runTauID(), crabFunctions.CrabController::submit(), and pkg.AbstractPkg::write().

std::string DTTTrigCalibration::digiLabel

private

Definition at line 66 of file DTTTrigCalibration.h.

bool DTTTrigCalibration::doSubtractT0

private

Definition at line 81 of file DTTTrigCalibration.h.

bool DTTTrigCalibration::findTMeanAndSigma

private

Definition at line 85 of file DTTTrigCalibration.h.

double DTTTrigCalibration::kFactor

private

Definition at line 87 of file DTTTrigCalibration.h.

int DTTTrigCalibration::maxDigiPerLayer

private

Definition at line 71 of file DTTTrigCalibration.h.

int DTTTrigCalibration::maxTDCCounts

private

Definition at line 69 of file DTTTrigCalibration.h.

TFile* DTTTrigCalibration::theFile

private

Definition at line 74 of file DTTTrigCalibration.h.

std::unique_ptr<DTTimeBoxFitter> DTTTrigCalibration::theFitter

private

Definition at line 90 of file DTTTrigCalibration.h.

std::map<DTSuperLayerId, TH1F*> DTTTrigCalibration::theHistoMap

private

Definition at line 77 of file DTTTrigCalibration.h.

std::map<DTLayerId, TH1F*> DTTTrigCalibration::theOccupancyMap

private

Definition at line 78 of file DTTTrigCalibration.h.

edm::ESGetToken<DTStatusFlag, DTStatusFlagRcd> DTTTrigCalibration::theStatusMapToken

private

Definition at line 93 of file DTTTrigCalibration.h.

std::unique_ptr<DTTTrigBaseSync> DTTTrigCalibration::theSync

private

Definition at line 92 of file DTTTrigCalibration.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation