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

#include <SiPixelPhase1RecHitsV.h>

Inheritance diagram for SiPixelPhase1RecHitsV:
SiPixelPhase1Base one::DQMEDAnalyzer< T > HistogramManagerHolder one::dqmimplementation::DQMBaseClass< T... >

Public Member Functions

void analyze (const edm::Event &, const edm::EventSetup &) override
 
 SiPixelPhase1RecHitsV (const edm::ParameterSet &conf)
 
- Public Member Functions inherited from SiPixelPhase1Base
void analyze (edm::Event const &e, edm::EventSetup const &eSetup) override=0
 
void bookHistograms (DQMStore::IBooker &iBooker, edm::Run const &run, edm::EventSetup const &iSetup) override
 
 SiPixelPhase1Base (const edm::ParameterSet &iConfig)
 
 ~SiPixelPhase1Base () override
 
- Public Member Functions inherited from one::DQMEDAnalyzer< T >
 DQMEDAnalyzer ()=default
 
 DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
 
 DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
 
 ~DQMEDAnalyzer () override=default
 
- Public Member Functions inherited from HistogramManagerHolder
 HistogramManagerHolder (const edm::ParameterSet &iConfig)
 

Private Types

enum  {
  IN_TIME_BUNCH, OUT_TIME_BUNCH, NSIMHITS, RECHIT_X,
  RECHIT_Y, RES_X, RES_Y, ERROR_X,
  ERROR_Y, PULL_X, PULL_Y
}
 

Private Attributes

edm::EDGetTokenT< SiPixelRecHitCollection > srcToken_
 
TrackerHitAssociator::Config trackerHitAssociatorConfig_
 

Additional Inherited Members

- Protected Types inherited from SiPixelPhase1Base
enum  { DCS }
 
- Protected Member Functions inherited from SiPixelPhase1Base
bool checktrigger (const edm::Event &iEvent, const edm::EventSetup &iSetup, const unsigned trgidx) const
 
- Protected Attributes inherited from HistogramManagerHolder
GeometryInterface geometryInterface
 
std::vector< HistogramManager > histo
 

Detailed Description

Definition at line 16 of file SiPixelPhase1RecHitsV.h.

Member Enumeration Documentation

anonymous enum
private
Enumerator
IN_TIME_BUNCH 
OUT_TIME_BUNCH 
NSIMHITS 
RECHIT_X 
RECHIT_Y 
RES_X 
RES_Y 
ERROR_X 
ERROR_Y 
PULL_X 
PULL_Y 

Definition at line 17 of file SiPixelPhase1RecHitsV.h.

17  {
18  IN_TIME_BUNCH,
19  OUT_TIME_BUNCH,
20  NSIMHITS,
21  RECHIT_X,
22  RECHIT_Y,
23  RES_X,
24  RES_Y,
25  ERROR_X,
26  ERROR_Y,
27  PULL_X,
28  PULL_Y,
29  };

Constructor & Destructor Documentation

SiPixelPhase1RecHitsV::SiPixelPhase1RecHitsV ( const edm::ParameterSet &  conf)
explicit

Definition at line 19 of file SiPixelPhase1RecHitsV.cc.

19  :
20  SiPixelPhase1Base(iConfig),
21  trackerHitAssociatorConfig_( iConfig, consumesCollector() ),
22  srcToken_ ( consumes<SiPixelRecHitCollection>(iConfig.getParameter<edm::InputTag>("src")) )
23 {}
SiPixelPhase1Base::SiPixelPhase1Base
SiPixelPhase1Base(const edm::ParameterSet &iConfig)
Definition: SiPixelPhase1Base.cc:14
SiPixelPhase1RecHitsV::srcToken_
edm::EDGetTokenT< SiPixelRecHitCollection > srcToken_
Definition: SiPixelPhase1RecHitsV.h:37
SiPixelPhase1RecHitsV::trackerHitAssociatorConfig_
TrackerHitAssociator::Config trackerHitAssociatorConfig_
Definition: SiPixelPhase1RecHitsV.h:36

