#include <DTLocalTriggerBaseTask.h>

Inheritance diagram for DTLocalTriggerBaseTask:

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

	~DTLocalTriggerBaseTask () 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	beginLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &context) override
	To reset the MEs. More...

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

void	endLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &context) override
	Perform trend plot operations. More...

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

void	bookHistos (DQMStore::IBooker &, const DTChamberId &chamb)
	Book the histograms. More...

void	bookHistos (DQMStore::IBooker &, int wh)
	Book the histograms. More...

void	runDDUAnalysis (edm::Handle< DTLocalTriggerCollection > &trigsDDU)
	Run analysis on ROS data. More...

void	runDDUvsTMAnalysis ()
	Run analysis on ROS data. More...

void	runTMAnalysis (std::vector< L1MuDTChambPhDigi > const phInTrigs, std::vector< L1MuDTChambPhDigi > const phOutTrigs, std::vector< L1MuDTChambThDigi > const *thTrigs)
	Run analysis on TM data. More...

void	setQLabels (MonitorElement *me, short int iaxis)
	Set Quality labels. More...

void	setQLabelsTheta (MonitorElement *me, short int iaxis)

std::string &	topFolder (std::string const &type)
	Get the Top folder (different between Physics and TP and TM/DDU) More...

Private Attributes
std::string	baseFolder [2]

int	bestAccRange

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

bool	detailedAnalysis

MonitorElement *	nEventMonitor

int	nEvents

int	nEventsInLS

int	nLumis

bool	processDDU

bool	processTM

int	targetBXDDU

int	targetBXTM

std::map< uint32_t, DTTPGCompareUnit >	theCompMapIn

std::map< uint32_t, DTTPGCompareUnit >	theCompMapOut

edm::ParameterSet	theParams

DTTrigGeomUtils *	theTrigGeomUtils

std::vector< std::string >	theTypes

edm::EDGetTokenT< L1MuDTChambPhContainer >	tm_phiIn_Token_

edm::EDGetTokenT< L1MuDTChambPhContainer >	tm_phiOut_Token_

edm::EDGetTokenT< L1MuDTChambThContainer >	tm_theta_Token_

bool	tpMode

std::map< uint32_t, DTTimeEvolutionHisto * >	trendHistos

edm::EDGetTokenT< DTLocalTriggerCollection >	trig_Token_

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

Friends
class	DTMonitorModule

Detailed Description

Definition at line 45 of file DTLocalTriggerBaseTask.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 63 of file DTLocalTriggerBaseTask.cc.

References bestAccRange, detailedAnalysis, edm::ParameterSet::getUntrackedParameter(), LogTrace, processDDU, processTM, targetBXDDU, targetBXTM, theParams, theTypes, tm_phiIn_Token_, tm_phiOut_Token_, tm_theta_Token_, topFolder(), tpMode, and trig_Token_.

                                                                         :
   nEvents(0), nLumis(0), theTrigGeomUtils(nullptr) {
 
   LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerBaseTask")
     << "[DTLocalTriggerBaseTask]: Constructor"<<endl;
 
   tpMode           = ps.getUntrackedParameter<bool>("testPulseMode");
   detailedAnalysis = ps.getUntrackedParameter<bool>("detailedAnalysis");
 
   targetBXTM  = ps.getUntrackedParameter<int>("targetBXTM");
   targetBXDDU  = ps.getUntrackedParameter<int>("targetBXDDU");
   bestAccRange = ps.getUntrackedParameter<int>("bestTrigAccRange");
 
   processTM   = ps.getUntrackedParameter<bool>("processTM");
   processDDU   = ps.getUntrackedParameter<bool>("processDDU");
 
   tm_phiIn_Token_ = consumes<L1MuDTChambPhContainer>(
       ps.getUntrackedParameter<InputTag>("inputTagTMphIn"));
   tm_phiOut_Token_ = consumes<L1MuDTChambPhContainer>(
       ps.getUntrackedParameter<InputTag>("inputTagTMphOut"));
   tm_theta_Token_ = consumes<L1MuDTChambThContainer>(
       ps.getUntrackedParameter<InputTag>("inputTagTMth"));
   trig_Token_ = consumes<DTLocalTriggerCollection>(
       ps.getUntrackedParameter<InputTag>("inputTagDDU"));
 
   if (processTM) theTypes.push_back("TM");
   if (processDDU) theTypes.push_back("DDU");
     
 
   if (tpMode) {
     topFolder("TM") = "DT/11-LocalTriggerTP-TM/";
     topFolder("DDU") = "DT/12-LocalTriggerTP-DDU/";
   } else {
     topFolder("TM") = "DT/03-LocalTrigger-TM/";
     topFolder("DDU") = "DT/04-LocalTrigger-DDU/";
   }
 
   theParams = ps;
 
 }

