DTLocalTriggerTask Class Reference

#include <DTLocalTriggerTask.h>

Inheritance diagram for DTLocalTriggerTask:

Public Member Functions
	DTLocalTriggerTask (const edm::ParameterSet &ps)
	Constructor. More...
	~DTLocalTriggerTask () override
	Destructor. More...
Public Member Functions inherited from one::DQMEDAnalyzer< edm::one::WatchLuminosityBlocks >
	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	bookBarrelHistos (DQMStore::IBooker &, std::string histoTag)
	Book the histograms. More...
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	Book the histograms. More...
void	bookHistos (DQMStore::IBooker &, const DTChamberId &dtCh, std::string folder, std::string histoTag)
void	bookWheelHistos (DQMStore::IBooker &, int wh, std::string histoTag)
	Book the histograms. More...
void	dqmBeginRun (const edm::Run &, const edm::EventSetup &) override
	Beginrun. More...
void	endLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &context) final
void	runDDUAnalysis (edm::Handle< DTLocalTriggerCollection > &trigsDDU)
	Run analysis on ROS data. More...
void	runDDUvsTMAnalysis (std::string &trigsrc)
	Run analysis on ROS data. More...
void	runSegmentAnalysis (edm::Handle< DTRecSegment4DCollection > &segments4D)
	Run analysis using DT 4D segments. More...
void	runTMAnalysis (std::vector< L1MuDTChambPhDigi > const *phTrigs, 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 (bool isTM)
	Get the Top folder (different between Physics and TP and TM/DDU) More...
void	triggerSource (const edm::Event &e)
	Get the L1A source. More...

Protected Attributes
const int	wheelArrayShift = 3

Private Attributes
std::string	baseFolderDDU
std::string	baseFolderTM
DTArr3int	dduphcode_best
DTArr3int	dduthcode_best
bool	detailedAnalysis
std::map< uint32_t, std::map< std::string, MonitorElement * > >	digiHistos
bool	doTMTheta
DTArr3PhDigi	iphbest
DTArr3LocalTrigger	iphbestddu
bool	isLocalRun
DTArr3ThDigi	ithbest
edm::EDGetTokenT< LTCDigiCollection >	ltcDigiCollectionToken_
DTArr3mapInt	mapDTTF
edm::ESHandle< DTGeometry >	muonGeom
int	nevents
edm::ParameterSet	parameters
DTArr3int	phcode_best
edm::EDGetTokenT< DTLocalTriggerCollection >	ros_Token_
edm::EDGetTokenT< DTRecSegment4DCollection >	seg_Token_
DTArr3int	thcode_best
MonitorElement *	tm_IDDataErrorPlot
edm::EDGetTokenT< L1MuDTChambPhContainer >	tm_Token_
edm::EDGetTokenT< L1MuDTChambThContainer >	tmTh_Token_
bool	tpMode
bool	track_ok [6][5][15]
DTTrigGeomUtils *	trigGeomUtils
std::string	trigsrc
bool	useDDU
bool	useSEG
bool	useTM
std::map< int, std::map< std::string, MonitorElement * > >	wheelHistos

Friends
class	DTMonitorModule

Detailed Description

Definition at line 53 of file DTLocalTriggerTask.h.

Constructor & Destructor Documentation

DTLocalTriggerTask::DTLocalTriggerTask

(

const edm::ParameterSet &

ps

)

Constructor.

Definition at line 35 of file DTLocalTriggerTask.cc.

References baseFolderDDU, baseFolderTM, detailedAnalysis, doTMTheta, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), isLocalRun, LogTrace, ltcDigiCollectionToken_, ros_Token_, seg_Token_, tm_Token_, tmTh_Token_, and tpMode.

                                                                :
  trigGeomUtils(nullptr),
  isLocalRun(ps.getUntrackedParameter<bool>("localrun", true))
 {
  if (!isLocalRun) {
    ltcDigiCollectionToken_ = consumes<LTCDigiCollection>(
        ps.getParameter<edm::InputTag>("ltcDigiCollectionTag"));
  }

  LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerTask") << "[DTLocalTriggerTask]: Constructor"<<endl;

  tpMode           = ps.getUntrackedParameter<bool>("testPulseMode", false);
  detailedAnalysis = ps.getUntrackedParameter<bool>("detailedAnalysis", false);
  doTMTheta       = ps.getUntrackedParameter<bool>("enableTMTheta", false);
  tm_Token_       = consumes<L1MuDTChambPhContainer>(
      ps.getUntrackedParameter<InputTag>("tm_label",InputTag("twinMuxStage2Digis:PhIn")));
  tmTh_Token_       = consumes<L1MuDTChambThContainer>(
      ps.getUntrackedParameter<edm::InputTag>("tmTh_label",InputTag("twinMuxStage2Digis:ThIn")));

  ros_Token_       = consumes<DTLocalTriggerCollection>(
      ps.getUntrackedParameter<InputTag>("ros_label", InputTag("dtunpacker")));
  seg_Token_       = consumes<DTRecSegment4DCollection>(
      ps.getUntrackedParameter<InputTag>("seg_label", InputTag("dt4DSegments")));

  if (tpMode) {
    baseFolderTM = "DT/11-LocalTriggerTP-TM/";
    baseFolderDDU = "DT/12-LocalTriggerTP-DDU/";
  }
  else {
    baseFolderTM = "DT/03-LocalTrigger-TM/";
    baseFolderDDU = "DT/04-LocalTrigger-DDU/";
  }

  parameters = ps;

}

DTLocalTriggerTask::~DTLocalTriggerTask ( )

override

Destructor.

Definition at line 73 of file DTLocalTriggerTask.cc.

References LogTrace, nevents, and trigGeomUtils.

                                        {

  LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerTask") << "[DTLocalTriggerTask]: analyzed " << nevents << " events" << endl;
  if (trigGeomUtils) { delete trigGeomUtils; }

}

Member Function Documentation

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

overrideprotected

Analyze.

Definition at line 249 of file DTLocalTriggerTask.cc.

References edm::Event::getByToken(), L1MuDTChambPhContainer::getContainer(), L1MuDTChambThContainer::getContainer(), edm::HandleBase::isValid(), nevents, ros_Token_, runDDUAnalysis(), runDDUvsTMAnalysis(), runSegmentAnalysis(), runTMAnalysis(), seg_Token_, tm_Token_, tmTh_Token_, tpMode, triggerSource(), trigsrc, useDDU, useSEG, and useTM.

                                                                         {
  if (!nevents){

    edm::Handle<L1MuDTChambPhContainer> l1DTTPGPh;
    e.getByToken(tm_Token_, l1DTTPGPh);
    edm::Handle<L1MuDTChambThContainer> l1DTTPGTh;
    e.getByToken(tmTh_Token_, l1DTTPGTh); 

    useTM = (l1DTTPGPh.isValid() || l1DTTPGTh.isValid()) && parameters.getUntrackedParameter<bool>("process_tm", true) ;

    Handle<DTLocalTriggerCollection> l1DDUTrigs;
    e.getByToken(ros_Token_,l1DDUTrigs);
    useDDU = l1DDUTrigs.isValid() && parameters.getUntrackedParameter<bool>("process_ddu", true) ;

    Handle<DTRecSegment4DCollection> all4DSegments;
    e.getByToken(seg_Token_, all4DSegments);
    useSEG = all4DSegments.isValid() && parameters.getUntrackedParameter<bool>("process_seg", true) ;

  }


  nevents++;

  triggerSource(e);

  if ( useTM ) {
    edm::Handle<L1MuDTChambPhContainer> l1DTTPGPh;
    e.getByToken(tm_Token_, l1DTTPGPh);
    vector<L1MuDTChambPhDigi> const*  l1PhTrig = l1DTTPGPh->getContainer();

    edm::Handle<L1MuDTChambThContainer> l1DTTPGTh;
    e.getByToken(tmTh_Token_, l1DTTPGTh);

    vector<L1MuDTChambThDigi> const*  l1ThTrig = l1DTTPGTh->getContainer();

    runTMAnalysis(l1PhTrig, l1ThTrig);
  }
  if ( useDDU ) {
    Handle<DTLocalTriggerCollection> l1DDUTrigs;
    e.getByToken(ros_Token_, l1DDUTrigs);

    runDDUAnalysis(l1DDUTrigs);
  }
  if ( !tpMode && useSEG ) {
    Handle<DTRecSegment4DCollection> segments4D;
    e.getByToken(seg_Token_, segments4D);

    runSegmentAnalysis(segments4D);
  }
  if ( !tpMode && useTM && useDDU ) {
    runDDUvsTMAnalysis(trigsrc);
  }


}

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

