CMS 3D CMS Logo
List of all members | Public Member Functions | Private Member Functions | Private Attributes
DTSegment2DQuality Class Reference

#include <DTSegment2DQuality.h>

Inheritance diagram for DTSegment2DQuality:
DQMGlobalEDAnalyzer< dtsegment2d::Histograms > edm::global::EDAnalyzer< edm::RunCache< dtsegment2d::Histograms >, Args... > edm::global::EDAnalyzerBase edm::EDConsumerBase

Public Member Functions

 DTSegment2DQuality (const edm::ParameterSet &pset)
 Constructor. More...
 
- Public Member Functions inherited from edm::global::EDAnalyzer< edm::RunCache< dtsegment2d::Histograms >, Args... >
 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 &, dtsegment2d::Histograms &) const override
 Book the DQM plots. More...
 
void dqmAnalyze (edm::Event const &, edm::EventSetup const &, dtsegment2d::Histograms const &) const override
 Perform the real analysis. More...
 

Private Attributes

bool debug_
 
edm::InputTag segment2DLabel_
 
edm::EDGetTokenT< DTRecSegment2DCollection > segment2DToken_
 
double sigmaResAngle_
 
double sigmaResPos_
 
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<typename ESProduct , Transition Tr = Transition::Event>
auto esConsumes (eventsetup::EventSetupRecordKey const &, ESInputTag const &tag)
 
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 2D DTSegments and plot resolution comparing reconstructed and simulated quantities

Author
S. Bolognesi and G. Cerminara - INFN Torino

Definition at line 33 of file DTSegment2DQuality.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 47 of file DTSegment2DQuality.cc.

References gather_cfg::cout, DTHitQualityUtils::debug, edm::ParameterSet::getParameter(), and edm::ParameterSet::getUntrackedParameter().

