#include <DTTTrigOffsetCalibration.h>

Inheritance diagram for DTTTrigOffsetCalibration:

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

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

	DTTTrigOffsetCalibration (const edm::ParameterSet &pset)

void	endJob () override

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

	~DTTTrigOffsetCalibration () override

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns >
	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

ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESResolverIndex > 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Types
typedef std::map< DTChamberId, std::vector< TH1F * > >	ChamberHistosMap

Private Member Functions
void	bookHistos (DTChamberId)

Private Attributes
const bool	doTTrigCorrection_

const edm::ESGetToken< DTGeometry, MuonGeometryRecord >	dtGeomToken_

TFile *	rootFile_

DTSegmentSelector *	select_

const std::string	theCalibChamber_

const edm::EDGetTokenT< DTRecSegment4DCollection >	theRecHits4DToken_

ChamberHistosMap	theT0SegHistoMap_

const DTTtrig *	tTrigMap_

const edm::ESGetToken< DTTtrig, DTTtrigRcd >	ttrigToken_

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

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

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

No description available.

Author: A. Vilela Pereira

Definition at line 31 of file DTTTrigOffsetCalibration.h.

Member Typedef Documentation

◆ ChamberHistosMap

typedef std::map<DTChamberId, std::vector<TH1F*> > DTTTrigOffsetCalibration::ChamberHistosMap

private

Definition at line 44 of file DTTTrigOffsetCalibration.h.

Constructor & Destructor Documentation

◆ DTTTrigOffsetCalibration()

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

Definition at line 36 of file DTTTrigOffsetCalibration.cc.

References edm::BeginRun, edm::EDConsumerBase::consumesCollector(), muonDTDigis_cfi::pset, rootFile_, CSCSkim_cfi::rootFileName, and select_.

     : theRecHits4DToken_(consumes<DTRecSegment4DCollection>(pset.getParameter<InputTag>("recHits4DLabel"))),
       doTTrigCorrection_(pset.getUntrackedParameter<bool>("doT0SegCorrection", false)),
       theCalibChamber_(pset.getUntrackedParameter<string>("calibChamber", "All")),
       ttrigToken_(
           esConsumes<edm::Transition::BeginRun>(edm::ESInputTag("", pset.getUntrackedParameter<string>("dbLabel")))),
       dtGeomToken_(esConsumes()) {
   LogVerbatim("Calibration") << "[DTTTrigOffsetCalibration] Constructor called!";
 
   edm::ConsumesCollector collector(consumesCollector());
   select_ = new DTSegmentSelector(pset, collector);
 
   // the root file which will contain the histos
   string rootFileName = pset.getUntrackedParameter<string>("rootFileName", "DTT0SegHistos.root");
   rootFile_ = new TFile(rootFileName.c_str(), "RECREATE");
   rootFile_->cd();
 }

◆ ~DTTTrigOffsetCalibration()

DTTTrigOffsetCalibration::~DTTTrigOffsetCalibration ( )

override

Definition at line 63 of file DTTTrigOffsetCalibration.cc.

References rootFile_.

                                                     {
   rootFile_->Close();
   LogVerbatim("Calibration") << "[DTTTrigOffsetCalibration] Destructor called!";
 }

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 68 of file DTTTrigOffsetCalibration.cc.