overrideprotected

To reset the MEs.

Definition at line 231 of file DTLocalTriggerTask.cc.

References digiHistos, trackerHits::histo, edm::LuminosityBlockBase::id(), LogTrace, edm::LuminosityBlockID::luminosityBlock(), and genParticles_cff::map.

                                                                                                       {

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

  if(lumiSeg.id().luminosityBlock()%parameters.getUntrackedParameter<int>("ResetCycle", 3) == 0) {
    for(map<uint32_t, map<string, MonitorElement*> > ::const_iterator histo = digiHistos.begin();
    histo != digiHistos.end();
    histo++) {
      for(map<string, MonitorElement*> ::const_iterator ht = (*histo).second.begin();
      ht != (*histo).second.end();
      ht++) {
    (*ht).second->Reset();
      }
    }
  }

}

void DTLocalTriggerTask::bookBarrelHistos ( DQMStore::IBooker &  , std::string  histoTag  )

protected

Book the histograms.

Definition at line 306 of file DTLocalTriggerTask.cc.

References DQMStore::IBooker::book1D(), MonitorElement::setAxisTitle(), DQMStore::IBooker::setCurrentFolder(), tm_IDDataErrorPlot, and topFolder().

Referenced by bookHistograms().

                                                                                    {

  bool isTM = histoTag.substr(0,2) == "TM";
  ibooker.setCurrentFolder(topFolder(isTM));
  if (histoTag == "TM_ErrorsChamberID") {
    tm_IDDataErrorPlot = ibooker.book1D(histoTag.c_str(),"TM Data ID Error",5,-2,3);
    tm_IDDataErrorPlot->setAxisTitle("wheel",1);
  }

  return;

}

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

overrideprotected

Book the histograms.

Definition at line 88 of file DTLocalTriggerTask.cc.

References bookBarrelHistos(), bookHistos(), bookWheelHistos(), detailedAnalysis, doTMTheta, LogTrace, trackingPlots::stat, and tpMode.

                                                                                                                     {

  LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerTask") << "[DTLocalTriggerTask]: bookHistograms" << endl;

  if(parameters.getUntrackedParameter<bool>("staticBooking", true)) {  // Static histo booking

    vector<string> trigSources;
    if(parameters.getUntrackedParameter<bool>("localrun", true)) {
      trigSources.push_back("");
    }
    else {
      trigSources.push_back("_DTonly");
      trigSources.push_back("_NoDT");
      trigSources.push_back("_DTalso");
    }
    vector<string>::const_iterator trigSrcIt  = trigSources.begin();
    vector<string>::const_iterator trigSrcEnd = trigSources.end();

    if(parameters.getUntrackedParameter<bool>("process_tm", true)) {
      bookBarrelHistos(ibooker, "TM_ErrorsChamberID");
    }

    if (tpMode) {
      for (int stat=1;stat<5;++stat){
    for (int wh=-2;wh<3;++wh){
      for (int sect=1;sect<13;++sect){
        DTChamberId dtChId(wh,stat,sect);

        if (parameters.getUntrackedParameter<bool>("process_tm", true)){ // TM data
          bookHistos(ibooker, dtChId,"LocalTriggerPhiIn","TM_BXvsQual"+(*trigSrcIt));
          bookHistos(ibooker, dtChId,"LocalTriggerPhiIn","TM_QualvsPhirad"+(*trigSrcIt));
        }

        if (parameters.getUntrackedParameter<bool>("process_ddu", true)){ // DDU data
          bookHistos(ibooker, dtChId,"LocalTriggerPhiIn","DDU_BXvsQual"+(*trigSrcIt));
        }

      }
    }
      } // end of loop
    }
    else {
      for (;trigSrcIt!=trigSrcEnd;++trigSrcIt){
    for (int wh=-2;wh<3;++wh){
      if (parameters.getUntrackedParameter<bool>("process_tm", true) &&
          parameters.getUntrackedParameter<bool>("process_ddu", true)){ // TM+DDU data
        bookWheelHistos(ibooker, wh,"COM_BXDiff"+(*trigSrcIt));
      }
      for (int sect=1;sect<13;++sect){
        for (int stat=1;stat<5;++stat){
          DTChamberId dtChId(wh,stat,sect);
          if (parameters.getUntrackedParameter<bool>("process_tm", true)){ // TM data

        bookHistos(ibooker, dtChId,"LocalTriggerPhiIn","TM_BXvsQual_In"+(*trigSrcIt));
        if (detailedAnalysis) {
          bookHistos(ibooker, dtChId,"LocalTriggerPhiIn","TM_QualvsPhirad_In"+(*trigSrcIt));
          bookHistos(ibooker, dtChId,"LocalTriggerPhiIn","TM_QualvsPhibend_In"+(*trigSrcIt));
        }
        bookHistos(ibooker, dtChId,"LocalTriggerPhiIn","TM_Flag1stvsQual_In"+(*trigSrcIt));
                bookHistos(ibooker, dtChId,"LocalTriggerPhiIn","TM_FlagUpDownvsQual_In"+(*trigSrcIt));

        bookHistos(ibooker, dtChId,"LocalTriggerPhiIn","TM_BestQual_In"+(*trigSrcIt));
        if (stat!=4 && doTMTheta){
          bookHistos(ibooker, dtChId,"LocalTriggerTheta","TM_PositionvsBX"+(*trigSrcIt));
              bookHistos(ibooker, dtChId,"LocalTriggerTheta","TM_QualityvsBX"+(*trigSrcIt));

              bookHistos(ibooker, dtChId,"LocalTriggerTheta","TM_PositionvsQual"+(*trigSrcIt));
                
          bookHistos(ibooker, dtChId,"LocalTriggerTheta","TM_ThetaBXvsQual"+(*trigSrcIt));
          bookHistos(ibooker, dtChId,"LocalTriggerTheta","TM_ThetaBestQual"+(*trigSrcIt));
        }

        if (parameters.getUntrackedParameter<bool>("process_seg", true)){ // TM + Segemnt
          bookHistos(ibooker, dtChId,"Segment","TM_PhitkvsPhitrig"+(*trigSrcIt));
          bookHistos(ibooker, dtChId,"Segment","TM_PhibtkvsPhibtrig"+(*trigSrcIt));
          bookHistos(ibooker, dtChId,"Segment","TM_PhiResidual"+(*trigSrcIt));
          bookHistos(ibooker, dtChId,"Segment","TM_PhiResidualvsLUTPhi"+(*trigSrcIt));
          bookHistos(ibooker, dtChId,"Segment","TM_PhibResidual"+(*trigSrcIt));
          bookHistos(ibooker, dtChId,"Segment","TM_HitstkvsQualtrig"+(*trigSrcIt));
          bookHistos(ibooker, dtChId,"Segment","TM_TrackPosvsAngle"+(*trigSrcIt));
          bookHistos(ibooker, dtChId,"Segment","TM_TrackPosvsAngleandTrig"+(*trigSrcIt));
          bookHistos(ibooker, dtChId,"Segment","TM_TrackPosvsAngleandTrigHHHL"+(*trigSrcIt));
          if(stat!=4){
            bookHistos(ibooker, dtChId,"Segment","TM_TrackThetaPosvsAngle"+(*trigSrcIt)); // theta view
            bookHistos(ibooker, dtChId,"Segment","TM_TrackThetaPosvsAngleandTrig"+(*trigSrcIt));
            bookHistos(ibooker, dtChId,"Segment","TM_TrackThetaPosvsAngleandTrigH"+(*trigSrcIt));     // TM theta quality not available!
          }
        }

          }

          if (parameters.getUntrackedParameter<bool>("process_ddu", true)){ // DDU data

        bookHistos(ibooker, dtChId,"LocalTriggerPhiIn","DDU_BXvsQual"+(*trigSrcIt));
        bookHistos(ibooker, dtChId,"LocalTriggerPhiIn","DDU_Flag1stvsQual"+(*trigSrcIt));
        bookHistos(ibooker, dtChId,"LocalTriggerPhiIn","DDU_BestQual"+(*trigSrcIt));
        if(stat!=4){                                                    // theta view
          bookHistos(ibooker, dtChId,"LocalTriggerTheta","DDU_ThetaBXvsQual"+(*trigSrcIt));
          bookHistos(ibooker, dtChId,"LocalTriggerTheta","DDU_ThetaBestQual"+(*trigSrcIt));
        }

        if (parameters.getUntrackedParameter<bool>("process_seg", true)){ // DDU + Segment
          bookHistos(ibooker, dtChId,"Segment","DDU_HitstkvsQualtrig"+(*trigSrcIt));
          bookHistos(ibooker, dtChId,"Segment","DDU_TrackPosvsAngle"+(*trigSrcIt));
          bookHistos(ibooker, dtChId,"Segment","DDU_TrackPosvsAngleandTrig"+(*trigSrcIt));
          bookHistos(ibooker, dtChId,"Segment","DDU_TrackPosvsAngleandTrigHHHL"+(*trigSrcIt));
          if(stat!=4){
            bookHistos(ibooker, dtChId,"Segment","DDU_TrackThetaPosvsAngle"+(*trigSrcIt)); // theta view
            bookHistos(ibooker, dtChId,"Segment","DDU_TrackThetaPosvsAngleandTrig"+(*trigSrcIt));
            bookHistos(ibooker, dtChId,"Segment","DDU_TrackThetaPosvsAngleandTrigH"+(*trigSrcIt));
          }
        }

          }

          if (parameters.getUntrackedParameter<bool>("process_tm", true) &&
          parameters.getUntrackedParameter<bool>("process_ddu", true)){ // TM+DDU data
        bookHistos(ibooker, dtChId,"LocalTriggerPhiIn","COM_QualDDUvsQualTM"+(*trigSrcIt));
          }

        }
      }
      for (int sect=13;sect<15;++sect){
        DTChamberId dtChId(wh,4,sect);
        if (parameters.getUntrackedParameter<bool>("process_tm", true) &&
        parameters.getUntrackedParameter<bool>("process_seg", true)){ // TM+SEG LUTs data
          bookHistos(ibooker, dtChId,"Segment","TM_PhitkvsPhitrig"+(*trigSrcIt));
          bookHistos(ibooker, dtChId,"Segment","TM_PhibtkvsPhibtrig"+(*trigSrcIt));
          bookHistos(ibooker, dtChId,"Segment","TM_PhiResidual"+(*trigSrcIt));
          bookHistos(ibooker, dtChId,"Segment","TM_PhiResidualvsLUTPhi"+(*trigSrcIt));
          bookHistos(ibooker, dtChId,"Segment","TM_PhibResidual"+(*trigSrcIt));
        }
      }
    }
      }// end of loop
    }

  }


}

