#include <DTLocalTriggerSynchTask.h>

Inheritance diagram for DTLocalTriggerSynchTask:

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

	~DTLocalTriggerSynchTask () 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 &event, const edm::EventSetup &context) override
	Analyze. More...

std::string &	baseDir ()

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	Book the histograms. More...

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

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

Protected Attributes
const int	wheelArrayShift = 3

Private Attributes
float	angleRange

std::string	baseDirectory

float	bxTime

edm::EDGetTokenT< DTLocalTriggerCollection >	ddu_Token_

int	fineDelay

float	minHitsPhi

edm::ESHandle< DTGeometry >	muonGeom

int	nBXHigh

int	nBXLow

int	nevents

DTArr3int	phCodeBestDDU

DTArr3int	phCodeBestTM

DTArr4int	phCodeBXTM

bool	processDDU

bool	rangeInBX

edm::EDGetTokenT< DTRecSegment4DCollection >	seg_Token_

DTArr3int	segHitBest

DTArr3int	thCodeBestDDU

MonitorElement *	tm_IDDataErrorPlot

edm::EDGetTokenT< L1MuDTChambPhContainer >	tm_Token_

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

std::unique_ptr< DTTTrigBaseSync >	tTrigSync

Friends
class	DTMonitorModule

Detailed Description

Definition at line 51 of file DTLocalTriggerSynchTask.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 43 of file DTLocalTriggerSynchTask.cc.

References angleRange, baseDirectory, bxTime, ddu_Token_, reco::get(), edm::ParameterSet::getParameter(), minHitsPhi, nBXHigh, nBXLow, processDDU, rangeInBX, seg_Token_, AlCaHLTBitMon_QueryRunRegistry::string, and tm_Token_.

                                                                           : nevents(0),
   tTrigSync{DTTTrigSyncFactory::get()->create(ps.getParameter<std::string>("tTrigMode"),
                                               ps.getParameter<edm::ParameterSet>("tTrigModeConfig"))}
 {
 
 
   edm::LogVerbatim ("DTLocalTriggerSynchTask")  << "[DTLocalTriggerSynchTask]: Constructor" << endl;
   tm_Token_  = consumes<L1MuDTChambPhContainer>(
       ps.getParameter<edm::InputTag>("TMInputTag"));
   seg_Token_  = consumes<DTRecSegment4DCollection>(
       ps.getParameter<edm::InputTag>("SEGInputTag"));
 
   processDDU = ps.getUntrackedParameter<bool>("processDDU",false);
 
   if (processDDU)
     ddu_Token_  = consumes<DTLocalTriggerCollection>(
          ps.getParameter<edm::InputTag>("DDUInputTag"));
 
   bxTime        = ps.getParameter<double>("bxTimeInterval");   // CB move this to static const or DB
   rangeInBX     = ps.getParameter<bool>("rangeWithinBX");
   nBXLow        = ps.getParameter<int>("nBXLow");
   nBXHigh       = ps.getParameter<int>("nBXHigh");
   angleRange    = ps.getParameter<double>("angleRange");
   minHitsPhi    = ps.getParameter<int>("minHitsPhi");
   baseDirectory = ps.getParameter<string>("baseDir");
 }

DTLocalTriggerSynchTask::~DTLocalTriggerSynchTask ( )

override

Destructor.

Definition at line 71 of file DTLocalTriggerSynchTask.cc.

References nevents.

                                                   {
 
   edm::LogVerbatim ("DTLocalTriggerSynchTask")  << "[DTLocalTriggerSynchTask]: analyzed " << nevents << " events" << endl;
 
 }

Member Function Documentation

void DTLocalTriggerSynchTask::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	context
	)

overrideprotected

Analyze.

Definition at line 106 of file DTLocalTriggerSynchTask.cc.

