SiPixelDigiModule Class Reference

#include <SiPixelDigiModule.h>

Public Types
typedef edm::DetSet< PixelDigi > ::const_iterator	DigiIterator

Public Member Functions
void	book (const edm::ParameterSet &iConfig, DQMStore::IBooker &iBooker, int type=0, bool twoD=true, bool hiRes=false, bool reducedSet=false, bool additInfo=false, bool isUpgrade=false)
	Book histograms. More...
int	fill (const edm::DetSetVector< PixelDigi > &input, MonitorElement *combBarrel, MonitorElement *chanBarrel, MonitorElement *chanBarrelL1, MonitorElement *chanBarrelL2, MonitorElement *chanBarrelL3, MonitorElement *chanBarrelL4, MonitorElement *combEndcap, const bool modon, const bool ladon, const bool layon, const bool phion, const bool bladeon, const bool diskon, const bool ringon, const bool twoD, const bool reducedSet, const bool twoDimModOn, const bool twoDimOnlyLayDisk, int &nDigisA, int &nDigisB, bool isUpgrade)
	Fill histograms. More...
std::pair< int, int >	getZeroLoEffROCs ()
void	resetRocMap ()
	SiPixelDigiModule ()
	Default constructor. More...
	SiPixelDigiModule (const uint32_t &id)
	Constructor with raw DetId. More...
	SiPixelDigiModule (const uint32_t &id, const int &ncols, const int &nrows)
	Constructor with raw DetId and sensor size. More...
	~SiPixelDigiModule ()
	Destructor. More...

Private Attributes
uint32_t	id_
MonitorElement *	meADC_
MonitorElement *	meADCBlade_
MonitorElement *	meADCDisk_
MonitorElement *	meADCLad_
MonitorElement *	meADCLay_
MonitorElement *	meADCPhi_
MonitorElement *	meADCRing_
MonitorElement *	meNDigis_
MonitorElement *	meNDigisBlade_
MonitorElement *	meNDigisDisk_
MonitorElement *	meNDigisLad_
MonitorElement *	meNDigisLay_
MonitorElement *	meNDigisPhi_
MonitorElement *	meNDigisRing_
MonitorElement *	mePixDigis_
MonitorElement *	mePixDigis_px_
MonitorElement *	mePixDigis_py_
MonitorElement *	mePixDigisDisk_
MonitorElement *	mePixDigisLad_
MonitorElement *	mePixDigisLad_px_
MonitorElement *	mePixDigisLad_py_
MonitorElement *	mePixDigisLay_
MonitorElement *	mePixDigisLay_px_
MonitorElement *	mePixDigisLay_py_
MonitorElement *	mePixDigisPhi_
MonitorElement *	mePixDigisPhi_px_
MonitorElement *	mePixDigisPhi_py_
MonitorElement *	mePixDigisRing_
MonitorElement *	mePixDigisRing_px_
MonitorElement *	mePixDigisRing_py_
MonitorElement *	mePixRocsDisk_ = NULL
MonitorElement *	mePixRocsLay_ = NULL
MonitorElement *	meZeroOccRocsDisk_ = NULL
MonitorElement *	meZeroOccRocsLay_ = NULL
int	ncols_
int	nrows_

Detailed Description

Definition at line 34 of file SiPixelDigiModule.h.

Member Typedef Documentation

typedef edm::DetSet<PixelDigi>::const_iterator SiPixelDigiModule::DigiIterator

Definition at line 47 of file SiPixelDigiModule.h.

Constructor & Destructor Documentation

SiPixelDigiModule::SiPixelDigiModule

(

)

Default constructor.

Framework.

Definition at line 26 of file SiPixelDigiModule.cc.

                                     : id_(0),
                     ncols_(416),
                     nrows_(160) 
{
}

SiPixelDigiModule::SiPixelDigiModule

(

const uint32_t &

id

)

Constructor with raw DetId.

Definition at line 32 of file SiPixelDigiModule.cc.

                                                       : 
  id_(id),
  ncols_(416),
  nrows_(160)
{ 
}

SiPixelDigiModule::SiPixelDigiModule	(	const uint32_t &	id,
		const int &	ncols,
		const int &	nrows
	)

Constructor with raw DetId and sensor size.

Definition at line 39 of file SiPixelDigiModule.cc.

                                                                                           : 
  id_(id),
  ncols_(ncols),
  nrows_(nrows)
{ 
}

SiPixelDigiModule::~SiPixelDigiModule

(

)

Destructor.

Definition at line 48 of file SiPixelDigiModule.cc.

{}

Member Function Documentation

void SiPixelDigiModule::book	(	const edm::ParameterSet &	iConfig,
		DQMStore::IBooker &	iBooker,
		int	type = 0,
		bool	twoD = true,
		bool	hiRes = false,
		bool	reducedSet = false,
		bool	additInfo = false,
		bool	isUpgrade = false
	)

Book histograms.

Definition at line 52 of file SiPixelDigiModule.cc.