void DTLocalTriggerTask::bookHistos ( DQMStore::IBooker &  , const DTChamberId &  dtCh, std::string  folder, std::string  histoTag  )

protected

Definition at line 319 of file DTLocalTriggerTask.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), digiHistos, printsummarytable::folder, fftjetimagerecorder_cfi::histoLabel, createfilelist::int, LogTrace, SiStripPI::max, simTrackMatching_cfi::maxBX, min(), simTrackMatching_cfi::minBX, pileupCalc::nbins, DTTrigGeomUtils::phiRange(), DetId::rawId(), mathSSE::return(), SimDataFormats::CaloAnalysis::sc, DTChamberId::sector(), DQMStore::IBooker::setCurrentFolder(), setQLabels(), setQLabelsTheta(), DTChamberId::station(), relativeConstraints::station, DTTrigGeomUtils::thetaRange(), topFolder(), trigGeomUtils, DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

Referenced by bookHistograms().

                                                                                                                      {

  int wh=dtCh.wheel();
  int sc=dtCh.sector();
  stringstream wheel; wheel << wh;
  stringstream station; station << dtCh.station();
  stringstream sector; sector << sc;

  double minBX=0;
  double maxBX=0;
  int  rangeBX=0;

  string histoType = histoTag.substr(3,histoTag.find("_",3)-3);
  bool isTM = histoTag.substr(0,2) == "TM";

  ibooker.setCurrentFolder(topFolder(isTM) + "Wheel" + wheel.str() +
            "/Sector" + sector.str() +
            "/Station" + station.str() + "/" + folder);

  string histoName = histoTag + "_W" + wheel.str() + "_Sec" + sector.str() + "_St" + station.str();

  LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerTask") << "[DTLocalTriggerTask]: booking " << topFolder(isTM) << "Wheel" << wheel.str()
                               << "/Sector" << sector.str()
                               << "/Station"<< station.str() << "/" << folder << "/" << histoName << endl;

  if (histoType.find("BX") != string::npos){
    if (histoTag.substr(0,2) == "TM"){
      minBX= parameters.getUntrackedParameter<int>("minBXTM",-2) - 0.5;
      maxBX= parameters.getUntrackedParameter<int>("maxBXTM",2) + 0.5;
    }
    else {
      minBX= parameters.getUntrackedParameter<int>("minBXDDU",0) - 0.5;
      maxBX= parameters.getUntrackedParameter<int>("maxBXDDU",20) + 0.5;
    }
    rangeBX = (int)(maxBX-minBX);
  }

  if ( folder == "LocalTriggerPhiIn") {

    if( histoType == "BXvsQual" ){
      (digiHistos[dtCh.rawId()])[histoTag] =
    ibooker.book2D(histoName,"BX vs trigger quality",7,-0.5,6.5,rangeBX,minBX,maxBX);
      setQLabels((digiHistos[dtCh.rawId()])[histoTag],1);
      return ;
    }
    if( histoType == "BestQual" ){
      (digiHistos[dtCh.rawId()])[histoTag] =
    ibooker.book1D(histoName,"Trigger quality of best primitives",7,-0.5,6.5);
      setQLabels((digiHistos[dtCh.rawId()])[histoTag],1);
      return ;
    }
    if( histoType == "QualvsPhirad" ){
      (digiHistos[dtCh.rawId()])[histoTag] =
    ibooker.book2D(histoName,"Trigger quality vs local position",100,-500.,500.,7,-0.5,6.5);
      setQLabels((digiHistos[dtCh.rawId()])[histoTag],2);
      return ;
    }
    if( histoType == "QualvsPhibend" ) {
      (digiHistos[dtCh.rawId()])[histoTag] =
    ibooker.book2D(histoName,"Trigger quality vs local direction",200,-40.,40.,7,-0.5,6.5);
      setQLabels((digiHistos[dtCh.rawId()])[histoTag],2);
      return ;
    }
    if( histoType == "Flag1stvsQual" ) {
      (digiHistos[dtCh.rawId()])[histoTag] =
    ibooker.book2D(histoName,"1st/2nd trig flag vs quality",7,-0.5,6.5,2,-0.5,1.5);
      setQLabels((digiHistos[dtCh.rawId()])[histoTag],1);
      return ;
    }
    if( histoType == "FlagUpDownvsQual" ) {
      (digiHistos[dtCh.rawId()])[histoTag] =
        ibooker.book2D(histoName,"Up/Down trig flag vs quality",7,-0.5,6.5,2,-0.5,1.5);
      setQLabels((digiHistos[dtCh.rawId()])[histoTag],1);
      return ;
    }

    if( histoType == "QualDDUvsQualTM" ){
      (digiHistos[dtCh.rawId()])[histoTag] =
    ibooker.book2D(histoName,"DDU quality vs TM quality",8,-1.5,6.5,8,-1.5,6.5);
      setQLabels((digiHistos[dtCh.rawId()])[histoTag],1);
      setQLabels((digiHistos[dtCh.rawId()])[histoTag],2);
      return ;
    }

  }
  else if ( folder == "LocalTriggerTheta")   {

    if( histoType == "PositionvsBX" ) {
      (digiHistos[dtCh.rawId()])[histoTag] =
    ibooker.book2D(histoName,"Theta trigger position vs BX",rangeBX,minBX,maxBX,7,-0.5,6.5);
      return ;
    }
    if( histoType == "PositionvsQual" ) {
      (digiHistos[dtCh.rawId()])[histoTag] =
    ibooker.book2D(histoName,"Theta trigger position vs quality",2,0.5,2.5,7,-0.5,6.5);
      setQLabelsTheta((digiHistos[dtCh.rawId()])[histoTag],1);
      return ;
    }
    if( histoType == "ThetaBXvsQual" ) {
      (digiHistos[dtCh.rawId()])[histoTag] =
    ibooker.book2D(histoName,"BX vs trigger quality",2,0.5,2.5,rangeBX,minBX,maxBX);
      setQLabelsTheta((digiHistos[dtCh.rawId()])[histoTag],1);
    }
    if( histoType == "ThetaBestQual" ){ //It does not make sense: 7 positions, how to define best?
      (digiHistos[dtCh.rawId()])[histoTag] =
    ibooker.book1D(histoName,"Trigger quality of best primitives (theta)",2,0.5,2.5); // 0 = not fired, 1 = L, 2 = H
      setQLabelsTheta((digiHistos[dtCh.rawId()])[histoTag],1);
      return ;
    }

  }
  else if ( folder == "Segment")   {

    if( histoType.find("TrackThetaPosvsAngle" ) == 0 ) {

      string histoLabel = "Position vs Angle (theta)";
      if (histoType.find("andTrigH") != string::npos) histoLabel += " for H triggers";
      else if (histoType.find("andTrig") != string::npos) histoLabel += " for triggers";

      float min,max;
      int nbins;
      trigGeomUtils->thetaRange(dtCh,min,max,nbins);
      (digiHistos[dtCh.rawId()])[histoTag] =
    ibooker.book2D(histoName,histoLabel,16,-40.,40.,nbins,min,max);
      return ;
    }
    if( histoType.find("TrackPosvsAngle") == 0 ){

      float min,max;
      int nbins;
      trigGeomUtils->phiRange(dtCh,min,max,nbins);

      string histoLabel = "Position vs Angle (phi)";
      if (histoType.find("andTrigHHHL")  != string::npos) histoLabel += " for HH/HL triggers";
      else if (histoType.find("andTrig") != string::npos) histoLabel += " for triggers";

      (digiHistos[dtCh.rawId()])[histoTag] =
    ibooker.book2D(histoName,histoLabel,16,-40.,40.,nbins,min,max);
      return ;
    }
    if( histoType == "PhitkvsPhitrig" ){
      (digiHistos[dtCh.rawId()])[histoTag] =
    ibooker.book2D(histoName,"Local position: segment vs trigger",100,-500.,500.,100,-500.,500.);
      return ;
    }
    if( histoType == "PhibtkvsPhibtrig" ){
      (digiHistos[dtCh.rawId()])[histoTag] =
    ibooker.book2D(histoName,"Local direction : segment vs trigger",200,-40.,40.,200,-40.,40.);
      return ;
    }
    if( histoType == "PhiResidual" ){
      (digiHistos[dtCh.rawId()])[histoTag] =
    ibooker.book1D(histoName,"Trigger local position - Segment local position (correlated triggers)",400,-10.,10.);
      return ;
    }
    if( histoType == "PhibResidual" ){
      (digiHistos[dtCh.rawId()])[histoTag] =
    ibooker.book1D(histoName,"Trigger local direction - Segment local direction (correlated triggers)",500,-10.,10.);
      return ;
    }
    if( histoType == "HitstkvsQualtrig" ){
      (digiHistos[dtCh.rawId()])[histoTag] =
    ibooker.book2D(histoName,"Segment hits (phi) vs trigger quality",7,-0.5,6.5,10,0.5,10.5);
      setQLabels((digiHistos[dtCh.rawId()])[histoTag],1);
      return ;
    }

  }

}