47  {
48  // get the debug parameter for verbose output
49  debug_ = pset.getUntrackedParameter<bool>("debug");
50  DTHitQualityUtils::debug = debug_;
51  // the name of the simhit collection
52  simHitLabel_ = pset.getUntrackedParameter<InputTag>("simHitLabel");
53  simHitToken_ = consumes<PSimHitContainer>(pset.getUntrackedParameter<InputTag>("simHitLabel"));
54  // the name of the 2D rec hit collection
55  segment2DLabel_ = pset.getUntrackedParameter<InputTag>("segment2DLabel");
56  segment2DToken_ = consumes<DTRecSegment2DCollection>(pset.getUntrackedParameter<InputTag>("segment2DLabel"));
57 
58  // sigma resolution on position
59  sigmaResPos_ = pset.getParameter<double>("sigmaResPos");
60  // sigma resolution on angle
61  sigmaResAngle_ = pset.getParameter<double>("sigmaResAngle");
62 
63  if (debug_) {
64  cout << "[DTSegment2DQuality] Constructor called " << endl;
65 }
66 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
DTSegment2DQuality::simHitToken_
edm::EDGetTokenT< edm::PSimHitContainer > simHitToken_
Definition: DTSegment2DQuality.h:49
DTHitQualityUtils::debug
std::atomic< bool > debug
Definition: DTHitQualityUtils.h:44
DTSegment2DQuality::segment2DToken_
edm::EDGetTokenT< DTRecSegment2DCollection > segment2DToken_
Definition: DTSegment2DQuality.h:50
DTSegment2DQuality::simHitLabel_
edm::InputTag simHitLabel_
Definition: DTSegment2DQuality.h:47
gather_cfg.cout
cout
Definition: gather_cfg.py:144
DTSegment2DQuality::segment2DLabel_
edm::InputTag segment2DLabel_
Definition: DTSegment2DQuality.h:48

Member Function Documentation

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

Book the DQM plots.

Implements DQMGlobalEDAnalyzer< dtsegment2d::Histograms >.

Definition at line 68 of file DTSegment2DQuality.cc.

References gather_cfg::cout, dtsegment2d::Histograms::h2DHitEff_RPhi, dtsegment2d::Histograms::h2DHitEff_RZ, dtsegment2d::Histograms::h2DHitEff_RZ_W0, dtsegment2d::Histograms::h2DHitEff_RZ_W1, dtsegment2d::Histograms::h2DHitEff_RZ_W2, dtsegment2d::Histograms::h2DHitRPhi, dtsegment2d::Histograms::h2DHitRZ, dtsegment2d::Histograms::h2DHitRZ_W0, dtsegment2d::Histograms::h2DHitRZ_W1, and dtsegment2d::Histograms::h2DHitRZ_W2.

68  {
69  histograms.h2DHitRPhi = std::make_unique<HRes2DHit> ("RPhi", booker, true, true);
70  histograms.h2DHitRZ = std::make_unique<HRes2DHit> ("RZ", booker, true, true);
71  histograms.h2DHitRZ_W0 = std::make_unique<HRes2DHit> ("RZ_W0", booker, true, true);
72  histograms.h2DHitRZ_W1 = std::make_unique<HRes2DHit> ("RZ_W1", booker, true, true);
73  histograms.h2DHitRZ_W2 = std::make_unique<HRes2DHit> ("RZ_W2", booker, true, true);
74 
75  histograms.h2DHitEff_RPhi = std::make_unique<HEff2DHit> ("RPhi", booker);
76  histograms.h2DHitEff_RZ = std::make_unique<HEff2DHit> ("RZ", booker);
77  histograms.h2DHitEff_RZ_W0 = std::make_unique<HEff2DHit> ("RZ_W0", booker);
78  histograms.h2DHitEff_RZ_W1 = std::make_unique<HEff2DHit> ("RZ_W1", booker);
79  histograms.h2DHitEff_RZ_W2 = std::make_unique<HEff2DHit> ("RZ_W2", booker);
80  if (debug_) {
81  cout << "[DTSegment2DQuality] hitsos created " << endl;
82  }
83 }
dtsegment2d::Histograms::h2DHitRZ_W2
std::unique_ptr< HRes2DHit > h2DHitRZ_W2
Definition: DTSegment2DQuality.cc:34
dtsegment2d::Histograms::h2DHitEff_RZ_W1
std::unique_ptr< HEff2DHit > h2DHitEff_RZ_W1
Definition: DTSegment2DQuality.cc:39
dtsegment2d::Histograms::h2DHitRZ_W0
std::unique_ptr< HRes2DHit > h2DHitRZ_W0
Definition: DTSegment2DQuality.cc:32
dtsegment2d::Histograms::h2DHitEff_RZ_W0
std::unique_ptr< HEff2DHit > h2DHitEff_RZ_W0
Definition: DTSegment2DQuality.cc:38
dtsegment2d::Histograms::h2DHitEff_RPhi
std::unique_ptr< HEff2DHit > h2DHitEff_RPhi
Definition: DTSegment2DQuality.cc:36
dtsegment2d::Histograms::h2DHitRPhi
std::unique_ptr< HRes2DHit > h2DHitRPhi
Definition: DTSegment2DQuality.cc:30
dtsegment2d::Histograms::h2DHitRZ_W1
std::unique_ptr< HRes2DHit > h2DHitRZ_W1
Definition: DTSegment2DQuality.cc:33
gather_cfg.cout
cout
Definition: gather_cfg.py:144
dtsegment2d::Histograms::h2DHitEff_RZ
std::unique_ptr< HEff2DHit > h2DHitEff_RZ
Definition: DTSegment2DQuality.cc:37
dtsegment2d::Histograms::h2DHitEff_RZ_W2
std::unique_ptr< HEff2DHit > h2DHitEff_RZ_W2
Definition: DTSegment2DQuality.cc:40
dtsegment2d::Histograms::h2DHitRZ
std::unique_ptr< HRes2DHit > h2DHitRZ
Definition: DTSegment2DQuality.cc:31
void DTSegment2DQuality::dqmAnalyze ( edm::Event const &  event,
edm::EventSetup const &  setup,
dtsegment2d::Histograms const &  histograms 
) const
overrideprivate

Perform the real analysis.

Definition at line 86 of file DTSegment2DQuality.cc.

References funct::abs(), gather_cfg::cout, DEFINE_FWK_MODULE, SoftLeptonByDistance_cfi::distance, PV3DBase< T, PVType, FrameType >::eta(), HRes2DHit::fill(), HEff2DHit::fill(), DTHitQualityUtils::findMuSimSegment(), DTHitQualityUtils::findMuSimSegmentDirAndPos(), DTHitQualityUtils::findSegmentAlphaAndBeta(), edm::EventSetup::get(), dtsegment2d::Histograms::h2DHitEff_RPhi, dtsegment2d::Histograms::h2DHitEff_RZ, dtsegment2d::Histograms::h2DHitEff_RZ_W0, dtsegment2d::Histograms::h2DHitEff_RZ_W1, dtsegment2d::Histograms::h2DHitEff_RZ_W2, dtsegment2d::Histograms::h2DHitRPhi, dtsegment2d::Histograms::h2DHitRZ, dtsegment2d::Histograms::h2DHitRZ_W0, dtsegment2d::Histograms::h2DHitRZ_W1, dtsegment2d::Histograms::h2DHitRZ_W2, edm::HandleBase::isValid(), DTRecSegment2D::localDirection(), DTRecSegment2D::localDirectionError(), DTRecSegment2D::localPosition(), DTRecSegment2D::localPositionError(), DTHitQualityUtils::mapMuSimHitsPerWire(), DTHitQualityUtils::mapSimHitsPerWire(), or, PV3DBase< T, PVType, FrameType >::phi(), rpcPointValidation_cfi::simHit, trackerHits::simHits, mathSSE::sqrt(), DTGeometry::superLayer(), GeomDet::toGlobal(), PV3DBase< T, PVType, FrameType >::x(), and LocalError::xx().

86  {
87  // Get the DT Geometry
88  ESHandle<DTGeometry> dtGeom;
89  setup.get<MuonGeometryRecord>().get(dtGeom);
90 
91  // Get the SimHit collection from the event
92  edm::Handle<PSimHitContainer> simHits;
93  event.getByToken(simHitToken_, simHits); // FIXME: second string to be removed
94 
95  // Map simHits by sl
96  map<DTSuperLayerId, PSimHitContainer > simHitsPerSl;
97  for (const auto & simHit : *simHits) {
98  // Create the id of the sl (the simHits in the DT known their wireId)
99  DTSuperLayerId slId = ((DTWireId(simHit.detUnitId())).layerId()).superlayerId();
100  // Fill the map
101  simHitsPerSl[slId].push_back(simHit);
102  }
103 
104  // Get the 2D rechits from the event
105  Handle<DTRecSegment2DCollection> segment2Ds;
106  event.getByToken(segment2DToken_, segment2Ds);
107 
108  if (not segment2Ds.isValid()) {
109  if (debug_) { cout << "[DTSegment2DQuality]**Warning: no 2DSegments with label: " << segment2DLabel_
110  << " in this event, skipping !" << endl;
111 }
112  return;
113  }
114 
115  // Loop over all superlayers containing a segment
116  DTRecSegment2DCollection::id_iterator slId;
117  for (slId = segment2Ds->id_begin();
118  slId != segment2Ds->id_end();
119  ++slId) {
120 
121  //------------------------- simHits ---------------------------//
122  // Get simHits of each superlayer
123  PSimHitContainer simHits = simHitsPerSl[(*slId)];
124 
125  // Map simhits per wire
126  map<DTWireId, PSimHitContainer> simHitsPerWire = DTHitQualityUtils::mapSimHitsPerWire(simHits);
127  map<DTWireId, const PSimHit*> muSimHitPerWire = DTHitQualityUtils::mapMuSimHitsPerWire(simHitsPerWire);
128  int nMuSimHit = muSimHitPerWire.size();
129  if (nMuSimHit == 0 or nMuSimHit == 1) {
130  if (debug_ and nMuSimHit == 1) {
131  cout << "[DTSegment2DQuality] Only " << nMuSimHit << " mu SimHit in this SL, skipping !" << endl;
132 }
133  continue; // If no or only one mu SimHit is found skip this SL
134  }
135  if (debug_) {
136  cout << "=== SL " << (*slId) << " has " << nMuSimHit << " SimHits" << endl;
137 }
138 
139  // Find outer and inner mu SimHit to build a segment
140  pair<const PSimHit*, const PSimHit*> inAndOutSimHit = DTHitQualityUtils::findMuSimSegment(muSimHitPerWire);
141  // Check that outermost and innermost SimHit are not the same
142  if (inAndOutSimHit.first == inAndOutSimHit.second ) {
143  cout << "[DTHitQualityUtils]***Warning: outermost and innermost SimHit are the same !" << endl;
144  continue;
145  }
146 
147  // Find direction and position of the sim Segment in SL RF
148  pair<LocalVector, LocalPoint> dirAndPosSimSegm = DTHitQualityUtils::findMuSimSegmentDirAndPos(inAndOutSimHit,
149  (*slId),&(*dtGeom));
150 
151  LocalVector simSegmLocalDir = dirAndPosSimSegm.first;
152  LocalPoint simSegmLocalPos = dirAndPosSimSegm.second;
153  if (debug_) {
154  cout << " Simulated segment: local direction " << simSegmLocalDir << endl
155  << " local position " << simSegmLocalPos << endl;
156 }
157  const DTSuperLayer* superLayer = dtGeom->superLayer(*slId);
158  GlobalPoint simSegmGlobalPos = superLayer->toGlobal(simSegmLocalPos);
159 
160  // Atan(x/z) angle and x position in SL RF
161  float angleSimSeg = DTHitQualityUtils::findSegmentAlphaAndBeta(simSegmLocalDir).first;
162  float posSimSeg = simSegmLocalPos.x();
163  // Position (in eta, phi coordinates) in the global RF
164  float etaSimSeg = simSegmGlobalPos.eta();
165  float phiSimSeg = simSegmGlobalPos.phi();
166 
167  //---------------------------- recHits --------------------------//
168  // Get the range of rechit for the corresponding slId
169  bool recHitFound = false;
170  DTRecSegment2DCollection::range range = segment2Ds->get(*slId);
171  int nsegm = distance(range.first, range.second);
172  if (debug_) {
173  cout << " Sl: " << *slId << " has " << nsegm
174  << " 2D segments" << endl;
175 }
176 
177  if (nsegm != 0) {
178  // Find the best RecHit: look for the 2D RecHit with the angle closest
179  // to that of segment made of SimHits.
180  // RecHits must have delta alpha and delta position within 5 sigma of
181  // the residual distribution (we are looking for residuals of segments
182  // usefull to the track fit) for efficency purpose
183  const DTRecSegment2D* bestRecHit = nullptr;
184  bool bestRecHitFound = false;
185  double deltaAlpha = 99999;
186 
187  // Loop over the recHits of this slId
188  for (DTRecSegment2DCollection::const_iterator segment2D = range.first;
189  segment2D != range.second;
190  ++segment2D) {
191  // Check the dimension
192  if ((*segment2D).dimension() != 2) {
193  if (debug_) { cout << "[DTSegment2DQuality]***Error: This is not 2D segment !!!" << endl;
194 }
195  abort();
196  }
197  // Segment Local Direction and position (in SL RF)
198  LocalVector recSegDirection = (*segment2D).localDirection();
199  LocalPoint recSegPosition = (*segment2D).localPosition();
200 
201  float recSegAlpha = DTHitQualityUtils::findSegmentAlphaAndBeta(recSegDirection).first;
202  if (debug_) {
203  cout << " RecSegment direction: " << recSegDirection << endl
204  << " position : " << recSegPosition << endl
205  << " alpha : " << recSegAlpha << endl;
206 }
207 
208  if (fabs(recSegAlpha - angleSimSeg) < deltaAlpha) {
209  deltaAlpha = fabs(recSegAlpha - angleSimSeg);
210  bestRecHit = &(*segment2D);
211  bestRecHitFound = true;
212  }
213  } // End of Loop over all 2D RecHits
214 
215  if (bestRecHitFound) {
216  // Best rechit direction and position in SL RF
217  LocalPoint bestRecHitLocalPos = bestRecHit->localPosition();
218  LocalVector bestRecHitLocalDir = bestRecHit->localDirection();
219 
220  LocalError bestRecHitLocalPosErr = bestRecHit->localPositionError();
221  LocalError bestRecHitLocalDirErr = bestRecHit->localDirectionError();
222 
223  float angleBestRHit = DTHitQualityUtils::findSegmentAlphaAndBeta(bestRecHitLocalDir).first;
224 
225  if (fabs(angleBestRHit - angleSimSeg) < 5*sigmaResAngle_ and
226  fabs(bestRecHitLocalPos.x() - posSimSeg) < 5*sigmaResPos_) {
227  recHitFound = true;
228  }
229 
230  // Fill Residual histos
231  HRes2DHit *hRes = nullptr;
232  if ((*slId).superlayer() == 1 or (*slId).superlayer() == 3) { // RPhi SL
233  hRes = histograms.h2DHitRPhi.get();
234  } else if ((*slId).superlayer() == 2) { // RZ SL
235  histograms.h2DHitRZ->fill(angleSimSeg,
236  angleBestRHit,
237  posSimSeg,
238  bestRecHitLocalPos.x(),
239  etaSimSeg,
240  phiSimSeg,
241  sqrt(bestRecHitLocalPosErr.xx()),
242  sqrt(bestRecHitLocalDirErr.xx())
243  );
244  if (abs((*slId).wheel()) == 0) {
245  hRes = histograms.h2DHitRZ_W0.get();
246  } else if (abs((*slId).wheel()) == 1) {
247  hRes = histograms.h2DHitRZ_W1.get();
248  } else if (abs((*slId).wheel()) == 2) {
249  hRes = histograms.h2DHitRZ_W2.get();
250 }
251  }
252  hRes->fill(angleSimSeg,
253  angleBestRHit,
254  posSimSeg,
255  bestRecHitLocalPos.x(),
256  etaSimSeg,
257  phiSimSeg,
258  sqrt(bestRecHitLocalPosErr.xx()),
259  sqrt(bestRecHitLocalDirErr.xx())
260  );
261  }
262  } // end of if (nsegm != 0)
263 
264  // Fill Efficiency plot
265  HEff2DHit *hEff = nullptr;
266  if ((*slId).superlayer() == 1 or (*slId).superlayer() == 3) { // RPhi SL
267  hEff = histograms.h2DHitEff_RPhi.get();
268  } else if ((*slId).superlayer() == 2) { // RZ SL
269  histograms.h2DHitEff_RZ->fill(etaSimSeg, phiSimSeg, posSimSeg, angleSimSeg, recHitFound);
270  if (abs((*slId).wheel()) == 0) {
271  hEff = histograms.h2DHitEff_RZ_W0.get();
272  } else if (abs((*slId).wheel()) == 1) {
273  hEff = histograms.h2DHitEff_RZ_W1.get();
274  } else if (abs((*slId).wheel()) == 2) {
275  hEff = histograms.h2DHitEff_RZ_W2.get();
276  }
277  }
278  hEff->fill(etaSimSeg, phiSimSeg, posSimSeg, angleSimSeg, recHitFound);
279  } // End of loop over superlayers
280 }
DTHitQualityUtils::findSegmentAlphaAndBeta
std::pair< double, double > findSegmentAlphaAndBeta(const LocalVector &direction)
DTRecSegment2D::localPositionError
LocalError localPositionError() const override
local position error in SL frame
Definition: DTRecSegment2D.cc:69
HRes2DHit::fill
void fill(float angleSimSegment, float angleRecSegment, float posSimSegment, float posRecSegment, float etaSimSegment, float phiSimSegment, float sigmaPos, float sigmaAngle)
Definition: Histograms.h:245
LocalError::xx
float xx() const
Definition: LocalError.h:24
rpcPointValidation_cfi.simHit
simHit
Definition: rpcPointValidation_cfi.py:24
edm::RangeMap< DTSuperLayerId, edm::OwnVector< DTSLRecSegment2D > >::range
std::pair< const_iterator, const_iterator > range
iterator range
Definition: RangeMap.h:50
DTSegment2DQuality::simHitToken_
edm::EDGetTokenT< edm::PSimHitContainer > simHitToken_
Definition: DTSegment2DQuality.h:49
DTHitQualityUtils::mapMuSimHitsPerWire
std::map< DTWireId, const PSimHit * > mapMuSimHitsPerWire(const std::map< DTWireId, edm::PSimHitContainer > &simHitWireMap)
Create a map between the Mu SimHits and corresponding MuBarWireId ;.
GeomDet::toGlobal
GlobalPoint toGlobal(const Local2DPoint &lp) const
Conversion to the global R.F. from the R.F. of the GeomDet.
Definition: GeomDet.h:54
HEff2DHit::fill
void fill(float etaSimSegm, float phiSimSegm, float posSimSegm, float angleSimSegm, bool fillRecHit)
Definition: Histograms.h:309
PV3DBase::phi
Geom::Phi< T > phi() const
Definition: PV3DBase.h:69
GeneralSetup.setup
def setup(process, global_tag, zero_tesla=False)
Definition: GeneralSetup.py:2
SoftLeptonByDistance_cfi.distance
distance
Definition: SoftLeptonByDistance_cfi.py:6
DTHitQualityUtils::mapSimHitsPerWire
std::map< DTWireId, edm::PSimHitContainer > mapSimHitsPerWire(const edm::PSimHitContainer &simhits)
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
or
The Signals That Services Can Subscribe To This is based on ActivityRegistry and is current per Services can connect to the signals distributed by the ActivityRegistry in order to monitor the activity of the application Each possible callback has some defined which we here list in angle e< void, edm::EventID const &, edm::Timestamp const & > We also list in braces which AR_WATCH_USING_METHOD_ is used for those or
Definition: Activities.doc:12
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
DTSegment2DQuality::segment2DToken_
edm::EDGetTokenT< DTRecSegment2DCollection > segment2DToken_
Definition: DTSegment2DQuality.h:50
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:74
DTRecSegment2D::localPosition
LocalPoint localPosition() const override
local position in SL frame
Definition: DTRecSegment2D.h:89
DTHitQualityUtils::findMuSimSegmentDirAndPos
std::pair< LocalVector, LocalPoint > findMuSimSegmentDirAndPos(const std::pair< const PSimHit *, const PSimHit * > &inAndOutSimHit, const DetId detId, const DTGeometry *muonGeom)
Find direction and position of a segment (in local RF) from outer and inner mu SimHit in the RF of ob...
PV3DBase::eta
T eta() const
Definition: PV3DBase.h:76
DTHitQualityUtils::findMuSimSegment
std::pair< const PSimHit *, const PSimHit * > findMuSimSegment(const std::map< DTWireId, const PSimHit * > &mapWireAndMuSimHit)
Find Innermost and outermost SimHit from Mu in a SL (they identify a simulated segment) ...
DTRecSegment2D::localDirection
LocalVector localDirection() const override
the local direction in SL frame
Definition: DTRecSegment2D.h:95
edm::OwnVector::const_iterator
Definition: OwnVector.h:41
edm::PSimHitContainer
std::vector< PSimHit > PSimHitContainer
Definition: PSimHitContainer.h:11
gather_cfg.cout
cout
Definition: gather_cfg.py:144
DTRecSegment2D::localDirectionError
LocalError localDirectionError() const override
the local direction error (xx,xy,yy) in SL frame: only xx is not 0.
Definition: DTRecSegment2D.cc:73
PV3DBase::x
T x() const
Definition: PV3DBase.h:62
DTSegment2DQuality::segment2DLabel_
edm::InputTag segment2DLabel_
Definition: DTSegment2DQuality.h:48
DTGeometry::superLayer
const DTSuperLayer * superLayer(const DTSuperLayerId &id) const
Return a DTSuperLayer given its id.
Definition: DTGeometry.cc:122

Member Data Documentation

bool DTSegment2DQuality::debug_
private

Definition at line 59 of file DTSegment2DQuality.h.

edm::InputTag DTSegment2DQuality::segment2DLabel_
private

Definition at line 48 of file DTSegment2DQuality.h.

edm::EDGetTokenT<DTRecSegment2DCollection> DTSegment2DQuality::segment2DToken_
private

Definition at line 50 of file DTSegment2DQuality.h.

double DTSegment2DQuality::sigmaResAngle_
private

Definition at line 56 of file DTSegment2DQuality.h.

double DTSegment2DQuality::sigmaResPos_
private

Definition at line 53 of file DTSegment2DQuality.h.

edm::InputTag DTSegment2DQuality::simHitLabel_
private

Definition at line 47 of file DTSegment2DQuality.h.

edm::EDGetTokenT<edm::PSimHitContainer> DTSegment2DQuality::simHitToken_
private

Definition at line 49 of file DTSegment2DQuality.h.

Generated for CMSSW Reference Manual by doxygen 1.8.11