---

DTEfficiencyTask Class Reference

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

• single cell efficiency All histos are produced per Layer.

More...

#include <DQM/DTMonitorModule/src/DTEfficiencyTask.h>

Inheritance diagram for DTEfficiencyTask:

List of all members.

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &setup)
void	beginJob (const edm::EventSetup &c)
	BeginJob.
void	beginLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &context)
	To reset the MEs.
	DTEfficiencyTask (const edm::ParameterSet &pset)
	Constructor.
void	endJob ()
	Endjob.
virtual	~DTEfficiencyTask ()
	Destructor.
Private Member Functions
void	bookHistos (DTLayerId lId, int fisrtWire, int lastWire)
void	fillHistos (DTLayerId lId, int firstWire, int lastWire, int missingWire, bool UnassHit)
void	fillHistos (DTLayerId lId, int firstWire, int lastWire, int numWire)
Private Attributes
bool	debug
std::map< DTLayerId, std::vector< MonitorElement * > >	histosPerL
edm::ParameterSet	parameters
DQMStore *	theDbe
std::string	theRecHitLabel
std::string	theRecHits4DLabel

---

Detailed Description

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

• single cell efficiency All histos are produced per Layer.

Date

2008/03/01 00:39:54

Revision

1.6

Author:

G. Mila - INFN Torino

Definition at line 31 of file DTEfficiencyTask.h.

---

Constructor & Destructor Documentation

DTEfficiencyTask::DTEfficiencyTask

(

const edm::ParameterSet &

pset

)

Constructor.

Definition at line 43 of file DTEfficiencyTask.cc.

References debug, edm::ParameterSet::getUntrackedParameter(), parameters, and theDbe.

                                                           {

  debug = pset.getUntrackedParameter<bool>("debug","false");

  // Get the DQM needed services
  theDbe = edm::Service<DQMStore>().operator->();
  theDbe->setCurrentFolder("DT/DTEfficiencyTask");

  parameters = pset;
}

DTEfficiencyTask::~DTEfficiencyTask

(

)

[virtual]

Destructor.

Definition at line 55 of file DTEfficiencyTask.cc.

                                   {
}  

---

Member Function Documentation

void DTEfficiencyTask::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)			[virtual]

Implements edm::EDAnalyzer.

Definition at line 89 of file DTEfficiencyTask.cc.

References edmNew::copy(), funct::cos(), GenMuonPlsPt100GeV_cfg::cout, debug, lat::endl(), fillHistos(), edm::EventSetup::get(), i, iter, DTChamber::layer(), DTWireId::layerId(), DTRangeMapAccessor::layersBySuperLayer(), range, DTRecSegment2D::specificRecHits(), DTSuperLayerId::superlayer(), DTLayerId::superlayerId(), DTChamber::superLayers(), theRecHitLabel, theRecHits4DLabel, PV3DBase< T, PVType, FrameType >::theta(), GeomDet::toGlobal(), GeomDet::toLocal(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                {


  if(debug)
    cout << "[DTEfficiencyTask] Analyze #Run: " << event.id().run()
         << " #Event: " << event.id().event() << endl;


  // Get the 4D segment collection from the event
  edm::Handle<DTRecSegment4DCollection> all4DSegments;
  event.getByLabel(theRecHits4DLabel, all4DSegments);

  // Get the rechit collection from the event
  Handle<DTRecHitCollection> dtRecHits;
  event.getByLabel(theRecHitLabel, dtRecHits);

  // Get the DT Geometry
  ESHandle<DTGeometry> dtGeom;
  setup.get<MuonGeometryRecord>().get(dtGeom);

  // Loop over all chambers containing a segment
  DTRecSegment4DCollection::id_iterator chamberId;
  for (chamberId = all4DSegments->id_begin();
       chamberId != all4DSegments->id_end();
       ++chamberId) {

    // Get the chamber
    const DTChamber* chamber = dtGeom->chamber(*chamberId); 

    // Get all 1D RecHits to be used for searches of hits not associated to segments and map them by wire
    const vector<const DTSuperLayer*> SLayers = chamber->superLayers();
    map<DTWireId, int> wireAnd1DRecHits;
    for(vector<const DTSuperLayer*>::const_iterator superlayer = SLayers.begin();
        superlayer != SLayers.end();
        superlayer++) {
        DTRecHitCollection::range  range = dtRecHits->get(DTRangeMapAccessor::layersBySuperLayer((*superlayer)->id()));
        // Loop over the rechits of this ChamberId
        for (DTRecHitCollection::const_iterator rechit = range.first;
             rechit!=range.second;
             ++rechit){
          wireAnd1DRecHits[(*rechit).wireId()] = (*rechit).wireId().wire();
        }
      }


    // Get the 4D segment range for the corresponding ChamerId
    DTRecSegment4DCollection::range  range = all4DSegments->get(*chamberId);
    int nsegm = distance(range.first, range.second);
    if(debug)
      cout << "   Chamber: " << *chamberId << " has " << nsegm
           << " 4D segments" << endl;

    
    // Loop over the rechits of this ChamerId
    for (DTRecSegment4DCollection::const_iterator segment4D = range.first;
         segment4D!=range.second;
         ++segment4D) {
      if(debug)
        cout << "   == RecSegment dimension: " << (*segment4D).dimension() << endl;

      // If Statio != 4 skip RecHits with dimension != 4
      // For the Station 4 consider 2D RecHits
      if((*chamberId).station() != 4 && (*segment4D).dimension() != 4) {
        if(debug)
          cout << "[DTEfficiencyTask]***Warning: RecSegment dimension is not 4 but "
               << (*segment4D).dimension() << ", skipping!" << endl;
        continue;
      } else if((*chamberId).station() == 4 && (*segment4D).dimension() != 2) {
        if(debug)
          cout << "[DTEfficiencyTask]***Warning: RecSegment dimension is not 2 but "
               << (*segment4D).dimension() << ", skipping!" << endl;
        continue;
      }

      vector<DTRecHit1D> recHits1D;   
      bool rPhi = false;
      bool rZ = false;

      // Get 1D RecHits and select only events with 7 or 8 hits in phi and 3 or 4 hits in theta (if any)
      const DTChamberRecSegment2D* phiSeg = (*segment4D).phiSegment();
      vector<DTRecHit1D> phiRecHits = phiSeg->specificRecHits();
      rPhi = true;
      if(phiRecHits.size() < 7 || phiRecHits.size() > 8 ) {
        if(debug)
          cout << "[DTEfficiencyTask] Phi segments has: " << phiRecHits.size()
               << " hits, skipping" << endl; // FIXME: info output
        continue;
      }
      copy(phiRecHits.begin(), phiRecHits.end(), back_inserter(recHits1D));
      const DTSLRecSegment2D* zSeg;
      if((*segment4D).dimension() == 4) {
        rZ = true;
        zSeg = (*segment4D).zSegment();
        vector<DTRecHit1D> zRecHits = zSeg->specificRecHits();
        if(zRecHits.size() < 3 || zRecHits.size() > 4 ) {
          if(debug)
            cout << "[DTEfficiencyTask] Theta segments has: " << zRecHits.size()
                 << " hits, skipping" << endl; // FIXME: info output
          continue;
        }
        copy(zRecHits.begin(), zRecHits.end(), back_inserter(recHits1D));
      }

      // Skip the segment if it has more than 1 hit on the same layer
      vector<DTWireId> wireMap; 
      for(vector<DTRecHit1D>::const_iterator recHit1D = recHits1D.begin();
          recHit1D != recHits1D.end();
          recHit1D++) {
        wireMap.push_back((*recHit1D).wireId());
      }

      bool hitsOnSameLayer = false;
      for(vector<DTWireId>::const_iterator channelId = wireMap.begin();
          channelId != wireMap.end(); channelId++) {
        vector<DTWireId>::const_iterator next = channelId;
        next++;
        for(vector<DTWireId>::const_iterator ite = next; ite != wireMap.end(); ite++) {
          if((*channelId).layerId() == (*ite).layerId()) {
            hitsOnSameLayer = true;
          }
        }
      }
      if(hitsOnSameLayer) {
        if(debug)
          cout << "[DTEfficiencyTask] This RecHit has 2 hits on the same layer, skipping!" << endl;
        continue;
      }


      // Select 2D segments with angle smaller than 45 deg
      LocalVector phiDirectionInChamber = (*phiSeg).localDirection();
      if(rPhi && fabs(phiDirectionInChamber.x()/phiDirectionInChamber.z()) > 1) {
        if(debug) {
          cout << "         RPhi segment has angle > 45 deg, skipping! " << endl;
          cout << "              Theta = " << phiDirectionInChamber.theta() << endl;
        }
        continue;
      }
      if(rZ) {
         LocalVector zDirectionInChamber = (*zSeg).localDirection();
         if(fabs(zDirectionInChamber.y()/zDirectionInChamber.z()) > 1) {
           if(debug){
             cout << "         RZ segment has angle > 45 deg, skipping! "  << endl;
             cout << "              Theta = " << zDirectionInChamber.theta() << endl;
           }
           continue;
         }
      }


      // Skip if the 4D segment has only 10 hits
      if(recHits1D.size() == 10) {
        if(debug)
          cout << "[DTEfficiencyTask] 4D Segment with only 10 hits, skipping!" << endl;
        continue;
      }


      // Analyse the case of 11 recHits for MB1,MB2,MB3 and of 7 recHits for MB4
      if((rPhi && recHits1D.size() == 7) || (rZ && recHits1D.size() == 11)) {

        if(debug) {
          if(rPhi && recHits1D.size() == 7)
            cout << "[DTEfficiencyTask] MB4 Segment with only 7 hits!" << endl;
          if(rZ && recHits1D.size() == 11)
             cout << "[DTEfficiencyTask] 4D Segment with only 11 hits!" << endl;
        }

        // Find the layer without RecHits ----------------------------------------
        const vector<const DTSuperLayer*> SupLayers = chamber->superLayers();
        map<DTLayerId, bool> layerMap; 
        map<DTWireId, float> wireAndPosInChamberAtLayerZ;
        // Loop over layers and wires to fill a map
        for(vector<const DTSuperLayer*>::const_iterator superlayer = SupLayers.begin();
            superlayer != SupLayers.end();
            superlayer++) {
          const vector<const DTLayer*> Layers = (*superlayer)->layers();
          for(vector<const DTLayer*>::const_iterator layer = Layers.begin();
              layer != Layers.end();
              layer++) {
            layerMap.insert(make_pair((*layer)->id(), false));
            const int firstWire = dtGeom->layer((*layer)->id())->specificTopology().firstChannel();
            const int lastWire = dtGeom->layer((*layer)->id())->specificTopology().lastChannel();
            for(int i=firstWire; i <= lastWire; i++) {
              DTWireId wireId((*layer)->id(), i);            
              float wireX = (*layer)->specificTopology().wirePosition(wireId.wire());
              LocalPoint wirePosInLay(wireX,0,0);
              GlobalPoint wirePosGlob = (*layer)->toGlobal(wirePosInLay);
              LocalPoint wirePosInChamber = chamber->toLocal(wirePosGlob); 
              if((*superlayer)->id().superlayer() == 1 || (*superlayer)->id().superlayer() == 3) {
                wireAndPosInChamberAtLayerZ.insert(make_pair(wireId, wirePosInChamber.x()));
              } else {
                wireAndPosInChamberAtLayerZ.insert(make_pair(wireId, wirePosInChamber.y()));
              }
            }
          }
        }

        // Loop over segment 1D RecHit
        map<DTLayerId, int> NumWireMap; 
        for(vector<DTRecHit1D>::const_iterator recHit = recHits1D.begin();
            recHit != recHits1D.end();
            recHit++) {
          layerMap[(*recHit).wireId().layerId()]= true;
          NumWireMap[(*recHit).wireId().layerId()]= (*recHit).wireId().wire();
        }

        DTLayerId missLayerId;
        //Loop over the map and find the layer without hits
        for(map<DTLayerId, bool>::const_iterator iter = layerMap.begin();
            iter != layerMap.end(); iter++) {
          if(!(*iter).second) missLayerId = (*iter).first;
        }
        if(debug)
          cout << "[DTEfficiencyTask] Layer without recHits is: " << missLayerId << endl;
        // -------------------------------------------------------



        const DTLayer* missLayer = chamber->layer(missLayerId);
        
        LocalPoint missLayerPosInChamber = chamber->toLocal(missLayer->toGlobal(LocalPoint(0,0,0)));
        
        // Segment position at Wire z in chamber local frame
        LocalPoint segPosAtZLayer = (*segment4D).localPosition()
          + (*segment4D).localDirection()*missLayerPosInChamber.z()/cos((*segment4D).localDirection().theta());
        
        DTWireId missWireId;

        // Find the id of the cell without hit ---------------------------------------------------
        for(map<DTWireId, float>::const_iterator wireAndPos = wireAndPosInChamberAtLayerZ.begin();
            wireAndPos != wireAndPosInChamberAtLayerZ.end();
            wireAndPos++) {
          DTWireId wireId = (*wireAndPos).first;
          if(wireId.layerId() == missLayerId) {
            if(missLayerId.superlayerId().superlayer() == 1 || missLayerId.superlayerId().superlayer() == 3 ) {
              if(fabs(segPosAtZLayer.x() - (*wireAndPos).second) < 2.1)
                missWireId = wireId;  
            } else {
              if(fabs(segPosAtZLayer.y() - (*wireAndPos).second) < 2.1)
                missWireId = wireId;
            }
          }
        }
        if(debug)
          cout << "[DTEfficiencyTask] Cell without hit is: " << missWireId << endl;
        // ----------------------------------------------------------


        bool foundUnAssRechit = false;

        // Look for unassociated hits in this cell
        if(wireAnd1DRecHits.find(missWireId) != wireAnd1DRecHits.end()) {
          if(debug)
            cout << "[DTEfficiencyTask] Unassociated Hit found!" << endl;
          foundUnAssRechit = true;
        }


        for(map<DTLayerId, bool>::const_iterator iter = layerMap.begin();
            iter != layerMap.end(); iter++) {
          if((*iter).second) 
            fillHistos((*iter).first, dtGeom->layer((*iter).first)->specificTopology().firstChannel(), dtGeom->layer((*iter).first)->specificTopology().lastChannel(), NumWireMap[(*iter).first]);
          else
            fillHistos((*iter).first, dtGeom->layer((*iter).first)->specificTopology().firstChannel(), dtGeom->layer((*iter).first)->specificTopology().lastChannel(), missWireId.wire(), foundUnAssRechit);
        }

      } // End of the loop for segment with 7 or 11 recHits
      
      if((rPhi && recHits1D.size() == 8) || (rZ && recHits1D.size() == 12)) {
        map<DTLayerId, int> NumWireMap; 
        DTLayerId LayerID;
        for(vector<DTRecHit1D>::const_iterator recHit = recHits1D.begin();
            recHit != recHits1D.end();
            recHit++) {
          LayerID = (*recHit).wireId().layerId();
          NumWireMap[LayerID]= (*recHit).wireId().wire();
        }
        for(map<DTLayerId, int>::const_iterator iter = NumWireMap.begin();
            iter != NumWireMap.end(); iter++) {
          fillHistos((*iter).first, dtGeom->layer((*iter).first)->specificTopology().firstChannel(), dtGeom->layer((*iter).first)->specificTopology().lastChannel(), NumWireMap[(*iter).first]);
        }
      }

    } // End of loop over the 4D segments inside a sigle chamber
  } // End of loop over all tha chamber with at least a 4D segment in the event
}

void DTEfficiencyTask::beginJob

(

const edm::EventSetup &

c

)

[virtual]

BeginJob.

Reimplemented from edm::EDAnalyzer.

Definition at line 59 of file DTEfficiencyTask.cc.

References edm::ParameterSet::getParameter(), parameters, theRecHitLabel, and theRecHits4DLabel.

                                                           {
  // the name of the 4D rec hits collection
  theRecHits4DLabel = parameters.getParameter<string>("recHits4DLabel");
  // the name of the rechits collection
  theRecHitLabel = parameters.getParameter<string>("recHitLabel");
}

void DTEfficiencyTask::beginLuminosityBlock	(	edm::LuminosityBlock const &	lumiSeg,
		edm::EventSetup const &	context
	)			[virtual]

To reset the MEs.

Reimplemented from edm::EDAnalyzer.

Definition at line 67 of file DTEfficiencyTask.cc.

References edm::ParameterSet::getUntrackedParameter(), histo, histosPerL, i, edm::LuminosityBlock::id(), edm::LuminosityBlockID::luminosityBlock(), parameters, and size.

                                                                                                     {

  
  if(lumiSeg.id().luminosityBlock()%parameters.getUntrackedParameter<int>("ResetCycle", 3) == 0) {
    for(map<DTLayerId, vector<MonitorElement*> > ::const_iterator histo = histosPerL.begin();
        histo != histosPerL.end();
        histo++) {
      int size = (*histo).second.size();
      for(int i=0; i<size; i++){
        (*histo).second[i]->Reset();
      }
    }
  }
  
}

void DTEfficiencyTask::bookHistos	(	DTLayerId	lId,
		int	fisrtWire,
		int	lastWire
	)			[private]

Definition at line 379 of file DTEfficiencyTask.cc.

References DQMStore::book1D(), DTSuperLayerId::chamberId(), GenMuonPlsPt100GeV_cfg::cout, debug, lat::endl(), combine::histos, histosPerL, DTLayerId::layer(), DTChamberId::sector(), DQMStore::setCurrentFolder(), DTChamberId::station(), DTSuperLayerId::superlayer(), DTLayerId::superlayerId(), theDbe, muonGeometry::wheel, and DTChamberId::wheel().

Referenced by fillHistos().

                                                                            {
  if(debug)
    cout << "   Booking histos for L: " << lId << endl;

  // Compose the chamber name
  stringstream wheel; wheel << lId.superlayerId().chamberId().wheel();  
  stringstream station; station << lId.superlayerId().chamberId().station();    
  stringstream sector; sector << lId.superlayerId().chamberId().sector();       
  stringstream superLayer; superLayer << lId.superlayerId().superlayer();       
  stringstream layer; layer << lId.layer();

  string lHistoName =
    "_W" + wheel.str() +
    "_St" + station.str() +
    "_Sec" + sector.str() +
    "_SL" + superLayer.str()+
    "_L" + layer.str();
  
  theDbe->setCurrentFolder("DT/DTEfficiencyTask/Wheel" + wheel.str() +
                           "/Station" + station.str() +
                           "/Sector" + sector.str() +
                           "/SuperLayer" +superLayer.str());
  // Create the monitor elements
  vector<MonitorElement *> histos;
  // histo for hits associated to the 4D reconstructed segment
  histos.push_back(theDbe->book1D("hEffOccupancy"+lHistoName, "4D segments recHits occupancy",lastWire-firstWire+1, firstWire-0.5, lastWire+0.5));
  // histo for hits not associated to the segment
  histos.push_back(theDbe->book1D("hEffUnassOccupancy"+lHistoName, "4D segments recHits and Hits not associated occupancy",lastWire-firstWire+1, firstWire-0.5, lastWire+0.5));
  // histo for cells associated to the 4D reconstructed segment
  histos.push_back(theDbe->book1D("hRecSegmOccupancy"+lHistoName, "4D segments cells occupancy",lastWire-firstWire+1, firstWire-0.5, lastWire+0.5));
  
  histosPerL[lId] = histos;
}

void DTEfficiencyTask::endJob

(

void

)

[virtual]

Endjob.

Reimplemented from edm::EDAnalyzer.

Definition at line 84 of file DTEfficiencyTask.cc.

References DQMStore::rmdir(), and theDbe.

                             {
  theDbe->rmdir("DT/DTEfficiencyTask");
}

void DTEfficiencyTask::fillHistos	(	DTLayerId	lId,
		int	firstWire,
		int	lastWire,
		int	missingWire,
		bool	UnassHit
	)			[private]

Definition at line 428 of file DTEfficiencyTask.cc.

References bookHistos(), combine::histos, and histosPerL.

                                                 {
 if(histosPerL.find(lId) == histosPerL.end()){
      bookHistos(lId, firstWire, lastWire);
  }
 vector<MonitorElement *> histos =  histosPerL[lId];
 if(unassHit) 
   histos[1]->Fill(missingWire);
 histos[2]->Fill(missingWire);
}

void DTEfficiencyTask::fillHistos	(	DTLayerId	lId,
		int	firstWire,
		int	lastWire,
		int	numWire
	)			[private]

Definition at line 415 of file DTEfficiencyTask.cc.

References bookHistos(), combine::histos, and histosPerL.

Referenced by analyze().

                                               {
  if(histosPerL.find(lId) == histosPerL.end()){
      bookHistos(lId, firstWire, lastWire);
  }
  vector<MonitorElement *> histos =  histosPerL[lId]; 
  histos[0]->Fill(numWire);
  histos[1]->Fill(numWire);
  histos[2]->Fill(numWire);
}

---

Member Data Documentation

bool DTEfficiencyTask::debug [private]

Definition at line 58 of file DTEfficiencyTask.h.

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

std::map<DTLayerId, std::vector<MonitorElement*> > DTEfficiencyTask::histosPerL [private]

Definition at line 75 of file DTEfficiencyTask.h.

Referenced by beginLuminosityBlock(), bookHistos(), and fillHistos().

edm::ParameterSet DTEfficiencyTask::parameters [private]

Definition at line 66 of file DTEfficiencyTask.h.

Referenced by beginJob(), beginLuminosityBlock(), and DTEfficiencyTask().

DQMStore* DTEfficiencyTask::theDbe [private]

Definition at line 55 of file DTEfficiencyTask.h.

Referenced by bookHistos(), DTEfficiencyTask(), and endJob().

std::string DTEfficiencyTask::theRecHitLabel [private]

Definition at line 64 of file DTEfficiencyTask.h.

Referenced by analyze(), and beginJob().

std::string DTEfficiencyTask::theRecHits4DLabel [private]

Definition at line 61 of file DTEfficiencyTask.h.

Referenced by analyze(), and beginJob().

---

The documentation for this class was generated from the following files:

• DQM/DTMonitorModule/src/DTEfficiencyTask.h
• DQM/DTMonitorModule/src/DTEfficiencyTask.cc

---