void DTLocalTriggerTask::bookWheelHistos ( DQMStore::IBooker &  , int  wh, std::string  histoTag  )

protected

Book the histograms.

Definition at line 490 of file DTLocalTriggerTask.cc.

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

Referenced by bookHistograms().

                                                                                           {

  stringstream wheel; wheel << wh;

  string histoType = histoTag.substr(3,histoTag.find("_",3)-3);
  bool isTM = histoTag.substr(0,2) == "TM";

  ibooker.setCurrentFolder(topFolder(isTM) + "Wheel" + wheel.str() + "/");

  string histoName = histoTag + "_W" + wheel.str();

  LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerTask") << "[DTLocalTriggerTask]: booking " << topFolder(isTM)
                               << "Wheel" << wheel.str() << "/" << histoName << endl;

  if( histoType.find("BXDiff") != string::npos ){
    MonitorElement *me = ibooker.bookProfile2D(histoName,"DDU-TM BX Difference",12,1,13,4,1,5,0.,20.);
    me->setAxisTitle("Sector",1);
    me->setAxisTitle("station",2);
    (wheelHistos[wh])[histoTag] = me;
    return ;
  }

}

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

overrideprotected

Beginrun.

Definition at line 80 of file DTLocalTriggerTask.cc.

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

                                                                                    {

  nevents = 0;
  context.get<MuonGeometryRecord>().get(muonGeom);
  trigGeomUtils = new DTTrigGeomUtils(muonGeom);

}

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

inlinefinalprotected

Definition at line 102 of file DTLocalTriggerTask.h.

References MillePedeFileConverter_cfg::e, and triggerSource().

