#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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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 >
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::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_FULL_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, edmPickEvents::event, HcalObjRepresent::Fill(), TrajectoryStateOnSurface::freeState(), combine::histos, mps_fire::i, reco::Track::innerDetId(), reco::TransientTrack::innermostMeasurementState(), createfilelist::int, TrajectoryStateOnSurface::isValid(), LogTrace, HLT_FULL_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_FULL_cff::track, 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_