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

#include <ShallowGainCalibration.h>

Inheritance diagram for ShallowGainCalibration:
edm::EDProducer edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

 ShallowGainCalibration (const edm::ParameterSet &)
 
- Public Member Functions inherited from edm::EDProducer
 EDProducer ()
 
ModuleDescription const & moduleDescription () const
 
virtual ~EDProducer ()
 
- Public Member Functions inherited from edm::ProducerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
std::vector< edm::ProductResolverIndex > const & indiciesForPutProducts (BranchType iBranchType) const
 
 ProducerBase ()
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription const &)> registrationCallback () const
 used by the fwk to register list of products More...
 
void resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)
 
virtual ~ProducerBase () noexcept(false)
 
- 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
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Member Functions

bool IsFarFromBorder (TrajectoryStateOnSurface *trajState, const uint32_t detid, const edm::EventSetup *iSetup)
 
void produce (edm::Event &, const edm::EventSetup &)
 
double thickness (DetId id)
 

Private Attributes

const edm::EDGetTokenT< TrajTrackAssociationCollection > association_token_
 
std::map< DetId, double > m_thicknessMap
 
const TrackerGeometry * m_tracker
 
std::string Prefix
 
std::string Suffix
 
const edm::EDGetTokenT< edm::View< reco::Track > > tracks_token_
 

Additional Inherited Members

- Public Types inherited from edm::EDProducer
typedef EDProducer ModuleType
 
- Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList TypeLabelList
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::EDProducer
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

Definition at line 59 of file ShallowGainCalibration.h.

Constructor & Destructor Documentation

ShallowGainCalibration::ShallowGainCalibration ( const edm::ParameterSet &  iConfig)
explicit

Definition at line 9 of file ShallowGainCalibration.cc.

References Prefix, produce(), and Suffix.