{}

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

protected

Run analysis on ROS data.

Definition at line 641 of file DTLocalTriggerTask.cc.

References dduphcode_best, dduthcode_best, digiHistos, mps_fire::i, iphbestddu, gen::k, jets_cff::quality, DetId::rawId(), tpMode, trigsrc, and wheelArrayShift.

Referenced by analyze().

                                                                                 {

  DTLocalTriggerCollection::DigiRangeIterator detUnitIt;

  for (int i=0;i<5;++i){
    for (int j=0;j<6;++j){
      for (int k=0;k<13;++k){
    dduphcode_best[j][i][k] = -1;
    dduthcode_best[j][i][k] = -1;
      }
    }
  }

  for (detUnitIt=trigsDDU->begin();
       detUnitIt!=trigsDDU->end();
       ++detUnitIt){

    const DTChamberId& id = (*detUnitIt).first;
    const DTLocalTriggerCollection::Range& range = (*detUnitIt).second;
    uint32_t indexCh = id.rawId();
    map<string, MonitorElement*> &innerME = digiHistos[indexCh];

    int wh = id.wheel();
    int sec = id.sector();
    int st = id.station();

    for (DTLocalTriggerCollection::const_iterator trigIt = range.first;
     trigIt!=range.second;
     ++trigIt){

      int bx = trigIt->bx();
      int quality = trigIt->quality();
      int thqual = trigIt->trTheta();
      int flag1st = trigIt->secondTrack() ? 1 : 0;

      // check if SC data exist: fill for any trigger
      if( quality>-1 && quality<7 ) {     // it is a phi trigger

    if(quality>dduphcode_best[wh+wheelArrayShift][st][sec]) {
      dduphcode_best[wh+wheelArrayShift][st][sec]=quality;
      iphbestddu[wh+wheelArrayShift][st][sec] = &(*trigIt);
    }

    if(tpMode) {
      innerME.find("DDU_BXvsQual"+trigsrc)->second->Fill(quality,bx-flag1st);     // SM BX vs Qual Phi view
    }
    else {
      innerME.find("DDU_BXvsQual"+trigsrc)->second->Fill(quality,bx-flag1st);     // SM BX vs Qual Phi view
      innerME.find("DDU_Flag1stvsQual"+trigsrc)->second->Fill(quality,flag1st); // SM Quality vs 1st/2nd track flag Phi view
    }
      }
      if( thqual>0 && !tpMode ) {  // it is a theta trigger

    if(thqual>dduthcode_best[wh+wheelArrayShift][st][sec] ) {
      dduthcode_best[wh+wheelArrayShift][st][sec]=thqual;
    }

    innerME.find("DDU_ThetaBXvsQual"+trigsrc)->second->Fill(thqual,bx);     // SM BX vs Qual Theta view
      }
    }

    // Fill Quality plots with best ddu triggers in phi & theta
    if (!tpMode) {
      if (dduphcode_best[wh+wheelArrayShift][st][sec]>-1 &&
      dduphcode_best[wh+wheelArrayShift][st][sec]<7){

    innerME.find("DDU_BestQual"+trigsrc)->second->Fill(dduphcode_best[wh+wheelArrayShift][st][sec]);  // Best Qual Trigger Phi view
      }
      if (dduthcode_best[wh+wheelArrayShift][st][sec]>0){
    innerME.find("DDU_ThetaBestQual"+trigsrc)->second->Fill(dduthcode_best[wh+wheelArrayShift][st][sec]); // Best Qual Trigger Theta view
      }
    }
  }
}

void DTLocalTriggerTask::runDDUvsTMAnalysis ( std::string &  trigsrc )

protected

Run analysis on ROS data.

Definition at line 899 of file DTLocalTriggerTask.cc.

References dduphcode_best, digiHistos, spr::find(), triggerObjects_cff::id, iphbest, iphbestddu, phcode_best, SimDataFormats::CaloAnalysis::sc, wheelArrayShift, and wheelHistos.

Referenced by analyze().

                                                          {

  string histoType ;
  string histoTag ;

  for (int st=1;st<5;++st){
    for (int wh=-2;wh<3;++wh){
      for (int sc=1;sc<13;++sc){
    if ( (phcode_best[wh+wheelArrayShift][st][sc]>-1 && phcode_best[wh+wheelArrayShift][st][sc]<7) ||
         (dduphcode_best[wh+wheelArrayShift][st][sc]>-1 && dduphcode_best[wh+wheelArrayShift][st][sc]<7) ){
      DTChamberId id(wh,st,sc);
      uint32_t indexCh = id.rawId();
      map<string, MonitorElement*> &innerME = digiHistos[indexCh];

      innerME.find("COM_QualDDUvsQualTM"+trigsrc)->second->Fill(phcode_best[wh+wheelArrayShift][st][sc],dduphcode_best[wh+wheelArrayShift][st][sc]);
      if ( (phcode_best[wh+wheelArrayShift][st][sc]>-1 && phcode_best[wh+wheelArrayShift][st][sc]<7) &&
           (dduphcode_best[wh+wheelArrayShift][st][sc]>-1 && dduphcode_best[wh+wheelArrayShift][st][sc]<7) ){
        int bxDDU = iphbestddu[wh+wheelArrayShift][st][sc]->bx() - iphbestddu[wh+wheelArrayShift][st][sc]->secondTrack();
        int bxTM = iphbest[wh+wheelArrayShift][st][sc]->bxNum() - iphbest[wh+wheelArrayShift][st][sc]->Ts2Tag();
        (wheelHistos[wh]).find("COM_BXDiff"+trigsrc)->second->Fill(sc,st,bxDDU-bxTM);
      }
    }
      }
    }
  }

}

void DTLocalTriggerTask::runSegmentAnalysis ( edm::Handle< DTRecSegment4DCollection > &  segments4D )

protected

Run analysis using DT 4D segments.

Definition at line 717 of file DTLocalTriggerTask.cc.

References DTTrigGeomUtils::computeSCCoordinates(), dduphcode_best, dduthcode_best, digiHistos, iphbest, phcode_best, DetId::rawId(), relativeConstraints::station, thcode_best, HiIsolationCommonParameters_cff::track, track_ok, DTTrigGeomUtils::trigDir(), trigGeomUtils, DTTrigGeomUtils::trigPos(), trigsrc, DTTrigGeomUtils::trigToSeg(), useDDU, useTM, and makeMuonMisalignmentScenario::wheel.

