CMS 3D CMS Logo
List of all members | Public Member Functions | Private Types | Private Member Functions | Private Attributes
DTVDriftSegmentCalibration Class Reference

#include <DTVDriftSegmentCalibration.h>

Inheritance diagram for DTVDriftSegmentCalibration:
edm::EDAnalyzer edm::EDConsumerBase

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
 
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::vector< ModuleDescription const * > &modules, 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 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>
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)
 
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 34 of file DTVDriftSegmentCalibration.cc.

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

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

◆ ~DTVDriftSegmentCalibration()

DTVDriftSegmentCalibration::~DTVDriftSegmentCalibration ( )
override

Definition at line 53 of file DTVDriftSegmentCalibration.cc.

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

References rootFile_.

Member Function Documentation

◆ analyze()

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

Implements edm::EDAnalyzer.

Definition at line 58 of file DTVDriftSegmentCalibration.cc.

58  {
59  rootFile_->cd();
60 
61  // Get the DT Geometry
62  ESHandle<DTGeometry> dtGeom;
63  eventSetup.get<MuonGeometryRecord>().get(dtGeom);
64 
65  // Get the rechit collection from the event
66  Handle<DTRecSegment4DCollection> all4DSegments;
67  event.getByLabel(theRecHits4DLabel_, all4DSegments);
68 
69  DTChamberId chosenChamberId;
70  if (theCalibChamber_ != "All") {
71  stringstream linestr;
72  int selWheel, selStation, selSector;
73  linestr << theCalibChamber_;
74  linestr >> selWheel >> selStation >> selSector;
75  chosenChamberId = DTChamberId(selWheel, selStation, selSector);
76  LogVerbatim("Calibration") << " Chosen chamber: " << chosenChamberId << endl;
77  }
78  // Loop over segments by chamber
79  DTRecSegment4DCollection::id_iterator chamberIdIt;
80  for (chamberIdIt = all4DSegments->id_begin(); chamberIdIt != all4DSegments->id_end(); ++chamberIdIt) {
81  // Calibrate just the chosen chamber/s
82  if ((theCalibChamber_ != "All") && ((*chamberIdIt) != chosenChamberId))
83  continue;
84 
85  // Book histos
86  if (theVDriftHistoMapTH1F_.find(*chamberIdIt) == theVDriftHistoMapTH1F_.end()) {
87  LogTrace("Calibration") << " Booking histos for Chamber: " << *chamberIdIt;
88  bookHistos(*chamberIdIt);
89  }
90 
91  // Get the chamber from the setup
92  const DTChamber* chamber = dtGeom->chamber(*chamberIdIt);
93  // Get the range for the corresponding ChamberId
94  DTRecSegment4DCollection::range range = all4DSegments->get((*chamberIdIt));
95  // Loop over the rechits of this DetUnit
96  for (DTRecSegment4DCollection::const_iterator segment = range.first; segment != range.second; ++segment) {
97  LogTrace("Calibration") << "Segment local pos (in chamber RF): " << (*segment).localPosition()
98  << "\nSegment global pos: " << chamber->toGlobal((*segment).localPosition());
99 
100  if (!(*select_)(*segment, event, eventSetup))
101  continue;
102 
103  // Fill v-drift values
104  if ((*segment).hasPhi()) {
105  //if( segment->phiSegment()->ist0Valid() ){
106  double segmentVDrift = segment->phiSegment()->vDrift();
107  if (segmentVDrift != 0.00) {
108  (theVDriftHistoMapTH1F_[*chamberIdIt])[0]->Fill(segmentVDrift);
109  (theVDriftHistoMapTH2F_[*chamberIdIt])[0]->Fill(segment->localPosition().x(), segmentVDrift);
110  (theVDriftHistoMapTH2F_[*chamberIdIt])[1]->Fill(segment->localPosition().y(), segmentVDrift);
111  }
112  }
113  // Probably not meaningful
114  if ((*segment).hasZed()) {
115  //if( segment->zSegment()->ist0Valid() ){
116  double segmentVDrift = segment->zSegment()->vDrift();
117  if (segmentVDrift != 0.00) {
118  (theVDriftHistoMapTH1F_[*chamberIdIt])[1]->Fill(segmentVDrift);
119  }
120  }
121  } // DTRecSegment4DCollection::const_iterator segment
122  } // DTRecSegment4DCollection::id_iterator chamberIdIt
123 } // DTVDriftSegmentCalibration::analyze

