#include <DTChamberEfficiency.h>

Inheritance diagram for DTChamberEfficiency:

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

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

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

void	beginStream (edm::StreamID id) final

	DQMEDAnalyzer ()

void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

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

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Protected Member Functions
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

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< Propagator >	propagator () const

MeasurementContainer	segQualityCut (const MeasurementContainer &seg_list) const

Private Attributes
bool	debug

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

edm::InputTag	labelRPCRecHits

edm::InputTag	thecscSegments

edm::InputTag	thedt4DSegments

Chi2MeasurementEstimator *	theEstimator

double	theMaxChi2

MuonDetLayerMeasurements *	theMeasurementExtractor

int	theMinNrec

std::string	theNavigationType

double	theNSigma

MuonServiceProxy *	theService

edm::InputTag	theTracksLabel_

edm::EDGetTokenT< reco::TrackCollection >	theTracksToken_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

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 49 of file DTChamberEfficiency.h.

Constructor & Destructor Documentation

◆ DTChamberEfficiency()

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

Definition at line 49 of file DTChamberEfficiency.cc.

                                                                  {
   // Get the debug parameter for verbose output
   debug = pSet.getUntrackedParameter<bool>("debug", false);
  
   LogVerbatim("DTDQM|DTMonitorModule|DTChamberEfficiency") << "DTChamberEfficiency: constructor called";
  
   // service parameters
   ParameterSet serviceParameters = pSet.getParameter<ParameterSet>("ServiceParameters");
   theService = new MuonServiceProxy(serviceParameters, consumesCollector());
  
   theTracksLabel_ = pSet.getParameter<InputTag>("TrackCollection");
   theTracksToken_ = consumes<reco::TrackCollection>(theTracksLabel_);
  
   theMaxChi2 = static_cast<unsigned int>(pSet.getParameter<double>("theMaxChi2"));
   theNSigma = pSet.getParameter<double>("theNSigma");
   theMinNrec = static_cast<int>(pSet.getParameter<double>("theMinNrec"));
  
   labelRPCRecHits = pSet.getParameter<InputTag>("theRPCRecHits");
  
   thedt4DSegments = pSet.getParameter<InputTag>("dt4DSegments");
   thecscSegments = pSet.getParameter<InputTag>("cscSegments");
  
   edm::ConsumesCollector iC = consumesCollector();
  
   theMeasurementExtractor = new MuonDetLayerMeasurements(
       thedt4DSegments, thecscSegments, labelRPCRecHits, InputTag(), InputTag(), iC, true, false, false, false);
  
   theNavigationType = pSet.getParameter<string>("NavigationType");
  
   theEstimator = new Chi2MeasurementEstimator(theMaxChi2, theNSigma);
 }

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.

◆ ~DTChamberEfficiency()

DTChamberEfficiency::~DTChamberEfficiency ( )

override

Definition at line 81 of file DTChamberEfficiency.cc.

                                           {
   LogTrace("DTDQM|DTMonitorModule|DTChamberEfficiency") << "DTChamberEfficiency: destructor called";
  
   // free memory
   delete theService;
   delete theMeasurementExtractor;
   delete theEstimator;
 }

References LogTrace.

Member Function Documentation