Referenced by analyze().

                                                                                       {

  DTRecSegment4DCollection::const_iterator track;

  // Find best tracks & good tracks
  memset(track_ok,false,450*sizeof(bool));

  DTRecSegment4DCollection::id_iterator chamberId;
  vector<const DTRecSegment4D*> best4DSegments;

  // Preliminary loop finds best 4D Segment and high quality ones
  for (chamberId = segments4D->id_begin(); chamberId != segments4D->id_end(); ++chamberId){

    DTRecSegment4DCollection::range  range = segments4D->get(*chamberId);
    const DTRecSegment4D* tmpBest=nullptr;
    int tmpdof = 0;
    int dof = 0;

    for ( track = range.first; track != range.second; ++track){

      if( (*track).hasPhi() ) {

    dof = (*track).phiSegment()->degreesOfFreedom();
    if ( dof>tmpdof ){
      tmpBest = &(*track);
      tmpdof = dof;

      int wheel = (*track).chamberId().wheel();
      int sector = (*track).chamberId().sector();
      int station = (*track).chamberId().station();
      if (sector==13){
        sector=4;
      }
      else if (sector==14){
        sector=10;
      }
      track_ok[wheel+3][station][sector] = (!track_ok[wheel+3][station][sector] && dof>=2);
    }

      }
    }
    if (tmpBest) best4DSegments.push_back(tmpBest);
  }

  vector<const DTRecSegment4D*>::const_iterator btrack;

  for ( btrack = best4DSegments.begin(); btrack != best4DSegments.end(); ++btrack ){

    if( (*btrack)->hasPhi() ) { // Phi component

      int wheel    = (*btrack)->chamberId().wheel();
      int station  = (*btrack)->chamberId().station();
      int sector   = (*btrack)->chamberId().sector();
      int scsector = 0;
      float x_track, y_track, x_angle, y_angle;
      trigGeomUtils->computeSCCoordinates((*btrack),scsector,x_track,x_angle,y_track,y_angle);
      int nHitsPhi = (*btrack)->phiSegment()->degreesOfFreedom()+2;

      DTChamberId dtChId(wheel,station,sector);      // get chamber for LUTs histograms (Sectors 1 to 14)
      uint32_t indexCh = dtChId.rawId();
      map<string, MonitorElement*> &innerMECh = digiHistos[indexCh];

      DTChamberId dtChIdSC = DTChamberId(wheel,station,scsector);  // get chamber for histograms SC granularity (sectors 1 to 12)
      indexCh = dtChIdSC.rawId();
      map<string, MonitorElement*> &innerME = digiHistos[indexCh];

      if (useDDU &&
      dduphcode_best[wheel+3][station][scsector] > -1 &&
      dduphcode_best[wheel+3][station][scsector] < 7 ) {

    // SM hits of the track vs quality of the trigger
    innerME.find("DDU_HitstkvsQualtrig"+trigsrc)->second->Fill(dduphcode_best[wheel+3][station][scsector],nHitsPhi);

      }

      if (useTM &&
      phcode_best[wheel+3][station][scsector] > -1 &&
      phcode_best[wheel+3][station][scsector] < 7 ) {

    innerME.find("TM_HitstkvsQualtrig"+trigsrc)->second->Fill(phcode_best[wheel+3][station][scsector],nHitsPhi);

    if (phcode_best[wheel+3][station][scsector]>3 && nHitsPhi>=7){

      float x_trigger     = trigGeomUtils->trigPos(iphbest[wheel+3][station][scsector]);
      float angle_trigger = trigGeomUtils->trigDir(iphbest[wheel+3][station][scsector]);
      trigGeomUtils->trigToSeg(station,x_trigger,x_angle);

      innerMECh.find("TM_PhitkvsPhitrig"+trigsrc)->second->Fill(x_trigger,x_track);
      innerMECh.find("TM_PhibtkvsPhibtrig"+trigsrc)->second->Fill(angle_trigger,x_angle);
      innerMECh.find("TM_PhiResidual"+trigsrc)->second->Fill(x_trigger-x_track);
      innerMECh.find("TM_PhibResidual"+trigsrc)->second->Fill(angle_trigger-x_angle);
    }
      }


      if (useTM) {

    // check for triggers elsewhere in the sector
    bool trigFlagTM =false;
    for (int ist=1; ist<5; ist++){
      if (ist!=station &&
          phcode_best[wheel+3][ist][scsector]>=2 &&
          phcode_best[wheel+3][ist][scsector]<7 &&
          track_ok[wheel+3][ist][scsector]==true){
        trigFlagTM = true;
        break;
      }
    }

    if (trigFlagTM && fabs(x_angle)<40. && nHitsPhi>=7){

      // position vs angle of track for reconstruced tracks (denom. for trigger efficiency)
      innerME.find("TM_TrackPosvsAngle"+trigsrc)->second->Fill(x_angle,x_track);
      if (phcode_best[wheel+3][station][scsector] >= 2 && phcode_best[wheel+3][station][scsector] < 7) {
        innerME.find("TM_TrackPosvsAngleandTrig"+trigsrc)->second->Fill(x_angle,x_track);
        if (phcode_best[wheel+3][station][scsector] > 4){  //HH & HL Triggers
          innerME.find("TM_TrackPosvsAngleandTrigHHHL"+trigsrc)->second->Fill(x_angle,x_track);
        }
      }

    }

    if ((*btrack)->hasZed() && trigFlagTM && fabs(y_angle)<40. && (*btrack)->zSegment()->degreesOfFreedom()>=1){

      // position va angle of track for reconstruced tracks (denom. for trigger efficiency) along theta direction
      innerME.find("TM_TrackThetaPosvsAngle"+trigsrc)->second->Fill(y_angle,y_track);
      if (thcode_best[wheel+3][station][scsector] > 0) {
        innerME.find("TM_TrackThetaPosvsAngleandTrig"+trigsrc)->second->Fill(y_angle,y_track);
              if (thcode_best[wheel+3][station][scsector] == 2) {
                innerME.find("TM_TrackThetaPosvsAngleandTrigH"+trigsrc)->second->Fill(y_angle,y_track);
              }
      }

    }
      }

      if (useDDU) {

    // check for triggers elsewhere in the sector
    bool trigFlagDDU =false;
    for (int ist=1; ist<5; ist++){
      if (ist!=station &&
          dduphcode_best[wheel+3][ist][scsector]>=2 &&
          dduphcode_best[wheel+3][ist][scsector]<7 &&
          track_ok[wheel+3][ist][scsector]==true){
        trigFlagDDU = true;
        break;
      }
    }

    if (trigFlagDDU && fabs(x_angle)<40. && nHitsPhi>=7){

      // position vs angle of track for reconstruced tracks (denom. for trigger efficiency)
      innerME.find("DDU_TrackPosvsAngle"+trigsrc)->second->Fill(x_angle,x_track);
      if (dduphcode_best[wheel+3][station][scsector] >= 2 && dduphcode_best[wheel+3][station][scsector] < 7) {
      innerME.find("DDU_TrackPosvsAngleandTrig"+trigsrc)->second->Fill(x_angle,x_track);
        if (dduphcode_best[wheel+3][station][scsector] > 4){  //HH & HL Triggers
          innerME.find("DDU_TrackPosvsAngleandTrigHHHL"+trigsrc)->second->Fill(x_angle,x_track);
        }
      }

    }

    if ((*btrack)->hasZed() && trigFlagDDU && fabs(y_angle)<40. && (*btrack)->zSegment()->degreesOfFreedom()>=1){

      // position va angle of track for reconstruced tracks (denom. for trigger efficiency) along theta direction
      innerME.find("DDU_TrackThetaPosvsAngle"+trigsrc)->second->Fill(y_angle,y_track);
      if (dduthcode_best[wheel+3][station][scsector] > 0) {
        innerME.find("DDU_TrackThetaPosvsAngleandTrig"+trigsrc)->second->Fill(y_angle,y_track);
        if (dduthcode_best[wheel+3][station][scsector] == 3) {
          innerME.find("DDU_TrackThetaPosvsAngleandTrigH"+trigsrc)->second->Fill(y_angle,y_track);
        }
      }

    }
      }
    }
  }

}

void DTLocalTriggerTask::runTMAnalysis ( std::vector< L1MuDTChambPhDigi > const *  phTrigs, std::vector< L1MuDTChambThDigi > const *  thTrigs  )

protected

Run analysis on TM data.

Definition at line 514 of file DTLocalTriggerTask.cc.