References bookHistos(), relativeConstraints::chamber, DTGeometry::chamber(), 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 49 of file DTVDriftSegmentCalibration.cc.

49 { TH1::SetDefaultSumw2(true); }

◆ beginRun()

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

Reimplemented from edm::EDAnalyzer.

Definition at line 51 of file DTVDriftSegmentCalibration.cc.

51 {}

◆ bookHistos()

void DTVDriftSegmentCalibration::bookHistos ( DTChamberId  chId)
private

Definition at line 150 of file DTVDriftSegmentCalibration.cc.

150  {
151  // Compose the chamber name
152  std::string wheel = std::to_string(chId.wheel());
153  std::string station = std::to_string(chId.station());
154  std::string sector = std::to_string(chId.sector());
155 
156  string chHistoName = "_W" + wheel + "_St" + station + "_Sec" + sector;
157 
158  vector<TH1F*> histosTH1F;
159  histosTH1F.push_back(
160  new TH1F(("hRPhiVDriftCorr" + chHistoName).c_str(), "v-drift corr. from Phi segments", 200, -0.4, 0.4));
161  if (chId.station() != 4)
162  histosTH1F.push_back(
163  new TH1F(("hRZVDriftCorr" + chHistoName).c_str(), "v-drift corr. from Z segments", 200, -0.4, 0.4));
164 
165  vector<TH2F*> histosTH2F;
166  histosTH2F.push_back(new TH2F(("hRPhiVDriftCorrVsSegmPosX" + chHistoName).c_str(),
167  "v-drift corr. vs. segment x position",
168  250,
169  -125.,
170  125.,
171  200,
172  -0.4,
173  0.4));
174  histosTH2F.push_back(new TH2F(("hRPhiVDriftCorrVsSegmPosY" + chHistoName).c_str(),
175  "v-drift corr. vs. segment y position",
176  250,
177  -125.,
178  125.,
179  200,
180  -0.4,
181  0.4));
182  //if(chId.station() != 4) ...
183 
184  theVDriftHistoMapTH1F_[chId] = histosTH1F;
185  theVDriftHistoMapTH2F_[chId] = histosTH2F;
186 }

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 125 of file DTVDriftSegmentCalibration.cc.

125  {
126  rootFile_->cd();
127 
128  LogVerbatim("Calibration") << "[DTVDriftSegmentCalibration] Writing histos to file!" << endl;
129 
130  for (ChamberHistosMapTH1F::const_iterator itChHistos = theVDriftHistoMapTH1F_.begin();
131  itChHistos != theVDriftHistoMapTH1F_.end();
132  ++itChHistos) {
133  vector<TH1F*>::const_iterator itHistTH1F = (*itChHistos).second.begin();
134  vector<TH1F*>::const_iterator itHistTH1F_end = (*itChHistos).second.end();
135  for (; itHistTH1F != itHistTH1F_end; ++itHistTH1F)
136  (*itHistTH1F)->Write();
137 
138  vector<TH2F*>::const_iterator itHistTH2F = theVDriftHistoMapTH2F_[(*itChHistos).first].begin();
139  vector<TH2F*>::const_iterator itHistTH2F_end = theVDriftHistoMapTH2F_[(*itChHistos).first].end();
140  for (; itHistTH2F != itHistTH2F_end; ++itHistTH2F)
141  (*itHistTH2F)->Write();
142  }
143 
144  /*if(writeVDriftDB_){
145  // ...
146  }*/
147 }

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