DTLocalTriggerBaseTask::~DTLocalTriggerBaseTask ( )

override

Destructor.

Definition at line 105 of file DTLocalTriggerBaseTask.cc.

References LogTrace, nEvents, and theTrigGeomUtils.

                                                 {
 
   LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerBaseTask")
     << "[DTLocalTriggerBaseTask]: analyzed " << nEvents << " events" << endl;
   if (theTrigGeomUtils) { delete theTrigGeomUtils; }
 
 }

Member Function Documentation

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

overrideprotected

Analyze.

Definition at line 177 of file DTLocalTriggerBaseTask.cc.

References MonitorElement::Fill(), edm::Event::getByToken(), L1MuDTChambPhContainer::getContainer(), L1MuDTChambThContainer::getContainer(), edm::ParameterSet::getUntrackedParameter(), edm::HandleBase::isValid(), nEventMonitor, nEvents, nEventsInLS, processDDU, processTM, runDDUAnalysis(), runDDUvsTMAnalysis(), runTMAnalysis(), theCompMapIn, theCompMapOut, theParams, tm_phiIn_Token_, tm_phiOut_Token_, tm_theta_Token_, and trig_Token_.

                                                                              {
 
   nEvents++;
   nEventsInLS++;
   nEventMonitor->Fill(nEvents);
 
   theCompMapIn.clear();
   theCompMapOut.clear();
 
   Handle<L1MuDTChambPhContainer> phiInTrigsTM;
   Handle<L1MuDTChambPhContainer> phiOutTrigsTM;
   Handle<L1MuDTChambThContainer> thetaTrigsTM;
   Handle<DTLocalTriggerCollection> trigsDDU;
 
   if (processTM) {
     InputTag inputTagTM = theParams.getUntrackedParameter<InputTag>("inputTagTM");
 
     e.getByToken(tm_phiIn_Token_, phiInTrigsTM);
     e.getByToken(tm_phiOut_Token_, phiOutTrigsTM);
     e.getByToken(tm_theta_Token_, thetaTrigsTM);
 
     if (phiInTrigsTM.isValid() && phiOutTrigsTM.isValid() && thetaTrigsTM.isValid()) {
       runTMAnalysis(phiInTrigsTM->getContainer(),phiOutTrigsTM->getContainer(),thetaTrigsTM->getContainer());
     } else {
       LogVerbatim("DTDQM|DTMonitorModule|DTLocalTriggerBaseTask")
     << "[DTLocalTriggerBaseTask]: one or more TM handles for Input Tag "
     << inputTagTM <<" not found!" << endl;
       return;
     }
   }
 
   if (processDDU) {
     InputTag inputTagDDU = theParams.getUntrackedParameter<InputTag>("inputTagDDU");
     e.getByToken(trig_Token_, trigsDDU);
 
     if (trigsDDU.isValid()) {
       runDDUAnalysis(trigsDDU);
     } else {
       LogVerbatim("DTDQM|DTMonitorModule|DTLocalTriggerBaseTask")
     << "[DTLocalTriggerBaseTask]: one or more DDU handles for Input Tag "
     << inputTagDDU <<" not found!" << endl;
       return;
     }
   }
 
   if (processTM && processDDU)
     runDDUvsTMAnalysis();
 
 }

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

overrideprotected

To reset the MEs.

Definition at line 128 of file DTLocalTriggerBaseTask.cc.

References chamberHistos, edm::ParameterSet::getUntrackedParameter(), LogTrace, nEventsInLS, nLumis, and theParams.

                                                                                                            {
 
   nEventsInLS=0;
   nLumis++;
   int resetCycle = theParams.getUntrackedParameter<int>("ResetCycle");
 
   LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerBaseTask")
     << "[DTLocalTriggerBaseTask]: Begin of LS transition" << endl;
 
   if( nLumis%resetCycle == 0 ) {
     map<uint32_t,map<string,MonitorElement*> >::const_iterator chambIt  = chamberHistos.begin();
     map<uint32_t,map<string,MonitorElement*> >::const_iterator chambEnd = chamberHistos.end();
     for(;chambIt!=chambEnd;++chambIt) {
       map<string,MonitorElement*>::const_iterator histoIt  = chambIt->second.begin();
       map<string,MonitorElement*>::const_iterator histoEnd = chambIt->second.end();
       for(;histoIt!=histoEnd;++histoIt) {
     histoIt->second->Reset();
       }
     }
   }
 
 }

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