References angle(), detailedAnalysis, digiHistos, doTMTheta, MonitorElement::Fill(), mps_fire::i, triggerObjects_cff::id, iphbest, ithbest, gen::k, phcode_best, DetId::rawId(), SimDataFormats::CaloAnalysis::sc, thcode_best, tm_IDDataErrorPlot, tpMode, DTTrigGeomUtils::trigDir(), trigGeomUtils, DTTrigGeomUtils::trigPos(), trigsrc, wheelArrayShift, and x.

Referenced by analyze().

                                                                   {

                               //exit(0);

  string histoType ;
  string histoTag ;


  // define best quality trigger segment (phi and theta)
  // in any station start from 1 and zero is kept empty
  for (int i=0;i<5;++i){
    for (int j=0;j<6;++j){
      for (int k=0;k<13;++k){
    phcode_best[j][i][k] = -1;
    thcode_best[j][i][k] = -1;
      }
    }
  }

  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; // SM The track finder goes from 0 to 11. I need them from 1 to 12 !!!!!
    int phst    = iph->stNum();
    int phbx    = iph->bxNum();
    int phcode  = iph->code();
    int phi1st  = iph->Ts2Tag();
    int updown  = iph->UpDownTag();

    // FIXME: workaround for TM data with station ID
    if(phst == 0) {
      tm_IDDataErrorPlot->Fill(phwheel);
      continue;
    }

    if(phcode>phcode_best[phwheel+3][phst][phsec] && phcode<7) {
      phcode_best[phwheel+3][phst][phsec]=phcode;
      iphbest[phwheel+3][phst][phsec] = &(*iph);
    }

    DTChamberId dtChId(phwheel,phst,phsec);

    float x     = trigGeomUtils->trigPos(&(*iph));
    float angle = trigGeomUtils->trigDir(&(*iph));
    uint32_t indexCh = dtChId.rawId();

    map<string, MonitorElement*> &innerME = digiHistos[indexCh];

    if (tpMode) {
      innerME.find("TM_BXvsQual"+trigsrc)->second->Fill(phcode,phbx-phi1st);    // SM BX vs Qual Phi view (1st tracks)
      innerME.find("TM_QualvsPhirad"+trigsrc)->second->Fill(x,phcode);          // SM Qual vs radial angle Phi view
    }
    else {
      innerME.find("TM_BXvsQual_In"+trigsrc)->second->Fill(phcode,phbx-phi1st);    // SM BX vs Qual Phi view (1st tracks)
      innerME.find("TM_Flag1stvsQual_In"+trigsrc)->second->Fill(phcode,phi1st);    // SM Qual 1st/2nd track flag Phi view
      innerME.find("TM_FlagUpDownvsQual_In"+trigsrc)->second->Fill(phcode,updown);    // SM Qual Up/Down track flag Phi view
      if (detailedAnalysis) {
    innerME.find("TM_QualvsPhirad_In"+trigsrc)->second->Fill(x,phcode);          // SM Qual vs radial angle Phi view
    innerME.find("TM_QualvsPhibend_In"+trigsrc)->second->Fill(angle,phcode);     // SM Qual vs bending Phi view
      }
    }

  }

  if (doTMTheta) {
    int thcode[7];
    vector<L1MuDTChambThDigi>::const_iterator ith  = thTrigs->begin();
    vector<L1MuDTChambThDigi>::const_iterator ithe = thTrigs->end();
    for(; ith != ithe; ++ith) {
      int thwheel = ith->whNum();
      int thsec   = ith->scNum() + 1; // SM The track finder goes from 0 to 11. I need them from 1 to 12 !!!!!
      int thst    = ith->stNum();
      int thbx    = ith->bxNum();

      for (int pos=0; pos<7; pos++) {
    thcode[pos] = ith->code(pos);

    if(thcode[pos]>thcode_best[thwheel+3][thst][thsec] ) {
      thcode_best[thwheel+3][thst][thsec]=thcode[pos];
      ithbest[thwheel+3][thst][thsec] = &(*ith);
    }
      }

      DTChamberId dtChId(thwheel,thst,thsec);
      uint32_t indexCh = dtChId.rawId();

      map<string, MonitorElement*> &innerME = digiHistos[indexCh];

      for (int pos=0; pos<7; pos++) { //SM fill position for non zero position bit in theta view
    if(thcode[pos]>0){ //Fired
      innerME.find("TM_PositionvsBX"+trigsrc)->second->Fill(thbx,pos);          // SM BX vs Position Theta view
          innerME.find("TM_PositionvsQual"+trigsrc)->second->Fill(thcode[pos],pos); //code = pos + qual; so 0, 1, 2 for 0, L, H resp.
          innerME.find("TM_ThetaBXvsQual"+trigsrc)->second->Fill(thcode[pos],thbx); //code = pos + qual; so 0, 1, 2 for 0, L, H resp.
    }
      }
    }
  }


  // Fill Quality plots with best TM triggers in phi & theta
  if (!tpMode) {
    for (int st=1;st<5;++st){
      for (int wh=-2;wh<3;++wh){
    for (int sc=1;sc<13;++sc){
      if (phcode_best[wh+wheelArrayShift][st][sc]>-1 && phcode_best[wh+wheelArrayShift][st][sc]<7){
        DTChamberId id(wh,st,sc);
        uint32_t indexCh = id.rawId();
        map<string, MonitorElement*> &innerME = digiHistos[indexCh];

        innerME.find("TM_BestQual_In"+trigsrc)->second->Fill(phcode_best[wh+wheelArrayShift][st][sc]);  // Best Qual Trigger Phi view
      }
          if (thcode_best[wh+wheelArrayShift][st][sc]>0 && thcode_best[wh+wheelArrayShift][st][sc]<3){
            DTChamberId id(wh,st,sc);
            uint32_t indexCh = id.rawId();
            map<string, MonitorElement*> &innerME = digiHistos[indexCh];

            innerME.find("TM_ThetaBestQual"+trigsrc)->second->Fill(thcode_best[wh+wheelArrayShift][st][sc]);  // Best Qual Trigger Theta view 
      }
    }
      }
    }
  }
}

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

protected

Set Quality labels.

Definition at line 927 of file DTLocalTriggerTask.cc.

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

