CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
ShallowClustersProducer.cc
Go to the documentation of this file.
1 #include "CalibTracker/SiStripCommon/interface/ShallowClustersProducer.h"
2 
3 #include "FWCore/Framework/interface/Event.h"
4 #include "FWCore/Framework/interface/ESHandle.h"
5 #include "FWCore/ParameterSet/interface/ParameterSet.h"
6 #include "RecoLocalTracker/SiStripClusterizer/interface/SiStripClusterInfo.h"
7 #include "DataFormats/SiStripDigi/interface/SiStripProcessedRawDigi.h"
8 #include "DataFormats/Common/interface/DetSetVectorNew.h"
9 #include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
10 #include "Geometry/Records/interface/TrackerTopologyRcd.h"
11 #include "boost/foreach.hpp"
12 
13 ShallowClustersProducer::ShallowClustersProducer(const edm::ParameterSet& iConfig)
14  : Prefix(iConfig.getParameter<std::string>("Prefix") )
15 {
16  produces <std::vector<unsigned> > ( Prefix + "number" );
17  produces <std::vector<unsigned> > ( Prefix + "width" );
18  produces <std::vector<float> > ( Prefix + "variance" );
19  produces <std::vector<float> > ( Prefix + "barystrip" );
20  produces <std::vector<float> > ( Prefix + "middlestrip" );
21  produces <std::vector<unsigned> > ( Prefix + "charge" );
22  produces <std::vector<float> > ( Prefix + "noise" );
23  produces <std::vector<float> > ( Prefix + "ston" );
24  produces <std::vector<unsigned> > ( Prefix + "seedstrip" );
25  produces <std::vector<unsigned> > ( Prefix + "seedindex" );
26  produces <std::vector<unsigned> > ( Prefix + "seedcharge" );
27  produces <std::vector<float> > ( Prefix + "seednoise" );
28  produces <std::vector<float> > ( Prefix + "seedgain" );
29  produces <std::vector<unsigned> > ( Prefix + "qualityisbad" );
30 
31  produces <std::vector<float> > ( Prefix + "rawchargeC" );
32  produces <std::vector<float> > ( Prefix + "rawchargeL" );
33  produces <std::vector<float> > ( Prefix + "rawchargeR" );
34  produces <std::vector<float> > ( Prefix + "rawchargeLL" );
35  produces <std::vector<float> > ( Prefix + "rawchargeRR" );
36  produces <std::vector<float> > ( Prefix + "eta" );
37  produces <std::vector<float> > ( Prefix + "foldedeta" );
38  produces <std::vector<float> > ( Prefix + "etaX" );
39  produces <std::vector<float> > ( Prefix + "etaasymm" );
40  produces <std::vector<float> > ( Prefix + "outsideasymm");
41  produces <std::vector<float> > ( Prefix + "neweta");
42  produces <std::vector<float> > ( Prefix + "newetaerr");
43 
44  produces <std::vector<unsigned> > ( Prefix + "detid" );
45  produces <std::vector<int> > ( Prefix + "subdetid" );
46  produces <std::vector<int> > ( Prefix + "module" );
47  produces <std::vector<int> > ( Prefix + "side" );
48  produces <std::vector<int> > ( Prefix + "layerwheel" );
49  produces <std::vector<int> > ( Prefix + "stringringrod" );
50  produces <std::vector<int> > ( Prefix + "petal" );
51  produces <std::vector<int> > ( Prefix + "stereo" );
52 
53  theClustersToken_ = consumes<edm::DetSetVector<SiStripCluster> > (iConfig.getParameter<edm::InputTag>("Clusters"));
54  theDigisToken_ = consumes<edm::DetSetVector<SiStripProcessedRawDigi> > (edm::InputTag("siStripProcessedRawDigis", ""));
55 }
56 
57 void ShallowClustersProducer::
58 produce(edm::Event& iEvent, const edm::EventSetup& iSetup) {
59  //Retrieve tracker topology from geometry
60  edm::ESHandle<TrackerTopology> tTopoHandle;
61  iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
62  const TrackerTopology* const tTopo = tTopoHandle.product();
63 
64  std::auto_ptr<std::vector<unsigned> > number ( new std::vector<unsigned>(7,0) );
65  std::auto_ptr<std::vector<unsigned> > width ( new std::vector<unsigned>() );
66  std::auto_ptr<std::vector<float> > variance ( new std::vector<float>() );
67  std::auto_ptr<std::vector<float> > barystrip ( new std::vector<float>() );
68  std::auto_ptr<std::vector<float> > middlestrip ( new std::vector<float>() );
69  std::auto_ptr<std::vector<unsigned> > charge ( new std::vector<unsigned>() );
70  std::auto_ptr<std::vector<float> > noise ( new std::vector<float>() );
71  std::auto_ptr<std::vector<float> > ston ( new std::vector<float>() );
72  std::auto_ptr<std::vector<unsigned> > seedstrip ( new std::vector<unsigned>() );
73  std::auto_ptr<std::vector<unsigned> > seedindex ( new std::vector<unsigned>() );
74  std::auto_ptr<std::vector<unsigned> > seedcharge ( new std::vector<unsigned>() );
75  std::auto_ptr<std::vector<float> > seednoise ( new std::vector<float>() );
76  std::auto_ptr<std::vector<float> > seedgain ( new std::vector<float>() );
77  std::auto_ptr<std::vector<unsigned> > qualityisbad ( new std::vector<unsigned>() );
78 
79  std::auto_ptr<std::vector<float> > rawchargeC ( new std::vector<float>() );
80  std::auto_ptr<std::vector<float> > rawchargeL ( new std::vector<float>() );
81  std::auto_ptr<std::vector<float> > rawchargeR ( new std::vector<float>() );
82  std::auto_ptr<std::vector<float> > rawchargeLL ( new std::vector<float>() );
83  std::auto_ptr<std::vector<float> > rawchargeRR ( new std::vector<float>() );
84  std::auto_ptr<std::vector<float> > etaX ( new std::vector<float>() );
85  std::auto_ptr<std::vector<float> > eta ( new std::vector<float>() );
86  std::auto_ptr<std::vector<float> > foldedeta ( new std::vector<float>() );
87  std::auto_ptr<std::vector<float> > etaasymm ( new std::vector<float>() );
88  std::auto_ptr<std::vector<float> > outsideasymm ( new std::vector<float>() );
89  std::auto_ptr<std::vector<float> > neweta ( new std::vector<float>() );
90  std::auto_ptr<std::vector<float> > newetaerr ( new std::vector<float>() );
91 
92  std::auto_ptr<std::vector<unsigned> > detid ( new std::vector<unsigned>() );
93  std::auto_ptr<std::vector<int> > subdetid ( new std::vector<int>() );
94  std::auto_ptr<std::vector<int> > side ( new std::vector<int>() );
95  std::auto_ptr<std::vector<int> > module ( new std::vector<int>() );
96  std::auto_ptr<std::vector<int> > layerwheel ( new std::vector<int>() );
97  std::auto_ptr<std::vector<int> > stringringrod ( new std::vector<int>() );
98  std::auto_ptr<std::vector<int> > petal ( new std::vector<int>() );
99  std::auto_ptr<std::vector<int> > stereo ( new std::vector<int>());
100 
101  edm::Handle<edmNew::DetSetVector<SiStripCluster> > clusters;
102  // iEvent.getByLabel(theClustersLabel, clusters);
103  iEvent.getByToken(theClustersToken_, clusters);
104 
105  edm::Handle<edm::DetSetVector<SiStripProcessedRawDigi> > rawProcessedDigis;
106  // iEvent.getByLabel("siStripProcessedRawDigis", "", rawProcessedDigis);
107  iEvent.getByToken(theDigisToken_,rawProcessedDigis);
108 
109  edmNew::DetSetVector<SiStripCluster>::const_iterator itClusters=clusters->begin();
110  for(;itClusters!=clusters->end();++itClusters){
111  uint32_t id = itClusters->id();
112  const moduleVars moduleV(id, tTopo);
113  for(edmNew::DetSet<SiStripCluster>::const_iterator cluster=itClusters->begin(); cluster!=itClusters->end();++cluster){
114 
115  const SiStripClusterInfo info(*cluster, iSetup, id);
116  const NearDigis digis = rawProcessedDigis.isValid() ? NearDigis(info, *rawProcessedDigis) : NearDigis(info);
117 
118  (number->at(0))++;
119  (number->at(moduleV.subdetid))++;
120  width->push_back( cluster->amplitudes().size() );
121  barystrip->push_back( cluster->barycenter() );
122  variance->push_back( info.variance() );
123  middlestrip->push_back( info.firstStrip() + info.width()/2.0 );
124  charge->push_back( info.charge() );
125  noise->push_back( info.noiseRescaledByGain() );
126  ston->push_back( info.signalOverNoise() );
127  seedstrip->push_back( info.maxStrip() );
128  seedindex->push_back( info.maxIndex() );
129  seedcharge->push_back( info.maxCharge() );
130  seednoise->push_back( info.stripNoisesRescaledByGain().at(info.maxIndex()) );
131  seedgain->push_back( info.stripGains().at(info.maxIndex()) );
132  qualityisbad->push_back( info.IsAnythingBad() );
133 
134  rawchargeC->push_back( digis.max );
135  rawchargeL->push_back( digis.left );
136  rawchargeR->push_back( digis.right );
137  rawchargeLL->push_back( digis.Lleft );
138  rawchargeRR->push_back( digis.Rright );
139  etaX->push_back( digis.etaX() );
140  eta->push_back( digis.eta() );
141  etaasymm->push_back( digis.etaasymm() );
142  outsideasymm->push_back( digis.outsideasymm() );
143  neweta->push_back( (digis.last-digis.first)/info.charge() );
144  newetaerr->push_back( (sqrt(digis.last+digis.first))/pow(info.charge(),1.5) );
145 
146  detid->push_back( id );
147  subdetid->push_back( moduleV.subdetid );
148  side->push_back( moduleV.side );
149  module->push_back( moduleV.module );
150  layerwheel->push_back( moduleV.layerwheel );
151  stringringrod->push_back( moduleV.stringringrod );
152  petal->push_back( moduleV.petal );
153  stereo->push_back( moduleV.stereo );
154  }
155  }
156 
157  iEvent.put( number, Prefix + "number" );
158  iEvent.put( width, Prefix + "width" );
159  iEvent.put( variance, Prefix + "variance" );
160  iEvent.put( barystrip, Prefix + "barystrip" );
161  iEvent.put( middlestrip, Prefix + "middlestrip" );
162  iEvent.put( charge, Prefix + "charge" );
163  iEvent.put( noise, Prefix + "noise" );
164  iEvent.put( ston, Prefix + "ston" );
165  iEvent.put( seedstrip, Prefix + "seedstrip" );
166  iEvent.put( seedindex, Prefix + "seedindex" );
167  iEvent.put( seedcharge, Prefix + "seedcharge" );
168  iEvent.put( seednoise, Prefix + "seednoise" );
169  iEvent.put( seedgain, Prefix + "seedgain" );
170  iEvent.put( qualityisbad, Prefix + "qualityisbad" );
171 
172  iEvent.put( rawchargeC, Prefix + "rawchargeC" );
173  iEvent.put( rawchargeL, Prefix + "rawchargeL" );
174  iEvent.put( rawchargeR, Prefix + "rawchargeR" );
175  iEvent.put( rawchargeLL, Prefix + "rawchargeLL" );
176  iEvent.put( rawchargeRR, Prefix + "rawchargeRR" );
177  iEvent.put( etaX, Prefix + "etaX" );
178  iEvent.put( eta, Prefix + "eta" );
179  iEvent.put( foldedeta, Prefix + "foldedeta" );
180  iEvent.put( etaasymm, Prefix + "etaasymm" );
181  iEvent.put( outsideasymm, Prefix + "outsideasymm" );
182  iEvent.put( neweta, Prefix + "neweta" );
183  iEvent.put( newetaerr, Prefix + "newetaerr" );
184 
185  iEvent.put( detid, Prefix + "detid" );
186  iEvent.put( subdetid, Prefix + "subdetid" );
187  iEvent.put( module, Prefix + "module" );
188  iEvent.put( side, Prefix + "side" );
189  iEvent.put( layerwheel, Prefix + "layerwheel" );
190  iEvent.put( stringringrod, Prefix + "stringringrod" );
191  iEvent.put( petal, Prefix + "petal" );
192  iEvent.put( stereo, Prefix + "stereo" );
193 
194 }
195 
196 ShallowClustersProducer::NearDigis::
197 NearDigis(const SiStripClusterInfo& info) {
198  max = info.maxCharge();
199  left = info.maxIndex() > uint16_t(0) ? info.stripCharges()[info.maxIndex()-1] : 0 ;
200  Lleft = info.maxIndex() > uint16_t(1) ? info.stripCharges()[info.maxIndex()-2] : 0 ;
201  right= unsigned(info.maxIndex()+1) < info.stripCharges().size() ? info.stripCharges()[info.maxIndex()+1] : 0 ;
202  Rright= unsigned(info.maxIndex()+2) < info.stripCharges().size() ? info.stripCharges()[info.maxIndex()+2] : 0 ;
203  first = info.stripCharges()[0];
204  last = info.stripCharges()[info.width()-1];
205 }
206 
207 ShallowClustersProducer::NearDigis::
208 NearDigis(const SiStripClusterInfo& info, const edm::DetSetVector<SiStripProcessedRawDigi>& rawProcessedDigis) {
209  edm::DetSetVector<SiStripProcessedRawDigi>::const_iterator digiframe = rawProcessedDigis.find(info.detId());
210  if( digiframe != rawProcessedDigis.end()) {
211  max = digiframe->data.at(info.maxStrip()).adc() ;
212  left = info.maxStrip() > uint16_t(0) ? digiframe->data.at(info.maxStrip()-1).adc() : 0 ;
213  Lleft = info.maxStrip() > uint16_t(1) ? digiframe->data.at(info.maxStrip()-2).adc() : 0 ;
214  right = unsigned(info.maxStrip()+1) < digiframe->data.size() ? digiframe->data.at(info.maxStrip()+1).adc() : 0 ;
215  Rright = unsigned(info.maxStrip()+2) < digiframe->data.size() ? digiframe->data.at(info.maxStrip()+2).adc() : 0 ;
216  first = digiframe->data.at(info.firstStrip()).adc();
217  last = digiframe->data.at(info.firstStrip()+info.width() - 1).adc();
218  } else {
219  *this = NearDigis(info);
220  }
221 }
222 
223 ShallowClustersProducer::moduleVars::
224 moduleVars(uint32_t detid, const TrackerTopology* tTopo) {
225  SiStripDetId subdet(detid);
226  subdetid = subdet.subDetector();
227  if( SiStripDetId::TIB == subdetid ) {
228 
229  module = tTopo->tibModule(detid);
230  side = tTopo->tibIsZMinusSide(detid)?-1:1;
231  layerwheel = tTopo->tibLayer(detid);
232  stringringrod = tTopo->tibString(detid);
233  stereo = tTopo->tibIsStereo(detid) ? 1 : 0;
234  } else
235  if( SiStripDetId::TID == subdetid ) {
236 
237  module = tTopo->tidModule(detid);
238  side = tTopo->tidIsZMinusSide(detid)?-1:1;
239  layerwheel = tTopo->tidWheel(detid);
240  stringringrod = tTopo->tidRing(detid);
241  stereo = tTopo->tidIsStereo(detid) ? 1 : 0;
242  } else
243  if( SiStripDetId::TOB == subdetid ) {
244 
245  module = tTopo->tobModule(detid);
246  side = tTopo->tobIsZMinusSide(detid)?-1:1;
247  layerwheel = tTopo->tobLayer(detid);
248  stringringrod = tTopo->tobRod(detid);
249  stereo = tTopo->tobIsStereo(detid) ? 1 : 0;
250  } else
251  if( SiStripDetId::TEC == subdetid ) {
252 
253  module = tTopo->tecModule(detid);
254  side = tTopo->tecIsZMinusSide(detid)?-1:1;
255  layerwheel = tTopo->tecWheel(detid);
256  stringringrod = tTopo->tecRing(detid);
257  petal = tTopo->tecPetalNumber(detid);
258  stereo = tTopo->tecIsStereo(detid) ? 1 : 0;
259  } else {
260  module = 0;
261  side = 0;
262  layerwheel=-1;
263  stringringrod = -1;
264  petal=-1;
265  }
266 }
ecalMGPA::adc
int adc(sample_type sample)
get the ADC sample (12 bits)
Definition: EcalMGPASample.h:11
SiStripClusterInfo::maxCharge
uint8_t maxCharge() const
Definition: SiStripClusterInfo.h:40
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
HLT_25ns14e33_v1_cff.noise
tuple noise
Definition: HLT_25ns14e33_v1_cff.py:2166
edmNew::DetSetVector::const_iterator
boost::transform_iterator< IterHelp, const_IdIter > const_iterator
Definition: DetSetVectorNew.h:124
edmNew::DetSetVector::end
const_iterator end(bool update=false) const
Definition: DetSetVectorNew.h:359
info
static const TGPicture * info(bool iBackgroundIsBlack)
Definition: FWCollectionSummaryWidget.cc:170
SiStripClusterInfo::firstStrip
uint16_t firstStrip() const
Definition: SiStripClusterInfo.h:28
create_public_lumi_plots.width
int width
Definition: create_public_lumi_plots.py:1008
edm::DetSetVector::find
iterator find(det_id_type id)
Definition: DetSetVector.h:294
TrackerTopology::tibLayer
unsigned int tibLayer(const DetId &id) const
Definition: TrackerTopology.h:142
TrackerTopology::tibString
unsigned int tibString(const DetId &id) const
Definition: TrackerTopology.h:344
TrackerTopology::tidRing
unsigned int tidRing(const DetId &id) const
Definition: TrackerTopology.h:223
SiStripClusterInfo::stripGains
std::vector< float > stripGains() const
Definition: SiStripClusterInfo.cc:77
SiStripClusterInfo::IsAnythingBad
bool IsAnythingBad() const
Definition: SiStripClusterInfo.cc:114
TrackerTopology::tobIsStereo
bool tobIsStereo(const DetId &id) const
Definition: TrackerTopology.h:261
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:457
SiStripClusterInfo::noiseRescaledByGain
float noiseRescaledByGain() const
Definition: SiStripClusterInfo.h:45
TrackerTopology::tecRing
unsigned int tecRing(const DetId &id) const
ring id
Definition: TrackerTopology.h:220
ShallowClustersProducer::ShallowClustersProducer
ShallowClustersProducer(const edm::ParameterSet &)
Definition: ShallowClustersProducer.cc:13
TrackerTopology::tidWheel
unsigned int tidWheel(const DetId &id) const
Definition: TrackerTopology.h:203
SiStripClusterInfo::variance
float variance() const
Definition: SiStripClusterInfo.cc:38
edmNew::DetSet::const_iterator
data_type const * const_iterator
Definition: DetSetNew.h:30
SiStripClusterInfo::maxIndex
uint16_t maxIndex() const
Definition: SiStripClusterInfo.h:41
edmNew::DetSetVector::id
id_type id(size_t cell) const
Definition: DetSetVectorNew.h:397
ShallowClustersProducer::moduleVars::moduleVars
moduleVars(uint32_t, const TrackerTopology *)
Definition: ShallowClustersProducer.cc:224
HLT_25ns14e33_v1_cff.InputTag
tuple InputTag
Definition: HLT_25ns14e33_v1_cff.py:45333
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
SiStripClusterInfo::stripCharges
auto stripCharges() const -> decltype(cluster() ->amplitudes())
Definition: SiStripClusterInfo.h:33
TrackerTopology::tecIsZMinusSide
bool tecIsZMinusSide(const DetId &id) const
Definition: TrackerTopology.h:258
TrackerTopology::tidIsStereo
bool tidIsStereo(const DetId &id) const
Definition: TrackerTopology.h:264
TrackerTopology::tidIsZMinusSide
bool tidIsZMinusSide(const DetId &id) const
Definition: TrackerTopology.h:255
TrackerTopology::tecIsStereo
bool tecIsStereo(const DetId &id) const
Definition: TrackerTopology.h:262
SiStripClusterInfo::signalOverNoise
float signalOverNoise() const
Definition: SiStripClusterInfo.h:47
HLT_25ns14e33_v1_cff.clusters
tuple clusters
Definition: HLT_25ns14e33_v1_cff.py:26184
iEvent
int iEvent
Definition: GenABIO.cc:230
SiStripClusterInfo::charge
uint16_t charge() const
Definition: SiStripClusterInfo.h:39
TrackerTopology::tibIsZMinusSide
bool tibIsZMinusSide(const DetId &id) const
Definition: TrackerTopology.h:252
edm::Event::put
OrphanHandle< PROD > put(std::auto_ptr< PROD > product)
Put a new product.
Definition: Event.h:115
EnsembleCalibrationLA_cfg.Prefix
tuple Prefix
Definition: EnsembleCalibrationLA_cfg.py:18
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
TrackerTopology::tidModule
unsigned int tidModule(const DetId &id) const
Definition: TrackerTopology.h:168
TrackerTopology::tobIsZMinusSide
bool tobIsZMinusSide(const DetId &id) const
Definition: TrackerTopology.h:249
SiStripClusterInfo::width
uint16_t width() const
Definition: SiStripClusterInfo.h:27
ShallowClustersProducer::theClustersToken_
edm::EDGetTokenT< edm::DetSetVector< SiStripCluster > > theClustersToken_
Definition: ShallowClustersProducer.h:44
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:75
ShallowClustersProducer::produce
void produce(edm::Event &, const edm::EventSetup &)
Definition: ShallowClustersProducer.cc:58
edm::DetSetVector::end
iterator end()
Return the off-the-end iterator.
Definition: DetSetVector.h:365
edm::DetSetVector::size
size_type size() const
Return the number of contained DetSets.
Definition: DetSetVector.h:287
TrackerTopology::tibModule
unsigned int tibModule(const DetId &id) const
Definition: TrackerTopology.h:165
TrackerTopology::tecModule
unsigned int tecModule(const DetId &id) const
Definition: TrackerTopology.h:162
SiStripDetId
Detector identifier class for the strip tracker.
Definition: SiStripDetId.h:17
SiStripDetId::subDetector
SubDetector subDetector() const
Definition: SiStripDetId.h:114
SiStripClusterInfo::maxStrip
uint16_t maxStrip() const
Definition: SiStripClusterInfo.h:30
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
TrackerTopology::tibIsStereo
bool tibIsStereo(const DetId &id) const
Definition: TrackerTopology.h:263
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
TrackerTopology::tobModule
unsigned int tobModule(const DetId &id) const
Definition: TrackerTopology.h:159
SiStripClusterInfo::detId
uint32_t detId() const
Definition: SiStripClusterInfo.h:26
ShallowClustersProducer::theDigisToken_
edm::EDGetTokenT< edm::DetSetVector< SiStripProcessedRawDigi > > theDigisToken_
Definition: ShallowClustersProducer.h:45
TrackerTopology::tecPetalNumber
unsigned int tecPetalNumber(const DetId &id) const
Definition: TrackerTopology.h:229
TrackerTopology::tobRod
unsigned int tobRod(const DetId &id) const
Definition: TrackerTopology.h:195
edm::DetSetVector::const_iterator
collection_type::const_iterator const_iterator
Definition: DetSetVector.h:108
TrackerTopology::tecWheel
unsigned int tecWheel(const DetId &id) const
Definition: TrackerTopology.h:200
module
Definition: vlib.h:208
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
edmNew::DetSetVector::begin
const_iterator begin(bool update=false) const
Definition: DetSetVectorNew.h:353
TrackerTopology::tobLayer
unsigned int tobLayer(const DetId &id) const
Definition: TrackerTopology.h:139
SiStripClusterInfo::stripNoisesRescaledByGain
std::vector< float > stripNoisesRescaledByGain() const
Definition: SiStripClusterInfo.cc:52
Generated for CMSSW Reference Manual by doxygen 1.8.5