DTRecHitQuality Class Reference

#include <DTRecHitQuality.h>

Inheritance diagram for DTRecHitQuality:

Public Member Functions
	DTRecHitQuality (const edm::ParameterSet &pset)
	Constructor. More...
Public Member Functions inherited from edm::global::EDAnalyzer< edm::RunCache< Histograms > >
	EDAnalyzer ()=default
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
bool	wantsStreamLuminosityBlocks () const final
bool	wantsStreamRuns () const final
Public Member Functions inherited from edm::global::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
	~EDAnalyzerBase () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	bookHistograms (DQMStore::ConcurrentBooker &, edm::Run const &, edm::EventSetup const &, Histograms &) const override
	Book the DQM plots. More...
template<typename type >
void	compute (const DTGeometry *dtGeom, const std::map< DTWireId, std::vector< PSimHit >> &simHitsPerWire, const std::map< DTWireId, std::vector< type >> &recHitsPerWire, Histograms const &histograms, int step) const
void	dqmAnalyze (edm::Event const &, edm::EventSetup const &, Histograms const &) const override
	Perform the real analysis. More...
template<typename type >
const type *	findBestRecHit (const DTLayer *layer, const DTWireId &wireId, const std::vector< type > &recHits, float simHitDist) const
std::map< DTWireId, std::vector< DTRecHit1DPair > >	map1DRecHitsPerWire (const DTRecHitCollection *dt1DRecHitPairs) const
std::map< DTWireId, std::vector< DTRecHit1D > >	map1DRecHitsPerWire (const DTRecSegment2DCollection *segment2Ds) const
std::map< DTWireId, std::vector< DTRecHit1D > >	map1DRecHitsPerWire (const DTRecSegment4DCollection *segment4Ds) const
float	recHitDistFromWire (const DTRecHit1DPair &hitPair, const DTLayer *layer) const
float	recHitDistFromWire (const DTRecHit1D &recHit, const DTLayer *layer) const
float	recHitPositionError (const DTRecHit1DPair &recHit) const
float	recHitPositionError (const DTRecHit1D &recHit) const
float	simHitDistFromFE (const DTLayer *layer, const DTWireId &wireId, const PSimHit &hit) const
float	simHitDistFromWire (const DTLayer *layer, const DTWireId &wireId, const PSimHit &hit) const
float	simHitImpactAngle (const DTLayer *layer, const DTWireId &wireId, const PSimHit &hit) const

Private Attributes
bool	debug_
bool	doall_
bool	doStep1_
bool	doStep2_
bool	doStep3_
bool	local_
edm::InputTag	recHitLabel_
edm::EDGetTokenT< DTRecHitCollection >	recHitToken_
edm::InputTag	segment2DLabel_
edm::EDGetTokenT< DTRecSegment2DCollection >	segment2DToken_
edm::InputTag	segment4DLabel_
edm::EDGetTokenT< DTRecSegment4DCollection >	segment4DToken_
edm::InputTag	simHitLabel_
edm::EDGetTokenT< edm::PSimHitContainer >	simHitToken_

Additional Inherited Members
Public Types inherited from edm::global::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::global::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
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

Basic analyzer class which accesses 1D DTRecHits and plot resolution comparing reconstructed and simulated quantities

Residual/pull plots are filled for the rechit with distance from wire closer to that of the muon simhit.

Efficiencies are defined as the fraction of muon simhits with a rechit in the same cell, for the given reconstruction step. Hence, for S2 and S3 the definition incorporate the segment reconstruction efficiency.

Author

G. Cerminara - INFN Torino

Definition at line 45 of file DTRecHitQuality.h.

Constructor & Destructor Documentation

DTRecHitQuality::DTRecHitQuality

(

const edm::ParameterSet &

pset

)

Constructor.

Definition at line 85 of file DTRecHitQuality.cc.

References edm::ParameterSet::getUntrackedParameter().

                                                         {
  // Get the debug parameter for verbose output
  debug_ = pset.getUntrackedParameter<bool>("debug");
  // the name of the simhit collection
  simHitLabel_ = pset.getUntrackedParameter<InputTag>("simHitLabel");
  simHitToken_ = consumes<PSimHitContainer>(pset.getUntrackedParameter<InputTag>("simHitLabel"));
  // the name of the 1D rec hit collection
  recHitLabel_ = pset.getUntrackedParameter<InputTag>("recHitLabel");
  recHitToken_ = consumes<DTRecHitCollection>(pset.getUntrackedParameter<InputTag>("recHitLabel"));
  // the name of the 2D rec hit collection
  segment2DLabel_ = pset.getUntrackedParameter<InputTag>("segment2DLabel");
  segment2DToken_ = consumes<DTRecSegment2DCollection>(pset.getUntrackedParameter<InputTag>("segment2DLabel"));
  // the name of the 4D rec hit collection
  segment4DLabel_ = pset.getUntrackedParameter<InputTag>("segment4DLabel");
  segment4DToken_ = consumes<DTRecSegment4DCollection>(pset.getUntrackedParameter<InputTag>("segment4DLabel"));
  // Switches for analysis at various steps
  doStep1_ = pset.getUntrackedParameter<bool>("doStep1", false);
  doStep2_ = pset.getUntrackedParameter<bool>("doStep2", false);
  doStep3_ = pset.getUntrackedParameter<bool>("doStep3", false);
  doall_ = pset.getUntrackedParameter<bool>("doall", false);
  local_ = pset.getUntrackedParameter<bool>("local", true);
}