overrideprivate

Definition at line 113 of file DTLocalTriggerBaseTask.cc.

References DQMStore::IBooker::bookFloat(), bookHistos(), nEventMonitor, processDDU, DQMStore::IBooker::setCurrentFolder(), and trackingPlots::stat.

Referenced by topFolder().

                                                                                                                          {
 
   ibooker.setCurrentFolder("DT/EventInfo/Counters");
   nEventMonitor = ibooker.bookFloat("nProcessedEventsTrigger");
   for (int wh=-2;wh<3;++wh){
     for (int stat=1;stat<5;++stat){
       for (int sect=1;sect<13;++sect){
     bookHistos(ibooker, DTChamberId(wh,stat,sect));
       }
     }
     if (processDDU) bookHistos(ibooker, wh);
   }
 }

void DTLocalTriggerBaseTask::bookHistos	(	DQMStore::IBooker &	ibooker,
		const DTChamberId &	chamb
	)

private

Book the histograms.

Definition at line 228 of file DTLocalTriggerBaseTask.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), chamberHistos, detailedAnalysis, edm::ParameterSet::getUntrackedParameter(), LogTrace, simTrackMatching_cfi::maxBX, simTrackMatching_cfi::minBX, DTTrigGeomUtils::phiRange(), processDDU, processTM, DetId::rawId(), DTChamberId::sector(), DQMStore::IBooker::setCurrentFolder(), setQLabels(), setQLabelsTheta(), DTChamberId::station(), relativeConstraints::station, theParams, theTrigGeomUtils, theTypes, topFolder(), tpMode, trendHistos, DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