References angleRange, bxTime, ddu_Token_, DTRecSegment2D::degreesOfFreedom(), dir, L1MuDTChambPhContainer::getContainer(), hfClusterShapes_cfi::hits, mps_fire::i, createfilelist::int, DTRecSegment2D::ist0Valid(), gen::k, minHitsPhi, nevents, phCodeBestDDU, phCodeBestTM, Geom::pi(), processDDU, jets_cff::quality, rangeInBX, DetId::rawId(), DTChamberId::sector(), seg_Token_, DTChamberId::station(), relativeConstraints::station, ntuplemaker::time, tm_Token_, HiIsolationCommonParameters_cff::track, triggerHistos, tTrigSync, DTChamberId::wheel(), makeMuonMisalignmentScenario::wheel, and wheelArrayShift.

                                                                                         {
 
   nevents++;
 
   for (int i=0;i<5;++i){
     for (int j=0;j<6;++j){
       for (int k=0;k<13;++k){
     phCodeBestTM[j][i][k] = -1;
     phCodeBestDDU[j][i][k] = -1;
       }
     }
   }
 
   // Get best TM triggers
   edm::Handle<L1MuDTChambPhContainer> l1DTTPGPh;
   event.getByToken(tm_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[0-11] -> DT[1-12] Sector Numbering
     int phst    = iph->stNum();
     int phcode  = iph->code();
 
     if(phcode>phCodeBestTM[phwheel+3][phst][phsec] && phcode<7) {
       phCodeBestTM[phwheel+3][phst][phsec]=phcode;
     }
 
   }
 
   // Get best DDU triggers
   if (processDDU){
     Handle<DTLocalTriggerCollection> trigsDDU;
     event.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;
 
       int wh = id.wheel();
       int sec = id.sector();
       int st = id.station();
 
       for (DTLocalTriggerCollection::const_iterator trigIt = range.first; trigIt!=range.second;++trigIt){
     
     int quality = trigIt->quality();
     
     if(quality>-1 && quality<7 &&
        quality>phCodeBestDDU[wh+wheelArrayShift][st][sec]) {
       phCodeBestDDU[wh+wheelArrayShift][st][sec]=quality;
     }
       }
     }
   }
 
   //Get best segments (highest number of phi hits)
   vector<const DTRecSegment4D*> bestSegments4D;
   Handle<DTRecSegment4DCollection> segments4D;
   event.getByToken(seg_Token_, segments4D);
   DTRecSegment4DCollection::const_iterator track;
   DTRecSegment4DCollection::id_iterator chambIdIt;
 
   for (chambIdIt = segments4D->id_begin(); chambIdIt != segments4D->id_end(); ++chambIdIt){
 
     DTRecSegment4DCollection::range  range = segments4D->get(*chambIdIt);
     const DTRecSegment4D* best=nullptr;
     int hitsBest = 0;
     int hits = 0;
 
     for ( track = range.first; track != range.second; ++track){
       if( (*track).hasPhi() ) {
     hits = (*track).phiSegment()->degreesOfFreedom()+2;
     if ( hits>hitsBest ){
       best = &(*track);
       hitsBest = hits;
     }
       }
     }
     if (best) {
       bestSegments4D.push_back(best);
     }
   }
 
 
   // Filling histos
   vector<const DTRecSegment4D*>::const_iterator bestSegIt  = bestSegments4D.begin();
   vector<const DTRecSegment4D*>::const_iterator bestSegEnd = bestSegments4D.end();
   for (; bestSegIt!=bestSegEnd; ++bestSegIt ){
 
     float dir = atan((*bestSegIt)->localDirection().x()/ (*bestSegIt)->localDirection().z())*180/Geom::pi(); // CB cerca un modo migliore x farlo
     const DTRecSegment2D* seg2D = (*bestSegIt)->phiSegment();
     int nHitsPhi = seg2D->degreesOfFreedom()+2;
     DTChamberId chambId = (*bestSegIt)->chamberId();
     map<string, MonitorElement*> &innerME = triggerHistos[chambId.rawId()];
 
     if (fabs(dir)<angleRange &&
     nHitsPhi>=minHitsPhi &&
     seg2D->ist0Valid()){
 
       float t0seg = (*bestSegIt)->phiSegment()->t0();
       float tTrig = (tTrigSync->offset(DTWireId(chambId,1,1,2)) + tTrigSync->offset(DTWireId(chambId,3,1,2)) )/2;
       float time = tTrig+t0seg;
       float htime = rangeInBX ? time-int(time/bxTime)*bxTime : time-int(tTrig/bxTime)*bxTime;
 
       int wheel   = chambId.wheel();
       int sector  = chambId.sector();
       int station = chambId.station();
       int scsector = sector>12 ? sector==13 ? 4 : 10 : sector;
 
       int qualTM = phCodeBestTM[wheel+3][station][scsector];
       int qualDDU = phCodeBestDDU[wheel+3][station][scsector];
 
       if (fabs(t0seg)>0.01) {
     innerME.find("SEG_TrackCrossingTime")->second->Fill(htime);
     if ( qualTM>=0 ) innerME.find("TM_TrackCrossingTimeAll")->second->Fill(htime);
     if ( qualTM==6 ) innerME.find("TM_TrackCrossingTimeHH")->second->Fill(htime);
     if ( processDDU && qualDDU>=0 ) innerME.find("DDU_TrackCrossingTimeAll")->second->Fill(htime);
     if ( processDDU && qualDDU==6 ) innerME.find("DDU_TrackCrossingTimeHH")->second->Fill(htime);
       }
 
     }
   }
 
 }

std::string& DTLocalTriggerSynchTask::baseDir ( )

inlineprotected

Definition at line 77 of file DTLocalTriggerSynchTask.h.

References baseDirectory.

Referenced by bookHistograms(), and bookHistos().

