#include <DTTriggerEfficiencyTask.h>

Inheritance diagram for DTTriggerEfficiencyTask:

Public Member Functions
	DTTriggerEfficiencyTask (const edm::ParameterSet &ps)
	Constructor. More...

	~DTTriggerEfficiencyTask () override
	Destructor. More...

Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override
	Analyze. More...

void	bookChamberHistos (DQMStore::IBooker &ibooker, const DTChamberId &dtCh, std::string histoTag, std::string folder="")
	Book chamber granularity histograms. More...

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

void	bookWheelHistos (DQMStore::IBooker &ibooker, int wheel, std::string histoTag, std::string folder="")
	Book wheel granularity histograms. More...

void	dqmBeginRun (const edm::Run &, const edm::EventSetup &) override
	BeginRun. More...

bool	hasRPCTriggers (const edm::Event &e)
	checks for RPC Triggers More...

std::string	topFolder (std::string source)
	return the top folder More...

Private Attributes
std::map< uint32_t, std::map< std::string, MonitorElement * > >	chamberHistos

bool	checkRPCtriggers

edm::EDGetTokenT< DTLocalTriggerCollection >	ddu_Token_

bool	detailedPlots

edm::EDGetTokenT< L1MuGMTReadoutCollection >	gmt_Token_

edm::InputTag	inputTagSEG

int	maxBXDDU

int	minBXDDU

edm::ESHandle< DTGeometry >	muonGeom

edm::EDGetTokenT< reco::MuonCollection >	muons_Token_

int	nevents
	To reset the MEs. More...

int	nMinHitsPhi

edm::ParameterSet	parameters

float	phiAccRange

bool	processDDU

std::vector< std::string >	processTags

bool	processTM

std::string	SegmArbitration

edm::EDGetTokenT< L1MuDTChambPhContainer >	tm_Token_

DTTrigGeomUtils *	trigGeomUtils

std::map< int, std::map< std::string, MonitorElement * > >	wheelHistos

Detailed Description

Definition at line 45 of file DTTriggerEfficiencyTask.h.

Constructor & Destructor Documentation

DTTriggerEfficiencyTask::DTTriggerEfficiencyTask ( const edm::ParameterSet & ps )

Constructor.

Definition at line 43 of file DTTriggerEfficiencyTask.cc.

References checkRPCtriggers, ddu_Token_, detailedPlots, gmt_Token_, inputTagSEG, LogTrace, maxBXDDU, minBXDDU, muons_Token_, nMinHitsPhi, phiAccRange, processDDU, processTags, processTM, SegmArbitration, AlCaHLTBitMon_QueryRunRegistry::string, and tm_Token_.

                                                                           : trigGeomUtils(nullptr) {
 
   LogTrace ("DTDQM|DTMonitorModule|DTTriggerEfficiencyTask")  << "[DTTriggerEfficiencyTask]: Constructor" << endl;
 
   parameters = ps;
 
   muons_Token_ = consumes<reco::MuonCollection>(
       parameters.getUntrackedParameter<edm::InputTag>("inputTagMuons"));
   tm_Token_   = consumes<L1MuDTChambPhContainer>(
       parameters.getUntrackedParameter<edm::InputTag>("inputTagTM"));
   ddu_Token_   = consumes<DTLocalTriggerCollection>(
       parameters.getUntrackedParameter<edm::InputTag>("inputTagDDU"));
   inputTagSEG  = parameters.getUntrackedParameter<edm::InputTag>("inputTagSEG");
   gmt_Token_   = consumes<L1MuGMTReadoutCollection>(
       parameters.getUntrackedParameter<edm::InputTag>("inputTagGMT"));
 
   SegmArbitration = parameters.getUntrackedParameter<std::string>("SegmArbitration");
 
   detailedPlots = parameters.getUntrackedParameter<bool>("detailedAnalysis");
   processTM = parameters.getUntrackedParameter<bool>("processTM");
   processDDU = parameters.getUntrackedParameter<bool>("processDDU");
   minBXDDU = parameters.getUntrackedParameter<int>("minBXDDU");
   maxBXDDU = parameters.getUntrackedParameter<int>("maxBXDDU");
 
   checkRPCtriggers = parameters.getUntrackedParameter<bool>("checkRPCtriggers");
   nMinHitsPhi = parameters.getUntrackedParameter<int>("nMinHitsPhi");
   phiAccRange = parameters.getUntrackedParameter<double>("phiAccRange");
 
   if (processTM) processTags.push_back("TM");
   if (processDDU) {processTags.push_back("DDU");
           ddu_Token_   = consumes<DTLocalTriggerCollection>(
       parameters.getUntrackedParameter<edm::InputTag>("inputTagDDU"));
     }
   if (!processTM && !processDDU) LogError ("DTDQM|DTMonitorModule|DTTriggerEfficiencyTask")  << "[DTTriggerEfficiencyTask]: Error, no trigger source (DDU or Twinmux) has been selected!!" <<endl;
 
 }

