74 bool isHalfModule =
false;
117 meX_ = iBooker.
book1D(hid,
"Cluster barycenter X (row #)",200,0.,200.);
121 meY_ = iBooker.
book1D(hid,
"Cluster barycenter Y (column #)",500,0.,500.);
125 meSizeX_ = iBooker.
book1D(hid,
"Cluster x-width (rows)",10,0.,10.);
129 meSizeY_ = iBooker.
book1D(hid,
"Cluster y-width (columns)",15,0.,15.);
147 delete theHistogramId;
151 if(barrel && type==7){
152 hid = src.
label() +
"_Barrel";
153 meSizeYvsEtaBarrel_= iBooker.
book2D(
"sizeYvsEta_" + hid,
"Cluster size along beamline vs. Cluster position #eta",60,-3.,3.,40,0.,40.);
157 if(type==1 && barrel){
160 char sladder[80]; sprintf(sladder,
"Ladder_%02i",DBladder);
161 hid = src.
label() +
"_" + sladder;
162 if(isHalfModule) hid +=
"H";
187 meXLad_ = iBooker.
book1D(
"x_" + hid,
"Cluster barycenter X (row #)",200,0.,200.);
190 meYLad_ = iBooker.
book1D(
"y_" + hid,
"Cluster barycenter Y (column #)",500,0.,500.);
196 meSizeYLad_ = iBooker.
book1D(
"sizeY_" + hid,
"Cluster y-width (columns)",15,0.,15.);
213 if(type==2 && barrel){
217 char slayer[80]; sprintf(slayer,
"Layer_%i",DBlayer);
218 hid = src.
label() +
"_" + slayer;
242 meXLay_ = iBooker.
book1D(
"x_" + hid,
"Cluster barycenter X (row #)",200,0.,200.);
245 meYLay_ = iBooker.
book1D(
"y_" + hid,
"Cluster barycenter Y (column #)",500,0.,500.);
251 meSizeYLay_ = iBooker.
book1D(
"sizeY_" + hid,
"Cluster y-width (columns)",15,0.,15.);
275 if(type==3 && barrel){
278 char smodule[80]; sprintf(smodule,
"Ring_%i",DBmodule);
279 hid = src.
label() +
"_" + smodule;
303 meXPhi_ = iBooker.
book1D(
"x_" + hid,
"Cluster barycenter X (row #)",200,0.,200.);
306 meYPhi_ = iBooker.
book1D(
"y_" + hid,
"Cluster barycenter Y (column #)",500,0.,500.);
312 meSizeYPhi_ = iBooker.
book1D(
"sizeY_" + hid,
"Cluster y-width (columns)",15,0.,15.);
337 if(type==4 && endcap){
341 char sblade[80]; sprintf(sblade,
"Blade_%02i",blade);
342 hid = src.
label() +
"_" + sblade;
366 meXBlade_ = iBooker.
book1D(
"x_" + hid,
"Cluster barycenter X (row #)",200,0.,200.);
369 meYBlade_ = iBooker.
book1D(
"y_" + hid,
"Cluster barycenter Y (column #)",500,0.,500.);
379 if(type==5 && endcap){
383 char sdisk[80]; sprintf(sdisk,
"Disk_%i",disk);
384 hid = src.
label() +
"_" + sdisk;
408 meXDisk_ = iBooker.
book1D(
"x_" + hid,
"Cluster barycenter X (row #)",200,0.,200.);
411 meYDisk_ = iBooker.
book1D(
"y_" + hid,
"Cluster barycenter Y (column #)",500,0.,500.);
422 if(type==6 && endcap){
428 char slab[80]; sprintf(slab,
"Panel_%i_Ring_%i",panel, module);
429 hid = src.
label() +
"_" + slab;
453 meXRing_ = iBooker.
book1D(
"x_" + hid,
"Cluster barycenter X (row #)",200,0.,200.);
456 meYRing_ = iBooker.
book1D(
"y_" + hid,
"Cluster barycenter Y (column #)",500,0.,500.);
483 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) {
489 unsigned int numberOfClusters = 0;
490 unsigned int numberOfFpixClusters = 0;
492 if( isearch != input.
end() ) {
497 for(di = isearch->
begin(); di != isearch->
end(); di++) {
499 if(endcap) numberOfFpixClusters++;
500 float charge = 0.001*(di->charge());
504 int sizeX = di->sizeX();
505 int sizeY = di->sizeY();
506 int minPixelRow = di->minPixelRow();
507 int maxPixelRow = di->maxPixelRow();
508 int minPixelCol = di->minPixelCol();
509 int maxPixelCol = di->maxPixelCol();
521 if (!(DBlayer > layers.size()) && (layers[DBlayer-1])) layers[DBlayer-1]->
Fill(clustgp.
z(),clustgp.
phi());
525 if (!(DBdisk > diskspz.size()) && (diskspz[DBdisk-1])) diskspz[DBdisk-1]->
Fill(clustgp.
x(),clustgp.
y());
527 if (!(DBdisk > disksmz.size()) && (disksmz[DBdisk-1])) disksmz[DBdisk-1]->
Fill(clustgp.
x(),clustgp.
y());
547 if(barrel && smileyon){
611 if(bladeon && endcap){
626 if(diskon && endcap){
642 if(ringon && endcap){
677 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_