SiPixelDigiModule Class Reference

#include <SiPixelDigiModule.h>

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

Public Member Functions
void	book (const edm::ParameterSet &iConfig, int type=0, bool twoD=true, bool hiRes=false, bool reducedSet=false, bool additInfo=false)
	Book histograms. More...
int	fill (const edm::DetSetVector< PixelDigi > &input, 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)
	Fill histograms. More...
	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_
MonitorElement *	mePixRocsLay_
MonitorElement *	meZeroOccRocsDisk_
MonitorElement *	meZeroOccRocsLay_
int	ncols_
int	nrows_

Detailed Description

Definition at line 33 of file SiPixelDigiModule.h.

Member Typedef Documentation

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

Definition at line 46 of file SiPixelDigiModule.h.

Constructor & Destructor Documentation

SiPixelDigiModule::SiPixelDigiModule

(

)

Default constructor.

Framework.

Definition at line 24 of file SiPixelDigiModule.cc.

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

SiPixelDigiModule::SiPixelDigiModule

(

const uint32_t &

id

)

Constructor with raw DetId.

Definition at line 30 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 37 of file SiPixelDigiModule.cc.

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

SiPixelDigiModule::~SiPixelDigiModule

(

)

Destructor.

Definition at line 46 of file SiPixelDigiModule.cc.

{}

Member Function Documentation

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

Book histograms.

Definition at line 50 of file SiPixelDigiModule.cc.