References Reference_intrackfit_cff::barrel, PixelEndcapNameUpgrade::bladeName(), PixelEndcapName::bladeName(), DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), PixelEndcapNameUpgrade::diskName(), PixelEndcapName::diskName(), Reference_intrackfit_cff::endcap, edm::ParameterSet::getParameter(), id_, PixelBarrelNameUpgrade::isHalfModule(), PixelBarrelName::isHalfModule(), edm::InputTag::label(), PixelBarrelNameUpgrade::ladderName(), PixelBarrelName::ladderName(), PixelBarrelNameUpgrade::layerName(), PixelBarrelName::layerName(), meADC_, meADCBlade_, meADCDisk_, meADCLad_, meADCLay_, meADCPhi_, meADCRing_, meNDigis_, meNDigisBlade_, meNDigisDisk_, meNDigisLad_, meNDigisLay_, meNDigisPhi_, meNDigisRing_, mePixDigis_, mePixDigis_px_, mePixDigis_py_, mePixDigisDisk_, mePixDigisLad_, mePixDigisLad_px_, mePixDigisLad_py_, mePixDigisLay_, mePixDigisLay_px_, mePixDigisLay_py_, mePixDigisPhi_, mePixDigisPhi_px_, mePixDigisPhi_py_, mePixDigisRing_, mePixDigisRing_px_, mePixDigisRing_py_, mePixRocsDisk_, mePixRocsLay_, meZeroOccRocsDisk_, meZeroOccRocsLay_, python.rootplot.argparse::module, PixelBarrelNameUpgrade::moduleName(), PixelBarrelName::moduleName(), ncols_, nrows_, PixelEndcapNameUpgrade::pannelName(), PixelEndcapName::pannelName(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, PixelEndcapNameUpgrade::plaquetteName(), PixelEndcapName::plaquetteName(), MonitorElement::setAxisTitle(), SiPixelHistogramId::setHistoId(), alcazmumu_cfi::src, AlCaHLTBitMon_QueryRunRegistry::string, and DetId::subdetId().

                                                                                                                                                                        {

  //isUpgrade = iConfig.getUntrackedParameter<bool>("isUpgrade");
    
  bool barrel = DetId(id_).subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel);
  bool endcap = DetId(id_).subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap);
  bool isHalfModule = false;
  if(barrel){
    if (!isUpgrade) {
    isHalfModule = PixelBarrelName(DetId(id_)).isHalfModule(); 
    } else if (isUpgrade) {
      isHalfModule = PixelBarrelNameUpgrade(DetId(id_)).isHalfModule(); 
    }
  }

  std::string hid;
  // Get collection name and instantiate Histo Id builder
  edm::InputTag src = iConfig.getParameter<edm::InputTag>( "src" );
  
  int nbinx=ncols_/2, nbiny=nrows_/2;
  std::string twodtitle           = "Number of Digis (1bin=four pixels)"; 
  std::string pxtitle             = "Number of Digis (1bin=two columns)";
  std::string pytitle             = "Number of Digis (1bin=two rows)";
  std::string twodroctitle        = "ROC Occupancy (1bin=one ROC)";
  std::string twodzeroOccroctitle = "Zero Occupancy ROC Map (1bin=one ROC) for ";
  if(hiRes){
    nbinx = ncols_;
    nbiny = nrows_;
    twodtitle    = "Number of Digis (1bin=one pixel)";
    pxtitle = "Number of Digis (1bin=one column)";
    pytitle = "Number of Digis (1bin=one row)";
  }
  if(type==0){
    SiPixelHistogramId* theHistogramId = new SiPixelHistogramId( src.label() );
    // Number of digis
    hid = theHistogramId->setHistoId("ndigis",id_);
    meNDigis_ = iBooker.book1D(hid,"Number of Digis",25,0.,25.);
    meNDigis_->setAxisTitle("Number of digis",1);
    // Charge in ADC counts
    hid = theHistogramId->setHistoId("adc",id_);
    meADC_ = iBooker.book1D(hid,"Digi charge",128,0.,256.);
    meADC_->setAxisTitle("ADC counts",1);
    if(!reducedSet)
    {
    if(twoD){
      if(additInfo){
    // 2D hit map
    hid = theHistogramId->setHistoId("hitmap",id_);
    mePixDigis_ = iBooker.book2D(hid,twodtitle,nbinx,0.,float(ncols_),nbiny,0.,float(nrows_));
    mePixDigis_->setAxisTitle("Columns",1);
    mePixDigis_->setAxisTitle("Rows",2);
    //std::cout << "During booking: type is "<< type << ", ID is "<< id_ << ", pwd for booking is " << theDMBE->pwd() << ", Plot name: " << hid << std::endl;
      }
    }
    else{
      // projections of 2D hit map
      hid = theHistogramId->setHistoId("hitmap",id_);
      mePixDigis_px_ = iBooker.book1D(hid+"_px",pxtitle,nbinx,0.,float(ncols_));
      mePixDigis_py_ = iBooker.book1D(hid+"_py",pytitle,nbiny,0.,float(nrows_));
      mePixDigis_px_->setAxisTitle("Columns",1);
      mePixDigis_py_->setAxisTitle("Rows",1);
    }
    }
    delete theHistogramId;

  }
  
  if(type==1 && barrel){
    uint32_t DBladder;
    if (!isUpgrade) { DBladder = PixelBarrelName(DetId(id_)).ladderName();}
    else { DBladder = PixelBarrelNameUpgrade(DetId(id_)).ladderName();}
    char sladder[80]; sprintf(sladder,"Ladder_%02i",DBladder);
    hid = src.label() + "_" + sladder;
    if(isHalfModule) hid += "H";
    else hid += "F";
    // Number of digis
    meNDigisLad_ = iBooker.book1D("ndigis_"+hid,"Number of Digis",25,0.,25.);
    meNDigisLad_->setAxisTitle("Number of digis",1);
    // Charge in ADC counts
    meADCLad_ = iBooker.book1D("adc_" + hid,"Digi charge",128,0.,256.);
    meADCLad_->setAxisTitle("ADC counts",1);
    if(!reducedSet)
    {
    if(twoD){
      // 2D hit map
      mePixDigisLad_ = iBooker.book2D("hitmap_"+hid,twodtitle,nbinx,0.,float(ncols_),nbiny,0.,float(nrows_));
      mePixDigisLad_->setAxisTitle("Columns",1);
      mePixDigisLad_->setAxisTitle("Rows",2);
      //std::cout << "During booking: type is "<< type << ", ID is "<< id_ << ", pwd for booking is " << theDMBE->pwd() << ", Plot name: " << hid << std::endl;
    }
    else{
      // projections of 2D hit map
      mePixDigisLad_px_ = iBooker.book1D("hitmap_"+hid+"_px",pxtitle,nbinx,0.,float(ncols_));
      mePixDigisLad_py_ = iBooker.book1D("hitmap_"+hid+"_py",pytitle,nbiny,0.,float(nrows_));
      mePixDigisLad_px_->setAxisTitle("Columns",1);
      mePixDigisLad_py_->setAxisTitle("Rows",1);    
    }
    }
  }
  if(type==2 && barrel){
    uint32_t DBlayer;
    if (!isUpgrade) { DBlayer = PixelBarrelName(DetId(id_)).layerName(); }
    else { DBlayer = PixelBarrelNameUpgrade(DetId(id_)).layerName(); }
    char slayer[80]; sprintf(slayer,"Layer_%i",DBlayer);
    hid = src.label() + "_" + slayer;
    if(!additInfo){
      // Number of digis
      meNDigisLay_ = iBooker.book1D("ndigis_"+hid,"Number of Digis",25,0.,25.);
      meNDigisLay_->setAxisTitle("Number of digis",1);
    // Charge in ADC counts
      meADCLay_ = iBooker.book1D("adc_" + hid,"Digi charge",128,0.,256.);
      meADCLay_->setAxisTitle("ADC counts",1);
    }
    if(!reducedSet){
      if(twoD || additInfo){
    // 2D hit map
    if(isHalfModule){
      mePixDigisLay_ = iBooker.book2D("hitmap_"+hid,twodtitle,nbinx,0.,float(ncols_),2*nbiny,0.,float(2*nrows_));
    }
    else{
      mePixDigisLay_ = iBooker.book2D("hitmap_"+hid,twodtitle,nbinx,0.,float(ncols_),nbiny,0.,float(nrows_));

    }
    mePixDigisLay_->setAxisTitle("Columns",1);
    mePixDigisLay_->setAxisTitle("Rows",2);
    
    //std::cout << "During booking: type is "<< type << ", ID is "<< id_ << ", pwd for booking is " << theDMBE->pwd() << ", Plot name: " << hid << std::endl;
    int yROCbins[3] = {18,30,42};
    mePixRocsLay_ = iBooker.book2D("rocmap_"+hid,twodroctitle,32,0.,32.,yROCbins[DBlayer-1],1.5,1.5+float(yROCbins[DBlayer-1]/2));
    mePixRocsLay_->setAxisTitle("ROCs per Module",1);
    mePixRocsLay_->setAxisTitle("ROCs per 1/2 Ladder",2);
    meZeroOccRocsLay_ = iBooker.book2D("zeroOccROC_map",twodzeroOccroctitle+hid,32,0.,32.,yROCbins[DBlayer-1],1.5,1.5+float(yROCbins[DBlayer-1]/2));
    meZeroOccRocsLay_->setAxisTitle("ROCs per Module",1);
    meZeroOccRocsLay_->setAxisTitle("ROCs per 1/2 Ladder",2);
      }
      if(!twoD && !additInfo){
    // projections of 2D hit map
    mePixDigisLay_px_ = iBooker.book1D("hitmap_"+hid+"_px",pxtitle,nbinx,0.,float(ncols_));
    if(isHalfModule){
      mePixDigisLay_py_ = iBooker.book1D("hitmap_"+hid+"_py",pytitle,2*nbiny,0.,float(2*nrows_));
    }
    else{
      mePixDigisLay_py_ = iBooker.book1D("hitmap_"+hid+"_py",pytitle,nbiny,0.,float(nrows_));
    }
    mePixDigisLay_px_->setAxisTitle("Columns",1);
    mePixDigisLay_py_->setAxisTitle("Rows",1);
      }
    }
  }
  if(type==3 && barrel){
    uint32_t DBmodule;
    if (!isUpgrade) { DBmodule = PixelBarrelName(DetId(id_)).moduleName(); }
    else { DBmodule = PixelBarrelNameUpgrade(DetId(id_)).moduleName(); }
    char smodule[80]; sprintf(smodule,"Ring_%i",DBmodule);
    hid = src.label() + "_" + smodule;
    // Number of digis
    meNDigisPhi_ = iBooker.book1D("ndigis_"+hid,"Number of Digis",25,0.,25.);
    meNDigisPhi_->setAxisTitle("Number of digis",1);
    // Charge in ADC counts
    meADCPhi_ = iBooker.book1D("adc_" + hid,"Digi charge",128,0.,256.);
    meADCPhi_->setAxisTitle("ADC counts",1);
    if(!reducedSet)
      {
    if(twoD){
      
      // 2D hit map
      if(isHalfModule){
        mePixDigisPhi_ = iBooker.book2D("hitmap_"+hid,twodtitle,nbinx,0.,float(ncols_),2*nbiny,0.,float(2*nrows_));
      }
      else {
        mePixDigisPhi_ = iBooker.book2D("hitmap_"+hid,twodtitle,nbinx,0.,float(ncols_),nbiny,0.,float(nrows_));
      }
      mePixDigisPhi_->setAxisTitle("Columns",1);
      mePixDigisPhi_->setAxisTitle("Rows",2);
      //std::cout << "During booking: type is "<< type << ", ID is "<< id_ << ", pwd for booking is " << theDMBE->pwd() << ", Plot name: " << hid << std::endl;
    }
    else{
      // projections of 2D hit map
      mePixDigisPhi_px_ = iBooker.book1D("hitmap_"+hid+"_px",pxtitle,nbinx,0.,float(ncols_));
      if(isHalfModule){
        mePixDigisPhi_py_ = iBooker.book1D("hitmap_"+hid+"_py",pytitle,2*nbiny,0.,float(2*nrows_));
      }
      else{
        mePixDigisPhi_py_ = iBooker.book1D("hitmap_"+hid+"_py",pytitle,nbiny,0.,float(nrows_));
      }
      mePixDigisPhi_px_->setAxisTitle("Columns",1);
      mePixDigisPhi_py_->setAxisTitle("Rows",1);
    }
      }
  }
  if(type==4 && endcap){
    uint32_t blade;
    if (!isUpgrade) { blade= PixelEndcapName(DetId(id_)).bladeName(); }
    else { blade= PixelEndcapNameUpgrade(DetId(id_)).bladeName(); }
    
    char sblade[80]; sprintf(sblade, "Blade_%02i",blade);
    hid = src.label() + "_" + sblade;
    // Number of digis
    meNDigisBlade_ = iBooker.book1D("ndigis_"+hid,"Number of Digis",25,0.,25.);
    meNDigisBlade_->setAxisTitle("Number of digis",1);
    // Charge in ADC counts
    meADCBlade_ = iBooker.book1D("adc_" + hid,"Digi charge",128,0.,256.);
    meADCBlade_->setAxisTitle("ADC counts",1);
  }
  if(type==5 && endcap){
    uint32_t disk;
    if (!isUpgrade) { disk = PixelEndcapName(DetId(id_)).diskName(); }
    else { disk = PixelEndcapNameUpgrade(DetId(id_)).diskName(); }
    
    char sdisk[80]; sprintf(sdisk, "Disk_%i",disk);
    hid = src.label() + "_" + sdisk;
    if(!additInfo){
      // Number of digis
      meNDigisDisk_ = iBooker.book1D("ndigis_"+hid,"Number of Digis",25,0.,25.);
      meNDigisDisk_->setAxisTitle("Number of digis",1);
      // Charge in ADC counts
      meADCDisk_ = iBooker.book1D("adc_" + hid,"Digi charge",128,0.,256.);
      meADCDisk_->setAxisTitle("ADC counts",1);
    }
    if(additInfo){
      mePixDigisDisk_ = iBooker.book2D("hitmap_"+hid,twodtitle,260,0.,260.,160,0.,160.);
      mePixDigisDisk_->setAxisTitle("Columns",1);
      mePixDigisDisk_->setAxisTitle("Rows",2);
      //ROC information in disks
      mePixRocsDisk_  = iBooker.book2D("rocmap_"+hid,twodroctitle,26,0.,26.,24,1.,13.);
      mePixRocsDisk_ ->setAxisTitle("ROCs per Module (2 Panels)",1);
      mePixRocsDisk_ ->setAxisTitle("Blade Number",2);
      meZeroOccRocsDisk_  = iBooker.book2D("zeroOccROC_map",twodzeroOccroctitle+hid,26,0.,26.,24,1.,13.);
      meZeroOccRocsDisk_ ->setAxisTitle("Zero-Occupancy ROCs per Module (2 Panels)",1);
      meZeroOccRocsDisk_ ->setAxisTitle("Blade Number",2);
    }
  }
  if(type==6 && endcap){
    uint32_t panel;
    uint32_t module;
    if (!isUpgrade) {
      panel= PixelEndcapName(DetId(id_)).pannelName();
      module= PixelEndcapName(DetId(id_)).plaquetteName();
    } else {
      panel= PixelEndcapNameUpgrade(DetId(id_)).pannelName();
      module= PixelEndcapNameUpgrade(DetId(id_)).plaquetteName();
    }
    
    char slab[80]; sprintf(slab, "Panel_%i_Ring_%i",panel, module);
    hid = src.label() + "_" + slab;
    // Number of digis
    meNDigisRing_ = iBooker.book1D("ndigis_"+hid,"Number of Digis",25,0.,25.);
    meNDigisRing_->setAxisTitle("Number of digis",1);
    // Charge in ADC counts
    meADCRing_ = iBooker.book1D("adc_" + hid,"Digi charge",128,0.,256.);
    meADCRing_->setAxisTitle("ADC counts",1);
    if(!reducedSet)
    {
    if(twoD){
      // 2D hit map
      mePixDigisRing_ = iBooker.book2D("hitmap_"+hid,twodtitle,nbinx,0.,float(ncols_),nbiny,0.,float(nrows_));
      mePixDigisRing_->setAxisTitle("Columns",1);
      mePixDigisRing_->setAxisTitle("Rows",2);
      //std::cout << "During booking: type is "<< type << ", ID is "<< id_ << ", pwd for booking is " << theDMBE->pwd() << ", Plot name: " << hid << std::endl;
    }
    else{
      // projections of 2D hit map
      mePixDigisRing_px_ = iBooker.book1D("hitmap_"+hid+"_px",pxtitle,nbinx,0.,float(ncols_));
      mePixDigisRing_py_ = iBooker.book1D("hitmap_"+hid+"_py",pytitle,nbiny,0.,float(nrows_));
      mePixDigisRing_px_->setAxisTitle("Columns",1);
      mePixDigisRing_py_->setAxisTitle("Rows",1);
    }
    }
  }
}