Referenced by bookHistograms(), and topFolder().

                                                                                           {
 
   uint32_t rawId = dtCh.rawId();
 
   stringstream wheel; wheel << dtCh.wheel();
   stringstream station; station << dtCh.station();
   stringstream sector; sector << dtCh.sector();
 
   map<string,int> minBX;
   map<string,int> maxBX;
 
   minBX["TM"] = theParams.getUntrackedParameter<int>("minBXTM");
   maxBX["TM"] = theParams.getUntrackedParameter<int>("maxBXTM");
   minBX["DDU"] = theParams.getUntrackedParameter<int>("minBXDDU");
   maxBX["DDU"] = theParams.getUntrackedParameter<int>("maxBXDDU");
 
   int nTimeBins  = theParams.getUntrackedParameter<int>("nTimeBins");
   int nLSTimeBin = theParams.getUntrackedParameter<int>("nLSTimeBin");
 
   string chTag = "_W" + wheel.str() + "_Sec" + sector.str() + "_St" + station.str();
   string labelInOut = "";
 
   vector<string>::const_iterator typeIt  = theTypes.begin();
   vector<string>::const_iterator typeEnd = theTypes.end();
 
   for (; typeIt!=typeEnd; ++typeIt) {
 
     LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerBaseTask")
       << "[DTLocalTriggerBaseTask]: booking histos for " << topFolder((*typeIt)) << "Wheel"
       << wheel.str() << "/Sector" << sector.str() << "/Station"<< station.str() << endl;
 
     for (int InOut=0; InOut<2;InOut++){
     // Book Phi View Related Plots
     
     if(InOut==0)   
     {ibooker.setCurrentFolder(topFolder(*typeIt) + "Wheel" + wheel.str() + "/Sector"
               + sector.str() + "/Station" + station.str() + "/LocalTriggerPhiIn");
     labelInOut = "_In";}
     else if (InOut==1){
     ibooker.setCurrentFolder(topFolder(*typeIt) + "Wheel" + wheel.str() + "/Sector"
                           + sector.str() + "/Station" + station.str() + "/LocalTriggerPhiOut");
     labelInOut = "_Out";}
 
 
     string histoTag = (*typeIt) + "_BXvsQual" + labelInOut;
     chamberHistos[rawId][histoTag] = ibooker.book2D(histoTag+chTag,"BX vs trigger quality",
        7,-0.5,6.5,(int)(maxBX[(*typeIt)]-minBX[*typeIt]+1),minBX[*typeIt]-.5,maxBX[*typeIt]+.5);
     setQLabels((chamberHistos[rawId])[histoTag],1);
 
     if (!tpMode) {
       histoTag = (*typeIt) + "_BestQual" + labelInOut;
       chamberHistos[rawId][histoTag] = ibooker.book1D(histoTag+chTag,
              "Trigger quality of best primitives",7,-0.5,6.5);
       setQLabels(chamberHistos[rawId][histoTag],1);
       
       histoTag = (*typeIt) + "_Flag1stvsQual" + labelInOut;
       chamberHistos[dtCh.rawId()][histoTag] = ibooker.book2D(histoTag+chTag,
               "1st/2nd trig flag vs quality",7,-0.5,6.5,2,-0.5,1.5);
       setQLabels(chamberHistos[rawId][histoTag],1);
 
       histoTag = (*typeIt) + "_FlagUpDownvsQual" + labelInOut;
       chamberHistos[dtCh.rawId()][histoTag] = ibooker.book2D(histoTag+chTag,
                   "Up/Down trig flag vs quality",7,-0.5,6.5,2,-0.5,1.5);
       setQLabels(chamberHistos[rawId][histoTag],1);
 
     }
 
     if (*typeIt=="TM") {
       float minPh, maxPh; int nBinsPh;
       theTrigGeomUtils->phiRange(dtCh,minPh,maxPh,nBinsPh);
 
       histoTag = (*typeIt) + "_QualvsPhirad" + labelInOut;
       chamberHistos[rawId][histoTag] = ibooker.book2D(histoTag+chTag,
            "Trigger quality vs local position",nBinsPh,minPh,maxPh,7,-0.5,6.5);
       setQLabels(chamberHistos[rawId][histoTag],2);
 
       if (detailedAnalysis && !tpMode) {
     histoTag = (*typeIt) + "_QualvsPhibend" + labelInOut;
     chamberHistos[rawId][histoTag] = ibooker.book2D(histoTag+chTag,
           "Trigger quality vs local direction",200,-40.,40.,7,-0.5,6.5);
     setQLabels((chamberHistos[dtCh.rawId()])[histoTag],2);
       }
     }
     } //InOut loop
 
     // Book Theta View Related Plots
     ibooker.setCurrentFolder(topFolder(*typeIt) + "Wheel" + wheel.str() + "/Sector"
              + sector.str() + "/Station" + station.str() + "/LocalTriggerTheta");
 
     string histoTag = "";
     if((*typeIt)=="TM" && dtCh.station()!=4) {
       histoTag = (*typeIt) + "_PositionvsBX";
       chamberHistos[rawId][histoTag] = ibooker.book2D(histoTag+chTag,"Theta trigger position vs BX",
                   (int)(maxBX[(*typeIt)]-minBX[*typeIt]+1),minBX[*typeIt]-.5,maxBX[*typeIt]+.5,7,-0.5,6.5);
       histoTag = (*typeIt) + "_PositionvsQual";
       chamberHistos[rawId][histoTag] = ibooker.book2D(histoTag+chTag,"Theta trigger position vs quality",
                               2,0.5,2.5,7,-0.5,6.5);
       setQLabelsTheta(chamberHistos[rawId][histoTag],1);
       histoTag = (*typeIt) + "_ThetaBXvsQual";
       chamberHistos[rawId][histoTag] = ibooker.book2D(histoTag+chTag,"BX vs trigger quality",
                               2,0.5,2.5,(int)(maxBX[(*typeIt)]-minBX[*typeIt]+1),minBX[*typeIt]-.5,maxBX[*typeIt]+.5);
       setQLabelsTheta(chamberHistos[rawId][histoTag],1);
 //      histoTag = (*typeIt) + "_ThetaBestQual";
 //      chamberHistos[rawId][histoTag] = ibooker.book1D(histoTag+chTag,
 //                              "Trigger quality of best primitives (theta)",2,0.5,2.5); // 0 = not fired, 1 = L, 2 = H
 //      setQLabelsTheta(chamberHistos[rawId][histoTag],1);
     } else {
         if(dtCh.station()!=4){    
       histoTag = (*typeIt) + "_ThetaBXvsQual";
       chamberHistos[rawId][histoTag] =  ibooker.book2D(histoTag+chTag,"BX vs trigger quality",7,-0.5,6.5,
                        (int)(maxBX[(*typeIt)]-minBX[*typeIt]+1),minBX[*typeIt]-.5,maxBX[*typeIt]+.5);
       setQLabels((chamberHistos[dtCh.rawId()])[histoTag],1);
 
       histoTag = (*typeIt) + "_ThetaBestQual";
       chamberHistos[rawId][histoTag] = ibooker.book1D(histoTag+chTag,
       "Trigger quality of best primitives (theta)",7,-0.5,6.5);
       setQLabels((chamberHistos[dtCh.rawId()])[histoTag],1);
     }
     }
 
   }
 
   if (processTM && processDDU) {
     // Book TM/DDU Comparison Plots
     // NOt needed In and Out comparison, only IN in DDU....
     ibooker.setCurrentFolder(topFolder("DDU") + "Wheel" + wheel.str() + "/Sector"
                + sector.str() + "/Station" + station.str() + "/LocalTriggerPhiIn");
 
 
     string histoTag = "COM_QualDDUvsQualTM";
     chamberHistos[rawId][histoTag] = ibooker.book2D(histoTag+chTag,
             "DDU quality vs TM quality",8,-1.5,6.5,8,-1.5,6.5);
     setQLabels((chamberHistos[rawId])[histoTag],1);
     setQLabels((chamberHistos[rawId])[histoTag],2);
 
     histoTag = "COM_MatchingTrend";
     trendHistos[rawId] = new DTTimeEvolutionHisto(ibooker,histoTag+chTag,
                           "Fraction of DDU-TM matches w.r.t. proc evts",
                           nTimeBins,nLSTimeBin,true,0);
   }
 
 }