◆ analyze()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 119 of file DTChamberEfficiency.cc.

                                                                                   {
   LogTrace("DTDQM|DTMonitorModule|DTChamberEfficiency")
       << "--- [DTChamberEfficiency] Event analysed #Run: " << event.id().run() << " #Event: " << event.id().event()
       << endl;
  
   theService->update(eventSetup);
   theMeasurementExtractor->setEvent(event);
  
   //Read tracks from event
   Handle<reco::TrackCollection> tracks;
   event.getByToken(theTracksToken_, tracks);
  
   if (tracks.isValid()) {  // check the validity of the collection
  
     const edm::ESHandle<GlobalTrackingGeometry>& globalTrackingGeometry = theService->trackingGeometry();
     const MagneticField* magneticField = theService->magneticField().product();
  
     //loop over the muons
     for (reco::TrackCollection::const_iterator track = tracks->begin(); track != tracks->end(); ++track) {
       reco::TransientTrack trans_track(*track, magneticField, globalTrackingGeometry);
       const int recHitsize = (int)trans_track.recHitsSize();
       if (recHitsize < theMinNrec)
         continue;
  
       // printout the DT rechits used by the track
       if (debug) {
         LogTrace("DTDQM|DTMonitorModule|DTChamberEfficiency") << "--- New track" << endl;
         set<DTChamberId> chAlrUsed;
         for (trackingRecHit_iterator rHit = trans_track.recHitsBegin(); rHit != trans_track.recHitsEnd(); ++rHit) {
           DetId rHitid = (*rHit)->geographicalId();
           if (!(rHitid.det() == DetId::Muon && rHitid.subdetId() == MuonSubdetId::DT))
             continue;
           DTChamberId wId(rHitid.rawId());
           if (chAlrUsed.find(wId) != chAlrUsed.end())
             continue;
           chAlrUsed.insert(wId);
           LogTrace("DTDQM|DTMonitorModule|DTChamberEfficiency") << "   " << wId << endl;
         }
       }
  
       // Get the layer on which the seed relies
       DetId id = trans_track.track().innerDetId();
       const DetLayer* initialLayer = theService->detLayerGeometry()->idToLayer(id);
  
       TrajectoryStateOnSurface init_fs = trans_track.innermostMeasurementState();
       const FreeTrajectoryState* init_fs_free = init_fs.freeState();
  
       //get the list of compatible layers
       vector<const DetLayer*> layer_list =
           compatibleLayers(*theService->muonNavigationSchool(), initialLayer, *init_fs_free, alongMomentum);
       vector<const DetLayer*> layer_list_2 =
           compatibleLayers(*theService->muonNavigationSchool(), initialLayer, *init_fs_free, oppositeToMomentum);
  
       layer_list.insert(layer_list.end(), layer_list_2.begin(), layer_list_2.end());
  
       set<DTChamberId> alreadyCheckedCh;
  
       //loop over the list of compatible layers
       for (int i = 0; i < (int)layer_list.size(); i++) {
         //propagate the track to the i-th layer
         TrajectoryStateOnSurface tsos = propagator()->propagate(init_fs, layer_list.at(i)->surface());
         if (!tsos.isValid())
           continue;
  
         //determine the chambers kinematically compatible with the track on the i-th layer
         vector<DetWithState> dss = layer_list.at(i)->compatibleDets(tsos, *propagator(), *theEstimator);
  
         if (dss.empty())
           continue;
  
         // get the first det (it's the most compatible)
         const DetWithState detWithState = dss.front();
         const DetId idDetLay = detWithState.first->geographicalId();
  
         // check if this is a DT and the track has the needed quality
         if (!chamberSelection(idDetLay, trans_track))
           continue;
  
         DTChamberId DTid = (DTChamberId)idDetLay;
  
         // check if the chamber has already been counted
         if (alreadyCheckedCh.find(DTid) != alreadyCheckedCh.end())
           continue;
         alreadyCheckedCh.insert(DTid);
  
         // get the compatible measurements
         MeasurementContainer detMeasurements_initial = theMeasurementExtractor->measurements(
             layer_list.at(i), detWithState.first, detWithState.second, *theEstimator, event);
         LogTrace("DTDQM|DTMonitorModule|DTChamberEfficiency")
             << "     chamber: " << DTid << " has: " << detMeasurements_initial.size() << " comp. meas." << endl;
  
         //we want to be more picky about the quality of the segments:
         //exclude the segments with less than 12 hits
         MeasurementContainer detMeasurements = segQualityCut(detMeasurements_initial);
  
         // get the histos for this chamber
         vector<MonitorElement*> histos = histosPerW[DTid.wheel() + 2];
         // fill them
         if (!detMeasurements_initial.empty())
           histos[0]->Fill(DTid.sector(), DTid.station(), 1.);
         if (!detMeasurements.empty())
           histos[1]->Fill(DTid.sector(), DTid.station(), 1.);
         histos[2]->Fill(DTid.sector(), DTid.station(), 1.);
       }
     }
   } else {
     LogInfo("DTDQM|DTMonitorModule|DTChamberEfficiency")
         << "[DTChamberEfficiency] Collection: " << theTracksLabel_ << " is not valid!" << endl;
   }
   return;
 }

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