References Reference_intrackfit_cff::barrel, PixelEndcapName::bladeName(), PixelEndcapName::diskName(), Reference_intrackfit_cff::endcap, edm::ParameterSet::getParameter(), id_, PixelBarrelName::isHalfModule(), edm::InputTag::label(), PixelBarrelName::ladderName(), 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_, PixelBarrelName::moduleName(), ncols_, nrows_, cppFunctionSkipper::operator, PixelEndcapName::pannelName(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, PixelEndcapName::plaquetteName(), MonitorElement::setAxisTitle(), SiPixelHistogramId::setHistoId(), alcazmumu_cfi::src, AlCaHLTBitMon_QueryRunRegistry::string, 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;
  if(barrel){
    isHalfModule = PixelBarrelName(DetId(id_)).isHalfModule(); 
  }

  std::string hid;
  // Get collection name and instantiate Histo Id builder
  edm::InputTag src = iConfig.getParameter<edm::InputTag>( "src" );
  

  // Get DQM interface
  DQMStore* theDMBE = edm::Service<DQMStore>().operator->();

  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_ = theDMBE->book1D(hid,"Number of Digis",25,0.,25.);
    meNDigis_->setAxisTitle("Number of digis",1);
    // Charge in ADC counts
    hid = theHistogramId->setHistoId("adc",id_);
    meADC_ = theDMBE->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_ = theDMBE->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_ = theDMBE->book1D(hid+"_px",pxtitle,nbinx,0.,float(ncols_));
      mePixDigis_py_ = theDMBE->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 = PixelBarrelName(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_ = theDMBE->book1D("ndigis_"+hid,"Number of Digis",25,0.,25.);
    meNDigisLad_->setAxisTitle("Number of digis",1);
    // Charge in ADC counts
    meADCLad_ = theDMBE->book1D("adc_" + hid,"Digi charge",128,0.,256.);
    meADCLad_->setAxisTitle("ADC counts",1);
    if(!reducedSet)
    {
    if(twoD){
      // 2D hit map
      mePixDigisLad_ = theDMBE->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_ = theDMBE->book1D("hitmap_"+hid+"_px",pxtitle,nbinx,0.,float(ncols_));
      mePixDigisLad_py_ = theDMBE->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 = PixelBarrelName(DetId(id_)).layerName();
    char slayer[80]; sprintf(slayer,"Layer_%i",DBlayer);
    hid = src.label() + "_" + slayer;
    if(!additInfo){
      // Number of digis
      meNDigisLay_ = theDMBE->book1D("ndigis_"+hid,"Number of Digis",25,0.,25.);
      meNDigisLay_->setAxisTitle("Number of digis",1);
    // Charge in ADC counts
      meADCLay_ = theDMBE->book1D("adc_" + hid,"Digi charge",128,0.,256.);
      meADCLay_->setAxisTitle("ADC counts",1);
    }
    if(!reducedSet){
      if(twoD || additInfo){
    // 2D hit map
    if(isHalfModule){
      mePixDigisLay_ = theDMBE->book2D("hitmap_"+hid,twodtitle,nbinx,0.,float(ncols_),2*nbiny,0.,float(2*nrows_));
    }
    else{
      mePixDigisLay_ = theDMBE->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_ = theDMBE->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_ = theDMBE->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_ = theDMBE->book1D("hitmap_"+hid+"_px",pxtitle,nbinx,0.,float(ncols_));
    if(isHalfModule){
      mePixDigisLay_py_ = theDMBE->book1D("hitmap_"+hid+"_py",pytitle,2*nbiny,0.,float(2*nrows_));
    }
    else{
      mePixDigisLay_py_ = theDMBE->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 = PixelBarrelName(DetId(id_)).moduleName();
    char smodule[80]; sprintf(smodule,"Ring_%i",DBmodule);
    hid = src.label() + "_" + smodule;
    // Number of digis
    meNDigisPhi_ = theDMBE->book1D("ndigis_"+hid,"Number of Digis",25,0.,25.);
    meNDigisPhi_->setAxisTitle("Number of digis",1);
    // Charge in ADC counts
    meADCPhi_ = theDMBE->book1D("adc_" + hid,"Digi charge",128,0.,256.);
    meADCPhi_->setAxisTitle("ADC counts",1);
    if(!reducedSet)
      {
    if(twoD){
      
      // 2D hit map
      if(isHalfModule){
        mePixDigisPhi_ = theDMBE->book2D("hitmap_"+hid,twodtitle,nbinx,0.,float(ncols_),2*nbiny,0.,float(2*nrows_));
      }
      else {
        mePixDigisPhi_ = theDMBE->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_ = theDMBE->book1D("hitmap_"+hid+"_px",pxtitle,nbinx,0.,float(ncols_));
      if(isHalfModule){
        mePixDigisPhi_py_ = theDMBE->book1D("hitmap_"+hid+"_py",pytitle,2*nbiny,0.,float(2*nrows_));
      }
      else{
        mePixDigisPhi_py_ = theDMBE->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= PixelEndcapName(DetId(id_)).bladeName();
    
    char sblade[80]; sprintf(sblade, "Blade_%02i",blade);
    hid = src.label() + "_" + sblade;
    // Number of digis
    meNDigisBlade_ = theDMBE->book1D("ndigis_"+hid,"Number of Digis",25,0.,25.);
    meNDigisBlade_->setAxisTitle("Number of digis",1);
    // Charge in ADC counts
    meADCBlade_ = theDMBE->book1D("adc_" + hid,"Digi charge",128,0.,256.);
    meADCBlade_->setAxisTitle("ADC counts",1);
  }
  if(type==5 && endcap){
    uint32_t disk = PixelEndcapName(DetId(id_)).diskName();
    
    char sdisk[80]; sprintf(sdisk, "Disk_%i",disk);
    hid = src.label() + "_" + sdisk;
    if(!additInfo){
      // Number of digis
      meNDigisDisk_ = theDMBE->book1D("ndigis_"+hid,"Number of Digis",25,0.,25.);
      meNDigisDisk_->setAxisTitle("Number of digis",1);
      // Charge in ADC counts
      meADCDisk_ = theDMBE->book1D("adc_" + hid,"Digi charge",128,0.,256.);
      meADCDisk_->setAxisTitle("ADC counts",1);
    }
    if(additInfo){
      mePixDigisDisk_ = theDMBE->book2D("hitmap_"+hid,twodtitle,260,0.,260.,160,0.,160.);
      mePixDigisDisk_->setAxisTitle("Columns",1);
      mePixDigisDisk_->setAxisTitle("Rows",2);
      //ROC information in disks
      mePixRocsDisk_  = theDMBE->book2D("rocmap_"+hid,twodroctitle,26,0.,26.,24,1.,13.);
      mePixRocsDisk_ ->setAxisTitle("ROCs per Module (2 Panels)",1);
      mePixRocsDisk_ ->setAxisTitle("Blade Number",2);
      meZeroOccRocsDisk_  = theDMBE->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= PixelEndcapName(DetId(id_)).pannelName();
    uint32_t module= PixelEndcapName(DetId(id_)).plaquetteName();
    char slab[80]; sprintf(slab, "Panel_%i_Ring_%i",panel, module);
    hid = src.label() + "_" + slab;
    // Number of digis
    meNDigisRing_ = theDMBE->book1D("ndigis_"+hid,"Number of Digis",25,0.,25.);
    meNDigisRing_->setAxisTitle("Number of digis",1);
    // Charge in ADC counts
    meADCRing_ = theDMBE->book1D("adc_" + hid,"Digi charge",128,0.,256.);
    meADCRing_->setAxisTitle("ADC counts",1);
    if(!reducedSet)
    {
    if(twoD){
      // 2D hit map
      mePixDigisRing_ = theDMBE->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_ = theDMBE->book1D("hitmap_"+hid+"_px",pxtitle,nbinx,0.,float(ncols_));
      mePixDigisRing_py_ = theDMBE->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,
		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
	)

Fill histograms.

Definition at line 304 of file SiPixelDigiModule.cc.

References ecalMGPA::adc(), Reference_intrackfit_cff::barrel, edm::DetSetVector< T >::begin(), PixelEndcapName::bladeName(), cuy::col, edm::DetSetVector< T >::end(), Reference_intrackfit_cff::endcap, MonitorElement::Fill(), HcalObjRepresent::Fill(), edm::DetSetVector< T >::find(), DQMStore::get(), i, id_, PixelBarrelName::isHalfModule(), PixelBarrelName::ladderName(), 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, PixelBarrelName::moduleName(), cppFunctionSkipper::operator, 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){
    isHalfModule = PixelBarrelName(DetId(id_)).isHalfModule(); 
    DBladder = PixelBarrelName(DetId(id_)).ladderName();
  }

  // Get DQM interface
  DQMStore* theDMBE = edm::Service<DQMStore>().operator->();
  //std::cout<<"id_ = "<<id_<<" , dmbe="<<theDMBE->pwd()<<std::endl;
  //std::cout<<"********************"<<std::endl;
  edm::DetSetVector<PixelDigi>::const_iterator isearch = input.find(id_); // search  digis of detid
  
  unsigned int numberOfDigisMod = 0;
  int numberOfDigis[8]; for(int i=0; i!=8; 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++;
      PixelBarrelName::Shell DBshell = PixelBarrelName(DetId(id_)).shell();
      int DBlayer  = PixelBarrelName(DetId(id_)).layerName();
      int 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]++;
      }
        }
      }
      if(modon){
    if(!reducedSet){
      if(twoD) {
        if(twoDimModOn) (mePixDigis_)->Fill((float)col,(float)row);
        //std::cout << "Col: " << col << ", Row: " << row << ", for DBlayer " << DBlayer << " and isladder " << DBladder << " and module " << PixelBarrelName(DetId(id_)).moduleName() << " and side is " << DBshell << std::endl;
        //std::cout<<"id_ = "<<id_<<" , dmbe="<<theDMBE->pwd()<<std::endl;                                                                                                                  
      }
      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);
        //Copying full 1/2 module to empty 1/2 module...
        //if(isHalfModule) mePixRocsLay_->Fill(rocx,rocy+0.5);
        //end of ROC filling...

        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= PixelEndcapName(DetId(id_)).pannelName();
      int DBblade= PixelEndcapName(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);
      //Now handle the 1/2 module cases by cloning those bins and filling...
      //if (DBpanel==1 && (DBmodule==1||DBmodule==4)){
      //rocy = rocy + 0.5;
      //mePixRocsDisk_->Fill(rocx,rocy);}
      //end ROC monitoring
    }
      }
      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){ 
      MonitorElement* me=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCOMB_Barrel");
      if(me) me->Fill((float)numberOfDigisMod);
      me=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_Barrel");
      if(me){ if(numberOfDigis[0]>0) me->Fill((float)numberOfDigis[0]); if(numberOfDigis[1]>0) me->Fill((float)numberOfDigis[1]); }
      me=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelL1");
      if(me){ if(numberOfDigis[2]>0) me->Fill((float)numberOfDigis[2]); }
      me=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelL2");
      if(me){ if(numberOfDigis[3]>0) me->Fill((float)numberOfDigis[3]); }
      me=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelL3");
      if(me){ if(numberOfDigis[4]>0) me->Fill((float)numberOfDigis[4]); }
    }else if(endcap){
      MonitorElement* me=theDMBE->get("Pixel/Endcap/ALLMODS_ndigisCOMB_Endcap");
      if(me) me->Fill((float)numberOfDigisMod);
    }
  }
  
  //std::cout<<"numberOfDigis for this module: "<<numberOfDigis<<std::endl;
  return numberOfDigisMod;
}

Member Data Documentation

uint32_t SiPixelDigiModule::id_

private

Definition at line 64 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meADC_

private

Definition at line 68 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meADCBlade_

private

Definition at line 96 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meADCDisk_

private

Definition at line 99 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meADCLad_

private

Definition at line 75 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meADCLay_

private

Definition at line 81 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meADCPhi_

private

Definition at line 89 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meADCRing_

private

Definition at line 105 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meNDigis_

private

Definition at line 67 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meNDigisBlade_

private

Definition at line 95 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meNDigisDisk_

private

Definition at line 98 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meNDigisLad_

private

Definition at line 74 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meNDigisLay_

private

Definition at line 80 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meNDigisPhi_

private

Definition at line 88 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meNDigisRing_

private

Definition at line 104 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigis_

private

Definition at line 69 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigis_px_

private

Definition at line 70 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigis_py_

private

Definition at line 71 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisDisk_

private

Definition at line 100 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisLad_

private

Definition at line 76 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisLad_px_

private

Definition at line 77 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisLad_py_

private

Definition at line 78 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisLay_

private

Definition at line 82 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisLay_px_

private

Definition at line 85 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisLay_py_

private

Definition at line 86 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisPhi_

private

Definition at line 90 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisPhi_px_

private

Definition at line 91 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisPhi_py_

private

Definition at line 92 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisRing_

private

Definition at line 106 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisRing_px_

private

Definition at line 107 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixDigisRing_py_

private

Definition at line 108 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixRocsDisk_

private

Definition at line 101 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::mePixRocsLay_

private

Definition at line 83 of file SiPixelDigiModule.h.

Referenced by book(), and fill().

MonitorElement* SiPixelDigiModule::meZeroOccRocsDisk_

private

Definition at line 102 of file SiPixelDigiModule.h.

Referenced by book().

MonitorElement* SiPixelDigiModule::meZeroOccRocsLay_

private

Definition at line 84 of file SiPixelDigiModule.h.

Referenced by book().

int SiPixelDigiModule::ncols_

private

Definition at line 65 of file SiPixelDigiModule.h.

Referenced by book().

int SiPixelDigiModule::nrows_

private

Definition at line 66 of file SiPixelDigiModule.h.

Referenced by book().