Member Function Documentation

void DTRecHitQuality::bookHistograms ( DQMStore::ConcurrentBooker &  booker, edm::Run const &  run, edm::EventSetup const &  setup, Histograms &  histograms  ) const

overrideprivatevirtual

Book the DQM plots.

Implements DQMGlobalEDAnalyzer< Histograms >.

Definition at line 108 of file DTRecHitQuality.cc.

References Histograms::hEff_S1RPhi, Histograms::hEff_S1RPhiWS, Histograms::hEff_S1RZ, Histograms::hEff_S1RZ_W0, Histograms::hEff_S1RZ_W1, Histograms::hEff_S1RZ_W2, Histograms::hEff_S1RZWS, Histograms::hEff_S2RPhi, Histograms::hEff_S2RZ, Histograms::hEff_S2RZ_W0, Histograms::hEff_S2RZ_W1, Histograms::hEff_S2RZ_W2, Histograms::hEff_S3RPhi, Histograms::hEff_S3RPhiWS, Histograms::hEff_S3RZ, Histograms::hEff_S3RZ_W0, Histograms::hEff_S3RZ_W1, Histograms::hEff_S3RZ_W2, Histograms::hEff_S3RZWS, Histograms::hRes_S1RPhi, Histograms::hRes_S1RPhi_W0, Histograms::hRes_S1RPhi_W1, Histograms::hRes_S1RPhi_W2, Histograms::hRes_S1RZ, Histograms::hRes_S1RZ_W0, Histograms::hRes_S1RZ_W1, Histograms::hRes_S1RZ_W2, Histograms::hRes_S2RPhi, Histograms::hRes_S2RPhi_W0, Histograms::hRes_S2RPhi_W1, Histograms::hRes_S2RPhi_W2, Histograms::hRes_S2RZ, Histograms::hRes_S2RZ_W0, Histograms::hRes_S2RZ_W1, Histograms::hRes_S2RZ_W2, Histograms::hRes_S3RPhi, Histograms::hRes_S3RPhi_W0, Histograms::hRes_S3RPhi_W1, Histograms::hRes_S3RPhi_W2, Histograms::hRes_S3RPhiWS, Histograms::hRes_S3RZ, Histograms::hRes_S3RZ_W0, Histograms::hRes_S3RZ_W1, Histograms::hRes_S3RZ_W2, Histograms::hRes_S3RZWS, alignCSCRings::s, and w.

                                                                                                                                                    {
  if (doall_ && doStep1_) {
    histograms.hRes_S1RPhi    = new HRes1DHit("S1RPhi", booker, true, local_);      // RecHits, 1. step, RPhi
    histograms.hRes_S1RPhi_W0 = new HRes1DHit("S1RPhi_W0", booker, true, local_);   // RecHits, 1. step, RZ, wheel 0
    histograms.hRes_S1RPhi_W1 = new HRes1DHit("S1RPhi_W1", booker, true, local_);   // RecHits, 1. step, RZ, wheel +-1
    histograms.hRes_S1RPhi_W2 = new HRes1DHit("S1RPhi_W2", booker, true, local_);   // RecHits, 1. step, RZ, wheel +-2
    histograms.hRes_S1RZ      = new HRes1DHit("S1RZ", booker, true, local_);        // RecHits, 1. step, RZ
    histograms.hRes_S1RZ_W0   = new HRes1DHit("S1RZ_W0", booker, true, local_);     // RecHits, 1. step, RZ, wheel 0
    histograms.hRes_S1RZ_W1   = new HRes1DHit("S1RZ_W1", booker, true, local_);     // RecHits, 1. step, RZ, wheel +-1
    histograms.hRes_S1RZ_W2   = new HRes1DHit("S1RZ_W2", booker, true, local_);     // RecHits, 1. step, RZ, wheel +-2
    histograms.hEff_S1RPhi    = new HEff1DHit("S1RPhi", booker);                    // RecHits, 1. step, RPhi
    histograms.hEff_S1RZ      = new HEff1DHit("S1RZ", booker);                      // RecHits, 1. step, RZ
    histograms.hEff_S1RZ_W0   = new HEff1DHit("S1RZ_W0", booker);                   // RecHits, 1. step, RZ, wheel 0
    histograms.hEff_S1RZ_W1   = new HEff1DHit("S1RZ_W1", booker);                   // RecHits, 1. step, RZ, wheel +-1
    histograms.hEff_S1RZ_W2   = new HEff1DHit("S1RZ_W2", booker);                   // RecHits, 1. step, RZ, wheel +-2
  }
  if (doall_ && doStep2_) {
    histograms.hRes_S2RPhi    = new HRes1DHit("S2RPhi", booker, true, local_);      // RecHits, 2. step, RPhi
    histograms.hRes_S2RPhi_W0 = new HRes1DHit("S2RPhi_W0", booker, true, local_);   // RecHits, 2. step, RPhi, wheel 0
    histograms.hRes_S2RPhi_W1 = new HRes1DHit("S2RPhi_W1", booker, true, local_);   // RecHits, 2. step, RPhi, wheel +-1
    histograms.hRes_S2RPhi_W2 = new HRes1DHit("S2RPhi_W2", booker, true, local_);   // RecHits, 2. step, RPhi, wheel +-2
    histograms.hRes_S2RZ      = new HRes1DHit("S2RZ", booker, true, local_);        // RecHits, 2. step, RZ
    histograms.hRes_S2RZ_W0   = new HRes1DHit("S2RZ_W0", booker, true, local_);     // RecHits, 2. step, RZ, wheel 0
    histograms.hRes_S2RZ_W1   = new HRes1DHit("S2RZ_W1", booker, true, local_);     // RecHits, 2. step, RZ, wheel +-1
    histograms.hRes_S2RZ_W2   = new HRes1DHit("S2RZ_W2", booker, true, local_);     // RecHits, 2. step, RZ, wheel +-2
    histograms.hEff_S2RPhi    = new HEff1DHit("S2RPhi", booker);                    // RecHits, 2. step, RPhi
    histograms.hEff_S2RZ_W0   = new HEff1DHit("S2RZ_W0", booker);                   // RecHits, 2. step, RZ, wheel 0
    histograms.hEff_S2RZ_W1   = new HEff1DHit("S2RZ_W1", booker);                   // RecHits, 2. step, RZ, wheel +-1
    histograms.hEff_S2RZ_W2   = new HEff1DHit("S2RZ_W2", booker);                   // RecHits, 2. step, RZ, wheel +-2
    histograms.hEff_S2RZ      = new HEff1DHit("S2RZ", booker);                      // RecHits, 2. step, RZ
  }
  if (doStep3_) {
    histograms.hRes_S3RPhi    = new HRes1DHit("S3RPhi", booker, doall_, local_);    // RecHits, 3. step, RPhi
    histograms.hRes_S3RPhi_W0 = new HRes1DHit("S3RPhi_W0", booker, doall_, local_); // RecHits, 3. step, RPhi, wheel 0
    histograms.hRes_S3RPhi_W1 = new HRes1DHit("S3RPhi_W1", booker, doall_, local_); // RecHits, 3. step, RPhi, wheel +-1
    histograms.hRes_S3RPhi_W2 = new HRes1DHit("S3RPhi_W2", booker, doall_, local_); // RecHits, 3. step, RPhi, wheel +-2
    histograms.hRes_S3RZ      = new HRes1DHit("S3RZ", booker, doall_, local_);      // RecHits, 3. step, RZ
    histograms.hRes_S3RZ_W0   = new HRes1DHit("S3RZ_W0", booker, doall_, local_);   // RecHits, 3. step, RZ, wheel 0
    histograms.hRes_S3RZ_W1   = new HRes1DHit("S3RZ_W1", booker, doall_, local_);   // RecHits, 3. step, RZ, wheel +-1
    histograms.hRes_S3RZ_W2   = new HRes1DHit("S3RZ_W2", booker, doall_, local_);   // RecHits, 3. step, RZ, wheel +-2

    if (local_) {
      // Plots with finer granularity, not to be included in DQM
      TString name1 = "RPhi_W";
      TString name2 = "RZ_W";
      for (long w = 0; w <= 2; ++w) {
        for (long s = 1;s <= 4; ++s) {
          histograms.hRes_S3RPhiWS[w][s-1] = new HRes1DHit(("S3"+name1+w+"_St"+s).Data(), booker, doall_, local_);
          histograms.hEff_S1RPhiWS[w][s-1] = new HEff1DHit(("S1"+name1+w+"_St"+s).Data(), booker);
          histograms.hEff_S3RPhiWS[w][s-1] = new HEff1DHit(("S3"+name1+w+"_St"+s).Data(), booker);
          if (s != 4) {
            histograms.hRes_S3RZWS[w][s-1] = new HRes1DHit(("S3"+name2+w+"_St"+s).Data(), booker, doall_, local_);
            histograms.hEff_S1RZWS[w][s-1] = new HEff1DHit(("S1"+name2+w+"_St"+s).Data(), booker);
            histograms.hEff_S3RZWS[w][s-1] = new HEff1DHit(("S3"+name2+w+"_St"+s).Data(), booker);
          }
        }
      }
    }

    if (doall_) {
      histograms.hEff_S3RPhi  = new HEff1DHit("S3RPhi", booker);                    // RecHits, 3. step, RPhi
      histograms.hEff_S3RZ    = new HEff1DHit("S3RZ", booker);                      // RecHits, 3. step, RZ
      histograms.hEff_S3RZ_W0 = new HEff1DHit("S3RZ_W0", booker);                   // RecHits, 3. step, RZ, wheel 0
      histograms.hEff_S3RZ_W1 = new HEff1DHit("S3RZ_W1", booker);                   // RecHits, 3. step, RZ, wheel +-1
      histograms.hEff_S3RZ_W2 = new HEff1DHit("S3RZ_W2", booker);                   // RecHits, 3. step, RZ, wheel +-2
    }
  }
}