DTTriggerEfficiencyTask::~DTTriggerEfficiencyTask ( )

override

Destructor.

Definition at line 81 of file DTTriggerEfficiencyTask.cc.

References LogTrace, and nevents.

                                                   {
 
   LogTrace ("DTDQM|DTMonitorModule|DTTriggerEfficiencyTask")  << "[DTTriggerEfficiencyTask]: analyzed " << nevents << " events" << endl;
 
 }

Member Function Documentation

void DTTriggerEfficiencyTask::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

overrideprotected

Analyze.

Definition at line 120 of file DTTriggerEfficiencyTask.cc.

References reco::MuonSegmentMatch::BelongsToTrackByCleaning, reco::MuonSegmentMatch::BelongsToTrackByDR, reco::MuonSegmentMatch::BestInChamberByDR, relativeConstraints::chamber, chamberHistos, checkRPCtriggers, DTTrigGeomUtils::computeSCCoordinates(), ddu_Token_, detailedPlots, MuonSubdetId::DT, edm::Event::getByToken(), L1MuDTChambPhContainer::getContainer(), hasRPCTriggers(), triggerObjects_cff::id, RPCpg::mu, extraflags_cff::muons, muons_Token_, nevents, nMinHitsPhi, phiAccRange, processDDU, processTags, jets_cff::quality, DetId::rawId(), SegmArbitration, relativeConstraints::station, tm_Token_, trigGeomUtils, makeMuonMisalignmentScenario::wheel, wheelHistos, x, xdir, y, and ydir.

