CMS 3D CMS Logo

List of all members | Public Member Functions | Protected Member Functions | Private Member Functions | Private Attributes
DTChamberEfficiency Class Reference

#include <DTChamberEfficiency.h>

Inheritance diagram for DTChamberEfficiency:
DQMEDAnalyzer edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, edm::Accumulator > edm::one::EDProducerBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

void analyze (const edm::Event &event, const edm::EventSetup &eventSetup) override
 
void dqmBeginRun (const edm::Run &, const edm::EventSetup &) override
 
 DTChamberEfficiency (const edm::ParameterSet &pset)
 
 ~DTChamberEfficiency () override
 
- Public Member Functions inherited from DQMEDAnalyzer
void accumulate (edm::Event const &event, edm::EventSetup const &setup) final
 
virtual void analyze (edm::Event const &, edm::EventSetup const &)
 
void beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
 
void beginRun (edm::Run const &run, edm::EventSetup const &setup) final
 
virtual void dqmBeginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
 
virtual void dqmBeginRun (edm::Run const &, edm::EventSetup const &)
 
 DQMEDAnalyzer ()
 
virtual void dqmEndLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
 
virtual void dqmEndRun (edm::Run const &, edm::EventSetup const &)
 
void endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final
 
void endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final
 
void endRun (edm::Run const &, edm::EventSetup const &) final
 
void endRunProduce (edm::Run &run, edm::EventSetup const &setup) final
 
virtual bool getCanSaveByLumi ()
 
- Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, edm::Accumulator >
 EDProducer ()=default
 
SerialTaskQueueglobalLuminosityBlocksQueue () final
 
SerialTaskQueueglobalRunsQueue () final
 
bool hasAbilityToProduceInBeginLumis () const final
 
bool hasAbilityToProduceInBeginRuns () const final
 
bool hasAbilityToProduceInEndLumis () const final
 
bool hasAbilityToProduceInEndRuns () const final
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
- Public Member Functions inherited from edm::one::EDProducerBase
 EDProducerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const
 
bool wantsStreamRuns () const
 
 ~EDProducerBase () override
 
- Public Member Functions inherited from edm::ProducerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
std::vector< edm::ProductResolverIndex > const & indiciesForPutProducts (BranchType iBranchType) const
 
 ProducerBase ()
 
std::vector< edm::ProductResolverIndex > const & putTokenIndexToProductResolverIndex () const
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription const &)> registrationCallback () const
 used by the fwk to register list of products More...
 
void resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
 
 ~ProducerBase () noexcept(false) override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () 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
 
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 const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBaseoperator= (EDConsumerBase &&)=default
 
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)
 

Protected Member Functions

void bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
 
- Protected Member Functions inherited from edm::ProducerBase
ProducesCollector producesCollector ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
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<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)
 

Private Member Functions

bool chamberSelection (const DetId &idDetLay, reco::TransientTrack &trans_track) const
 
std::vector< const DetLayer * > compatibleLayers (const NavigationSchool &navigationSchool, const DetLayer *initialLayer, const FreeTrajectoryState &fts, PropagationDirection propDir)
 
edm::ESHandle< Propagatorpropagator () const
 
MeasurementContainer segQualityCut (const MeasurementContainer &seg_list) const
 

Private Attributes

bool debug
 
edm::ESHandle< DTGeometrydtGeom
 
std::vector< std::vector< MonitorElement * > > histosPerW
 
edm::InputTag labelRPCRecHits
 
edm::ESHandle< MagneticFieldmagfield
 
edm::InputTag thecscSegments
 
edm::InputTag thedt4DSegments
 
Chi2MeasurementEstimatortheEstimator
 
double theMaxChi2
 
MuonDetLayerMeasurementstheMeasurementExtractor
 
int theMinNrec
 
std::string theNavigationType
 
double theNSigma
 
MuonServiceProxytheService
 
edm::ESHandle< GlobalTrackingGeometrytheTrackingGeometry
 
edm::InputTag theTracksLabel_
 
edm::EDGetTokenT< reco::TrackCollectiontheTracksToken_
 

Additional Inherited Members

- Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore DQMStore
 
typedef dqm::reco::MonitorElement MonitorElement
 
- Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase ModuleType
 
- Public Types inherited from edm::ProducerBase
using ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
 
typedef ProductRegistryHelper::TypeLabelList TypeLabelList
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMTokenlumiToken_
 
edm::EDPutTokenT< DQMTokenrunToken_
 

Detailed Description

Description:

This class provides the histograms for the calculation of the efficiency of muons reconstruction in the DTs. It is applicable both in presence or absence of a magnetic field. Histos are 2D Sector vs Chamber plots for each wheel

Author
: Mario Pelliccioni - INFN Torino pelli.nosp@m.cci@.nosp@m.cern..nosp@m.ch $date : 05/12/2008 16:51:04 CET $

Modification:

Definition at line 55 of file DTChamberEfficiency.h.

Constructor & Destructor Documentation

DTChamberEfficiency::DTChamberEfficiency ( const edm::ParameterSet pset)

Definition at line 53 of file DTChamberEfficiency.cc.

References Chi2MeasurementEstimator_cfi::Chi2MeasurementEstimator, debug, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HLT_2018_cff::InputTag, MuonServiceProxy_cff::MuonServiceProxy, dtChamberEfficiency_cfi::theMaxChi2, dtChamberEfficiency_cfi::theMinNrec, and createJobs::theNSigma.

53  {
54  // Get the debug parameter for verbose output
55  debug = pSet.getUntrackedParameter<bool>("debug", false);
56 
57  LogVerbatim("DTDQM|DTMonitorModule|DTChamberEfficiency") << "DTChamberEfficiency: constructor called";
58 
59  // service parameters
60  ParameterSet serviceParameters = pSet.getParameter<ParameterSet>("ServiceParameters");
61  theService = new MuonServiceProxy(serviceParameters);
62 
63  theTracksLabel_ = pSet.getParameter<InputTag>("TrackCollection");
64  theTracksToken_ = consumes<reco::TrackCollection>(theTracksLabel_);
65 
66  theMaxChi2 = static_cast<unsigned int>(pSet.getParameter<double>("theMaxChi2"));
67  theNSigma = pSet.getParameter<double>("theNSigma");
68  theMinNrec = static_cast<int>(pSet.getParameter<double>("theMinNrec"));
69 
70  labelRPCRecHits = pSet.getParameter<InputTag>("theRPCRecHits");
71 
72  thedt4DSegments = pSet.getParameter<InputTag>("dt4DSegments");
73  thecscSegments = pSet.getParameter<InputTag>("cscSegments");
74 
76 
78  thedt4DSegments, thecscSegments, labelRPCRecHits, InputTag(), InputTag(), iC, true, false, false, false);
79 
80  theNavigationType = pSet.getParameter<string>("NavigationType");
81 
83 }
MuonDetLayerMeasurements * theMeasurementExtractor
T getParameter(std::string const &) const
MuonServiceProxy * theService
edm::InputTag thecscSegments
edm::InputTag labelRPCRecHits
edm::EDGetTokenT< reco::TrackCollection > theTracksToken_
edm::InputTag thedt4DSegments
ConsumesCollector consumesCollector()
Use a ConsumesCollector to gather consumes information from helper functions.
edm::InputTag theTracksLabel_
Chi2MeasurementEstimator * theEstimator
DTChamberEfficiency::~DTChamberEfficiency ( )
override

Definition at line 85 of file DTChamberEfficiency.cc.

References LogTrace.

85  {
86  LogTrace("DTDQM|DTMonitorModule|DTChamberEfficiency") << "DTChamberEfficiency: destructor called";
87 
88  // free memory
89  delete theService;
91  delete theEstimator;
92 }
MuonDetLayerMeasurements * theMeasurementExtractor
MuonServiceProxy * theService
#define LogTrace(id)
Chi2MeasurementEstimator * theEstimator

Member Function Documentation

void DTChamberEfficiency::analyze ( const edm::Event event,
const edm::EventSetup eventSetup 
)
override

Definition at line 133 of file DTChamberEfficiency.cc.