template<typename type >

void DTRecHitQuality::compute ( const DTGeometry *  dtGeom, const std::map< DTWireId, std::vector< PSimHit >> &  simHitsPerWire, const std::map< DTWireId, std::vector< type >> &  recHitsPerWire, Histograms const &  histograms, int  step  ) const

private

Definition at line 397 of file DTRecHitQuality.cc.

References funct::abs(), gather_cfg::cout, PV3DBase< T, PVType, FrameType >::eta(), HRes1DHit::fill(), HEff1DHit::fill(), DTHitQualityUtils::findMuSimHit(), Histograms::hEff_S1RPhi, Histograms::hEff_S1RPhiWS, Histograms::hEff_S1RZ, Histograms::hEff_S1RZ_W0, Histograms::hEff_S1RZ_W1, Histograms::hEff_S1RZ_W2, Histograms::hEff_S1RZWS, Histograms::hEff_S2RPhi, Histograms::hEff_S2RZ, Histograms::hEff_S2RZ_W0, Histograms::hEff_S2RZ_W1, Histograms::hEff_S2RZ_W2, Histograms::hEff_S3RPhi, Histograms::hEff_S3RPhiWS, Histograms::hEff_S3RZ, Histograms::hEff_S3RZ_W0, Histograms::hEff_S3RZ_W1, Histograms::hEff_S3RZ_W2, Histograms::hEff_S3RZWS, Histograms::hRes_S1RPhi, Histograms::hRes_S1RPhi_W0, Histograms::hRes_S1RPhi_W1, Histograms::hRes_S1RPhi_W2, Histograms::hRes_S1RZ, Histograms::hRes_S1RZ_W0, Histograms::hRes_S1RZ_W1, Histograms::hRes_S1RZ_W2, Histograms::hRes_S2RPhi, Histograms::hRes_S2RPhi_W0, Histograms::hRes_S2RPhi_W1, Histograms::hRes_S2RPhi_W2, Histograms::hRes_S2RZ, Histograms::hRes_S2RZ_W0, Histograms::hRes_S2RZ_W1, Histograms::hRes_S2RZ_W2, Histograms::hRes_S3RPhi, Histograms::hRes_S3RPhi_W0, Histograms::hRes_S3RPhi_W1, Histograms::hRes_S3RPhi_W2, Histograms::hRes_S3RPhiWS, Histograms::hRes_S3RZ, Histograms::hRes_S3RZ_W0, Histograms::hRes_S3RZ_W1, Histograms::hRes_S3RZ_W2, Histograms::hRes_S3RZWS, DTGeometry::layer(), PSimHit::localPosition(), PV3DBase< T, PVType, FrameType >::phi(), DTChamberId::station(), DTSuperLayerId::superLayer(), GeomDet::toGlobal(), DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

                                                                            {
  // Loop over cells with a muon SimHit
  for (const auto & wireAndSHits : simHitsPerWire) {
    DTWireId wireId = wireAndSHits.first;
    int wheel = wireId.wheel();
    int sl = wireId.superLayer();

    vector<PSimHit> simHitsInCell = wireAndSHits.second;

    // Get the layer
    const DTLayer* layer = dtGeom->layer(wireId);

    // Look for a mu hit in the cell
    const PSimHit* muSimHit = DTHitQualityUtils::findMuSimHit(simHitsInCell);
    if (muSimHit == nullptr) {
      if (debug_) {
        cout << "   No mu SimHit in channel: " << wireId << ", skipping! " << endl;
}
      continue; // Skip this cell
    }

    // Find the distance of the simhit from the wire
    float simHitWireDist = simHitDistFromWire(layer, wireId, *muSimHit);
    // Skip simhits out of the cell
    if (simHitWireDist>2.1) {
      if (debug_) {
        cout << "  [DTRecHitQuality]###Warning: The mu SimHit in out of the cell, skipping!" << endl;
}
      continue; // Skip this cell
    }
    GlobalPoint simHitGlobalPos = layer->toGlobal(muSimHit->localPosition());

    // find SH impact angle
    float simHitTheta = simHitImpactAngle(layer, wireId, *muSimHit);

    // find SH distance from FE
    float simHitFEDist = simHitDistFromFE(layer, wireId, *muSimHit);

    bool recHitReconstructed = false;

    // Look for RecHits in the same cell
    if (recHitsPerWire.find(wireId) == recHitsPerWire.end()) {
      // No RecHit found in this cell
      if (debug_) {
        cout << "   No RecHit found at Step: " << step << " in cell: " << wireId << endl;
}
    } else {
      recHitReconstructed = true;
      // vector<type> recHits = (*wireAndRecHits).second;
      const vector<type>& recHits = recHitsPerWire.at(wireId);
      if (debug_) {
        cout << "   " << recHits.size() << " RecHits, Step " << step << " in channel: " << wireId << endl;
}

      // Find the best RecHit
      const type* theBestRecHit = findBestRecHit(layer, wireId, recHits, simHitWireDist);


      float recHitWireDist =  recHitDistFromWire(*theBestRecHit, layer);
      if (debug_) {
        cout << "    SimHit distance from wire: " << simHitWireDist << endl
             << "    SimHit distance from FE:   " << simHitFEDist << endl
             << "    SimHit angle in layer RF:  " << simHitTheta << endl
             << "    RecHit distance from wire: " << recHitWireDist << endl;
}
      float recHitErr = recHitPositionError(*theBestRecHit);
      HRes1DHit *hRes = nullptr;
      HRes1DHit *hResTot = nullptr;

      // Mirror angle in phi so that + and - wheels can be plotted together
      if (mirrorMinusWheels && wheel<0 && sl!= 2) {
        simHitTheta *= -1.;
        // Note: local X, if used, would have to be mirrored as well
      }

      // Fill residuals and pulls
      // Select the histo to be filled
      if (step == 1) {
        // Step 1
        if (sl != 2) {
          hResTot = histograms.hRes_S1RPhi;
          if (wheel == 0) {
            hRes = histograms.hRes_S1RPhi_W0;
}
          if (abs(wheel) == 1) {
            hRes = histograms.hRes_S1RPhi_W1;
}
          if (abs(wheel) == 2) {
            hRes = histograms.hRes_S1RPhi_W2;
}
        } else {
          hResTot = histograms.hRes_S1RZ;
          if (wheel == 0) {
            hRes = histograms.hRes_S1RZ_W0;
}
          if (abs(wheel) == 1) {
            hRes = histograms.hRes_S1RZ_W1;
}
          if (abs(wheel) == 2) {
            hRes = histograms.hRes_S1RZ_W2;
}
        }

      } else if (step == 2) {
        // Step 2
        if (sl != 2) {
          hRes = histograms.hRes_S2RPhi;
          if (wheel == 0) {
            hRes = histograms.hRes_S2RPhi_W0;
}
          if (abs(wheel) == 1) {
            hRes = histograms.hRes_S2RPhi_W1;
}
          if (abs(wheel) == 2) {
            hRes = histograms.hRes_S2RPhi_W2;
}
        } else {
          hResTot = histograms.hRes_S2RZ;
          if (wheel == 0) {
            hRes = histograms.hRes_S2RZ_W0;
}
          if (abs(wheel) == 1) {
            hRes = histograms.hRes_S2RZ_W1;
}
          if (abs(wheel) == 2) {
            hRes = histograms.hRes_S2RZ_W2;
}
        }

      } else if (step == 3) {
        // Step 3
        if (sl != 2) {
          hResTot = histograms.hRes_S3RPhi;
          if (wheel == 0) {
            hRes = histograms.hRes_S3RPhi_W0;
}
          if (abs(wheel) == 1) {
            hRes = histograms.hRes_S3RPhi_W1;
}
          if (abs(wheel) == 2) {
            hRes = histograms.hRes_S3RPhi_W2;
}
          if (local_) {
            histograms.hRes_S3RPhiWS[abs(wheel)][wireId.station()-1]->fill(simHitWireDist, simHitTheta, simHitFEDist, recHitWireDist, simHitGlobalPos.eta(), simHitGlobalPos.phi(), recHitErr, wireId.station());
}
        } else {
          hResTot = histograms.hRes_S3RZ;
          if (wheel == 0) {
            hRes = histograms.hRes_S3RZ_W0;
}
          if (abs(wheel) == 1) {
            hRes = histograms.hRes_S3RZ_W1;
}
          if (abs(wheel) == 2) {
            hRes = histograms.hRes_S3RZ_W2;
}
          if (local_) {
            histograms.hRes_S3RZWS[abs(wheel)][wireId.station()-1]->fill(simHitWireDist, simHitTheta, simHitFEDist, recHitWireDist, simHitGlobalPos.eta(), simHitGlobalPos.phi(), recHitErr, wireId.station());
}
        }
      }

      // Fill
      hRes->fill(simHitWireDist, simHitTheta, simHitFEDist, recHitWireDist, simHitGlobalPos.eta(),
                 simHitGlobalPos.phi(), recHitErr, wireId.station());
      if (hResTot != nullptr) {
        hResTot->fill(simHitWireDist, simHitTheta, simHitFEDist, recHitWireDist, simHitGlobalPos.eta(),
                      simHitGlobalPos.phi(), recHitErr, wireId.station());
}
    }

    // Fill Efficiencies
    if (doall_) {
      HEff1DHit *hEff = nullptr;
      HEff1DHit *hEffTot = nullptr;
      if (step == 1) {
        // Step 1
        if (sl != 2) {
          hEff = histograms.hEff_S1RPhi;
          if (local_) { histograms.hEff_S1RPhiWS[abs(wheel)][wireId.station()-1]->fill(simHitWireDist, simHitGlobalPos.eta(), simHitGlobalPos.phi(), recHitReconstructed);
}
        } else {
          hEffTot = histograms.hEff_S1RZ;
          if (wheel == 0) {
            hEff = histograms.hEff_S1RZ_W0;
}
          if (abs(wheel) == 1) {
            hEff = histograms.hEff_S1RZ_W1;
}
          if (abs(wheel) == 2) {
            hEff = histograms.hEff_S1RZ_W2;
}
          if (local_) { histograms.hEff_S1RZWS[abs(wheel)][wireId.station()-1]->fill(simHitWireDist, simHitGlobalPos.eta(), simHitGlobalPos.phi(), recHitReconstructed);
}
        }

      } else if (step == 2) {
        // Step 2
        if (sl != 2) {
          hEff = histograms.hEff_S2RPhi;
        } else {
          hEffTot = histograms.hEff_S2RZ;
          if (wheel == 0) {
            hEff = histograms.hEff_S2RZ_W0;
}
          if (abs(wheel) == 1) {
            hEff = histograms.hEff_S2RZ_W1;
}
          if (abs(wheel) == 2) {
            hEff = histograms.hEff_S2RZ_W2;
}
        }

      } else if (step == 3) {
        // Step 3
        if (sl != 2) {
          hEff = histograms.hEff_S3RPhi;
          if (local_) { histograms.hEff_S3RPhiWS[abs(wheel)][wireId.station()-1]->fill(simHitWireDist, simHitGlobalPos.eta(), simHitGlobalPos.phi(), recHitReconstructed);
}
        } else {
          hEffTot = histograms.hEff_S3RZ;
          if (wheel == 0) {
            hEff = histograms.hEff_S3RZ_W0;
}
          if (abs(wheel) == 1) {
            hEff = histograms.hEff_S3RZ_W1;
}
          if (abs(wheel) == 2) {
            hEff = histograms.hEff_S3RZ_W2;
}
          if (local_) { histograms.hEff_S3RZWS[abs(wheel)][wireId.station()-1]->fill(simHitWireDist, simHitGlobalPos.eta(), simHitGlobalPos.phi(), recHitReconstructed);
}
        }

      }
      // Fill
      hEff->fill(simHitWireDist, simHitGlobalPos.eta(), simHitGlobalPos.phi(), recHitReconstructed);
      if (hEffTot != nullptr) {
        hEffTot->fill(simHitWireDist, simHitGlobalPos.eta(), simHitGlobalPos.phi(), recHitReconstructed);
}
    }
  }
}