Referenced by topFolder().

                                                                               {
 
   nevents++;
 
   if (checkRPCtriggers){ //For pre-2016 Era compatibility
     if (!hasRPCTriggers(e)) { return; }
     }
   map<DTChamberId,const L1MuDTChambPhDigi*> phBestTM;
   map<DTChamberId,const DTLocalTrigger*>    phBestDDU;
   // Getting best TM Stuff
   edm::Handle<L1MuDTChambPhContainer> l1DTTPGPh;
   e.getByToken(tm_Token_, l1DTTPGPh);
   vector<L1MuDTChambPhDigi> const*  phTrigs = l1DTTPGPh->getContainer();
   //empty from dttfDigis, needs emulator working?
   vector<L1MuDTChambPhDigi>::const_iterator iph  = phTrigs->begin();
   vector<L1MuDTChambPhDigi>::const_iterator iphe = phTrigs->end();
   for(; iph !=iphe ; ++iph) {
 
     int phwheel = iph->whNum();
     int phsec   = iph->scNum() + 1; // DTTF numbering [0:11] -> DT numbering [1:12]
     int phst    = iph->stNum();
     int phcode  = iph->code();
 
     DTChamberId chId(phwheel,phst,phsec);
 
     if( phcode < 7 && (phBestTM.find(chId) == phBestTM.end() ||
           phcode>phBestTM[chId]->code()) ) phBestTM[chId] = &(*iph);
   }
 
   //Getting Best DDU Stuff
   if (processDDU){
   Handle<DTLocalTriggerCollection> trigsDDU;
   e.getByToken(ddu_Token_, trigsDDU);
   DTLocalTriggerCollection::DigiRangeIterator detUnitIt;
 
   for (detUnitIt=trigsDDU->begin();detUnitIt!=trigsDDU->end();++detUnitIt){
 
     const DTChamberId& id = (*detUnitIt).first;
     const DTLocalTriggerCollection::Range& range = (*detUnitIt).second;
 
     DTLocalTriggerCollection::const_iterator trigIt  = range.first;
     DTLocalTriggerCollection::const_iterator trigEnd = range.second;
     for (; trigIt!= trigEnd;++trigIt){
       int quality = trigIt->quality();
       if(quality>-1 && quality<7 &&
           (phBestDDU.find(id) == phBestDDU.end() ||
            quality>phBestDDU[id]->quality()) ) phBestDDU[id] = &(*trigIt);
     }
 
   }
   }//processDDU
 
   //Getting Best Segments
   vector<const DTRecSegment4D*> best4DSegments;
 
   Handle<reco::MuonCollection> muons;
   e.getByToken(muons_Token_, muons);
   reco::MuonCollection::const_iterator mu;
 
   for( mu = muons->begin(); mu != muons->end(); ++mu ) {
 
     // Make sure that is standalone muon
     if( !((*mu).isStandAloneMuon()) ) {continue;}
 
     // Get the chambers compatible with the muon
     const vector<reco::MuonChamberMatch> matchedChambers = (*mu).matches();
     vector<reco::MuonChamberMatch>::const_iterator chamber;
 
     for( chamber = matchedChambers.begin(); chamber != matchedChambers.end(); ++chamber ) {
 
       // look only in DTs
       if( chamber->detector() != MuonSubdetId::DT ) {continue;}
 
       // Get the matched segments in the chamber
       const vector<reco::MuonSegmentMatch> matchedSegments = (*chamber).segmentMatches;
       vector<reco::MuonSegmentMatch>::const_iterator segment;
 
       for( segment = matchedSegments.begin(); segment != matchedSegments.end(); ++segment ) {
 
         edm::Ref<DTRecSegment4DCollection> dtSegment = segment->dtSegmentRef;
 
         // Segment Arbitration
         if( SegmArbitration == "SegmentArbitration"
             && !((*segment).isMask(reco::MuonSegmentMatch::BestInChamberByDR)) ) {continue;}
 
         if( SegmArbitration == "SegmentAndTrackArbitration"
             && (!((*segment).isMask(reco::MuonSegmentMatch::BestInChamberByDR)) ||
               !((*segment).isMask(reco::MuonSegmentMatch::BelongsToTrackByDR))) ) {continue;}
 
         if( SegmArbitration == "SegmentAndTrackArbitrationCleaned"
             && (!((*segment).isMask(reco::MuonSegmentMatch::BestInChamberByDR))  ||
               !((*segment).isMask(reco::MuonSegmentMatch::BelongsToTrackByDR)) ||
               !((*segment).isMask(reco::MuonSegmentMatch::BelongsToTrackByCleaning))) ) {continue;}
 
 
         if( (*dtSegment).hasPhi() ) {
           best4DSegments.push_back(&(*dtSegment));
         }
 
       }// end loop on matched segments
     }// end loop on compatible chambers
   }// end loop on muons
 
   // Plot filling
   vector<const DTRecSegment4D*>::const_iterator btrack;
   for ( btrack = best4DSegments.begin(); btrack != best4DSegments.end(); ++btrack ){
     int wheel    = (*btrack)->chamberId().wheel();
     int station  = (*btrack)->chamberId().station();
     int scsector = 0;
     float x, xdir, y, ydir;
     trigGeomUtils->computeSCCoordinates((*btrack),scsector,x,xdir,y,ydir);
     int nHitsPhi = (*btrack)->phiSegment()->degreesOfFreedom()+2;
     DTChamberId dtChId(wheel,station,scsector);
     uint32_t indexCh = dtChId.rawId();
     map<string, MonitorElement*> &innerChME = chamberHistos[indexCh];
     map<string, MonitorElement*> &innerWhME = wheelHistos[wheel];
 
     if (fabs(xdir)<phiAccRange && nHitsPhi>=nMinHitsPhi){
       vector<string>::const_iterator tagIt  = processTags.begin();
       vector<string>::const_iterator tagEnd = processTags.end();
       for (; tagIt!=tagEnd; ++tagIt) {
     int qual = (*tagIt) == "TM" ?
           phBestTM.find(dtChId) != phBestTM.end() ? phBestTM[dtChId]->code() : -1 :
           phBestDDU.find(dtChId) != phBestDDU.end() ? phBestDDU[dtChId]->quality() : -1;
         innerWhME.find((*tagIt) + "_TrigEffDenum")->second->Fill(scsector,station);
 
         if ( qual>=0 && qual<7 ) {
           innerWhME.find((*tagIt) + "_TrigEffNum")->second->Fill(scsector,station);
           if ( qual>=4 ) {
             innerWhME.find((*tagIt) + "_TrigEffCorrNum")->second->Fill(scsector,station);
           }
         }
         if (detailedPlots) {
           innerChME.find((*tagIt) + "_TrackPosvsAngle")->second->Fill(xdir,x);
           if ( qual>=0 && qual<7 ) {
             innerChME.find((*tagIt) + "_TrackPosvsAngleAnyQual")->second->Fill(xdir,x);
             if ( qual>=4 ) {
               innerChME.find((*tagIt) + "_TrackPosvsAngleCorr")->second->Fill(xdir,x);
             }
           }
         }
       }
     }
   }
 
 }