77 { return baseDirectory; }

DTLocalTriggerSynchTask::baseDirectory

std::string baseDirectory

Definition: DTLocalTriggerSynchTask.h:101

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

overrideprotected

Book the histograms.

Definition at line 77 of file DTLocalTriggerSynchTask.cc.

References baseDir(), DQMStore::IBooker::bookFloat(), bookHistos(), bxTime, DTGeometry::chambers(), MonitorElement::Fill(), muonGeom, DQMStore::IBooker::setCurrentFolder(), triggerHistos, and tTrigSync.

                                                                                                                           {
 
   edm::LogVerbatim ("DTLocalTriggerSynchTask") <<"[DTLocalTriggerSynchTask]: Book Histograms"<<endl;
 
   ibooker.setCurrentFolder(baseDir());
   ibooker.bookFloat("BXTimeSpacing")->Fill(bxTime);
 
   tTrigSync->setES(context);
 
 
   std::vector<const DTChamber*>::const_iterator chambIt  = muonGeom->chambers().begin();
   std::vector<const DTChamber*>::const_iterator chambEnd = muonGeom->chambers().end();
 
   for (; chambIt!=chambEnd; ++chambIt) {
     bookHistos(ibooker,(*chambIt)->id());
     triggerHistos[(*chambIt)->id().rawId()]["tTrig_SL1"]->Fill(tTrigSync->offset(DTWireId((*chambIt)->id(),1,1,2)));
     triggerHistos[(*chambIt)->id().rawId()]["tTrig_SL3"]->Fill(tTrigSync->offset(DTWireId((*chambIt)->id(),3,1,2)));
   }
 
 
 }

void DTLocalTriggerSynchTask::bookHistos	(	DQMStore::IBooker &	ibooker,
		const DTChamberId &	dtCh
	)

protected

Book the histograms.

Definition at line 236 of file DTLocalTriggerSynchTask.cc.

References baseDir(), DQMStore::IBooker::book1D(), DQMStore::IBooker::bookFloat(), bxTime, SiStripPI::max, min(), pileupCalc::nbins, nBXHigh, nBXLow, processDDU, rangeInBX, DetId::rawId(), DTChamberId::sector(), DQMStore::IBooker::setCurrentFolder(), DTChamberId::station(), relativeConstraints::station, triggerHistos, DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

Referenced by bookHistograms().

                                                                                              {
 
   stringstream wheel; wheel << dtChId.wheel();
   stringstream station; station << dtChId.station();
   stringstream sector; sector << dtChId.sector();
   uint32_t chRawId = dtChId.rawId();
 
   ibooker.setCurrentFolder(baseDir() + "/Wheel" + wheel.str() + "/Sector" + sector.str() + "/Station" + station.str() );
 
   std::vector<string> histoTags = { "SEG_TrackCrossingTime", "TM_TrackCrossingTimeAll", "TM_TrackCrossingTimeHH" };
 
   if (processDDU) {
     histoTags.push_back("DDU_TrackCrossingTimeAll");
     histoTags.push_back("DDU_TrackCrossingTimeHH");
   }
 
   float min = rangeInBX ?      0 : nBXLow*bxTime;
   float max = rangeInBX ? bxTime : nBXHigh*bxTime;
   int nbins = static_cast<int>(ceil( rangeInBX ? bxTime : (nBXHigh-nBXLow)*bxTime));
 
   for (const auto & histoTag : histoTags) {
     string histoName = histoTag + (rangeInBX ? "InBX" : "") + "_W" + wheel.str() + "_Sec" + sector.str() + "_St" + station.str();
     edm::LogVerbatim ("DTLocalTriggerSynchTask") << "[DTLocalTriggerSynchTask]: booking "
                          << baseDir() + "/Wheel" << wheel.str()
                         << "/Sector" << sector.str()
                         << "/Station"<< station.str()
                         << "/" << histoName << endl;
 
     triggerHistos[chRawId][histoTag] = ibooker.book1D(histoName.c_str(),"Track time distribution",nbins,min,max);
   }
 
   string floatTag[2] = { "tTrig_SL1", "tTrig_SL3" };
 
   for (int iFloat=0;iFloat<2;++iFloat) {
     string floatName = floatTag[iFloat] + "_W" + wheel.str() + "_Sec" + sector.str() + "_St" + station.str();
     triggerHistos[chRawId][floatTag[iFloat]] = ibooker.bookFloat(floatName);
   }
 
 }

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

overrideprotected

Beginrun.

Definition at line 100 of file DTLocalTriggerSynchTask.cc.

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

                                                                                    {
 
   context.get<MuonGeometryRecord>().get(muonGeom);
 }

Friends And Related Function Documentation

friend class DTMonitorModule

friend

Definition at line 53 of file DTLocalTriggerSynchTask.h.

Member Data Documentation