void DTRecHitQuality::dqmAnalyze ( edm::Event const &  event, edm::EventSetup const &  setup, Histograms const &  histograms  ) const

overrideprivate

Perform the real analysis.

Definition at line 178 of file DTRecHitQuality.cc.

References bookConverter::compute(), gather_cfg::cout, edm::EventSetup::get(), ecalSelectiveReadoutValidation_cfi::histograms, edm::HandleBase::isValid(), DTHitQualityUtils::mapSimHitsPerWire(), edm::Handle< T >::product(), edm::ESHandle< T >::product(), and trackerHits::simHits.

                                                                                                                      {
  if (debug_) {
    cout << "--- [DTRecHitQuality] Analysing Event: #Run: " << event.id().run()
         << " #Event: " << event.id().event() << endl;
}

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

  // Get the SimHit collection from the event
  Handle<PSimHitContainer> simHits;
  event.getByToken(simHitToken_, simHits);

  // Map simhits per wire
  map<DTWireId, PSimHitContainer > simHitsPerWire =
    DTHitQualityUtils::mapSimHitsPerWire(*(simHits.product()));

  //=======================================================================================
  // RecHit analysis at Step 1
  if (doStep1_ && doall_) {
    if (debug_) {
      cout << "  -- DTRecHit S1: begin analysis:" << endl;
}
    // Get the rechit collection from the event
    Handle<DTRecHitCollection> dtRecHits;
    event.getByToken(recHitToken_, dtRecHits);

    if (!dtRecHits.isValid()) {
      if (debug_) { cout << "[DTRecHitQuality]**Warning: no 1DRechits with label: " << recHitLabel_ << " in this event, skipping!" << endl;
}
      return;
    }

    // Map rechits per wire
    auto const& recHitsPerWire = map1DRecHitsPerWire(dtRecHits.product());
    compute(dtGeom.product(), simHitsPerWire, recHitsPerWire, histograms, 1);
  }

  //=======================================================================================
  // RecHit analysis at Step 2
  if (doStep2_ && doall_) {
    if (debug_) {
      cout << "  -- DTRecHit S2: begin analysis:" << endl;
}

    // Get the 2D rechits from the event
    Handle<DTRecSegment2DCollection> segment2Ds;
    event.getByToken(segment2DToken_, segment2Ds);

    if (!segment2Ds.isValid()) {
      if (debug_) {
        cout << "[DTRecHitQuality]**Warning: no 2DSegments with label: " << segment2DLabel_
             << " in this event, skipping!" << endl;
}

    }
    else{
      // Map rechits per wire
      auto const& recHitsPerWire = map1DRecHitsPerWire(segment2Ds.product());
      compute(dtGeom.product(), simHitsPerWire, recHitsPerWire, histograms, 2);
    }
  }

  //=======================================================================================
  // RecHit analysis at Step 3
  if (doStep3_) {
    if (debug_) {
      cout << "  -- DTRecHit S3: begin analysis:" << endl;
}

    // Get the 4D rechits from the event
    Handle<DTRecSegment4DCollection> segment4Ds;
    event.getByToken(segment4DToken_, segment4Ds);

    if (!segment4Ds.isValid()) {
      if (debug_) { cout << "[DTRecHitQuality]**Warning: no 4D Segments with label: " << segment4DLabel_
        << " in this event, skipping!" << endl;
}
      return;
    }

    // Map rechits per wire
    auto const& recHitsPerWire = map1DRecHitsPerWire(segment4Ds.product());
    compute(dtGeom.product(), simHitsPerWire, recHitsPerWire, histograms, 3);
  }
}