10  : tracks_token_( consumes< edm::View<reco::Track> >(iConfig.getParameter<edm::InputTag>("Tracks")) ),
11  association_token_( consumes< TrajTrackAssociationCollection >(iConfig.getParameter<edm::InputTag>("Tracks")) ),
12  Suffix ( iConfig.getParameter<std::string>("Suffix") ),
13  Prefix ( iConfig.getParameter<std::string>("Prefix") )
14 {
15  produces <std::vector<int> > ( Prefix + "trackindex" + Suffix );
16  produces <std::vector<unsigned int> > ( Prefix + "rawid" + Suffix );
17  produces <std::vector<double> > ( Prefix + "localdirx" + Suffix );
18  produces <std::vector<double> > ( Prefix + "localdiry" + Suffix );
19  produces <std::vector<double> > ( Prefix + "localdirz" + Suffix );
20  produces <std::vector<unsigned short> > ( Prefix + "firststrip" + Suffix );
21  produces <std::vector<unsigned short> > ( Prefix + "nstrips" + Suffix );
22  produces <std::vector<bool> > ( Prefix + "saturation" + Suffix );
23  produces <std::vector<bool> > ( Prefix + "overlapping" + Suffix );
24  produces <std::vector<bool> > ( Prefix + "farfromedge" + Suffix );
25  produces <std::vector<unsigned int> > ( Prefix + "charge" + Suffix );
26  produces <std::vector<double> > ( Prefix + "path" + Suffix );
27  produces <std::vector<double> > ( Prefix + "chargeoverpath" + Suffix );
28  produces <std::vector<unsigned char> > ( Prefix + "amplitude" + Suffix );
29  produces <std::vector<double> > ( Prefix + "gainused" + Suffix );
30  produces <std::vector<double> > ( Prefix + "gainusedTick" + Suffix );
31 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
edm::EDConsumerBase::consumes
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
Definition: EDConsumerBase.h:104
ShallowGainCalibration::tracks_token_
const edm::EDGetTokenT< edm::View< reco::Track > > tracks_token_
Definition: ShallowGainCalibration.h:63
ShallowGainCalibration::association_token_
const edm::EDGetTokenT< TrajTrackAssociationCollection > association_token_
Definition: ShallowGainCalibration.h:64

Member Function Documentation

bool ShallowGainCalibration::IsFarFromBorder ( TrajectoryStateOnSurface *  trajState,
const uint32_t  detid,
const edm::EventSetup *  iSetup 
)
private

Definition at line 255 of file ShallowGainCalibration.cc.

References Surface::bounds(), gather_cfg::cout, edm::EventSetup::get(), TrackerGeometry::idToDetUnit(), Bounds::length(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localPosition(), metProducer_cfi::parameters, LocalTrajectoryError::positionError(), GeomDet::surface(), PV3DBase< T, PVType, FrameType >::y(), and LocalError::yy().

Referenced by produce().

256 {
257  edm::ESHandle<TrackerGeometry> tkGeom; iSetup->get<TrackerDigiGeometryRecord>().get( tkGeom );
258 
259  LocalPoint HitLocalPos = trajState->localPosition();
260  LocalError HitLocalError = trajState->localError().positionError() ;
261 
262  const GeomDetUnit* it = tkGeom->idToDetUnit(DetId(detid));
263  if (dynamic_cast<const StripGeomDetUnit*>(it)==0 && dynamic_cast<const PixelGeomDetUnit*>(it)==0) {
264  std::cout << "this detID doesn't seem to belong to the Tracker" << std::endl;
265  return false;
266  }
267 
268  const BoundPlane plane = it->surface();
269  const TrapezoidalPlaneBounds* trapezoidalBounds( dynamic_cast<const TrapezoidalPlaneBounds*>(&(plane.bounds())));
270  const RectangularPlaneBounds* rectangularBounds( dynamic_cast<const RectangularPlaneBounds*>(&(plane.bounds())));
271 
272  double DistFromBorder = 1.0;
273  double HalfLength = it->surface().bounds().length() /2.0;
274 
275  if(trapezoidalBounds)
276  {
277  std::array<const float, 4> const & parameters = (*trapezoidalBounds).parameters();
278  HalfLength = parameters[3];
279  }else if(rectangularBounds){
280  HalfLength = it->surface().bounds().length() /2.0;
281  }else{return false;}
282 
283  if (fabs(HitLocalPos.y())+HitLocalError.yy() >= (HalfLength - DistFromBorder) ) return false;
284 
285  return true;
286 }
Bounds::length
virtual float length() const =0
metProducer_cfi.parameters
parameters
Definition: metProducer_cfi.py:59
PV3DBase::y
T y() const
Definition: PV3DBase.h:63
Surface::bounds
const Bounds & bounds() const
Definition: Surface.h:120
GeomDet::surface
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:42
LocalTrajectoryError::positionError
LocalError positionError() const
Definition: LocalTrajectoryError.h:83
TrackerGeometry::idToDetUnit
const TrackerGeomDet * idToDetUnit(DetId) const override
Return the pointer to the GeomDetUnit corresponding to a given DetId.
Definition: TrackerGeometry.cc:202
LocalError::yy
float yy() const
Definition: LocalError.h:26
TrajectoryStateOnSurface::localError
const LocalTrajectoryError & localError() const
Definition: TrajectoryStateOnSurface.h:119
DetId
Definition: DetId.h:18
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
gather_cfg.cout
cout
Definition: gather_cfg.py:145
void ShallowGainCalibration::produce ( edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
)
private

Definition at line 34 of file ShallowGainCalibration.cc.

References a, CustomPhysics_cfi::amplitude, SiStripCluster::amplitudes(), TrackValidation_cff::association, association_token_, genericTrackCandidates_cff::associations, edm::AssociationMap< Tag >::begin(), ALCARECOTkAlJpsiMuMu_cff::charge, SiStripRecHit1D::cluster(), SiStripRecHit2D::cluster(), SiPixelRecHit::cluster(), edm::AssociationMap< edm::OneToOne< std::vector< Trajectory >, reco::TrackCollection, unsigned short > >::const_iterator, edm::AssociationMap< Tag >::end(), shallow::findTrackIndex(), SiStripCluster::firstStrip(), TrackingRecHit::geographicalId(), edm::EventSetup::get(), SiStripGain::getApvGain(), edm::Event::getByToken(), SiStripGain::getRange(), IsFarFromBorder(), edm::ESHandleBase::isValid(), TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localDirection(), m_tracker, PV3DBase< T, PVType, FrameType >::mag(), Trajectory::measurements(), SiPixelCluster::minPixelCol(), SiPixelCluster::minPixelRow(), SiStripMatchedRecHit2D::monoCluster(), SiStripMatchedRecHit2D::monoId(), eostools::move(), NULL, Overlapping, callgraph::path, SiPixelCluster::pixelADC(), Prefix, edm::Event::put(), DetId::rawId(), SiStripMatchedRecHit2D::stereoCluster(), SiStripMatchedRecHit2D::stereoId(), Suffix, thickness(), HiIsolationCommonParameters_cff::track, l1t::tracks, tracks_token_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), HTMLExport.HTMLExportStatic::export(), and ShallowGainCalibration().

34  {
35  auto trackindex = std::make_unique<std::vector<int>> ();
36  auto rawid = std::make_unique<std::vector<unsigned int>> ();
37  auto localdirx = std::make_unique<std::vector<double>> ();
38  auto localdiry = std::make_unique<std::vector<double>> ();
39  auto localdirz = std::make_unique<std::vector<double>> ();
40  auto firststrip = std::make_unique<std::vector<unsigned short>> ();
41  auto nstrips = std::make_unique<std::vector<unsigned short>> ();
42  auto saturation = std::make_unique<std::vector<bool>> ();
43  auto overlapping = std::make_unique<std::vector<bool>> ();
44  auto farfromedge = std::make_unique<std::vector<bool>> ();
45  auto charge = std::make_unique<std::vector<unsigned int>> ();
46  auto path = std::make_unique<std::vector<double>> ();
47  auto chargeoverpath= std::make_unique<std::vector<double>> ();
48  auto amplitude = std::make_unique<std::vector<unsigned char>> ();
49  auto gainused = std::make_unique<std::vector<double>> ();
50  auto gainusedTick = std::make_unique<std::vector<double>> ();
51 
52  edm::ESHandle<TrackerGeometry> theTrackerGeometry; iSetup.get<TrackerDigiGeometryRecord>().get( theTrackerGeometry );
53  m_tracker=&(* theTrackerGeometry );
54  edm::ESHandle<SiStripGain> gainHandle; iSetup.get<SiStripGainRcd>().get(gainHandle);
55  edm::Handle<edm::View<reco::Track> > tracks; iEvent.getByToken(tracks_token_, tracks);
56  edm::Handle<TrajTrackAssociationCollection> associations; iEvent.getByToken(association_token_, associations);
57 
58  for( TrajTrackAssociationCollection::const_iterator association = associations->begin(); association != associations->end(); association++) {
59  const Trajectory* traj = association->key.get();
60  const reco::Track* track = association->val.get();
61 
62  vector<TrajectoryMeasurement> measurements = traj->measurements();
63  for(vector<TrajectoryMeasurement>::const_iterator measurement_it = measurements.begin(); measurement_it!=measurements.end(); measurement_it++){
64  TrajectoryStateOnSurface trajState = measurement_it->updatedState();
65  if( !trajState.isValid() ) continue;
66 
67  const TrackingRecHit* hit = (*measurement_it->recHit()).hit();
68  const SiStripRecHit1D* sistripsimple1dhit = dynamic_cast<const SiStripRecHit1D*>(hit);
69  const SiStripRecHit2D* sistripsimplehit = dynamic_cast<const SiStripRecHit2D*>(hit);
70  const SiStripMatchedRecHit2D* sistripmatchedhit = dynamic_cast<const SiStripMatchedRecHit2D*>(hit);
71  const SiPixelRecHit* sipixelhit = dynamic_cast<const SiPixelRecHit*>(hit);
72 
73  const SiPixelCluster* PixelCluster = NULL;
74  const SiStripCluster* StripCluster = NULL;
75  uint32_t DetId = 0;
76 
77  for(unsigned int h=0;h<2;h++){
78  if(!sistripmatchedhit && h==1){
79  continue;
80  }else if(sistripmatchedhit && h==0){
81  StripCluster = &sistripmatchedhit->monoCluster();
82  DetId = sistripmatchedhit->monoId();
83  }else if(sistripmatchedhit && h==1){
84  StripCluster = &sistripmatchedhit->stereoCluster();;
85  DetId = sistripmatchedhit->stereoId();
86  }else if(sistripsimplehit){
87  StripCluster = (sistripsimplehit->cluster()).get();
88  DetId = sistripsimplehit->geographicalId().rawId();
89  }else if(sistripsimple1dhit){
90  StripCluster = (sistripsimple1dhit->cluster()).get();
91  DetId = sistripsimple1dhit->geographicalId().rawId();
92  }else if(sipixelhit){
93  PixelCluster = (sipixelhit->cluster()).get();
94  DetId = sipixelhit->geographicalId().rawId();
95  }else{
96  continue;
97  }
98 
99  LocalVector trackDirection = trajState.localDirection();
100  double cosine = trackDirection.z()/trackDirection.mag();
101  bool Saturation = false;
102  bool Overlapping = false;
103  unsigned int Charge = 0;
104  double Path = (10.0*thickness(DetId))/fabs(cosine);
105  double PrevGain = -1;
106  double PrevGainTick = -1;
107  int FirstStrip = 0;
108  int NStrips = 0;
109 
110  if(StripCluster){
111  const auto & Ampls = StripCluster->amplitudes();
112  FirstStrip = StripCluster->firstStrip();
113  NStrips = Ampls.size();
114  int APVId = FirstStrip/128;
115 
116 
117  if(gainHandle.isValid()){
118  PrevGain = gainHandle->getApvGain(APVId,gainHandle->getRange(DetId, 1),1);
119  PrevGainTick = gainHandle->getApvGain(APVId,gainHandle->getRange(DetId, 0),1);
120  }
121 
122  for(unsigned int a=0;a<Ampls.size();a++){
123  Charge+=Ampls[a];
124  if(Ampls[a] >=254)Saturation =true;
125  amplitude->push_back( Ampls[a] );
126  }
127 
128  if(FirstStrip==0 )Overlapping=true;
129  if(FirstStrip==128 )Overlapping=true;
130  if(FirstStrip==256 )Overlapping=true;
131  if(FirstStrip==384 )Overlapping=true;
132  if(FirstStrip==512 )Overlapping=true;
133  if(FirstStrip==640 )Overlapping=true;
134 
135  if(FirstStrip<=127 && FirstStrip+Ampls.size()>127)Overlapping=true;
136  if(FirstStrip<=255 && FirstStrip+Ampls.size()>255)Overlapping=true;
137  if(FirstStrip<=383 && FirstStrip+Ampls.size()>383)Overlapping=true;
138  if(FirstStrip<=511 && FirstStrip+Ampls.size()>511)Overlapping=true;
139  if(FirstStrip<=639 && FirstStrip+Ampls.size()>639)Overlapping=true;
140 
141  if(FirstStrip+Ampls.size()==127 )Overlapping=true;
142  if(FirstStrip+Ampls.size()==255 )Overlapping=true;
143  if(FirstStrip+Ampls.size()==383 )Overlapping=true;
144  if(FirstStrip+Ampls.size()==511 )Overlapping=true;
145  if(FirstStrip+Ampls.size()==639 )Overlapping=true;
146  if(FirstStrip+Ampls.size()==767 )Overlapping=true;
147  }else if(PixelCluster){
148  const auto& Ampls = PixelCluster->pixelADC();
149  int FirstRow = PixelCluster->minPixelRow();
150  int FirstCol = PixelCluster->minPixelCol();
151  FirstStrip = ((FirstRow/80)<<3 | (FirstCol/52)) * 128; //Hack to save the APVId
152  NStrips = 0;
153  Saturation = false;
154  Overlapping = false;
155 
156  for(unsigned int a=0;a<Ampls.size();a++){
157  Charge+=Ampls[a];
158  if(Ampls[a] >=254)Saturation =true;
159  }
160  }
161  double ChargeOverPath = (double)Charge / Path ;
162 
163  trackindex ->push_back( shallow::findTrackIndex(tracks, track) );
164  rawid ->push_back( DetId );
165  localdirx ->push_back( trackDirection.x() );
166  localdiry ->push_back( trackDirection.y() );
167  localdirz ->push_back( trackDirection.z() );
168  firststrip ->push_back( FirstStrip );
169  nstrips ->push_back( NStrips );
170  saturation ->push_back( Saturation );
171  overlapping ->push_back( Overlapping );
172  farfromedge ->push_back( StripCluster ? IsFarFromBorder(&trajState,DetId, &iSetup) : true );
173  charge ->push_back( Charge );
174  path ->push_back( Path );
175  chargeoverpath->push_back( ChargeOverPath );
176  gainused ->push_back( PrevGain );
177  gainusedTick ->push_back( PrevGainTick );
178  }
179  }
180  }
181 
182  iEvent.put(std::move(trackindex), Prefix + "trackindex" + Suffix );
183  iEvent.put(std::move(rawid ), Prefix + "rawid" + Suffix );
184  iEvent.put(std::move(localdirx ), Prefix + "localdirx" + Suffix );
185  iEvent.put(std::move(localdiry ), Prefix + "localdiry" + Suffix );
186  iEvent.put(std::move(localdirz ), Prefix + "localdirz" + Suffix );
187  iEvent.put(std::move(firststrip), Prefix + "firststrip" + Suffix );
188  iEvent.put(std::move(nstrips), Prefix + "nstrips" + Suffix );
189  iEvent.put(std::move(saturation), Prefix + "saturation" + Suffix );
190  iEvent.put(std::move(overlapping), Prefix + "overlapping" + Suffix );
191  iEvent.put(std::move(farfromedge), Prefix + "farfromedge" + Suffix );
192  iEvent.put(std::move(charge), Prefix + "charge" + Suffix );
193  iEvent.put(std::move(path), Prefix + "path" + Suffix );
194  iEvent.put(std::move(chargeoverpath),Prefix + "chargeoverpath"+ Suffix );
195  iEvent.put(std::move(amplitude), Prefix + "amplitude" + Suffix );
196  iEvent.put(std::move(gainused), Prefix + "gainused" + Suffix );
197  iEvent.put(std::move(gainusedTick), Prefix + "gainusedTick" + Suffix );
198 }
SiStripRecHit1D::cluster
ClusterRef cluster() const
Definition: SiStripRecHit1D.h:26
SiPixelCluster::minPixelCol
int minPixelCol() const
Definition: SiPixelCluster.h:140
edm::Event::put
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:122
SiStripMatchedRecHit2D::stereoId
unsigned int stereoId() const
Definition: SiStripMatchedRecHit2D.h:23
edm::AssociationMap::end
const_iterator end() const
last iterator over the map (read only)
Definition: AssociationMap.h:160
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:460
shallow::findTrackIndex
int findTrackIndex(const edm::Handle< edm::View< reco::Track > > &h, const reco::Track *t)
Definition: ShallowTools.cc:28
SiStripMatchedRecHit2D::monoCluster
SiStripCluster const & monoCluster() const
Definition: SiStripMatchedRecHit2D.h:39
TrajectoryStateOnSurface::localDirection
LocalVector localDirection() const
Definition: TrajectoryStateOnSurface.h:116
PV3DBase::y
T y() const
Definition: PV3DBase.h:63
genericTrackCandidates_cff.associations
associations
Definition: genericTrackCandidates_cff.py:45
NULL
#define NULL
Definition: scimark2.h:8
SiStripCluster::firstStrip
uint16_t firstStrip() const
Definition: SiStripCluster.h:46
CustomPhysics_cfi.amplitude
amplitude
Definition: CustomPhysics_cfi.py:12
ShallowGainCalibration::IsFarFromBorder
bool IsFarFromBorder(TrajectoryStateOnSurface *trajState, const uint32_t detid, const edm::EventSetup *iSetup)
Definition: ShallowGainCalibration.cc:255
SiStripGain::getApvGain
static float getApvGain(const uint16_t &apv, const SiStripApvGain::Range &range)
Definition: SiStripGain.h:69
DetId::rawId
uint32_t rawId() const
get the raw id
Definition: DetId.h:43
Trajectory::measurements
DataContainer const & measurements() const
Definition: Trajectory.h:196
PV3DBase::mag
T mag() const
Definition: PV3DBase.h:67
SiPixelCluster::minPixelRow
int minPixelRow() const
Definition: SiPixelCluster.h:138
edm::Path
Definition: Path.h:44
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
SiStripRecHit2D::cluster
ClusterRef cluster() const
Definition: SiStripRecHit2D.h:29
ALCARECOTkAlJpsiMuMu_cff.charge
charge
Definition: ALCARECOTkAlJpsiMuMu_cff.py:47
ShallowGainCalibration::tracks_token_
const edm::EDGetTokenT< edm::View< reco::Track > > tracks_token_
Definition: ShallowGainCalibration.h:63
SiPixelCluster::pixelADC
const std::vector< uint16_t > & pixelADC() const
Definition: SiPixelCluster.h:145
DetId
Definition: DetId.h:18
ShallowGainCalibration::m_tracker
const TrackerGeometry * m_tracker
Definition: ShallowGainCalibration.h:75
HiIsolationCommonParameters_cff.track
track
Definition: HiIsolationCommonParameters_cff.py:9
SiPixelRecHit::cluster
ClusterRef cluster() const
Definition: SiPixelRecHit.h:49
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
SiStripMatchedRecHit2D::stereoCluster
SiStripCluster const & stereoCluster() const
Definition: SiStripMatchedRecHit2D.h:36
ShallowGainCalibration::association_token_
const edm::EDGetTokenT< TrajTrackAssociationCollection > association_token_
Definition: ShallowGainCalibration.h:64
SiPixelCluster
Pixel cluster – collection of neighboring pixels above threshold.
Definition: SiPixelCluster.h:27
a
double a
Definition: hdecay.h:121
TrackValidation_cff.association
association
Definition: TrackValidation_cff.py:268
SiStripMatchedRecHit2D::monoId
unsigned int monoId() const
Definition: SiStripMatchedRecHit2D.h:24
edm::AssociationMap::begin
const_iterator begin() const
first iterator over the map (read only)
Definition: AssociationMap.h:158
TrackingRecHit::geographicalId
DetId geographicalId() const
Definition: TrackingRecHit.h:116
edm::ESHandleBase::isValid
bool isValid() const
Definition: ESHandle.h:47
PV3DBase::x
T x() const
Definition: PV3DBase.h:62
SiStripCluster::amplitudes
const std::vector< uint8_t > & amplitudes() const
Definition: SiStripCluster.h:59
eostools.move
def move(src, dest)
Definition: eostools.py:510
SiPixelRecHit
Our base class.
Definition: SiPixelRecHit.h:23
SiStripGain::getRange
const SiStripApvGain::Range getRange(uint32_t detID) const
Definition: SiStripGain.h:66
double ShallowGainCalibration::thickness ( DetId  id)
private

Definition at line 289 of file ShallowGainCalibration.cc.

References hcalTTPDigis_cfi::id, TrackerGeometry::idToDetUnit(), fastTrackerRecHitType::isPixel(), m_thicknessMap, and m_tracker.

Referenced by produce().

290 {
291  map<DetId,double>::iterator th=m_thicknessMap.find(id);
292  if(th!=m_thicknessMap.end())
293  return (*th).second;
294  else {
295  double detThickness=1.;
296  //compute thickness normalization
297  const GeomDetUnit* it = m_tracker->idToDetUnit(DetId(id));
298  bool isPixel = dynamic_cast<const PixelGeomDetUnit*>(it)!=0;
299  bool isStrip = dynamic_cast<const StripGeomDetUnit*>(it)!=0;
300  if (!isPixel && ! isStrip) {
301  //FIXME throw exception
302  edm::LogWarning("DeDxHitsProducer") << "\t\t this detID doesn't seem to belong to the Tracker";
303  detThickness = 1.;
304  }else{
305  detThickness = it->surface().bounds().thickness();
306  }
307 
308  m_thicknessMap[id]=detThickness;//computed value
309  return detThickness;
310  }
311 
312 }
ShallowGainCalibration::m_thicknessMap
std::map< DetId, double > m_thicknessMap
Definition: ShallowGainCalibration.h:76
TrackerGeometry::idToDetUnit
const TrackerGeomDet * idToDetUnit(DetId) const override
Return the pointer to the GeomDetUnit corresponding to a given DetId.
Definition: TrackerGeometry.cc:202
DetId
Definition: DetId.h:18
ShallowGainCalibration::m_tracker
const TrackerGeometry * m_tracker
Definition: ShallowGainCalibration.h:75
fastTrackerRecHitType::isPixel
bool isPixel(HitType hitType)
Definition: FastTrackerRecHit.h:32
hcalTTPDigis_cfi.id
id
Definition: hcalTTPDigis_cfi.py:7

Member Data Documentation

const edm::EDGetTokenT< TrajTrackAssociationCollection > ShallowGainCalibration::association_token_
private

Definition at line 64 of file ShallowGainCalibration.h.

Referenced by produce().

std::map<DetId,double> ShallowGainCalibration::m_thicknessMap
private

Definition at line 76 of file ShallowGainCalibration.h.

Referenced by thickness().

const TrackerGeometry* ShallowGainCalibration::m_tracker
private

Definition at line 75 of file ShallowGainCalibration.h.

Referenced by produce(), and thickness().

std::string ShallowGainCalibration::Prefix
private

Definition at line 67 of file ShallowGainCalibration.h.

Referenced by produce(), and ShallowGainCalibration().

std::string ShallowGainCalibration::Suffix
private

Definition at line 66 of file ShallowGainCalibration.h.

Referenced by produce(), and ShallowGainCalibration().

const edm::EDGetTokenT< edm::View<reco::Track> > ShallowGainCalibration::tracks_token_
private

Definition at line 63 of file ShallowGainCalibration.h.

Referenced by produce().

Generated for CMSSW Reference Manual by doxygen 1.8.11