References alongMomentum, debug, DetId::det(), MuonSubdetId::DT, event(), HcalObjRepresent::Fill(), TrajectoryStateOnSurface::freeState(), combine::histos, mps_fire::i, createfilelist::int, TrajectoryStateOnSurface::isValid(), edm::HandleBase::isValid(), LogTrace, volumeBasedMagneticField_1103l_cfi::magfield, DetId::Muon, oppositeToMomentum, TrackCandidateProducer_cfi::propagator, DetId::rawId(), reco::TransientTrack::recHitsSize(), DTChamberId::sector(), DTChamberId::station(), DetId::subdetId(), TrajectoryStateOnSurface::surface(), dtChamberEfficiency_cfi::theMinNrec, HLT_2018_cff::track, PDWG_EXOHSCP_cff::tracks, and DTChamberId::wheel().

133  {
134  LogTrace("DTDQM|DTMonitorModule|DTChamberEfficiency")
135  << "--- [DTChamberEfficiency] Event analysed #Run: " << event.id().run() << " #Event: " << event.id().event()
136  << endl;
137 
138  theService->update(eventSetup);
140 
141  //Read tracks from event
143  event.getByToken(theTracksToken_, tracks);
144 
145  if (tracks.isValid()) { // check the validity of the collection
146 
147  //loop over the muons
148  for (reco::TrackCollection::const_iterator track = tracks->begin(); track != tracks->end(); ++track) {
150  const int recHitsize = (int)trans_track.recHitsSize();
151  if (recHitsize < theMinNrec)
152  continue;
153 
154  // printout the DT rechits used by the track
155  if (debug) {
156  LogTrace("DTDQM|DTMonitorModule|DTChamberEfficiency") << "--- New track" << endl;
157  set<DTChamberId> chAlrUsed;
158  for (trackingRecHit_iterator rHit = trans_track.recHitsBegin(); rHit != trans_track.recHitsEnd(); ++rHit) {
159  DetId rHitid = (*rHit)->geographicalId();
160  if (!(rHitid.det() == DetId::Muon && rHitid.subdetId() == MuonSubdetId::DT))
161  continue;
162  DTChamberId wId(rHitid.rawId());
163  if (chAlrUsed.find(wId) != chAlrUsed.end())
164  continue;
165  chAlrUsed.insert(wId);
166  LogTrace("DTDQM|DTMonitorModule|DTChamberEfficiency") << " " << wId << endl;
167  }
168  }
169 
170  // Get the layer on which the seed relies
171  DetId id = trans_track.track().innerDetId();
172  const DetLayer* initialLayer = theService->detLayerGeometry()->idToLayer(id);
173 
174  TrajectoryStateOnSurface init_fs = trans_track.innermostMeasurementState();
175  const FreeTrajectoryState* init_fs_free = init_fs.freeState();
176 
177  //get the list of compatible layers
178  vector<const DetLayer*> layer_list =
179  compatibleLayers(*theService->muonNavigationSchool(), initialLayer, *init_fs_free, alongMomentum);
180  vector<const DetLayer*> layer_list_2 =
181  compatibleLayers(*theService->muonNavigationSchool(), initialLayer, *init_fs_free, oppositeToMomentum);
182 
183  layer_list.insert(layer_list.end(), layer_list_2.begin(), layer_list_2.end());
184 
185  set<DTChamberId> alreadyCheckedCh;
186 
187  //loop over the list of compatible layers
188  for (int i = 0; i < (int)layer_list.size(); i++) {
189  //propagate the track to the i-th layer
190  TrajectoryStateOnSurface tsos = propagator()->propagate(init_fs, layer_list.at(i)->surface());
191  if (!tsos.isValid())
192  continue;
193 
194  //determine the chambers kinematically compatible with the track on the i-th layer
195  vector<DetWithState> dss = layer_list.at(i)->compatibleDets(tsos, *propagator(), *theEstimator);
196 
197  if (dss.empty())
198  continue;
199 
200  // get the first det (it's the most compatible)
201  const DetWithState detWithState = dss.front();
202  const DetId idDetLay = detWithState.first->geographicalId();
203 
204  // check if this is a DT and the track has the needed quality
205  if (!chamberSelection(idDetLay, trans_track))
206  continue;
207 
208  DTChamberId DTid = (DTChamberId)idDetLay;
209 
210  // check if the chamber has already been counted
211  if (alreadyCheckedCh.find(DTid) != alreadyCheckedCh.end())
212  continue;
213  alreadyCheckedCh.insert(DTid);
214 
215  // get the compatible measurements
216  MeasurementContainer detMeasurements_initial = theMeasurementExtractor->measurements(
217  layer_list.at(i), detWithState.first, detWithState.second, *theEstimator, event);
218  LogTrace("DTDQM|DTMonitorModule|DTChamberEfficiency")
219  << " chamber: " << DTid << " has: " << detMeasurements_initial.size() << " comp. meas." << endl;
220 
221  //we want to be more picky about the quality of the segments:
222  //exclude the segments with less than 12 hits
223  MeasurementContainer detMeasurements = segQualityCut(detMeasurements_initial);
224 
225  // get the histos for this chamber
226  vector<MonitorElement*> histos = histosPerW[DTid.wheel() + 2];
227  // fill them
228  if (!detMeasurements_initial.empty())
229  histos[0]->Fill(DTid.sector(), DTid.station(), 1.);
230  if (!detMeasurements.empty())
231  histos[1]->Fill(DTid.sector(), DTid.station(), 1.);
232  histos[2]->Fill(DTid.sector(), DTid.station(), 1.);
233  }
234  }
235  } else {
236  LogInfo("DTDQM|DTMonitorModule|DTChamberEfficiency")
237  << "[DTChamberEfficiency] Collection: " << theTracksLabel_ << " is not valid!" << endl;
238  }
239  return;
240 }
MuonDetLayerMeasurements * theMeasurementExtractor
edm::ESHandle< GlobalTrackingGeometry > theTrackingGeometry
MuonServiceProxy * theService
MeasurementContainer segQualityCut(const MeasurementContainer &seg_list) const
std::vector< std::vector< MonitorElement * > > histosPerW
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
size_t recHitsSize() const
number of RecHits
edm::EDGetTokenT< reco::TrackCollection > theTracksToken_
const SurfaceType & surface() const
bool chamberSelection(const DetId &idDetLay, reco::TransientTrack &trans_track) const
edm::ESHandle< Propagator > propagator() const
std::pair< const GeomDet *, TrajectoryStateOnSurface > DetWithState
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
FreeTrajectoryState const * freeState(bool withErrors=true) const
constexpr int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:48
std::vector< const DetLayer * > compatibleLayers(const NavigationSchool &navigationSchool, const DetLayer *initialLayer, const FreeTrajectoryState &fts, PropagationDirection propDir)
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
bool isValid() const
Definition: HandleBase.h:70
MeasurementContainer measurements(const DetLayer *layer, const GeomDet *det, const TrajectoryStateOnSurface &stateOnDet, const MeasurementEstimator &est, const edm::Event &iEvent)
edm::ESHandle< MagneticField > magfield
#define LogTrace(id)
Definition: DetId.h:17
edm::InputTag theTracksLabel_
std::vector< TrajectoryMeasurement > MeasurementContainer
TrajectoryStateOnSurface propagate(STA const &state, SUR const &surface) const
Definition: Propagator.h:50
histos
Definition: combine.py:4
void setEvent(const edm::Event &)
set event
int sector() const
Definition: DTChamberId.h:49
Chi2MeasurementEstimator * theEstimator
static constexpr int DT
Definition: MuonSubdetId.h:11
int station() const
Return the station number.
Definition: DTChamberId.h:42
int wheel() const
Return the wheel number.
Definition: DTChamberId.h:39
T const * product() const
Definition: ESHandle.h:86
constexpr Detector det() const
get the detector field from this detid
Definition: DetId.h:46
void DTChamberEfficiency::bookHistograms ( DQMStore::IBooker ibooker,
edm::Run const &  iRun,
edm::EventSetup const &  context 
)
overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 104 of file DTChamberEfficiency.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book2D(), combine::histos, LogTrace, dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), and makeMuonMisalignmentScenario::wheel.