template<typename type >

const type * DTRecHitQuality::findBestRecHit ( const DTLayer *  layer, const DTWireId &  wireId, const std::vector< type > &  recHits, float  simHitDist  ) const

private

Definition at line 358 of file DTRecHitQuality.cc.

References rpcPointValidation_cfi::recHit.

                                                              {
  float res = 99999;
  const type* theBestRecHit = nullptr;
  // Loop over RecHits within the cell
  for (auto recHit = recHits.begin();
      recHit != recHits.end();
      recHit++) {
    float distTmp = recHitDistFromWire(*recHit, layer);
    if (fabs(distTmp-simHitDist) < res) {
      res = fabs(distTmp-simHitDist);
      theBestRecHit = &(*recHit);
    }
  } // End of loop over RecHits within the cell

  return theBestRecHit;
}

map< DTWireId, vector< DTRecHit1DPair > > DTRecHitQuality::map1DRecHitsPerWire ( const DTRecHitCollection *  dt1DRecHitPairs ) const

private

Definition at line 268 of file DTRecHitQuality.cc.

                                                                                    {
  map<DTWireId, vector<DTRecHit1DPair>> ret;

  for (const auto & dt1DRecHitPair : *dt1DRecHitPairs) {
    ret[dt1DRecHitPair.wireId()].push_back(dt1DRecHitPair);
  }

  return ret;
}