void DTLocalTriggerBaseTask::bookHistos	(	DQMStore::IBooker &	ibooker,
		int	wh
	)

private

Book the histograms.

Definition at line 371 of file DTLocalTriggerBaseTask.cc.

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

                                                                          {
 
   stringstream wheel; wheel << wh;
   ibooker.setCurrentFolder(topFolder("DDU") + "Wheel" + wheel.str() + "/");
   string whTag = "_W" + wheel.str();
 
   LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerBaseTask")
     << "[DTLocalTriggerBaseTask]: booking wheel histos for "
     << topFolder("DDU") << "Wheel" << wh << endl;
 
   string histoTag = "COM_BXDiff";
   MonitorElement *me = ibooker.bookProfile2D(histoTag+whTag,
         "DDU-TM BX Difference",12,1,13,4,1,5,0.,20.);
   me->setAxisTitle("Sector",1);
   me->setAxisTitle("station",2);
   wheelHistos[wh][histoTag] = me;
 
 }

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

overrideprotected

Beginrun.

Definition at line 165 of file DTLocalTriggerBaseTask.cc.

References edm::EventSetup::get(), LogTrace, and theTrigGeomUtils.

                                                                                         {
 
   LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerBaseTask")
     << "[DTLocalTriggerBaseTask]: BeginRun" << endl;
 
   ESHandle<DTGeometry> theGeom;
   context.get<MuonGeometryRecord>().get(theGeom);
   theTrigGeomUtils = new DTTrigGeomUtils(theGeom);
 
 }

void DTLocalTriggerBaseTask::endLuminosityBlock	(	const edm::LuminosityBlock &	lumiSeg,
		const edm::EventSetup &	context
	)

overrideprotected

Perform trend plot operations.

Definition at line 151 of file DTLocalTriggerBaseTask.cc.

References LogTrace, edm::LuminosityBlockBase::luminosityBlock(), nEventsInLS, and trendHistos.

                                                                                                          {
 
   LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerBaseTask")
     << "[DTLocalTriggerBaseTask]: End of LS transition" << endl;
 
   map<uint32_t,DTTimeEvolutionHisto* >::const_iterator chambIt  = trendHistos.begin();
   map<uint32_t,DTTimeEvolutionHisto* >::const_iterator chambEnd = trendHistos.end();
   for(;chambIt!=chambEnd;++chambIt) {
     chambIt->second->updateTimeSlot(lumiSeg.luminosityBlock(), nEventsInLS);
   }
 
 }

void DTLocalTriggerBaseTask::runDDUAnalysis ( edm::Handle< DTLocalTriggerCollection > & trigsDDU )

private

Run analysis on ROS data.

Definition at line 528 of file DTLocalTriggerBaseTask.cc.

References funct::abs(), bestAccRange, chamberHistos, DetId::rawId(), targetBXDDU, theCompMapIn, and tpMode.