Member Function Documentation

void SiPixelPhase1RecHitsV::analyze ( const edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
)
override

Definition at line 25 of file SiPixelPhase1RecHitsV.cc.

References funct::abs(), TrackerHitAssociator::associateHit(), edmNew::DetSetVector< T >::begin(), EncodedEventId::bunchCrossing(), DEFINE_FWK_MODULE, edmNew::DetSetVector< T >::end(), PSimHit::entryPoint(), ERROR_X, ERROR_Y, EncodedEventId::event(), PSimHit::eventId(), PSimHit::exitPoint(), edm::Event::getByToken(), HistogramManagerHolder::histo, IN_TIME_BUNCH, input, edm::HandleBase::isValid(), funct::m, NSIMHITS, OUT_TIME_BUNCH, PULL_X, PULL_Y, RECHIT_X, RECHIT_Y, RES_X, RES_Y, rpcPointValidation_cfi::simHit, mathSSE::sqrt(), srcToken_, trackerHitAssociatorConfig_, PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), PV3DBase< T, PVType, FrameType >::y(), and LocalError::yy().

25  {
26  edm::Handle<SiPixelRecHitCollection> input;
27  iEvent.getByToken(srcToken_, input);
28  if (!input.isValid()) return;
29 
30  TrackerHitAssociator associate(iEvent, trackerHitAssociatorConfig_);
31 
32  SiPixelRecHitCollection::const_iterator it;
33  for (it = input->begin(); it != input->end(); ++it) {
34  auto id = DetId(it->detId());
35 
36  for(SiPixelRecHit const& rechit : *it) {
37  std::vector<PSimHit> associateSimHit;
38  associateSimHit = associate.associateHit(rechit);
39  std::vector<PSimHit>::const_iterator closestIt = associateSimHit.begin();
40 
41  LocalPoint lp = rechit.localPosition();
42  float rechit_x = lp.x();
43  float rechit_y = lp.y();
44 
45  LocalError lerr = rechit.localPositionError();
46  float lerr_x = sqrt(lerr.xx());
47  float lerr_y = sqrt(lerr.yy());
48 
49  // loop over associated sim hits and find the closest
50  if ( !associateSimHit.empty() ) {
51  float closestSimHit = 9999.9;
52 
53  for (std::vector<PSimHit>::const_iterator m = associateSimHit.begin(); m < associateSimHit.end(); m++) {
54  float sim_x1 ( (*m).entryPoint().x() ), sim_x2 ( (*m).exitPoint().x() ), sim_xpos ( 0.5*(sim_x1+sim_x2) );
55  float sim_y1 ( (*m).entryPoint().y() ), sim_y2 ( (*m).exitPoint().y() ), sim_ypos ( 0.5*(sim_y1+sim_y2) );
56 
57  float xres ( std::abs(sim_xpos - rechit_x) ), yres ( std::abs(sim_ypos - rechit_y) );
58  float dist = std::sqrt( xres*xres + yres*yres );
59 
60  if ( dist < closestSimHit ) {
61  closestSimHit = dist;
62  closestIt = m;
63  }
64  }
65  }
66 
67  // Sim Hit stuff
68 
69  if ( !associateSimHit.empty() ) {
70 
71  const PSimHit& simHit = *closestIt;
72 
73  int bunch = simHit.eventId().bunchCrossing();
74  int event = simHit.eventId().event();
75 
76  float sim_x1 ( simHit.entryPoint().x() ), sim_x2 ( simHit.exitPoint().x() ), sim_xpos ( 0.5*(sim_x1 + sim_x2) );
77  float sim_y1 ( simHit.entryPoint().y() ), sim_y2 ( simHit.exitPoint().y() ), sim_ypos ( 0.5*(sim_y1 + sim_y2) );
78 
79  float res_x = (rechit_x - sim_xpos) * 10000.0;
80  float res_y = (rechit_y - sim_ypos) * 10000.0;
81 
82  float pull_x = ( rechit_x - sim_xpos ) / lerr_x;
83  float pull_y = ( rechit_y - sim_ypos ) / lerr_y;
84 
85  // Now Plotting stuff
86  if ( bunch == 0 ) histo[IN_TIME_BUNCH].fill(bunch, id, &iEvent);
87  if ( bunch != 0 ) histo[OUT_TIME_BUNCH].fill(bunch, id, &iEvent);
88 
89  histo[NSIMHITS].fill(event, id, &iEvent);
90 
91  histo[RECHIT_X].fill(rechit_x, id, &iEvent);
92  histo[RECHIT_Y].fill(rechit_y, id, &iEvent);
93 
94  histo[RES_X].fill(res_x, id, &iEvent);
95  histo[RES_Y].fill(res_y, id, &iEvent);
96 
97  histo[ERROR_X].fill(lerr_x, id, &iEvent);
98  histo[ERROR_Y].fill(lerr_y, id, &iEvent);
99 
100  histo[PULL_X].fill(pull_x, id, &iEvent);
101  histo[PULL_Y].fill(pull_y, id, &iEvent);
102  }
103  }
104  }
105 }
LocalError::xx
float xx() const
Definition: LocalError.h:24
rpcPointValidation_cfi.simHit
simHit
Definition: rpcPointValidation_cfi.py:24
edmNew::DetSetVector::const_iterator
boost::transform_iterator< IterHelp, const_IdIter > const_iterator
Definition: DetSetVectorNew.h:220
edmNew::DetSetVector::end
const_iterator end(bool update=false) const
Definition: DetSetVectorNew.h:592
EncodedEventId::event
int event() const
get the contents of the subdetector field (should be protected?)
Definition: EncodedEventId.h:30
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:517
PV3DBase::y
T y() const
Definition: PV3DBase.h:63
input
static std::string const input
Definition: EdmProvDump.cc:48
PSimHit::exitPoint
Local3DPoint exitPoint() const
Exit point in the local Det frame.
Definition: PSimHit.h:46
LocalError::yy
float yy() const
Definition: LocalError.h:26
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
EncodedEventId::bunchCrossing
int bunchCrossing() const
get the detector field from this detid
Definition: EncodedEventId.h:25
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
SiPixelPhase1RecHitsV::srcToken_
edm::EDGetTokenT< SiPixelRecHitCollection > srcToken_
Definition: SiPixelPhase1RecHitsV.h:37
PSimHit::eventId
EncodedEventId eventId() const
Definition: PSimHit.h:108
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:74
SiPixelPhase1RecHitsV::trackerHitAssociatorConfig_
TrackerHitAssociator::Config trackerHitAssociatorConfig_
Definition: SiPixelPhase1RecHitsV.h:36
DetId
Definition: DetId.h:18
HistogramManagerHolder::histo
std::vector< HistogramManager > histo
Definition: SiPixelPhase1Base.h:33
PV3DBase::x
T x() const
Definition: PV3DBase.h:62
PSimHit::entryPoint
Local3DPoint entryPoint() const
Entry point in the local Det frame.
Definition: PSimHit.h:43
edmNew::DetSetVector::begin
const_iterator begin(bool update=false) const
Definition: DetSetVectorNew.h:586
event
Definition: event.py:1
SiPixelRecHit
Our base class.
Definition: SiPixelRecHit.h:23

Member Data Documentation

edm::EDGetTokenT<SiPixelRecHitCollection> SiPixelPhase1RecHitsV::srcToken_
private

Definition at line 37 of file SiPixelPhase1RecHitsV.h.

Referenced by analyze().

TrackerHitAssociator::Config SiPixelPhase1RecHitsV::trackerHitAssociatorConfig_
private

Definition at line 36 of file SiPixelPhase1RecHitsV.h.

Referenced by analyze().

Generated for CMSSW Reference Manual by doxygen 1.8.11