◆ bookHistograms()

void DTChamberEfficiency::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	iRun,
		edm::EventSetup const &
	)

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 92 of file DTChamberEfficiency.cc.

                                                                                                            {
   LogTrace("DTDQM|DTMonitorModule|DTChamberEfficiency") << "DTChamberEfficiency: booking histos";
  
   // Create the monitor elements
   ibooker.setCurrentFolder("DT/05-ChamberEff/Task");
  
   for (int wheel = -2; wheel <= 2; wheel++) {
     vector<MonitorElement*> histos;
  
     stringstream wheel_str;
     wheel_str << wheel;
  
     histos.push_back(ibooker.book2D(
         "hCountSectVsChamb_All_W" + wheel_str.str(), "Countings for wheel " + wheel_str.str(), 14, 1., 15., 4, 1., 5.));
  
     histos.push_back(ibooker.book2D(
         "hCountSectVsChamb_Qual_W" + wheel_str.str(), "Countings for wheel " + wheel_str.str(), 14, 1., 15., 4, 1., 5.));
  
     histos.push_back(ibooker.book2D(
         "hExtrapSectVsChamb_W" + wheel_str.str(), "Extrapolations for wheel " + wheel_str.str(), 14, 1., 15., 4, 1., 5.));
  
     histosPerW.push_back(histos);
   }
  
   return;
 }

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

◆ chamberSelection()

bool DTChamberEfficiency::chamberSelection	(	const DetId &	idDetLay,
		reco::TransientTrack &	trans_track
	)		const

private

Definition at line 231 of file DTChamberEfficiency.cc.

                                                                                                        {
   //check that we have a muon and that is a DT detector
   if (!(idDetLay.det() == DetId::Muon && idDetLay.subdetId() == MuonSubdetId::DT))
     return false;
  
   if (trans_track.recHitsSize() == 2)
     if (trans_track.recHit(0)->geographicalId() == idDetLay || trans_track.recHit(1)->geographicalId() == idDetLay)
       return false;
  
   return true;
 }

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

◆ compatibleLayers()

vector< const DetLayer * > DTChamberEfficiency::compatibleLayers	(	const NavigationSchool &	navigationSchool,
		const DetLayer *	initialLayer,
		const FreeTrajectoryState &	fts,
		PropagationDirection	propDir
	)

private

Definition at line 270 of file DTChamberEfficiency.cc.

                                                                                             {
   vector<const DetLayer*> detLayers;
  
   if (theNavigationType == "Standard") {
     // ask for compatible layers
     detLayers = navigationSchool.compatibleLayers(*initialLayer, fts, propDir);
     // I have to fit by hand the first layer until the seedTSOS is defined on the first rechit layer
     // In fact the first layer is not returned by initialLayer->compatibleLayers.
  
     detLayers.insert(detLayers.begin(), initialLayer);
  
   } else if (theNavigationType == "Direct") {
     DirectMuonNavigation navigation(ESHandle<MuonDetLayerGeometry>(&*theService->detLayerGeometry()));
     detLayers = navigation.compatibleLayers(fts, propDir);
   } else
     LogError("DTDQM|DTMonitorModule|DTChamberEfficiency") << "No Properly Navigation Selected!!" << endl;
  
   return detLayers;
 }

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

◆ dqmBeginRun()

void DTChamberEfficiency::dqmBeginRun	(	const edm::Run &	,
		const edm::EventSetup &
	)

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 90 of file DTChamberEfficiency.cc.

90 {}

◆ propagator()

ESHandle< Propagator > DTChamberEfficiency::propagator ( ) const

inlineprivate