106  {
107  LogTrace("DTDQM|DTMonitorModule|DTChamberEfficiency") << "DTChamberEfficiency: booking histos";
108 
109  // Create the monitor elements
110  ibooker.setCurrentFolder("DT/05-ChamberEff/Task");
111 
112  for (int wheel = -2; wheel <= 2; wheel++) {
113  vector<MonitorElement*> histos;
114 
115  stringstream wheel_str;
116  wheel_str << wheel;
117 
118  histos.push_back(ibooker.book2D(
119  "hCountSectVsChamb_All_W" + wheel_str.str(), "Countings for wheel " + wheel_str.str(), 14, 1., 15., 4, 1., 5.));
120 
121  histos.push_back(ibooker.book2D(
122  "hCountSectVsChamb_Qual_W" + wheel_str.str(), "Countings for wheel " + wheel_str.str(), 14, 1., 15., 4, 1., 5.));
123 
124  histos.push_back(ibooker.book2D(
125  "hExtrapSectVsChamb_W" + wheel_str.str(), "Extrapolations for wheel " + wheel_str.str(), 14, 1., 15., 4, 1., 5.));
126 
127  histosPerW.push_back(histos);
128  }
129 
130  return;
131 }
std::vector< std::vector< MonitorElement * > > histosPerW
#define LogTrace(id)
histos
Definition: combine.py:4
bool DTChamberEfficiency::chamberSelection ( const DetId idDetLay,
reco::TransientTrack trans_track 
) const
private