FastTimerService_cff.range
range
Definition: FastTimerService_cff.py:34
DTVDriftSegmentCalibration::rootFile_
TFile * rootFile_
Definition: DTVDriftSegmentCalibration.h:45
DTVDriftSegmentCalibration::bookHistos
void bookHistos(DTChamberId)
Definition: DTVDriftSegmentCalibration.cc:150
DTVDriftSegmentCalibration::theVDriftHistoMapTH2F_
ChamberHistosMapTH2F theVDriftHistoMapTH2F_
Definition: DTVDriftSegmentCalibration.h:47
relativeConstraints.station
station
Definition: relativeConstraints.py:67
DTChamber
Definition: DTChamber.h:24
DTVDriftSegmentCalibration::theCalibChamber_
std::string theCalibChamber_
Definition: DTVDriftSegmentCalibration.h:43
edm::EDConsumerBase::consumesCollector
ConsumesCollector consumesCollector()
Use a ConsumesCollector to gather consumes information from helper functions.
Definition: EDConsumerBase.cc:71
edm::Handle< DTRecSegment4DCollection >
edm::RangeMap::id_iterator
identifier iterator
Definition: RangeMap.h:130
DTGeometry::chamber
const DTChamber * chamber(const DTChamberId &id) const
Return a DTChamber given its id.
Definition: DTGeometry.cc:90
edm::EventSetup::get
T get() const
Definition: EventSetup.h:73
DTVDriftSegmentCalibration::theVDriftHistoMapTH1F_
ChamberHistosMapTH1F theVDriftHistoMapTH1F_
Definition: DTVDriftSegmentCalibration.h:46
edm::ESHandle< DTGeometry >
DTSegmentSelector
Definition: DTSegmentSelector.h:24
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
edm::RangeMap::const_iterator
C::const_iterator const_iterator
constant access iterator type
Definition: RangeMap.h:43
makeMuonMisalignmentScenario.wheel
wheel
Definition: makeMuonMisalignmentScenario.py:319
edm::LogVerbatim
Definition: MessageLogger.h:297
get
#define get
DTVDriftSegmentCalibration::theRecHits4DLabel_
edm::InputTag theRecHits4DLabel_
Definition: DTVDriftSegmentCalibration.h:41
HcalObjRepresent::Fill
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:1053
DTChamberId::sector
int sector() const
Definition: DTChamberId.h:49
edm::RangeMap::range
std::pair< const_iterator, const_iterator > range
iterator range
Definition: RangeMap.h:50
relativeConstraints.chamber
chamber
Definition: relativeConstraints.py:53
CSCSkim_cfi.rootFileName
rootFileName
Definition: CSCSkim_cfi.py:9
DTVDriftSegmentCalibration::select_
DTSegmentSelector * select_
Definition: DTVDriftSegmentCalibration.h:39
DTChamberId
Definition: DTChamberId.h:14
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:671
event
How EventSelector::AcceptEvent() decides whether to accept an event for output otherwise it is excluding the probing of A single or multiple positive and the trigger will pass if any such matching triggers are PASS or EXCEPTION[A criterion thatmatches no triggers at all is detected and causes a throw.] A single negative with an expectation of appropriate bit checking in the decision and the trigger will pass if any such matching triggers are FAIL or EXCEPTION A wildcarded negative criterion that matches more than one trigger in the trigger but the state exists so we define the behavior If all triggers are the negative crieriion will lead to accepting the event(this again matches the behavior of "!*" before the partial wildcard feature was incorporated). The per-event "cost" of each negative criterion with multiple relevant triggers is about the same as ! *was in the past
MuonGeometryRecord
Definition: MuonGeometryRecord.h:34
DTChamberId::wheel
int wheel() const
Return the wheel number.
Definition: DTChamberId.h:39
edm::InputTag
Definition: InputTag.h:15
edm::ConsumesCollector
Definition: ConsumesCollector.h:39
DTChamberId::station
int station() const
Return the station number.
Definition: DTChamberId.h:42
muonDTDigis_cfi.pset
pset
Definition: muonDTDigis_cfi.py:27
Generated for CMSSW Reference Manual by doxygen 1.8.16