#include <DTChamberEfficiencyTask.h>

Inheritance diagram for DTChamberEfficiencyTask:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &setup) override

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

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

	DTChamberEfficiencyTask (const edm::ParameterSet &pset)
	Constructor. More...

	~DTChamberEfficiencyTask () 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	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

Private Member Functions
void	bookHistos (DQMStore::IBooker &ibooker, DTChamberId chId)

const DTRecSegment4D &	getBestSegment (const DTRecSegment4DCollection::range &segs) const

const DTRecSegment4D *	getBestSegment (const DTRecSegment4D s1, const DTRecSegment4D s2) const

LocalPoint	interpolate (const DTRecSegment4D &seg1, const DTRecSegment4D &seg3, const DTChamberId &MB2) const

bool	isGoodSegment (const DTRecSegment4D &seg) const

Private Attributes
bool	debug

bool	detailedAnalysis

edm::ESHandle< DTGeometry >	dtGeom

std::map< DTChamberId, std::vector< MonitorElement * > >	histosPerCh

bool	onlineMonitor

edm::ParameterSet	parameters

edm::EDGetTokenT< DTRecSegment4DCollection >	recHits4DToken_

edm::Handle< DTRecSegment4DCollection >	segs

double	theMinChi2NormSegment

double	theMinCloseDist

unsigned int	theMinHitsSegment

Detailed Description

DQM Analysis of 4D DT segments, it produces plots about:

single chamber efficiency All histos are produced per Chamber

Class based on the code written by S. Lacaprara : RecoLocalMuon / DTSegment / test / DTEffAnalyzer.h

Author: G. Mila - INFN Torino

Definition at line 43 of file DTChamberEfficiencyTask.h.

Constructor & Destructor Documentation

DTChamberEfficiencyTask::DTChamberEfficiencyTask ( const edm::ParameterSet & pset )

Constructor.

Definition at line 32 of file DTChamberEfficiencyTask.cc.

References debug, edm::ParameterSet::getUntrackedParameter(), and muonDTDigis_cfi::pset.

                                                                          {
 
   debug = pset.getUntrackedParameter<bool>("debug",false);
 
   edm::LogVerbatim ("DTDQM|DTMonitorModule|DTChamberEfficiencyTask") << "[DTChamberEfficiencyTask] Constructor called!";
 
   parameters = pset;
 
   // the name of the 4D rec hits collection
   recHits4DToken_ = consumes<DTRecSegment4DCollection>(
       edm::InputTag(parameters.getParameter<string>("recHits4DLabel")));
 
   // parameters to use for the segment quality check
   theMinHitsSegment = static_cast<unsigned int>(parameters.getParameter<int>("minHitsSegment"));
   theMinChi2NormSegment = parameters.getParameter<double>("minChi2NormSegment");
   // parameter to use for the exstrapolated segment check
   theMinCloseDist = parameters.getParameter<double>("minCloseDist");
 
   // the running modality
   onlineMonitor = parameters.getUntrackedParameter<bool>("onlineMonitor");
 
   // the analysis mode
   detailedAnalysis = parameters.getUntrackedParameter<bool>("detailedAnalysis");
 
 }

DTChamberEfficiencyTask::~DTChamberEfficiencyTask ( )

override

Destructor.

Definition at line 59 of file DTChamberEfficiencyTask.cc.

                                                  {
 
   edm::LogVerbatim ("DTDQM|DTMonitorModule|DTChamberEfficiencyTask") << "[DTChamberEfficiencyTask] Destructor called!";
 }

Member Function Documentation

void DTChamberEfficiencyTask::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

override

Definition at line 153 of file DTChamberEfficiencyTask.cc.

