50 src_( conf_.getParameter<edm::
InputTag>(
"src" ) ),
51 saveFile( conf_.getUntrackedParameter<bool>(
"saveFile",
false) ),
52 isPIB( conf_.getUntrackedParameter<bool>(
"isPIB",
false) ),
53 slowDown( conf_.getUntrackedParameter<bool>(
"slowDown",
false) ),
54 modOn( conf_.getUntrackedParameter<bool>(
"modOn",
true) ),
55 twoDimOn( conf_.getUntrackedParameter<bool>(
"twoDimOn",
true) ),
56 twoDimModOn( conf_.getUntrackedParameter<bool>(
"twoDimModOn",
true) ),
57 twoDimOnlyLayDisk( conf_.getUntrackedParameter<bool>(
"twoDimOnlyLayDisk",
false) ),
58 hiRes( conf_.getUntrackedParameter<bool>(
"hiRes",
false) ),
59 reducedSet( conf_.getUntrackedParameter<bool>(
"reducedSet",
false) ),
60 ladOn( conf_.getUntrackedParameter<bool>(
"ladOn",
false) ),
61 layOn( conf_.getUntrackedParameter<bool>(
"layOn",
false) ),
62 phiOn( conf_.getUntrackedParameter<bool>(
"phiOn",
false) ),
63 ringOn( conf_.getUntrackedParameter<bool>(
"ringOn",
false) ),
64 bladeOn( conf_.getUntrackedParameter<bool>(
"bladeOn",
false) ),
65 diskOn( conf_.getUntrackedParameter<bool>(
"diskOn",
false) ),
66 bigEventSize( conf_.getUntrackedParameter<int>(
"bigEventSize",1000) ),
67 isUpgrade( conf_.getUntrackedParameter<bool>(
"isUpgrade",
false) ),
82 if (
isUpgrade) { nTOTmodules=1856; }
else { nTOTmodules=1440; }
83 while(!infile.eof()&&nModsInFile<nTOTmodules) {
89 LogInfo (
"PixelDQM") <<
"SiPixelDigiSource::SiPixelDigiSource: Got DQM BackEnd interface"<<endl;
97 LogInfo (
"PixelDQM") <<
"SiPixelDigiSource::~SiPixelDigiSource: Destructor"<<endl;
129 for(
int i=0;
i!=40;
i++){
131 float averageOcc = 0.;
133 if(averageBPIXFed>0.) averageOcc =
nDigisPerFed[
i]/averageBPIXFed;
135 if(averageFPIXFed>0.) averageOcc =
nDigisPerFed[
i]/averageFPIXFed;
141 if (thisls % 10 == 0)
153 LogInfo (
"PixelDQM") <<
" SiPixelDigiSource::beginJob - Initialisation ... " << std::endl;
159 LogInfo (
"PixelDQM") <<
"2DIM IS " <<
twoDimOn <<
" and set to high resolution? " <<
hiRes <<
"\n";
245 int nEventDigis = 0;
int nActiveModules = 0;
247 std::map<uint32_t,SiPixelDigiModule*>::iterator struct_iter;
251 int NzeroROCs[2] = {0,-672};
252 int NloEffROCs[2] = {0,-672};
254 int numberOfDigisMod = (*struct_iter).second->fill(*input, iSetup,
265 if(numberOfDigisMod>0){
266 nEventDigis = nEventDigis + numberOfDigisMod;
279 for(
int i=0;
i!=nBPiXmodules; ++
i){
280 if((*struct_iter).first ==
I_detId[
i]){
282 int index1 = 0;
int index2 = 0;
301 if(panel==1){
if(module==1)
nDM1P1M1+=numberOfDigisMod;
302 else if(module==2)
nDM1P1M2+=numberOfDigisMod;
303 else if(module==3)
nDM1P1M3+=numberOfDigisMod;
304 else if(module==4)
nDM1P1M4+=numberOfDigisMod;}
305 else if(panel==2){
if(module==1)
nDM1P2M1+=numberOfDigisMod;
306 else if(module==2)
nDM1P2M2+=numberOfDigisMod;
307 else if(module==3)
nDM1P2M3+=numberOfDigisMod; }
308 if(blade<13 && blade>0 && (panel==1 || panel==2)) iter = i+2*(blade-1)+(panel-1);
311 if(panel==1){
if(module==1)
nDM2P1M1+=numberOfDigisMod;
312 else if(module==2)
nDM2P1M2+=numberOfDigisMod;
313 else if(module==3)
nDM2P1M3+=numberOfDigisMod;
314 else if(module==4)
nDM2P1M4+=numberOfDigisMod;}
315 else if(panel==2){
if(module==1)
nDM2P2M1+=numberOfDigisMod;
316 else if(module==2)
nDM2P2M2+=numberOfDigisMod;
317 else if(module==3)
nDM2P2M3+=numberOfDigisMod; }
318 if(blade<13 && blade>0 && (panel==1 || panel==2)) iter = i+2*(blade-1)+(panel-1);
323 if(panel==1){
if(module==1)
nDM1P1M1+=numberOfDigisMod;
324 else if(module==2)
nDM1P1M2+=numberOfDigisMod;
325 else if(module==3)
nDM1P1M3+=numberOfDigisMod;
326 else if(module==4)
nDM1P1M4+=numberOfDigisMod;}
327 else if(panel==2){
if(module==1)
nDM1P2M1+=numberOfDigisMod;
328 else if(module==2)
nDM1P2M2+=numberOfDigisMod;
329 else if(module==3)
nDM1P2M3+=numberOfDigisMod; }
330 if(blade<13 && blade>0 && (panel==1 || panel==2)) iter = i+2*(blade-1)+(panel-1);
333 if(panel==1){
if(module==1)
nDM2P1M1+=numberOfDigisMod;
334 else if(module==2)
nDM2P1M2+=numberOfDigisMod;
335 else if(module==3)
nDM2P1M3+=numberOfDigisMod;
336 else if(module==4)
nDM2P1M4+=numberOfDigisMod;}
337 else if(panel==2){
if(module==1)
nDM2P2M1+=numberOfDigisMod;
338 else if(module==2)
nDM2P2M2+=numberOfDigisMod;
339 else if(module==3)
nDM2P2M3+=numberOfDigisMod; }
340 if(blade<13 && blade>0 && (panel==1 || panel==2)) iter = i+2*(blade-1)+(panel-1);
345 if(panel==1){
if(module==1)
nDP1P1M1+=numberOfDigisMod;
346 else if(module==2)
nDP1P1M2+=numberOfDigisMod;
347 else if(module==3)
nDP1P1M3+=numberOfDigisMod;
348 else if(module==4)
nDP1P1M4+=numberOfDigisMod;}
349 else if(panel==2){
if(module==1)
nDP1P2M1+=numberOfDigisMod;
350 else if(module==2)
nDP1P2M2+=numberOfDigisMod;
351 else if(module==3)
nDP1P2M3+=numberOfDigisMod; }
352 if(blade<13 && blade>0 && (panel==1 || panel==2)) iter = i+2*(blade-1)+(panel-1);
355 if(panel==1){
if(module==1)
nDP2P1M1+=numberOfDigisMod;
356 else if(module==2)
nDP2P1M2+=numberOfDigisMod;
357 else if(module==3)
nDP2P1M3+=numberOfDigisMod;
358 else if(module==4)
nDP2P1M4+=numberOfDigisMod;}
359 else if(panel==2){
if(module==1)
nDP2P2M1+=numberOfDigisMod;
360 else if(module==2)
nDP2P2M2+=numberOfDigisMod;
361 else if(module==3)
nDP2P2M3+=numberOfDigisMod; }
362 if(blade<13 && blade>0 && (panel==1 || panel==2)) iter = i+2*(blade-1)+(panel-1);
367 if(panel==1){
if(module==1)
nDP1P1M1+=numberOfDigisMod;
368 else if(module==2)
nDP1P1M2+=numberOfDigisMod;
369 else if(module==3)
nDP1P1M3+=numberOfDigisMod;
370 else if(module==4)
nDP1P1M4+=numberOfDigisMod;}
371 else if(panel==2){
if(module==1)
nDP1P2M1+=numberOfDigisMod;
372 else if(module==2)
nDP1P2M2+=numberOfDigisMod;
373 else if(module==3)
nDP1P2M3+=numberOfDigisMod; }
374 if(blade<13 && blade>0 && (panel==1 || panel==2)) iter = i+2*(blade-1)+(panel-1);
377 if(panel==1){
if(module==1)
nDP2P1M1+=numberOfDigisMod;
378 else if(module==2)
nDP2P1M2+=numberOfDigisMod;
379 else if(module==3)
nDP2P1M3+=numberOfDigisMod;
380 else if(module==4)
nDP2P1M4+=numberOfDigisMod;}
381 else if(panel==2){
if(module==1)
nDP2P2M1+=numberOfDigisMod;
382 else if(module==2)
nDP2P2M2+=numberOfDigisMod;
383 else if(module==3)
nDP2P2M3+=numberOfDigisMod; }
384 if(blade<13 && blade>0 && (panel==1 || panel==2)) iter = i+2*(blade-1)+(panel-1);
389 for(
int i=nBPiXmodules; i!=nTOTmodules; i++){
390 if((*struct_iter).first ==
I_detId[i]){
408 if(panel==1){
if(module==1)
nDM1P1M1+=numberOfDigisMod; }
409 else if(panel==2){
if(module==1)
nDM1P2M1+=numberOfDigisMod; }
410 if(blade<12 && blade>0 && (panel==1 || panel==2)) iter = i+2*(blade-1)+(panel-1);
413 if(panel==1){
if(module==1)
nDM2P1M1+=numberOfDigisMod; }
414 else if(panel==2){
if(module==1)
nDM2P2M1+=numberOfDigisMod; }
415 if(blade<12 && blade>0 && (panel==1 || panel==2)) iter = i+2*(blade-1)+(panel-1);
418 if(panel==1){
if(module==1)
nDM3P1M1+=numberOfDigisMod; }
419 else if(panel==2){
if(module==1)
nDM3P2M1+=numberOfDigisMod; }
420 if(blade<12 && blade>0 && (panel==1 || panel==2)) iter = i+2*(blade-1)+(panel-1);
425 if(panel==1){
if(module==1)
nDM1P1M1+=numberOfDigisMod; }
426 else if(panel==2){
if(module==1)
nDM1P2M1+=numberOfDigisMod; }
427 if(blade<18 && blade>0 && (panel==1 || panel==2)) iter = i+2*(blade-1)+(panel-1);
430 if(panel==1){
if(module==1)
nDM2P1M1+=numberOfDigisMod; }
431 else if(panel==2){
if(module==1)
nDM2P2M1+=numberOfDigisMod; }
432 if(blade<18 && blade>0 && (panel==1 || panel==2)) iter = i+2*(blade-1)+(panel-1);
435 if(panel==1){
if(module==1)
nDM3P1M1+=numberOfDigisMod; }
436 else if(panel==2){
if(module==1)
nDM3P2M1+=numberOfDigisMod; }
437 if(blade<18 && blade>0 && (panel==1 || panel==2)) iter = i+2*(blade-1)+(panel-1);
442 if(panel==1){
if(module==1)
nDP1P1M1+=numberOfDigisMod; }
443 else if(panel==2){
if(module==1)
nDP1P2M1+=numberOfDigisMod; }
444 if(blade<12 && blade>0 && (panel==1 || panel==2)) iter = i+2*(blade-1)+(panel-1);
447 if(panel==1){
if(module==1)
nDP2P1M1+=numberOfDigisMod; }
448 else if(panel==2){
if(module==1)
nDP2P2M1+=numberOfDigisMod; }
449 if(blade<12 && blade>0 && (panel==1 || panel==2)) iter = i+2*(blade-1)+(panel-1);
452 if(panel==1){
if(module==1)
nDP3P1M1+=numberOfDigisMod; }
453 else if(panel==2){
if(module==1)
nDP3P2M1+=numberOfDigisMod; }
454 if(blade<12 && blade>0 && (panel==1 || panel==2)) iter = i+2*(blade-1)+(panel-1);
459 if(panel==1){
if(module==1)
nDP1P1M1+=numberOfDigisMod; }
460 else if(panel==2){
if(module==1)
nDP1P2M1+=numberOfDigisMod; }
461 if(blade<18 && blade>0 && (panel==1 || panel==2)) iter = i+2*(blade-1)+(panel-1);
464 if(panel==1){
if(module==1)
nDP2P1M1+=numberOfDigisMod; }
465 else if(panel==2){
if(module==1)
nDP2P2M1+=numberOfDigisMod; }
466 if(blade<18 && blade>0 && (panel==1 || panel==2)) iter = i+2*(blade-1)+(panel-1);
469 if(panel==1){
if(module==1)
nDP3P1M1+=numberOfDigisMod; }
470 else if(panel==2){
if(module==1)
nDP3P2M1+=numberOfDigisMod; }
471 if(blade<18 && blade>0 && (panel==1 || panel==2)) iter = i+2*(blade-1)+(panel-1);
475 for(
int i=nBPiXmodules; i!=nTOTmodules; i++){
476 if((*struct_iter).first ==
I_detId[i]){
485 std::pair<int,int> tempPair = (*struct_iter).second->getZeroLoEffROCs();
487 NzeroROCs[0] += tempPair.first;
488 NloEffROCs[0] += tempPair.second;
491 NzeroROCs[1] += tempPair.first;
492 NloEffROCs[1] += tempPair.second;
497 if (lumiSection%10> 2){
498 for (
int i =0;
i < 2; ++
i) NloEffROCs[
i] = NloEffROCs[
i] - NzeroROCs[
i];
565 if(nActiveModules>=4){
583 LogInfo (
"PixelDQM") <<
" SiPixelDigiSource::buildStructure" ;
587 LogVerbatim (
"PixelDQM") <<
" *** Geometry node for TrackerGeom is "<<&(*pDD)<<std::endl;
588 LogVerbatim (
"PixelDQM") <<
" *** I have " << pDD->dets().size() <<
" detectors"<<std::endl;
589 LogVerbatim (
"PixelDQM") <<
" *** I have " << pDD->detTypes().size() <<
" types"<<std::endl;
591 for(TrackerGeometry::DetContainer::const_iterator it = pDD->dets().begin(); it != pDD->dets().end(); it++){
593 if(dynamic_cast<PixelGeomDetUnit const *>((*it))!=0){
595 DetId detId = (*it)->geographicalId();
596 const GeomDetUnit * geoUnit = pDD->idToDetUnit( detId );
603 LogDebug (
"PixelDQM") <<
" ---> Adding Barrel Module " << detId.
rawId() << endl;
604 uint32_t
id = detId();
611 LogDebug (
"PixelDQM") <<
" ---> Adding Endcap Module " << detId.
rawId() << endl;
612 uint32_t
id = detId();
622 char sside[80]; sprintf(sside,
"HalfCylinder_%i",side);
623 char sdisk[80]; sprintf(sdisk,
"Disk_%i",disk);
624 char sblade[80]; sprintf(sblade,
"Blade_%02i",blade);
625 char spanel[80]; sprintf(spanel,
"Panel_%i",panel);
626 char smodule[80];sprintf(smodule,
"Module_%i",module);
629 bool mask = side_str.find(
"HalfCylinder_1")!=string::npos||
630 side_str.find(
"HalfCylinder_2")!=string::npos||
631 side_str.find(
"HalfCylinder_4")!=string::npos||
632 disk_str.find(
"Disk_2")!=string::npos;
635 if(
isPIB && mask)
continue;
639 LogDebug (
"PixelDQM") <<
" ---> Adding Endcap Module " << detId.
rawId() << endl;
640 uint32_t
id = detId();
650 char sside[80]; sprintf(sside,
"HalfCylinder_%i",side);
651 char sdisk[80]; sprintf(sdisk,
"Disk_%i",disk);
652 char sblade[80]; sprintf(sblade,
"Blade_%02i",blade);
653 char spanel[80]; sprintf(spanel,
"Panel_%i",panel);
654 char smodule[80];sprintf(smodule,
"Module_%i",module);
657 bool mask = side_str.find(
"HalfCylinder_1")!=string::npos||
658 side_str.find(
"HalfCylinder_2")!=string::npos||
659 side_str.find(
"HalfCylinder_4")!=string::npos||
660 disk_str.find(
"Disk_2")!=string::npos;
663 if(
isPIB && mask)
continue;
679 char title[80]; sprintf(title,
"Rate of events with >%i digis;LumiSection;Rate [Hz]",
bigEventSize);
681 char title1[80]; sprintf(title1,
"Pixel events vs. BX;BX;# events");
683 char title2[80]; sprintf(title2,
"Rate of Pixel events;LumiSection;Rate [Hz]");
685 char title3[80]; sprintf(title3,
"Number of Zero-Occupancy Barrel ROCs;LumiSection;N_{ZERO-OCCUPANCY} Barrel ROCs");
687 char title4[80]; sprintf(title4,
"Number of Low-Efficiency Barrel ROCs;LumiSection;N_{LO EFF} Barrel ROCs");
689 char title5[80]; sprintf(title5,
"Number of Zero-Occupancy Endcap ROCs;LumiSection;N_{ZERO-OCCUPANCY} Endcap ROCs");
691 char title6[80]; sprintf(title6,
"Number of Low-Efficiency Endcap ROCs;LumiSection;N_{LO EFF} Endcap ROCs");
693 char title7[80]; sprintf(title7,
"Average digi occupancy per FED;FED;NDigis/<NDigis>");
697 char title4[80]; sprintf(title4,
"FED Digi Occupancy (NDigis/<NDigis>) vs LumiSections;Lumi Section;FED");
700 std::map<uint32_t,SiPixelDigiModule*>::iterator struct_iter;
712 <<
"[SiPixelDigiSource::bookMEs] Creation of DQM folder failed";
719 LogDebug (
"PixelDQM") <<
"PROBLEM WITH LADDER-FOLDER\n";
727 LogDebug (
"PixelDQM") <<
"PROBLEM WITH LAYER-FOLDER\n";
735 LogDebug (
"PixelDQM") <<
"PROBLEM WITH PHI-FOLDER\n";
742 LogDebug (
"PixelDQM") <<
"PROBLEM WITH BLADE-FOLDER\n";
749 LogDebug (
"PixelDQM") <<
"PROBLEM WITH DISK-FOLDER\n";
756 LogDebug (
"PixelDQM") <<
"PROBLEM WITH RING-FOLDER\n";
765 std::stringstream ss1, ss2;
768 ss1.str(
std::string()); ss1 <<
"ALLMODS_ndigisCHAN_BarrelL" <<
i;
769 ss2.str(
std::string()); ss2 <<
"Number of Digis L" <<
i;
852 ss1.str(
std::string()); ss1 <<
"ALLMODS_ndigisCHAN_EndcapDp" <<
i;
853 ss2.str(
std::string()); ss2 <<
"Number of Digis Disk p" <<
i;
859 ss1.str(
std::string()); ss1 <<
"ALLMODS_ndigisCHAN_EndcapDm" <<
i;
860 ss2.str(
std::string()); ss2 <<
"Number of Digis Disk m" <<
i;
int plaquetteName() const
plaquetteId (in pannel)
LuminosityBlockID id() const
T getParameter(std::string const &) const
MonitorElement * meNDigisCHANBarrelCh19_
MonitorElement * meNDigisCHANBarrelCh35_
MonitorElement * meNDigisCHANBarrelCh2_
MonitorElement * meNDigisCHANBarrelCh30_
void setBinContent(int binx, double content)
set content of bin (1-D)
MonitorElement * meNDigisCHANBarrelCh31_
MonitorElement * meNDigisCHANBarrelCh8_
MonitorElement * bookProfile(Args &&...args)
MonitorElement * meNDigisCHANBarrelCh36_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
#define DEFINE_FWK_MODULE(type)
MonitorElement * loOccROCsEndcap
MonitorElement * meNDigisCHANBarrelCh17_
MonitorElement * meNDigisCHANBarrelCh25_
MonitorElement * meNDigisCHANBarrelCh6_
MonitorElement * meNDigisCOMBBarrel_
int bunchCrossing() const
std::string topFolderName_
edm::LuminosityBlockNumber_t luminosityBlock() const
MonitorElement * meNDigisCOMBEndcap_
MonitorElement * meNDigisCHANBarrelCh14_
MonitorElement * meNDigisCHANBarrelCh12_
MonitorElement * meNDigisCHANBarrelCh29_
static std::string const input
edm::EDGetTokenT< edm::DetSetVector< PixelDigi > > srcToken_
uint32_t rawId() const
get the raw id
MonitorElement * meNDigisCHANBarrelCh27_
int bladeName() const
blade id
virtual void bookMEs(DQMStore::IBooker &, const edm::EventSetup &iSetup)
SiPixelDigiSource(const edm::ParameterSet &conf)
MonitorElement * meNDigisCHANBarrelCh3_
MonitorElement * meNDigisCHANBarrelCh9_
MonitorElement * meNDigisCHANBarrelCh4_
MonitorElement * meNDigisCHANBarrelCh34_
MonitorElement * meNDigisCHANBarrelCh26_
bool setModuleFolder(const uint32_t &rawdetid=0, int type=0, bool isUpgrade=false)
Set folder name for a module or plaquette.
MonitorElement * noOccROCsBarrel
MonitorElement * meNDigisCHANEndcap_
MonitorElement * book1D(Args &&...args)
MonitorElement * meNDigisCHANBarrelCh28_
MonitorElement * noOccROCsEndcap
virtual void analyze(const edm::Event &, const edm::EventSetup &)
MonitorElement * meNDigisCHANBarrel_
std::vector< MonitorElement * > meNDigisCHANEndcapDms_
virtual void dqmBeginRun(const edm::Run &, edm::EventSetup const &)
MonitorElement * meNDigisCHANBarrelCh33_
MonitorElement * pixEventRate
MonitorElement * meNDigisCHANBarrelCh21_
unsigned int I_detId[1856]
int subdetId() const
get the contents of the subdetector field (not cast into any detector's numbering enum) ...
MonitorElement * meNDigisCHANBarrelCh32_
MonitorElement * meNDigisCHANBarrelCh10_
MonitorElement * loOccROCsBarrel
MonitorElement * meNDigisCHANBarrelCh1_
MonitorElement * averageDigiOccupancy
void setCurrentFolder(const std::string &fullpath)
MonitorElement * meNDigisCHANBarrelCh5_
MonitorElement * book2D(Args &&...args)
virtual void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
MonitorElement * meNDigisCHANBarrelCh7_
MonitorElement * meNDigisCHANBarrelCh15_
MonitorElement * meNDigisCHANBarrelCh13_
virtual const PixelTopology & specificTopology() const
Returns a reference to the pixel proxy topology.
MonitorElement * meNDigisCHANBarrelCh23_
int layerName() const
layer id
T const * product() const
LuminosityBlockNumber_t luminosityBlock() const
MonitorElement * meNDigisCHANBarrelCh18_
MonitorElement * bigEventRate
virtual void beginLuminosityBlock(edm::LuminosityBlock const &, edm::EventSetup const &) override
MonitorElement * meNDigisCHANBarrelCh20_
MonitorElement * avgfedDigiOccvsLumi
int pannelName() const
pannel id
MonitorElement * meNDigisCHANBarrelCh22_
int diskName() const
disk id
std::vector< MonitorElement * > meNDigisCHANEndcapDps_
MonitorElement * pixEvtsPerBX
virtual void endLuminosityBlock(edm::LuminosityBlock const &, edm::EventSetup const &) override
volatile std::atomic< bool > shutdown_flag false
std::vector< MonitorElement * > meNDigisCHANBarrelLs_
std::map< uint32_t, SiPixelDigiModule * > thePixelStructure
void setLumiFlag(void)
this ME is meant to be stored for each luminosity section
HalfCylinder halfCylinder() const
MonitorElement * meNDigisCHANBarrelCh16_
MonitorElement * meNDigisCHANBarrelCh11_
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
virtual void buildStructure(edm::EventSetup const &)
MonitorElement * meNDigisCHANBarrelCh24_