Referenced by analyze().

                                                                                      {
 
   DTLocalTriggerCollection::DigiRangeIterator detUnitIt  = trigsDDU->begin();
   DTLocalTriggerCollection::DigiRangeIterator detUnitEnd = trigsDDU->end();
 
   for (; detUnitIt!=detUnitEnd; ++detUnitIt){
 
     const DTChamberId& chId = (*detUnitIt).first;
     uint32_t rawId = chId.rawId();
 
     const DTLocalTriggerCollection::Range& range = (*detUnitIt).second;
     DTLocalTriggerCollection::const_iterator trigIt = range.first;
     map<string, MonitorElement*> &innerME = chamberHistos[rawId];
 
     int bestQualTheta = -1;
 
     for (; trigIt!=range.second; ++trigIt){
 
       int qualPhi   = trigIt->quality();
       int qualTheta = trigIt->trTheta();
       int flag1st   = trigIt->secondTrack() ? 1 : 0;
       int bx = trigIt->bx();
       int bxPhi = bx - flag1st; // phi BX assign is different for 1st & 2nd tracks
 
       if( qualPhi>-1 && qualPhi<7 ) { // it is a phi trigger
     if (abs(bx-targetBXDDU) <= bestAccRange &&
         theCompMapIn[rawId].qualDDU()<= qualPhi)
       theCompMapIn[rawId].setDDU(qualPhi,bxPhi);
     if(tpMode) {
       innerME["DDU_BXvsQual"]->Fill(qualPhi,bxPhi); // SM BX vs Qual Phi view
     } else {
       innerME["DDU_BXvsQual"]->Fill(qualPhi,bxPhi); // SM BX vs Qual Phi view
       innerME["DDU_Flag1stvsQual"]->Fill(qualPhi,flag1st); // SM Quality vs 1st/2nd track flag Phi view
     }
       }
 
       if( qualTheta>0 && !tpMode ){// it is a theta trigger & is not TP
     if (qualTheta > bestQualTheta){
       bestQualTheta = qualTheta;
     }
     innerME["DDU_ThetaBXvsQual"]->Fill(qualTheta,bx); // SM BX vs Qual Theta view
       }
     }
 
     // Fill Quality plots with best ddu triggers
     if (!tpMode && theCompMapIn.find(rawId)!= theCompMapIn.end()) {
       int bestQualPhi = theCompMapIn[rawId].qualDDU();
       if (bestQualPhi>-1)
     innerME["DDU_BestQual"]->Fill(bestQualPhi); // SM Best Qual Trigger Phi view
       if(bestQualTheta>0) {
     innerME["DDU_ThetaBestQual"]->Fill(bestQualTheta); // SM Best Qual Trigger Theta view
       }
     }
   }
 
 }

void DTLocalTriggerBaseTask::runDDUvsTMAnalysis ( )

private

Run analysis on ROS data.

Definition at line 586 of file DTLocalTriggerBaseTask.cc.

References chamberHistos, DTTPGCompareUnit::deltaBX(), DTTPGCompareUnit::hasBoth(), DTTPGCompareUnit::hasOne(), DTTPGCompareUnit::hasSameQual(), DTTPGCompareUnit::qualDDU(), DTTPGCompareUnit::qualTM(), DTChamberId::sector(), DTChamberId::station(), theCompMapIn, trendHistos, DTChamberId::wheel(), and wheelHistos.

Referenced by analyze().

                                                {
 
   map<uint32_t,DTTPGCompareUnit>::const_iterator compIt  = theCompMapIn.begin();
   map<uint32_t,DTTPGCompareUnit>::const_iterator compEnd = theCompMapIn.end();
 
   for (; compIt!=compEnd; ++compIt) {
 
     uint32_t rawId = compIt->first;
     DTChamberId chId(rawId);
     map<string, MonitorElement*> &innerME = chamberHistos[rawId];
 
     const DTTPGCompareUnit & compUnit = compIt->second;
     if ( compUnit.hasOne() ){
       innerME["COM_QualDDUvsQualTM"]->Fill(compUnit.qualTM(),compUnit.qualDDU());
     }
     if ( compUnit.hasBoth() ){
       wheelHistos[chId.wheel()]["COM_BXDiff"]->Fill(chId.sector(),chId.station(),compUnit.deltaBX());
       if (  compUnit.hasSameQual() ) {
     trendHistos[rawId]->accumulateValueTimeSlot(1);
       }
     }
   }
 
 }

void DTLocalTriggerBaseTask::runTMAnalysis	(	std::vector< L1MuDTChambPhDigi > const *	phInTrigs,
		std::vector< L1MuDTChambPhDigi > const *	phOutTrigs,
		std::vector< L1MuDTChambThDigi > const *	thTrigs
	)

private

Run analysis on TM data.

Definition at line 391 of file DTLocalTriggerBaseTask.cc.