float DTLocalTriggerSynchTask::angleRange

private

Definition at line 96 of file DTLocalTriggerSynchTask.h.

Referenced by analyze(), and DTLocalTriggerSynchTask().

std::string DTLocalTriggerSynchTask::baseDirectory

private

Definition at line 101 of file DTLocalTriggerSynchTask.h.

Referenced by baseDir(), and DTLocalTriggerSynchTask().

float DTLocalTriggerSynchTask::bxTime

private

Definition at line 91 of file DTLocalTriggerSynchTask.h.

Referenced by analyze(), bookHistograms(), bookHistos(), and DTLocalTriggerSynchTask().

edm::EDGetTokenT<DTLocalTriggerCollection> DTLocalTriggerSynchTask::ddu_Token_

private

Definition at line 108 of file DTLocalTriggerSynchTask.h.

Referenced by analyze(), and DTLocalTriggerSynchTask().

int DTLocalTriggerSynchTask::fineDelay

private

Definition at line 98 of file DTLocalTriggerSynchTask.h.

float DTLocalTriggerSynchTask::minHitsPhi

private

Definition at line 97 of file DTLocalTriggerSynchTask.h.

Referenced by analyze(), and DTLocalTriggerSynchTask().

edm::ESHandle<DTGeometry> DTLocalTriggerSynchTask::muonGeom

private

Definition at line 103 of file DTLocalTriggerSynchTask.h.

Referenced by bookHistograms(), and dqmBeginRun().

int DTLocalTriggerSynchTask::nBXHigh

private

Definition at line 95 of file DTLocalTriggerSynchTask.h.

Referenced by bookHistos(), and DTLocalTriggerSynchTask().

int DTLocalTriggerSynchTask::nBXLow

private

Definition at line 94 of file DTLocalTriggerSynchTask.h.

Referenced by bookHistos(), and DTLocalTriggerSynchTask().

int DTLocalTriggerSynchTask::nevents

private

Definition at line 83 of file DTLocalTriggerSynchTask.h.

Referenced by analyze(), and ~DTLocalTriggerSynchTask().

DTArr3int DTLocalTriggerSynchTask::phCodeBestDDU

private

Definition at line 87 of file DTLocalTriggerSynchTask.h.

Referenced by analyze().

DTArr3int DTLocalTriggerSynchTask::phCodeBestTM

private

Definition at line 85 of file DTLocalTriggerSynchTask.h.

Referenced by analyze().

DTArr4int DTLocalTriggerSynchTask::phCodeBXTM

private

Definition at line 86 of file DTLocalTriggerSynchTask.h.

bool DTLocalTriggerSynchTask::processDDU

private

Definition at line 93 of file DTLocalTriggerSynchTask.h.

Referenced by analyze(), bookHistos(), and DTLocalTriggerSynchTask().

bool DTLocalTriggerSynchTask::rangeInBX

private

Definition at line 92 of file DTLocalTriggerSynchTask.h.

Referenced by analyze(), bookHistos(), and DTLocalTriggerSynchTask().

edm::EDGetTokenT<DTRecSegment4DCollection> DTLocalTriggerSynchTask::seg_Token_

private

Definition at line 109 of file DTLocalTriggerSynchTask.h.

Referenced by analyze(), and DTLocalTriggerSynchTask().

DTArr3int DTLocalTriggerSynchTask::segHitBest

private

Definition at line 89 of file DTLocalTriggerSynchTask.h.

DTArr3int DTLocalTriggerSynchTask::thCodeBestDDU

private

Definition at line 88 of file DTLocalTriggerSynchTask.h.

MonitorElement* DTLocalTriggerSynchTask::tm_IDDataErrorPlot

private

Definition at line 105 of file DTLocalTriggerSynchTask.h.

edm::EDGetTokenT<L1MuDTChambPhContainer> DTLocalTriggerSynchTask::tm_Token_

private

Definition at line 107 of file DTLocalTriggerSynchTask.h.

Referenced by analyze(), and DTLocalTriggerSynchTask().

std::map<uint32_t, std::map<std::string, MonitorElement*> > DTLocalTriggerSynchTask::triggerHistos

private

Definition at line 104 of file DTLocalTriggerSynchTask.h.

Referenced by analyze(), bookHistograms(), and bookHistos().

std::unique_ptr<DTTTrigBaseSync> DTLocalTriggerSynchTask::tTrigSync

private

Definition at line 99 of file DTLocalTriggerSynchTask.h.

Referenced by analyze(), and bookHistograms().

const int DTLocalTriggerSynchTask::wheelArrayShift = 3

protected

Definition at line 79 of file DTLocalTriggerSynchTask.h.

Referenced by analyze().

Public Member Functions

Protected Member Functions

Protected Attributes

Private Attributes

Friends

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Friends And Related Function Documentation

Member Data Documentation