References DTRecSegment4D::hasPhi(), DTRecSegment4D::hasZed(), plotFactory::histos, HiCaloJetParameters_cff::interpolate, DTRecSegment4D::localPosition(), mag(), DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, DTChamberId::wheel(), makeMuonMisalignmentScenario::wheel, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

                                                                                        {
 
   edm::LogVerbatim ("DTDQM|DTMonitorModule|DTChamberEfficiencyTask") << "[DTChamberEfficiencyTask] Analyze #Run: " << event.id().run()
      << " #Event: " << event.id().event();
 
   // Get the 4D rechit collection from the event
   event.getByToken(recHits4DToken_, segs);
 
   int bottom=0, top=0;
 
 
   // Loop over all the chambers
   vector<const DTChamber*>::const_iterator ch_it = dtGeom->chambers().begin();
   vector<const DTChamber*>::const_iterator ch_end = dtGeom->chambers().end();
   for (; ch_it != ch_end; ++ch_it) {
 
     DTChamberId ch = (*ch_it)->id();
     int wheel =  ch.wheel();
     int sector = ch.sector();
     int station = ch.station();
 
 
     DTChamberId MidId(wheel, station, sector);
 
     // get efficiency for MB1 using MB2 and MB3
     if( station == 1 ) {
       bottom = 2;
       top = 3;
     }
 
     // get efficiency for MB2 using MB1 and MB3
     if( station == 2 ) {
       bottom = 1;
       top = 3;
     }
 
     // get efficiency for MB2 using MB2 and MB4
     if( station == 3 ) {
       bottom = 2;
       top = 4;
     }
 
     // get efficiency for MB4 using MB2 and MB3
     if( station == 4 ) {
       bottom = 2;
       top = 3;
     }
 
     // Select events with (good) segments in Bot and Top
     DTChamberId BotId(wheel, bottom, sector);
     DTChamberId TopId(wheel, top, sector);
 
     // Get segments in the bottom chambers (if any)
     DTRecSegment4DCollection::range segsBot= segs->get(BotId);
     int nSegsBot=segsBot.second-segsBot.first;
     // check if any segments is there
     if (nSegsBot==0) continue;
 
     vector<MonitorElement *> histos =  histosPerCh[MidId];
 
     // Get segments in the top chambers (if any)
     DTRecSegment4DCollection::range segsTop= segs->get(TopId);
     int nSegsTop=segsTop.second-segsTop.first;
 
     // Select one segment for the bottom chamber
     const DTRecSegment4D& bestBotSeg= getBestSegment(segsBot);
 
     // Select one segment for the top chamber
     DTRecSegment4D* pBestTopSeg=nullptr;
     if (nSegsTop>0)
       pBestTopSeg = const_cast<DTRecSegment4D*>(&getBestSegment(segsTop));
     //if top chamber is MB4 sector 10, consider also sector 14
     if (TopId.station() == 4 && TopId.sector() == 10) {
       DTChamberId TopId14(wheel, top, 14);
       DTRecSegment4DCollection::range segsTop14= segs->get(TopId14);
       int nSegsTop14=segsTop14.second-segsTop14.first;
       nSegsTop+=nSegsTop14;
       if (nSegsTop14) {
     DTRecSegment4D* pBestTopSeg14 = const_cast<DTRecSegment4D*>(&getBestSegment(segsTop14));
     // get best between sector 10 and 14
     pBestTopSeg = const_cast<DTRecSegment4D*>(getBestSegment(pBestTopSeg, pBestTopSeg14));
       }
     }
     if (!pBestTopSeg) continue;
     const DTRecSegment4D& bestTopSeg= *pBestTopSeg;
 
     DTRecSegment4DCollection::range segsMid= segs->get(MidId);
     int nSegsMid=segsMid.second-segsMid.first;
 
     if(detailedAnalysis){
       // very trivial efficiency, just count segments
       histos[3]->Fill(0);
       if (nSegsMid>0) histos[3]->Fill(1);
     }
 
     // get position at Mid by interpolating the position (not direction) of best
     // segment in Bot and Top to Mid surface
     LocalPoint posAtMid = interpolate(bestBotSeg, bestTopSeg, MidId);
 
     // is best segment good enough?
     if (isGoodSegment(bestBotSeg) && isGoodSegment(bestTopSeg)) {
       if(detailedAnalysis)
     histos[4]->Fill(posAtMid.x(),posAtMid.y());
       //check if interpolation fall inside middle chamber
       if ((dtGeom->chamber(MidId))->surface().bounds().inside(posAtMid)) {
     histos[0]->Fill(posAtMid.x(),posAtMid.y());
     if (nSegsMid>0) {
 
       if(detailedAnalysis){
         histos[3]->Fill(2);
         histos[5]->Fill(posAtMid.x(),posAtMid.y());
       }
 
       const DTRecSegment4D& bestMidSeg= getBestSegment(segsMid);
       // check if middle segments is good enough
       if (isGoodSegment(bestMidSeg)) {
 
         if(detailedAnalysis)
           histos[6]->Fill(posAtMid.x(),posAtMid.y());
         LocalPoint midSegPos=bestMidSeg.localPosition();
 
         // check if middle segments is also close enough
         double dist;
         if (bestMidSeg.hasPhi()) {
           if (bestTopSeg.hasZed() && bestBotSeg.hasZed() && bestMidSeg.hasZed()) {
         dist = (midSegPos-posAtMid).mag();
           } else {
             dist = fabs((midSegPos-posAtMid).x());
           }
         } else {
           dist = fabs((midSegPos-posAtMid).y());
         }
         if (dist < theMinCloseDist ) {
           histos[1]->Fill(posAtMid.x(),posAtMid.y());
         }
         if(detailedAnalysis)
           histos[2]->Fill(dist);
       }
     }
       }
     }
   }// loop over stations
 
 }

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