void DTTriggerEfficiencyTask::bookChamberHistos	(	DQMStore::IBooker &	ibooker,
		const DTChamberId &	dtCh,
		std::string	histoTag,
		std::string	folder = `""`
	)

protected

Book chamber granularity histograms.

Definition at line 296 of file DTTriggerEfficiencyTask.cc.

References DQMStore::IBooker::book2D(), chamberHistos, printsummarytable::folder, fftjetimagerecorder_cfi::histoLabel, LogTrace, SiStripPI::max, min(), pileupCalc::nbins, DTTrigGeomUtils::phiRange(), DetId::rawId(), SimDataFormats::CaloAnalysis::sc, DTChamberId::sector(), DQMStore::IBooker::setCurrentFolder(), DTChamberId::station(), relativeConstraints::station, topFolder(), trigGeomUtils, DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

Referenced by bookHistograms().

                                                                                    {
 
   int wh = dtCh.wheel();
   int sc = dtCh.sector();
   int st = dtCh.station();
   stringstream wheel; wheel << wh;
   stringstream station; station << st;
   stringstream sector; sector << sc;
 
   string hwFolder      = topFolder(histoType);
   string bookingFolder = hwFolder + "Wheel" + wheel.str() + "/Sector" + sector.str() + "/Station" + station.str() + "/" + folder;
   string histoTag      = "_W" + wheel.str() + "_Sec" + sector.str() + "_St" + station.str();
 
   ibooker.setCurrentFolder(bookingFolder);
 
   LogTrace ("DTDQM|DTMonitorModule|DTTriggerEfficiencyTask")
     << "[DTTriggerEfficiencyTask]: booking histos in " << bookingFolder << endl;
 
   float min, max;
   int nbins;
   trigGeomUtils->phiRange(dtCh,min,max,nbins,20);
 
   string histoName = histoType + "_TrackPosvsAngle" +  histoTag;
   string histoLabel = "Position vs Angle (phi)";
 
   (chamberHistos[dtCh.rawId()])[histoType + "_TrackPosvsAngle"] =
     ibooker.book2D(histoName,histoLabel,12,-30.,30.,nbins,min,max);
 
   histoName = histoType + "_TrackPosvsAngleAnyQual" +  histoTag;
   histoLabel = "Position vs Angle (phi) for any qual triggers";
 
   (chamberHistos[dtCh.rawId()])[histoType + "_TrackPosvsAngleAnyQual"] =
     ibooker.book2D(histoName,histoLabel,12,-30.,30.,nbins,min,max);
 
   histoName = histoType + "_TrackPosvsAngleCorr" +  histoTag;
   histoLabel = "Position vs Angle (phi) for correlated triggers";
 
   (chamberHistos[dtCh.rawId()])[histoType + "_TrackPosvsAngleCorr"] =
     ibooker.book2D(histoName,histoLabel,12,-30.,30.,nbins,min,max);
 
 }

void DTTriggerEfficiencyTask::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	run,
		edm::EventSetup const &	context
	)

overrideprotected

Definition at line 95 of file DTTriggerEfficiencyTask.cc.

References bookChamberHistos(), bookWheelHistos(), detailedPlots, LogTrace, nevents, processTags, and trackingPlots::stat.

                                                                             {
 
   LogTrace ("DTDQM|DTMonitorModule|DTTriggerEfficiencyTask") << "[DTTriggerEfficiencyTask]: bookHistograms" << endl;
 
   nevents = 0;
   for (int wh=-2;wh<=2;++wh){
     vector<string>::const_iterator tagIt  = processTags.begin();
     vector<string>::const_iterator tagEnd = processTags.end();
     for (; tagIt!=tagEnd; ++tagIt) {
 
       bookWheelHistos(ibooker,wh,(*tagIt),"Task");
       if (detailedPlots) {
         for (int stat=1;stat<=4;++stat){
           for (int sect=1;sect<=12;++sect){
             bookChamberHistos(ibooker,DTChamberId(wh,stat,sect),(*tagIt),"Segment");
           }
         }
       }
     }
   }
 }

void DTTriggerEfficiencyTask::bookWheelHistos	(	DQMStore::IBooker &	ibooker,
		int	wheel,
		std::string	histoTag,
		std::string	folder = `""`
	)

protected

Book wheel granularity histograms.

