DTVDriftSegmentCalibration Class Reference

#include <DTVDriftSegmentCalibration.h>

Inheritance diagram for DTVDriftSegmentCalibration:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &eventSetup) override
void	beginJob () override
void	beginRun (const edm::Run &run, const edm::EventSetup &setup) override
	DTVDriftSegmentCalibration (const edm::ParameterSet &pset)
void	endJob () override
	~DTVDriftSegmentCalibration () override
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
	~EDAnalyzer () 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 &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
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 &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Types
typedef std::map< DTChamberId, std::vector< TH1F * > >	ChamberHistosMapTH1F
typedef std::map< DTChamberId, std::vector< TH2F * > >	ChamberHistosMapTH2F

Private Member Functions
void	bookHistos (DTChamberId)

Private Attributes
TFile *	rootFile_
DTSegmentSelector *	select_
std::string	theCalibChamber_
edm::InputTag	theRecHits4DLabel_
ChamberHistosMapTH1F	theVDriftHistoMapTH1F_
ChamberHistosMapTH2F	theVDriftHistoMapTH2F_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsInputProcessBlocks ()
static bool	wantsProcessBlocks ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
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...
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Produces histograms from v-drift computation in segment fit to be used for v-drift calibration

Author

A. Vilela Pereira

Definition at line 22 of file DTVDriftSegmentCalibration.h.

Member Typedef Documentation

ChamberHistosMapTH1F

typedef std::map<DTChamberId, std::vector<TH1F*> > DTVDriftSegmentCalibration::ChamberHistosMapTH1F

private

Definition at line 35 of file DTVDriftSegmentCalibration.h.

ChamberHistosMapTH2F

typedef std::map<DTChamberId, std::vector<TH2F*> > DTVDriftSegmentCalibration::ChamberHistosMapTH2F

private

Definition at line 36 of file DTVDriftSegmentCalibration.h.

Constructor & Destructor Documentation

DTVDriftSegmentCalibration()

DTVDriftSegmentCalibration::DTVDriftSegmentCalibration

(

const edm::ParameterSet &

pset

)

Definition at line 33 of file DTVDriftSegmentCalibration.cc.

    : theRecHits4DLabel_(pset.getParameter<InputTag>("recHits4DLabel")),
      //writeVDriftDB_(pset.getUntrackedParameter<bool>("writeVDriftDB", false)),
      theCalibChamber_(pset.getUntrackedParameter<string>("calibChamber", "All")) {
  LogVerbatim("Calibration") << "[DTVDriftSegmentCalibration] Constructor called!";
 
  edm::ConsumesCollector collector(consumesCollector());
  select_ = new DTSegmentSelector(pset, collector);
  consumes<DTRecSegment4DCollection>(edm::InputTag(theRecHits4DLabel_));
  // the root file which will contain the histos
  string rootFileName = pset.getUntrackedParameter<string>("rootFileName", "DTVDriftHistos.root");
  rootFile_ = new TFile(rootFileName.c_str(), "RECREATE");
  rootFile_->cd();
}

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

~DTVDriftSegmentCalibration()

DTVDriftSegmentCalibration::~DTVDriftSegmentCalibration ( )

override

Definition at line 52 of file DTVDriftSegmentCalibration.cc.

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

References rootFile_.

Member Function Documentation

analyze()

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

overridevirtual

Implements edm::EDAnalyzer.

Definition at line 57 of file DTVDriftSegmentCalibration.cc.

                                                                                         {
  rootFile_->cd();
 
  // Get the DT Geometry
  ESHandle<DTGeometry> dtGeom;
  eventSetup.get<MuonGeometryRecord>().get(dtGeom);
 
  // Get the rechit collection from the event
  Handle<DTRecSegment4DCollection> all4DSegments;
  event.getByLabel(theRecHits4DLabel_, all4DSegments);
 
  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;
  }
  // Loop over segments by chamber
  DTRecSegment4DCollection::id_iterator chamberIdIt;
  for (chamberIdIt = all4DSegments->id_begin(); chamberIdIt != all4DSegments->id_end(); ++chamberIdIt) {
    // Calibrate just the chosen chamber/s
    if ((theCalibChamber_ != "All") && ((*chamberIdIt) != chosenChamberId))
      continue;
 
    // Book histos
    if (theVDriftHistoMapTH1F_.find(*chamberIdIt) == theVDriftHistoMapTH1F_.end()) {
      LogTrace("Calibration") << "   Booking histos for Chamber: " << *chamberIdIt;
      bookHistos(*chamberIdIt);
    }
 
    // Get the chamber from the setup
    const DTChamber* chamber = dtGeom->chamber(*chamberIdIt);
    // 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 v-drift values
      if ((*segment).hasPhi()) {
        //if( segment->phiSegment()->ist0Valid() ){
        double segmentVDrift = segment->phiSegment()->vDrift();
        if (segmentVDrift != 0.00) {
          (theVDriftHistoMapTH1F_[*chamberIdIt])[0]->Fill(segmentVDrift);
          (theVDriftHistoMapTH2F_[*chamberIdIt])[0]->Fill(segment->localPosition().x(), segmentVDrift);
          (theVDriftHistoMapTH2F_[*chamberIdIt])[1]->Fill(segment->localPosition().y(), segmentVDrift);
        }
      }
      // Probably not meaningful
      if ((*segment).hasZed()) {
        //if( segment->zSegment()->ist0Valid() ){
        double segmentVDrift = segment->zSegment()->vDrift();
        if (segmentVDrift != 0.00) {
          (theVDriftHistoMapTH1F_[*chamberIdIt])[1]->Fill(segmentVDrift);
        }
      }
    }  // DTRecSegment4DCollection::const_iterator segment
  }    // DTRecSegment4DCollection::id_iterator chamberIdIt
}  // DTVDriftSegmentCalibration::analyze