map< DTWireId, vector< DTRecHit1D > > DTRecHitQuality::map1DRecHitsPerWire ( const DTRecSegment2DCollection *  segment2Ds ) const

private

Definition at line 281 of file DTRecHitQuality.cc.

                                                                                     {
  map<DTWireId, vector<DTRecHit1D>> ret;

  // Loop over all 2D segments
  for (const auto & segment2D : *segment2Ds) {
    vector<DTRecHit1D> component1DHits = segment2D.specificRecHits();
    // Loop over all component 1D hits
    for (auto & component1DHit : component1DHits) {
      ret[component1DHit.wireId()].push_back(component1DHit);
    }
  }
  return ret;
}

map< DTWireId, std::vector< DTRecHit1D > > DTRecHitQuality::map1DRecHitsPerWire ( const DTRecSegment4DCollection *  segment4Ds ) const

private

Definition at line 297 of file DTRecHitQuality.cc.

References hfClusterShapes_cfi::hits.

                                                                                     {
  map<DTWireId, vector<DTRecHit1D>> ret;
  // Loop over all 4D segments
  for (const auto & segment4D : *segment4Ds) {
    // Get component 2D segments
    vector<const TrackingRecHit*> segment2Ds = segment4D.recHits();
    // Loop over 2D segments:
    for (auto & segment2D : segment2Ds) {
      // Get 1D component rechits
      vector<const TrackingRecHit*> hits = segment2D->recHits();
      // Loop over them
      for (auto & hit : hits) {
        const auto* hit1D = dynamic_cast<const DTRecHit1D*>(hit);
        ret[hit1D->wireId()].push_back(*hit1D);
      }
    }
  }

  return ret;
}

