#include <DTTriggerEfficiencyTask.h>

Inheritance diagram for DTTriggerEfficiencyTask:

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

virtual	~DTTriggerEfficiencyTask ()
	Destructor. More...

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

void	beginJob ()

void	beginLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &context)
	To reset the MEs. More...

void	beginRun (const edm::Run &run, const edm::EventSetup &context)
	BeginRun. More...

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

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

void	endJob (void)
	EndJob. More...

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

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

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 Attributes
std::map< uint32_t, std::map < std::string, MonitorElement * > >	chamberHistos

DQMStore *	dbe

edm::EDGetTokenT < L1MuDTChambPhContainer >	dcc_Token_

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

int	nMinHitsPhi

edm::ParameterSet	parameters

float	phiAccRange

bool	processDCC

bool	processDDU

std::vector< std::string >	processTags

std::string	SegmArbitration

DTTrigGeomUtils *	trigGeomUtils

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

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Detailed Description

Definition at line 44 of file DTTriggerEfficiencyTask.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 43 of file DTTriggerEfficiencyTask.cc.

References dbe, dcc_Token_, ddu_Token_, detailedPlots, edm::ParameterSet::getUntrackedParameter(), gmt_Token_, inputTagSEG, LogTrace, maxBXDDU, minBXDDU, muons_Token_, nMinHitsPhi, cppFunctionSkipper::operator, parameters, phiAccRange, processDCC, processDDU, processTags, SegmArbitration, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                           : trigGeomUtils(0) {
 
   LogTrace ("DTDQM|DTMonitorModule|DTTriggerEfficiencyTask")  << "[DTTriggerEfficiencyTask]: Constructor" << endl;
 
   parameters = ps;
   dbe = edm::Service<DQMStore>().operator->();
 
   muons_Token_ = consumes<reco::MuonCollection>(
       parameters.getUntrackedParameter<edm::InputTag>("inputTagMuons"));
   dcc_Token_   = consumes<L1MuDTChambPhContainer>(
       parameters.getUntrackedParameter<edm::InputTag>("inputTagDCC"));
   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");
   processDCC = parameters.getUntrackedParameter<bool>("processDCC");
   processDDU = parameters.getUntrackedParameter<bool>("processDDU");
   minBXDDU = parameters.getUntrackedParameter<int>("minBXDDU");
   maxBXDDU = parameters.getUntrackedParameter<int>("maxBXDDU");
 
   nMinHitsPhi = parameters.getUntrackedParameter<int>("nMinHitsPhi");
   phiAccRange = parameters.getUntrackedParameter<double>("phiAccRange");
 
   if (processDCC) processTags.push_back("DCC");
   if (processDDU) processTags.push_back("DDU");
 
 }

DTTriggerEfficiencyTask::~DTTriggerEfficiencyTask ( )

virtual

Destructor.

Definition at line 77 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
	)

protectedvirtual

Analyze.

Implements edm::EDAnalyzer.

Definition at line 130 of file DTTriggerEfficiencyTask.cc.