Definition at line 339 of file DTTriggerEfficiencyTask.cc.

References DQMStore::IBooker::book2D(), LogTrace, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::IBooker::setCurrentFolder(), topFolder(), makeMuonMisalignmentScenario::wheel, and wheelHistos.

Referenced by bookHistograms().

                                                                {
 
   stringstream wh; wh << wheel;
   string basedir;
   if (hTag.find("Summary") != string::npos ) {
     basedir = topFolder(hTag);   //Book summary histo outside folder directory
   } else {
     basedir = topFolder(hTag) + folder + "/" ;
   }
 
   ibooker.setCurrentFolder(basedir);
 
   string hTagName = "_W" + wh.str();
 
   LogTrace("DTDQM|DTMonitorModule|DTTriggerEfficiencyTask")
     << "[DTTriggerEfficiencyTask]: booking histos in "<< basedir << endl;
 
   string hName = hTag + "_TrigEffDenum" + hTagName;
 
   MonitorElement* me = ibooker.book2D(hName.c_str(),hName.c_str(),12,1,13,4,1,5);
 
   me->setBinLabel(1,"MB1",2);
   me->setBinLabel(2,"MB2",2);
   me->setBinLabel(3,"MB3",2);
   me->setBinLabel(4,"MB4",2);
   me->setAxisTitle("Sector",1);
 
   wheelHistos[wheel][hTag + "_TrigEffDenum"] = me;
 
   hName = hTag + "_TrigEffNum" + hTagName;
   me = ibooker.book2D(hName.c_str(),hName.c_str(),12,1,13,4,1,5);
 
   me->setBinLabel(1,"MB1",2);
   me->setBinLabel(2,"MB2",2);
   me->setBinLabel(3,"MB3",2);
   me->setBinLabel(4,"MB4",2);
   me->setAxisTitle("Sector",1);
 
   wheelHistos[wheel][hTag + "_TrigEffNum"] = me;
 
   hName = hTag + "_TrigEffCorrNum" + hTagName;
   me = ibooker.book2D(hName.c_str(),hName.c_str(),12,1,13,4,1,5);
 
   me->setBinLabel(1,"MB1",2);
   me->setBinLabel(2,"MB2",2);
   me->setBinLabel(3,"MB3",2);
   me->setBinLabel(4,"MB4",2);
   me->setAxisTitle("Sector",1);
 
   wheelHistos[wheel][hTag + "_TrigEffCorrNum"] = me;
 
   return;
 }

void DTTriggerEfficiencyTask::dqmBeginRun	(	const edm::Run &	run,
		const edm::EventSetup &	context
	)

overrideprotected

BeginRun.

Definition at line 87 of file DTTriggerEfficiencyTask.cc.

References edm::EventSetup::get(), muonGeom, and trigGeomUtils.

                                                                                          {
 
   // Get the geometry
   context.get<MuonGeometryRecord>().get(muonGeom);
   trigGeomUtils = new DTTrigGeomUtils(muonGeom);
 
 }

bool DTTriggerEfficiencyTask::hasRPCTriggers ( const edm::Event & e )

protected

checks for RPC Triggers

Definition at line 267 of file DTTriggerEfficiencyTask.cc.

References edm::Event::getByToken(), L1MuGMTReadoutCollection::getRecords(), gmt_Token_, and jets_cff::quality.

Referenced by analyze().

                                                               {
 
   edm::Handle<L1MuGMTReadoutCollection> gmtrc;
   e.getByToken(gmt_Token_, gmtrc);
 
   std::vector<L1MuGMTReadoutRecord> gmt_records = gmtrc->getRecords();
   std::vector<L1MuGMTReadoutRecord>::const_iterator igmtrr = gmt_records.begin();
   std::vector<L1MuGMTReadoutRecord>::const_iterator egmtrr = gmt_records.end();
   for(; igmtrr!=egmtrr; igmtrr++) {
 
     std::vector<L1MuGMTExtendedCand> candsGMT = igmtrr->getGMTCands();
     std::vector<L1MuGMTExtendedCand>::const_iterator candGMTIt   = candsGMT.begin();
     std::vector<L1MuGMTExtendedCand>::const_iterator candGMTEnd  = candsGMT.end();
 
     for(; candGMTIt!=candGMTEnd; ++candGMTIt){
       if(!candGMTIt->empty()) {
         int quality = candGMTIt->quality();
         if(candGMTIt->bx()==0 &&
             (quality == 5 || quality == 7)){
           return true;
         }
       }
     }
   }
 
   return false;
 
 }