#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 &ev, edm::EventSetup const &es) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

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

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

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

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

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

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

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () 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 () 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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

EDConsumerBase &	operator= (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)

virtual	~EDConsumerBase () noexcept(false)

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

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

MeasurementContainer	segQualityCut (const MeasurementContainer &seg_list) const

Private Attributes
bool	debug

edm::ESHandle< DTGeometry >	dtGeom

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

edm::InputTag	labelRPCRecHits

edm::ESHandle< MagneticField >	magfield

edm::InputTag	thecscSegments

edm::InputTag	thedt4DSegments

Chi2MeasurementEstimator *	theEstimator

double	theMaxChi2

MuonDetLayerMeasurements *	theMeasurementExtractor

int	theMinNrec

std::string	theNavigationType

double	theNSigma

MuonServiceProxy *	theService

edm::ESHandle< GlobalTrackingGeometry >	theTrackingGeometry

edm::InputTag	theTracksLabel_

edm::EDGetTokenT< reco::TrackCollection >	theTracksToken_

Additional Inherited Members
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< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

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

Constructor & Destructor Documentation

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

Definition at line 55 of file DTChamberEfficiency.cc.

References Chi2MeasurementEstimator_cfi::Chi2MeasurementEstimator, debug, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), MuonServiceProxy_cff::MuonServiceProxy, and createJobs::theNSigma.

 {
   // 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);
 
   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);
 }

DTChamberEfficiency::~DTChamberEfficiency ( )

override

Definition at line 90 of file DTChamberEfficiency.cc.

References LogTrace.

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

Member Function Documentation

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

override

Definition at line 143 of file DTChamberEfficiency.cc.

References alongMomentum, debug, DetId::det(), MuonSubdetId::DT, event(), HcalObjRepresent::Fill(), TrajectoryStateOnSurface::freeState(), plotFactory::histos, mps_fire::i, createfilelist::int, TrajectoryStateOnSurface::isValid(), edm::HandleBase::isValid(), LogTrace, DetId::Muon, oppositeToMomentum, PhotonConversionTrajectorySeedProducerFromQuadruplets_cfi::propagator, DetId::rawId(), reco::TransientTrack::recHitsSize(), DTChamberId::sector(), DTChamberId::station(), DetId::subdetId(), TrajectoryStateOnSurface::surface(), HiIsolationCommonParameters_cff::track, l1t::tracks, and DTChamberId::wheel().

 {
 
   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
 
     //loop over the muons
     for(reco::TrackCollection::const_iterator track = tracks->begin(); track!=tracks->end(); ++track) {
 
 
       reco::TransientTrack trans_track(*track,magfield.product(),theTrackingGeometry);
       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;
 }

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 113 of file DTChamberEfficiency.cc.

References DQMStore::IBooker::book2D(), plotFactory::histos, LogTrace, DQMStore::IBooker::setCurrentFolder(), and makeMuonMisalignmentScenario::wheel.

                                                                                                                       {
 
   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;
 }

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

private

Definition at line 253 of file DTChamberEfficiency.cc.

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

 {
   //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;
 }

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

private

Definition at line 293 of file DTChamberEfficiency.cc.

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

 {
 
 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;
 }

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

override

Definition at line 101 of file DTChamberEfficiency.cc.

References edm::EventSetup::get().

 {
   // Get the DT Geometry
   setup.get<MuonGeometryRecord>().get(dtGeom);
 
   setup.get<IdealMagneticFieldRecord>().get(magfield);
   setup.get<GlobalTrackingGeometryRecord>().get(theTrackingGeometry);
 
   return;
 }

ESHandle< Propagator > DTChamberEfficiency::propagator ( ) const

inlineprivate

Definition at line 319 of file DTChamberEfficiency.cc.

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

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

private

Definition at line 265 of file DTChamberEfficiency.cc.

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

 {
 
   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;
 }

Member Data Documentation

bool DTChamberEfficiency::debug

private

Definition at line 89 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(), and pkg.AbstractPkg::write().

edm::ESHandle<DTGeometry> DTChamberEfficiency::dtGeom

private

Definition at line 104 of file DTChamberEfficiency.h.

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

private

Definition at line 113 of file DTChamberEfficiency.h.

edm::InputTag DTChamberEfficiency::labelRPCRecHits

private

Definition at line 94 of file DTChamberEfficiency.h.

edm::ESHandle<MagneticField> DTChamberEfficiency::magfield

private

Definition at line 110 of file DTChamberEfficiency.h.

edm::InputTag DTChamberEfficiency::thecscSegments

private

Definition at line 96 of file DTChamberEfficiency.h.

edm::InputTag DTChamberEfficiency::thedt4DSegments

private

Definition at line 95 of file DTChamberEfficiency.h.

Chi2MeasurementEstimator* DTChamberEfficiency::theEstimator

private

Definition at line 108 of file DTChamberEfficiency.h.

double DTChamberEfficiency::theMaxChi2

private

Definition at line 98 of file DTChamberEfficiency.h.

MuonDetLayerMeasurements* DTChamberEfficiency::theMeasurementExtractor

private

Definition at line 107 of file DTChamberEfficiency.h.

int DTChamberEfficiency::theMinNrec

private

Definition at line 100 of file DTChamberEfficiency.h.

std::string DTChamberEfficiency::theNavigationType

private

Definition at line 102 of file DTChamberEfficiency.h.

double DTChamberEfficiency::theNSigma

private

Definition at line 99 of file DTChamberEfficiency.h.

MuonServiceProxy* DTChamberEfficiency::theService

private

Definition at line 106 of file DTChamberEfficiency.h.

edm::ESHandle<GlobalTrackingGeometry> DTChamberEfficiency::theTrackingGeometry

private

Definition at line 111 of file DTChamberEfficiency.h.

edm::InputTag DTChamberEfficiency::theTracksLabel_

private

Definition at line 91 of file DTChamberEfficiency.h.

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

private

Definition at line 92 of file DTChamberEfficiency.h.

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation