DTEfficiencyTask Class Reference

#include <DTEfficiencyTask.h>

Inheritance diagram for DTEfficiencyTask:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &setup)
void	beginJob ()
	BeginJob. More...
void	beginLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &context)
	To reset the MEs. More...
	DTEfficiencyTask (const edm::ParameterSet &pset)
	Constructor. More...
void	endJob ()
	Endjob. More...
virtual	~DTEfficiencyTask ()
	Destructor. More...
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndex	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Member Functions
void	bookHistos (DTLayerId lId, int fisrtWire, int lastWire)
void	fillHistos (DTLayerId lId, int firstWire, int lastWire, int numWire)
void	fillHistos (DTLayerId lId, int firstWire, int lastWire, int missingWire, bool UnassHit)

Private Attributes
bool	debug
std::map< DTLayerId, std::vector< MonitorElement * > >	histosPerL
edm::ParameterSet	parameters
DQMStore *	theDbe
std::string	theRecHitLabel
std::string	theRecHits4DLabel

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
CurrentProcessingContext const *	currentContext () const
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

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

• single cell efficiency All histos are produced per Layer

Date:

2010/01/05 10:14:40

Revision:

1.7

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(), cppFunctionSkipper::operator, and Parameters::parameters.

                                                           {

  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 filterCSVwithJSON::copy, funct::cos(), gather_cfg::cout, debug, edm::EventSetup::get(), i, DTChamber::layer(), DTWireId::layerId(), DTRangeMapAccessor::layersBySuperLayer(), DTRecSegment2D::specificRecHits(), DTSuperLayerId::superlayer(), DTLayerId::superlayerId(), DTChamber::superLayers(), 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 = 0;
      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 <= 0; 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 ( void  )

virtual

BeginJob.

Reimplemented from edm::EDAnalyzer.

Definition at line 59 of file DTEfficiencyTask.cc.

References Parameters::parameters.

                               {
  // 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 timingPdfMaker::histo, i, edm::LuminosityBlockBase::id(), edm::LuminosityBlockID::luminosityBlock(), Association::map, Parameters::parameters, and findQualityFiles::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 DTSuperLayerId::chamberId(), gather_cfg::cout, debug, mergeVDriftHistosByStation::histos, DTLayerId::layer(), DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, DTSuperLayerId::superlayer(), DTLayerId::superlayerId(), and DTChamberId::wheel().

                                                                            {
  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.

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

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

private

Definition at line 415 of file DTEfficiencyTask.cc.

References bookHistos(), and mergeVDriftHistosByStation::histos.

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

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

private

Definition at line 428 of file DTEfficiencyTask.cc.

References bookHistos(), and mergeVDriftHistosByStation::histos.

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

Member Data Documentation

bool DTEfficiencyTask::debug

private

Definition at line 58 of file DTEfficiencyTask.h.

Referenced by pkg.AbstractPkg::generate(), pkg.AbstractPkg::get_kwds(), and pkg.AbstractPkg::write().

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

private

Definition at line 75 of file DTEfficiencyTask.h.

edm::ParameterSet DTEfficiencyTask::parameters

private

Definition at line 66 of file DTEfficiencyTask.h.

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

DQMStore* DTEfficiencyTask::theDbe

private

Definition at line 55 of file DTEfficiencyTask.h.

std::string DTEfficiencyTask::theRecHitLabel

private

Definition at line 64 of file DTEfficiencyTask.h.

std::string DTEfficiencyTask::theRecHits4DLabel

private

Definition at line 61 of file DTEfficiencyTask.h.