int SiPixelDigiModule::fill	(	const edm::DetSetVector< PixelDigi > &	input,
		MonitorElement *	combBarrel,
		MonitorElement *	chanBarrel,
		MonitorElement *	chanBarrelL1,
		MonitorElement *	chanBarrelL2,
		MonitorElement *	chanBarrelL3,
		MonitorElement *	chanBarrelL4,
		MonitorElement *	combEndcap,
		const bool	modon,
		const bool	ladon,
		const bool	layon,
		const bool	phion,
		const bool	bladeon,
		const bool	diskon,
		const bool	ringon,
		const bool	twoD,
		const bool	reducedSet,
		const bool	twoDimModOn,
		const bool	twoDimOnlyLayDisk,
		int &	nDigisA,
		int &	nDigisB,
		bool	isUpgrade
	)

Fill histograms.

Definition at line 327 of file SiPixelDigiModule.cc.

References ecalMGPA::adc(), Reference_intrackfit_cff::barrel, edm::DetSetVector< T >::begin(), PixelEndcapNameUpgrade::bladeName(), PixelEndcapName::bladeName(), cuy::col, edm::DetSetVector< T >::end(), Reference_intrackfit_cff::endcap, MonitorElement::Fill(), HcalObjRepresent::Fill(), edm::DetSetVector< T >::find(), i, id_, PixelBarrelNameUpgrade::isHalfModule(), PixelBarrelName::isHalfModule(), PixelBarrelNameUpgrade::ladderName(), PixelBarrelName::ladderName(), PixelBarrelNameUpgrade::layerName(), PixelBarrelName::layerName(), meADC_, meADCBlade_, meADCDisk_, meADCLad_, meADCLay_, meADCPhi_, meADCRing_, meNDigis_, meNDigisBlade_, meNDigisDisk_, meNDigisLad_, meNDigisLay_, meNDigisPhi_, meNDigisRing_, mePixDigis_, mePixDigis_px_, mePixDigis_py_, mePixDigisDisk_, mePixDigisLad_, mePixDigisLad_px_, mePixDigisLad_py_, mePixDigisLay_, mePixDigisLay_px_, mePixDigisLay_py_, mePixDigisPhi_, mePixDigisPhi_px_, mePixDigisPhi_py_, mePixDigisRing_, mePixDigisRing_px_, mePixDigisRing_py_, mePixRocsDisk_, mePixRocsLay_, PixelBarrelName::mI, PixelBarrelName::mO, PixelBarrelNameUpgrade::moduleName(), PixelBarrelName::moduleName(), PixelEndcapNameUpgrade::pannelName(), PixelEndcapName::pannelName(), PixelBarrelName::pI, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, PixelBarrelName::pO, PixelBarrelName::shell(), and DetId::subdetId().

                                                            {
  bool barrel = DetId(id_).subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel);
  bool endcap = DetId(id_).subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap);
  bool isHalfModule = false;
  uint32_t DBladder = 0;
  if(barrel && !isUpgrade){
    isHalfModule = PixelBarrelName(DetId(id_)).isHalfModule(); 
    DBladder = PixelBarrelName(DetId(id_)).ladderName();
  } else if (barrel && isUpgrade) {
    isHalfModule = PixelBarrelNameUpgrade(DetId(id_)).isHalfModule(); 
    DBladder = PixelBarrelNameUpgrade(DetId(id_)).ladderName();
  }

  edm::DetSetVector<PixelDigi>::const_iterator isearch = input.find(id_); // search  digis of detid
  
  unsigned int numberOfDigisMod = 0;
  int msize;
  if (isUpgrade) {msize=10;} else {msize=8;}
  int numberOfDigis[msize]; for(int i=0; i!=msize; i++) numberOfDigis[i]=0; 
  nDigisA=0; nDigisB=0;
  if( isearch != input.end() ) {  // Not an empty iterator
    
    // Look at digis now
    edm::DetSet<PixelDigi>::const_iterator  di;
    for(di = isearch->data.begin(); di != isearch->data.end(); di++) {
      int adc = di->adc();    // charge
      int col = di->column(); // column 
      int row = di->row();    // row
      numberOfDigisMod++;
      
      int DBlayer = 0;
      int DBmodule =0;
      
      if (!isUpgrade) {
      PixelBarrelName::Shell DBshell = PixelBarrelName(DetId(id_)).shell();
        DBlayer  = PixelBarrelName(DetId(id_)).layerName();
        DBmodule = PixelBarrelName(DetId(id_)).moduleName();
      if(barrel){
        if(isHalfModule){
          if(DBshell==PixelBarrelName::pI||DBshell==PixelBarrelName::pO){
        numberOfDigis[0]++; nDigisA++;
        if(DBlayer==1) numberOfDigis[2]++;
        if(DBlayer==2) numberOfDigis[3]++;
        if(DBlayer==3) numberOfDigis[4]++;
      }
          if(DBshell==PixelBarrelName::mI||DBshell==PixelBarrelName::mO){
        numberOfDigis[1]++; nDigisB++;
        if(DBlayer==1) numberOfDigis[5]++;
        if(DBlayer==2) numberOfDigis[6]++;
        if(DBlayer==3) numberOfDigis[7]++;
      }
        }else{
          if(row<80){
        numberOfDigis[0]++; nDigisA++;
        if(DBlayer==1) numberOfDigis[2]++;
        if(DBlayer==2) numberOfDigis[3]++;
        if(DBlayer==3) numberOfDigis[4]++;
      }else{ 
        numberOfDigis[1]++; nDigisB++;
        if(DBlayer==1) numberOfDigis[5]++;
        if(DBlayer==2) numberOfDigis[6]++;
        if(DBlayer==3) numberOfDigis[7]++;
      }
        }
      }
      } else if (isUpgrade) {
        //PixelBarrelNameUpgrade::Shell DBshell = PixelBarrelNameUpgrade(DetId(id_)).shell();
        DBlayer  = PixelBarrelNameUpgrade(DetId(id_)).layerName();
        DBmodule = PixelBarrelNameUpgrade(DetId(id_)).moduleName();
    if(barrel){
      if(row<80){
        numberOfDigis[0]++; nDigisA++;
        if(DBlayer==1) numberOfDigis[2]++;
        if(DBlayer==2) numberOfDigis[3]++;
        if(DBlayer==3) numberOfDigis[4]++;
        if(DBlayer==4) numberOfDigis[5]++;
      }else{ 
        numberOfDigis[1]++; nDigisB++;
        if(DBlayer==1) numberOfDigis[6]++;
        if(DBlayer==2) numberOfDigis[7]++;
        if(DBlayer==3) numberOfDigis[8]++;
        if(DBlayer==4) numberOfDigis[9]++;
      }
    }
      }
      
      if(modon){
    if(!reducedSet){
      if(twoD) {
        if(twoDimModOn) (mePixDigis_)->Fill((float)col,(float)row);
      }
      else {
        (mePixDigis_px_)->Fill((float)col);
        (mePixDigis_py_)->Fill((float)row);
      }
    }
    (meADC_)->Fill((float)adc);
      }
      if(ladon && barrel){
    (meADCLad_)->Fill((float)adc);
    if(!reducedSet){
      if(twoD) (mePixDigisLad_)->Fill((float)col,(float)row);
      else {
      (mePixDigisLad_px_)->Fill((float)col);
      (mePixDigisLad_py_)->Fill((float)row);
      }
    }
      }
      if((layon || twoDimOnlyLayDisk) && barrel){
    if(!twoDimOnlyLayDisk) (meADCLay_)->Fill((float)adc);
    if(!reducedSet){
      if((layon && twoD) || twoDimOnlyLayDisk){
        //ROC histos...
        float rocx = (float)col/52. + 8.0*float(DBmodule-1);
        float rocy = (float)row/160.+float(DBladder);
        //Shift 1st ladder (half modules) up by 1 bin
        if(DBladder==1) rocy = rocy + 0.5;
        mePixRocsLay_->Fill(rocx,rocy);

        if(isHalfModule && DBladder==1){
          (mePixDigisLay_)->Fill((float)col,(float)row+80);
        }
        else (mePixDigisLay_)->Fill((float)col,(float)row);
      }
      if((layon && !twoD) && !twoDimOnlyLayDisk){
        (mePixDigisLay_px_)->Fill((float)col);
        if(isHalfModule && DBladder==1) {
          (mePixDigisLay_py_)->Fill((float)row+80);
        }
        else (mePixDigisLay_py_)->Fill((float)row);
      }
    }
      }
      if(phion && barrel){
    (meADCPhi_)->Fill((float)adc);
    if(!reducedSet)
    {
    if(twoD){
      if(isHalfModule && DBladder==1){
        (mePixDigisPhi_)->Fill((float)col,(float)row+80);
      }
      else (mePixDigisPhi_)->Fill((float)col,(float)row);
    }
    else {
      (mePixDigisPhi_px_)->Fill((float)col);
      if(isHalfModule && DBladder==1) {
        (mePixDigisPhi_py_)->Fill((float)row+80);
      }
      else (mePixDigisPhi_py_)->Fill((float)row);
    }
    }
      }
      if(bladeon && endcap){
    (meADCBlade_)->Fill((float)adc);
      }

      if((diskon || twoDimOnlyLayDisk) && endcap){
    if(!twoDimOnlyLayDisk) (meADCDisk_)->Fill((float)adc);
    if(twoDimOnlyLayDisk){
      (mePixDigisDisk_)->Fill((float)col,(float)row);
      //ROC monitoring
      int DBpanel;
      int DBblade;
      if (!isUpgrade) {
        DBpanel= PixelEndcapName(DetId(id_)).pannelName();
        DBblade= PixelEndcapName(DetId(id_)).bladeName();
      } else {
        DBpanel= PixelEndcapNameUpgrade(DetId(id_)).pannelName();
        DBblade= PixelEndcapNameUpgrade(DetId(id_)).bladeName();
      }
      float offx = 0.;
      //This crazy offset takes into account the roc and module fpix configuration
      for (int i = DBpanel; i < DBmodule; ++i) {offx = offx + float(5+DBpanel-i);}
      float rocx = (float)col/52. + offx + 14.0*float(DBpanel-1);
      float rocy = (float)row/160.+float(DBblade);
      mePixRocsDisk_->Fill(rocx,rocy);
    }
      }
      if(ringon && endcap){
    (meADCRing_)->Fill((float)adc);
    if(!reducedSet)
    {
    if(twoD) (mePixDigisRing_)->Fill((float)col,(float)row);
    else {
      (mePixDigisRing_px_)->Fill((float)col);
      (mePixDigisRing_py_)->Fill((float)row);
    }
    }
      }
    }
    if(modon) (meNDigis_)->Fill((float)numberOfDigisMod);
    if(ladon && barrel) (meNDigisLad_)->Fill((float)numberOfDigisMod);
    if(layon && barrel && !twoDimOnlyLayDisk) (meNDigisLay_)->Fill((float)numberOfDigisMod);
    if(phion && barrel) (meNDigisPhi_)->Fill((float)numberOfDigisMod);
    if(bladeon && endcap) (meNDigisBlade_)->Fill((float)numberOfDigisMod);
    if(diskon && endcap && !twoDimOnlyLayDisk) (meNDigisDisk_)->Fill((float)numberOfDigisMod);
    if(ringon && endcap) (meNDigisRing_)->Fill((float)numberOfDigisMod);
    if(barrel){ 
      if(combBarrel) combBarrel->Fill((float)numberOfDigisMod);
      if(chanBarrel){ if(numberOfDigis[0]>0) chanBarrel->Fill((float)numberOfDigis[0]); if(numberOfDigis[1]>0) chanBarrel->Fill((float)numberOfDigis[1]); }
      if(chanBarrelL1){ if(numberOfDigis[2]>0) chanBarrelL1->Fill((float)numberOfDigis[2]); }
      if(chanBarrelL2){ if(numberOfDigis[3]>0) chanBarrelL2->Fill((float)numberOfDigis[3]); }
      if(chanBarrelL3){ if(numberOfDigis[4]>0) chanBarrelL3->Fill((float)numberOfDigis[4]); }
      if(chanBarrelL4){ if(numberOfDigis[5]>0) chanBarrelL4->Fill((float)numberOfDigis[5]); }
    }else if(endcap){
      if(combEndcap) combEndcap->Fill((float)numberOfDigisMod);
    }
  }
  
  //std::cout<<"numberOfDigis for this module: "<<numberOfDigis<<std::endl;
  return numberOfDigisMod;
}