Definition at line 242 of file DTChamberEfficiency.cc.

References DetId::det(), MuonSubdetId::DT, DetId::Muon, reco::TransientTrack::recHit(), reco::TransientTrack::recHitsSize(), and DetId::subdetId().

242  {
243  //check that we have a muon and that is a DT detector
244  if (!(idDetLay.det() == DetId::Muon && idDetLay.subdetId() == MuonSubdetId::DT))
245  return false;
246 
247  if (trans_track.recHitsSize() == 2)
248  if (trans_track.recHit(0)->geographicalId() == idDetLay || trans_track.recHit(1)->geographicalId() == idDetLay)
249  return false;
250 
251  return true;
252 }
size_t recHitsSize() const
number of RecHits
constexpr int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:48
TrackingRecHitRef recHit(size_t i) const
get n-th recHit
static constexpr int DT
Definition: MuonSubdetId.h:11
constexpr Detector det() const
get the detector field from this detid
Definition: DetId.h:46
vector< const DetLayer * > DTChamberEfficiency::compatibleLayers ( const NavigationSchool navigationSchool,
const DetLayer initialLayer,
const FreeTrajectoryState fts,
PropagationDirection  propDir 
)
private

Definition at line 281 of file DTChamberEfficiency.cc.

References DirectMuonNavigation::compatibleLayers(), and NavigationSchool::compatibleLayers().

284  {
285  vector<const DetLayer*> detLayers;
286 
287  if (theNavigationType == "Standard") {
288  // ask for compatible layers
289  detLayers = navigationSchool.compatibleLayers(*initialLayer, fts, propDir);
290  // I have to fit by hand the first layer until the seedTSOS is defined on the first rechit layer
291  // In fact the first layer is not returned by initialLayer->compatibleLayers.
292 
293  detLayers.insert(detLayers.begin(), initialLayer);
294 
295  } else if (theNavigationType == "Direct") {
296  DirectMuonNavigation navigation(ESHandle<MuonDetLayerGeometry>(&*theService->detLayerGeometry()));
297  detLayers = navigation.compatibleLayers(fts, propDir);
298  } else
299  LogError("DTDQM|DTMonitorModule|DTChamberEfficiency") << "No Properly Navigation Selected!!" << endl;
300 
301  return detLayers;
302 }
MuonServiceProxy * theService
std::vector< const DetLayer * > compatibleLayers(const DetLayer &detLayer, Args &&...args) const
Returns all layers compatible.
void DTChamberEfficiency::dqmBeginRun ( const edm::Run run,
const edm::EventSetup setup 
)
override

Definition at line 94 of file DTChamberEfficiency.cc.

References edm::EventSetup::get(), and volumeBasedMagneticField_1103l_cfi::magfield.

94  {
95  // Get the DT Geometry
96  setup.get<MuonGeometryRecord>().get(dtGeom);
97 
100 
101  return;
102 }
edm::ESHandle< GlobalTrackingGeometry > theTrackingGeometry
edm::ESHandle< DTGeometry > dtGeom
edm::ESHandle< MagneticField > magfield
T get() const
Definition: EventSetup.h:73
ESHandle< Propagator > DTChamberEfficiency::propagator ( ) const
inlineprivate