References bookHistos(), relativeConstraints::chamber, DTGeometry::chamber(), dtGeomToken_, edmPickEvents::event, options_cfi::eventSetup, HcalObjRepresent::Fill(), LogTrace, FastTimerService_cff::range, rootFile_, select_, theCalibChamber_, theRecHits4DToken_, and theT0SegHistoMap_.

                                                                                        {
   rootFile_->cd();
   DTChamberId chosenChamberId;
 
   if (theCalibChamber_ != "All") {
     stringstream linestr;
     int selWheel, selStation, selSector;
     linestr << theCalibChamber_;
     linestr >> selWheel >> selStation >> selSector;
     chosenChamberId = DTChamberId(selWheel, selStation, selSector);
     LogVerbatim("Calibration") << " chosen chamber " << chosenChamberId << endl;
   }
 
   // Get the DT Geometry
   ESHandle<DTGeometry> dtGeom;
   dtGeom = eventSetup.getHandle(dtGeomToken_);
 
   // Get the rechit collection from the event
   const Handle<DTRecSegment4DCollection>& all4DSegments = event.getHandle(theRecHits4DToken_);
 
   // Loop over segments by chamber
   DTRecSegment4DCollection::id_iterator chamberIdIt;
   for (chamberIdIt = all4DSegments->id_begin(); chamberIdIt != all4DSegments->id_end(); ++chamberIdIt) {
     // Get the chamber from the setup
     const DTChamber* chamber = dtGeom->chamber(*chamberIdIt);
     LogTrace("Calibration") << "Chamber Id: " << *chamberIdIt;
 
     // Book histos
     if (theT0SegHistoMap_.find(*chamberIdIt) == theT0SegHistoMap_.end()) {
       bookHistos(*chamberIdIt);
     }
 
     // Calibrate just the chosen chamber/s
     if ((theCalibChamber_ != "All") && ((*chamberIdIt) != chosenChamberId))
       continue;
 
     // Get the range for the corresponding ChamberId
     DTRecSegment4DCollection::range range = all4DSegments->get((*chamberIdIt));
 
     // Loop over the rechits of this DetUnit
     for (DTRecSegment4DCollection::const_iterator segment = range.first; segment != range.second; ++segment) {
       LogTrace("Calibration") << "Segment local pos (in chamber RF): " << (*segment).localPosition()
                               << "\nSegment global pos: " << chamber->toGlobal((*segment).localPosition());
 
       if (!(*select_)(*segment, event, eventSetup))
         continue;
 
       // Fill t0-seg values
       if ((*segment).hasPhi()) {
         //if( segment->phiSegment()->ist0Valid() ){
         if ((segment->phiSegment()->t0()) != 0.00) {
           (theT0SegHistoMap_[*chamberIdIt])[0]->Fill(segment->phiSegment()->t0());
         }
       }
       if ((*segment).hasZed()) {
         //if( segment->zSegment()->ist0Valid() ){
         if ((segment->zSegment()->t0()) != 0.00) {
           (theT0SegHistoMap_[*chamberIdIt])[1]->Fill(segment->zSegment()->t0());
         }
       }
     }  // DTRecSegment4DCollection::const_iterator segment
   }    // DTRecSegment4DCollection::id_iterator chamberIdIt
 }  // DTTTrigOffsetCalibration::analyze

◆ beginRun()

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

override

Definition at line 54 of file DTTTrigOffsetCalibration.cc.

References doTTrigCorrection_, singleTopDQM_cfi::setup, dttriganalyzer_cfi::tTrig, tTrigMap_, and ttrigToken_.

                                                                                      {
   if (doTTrigCorrection_) {
     ESHandle<DTTtrig> tTrig;
     tTrig = setup.getHandle(ttrigToken_);
     tTrigMap_ = &setup.getData(ttrigToken_);
     LogVerbatim("Calibration") << "[DTTTrigOffsetCalibration]: TTrig version: " << tTrig->version() << endl;
   }
 }

◆ bookHistos()

void DTTTrigOffsetCalibration::bookHistos ( DTChamberId chId )

private

Definition at line 184 of file DTTTrigOffsetCalibration.cc.

References combine::histos, LogTrace, nano_mu_digi_cff::sector, DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, AlCaHLTBitMon_QueryRunRegistry::string, theT0SegHistoMap_, to_string(), DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

Referenced by analyze().

                                                           {
   LogTrace("Calibration") << "   Booking histos for Chamber: " << chId;
 
   // Compose the chamber name
   std::string wheel = std::to_string(chId.wheel());
   std::string station = std::to_string(chId.station());
   std::string sector = std::to_string(chId.sector());
 
   string chHistoName = "_W" + wheel + "_St" + station + "_Sec" + sector;
 
   vector<TH1F*> histos;
   // Note the order matters
   histos.push_back(new TH1F(("hRPhiSegT0" + chHistoName).c_str(), "t0 from Phi segments", 500, -60., 60.));
   if (chId.station() != 4)
     histos.push_back(new TH1F(("hRZSegT0" + chHistoName).c_str(), "t0 from Z segments", 500, -60., 60.));
 
   theT0SegHistoMap_[chId] = histos;
 }

◆ endJob()

void DTTTrigOffsetCalibration::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 132 of file DTTTrigOffsetCalibration.cc.

References doTTrigCorrection_, DTTtrig::get(), dttriganalyzer_cfi::kFactor, DTTimeUnits::ns, rootFile_, DTChamberId::station(), theT0SegHistoMap_, dttriganalyzer_cfi::tTrig, tTrigMap_, and DTCalibDBUtils::writeToDB().