References funct::abs(), bestAccRange, chamberHistos, detailedAnalysis, dir, DetId::rawId(), targetBXTM, theCompMapIn, theCompMapOut, theTrigGeomUtils, tpMode, DTTrigGeomUtils::trigDir(), and DTTrigGeomUtils::trigPos().

Referenced by analyze().

 {
   vector<L1MuDTChambPhDigi>::const_iterator iph  = phInTrigs->begin();
   vector<L1MuDTChambPhDigi>::const_iterator iphe = phInTrigs->end();
 
   for(; iph !=iphe ; ++iph) {
 
     int wh    = iph->whNum();
     int sec   = iph->scNum() + 1; // DTTF->DT Convention
     int st    = iph->stNum();
     int qual  = iph->code();
     int is1st = iph->Ts2Tag() ? 1 : 0;
     int bx    = iph->bxNum() - is1st;
     int updown = iph->UpDownTag();
 
     if (qual <0 || qual>6) continue; // Check that quality is in a valid range
 
     DTChamberId dtChId(wh,st,sec);
     uint32_t rawId = dtChId.rawId();
 
     float pos = theTrigGeomUtils->trigPos(&(*iph));
     float dir = theTrigGeomUtils->trigDir(&(*iph));
 
     if (abs(bx-targetBXTM)<= bestAccRange &&
     theCompMapIn[rawId].qualTM() <= qual)
       theCompMapIn[rawId].setTM(qual,bx);
 
     map<string, MonitorElement*> &innerME = chamberHistos[rawId];
     if (tpMode) {
       innerME["TM_BXvsQual_In"]->Fill(qual,bx);      // SM BX vs Qual Phi view (1st tracks)
       innerME["TM_QualvsPhirad_In"]->Fill(pos,qual); // SM Qual vs radial angle Phi view
     } else {
       innerME["TM_BXvsQual_In"]->Fill(qual,bx);         // SM BX vs Qual Phi view (1st tracks)
       innerME["TM_Flag1stvsQual_In"]->Fill(qual,is1st); // SM Qual 1st/2nd track flag Phi view
       innerME["TM_FlagUpDownvsQual_In"]->Fill(qual,updown); // SM Qual Up/Down track flag Phi view
       if (!is1st) innerME["TM_QualvsPhirad_In"]->Fill(pos,qual);  // SM Qual vs radial angle Phi view ONLY for 1st tracks
       if (detailedAnalysis) {
     innerME["TM_QualvsPhibend_In"]->Fill(dir,qual); // SM Qual vs bending Phi view
       }
     }
 
   }
 
   iph  = phOutTrigs->begin();
   iphe = phOutTrigs->end();
 
   for(; iph !=iphe ; ++iph) {
 
     int wh    = iph->whNum();
     int sec   = iph->scNum() + 1; // DTTF->DT Convention
     int st    = iph->stNum();
     int qual  = iph->code();
     int is1st = iph->Ts2Tag() ? 1 : 0;
     int bx    = iph->bxNum() - is1st;
     int updown = iph->UpDownTag();
     if (qual <0 || qual>6) continue; // Check that quality is in a valid range
 
     DTChamberId dtChId(wh,st,sec);
     uint32_t rawId = dtChId.rawId();
 
     float pos = theTrigGeomUtils->trigPos(&(*iph));
     float dir = theTrigGeomUtils->trigDir(&(*iph));
 
     if (abs(bx-targetBXTM)<= bestAccRange &&
         theCompMapOut[rawId].qualTM() <= qual)
       theCompMapOut[rawId].setTM(qual,bx);
 
     map<string, MonitorElement*> &innerME = chamberHistos[rawId];
     if (tpMode) {
       innerME["TM_BXvsQual_Out"]->Fill(qual,bx);      // SM BX vs Qual Phi view (1st tracks)
       innerME["TM_QualvsPhirad_Out"]->Fill(pos,qual); // SM Qual vs radial angle Phi view
     } else {
       innerME["TM_BXvsQual_Out"]->Fill(qual,bx);         // SM BX vs Qual Phi view (1st tracks)
       innerME["TM_Flag1stvsQual_Out"]->Fill(qual,is1st); // SM Qual 1st/2nd track flag Phi view
       innerME["TM_FlagUpDownvsQual_Out"]->Fill(qual,updown); // SM Qual Up/Down track flag Phi view
 
       if (!is1st) innerME["TM_QualvsPhirad_Out"]->Fill(pos,qual);  // SM Qual vs radial angle Phi view ONLY for 1st tracks
       if (detailedAnalysis) {
         innerME["TM_QualvsPhibend_Out"]->Fill(dir,qual); // SM Qual vs bending Phi view
       }
     }
 
   }
 
   vector<L1MuDTChambThDigi>::const_iterator ith  = thTrigs->begin();
   vector<L1MuDTChambThDigi>::const_iterator ithe = thTrigs->end();
 
   for(; ith != ithe; ++ith) {
     int wh  = ith->whNum();
     int sec = ith->scNum() + 1; // DTTF -> DT Convention
     int st  = ith->stNum();
     int bx  = ith->bxNum();
 
     int thcode[7];
 
     for (int pos=0; pos<7; pos++)
       thcode[pos] = ith->code(pos);
 
     DTChamberId dtChId(wh,st,sec);
     uint32_t rawId = dtChId.rawId();
 
     map<string, MonitorElement*> &innerME = chamberHistos[rawId];
 
     for (int pos=0; pos<7; pos++)
       if (thcode[pos]>0) {//Fired
     innerME["TM_PositionvsBX"]->Fill(bx,pos); // SM BX vs Position Theta view
         innerME["TM_PositionvsQual"]->Fill(thcode[pos],pos); //code = pos + qual; so 0, 1, 2 for 0, L, H resp.
         innerME["TM_ThetaBXvsQual"]->Fill(thcode[pos],bx); //code = pos + qual; so 0, 1, 2 for 0, L, H resp.
       }
    }
   // Fill Quality plots with best TM triggers (phi view In)
   if (!tpMode) {
     map<uint32_t,DTTPGCompareUnit>::const_iterator compIt  = theCompMapIn.begin();
     map<uint32_t,DTTPGCompareUnit>::const_iterator compEnd = theCompMapIn.end();
     for (; compIt!=compEnd; ++compIt) {
       int bestQual = compIt->second.qualTM();
       if (bestQual > -1)
     chamberHistos[compIt->first]["TM_BestQual_In"]->Fill(bestQual);  // SM Best Qual Trigger Phi view
     }
   }
 
   // Fill Quality plots with best TM triggers (phi view Out)
        if (!tpMode) {
     map<uint32_t,DTTPGCompareUnit>::const_iterator compIt  = theCompMapOut.begin();
     map<uint32_t,DTTPGCompareUnit>::const_iterator compEnd = theCompMapOut.end();
     for (; compIt!=compEnd; ++compIt) {
       int bestQual = compIt->second.qualTM();
       if (bestQual > -1)
         chamberHistos[compIt->first]["TM_BestQual_Out"]->Fill(bestQual);  // SM Best Qual Trigger Phi view
     }
   }
 
 }