Definition at line 304 of file DTChamberEfficiency.cc.

304  {
305  return theService->propagator("SteppingHelixPropagatorAny");
306 }
MuonServiceProxy * theService
MeasurementContainer DTChamberEfficiency::segQualityCut ( const MeasurementContainer seg_list) const
private

Definition at line 254 of file DTChamberEfficiency.cc.

References createfilelist::int, mps_fire::result, and DTChamberId::station().

254  {
256 
257  for (MeasurementContainer::const_iterator mescont_Itr = seg_list.begin(); mescont_Itr != seg_list.end();
258  ++mescont_Itr) {
259  //get the rechits of the segment
260  TransientTrackingRecHit::ConstRecHitContainer recHit_list = mescont_Itr->recHit()->transientHits();
261 
262  //loop over the rechits and get the number of hits
263  int nhit_seg(0);
264  for (TransientTrackingRecHit::ConstRecHitContainer::const_iterator recList_Itr = recHit_list.begin();
265  recList_Itr != recHit_list.end();
266  ++recList_Itr) {
267  nhit_seg += (int)(*recList_Itr)->transientHits().size();
268  }
269 
270  DTChamberId tmpId = (DTChamberId)mescont_Itr->recHit()->hit()->geographicalId();
271 
272  if (tmpId.station() < 4 && nhit_seg >= 12)
273  result.push_back(*mescont_Itr);
274  if (tmpId.station() == 4 && nhit_seg >= 8)
275  result.push_back(*mescont_Itr);
276  }
277 
278  return result;
279 }
std::vector< ConstRecHitPointer > ConstRecHitContainer
std::vector< TrajectoryMeasurement > MeasurementContainer
int station() const
Return the station number.
Definition: DTChamberId.h:42

Member Data Documentation

bool DTChamberEfficiency::debug
private
edm::ESHandle<DTGeometry> DTChamberEfficiency::dtGeom
private

Definition at line 98 of file DTChamberEfficiency.h.

std::vector<std::vector<MonitorElement*> > DTChamberEfficiency::histosPerW
private

Definition at line 107 of file DTChamberEfficiency.h.

edm::InputTag DTChamberEfficiency::labelRPCRecHits
private

Definition at line 88 of file DTChamberEfficiency.h.

edm::ESHandle<MagneticField> DTChamberEfficiency::magfield
private

Definition at line 104 of file DTChamberEfficiency.h.

edm::InputTag DTChamberEfficiency::thecscSegments
private

Definition at line 90 of file DTChamberEfficiency.h.

edm::InputTag DTChamberEfficiency::thedt4DSegments
private

Definition at line 89 of file DTChamberEfficiency.h.

Chi2MeasurementEstimator* DTChamberEfficiency::theEstimator
private

Definition at line 102 of file DTChamberEfficiency.h.

double DTChamberEfficiency::theMaxChi2
private

Definition at line 92 of file DTChamberEfficiency.h.

MuonDetLayerMeasurements* DTChamberEfficiency::theMeasurementExtractor
private

Definition at line 101 of file DTChamberEfficiency.h.

int DTChamberEfficiency::theMinNrec
private

Definition at line 94 of file DTChamberEfficiency.h.

std::string DTChamberEfficiency::theNavigationType
private

Definition at line 96 of file DTChamberEfficiency.h.

double DTChamberEfficiency::theNSigma
private

Definition at line 93 of file DTChamberEfficiency.h.

MuonServiceProxy* DTChamberEfficiency::theService
private

Definition at line 100 of file DTChamberEfficiency.h.

edm::ESHandle<GlobalTrackingGeometry> DTChamberEfficiency::theTrackingGeometry
private

Definition at line 105 of file DTChamberEfficiency.h.

edm::InputTag DTChamberEfficiency::theTracksLabel_
private

Definition at line 85 of file DTChamberEfficiency.h.

edm::EDGetTokenT<reco::TrackCollection> DTChamberEfficiency::theTracksToken_
private

Definition at line 86 of file DTChamberEfficiency.h.