Definition at line 293 of file DTChamberEfficiency.cc.

                                                                   {
   return theService->propagator("SteppingHelixPropagatorAny");
 }

◆ segQualityCut()

MeasurementContainer DTChamberEfficiency::segQualityCut ( const MeasurementContainer & seg_list ) const

private

Definition at line 243 of file DTChamberEfficiency.cc.

                                                                                                   {
   MeasurementContainer result;
  
   for (MeasurementContainer::const_iterator mescont_Itr = seg_list.begin(); mescont_Itr != seg_list.end();
        ++mescont_Itr) {
     //get the rechits of the segment
     TransientTrackingRecHit::ConstRecHitContainer recHit_list = mescont_Itr->recHit()->transientHits();
  
     //loop over the rechits and get the number of hits
     int nhit_seg(0);
     for (TransientTrackingRecHit::ConstRecHitContainer::const_iterator recList_Itr = recHit_list.begin();
          recList_Itr != recHit_list.end();
          ++recList_Itr) {
       nhit_seg += (int)(*recList_Itr)->transientHits().size();
     }
  
     DTChamberId tmpId = (DTChamberId)mescont_Itr->recHit()->hit()->geographicalId();
  
     if (tmpId.station() < 4 && nhit_seg >= 12)
       result.push_back(*mescont_Itr);
     if (tmpId.station() == 4 && nhit_seg >= 8)
       result.push_back(*mescont_Itr);
   }
  
   return result;
 }

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

Member Data Documentation

◆ debug

bool DTChamberEfficiency::debug

private

Definition at line 77 of file DTChamberEfficiency.h.

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

◆ histosPerW

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

private

Definition at line 96 of file DTChamberEfficiency.h.

◆ labelRPCRecHits

edm::InputTag DTChamberEfficiency::labelRPCRecHits

private

Definition at line 82 of file DTChamberEfficiency.h.

◆ thecscSegments

edm::InputTag DTChamberEfficiency::thecscSegments

private

Definition at line 84 of file DTChamberEfficiency.h.

◆ thedt4DSegments

edm::InputTag DTChamberEfficiency::thedt4DSegments

private

Definition at line 83 of file DTChamberEfficiency.h.

◆ theEstimator

Chi2MeasurementEstimator* DTChamberEfficiency::theEstimator

private

Definition at line 94 of file DTChamberEfficiency.h.

◆ theMaxChi2

double DTChamberEfficiency::theMaxChi2

private

Definition at line 86 of file DTChamberEfficiency.h.

◆ theMeasurementExtractor

MuonDetLayerMeasurements* DTChamberEfficiency::theMeasurementExtractor

private

Definition at line 93 of file DTChamberEfficiency.h.

◆ theMinNrec

int DTChamberEfficiency::theMinNrec

private

Definition at line 88 of file DTChamberEfficiency.h.

◆ theNavigationType

std::string DTChamberEfficiency::theNavigationType

private

Definition at line 90 of file DTChamberEfficiency.h.

◆ theNSigma

double DTChamberEfficiency::theNSigma

private

Definition at line 87 of file DTChamberEfficiency.h.

◆ theService

MuonServiceProxy* DTChamberEfficiency::theService

private

Definition at line 92 of file DTChamberEfficiency.h.

◆ theTracksLabel_

edm::InputTag DTChamberEfficiency::theTracksLabel_

private

Definition at line 79 of file DTChamberEfficiency.h.

◆ theTracksToken_

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

private

Definition at line 80 of file DTChamberEfficiency.h.

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ DTChamberEfficiency()

◆ ~DTChamberEfficiency()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ chamberSelection()

◆ compatibleLayers()

◆ dqmBeginRun()

◆ propagator()

◆ segQualityCut()

Member Data Documentation

◆ debug

◆ histosPerW

◆ labelRPCRecHits

◆ thecscSegments

◆ thedt4DSegments

◆ theEstimator

◆ theMaxChi2

◆ theMeasurementExtractor

◆ theMinNrec

◆ theNavigationType

◆ theNSigma

◆ theService

◆ theTracksLabel_

◆ theTracksToken_