override

To reset the MEs.

Definition at line 64 of file DTChamberEfficiencyTask.cc.

References trackerHits::histo, mps_fire::i, edm::LuminosityBlockBase::id(), edm::LuminosityBlockID::luminosityBlock(), genParticles_cff::map, and findQualityFiles::size.

                                                                                                             {
 
   edm::LogVerbatim ("DTDQM|DTMonitorModule|DTChamberEfficiencyTask")<<"[DTChamberEfficiencyTask]: Begin of LS transition";
 
   if(lumiSeg.id().luminosityBlock()%parameters.getUntrackedParameter<int>("ResetCycle", 3) == 0 && onlineMonitor) {
     for(map<DTChamberId, vector<MonitorElement*> > ::const_iterator histo = histosPerCh.begin();
     histo != histosPerCh.end();
     histo++) {
       int size = (*histo).second.size();
       for(int i=0; i<size; i++){
     (*histo).second[i]->Reset();
       }
     }
   }
 
 }

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

overrideprotected

Definition at line 89 of file DTChamberEfficiencyTask.cc.

References bookHistos(), and DQMStore::IBooker::setCurrentFolder().

                                                                                                                           {
 
   ibooker.setCurrentFolder("DT/DTChamberEfficiencyTask");
 
   // Loop over all the chambers
   vector<const DTChamber*>::const_iterator ch_it = dtGeom->chambers().begin();
   vector<const DTChamber*>::const_iterator ch_end = dtGeom->chambers().end();
   for (; ch_it != ch_end; ++ch_it) {
     // histo booking
     bookHistos(ibooker, (*ch_it)->id());
   }
 
 }

void DTChamberEfficiencyTask::bookHistos	(	DQMStore::IBooker &	ibooker,
		DTChamberId	chId
	)

private

Definition at line 104 of file DTChamberEfficiencyTask.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), plotFactory::histos, DTChamberId::sector(), DQMStore::IBooker::setCurrentFolder(), DTChamberId::station(), relativeConstraints::station, DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

                                                                                     {
 
   edm::LogVerbatim ("DTDQM|DTMonitorModule|DTChamberEfficiencyTask") << "   Booking histos for CH : " << chId;
 
   // Compose the chamber name
   stringstream wheel; wheel << chId.wheel();
   stringstream station; station << chId.station();
   stringstream sector; sector << chId.sector();
 
   string HistoName =
     "_W" + wheel.str() +
     "_St" + station.str() +
     "_Sec" + sector.str();
 
   ibooker.setCurrentFolder("DT/01-DTChamberEfficiency/Task/Wheel" + wheel.str() +
                "/Sector" + sector.str() +
                            "/Station" + station.str());
 
   // Create the monitor elements
   vector<MonitorElement *> histos;
 
   //efficiency selection cuts
   // a- number of segments of the top chamber > 0 && number of segments of the bottom chamber > 0
   // b- number of segments of the middle chamber > 0
   // c- check of the top and bottom segment quality
   // d- check if interpolation falls inside the middle chamber
   // e- check of the middle segment quality
   // f- check if the distance between the reconstructed and the exstrapolated segments is ok
 
 
   // histo for efficiency with cuts a-/c-/d-
   histos.push_back(ibooker.book2D("hEffGoodSegVsPosDen"+HistoName,"Eff vs local position (good) ",25,-250.,250., 25,-250.,250.));
   // histo for efficiency with cuts a-/b-/c-/d-/e-/f-
   histos.push_back(ibooker.book2D("hEffGoodCloseSegVsPosNum"+HistoName, "Eff vs local position (good and close segs) ", 25,-250.,250., 25,-250.,250.));
   if(detailedAnalysis){
     histos.push_back(ibooker.book1D("hDistSegFromExtrap"+HistoName, "Distance segments from extrap position ",200,0.,200.));
     // histo for efficiency from segment counting
     histos.push_back(ibooker.book1D("hNaiveEffSeg"+HistoName, "Naive eff ",10,0.,10.));
     // histo for efficiency with cuts a-/c-
   histos.push_back(ibooker.book2D("hEffSegVsPosDen"+HistoName,"Eff vs local position (all) ",25,-250.,250., 25,-250.,250.));
     // histo for efficiency with cuts a-/b-/c-/d-
     histos.push_back(ibooker.book2D("hEffSegVsPosNum"+HistoName, "Eff vs local position ",25,-250.,250., 25,-250.,250.));
     // histo for efficiency with cuts a-/b-/c-/d-/e-
     histos.push_back(ibooker.book2D("hEffGoodSegVsPosNum"+HistoName, "Eff vs local position (good segs) ", 25,-250.,250., 25,-250.,250.));
   }
   histosPerCh[chId] = histos;
 }