void DTLocalTriggerBaseTask::setQLabels	(	MonitorElement *	me,
		short int	iaxis
	)

private

Set Quality labels.

Definition at line 612 of file DTLocalTriggerBaseTask.cc.

References MonitorElement::getTH1(), trackerHits::histo, mps_fire::i, and tablePrinter::labels.

Referenced by bookHistos(), and topFolder().

                                                                           {
 
   TH1* histo = me->getTH1();
   if (!histo) return;
 
   TAxis* axis=nullptr;
   if (iaxis==1) {
     axis=histo->GetXaxis();
   }
   else if(iaxis==2) {
     axis=histo->GetYaxis();
   }
   if (!axis) return;
 
   string labels[7] = {"LI","LO","HI","HO","LL","HL","HH"};
   int istart = axis->GetXmin()<-1 ? 2 : 1;
   for (int i=0;i<7;i++) {
     axis->SetBinLabel(i+istart,labels[i].c_str());
   }
 
 }

void DTLocalTriggerBaseTask::setQLabelsTheta	(	MonitorElement *	me,
		short int	iaxis
	)

private

Definition at line 634 of file DTLocalTriggerBaseTask.cc.

References MonitorElement::getTH1(), trackerHits::histo, mps_fire::i, and tablePrinter::labels.

Referenced by bookHistos(), and topFolder().

                                                                                {
 
   TH1* histo = me->getTH1();
   if (!histo) return;
 
   TAxis* axis=nullptr;
   if (iaxis==1) {
     axis=histo->GetXaxis();
   }
   else if(iaxis==2) {
     axis=histo->GetYaxis();
   }
   if (!axis) return;
 
   string labels[2] = {"L","H"};
   int istart = axis->GetXmin()<-1 ? 2 : 1;
   for (int i=0;i<2;i++) {
     axis->SetBinLabel(i+istart,labels[i].c_str());
   }
 
 }