Referenced by bookHistos().

                                                                      {

  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 DTLocalTriggerTask::setQLabelsTheta ( MonitorElement *  me, short int  iaxis  )

protected

Definition at line 949 of file DTLocalTriggerTask.cc.

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

Referenced by bookHistos().

                                                                            {

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

}

std::string& DTLocalTriggerTask::topFolder ( bool  isTM )

inlineprotected

Get the Top folder (different between Physics and TP and TM/DDU)

Definition at line 108 of file DTLocalTriggerTask.h.

References baseFolderDDU, and baseFolderTM.

Referenced by bookBarrelHistos(), bookHistos(), and bookWheelHistos().

{ return isTM ? baseFolderTM : baseFolderDDU; }

void DTLocalTriggerTask::triggerSource ( const edm::Event &  e )

protected

Get the L1A source.

Definition at line 973 of file DTLocalTriggerTask.cc.

References edm::Event::getByToken(), mps_fire::i, isLocalRun, ltcDigiCollectionToken_, and trigsrc.

Referenced by analyze(), and endLuminosityBlock().

                                                        {


  if (!isLocalRun){

    Handle<LTCDigiCollection> ltcdigis;
    e.getByToken(ltcDigiCollectionToken_, ltcdigis);

    for (std::vector<LTCDigi>::const_iterator ltc_it = ltcdigis->begin(); ltc_it != ltcdigis->end(); ltc_it++){

      size_t otherTriggerSum=0;
      for (size_t i = 1; i < 6; i++) {
    otherTriggerSum += size_t((*ltc_it).HasTriggered(i));
      }
      if ((*ltc_it).HasTriggered(0) && otherTriggerSum == 0)
    trigsrc = "_DTonly";
      else if (!(*ltc_it).HasTriggered(0))
    trigsrc = "_NoDT";
      else if ((*ltc_it).HasTriggered(0) && otherTriggerSum > 0)
    trigsrc = "_DTalso";

    }
    return;
  }

  trigsrc = "";
  return;

}

Friends And Related Function Documentation

friend class DTMonitorModule

friend

Definition at line 55 of file DTLocalTriggerTask.h.

Member Data Documentation

std::string DTLocalTriggerTask::baseFolderDDU

private

Definition at line 125 of file DTLocalTriggerTask.h.

Referenced by DTLocalTriggerTask(), and topFolder().

std::string DTLocalTriggerTask::baseFolderTM

private

Definition at line 124 of file DTLocalTriggerTask.h.

Referenced by DTLocalTriggerTask(), and topFolder().

DTArr3int DTLocalTriggerTask::dduphcode_best

private

Definition at line 131 of file DTLocalTriggerTask.h.

Referenced by runDDUAnalysis(), runDDUvsTMAnalysis(), and runSegmentAnalysis().

DTArr3int DTLocalTriggerTask::dduthcode_best

private

Definition at line 133 of file DTLocalTriggerTask.h.

Referenced by runDDUAnalysis(), and runSegmentAnalysis().

bool DTLocalTriggerTask::detailedAnalysis

private

Definition at line 127 of file DTLocalTriggerTask.h.

Referenced by bookHistograms(), DTLocalTriggerTask(), and runTMAnalysis().

std::map<uint32_t, std::map<std::string, MonitorElement*> > DTLocalTriggerTask::digiHistos

private

Definition at line 143 of file DTLocalTriggerTask.h.

Referenced by beginLuminosityBlock(), bookHistos(), runDDUAnalysis(), runDDUvsTMAnalysis(), runSegmentAnalysis(), and runTMAnalysis().

bool DTLocalTriggerTask::doTMTheta

private

Definition at line 126 of file DTLocalTriggerTask.h.

Referenced by bookHistograms(), DTLocalTriggerTask(), and runTMAnalysis().

DTArr3PhDigi DTLocalTriggerTask::iphbest

private

Definition at line 135 of file DTLocalTriggerTask.h.

Referenced by runDDUvsTMAnalysis(), runSegmentAnalysis(), and runTMAnalysis().

DTArr3LocalTrigger DTLocalTriggerTask::iphbestddu

private

Definition at line 136 of file DTLocalTriggerTask.h.

Referenced by runDDUAnalysis(), and runDDUvsTMAnalysis().

bool DTLocalTriggerTask::isLocalRun

private

Definition at line 148 of file DTLocalTriggerTask.h.

Referenced by DTLocalTriggerTask(), and triggerSource().

DTArr3ThDigi DTLocalTriggerTask::ithbest

private

Definition at line 137 of file DTLocalTriggerTask.h.

Referenced by runTMAnalysis().

edm::EDGetTokenT<LTCDigiCollection> DTLocalTriggerTask::ltcDigiCollectionToken_

private

Definition at line 118 of file DTLocalTriggerTask.h.

Referenced by DTLocalTriggerTask(), and triggerSource().

DTArr3mapInt DTLocalTriggerTask::mapDTTF

private

Definition at line 134 of file DTLocalTriggerTask.h.

edm::ESHandle<DTGeometry> DTLocalTriggerTask::muonGeom

private

Definition at line 141 of file DTLocalTriggerTask.h.

Referenced by dqmBeginRun().

int DTLocalTriggerTask::nevents

private

Definition at line 122 of file DTLocalTriggerTask.h.

Referenced by analyze(), dqmBeginRun(), and ~DTLocalTriggerTask().

edm::ParameterSet DTLocalTriggerTask::parameters

private

Definition at line 140 of file DTLocalTriggerTask.h.

Referenced by Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputTags(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::recursePSetProperties().

DTArr3int DTLocalTriggerTask::phcode_best

private

Definition at line 130 of file DTLocalTriggerTask.h.

Referenced by runDDUvsTMAnalysis(), runSegmentAnalysis(), and runTMAnalysis().

edm::EDGetTokenT<DTLocalTriggerCollection> DTLocalTriggerTask::ros_Token_

private

Definition at line 116 of file DTLocalTriggerTask.h.

Referenced by analyze(), and DTLocalTriggerTask().

edm::EDGetTokenT<DTRecSegment4DCollection> DTLocalTriggerTask::seg_Token_

private

Definition at line 117 of file DTLocalTriggerTask.h.

Referenced by analyze(), and DTLocalTriggerTask().

DTArr3int DTLocalTriggerTask::thcode_best

private

Definition at line 132 of file DTLocalTriggerTask.h.

Referenced by runSegmentAnalysis(), and runTMAnalysis().

MonitorElement* DTLocalTriggerTask::tm_IDDataErrorPlot

private

Definition at line 146 of file DTLocalTriggerTask.h.

Referenced by bookBarrelHistos(), and runTMAnalysis().

edm::EDGetTokenT<L1MuDTChambPhContainer> DTLocalTriggerTask::tm_Token_

private

Definition at line 114 of file DTLocalTriggerTask.h.

Referenced by analyze(), and DTLocalTriggerTask().

edm::EDGetTokenT<L1MuDTChambThContainer> DTLocalTriggerTask::tmTh_Token_

private

Definition at line 115 of file DTLocalTriggerTask.h.

Referenced by analyze(), and DTLocalTriggerTask().

bool DTLocalTriggerTask::tpMode

private

Definition at line 123 of file DTLocalTriggerTask.h.

Referenced by analyze(), bookHistograms(), DTLocalTriggerTask(), runDDUAnalysis(), and runTMAnalysis().

bool DTLocalTriggerTask::track_ok[6][5][15]

private

Definition at line 138 of file DTLocalTriggerTask.h.

Referenced by runSegmentAnalysis().

DTTrigGeomUtils* DTLocalTriggerTask::trigGeomUtils

private

Definition at line 142 of file DTLocalTriggerTask.h.

Referenced by bookHistos(), dqmBeginRun(), runSegmentAnalysis(), runTMAnalysis(), and ~DTLocalTriggerTask().

std::string DTLocalTriggerTask::trigsrc

private

Definition at line 121 of file DTLocalTriggerTask.h.

Referenced by analyze(), runDDUAnalysis(), runSegmentAnalysis(), runTMAnalysis(), and triggerSource().

bool DTLocalTriggerTask::useDDU

private

Definition at line 120 of file DTLocalTriggerTask.h.

Referenced by analyze(), and runSegmentAnalysis().

bool DTLocalTriggerTask::useSEG

private

Definition at line 120 of file DTLocalTriggerTask.h.

Referenced by analyze().

bool DTLocalTriggerTask::useTM

private

Definition at line 120 of file DTLocalTriggerTask.h.

Referenced by analyze(), and runSegmentAnalysis().

const int DTLocalTriggerTask::wheelArrayShift = 3

protected

Definition at line 110 of file DTLocalTriggerTask.h.

Referenced by runDDUAnalysis(), runDDUvsTMAnalysis(), and runTMAnalysis().

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

private

Definition at line 144 of file DTLocalTriggerTask.h.

Referenced by bookWheelHistos(), and runDDUvsTMAnalysis().