References bookHistos(), relativeConstraints::chamber, DTGeometry::chamber(), edmPickEvents::event, HcalObjRepresent::Fill(), edm::EventSetup::get(), get, LogTrace, FastTimerService_cff::range, rootFile_, select_, theCalibChamber_, theRecHits4DLabel_, theVDriftHistoMapTH1F_, and theVDriftHistoMapTH2F_.

beginJob()

void DTVDriftSegmentCalibration::beginJob ( void  )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 48 of file DTVDriftSegmentCalibration.cc.

{ TH1::SetDefaultSumw2(true); }

beginRun()

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

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 50 of file DTVDriftSegmentCalibration.cc.

{}

bookHistos()

void DTVDriftSegmentCalibration::bookHistos ( DTChamberId  chId )

private

Definition at line 149 of file DTVDriftSegmentCalibration.cc.

                                                            {
  // 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*> histosTH1F;
  histosTH1F.push_back(
      new TH1F(("hRPhiVDriftCorr" + chHistoName).c_str(), "v-drift corr. from Phi segments", 200, -0.4, 0.4));
  if (chId.station() != 4)
    histosTH1F.push_back(
        new TH1F(("hRZVDriftCorr" + chHistoName).c_str(), "v-drift corr. from Z segments", 200, -0.4, 0.4));
 
  vector<TH2F*> histosTH2F;
  histosTH2F.push_back(new TH2F(("hRPhiVDriftCorrVsSegmPosX" + chHistoName).c_str(),
                                "v-drift corr. vs. segment x position",
                                250,
                                -125.,
                                125.,
                                200,
                                -0.4,
                                0.4));
  histosTH2F.push_back(new TH2F(("hRPhiVDriftCorrVsSegmPosY" + chHistoName).c_str(),
                                "v-drift corr. vs. segment y position",
                                250,
                                -125.,
                                125.,
                                200,
                                -0.4,
                                0.4));
  //if(chId.station() != 4) ...
 
  theVDriftHistoMapTH1F_[chId] = histosTH1F;
  theVDriftHistoMapTH2F_[chId] = histosTH2F;
}

References DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, AlCaHLTBitMon_QueryRunRegistry::string, theVDriftHistoMapTH1F_, theVDriftHistoMapTH2F_, DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

Referenced by analyze().

endJob()

void DTVDriftSegmentCalibration::endJob ( void  )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 124 of file DTVDriftSegmentCalibration.cc.

                                        {
  rootFile_->cd();
 
  LogVerbatim("Calibration") << "[DTVDriftSegmentCalibration] Writing histos to file!" << endl;
 
  for (ChamberHistosMapTH1F::const_iterator itChHistos = theVDriftHistoMapTH1F_.begin();
       itChHistos != theVDriftHistoMapTH1F_.end();
       ++itChHistos) {
    vector<TH1F*>::const_iterator itHistTH1F = (*itChHistos).second.begin();
    vector<TH1F*>::const_iterator itHistTH1F_end = (*itChHistos).second.end();
    for (; itHistTH1F != itHistTH1F_end; ++itHistTH1F)
      (*itHistTH1F)->Write();
 
    vector<TH2F*>::const_iterator itHistTH2F = theVDriftHistoMapTH2F_[(*itChHistos).first].begin();
    vector<TH2F*>::const_iterator itHistTH2F_end = theVDriftHistoMapTH2F_[(*itChHistos).first].end();
    for (; itHistTH2F != itHistTH2F_end; ++itHistTH2F)
      (*itHistTH2F)->Write();
  }
 
  /*if(writeVDriftDB_){
     // ...
  }*/
}

References rootFile_, theVDriftHistoMapTH1F_, and theVDriftHistoMapTH2F_.

Referenced by o2olib.O2ORunMgr::executeJob().

Member Data Documentation

rootFile_

TFile* DTVDriftSegmentCalibration::rootFile_

private

Definition at line 45 of file DTVDriftSegmentCalibration.h.

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

select_

DTSegmentSelector* DTVDriftSegmentCalibration::select_

private

Definition at line 39 of file DTVDriftSegmentCalibration.h.

Referenced by analyze(), and DTVDriftSegmentCalibration().

theCalibChamber_

std::string DTVDriftSegmentCalibration::theCalibChamber_

private

Definition at line 43 of file DTVDriftSegmentCalibration.h.

Referenced by analyze().

theRecHits4DLabel_

edm::InputTag DTVDriftSegmentCalibration::theRecHits4DLabel_

private

Definition at line 41 of file DTVDriftSegmentCalibration.h.

Referenced by analyze(), and DTVDriftSegmentCalibration().

theVDriftHistoMapTH1F_

ChamberHistosMapTH1F DTVDriftSegmentCalibration::theVDriftHistoMapTH1F_

private

Definition at line 46 of file DTVDriftSegmentCalibration.h.

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

theVDriftHistoMapTH2F_

ChamberHistosMapTH2F DTVDriftSegmentCalibration::theVDriftHistoMapTH2F_

private

Definition at line 47 of file DTVDriftSegmentCalibration.h.

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