void DTChamberEfficiencyTask::dqmBeginRun	(	const edm::Run &	run,
		const edm::EventSetup &	setup
	)

override

BeginRun.

Definition at line 82 of file DTChamberEfficiencyTask.cc.

References edm::EventSetup::get().

                                                                                       {
 
   // Get the DT Geometry
   setup.get<MuonGeometryRecord>().get(dtGeom);
 
 }

const DTRecSegment4D & DTChamberEfficiencyTask::getBestSegment ( const DTRecSegment4DCollection::range & segs ) const

private

Definition at line 302 of file DTChamberEfficiencyTask.cc.

                                                                                                             {
   DTRecSegment4DCollection::const_iterator bestIter;
   unsigned int nHitBest=0;
   double chi2Best=99999.;
   for (DTRecSegment4DCollection::const_iterator seg=segs.first ;
        seg!=segs.second ; ++seg ) {
     unsigned int nHits= ((*seg).hasPhi() ? (*seg).phiSegment()->recHits().size() : 0 ) ;
     nHits+= ((*seg).hasZed() ?  (*seg).zSegment()->recHits().size() : 0 );
 
     if (nHits==nHitBest) {
       if ((*seg).chi2()/(*seg).degreesOfFreedom() < chi2Best ) {
         chi2Best=(*seg).chi2()/(*seg).degreesOfFreedom();
         bestIter = seg;
       }
     }
     else if (nHits>nHitBest) {
       nHitBest=nHits;
       bestIter = seg;
     }
   }
   return *bestIter;
 }

const DTRecSegment4D * DTChamberEfficiencyTask::getBestSegment	(	const DTRecSegment4D *	s1,
		const DTRecSegment4D *	s2
	)		const

private

Definition at line 325 of file DTChamberEfficiencyTask.cc.

References DTRecSegment4D::chi2(), DTRecSegment4D::degreesOfFreedom(), DTRecSegment4D::hasPhi(), DTRecSegment4D::hasZed(), DTRecSegment4D::phiSegment(), DTRecSegment2D::recHits(), indexGen::s2, and DTRecSegment4D::zSegment().

                                                                                    {
 
   if (!s1) return s2;
   if (!s2) return s1;
   unsigned int nHits1= (s1->hasPhi() ? s1->phiSegment()->recHits().size() : 0 ) ;
   nHits1+= (s1->hasZed() ?  s1->zSegment()->recHits().size() : 0 );
 
   unsigned int nHits2= (s2->hasPhi() ? s2->phiSegment()->recHits().size() : 0 ) ;
   nHits2+= (s2->hasZed() ?  s2->zSegment()->recHits().size() : 0 );
 
   if (nHits1==nHits2) {
     if (s1->chi2()/s1->degreesOfFreedom() < s2->chi2()/s2->degreesOfFreedom() )
       return s1;
     else
       return s2;
   }
   else if (nHits1>nHits2) return s1;
   return s2;
 }

LocalPoint DTChamberEfficiencyTask::interpolate	(	const DTRecSegment4D &	seg1,
		const DTRecSegment4D &	seg3,
		const DTChamberId &	MB2
	)		const

private

Definition at line 347 of file DTChamberEfficiencyTask.cc.

