CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
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)
 
 ProducerBase ()
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription
const &)> 		registrationCallback () const
 used by the fwk to register list of products More...
 
virtual ~ProducerBase ()
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
 EDConsumerBase ()
 
ProductHolderIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
std::vector
< ProductHolderIndexAndSkipBit >
const & 		itemsToGetFromEvent () const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
bool registeredToConsume (ProductHolderIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
 
virtual ~EDConsumerBase ()
 

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

std::map< DetId, double > m_thicknessMap
 
const TrackerGeometry * m_tracker
 
std::string Prefix
 
std::string Suffix
 
edm::InputTag theTracksLabel
 

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 12 of file ShallowGainCalibration.cc.

References Prefix, and Suffix.

13  : theTracksLabel( iConfig.getParameter<edm::InputTag>("Tracks") ),
14  Suffix ( iConfig.getParameter<std::string>("Suffix") ),
15  Prefix ( iConfig.getParameter<std::string>("Prefix") )
16 {
17  produces <std::vector<int> > ( Prefix + "trackindex" + Suffix );
18  produces <std::vector<unsigned int> > ( Prefix + "rawid" + Suffix );
19  produces <std::vector<float> > ( Prefix + "localdirx" + Suffix );
20  produces <std::vector<float> > ( Prefix + "localdiry" + Suffix );
21  produces <std::vector<float> > ( Prefix + "localdirz" + Suffix );
22  produces <std::vector<unsigned short> > ( Prefix + "firststrip" + Suffix );
23  produces <std::vector<unsigned short> > ( Prefix + "nstrips" + Suffix );
24  produces <std::vector<bool> > ( Prefix + "saturation" + Suffix );
25  produces <std::vector<bool> > ( Prefix + "overlapping" + Suffix );
26  produces <std::vector<bool> > ( Prefix + "farfromedge" + Suffix );
27  produces <std::vector<unsigned int> > ( Prefix + "charge" + Suffix );
28  produces <std::vector<float> > ( Prefix + "path" + Suffix );
29  produces <std::vector<float> > ( Prefix + "chargeoverpath" + Suffix );
30  produces <std::vector<unsigned char> > ( Prefix + "amplitude" + Suffix );
31  produces <std::vector<double> > ( Prefix + "gainused" + Suffix );
32  produces <std::vector<double> > ( Prefix + "gainusedTick" + Suffix );
33 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255

Member Function Documentation

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

Definition at line 257 of file ShallowGainCalibration.cc.

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

Referenced by produce().

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

Implements edm::EDProducer.

Definition at line 36 of file ShallowGainCalibration.cc.

References a, SiStripCluster::amplitudes(), RecoTauCleanerPlugins::charge, SiStripRecHit1D::cluster(), SiStripRecHit2D::cluster(), SiPixelRecHit::cluster(), edm::AssociationMap< edm::OneToOne< std::vector< Trajectory >, reco::TrackCollection, unsigned short > >::const_iterator, shallow::findTrackIndex(), SiStripCluster::firstStrip(), TrackingRecHit::geographicalId(), edm::EventSetup::get(), edm::Event::getByLabel(), h, 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(), NULL, fed_dqm_sourceclient-live_cfg::path, SiPixelCluster::pixelADC(), Prefix, edm::Event::put(), DetId::rawId(), SiStripMatchedRecHit2D::stereoCluster(), SiStripMatchedRecHit2D::stereoId(), Suffix, theTracksLabel, thickness(), testEve_cfg::tracks, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

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

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

Definition at line 291 of file ShallowGainCalibration.cc.

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

Referenced by produce().

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

Member Data Documentation

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

Definition at line 74 of file ShallowGainCalibration.h.

Referenced by thickness().

const TrackerGeometry* ShallowGainCalibration::m_tracker
private

Definition at line 73 of file ShallowGainCalibration.h.

Referenced by produce(), and thickness().

std::string ShallowGainCalibration::Prefix
private

Definition at line 65 of file ShallowGainCalibration.h.

Referenced by produce(), and ShallowGainCalibration().

std::string ShallowGainCalibration::Suffix
private

Definition at line 64 of file ShallowGainCalibration.h.

Referenced by produce(), and ShallowGainCalibration().

edm::InputTag ShallowGainCalibration::theTracksLabel
private

Definition at line 63 of file ShallowGainCalibration.h.

Referenced by produce().

Generated for CMSSW Reference Manual by doxygen 1.8.5