std::pair< int, int > SiPixelDigiModule::getZeroLoEffROCs

(

)

Definition at line 552 of file SiPixelDigiModule.cc.

References MonitorElement::Fill(), MonitorElement::getBinContent(), MonitorElement::getEntries(), MonitorElement::getNbinsX(), MonitorElement::getNbinsY(), i, j, mePixRocsDisk_, mePixRocsLay_, meZeroOccRocsDisk_, and meZeroOccRocsLay_.

                                                    {
  int nZeroROC = 0;
  int nLoEffROC = 0;
  float SF = 1.0;
  if (mePixRocsDisk_ && meZeroOccRocsDisk_){
    if (mePixRocsDisk_->getEntries() > 0) SF = float(mePixRocsDisk_->getNbinsX()*mePixRocsDisk_->getNbinsY()/mePixRocsDisk_->getEntries());
    for (int i = 1; i < mePixRocsDisk_->getNbinsX(); ++i){
      for (int j = 1; j < mePixRocsDisk_->getNbinsY(); ++j){
    float localX = float(i) - 0.5;
    float localY = float(j)/2.0 + 0.75;
    if (mePixRocsDisk_->getBinContent(i,j)    <  1 ) {nZeroROC++; meZeroOccRocsDisk_->Fill(localX,localY);}
    if (mePixRocsDisk_->getBinContent(i,j)*SF < 0.25){nLoEffROC++;}
      }
    }
    return std::pair<int,int>(nZeroROC,nLoEffROC);
  }
  if (mePixRocsLay_ && meZeroOccRocsLay_){
    if (mePixRocsLay_->getEntries() > 0) SF = float(mePixRocsLay_->getNbinsX()*mePixRocsLay_->getNbinsY()/mePixRocsLay_->getEntries());
    for (int i = 1; i < mePixRocsLay_->getNbinsX(); ++i){
      for (int j = 1; j < mePixRocsLay_->getNbinsY(); ++j){
    float localX = float(i) - 0.5;
    float localY = float(j)/2.0 + 1.25;
    if (mePixRocsLay_->getBinContent(i,j)    <  1 ) {nZeroROC++; meZeroOccRocsLay_->Fill(localX,localY);}
    if (mePixRocsLay_->getBinContent(i,j)*SF < 0.25){nLoEffROC++;}
      }
    }
    return std::pair<int,int>(nZeroROC,nLoEffROC);
  }
  return std::pair<int,int>(0,0);
}