References DTRecSegment4D::chamberId(), dir, DTRecSegment4D::hasPhi(), DTRecSegment4D::hasZed(), DTRecSegment4D::localPosition(), toLocal(), csvLumiCalc::unit, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                       {
   // Get GlobalPoition of Seg in MB1
   GlobalPoint gpos1=(dtGeom->chamber(seg1.chamberId()))->toGlobal(seg1.localPosition());
 
   // Get GlobalPoition of Seg in MB3
   GlobalPoint gpos3=(dtGeom->chamber(seg3.chamberId()))->toGlobal(seg3.localPosition());
 
   // interpolate
   // get all in MB2 frame
   LocalPoint pos1=(dtGeom->chamber(id2))->toLocal(gpos1);
   LocalPoint pos3=(dtGeom->chamber(id2))->toLocal(gpos3);
 
   // case 1: 1 and 3 has both projection. No problem
 
   // case 2: one projection is missing for one of the segments. Keep the other's segment position
   if (!seg1.hasZed()) pos1=LocalPoint(pos1.x(),pos3.y(),pos1.z());
   if (!seg3.hasZed()) pos3=LocalPoint(pos3.x(),pos1.y(),pos3.z());
 
   if (!seg1.hasPhi()) pos1=LocalPoint(pos3.x(),pos1.y(),pos1.z());
   if (!seg3.hasPhi()) pos3=LocalPoint(pos1.x(),pos3.y(),pos3.z());
 
   // direction
   LocalVector dir = (pos3-pos1).unit(); // z points inward!
   LocalPoint pos2 = pos1+dir*pos1.z()/(-dir.z());
 
   return pos2;
 }

bool DTChamberEfficiencyTask::isGoodSegment ( const DTRecSegment4D & seg ) const

private

Definition at line 378 of file DTChamberEfficiencyTask.cc.

References DTRecSegment4D::chamberId(), DTRecSegment4D::chi2(), DTRecSegment4D::degreesOfFreedom(), DTRecSegment4D::hasPhi(), DTRecSegment4D::hasZed(), DTRecSegment4D::phiSegment(), DTRecSegment2D::recHits(), DTChamberId::station(), and DTRecSegment4D::zSegment().

                                                                            {
   if (seg.chamberId().station()!=4 && !seg.hasZed()) return false;
   unsigned int nHits= (seg.hasPhi() ? seg.phiSegment()->recHits().size() : 0 ) ;
   nHits+= (seg.hasZed() ?  seg.zSegment()->recHits().size() : 0 );
   return ( nHits >= theMinHitsSegment &&
        seg.chi2()/seg.degreesOfFreedom() < theMinChi2NormSegment );
 }

Member Data Documentation

bool DTChamberEfficiencyTask::debug

private

Definition at line 77 of file DTChamberEfficiencyTask.h.

Referenced by util.rrapi.RRApi::dprint(), rrapi.RRApi::dprint(), pkg.AbstractPkg::generate(), util.rrapi.RRApi::get(), rrapi.RRApi::get(), pkg.AbstractPkg::get_kwds(), and pkg.AbstractPkg::write().

bool DTChamberEfficiencyTask::detailedAnalysis

private

Definition at line 81 of file DTChamberEfficiencyTask.h.

edm::ESHandle<DTGeometry> DTChamberEfficiencyTask::dtGeom

private

Definition at line 94 of file DTChamberEfficiencyTask.h.

std::map<DTChamberId, std::vector<MonitorElement*> > DTChamberEfficiencyTask::histosPerCh

private

Definition at line 88 of file DTChamberEfficiencyTask.h.

bool DTChamberEfficiencyTask::onlineMonitor

private

Definition at line 79 of file DTChamberEfficiencyTask.h.

edm::ParameterSet DTChamberEfficiencyTask::parameters

private

Definition at line 86 of file DTChamberEfficiencyTask.h.

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

edm::EDGetTokenT<DTRecSegment4DCollection> DTChamberEfficiencyTask::recHits4DToken_

private

Definition at line 84 of file DTChamberEfficiencyTask.h.

edm::Handle<DTRecSegment4DCollection> DTChamberEfficiencyTask::segs

private

Definition at line 95 of file DTChamberEfficiencyTask.h.

double DTChamberEfficiencyTask::theMinChi2NormSegment

private

Definition at line 91 of file DTChamberEfficiencyTask.h.

double DTChamberEfficiencyTask::theMinCloseDist

private

Definition at line 92 of file DTChamberEfficiencyTask.h.

unsigned int DTChamberEfficiencyTask::theMinHitsSegment

private

Definition at line 90 of file DTChamberEfficiencyTask.h.

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation