75 bool isHalfModule =
false;
118 meX_ = iBooker.
book1D(hid,
"Cluster barycenter X (row #)",200,0.,200.);
122 meY_ = iBooker.
book1D(hid,
"Cluster barycenter Y (column #)",500,0.,500.);
126 meSizeX_ = iBooker.
book1D(hid,
"Cluster x-width (rows)",10,0.,10.);
130 meSizeY_ = iBooker.
book1D(hid,
"Cluster y-width (columns)",15,0.,15.);
148 delete theHistogramId;
152 if(barrel && type==7){
153 hid = src.
label() +
"_Barrel";
154 meSizeYvsEtaBarrel_= iBooker.
book2D(
"sizeYvsEta_" + hid,
"Cluster size along beamline vs. Cluster position #eta",60,-3.,3.,40,0.,40.);
158 if(type==1 && barrel){
161 char sladder[80]; sprintf(sladder,
"Ladder_%02i",DBladder);
162 hid = src.
label() +
"_" + sladder;
163 if(isHalfModule) hid +=
"H";
188 meXLad_ = iBooker.
book1D(
"x_" + hid,
"Cluster barycenter X (row #)",200,0.,200.);
191 meYLad_ = iBooker.
book1D(
"y_" + hid,
"Cluster barycenter Y (column #)",500,0.,500.);
197 meSizeYLad_ = iBooker.
book1D(
"sizeY_" + hid,
"Cluster y-width (columns)",15,0.,15.);
214 if(type==2 && barrel){
218 char slayer[80]; sprintf(slayer,
"Layer_%i",DBlayer);
219 hid = src.
label() +
"_" + slayer;
243 meXLay_ = iBooker.
book1D(
"x_" + hid,
"Cluster barycenter X (row #)",200,0.,200.);
246 meYLay_ = iBooker.
book1D(
"y_" + hid,
"Cluster barycenter Y (column #)",500,0.,500.);
252 meSizeYLay_ = iBooker.
book1D(
"sizeY_" + hid,
"Cluster y-width (columns)",15,0.,15.);
276 if(type==3 && barrel){
279 char smodule[80]; sprintf(smodule,
"Ring_%i",DBmodule);
280 hid = src.
label() +
"_" + smodule;
304 meXPhi_ = iBooker.
book1D(
"x_" + hid,
"Cluster barycenter X (row #)",200,0.,200.);
307 meYPhi_ = iBooker.
book1D(
"y_" + hid,
"Cluster barycenter Y (column #)",500,0.,500.);
313 meSizeYPhi_ = iBooker.
book1D(
"sizeY_" + hid,
"Cluster y-width (columns)",15,0.,15.);
338 if(type==4 && endcap){
342 char sblade[80]; sprintf(sblade,
"Blade_%02i",blade);
343 hid = src.
label() +
"_" + sblade;
367 meXBlade_ = iBooker.
book1D(
"x_" + hid,
"Cluster barycenter X (row #)",200,0.,200.);
370 meYBlade_ = iBooker.
book1D(
"y_" + hid,
"Cluster barycenter Y (column #)",500,0.,500.);
380 if(type==5 && endcap){
384 char sdisk[80]; sprintf(sdisk,
"Disk_%i",disk);
385 hid = src.
label() +
"_" + sdisk;
409 meXDisk_ = iBooker.
book1D(
"x_" + hid,
"Cluster barycenter X (row #)",200,0.,200.);
412 meYDisk_ = iBooker.
book1D(
"y_" + hid,
"Cluster barycenter Y (column #)",500,0.,500.);
423 if(type==6 && endcap){
429 char slab[80]; sprintf(slab,
"Panel_%i_Ring_%i",panel, module);
430 hid = src.
label() +
"_" + slab;
454 meXRing_ = iBooker.
book1D(
"x_" + hid,
"Cluster barycenter X (row #)",200,0.,200.);
457 meYRing_ = iBooker.
book1D(
"y_" + hid,
"Cluster barycenter Y (column #)",500,0.,500.);
484 int SiPixelClusterModule::fill(
const edmNew::DetSetVector<SiPixelCluster>&
input,
const TrackerGeometry*
tracker,std::vector<MonitorElement*>&
layers,std::vector<MonitorElement*>& diskspz,std::vector<MonitorElement*>& disksmz,
bool modon,
bool ladon,
bool layon,
bool phion,
bool bladeon,
bool diskon,
bool ringon,
bool twoD,
bool reducedSet,
bool smileyon,
bool isUpgrade) {
490 unsigned int numberOfClusters = 0;
491 unsigned int numberOfFpixClusters = 0;
493 if( isearch != input.
end() ) {
498 for(di = isearch->
begin(); di != isearch->
end(); di++) {
500 if(endcap) numberOfFpixClusters++;
501 float charge = 0.001*(di->charge());
505 int sizeX = di->sizeX();
506 int sizeY = di->sizeY();
507 int minPixelRow = di->minPixelRow();
508 int maxPixelRow = di->maxPixelRow();
509 int minPixelCol = di->minPixelCol();
510 int maxPixelCol = di->maxPixelCol();
522 if (!(DBlayer > layers.size()) && (layers[DBlayer-1])) layers[DBlayer-1]->
Fill(clustgp.
z(),clustgp.
phi());
526 if (!(DBdisk > diskspz.size()) && (diskspz[DBdisk-1])) diskspz[DBdisk-1]->
Fill(clustgp.
x(),clustgp.
y());
528 if (!(DBdisk > disksmz.size()) && (disksmz[DBdisk-1])) disksmz[DBdisk-1]->
Fill(clustgp.
x(),clustgp.
y());
548 if(barrel && smileyon){
612 if(bladeon && endcap){
627 if(diskon && endcap){
643 if(ringon && endcap){
678 return numberOfFpixClusters;
int plaquetteName() const
plaquetteId (in pannel)
MonitorElement * meMinColLay_
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
T getParameter(std::string const &) const
MonitorElement * meMinRowRing_
boost::transform_iterator< IterHelp, const_IdIter > const_iterator
const_iterator end(bool update=false) const
MonitorElement * meMaxColLad_
MonitorElement * meSizePhi_
MonitorElement * meCharge_
MonitorElement * mePixClustersLad_px_
std::vector< LayerSetAndLayers > layers(const SeedingLayerSetsHits &sets)
int fill(const edmNew::DetSetVector< SiPixelCluster > &input, const TrackerGeometry *tracker, std::vector< MonitorElement * > &layers, std::vector< MonitorElement * > &diskspz, std::vector< MonitorElement * > &disksmz, bool modon=true, bool ladon=false, bool layon=false, bool phion=false, bool bladeon=false, bool diskon=false, bool ringon=false, bool twoD=true, bool reducedSet=false, bool smileyon=false, bool isUpgrade=false)
Fill histograms.
MonitorElement * meMaxCol_
virtual LocalPoint localPosition(const MeasurementPoint &) const =0
MonitorElement * meMaxColPhi_
MonitorElement * meNClustersLad_
MonitorElement * meYBlade_
int moduleName() const
module id (index in z)
MonitorElement * meMinRowLad_
MonitorElement * meSizeYLay_
const TrackerTopology * pTT
MonitorElement * meSizeYLad_
MonitorElement * meSizeXLad_
MonitorElement * meSizeXLay_
MonitorElement * meMinRow_
Geom::Phi< T > phi() const
MonitorElement * meXDisk_
MonitorElement * meMaxRowRing_
MonitorElement * mePixClustersLay_
MonitorElement * meSizeYBlade_
SiPixelClusterModule()
Default constructor.
data_type const * const_iterator
MonitorElement * meMinRowBlade_
MonitorElement * meNClustersPhi_
MonitorElement * meChargePhi_
MonitorElement * meMinColPhi_
MonitorElement * meNClustersBlade_
MonitorElement * mePixClustersPhi_py_
MonitorElement * meMinRowDisk_
const Plane & surface() const
The nominal surface of the GeomDet.
MonitorElement * mePixClustersRing_px_
MonitorElement * meMaxRowLad_
static std::string const input
MonitorElement * meMinRowLay_
MonitorElement * meMinColBlade_
void book(const edm::ParameterSet &iConfig, const edm::EventSetup &iSetup, DQMStore::IBooker &iBooker, int type=0, bool twoD=true, bool reducedSet=false, bool isUpgrade=false)
Book histograms.
MonitorElement * meXRing_
MonitorElement * meSizeYDisk_
MonitorElement * meSizeDisk_
MonitorElement * meMaxColRing_
MonitorElement * meMinColLad_
int bladeName() const
blade id
MonitorElement * meMaxRowBlade_
MonitorElement * meChargeRing_
Measurement2DPoint MeasurementPoint
Measurement points are two-dimensional by default.
MonitorElement * mePixClustersLad_
bool isHalfModule() const
full or half module
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
~SiPixelClusterModule()
Destructor.
MonitorElement * meMaxRowPhi_
MonitorElement * meMinRowPhi_
MonitorElement * meMinColDisk_
MonitorElement * meMinCol_
MonitorElement * mePixClusters_py_
MonitorElement * book1D(Args &&...args)
MonitorElement * mePixClustersLay_px_
MonitorElement * meNClustersDisk_
MonitorElement * meSizeYRing_
MonitorElement * mePixClustersRing_
MonitorElement * meSizeX_
MonitorElement * meSizeY_
MonitorElement * meMaxRowLay_
MonitorElement * meChargeDisk_
MonitorElement * meMaxColBlade_
MonitorElement * mePixClustersPhi_
int subdetId() const
get the contents of the subdetector field (not cast into any detector's numbering enum) ...
MonitorElement * meChargeBlade_
MonitorElement * mePixClustersLay_py_
MonitorElement * meMinColRing_
MonitorElement * meNClustersLay_
MonitorElement * meMaxColLay_
std::string setHistoId(std::string variable, uint32_t &rawId)
Set Histogram Id.
MonitorElement * meMaxRowDisk_
int ladderName() const
ladder id (index in phi)
MonitorElement * meYRing_
MonitorElement * book2D(Args &&...args)
MonitorElement * meNClustersRing_
MonitorElement * meSizeYvsEtaBarrel_
MonitorElement * mePixClustersLad_py_
virtual const PixelTopology & specificTopology() const
Returns a reference to the pixel proxy topology.
int layerName() const
layer id
T const * product() const
MonitorElement * meYDisk_
MonitorElement * meSizeRing_
MonitorElement * mePixClusters_px_
const_iterator find(id_type i, bool update=false) const
MonitorElement * meSizeXPhi_
MonitorElement * meSizeXRing_
MonitorElement * mePixClusters_
MonitorElement * meSizeXBlade_
int pannelName() const
pannel id
MonitorElement * meMaxRow_
MonitorElement * meChargeLad_
MonitorElement * meSizeLad_
MonitorElement * meSizeBlade_
MonitorElement * mePixClustersRing_py_
int diskName() const
disk id
MonitorElement * meSizeXDisk_
MonitorElement * meSizeYPhi_
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
MonitorElement * mePixClustersPhi_px_
MonitorElement * meXBlade_
MonitorElement * meChargeLay_
tuple size
Write out results.
const_iterator begin(bool update=false) const
MonitorElement * meNClusters_
virtual const TrackerGeomDet * idToDet(DetId) const
MonitorElement * meSizeLay_
MonitorElement * meMaxColDisk_