Referenced by o2olib.O2ORunMgr::executeJob().

                                       {
   rootFile_->cd();
 
   LogVerbatim("Calibration") << "[DTTTrigOffsetCalibration] Writing histos to file!" << endl;
 
   for (ChamberHistosMap::const_iterator itChHistos = theT0SegHistoMap_.begin(); itChHistos != theT0SegHistoMap_.end();
        ++itChHistos) {
     for (vector<TH1F*>::const_iterator itHist = (*itChHistos).second.begin(); itHist != (*itChHistos).second.end();
          ++itHist)
       (*itHist)->Write();
   }
 
   if (doTTrigCorrection_) {
     // Create the object to be written to DB
     DTTtrig* tTrig = new DTTtrig();
 
     for (ChamberHistosMap::const_iterator itChHistos = theT0SegHistoMap_.begin(); itChHistos != theT0SegHistoMap_.end();
          ++itChHistos) {
       DTChamberId chId = itChHistos->first;
       // Get SuperLayerId's for each ChamberId
       vector<DTSuperLayerId> slIds;
       slIds.push_back(DTSuperLayerId(chId, 1));
       slIds.push_back(DTSuperLayerId(chId, 3));
       if (chId.station() != 4)
         slIds.push_back(DTSuperLayerId(chId, 2));
 
       for (vector<DTSuperLayerId>::const_iterator itSl = slIds.begin(); itSl != slIds.end(); ++itSl) {
         // Get old values from DB
         float ttrigMean = 0;
         float ttrigSigma = 0;
         float kFactor = 0;
         tTrigMap_->get(*itSl, ttrigMean, ttrigSigma, kFactor, DTTimeUnits::ns);
         //FIXME: verify if values make sense
         // Set new values
         float ttrigMeanNew = ttrigMean;
         float ttrigSigmaNew = ttrigSigma;
         float t0SegMean =
             (itSl->superLayer() != 2) ? itChHistos->second[0]->GetMean() : itChHistos->second[1]->GetMean();
 
         float kFactorNew = (kFactor * ttrigSigma + t0SegMean) / ttrigSigma;
 
         tTrig->set(*itSl, ttrigMeanNew, ttrigSigmaNew, kFactorNew, DTTimeUnits::ns);
       }
     }
     LogVerbatim("Calibration") << "[DTTTrigOffsetCalibration] Writing ttrig object to DB!" << endl;
     // Write the object to DB
     string tTrigRecord = "DTTtrigRcd";
     DTCalibDBUtils::writeToDB(tTrigRecord, tTrig);
   }
 }

◆ endRun()

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

inlineoverride

Definition at line 40 of file DTTTrigOffsetCalibration.h.

40 {};

Member Data Documentation

◆ doTTrigCorrection_

const bool DTTTrigOffsetCalibration::doTTrigCorrection_

private

Definition at line 50 of file DTTTrigOffsetCalibration.h.

Referenced by beginRun(), and endJob().

◆ dtGeomToken_

const edm::ESGetToken<DTGeometry, MuonGeometryRecord> DTTTrigOffsetCalibration::dtGeomToken_

private

Definition at line 58 of file DTTTrigOffsetCalibration.h.

Referenced by analyze().

◆ rootFile_

TFile* DTTTrigOffsetCalibration::rootFile_

private

Definition at line 53 of file DTTTrigOffsetCalibration.h.

Referenced by analyze(), DTTTrigOffsetCalibration(), endJob(), and ~DTTTrigOffsetCalibration().

◆ select_

DTSegmentSelector* DTTTrigOffsetCalibration::select_

private

Definition at line 47 of file DTTTrigOffsetCalibration.h.

Referenced by analyze(), and DTTTrigOffsetCalibration().

◆ theCalibChamber_

const std::string DTTTrigOffsetCalibration::theCalibChamber_

private

Definition at line 51 of file DTTTrigOffsetCalibration.h.

Referenced by analyze().

◆ theRecHits4DToken_

const edm::EDGetTokenT<DTRecSegment4DCollection> DTTTrigOffsetCalibration::theRecHits4DToken_

private

Definition at line 49 of file DTTTrigOffsetCalibration.h.

Referenced by analyze().

◆ theT0SegHistoMap_

ChamberHistosMap DTTTrigOffsetCalibration::theT0SegHistoMap_

private

Definition at line 55 of file DTTTrigOffsetCalibration.h.

Referenced by analyze(), bookHistos(), and endJob().

◆ tTrigMap_

const DTTtrig* DTTTrigOffsetCalibration::tTrigMap_

private

Definition at line 54 of file DTTTrigOffsetCalibration.h.

Referenced by beginRun(), and endJob().

◆ ttrigToken_

const edm::ESGetToken<DTTtrig, DTTtrigRcd> DTTTrigOffsetCalibration::ttrigToken_

private

Definition at line 57 of file DTTTrigOffsetCalibration.h.

Referenced by beginRun().

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ ChamberHistosMap

Constructor & Destructor Documentation

◆ DTTTrigOffsetCalibration()

◆ ~DTTTrigOffsetCalibration()

Member Function Documentation

◆ analyze()

◆ beginRun()

◆ bookHistos()

◆ endJob()

◆ endRun()

Member Data Documentation

◆ doTTrigCorrection_

◆ dtGeomToken_

◆ rootFile_

◆ select_

◆ theCalibChamber_

◆ theRecHits4DToken_

◆ theT0SegHistoMap_

◆ tTrigMap_

◆ ttrigToken_