float DTRecHitQuality::recHitDistFromWire ( const DTRecHit1DPair &  hitPair, const DTLayer *  layer  ) const

private

Definition at line 381 of file DTRecHitQuality.cc.

References DTEnums::Left, DTRecHit1DPair::localPosition(), DTEnums::Right, and PV3DBase< T, PVType, FrameType >::x().

                                                                                             {
  // Compute the rechit distance from wire
  return fabs(hitPair.localPosition(DTEnums::Left).x() -
              hitPair.localPosition(DTEnums::Right).x())/2.;
}

float DTRecHitQuality::recHitDistFromWire ( const DTRecHit1D &  recHit, const DTLayer *  layer  ) const

private

Definition at line 391 of file DTRecHitQuality.cc.

References DTRecHit1D::localPosition(), DTLayer::specificTopology(), DTWireId::wire(), DTRecHit1D::wireId(), DTTopology::wirePosition(), and PV3DBase< T, PVType, FrameType >::x().

                                                                                        {
  return fabs(recHit.localPosition().x() - layer->specificTopology().wirePosition(recHit.wireId().wire()));
}

float DTRecHitQuality::recHitPositionError ( const DTRecHit1DPair &  recHit ) const

private

Definition at line 645 of file DTRecHitQuality.cc.

References DTEnums::Left, DTRecHit1DPair::localPositionError(), mathSSE::sqrt(), and LocalError::xx().

                                                                             {
  return sqrt(recHit.localPositionError(DTEnums::Left).xx());
}