void SiPixelDigiModule::resetRocMap

(

)

Definition at line 546 of file SiPixelDigiModule.cc.

References mePixRocsDisk_, mePixRocsLay_, and MonitorElement::Reset().

                                   {
  if (mePixRocsDisk_) mePixRocsDisk_->Reset();
  if (mePixRocsLay_) mePixRocsLay_->Reset();
}

Member Data Documentation

uint32_t SiPixelDigiModule::id_

private

Definition at line 67 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meADC_

private

Definition at line 71 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meADCBlade_

private

Definition at line 99 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meADCDisk_

private

Definition at line 102 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meADCLad_

private

Definition at line 78 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meADCLay_

private

Definition at line 84 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meADCPhi_

private

Definition at line 92 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meADCRing_

private

Definition at line 108 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meNDigis_

private

Definition at line 70 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meNDigisBlade_

private

Definition at line 98 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meNDigisDisk_

private

Definition at line 101 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meNDigisLad_

private

Definition at line 77 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meNDigisLay_

private

Definition at line 83 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meNDigisPhi_

private

Definition at line 91 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meNDigisRing_

private

Definition at line 107 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigis_

private

Definition at line 72 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigis_px_

private

Definition at line 73 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigis_py_

private

Definition at line 74 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisDisk_

private