References reco::MuonSegmentMatch::BelongsToTrackByCleaning, reco::MuonSegmentMatch::BelongsToTrackByDR, reco::MuonSegmentMatch::BestInChamberByDR, chamberHistos, DTTrigGeomUtils::computeSCCoordinates(), dcc_Token_, ddu_Token_, detailedPlots, MuonSubdetId::DT, edm::Event::getByToken(), hasRPCTriggers(), RPCpg::mu, patZpeak::muons, muons_Token_, nevents, nMinHitsPhi, phiAccRange, processTags, DetId::rawId(), SegmArbitration, relativeConstraints::station, trigGeomUtils, wheelHistos, x, xdir, detailsBasic3DVector::y, and ydir.

                                                                               {
 
   nevents++;
 
   if (!hasRPCTriggers(e)) { return; }
 
   map<DTChamberId,const L1MuDTChambPhDigi*> phBestDCC;
   map<DTChamberId,const DTLocalTrigger*>    phBestDDU;
 
   // Getting best DCC Stuff
   edm::Handle<L1MuDTChambPhContainer> l1DTTPGPh;
   e.getByToken(dcc_Token_, l1DTTPGPh);
   vector<L1MuDTChambPhDigi> const*  phTrigs = l1DTTPGPh->getContainer();
 
   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 && (phBestDCC.find(chId) == phBestDCC.end() ||
           phcode>phBestDCC[chId]->code()) ) phBestDCC[chId] = &(*iph);
   }
 
   //Getting Best DDU Stuff
   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);
     }
 
   }
 
   //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) == "DCC" ?
           phBestDCC.find(dtChId) != phBestDCC.end() ? phBestDCC[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::beginJob ( void )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 84 of file DTTriggerEfficiencyTask.cc.

References LogTrace.

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

void DTTriggerEfficiencyTask::beginLuminosityBlock	(	const edm::LuminosityBlock &	lumiSeg,
		const edm::EventSetup &	context
	)

protectedvirtual

To reset the MEs.

Reimplemented from edm::EDAnalyzer.

Definition at line 116 of file DTTriggerEfficiencyTask.cc.

References LogTrace.

                                                                                                             {
 
   LogTrace ("DTDQM|DTMonitorModule|DTTriggerEfficiencyTask") <<"[DTTriggerEfficiencyTask]: Begin of LS transition"<<endl;
 
 }

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

protectedvirtual

BeginRun.

Reimplemented from edm::EDAnalyzer.

Definition at line 90 of file DTTriggerEfficiencyTask.cc.

References bookChamberHistos(), bookWheelHistos(), detailedPlots, DTChamberId, edm::EventSetup::get(), LogTrace, muonGeom, nevents, processTags, and trigGeomUtils.

                                                                                      {
 
   LogTrace ("DTDQM|DTMonitorModule|DTTriggerEfficiencyTask") << "[DTTriggerEfficiencyTask]: BeginRun" << 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(wh,(*tagIt),"Task");
       if (detailedPlots) {
         for (int stat=1;stat<=4;++stat){
           for (int sect=1;sect<=12;++sect){
             bookChamberHistos(DTChamberId(wh,stat,sect),(*tagIt),"Segment");
           }
         }
       }
     }
   }
 
   context.get<MuonGeometryRecord>().get(muonGeom);
   trigGeomUtils = new DTTrigGeomUtils(muonGeom);
 
 }

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

protected

Book chamber granularity histograms.

Definition at line 308 of file DTTriggerEfficiencyTask.cc.

References DQMStore::book2D(), chamberHistos, dbe, LogTrace, max(), bookConverter::min, pileupCalc::nbins, DTTrigGeomUtils::phiRange(), DetId::rawId(), DTChamberId::sector(), DQMStore::setCurrentFolder(), DTChamberId::station(), relativeConstraints::station, topFolder(), trigGeomUtils, and DTChamberId::wheel().

Referenced by beginRun().

                                                                                                         {
 
   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();
 
   dbe->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"] =
     dbe->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"] =
     dbe->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"] =
     dbe->book2D(histoName,histoLabel,12,-30.,30.,nbins,min,max);
 
 }

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

protected

Book wheel granularity histograms.

Definition at line 347 of file DTTriggerEfficiencyTask.cc.

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

Referenced by beginRun().

                                                                                  {
 
   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 + "/" ;
   }
 
   dbe->setCurrentFolder(basedir);
 
   string hTagName = "_W" + wh.str();
 
   LogTrace("DTDQM|DTMonitorModule|DTTriggerEfficiencyTask")
     << "[DTTriggerEfficiencyTask]: booking histos in "<< basedir << endl;
 
   string hName = hTag + "_TrigEffDenum" + hTagName;
   MonitorElement* me = dbe->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 = dbe->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 = dbe->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::endJob ( void )

protectedvirtual

EndJob.

Reimplemented from edm::EDAnalyzer.

Definition at line 123 of file DTTriggerEfficiencyTask.cc.

References LogTrace, and nevents.

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

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

protected

checks for RPC Triggers

Definition at line 279 of file DTTriggerEfficiencyTask.cc.

References edm::Event::getByToken(), and gmt_Token_.

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