float DTRecHitQuality::recHitPositionError ( const DTRecHit1D &  recHit ) const

private

Definition at line 650 of file DTRecHitQuality.cc.

References DEFINE_FWK_MODULE, DTRecHit1D::localPositionError(), mathSSE::sqrt(), and LocalError::xx().

                                                                         {
  return sqrt(recHit.localPositionError().xx());
}

float DTRecHitQuality::simHitDistFromFE ( const DTLayer *  layer, const DTWireId &  wireId, const PSimHit &  hit  ) const

private

Definition at line 342 of file DTRecHitQuality.cc.

References DTTopology::cellLenght(), PSimHit::entryPoint(), PSimHit::exitPoint(), DTLayer::specificTopology(), and PV3DBase< T, PVType, FrameType >::y().

                                                                  {
  LocalPoint entryP = hit.entryPoint();
  LocalPoint exitP = hit.exitPoint();
  float wireLenght = layer->specificTopology().cellLenght();
  // FIXME: should take only wireLenght/2.;
  // moreover, pos+cellLenght/2. is shorter than the distance from FE.
  // In fact it would make more sense to make plots vs y.
  return (entryP.y()+exitP.y())/2.+wireLenght;
}

float DTRecHitQuality::simHitDistFromWire ( const DTLayer *  layer, const DTWireId &  wireId, const PSimHit &  hit  ) const

private

Definition at line 319 of file DTRecHitQuality.cc.

References PSimHit::entryPoint(), PSimHit::exitPoint(), DTLayer::specificTopology(), DTWireId::wire(), DTTopology::wirePosition(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::z().

                                                                    {
  float xwire = layer->specificTopology().wirePosition(wireId.wire());
  LocalPoint entryP = hit.entryPoint();
  LocalPoint exitP = hit.exitPoint();
  float xEntry = entryP.x() - xwire;
  float xExit  = exitP.x() - xwire;

  return fabs(xEntry - (entryP.z()*(xExit-xEntry))/(exitP.z()-entryP.z())); // FIXME: check...
}

float DTRecHitQuality::simHitImpactAngle ( const DTLayer *  layer, const DTWireId &  wireId, const PSimHit &  hit  ) const

private

Definition at line 332 of file DTRecHitQuality.cc.

References PSimHit::entryPoint(), PSimHit::exitPoint(), theta(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::z().

                                                                   {
  LocalPoint entryP = hit.entryPoint();
  LocalPoint exitP = hit.exitPoint();
  float theta =(exitP.x()-entryP.x())/(exitP.z()-entryP.z());
  return atan(theta);
}

Member Data Documentation

bool DTRecHitQuality::debug_

private

Definition at line 60 of file DTRecHitQuality.h.

bool DTRecHitQuality::doall_

private

Definition at line 78 of file DTRecHitQuality.h.

bool DTRecHitQuality::doStep1_

private

Definition at line 74 of file DTRecHitQuality.h.

bool DTRecHitQuality::doStep2_

private

Definition at line 75 of file DTRecHitQuality.h.

bool DTRecHitQuality::doStep3_

private

Definition at line 76 of file DTRecHitQuality.h.

bool DTRecHitQuality::local_

private

Definition at line 77 of file DTRecHitQuality.h.

edm::InputTag DTRecHitQuality::recHitLabel_

private

Definition at line 69 of file DTRecHitQuality.h.

edm::EDGetTokenT<DTRecHitCollection> DTRecHitQuality::recHitToken_

private

Definition at line 64 of file DTRecHitQuality.h.

edm::InputTag DTRecHitQuality::segment2DLabel_

private

Definition at line 70 of file DTRecHitQuality.h.

edm::EDGetTokenT<DTRecSegment2DCollection> DTRecHitQuality::segment2DToken_

private

Definition at line 65 of file DTRecHitQuality.h.

edm::InputTag DTRecHitQuality::segment4DLabel_

private

Definition at line 71 of file DTRecHitQuality.h.

edm::EDGetTokenT<DTRecSegment4DCollection> DTRecHitQuality::segment4DToken_

private

Definition at line 66 of file DTRecHitQuality.h.

edm::InputTag DTRecHitQuality::simHitLabel_

private

Definition at line 66 of file DTRecHitQuality.h.

edm::EDGetTokenT<edm::PSimHitContainer> DTRecHitQuality::simHitToken_

private

Definition at line 63 of file DTRecHitQuality.h.