Definition at line 103 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisLad_

private

Definition at line 79 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisLad_px_

private

Definition at line 80 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisLad_py_

private

Definition at line 81 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisLay_

private

Definition at line 85 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisLay_px_

private

Definition at line 88 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisLay_py_

private

Definition at line 89 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisPhi_

private

Definition at line 93 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisPhi_px_

private

Definition at line 94 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisPhi_py_

private

Definition at line 95 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisRing_

private

Definition at line 109 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisRing_px_

private

Definition at line 110 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisRing_py_

private

Definition at line 111 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixRocsDisk_ = NULL

private

Definition at line 104 of file SiPixelDigiModule.h.

Referenced by book(), fill(), getZeroLoEffROCs(), and resetRocMap().

MonitorElement* SiPixelDigiModule::mePixRocsLay_ = NULL

private

Definition at line 86 of file SiPixelDigiModule.h.

Referenced by book(), fill(), getZeroLoEffROCs(), and resetRocMap().

MonitorElement* SiPixelDigiModule::meZeroOccRocsDisk_ = NULL

private

Definition at line 105 of file SiPixelDigiModule.h.

Referenced by book(), and getZeroLoEffROCs().

MonitorElement* SiPixelDigiModule::meZeroOccRocsLay_ = NULL

private

Definition at line 87 of file SiPixelDigiModule.h.

Referenced by book(), and getZeroLoEffROCs().

int SiPixelDigiModule::ncols_

private

Definition at line 68 of file SiPixelDigiModule.h.

Referenced by book().

int SiPixelDigiModule::nrows_

private

Definition at line 69 of file SiPixelDigiModule.h.

Referenced by book().