SiPixelDQMRocLevelAnalyzer Class Reference

#include <DQM/SiPixelDQMRocLevelAnalyzer/src/SiPixelDQMRocLevelAnalyzer.cc>

Inheritance diagram for SiPixelDQMRocLevelAnalyzer:

Public Member Functions
	SiPixelDQMRocLevelAnalyzer (const edm::ParameterSet &)
	~SiPixelDQMRocLevelAnalyzer ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob ()
virtual void	endJob ()
void	FillRocLevelHistos (TH1F *hrocdep, TH1F *hdist, std::vector< double > &vecx, std::vector< double > &vecy)
void	RocSummary (std::string tagname)
void	RocSumOneModule (int maxr, int maxc, MonitorElement *const &me, std::vector< double > &vecCN, std::vector< double > &vecMean, std::vector< double > &vecSD, int &chipNumber)

Private Attributes
TH1F *	bhGainMean
TH1F *	bhGainMean_dist
TH1F *	bhGainRMS
TH1F *	bhGainRMS_dist
TH1F *	bhNoiseMean
TH1F *	bhNoiseMean_dist
TH1F *	bhNoiseRMS
TH1F *	bhNoiseRMS_dist
TH1F *	bhPedestalMean
TH1F *	bhPedestalMean_dist
TH1F *	bhPedestalRMS
TH1F *	bhPedestalRMS_dist
TH1F *	bhPixelAlive
TH1F *	bhPixelAlive_dist
TH1F *	bhThresholdMean
TH1F *	bhThresholdMean_dist
TH1F *	bhThresholdRMS
TH1F *	bhThresholdRMS_dist
bool	bPixelAlive
bool	bRS
edm::ParameterSet	conf_
DQMStore *	dbe
TH1F *	ehGainMean
TH1F *	ehGainMean_dist
TH1F *	ehGainRMS
TH1F *	ehGainRMS_dist
TH1F *	ehNoiseMean
TH1F *	ehNoiseMean_dist
TH1F *	ehNoiseRMS
TH1F *	ehNoiseRMS_dist
TH1F *	ehPedestalMean
TH1F *	ehPedestalMean_dist
TH1F *	ehPedestalRMS
TH1F *	ehPedestalRMS_dist
TH1F *	ehPixelAlive
TH1F *	ehPixelAlive_dist
TH1F *	ehThresholdMean
TH1F *	ehThresholdMean_dist
TH1F *	ehThresholdRMS
TH1F *	ehThresholdRMS_dist
bool	fRS
edm::Service< TFileService >	fs_
std::vector< MonitorElement * >	mes
std::vector< double >	vbpixCN
std::vector< double >	vbpixM
std::vector< double >	vbpixSD
std::vector< double >	vfpixCN
std::vector< double >	vfpixM
std::vector< double >	vfpixSD

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 47 of file SiPixelDQMRocLevelAnalyzer.h.

Constructor & Destructor Documentation

SiPixelDQMRocLevelAnalyzer::SiPixelDQMRocLevelAnalyzer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 3 of file SiPixelDQMRocLevelAnalyzer.cc.

                                                                                    :conf_(iConfig)
{


}

SiPixelDQMRocLevelAnalyzer::~SiPixelDQMRocLevelAnalyzer

(

)

Definition at line 10 of file SiPixelDQMRocLevelAnalyzer.cc.

{


}

Member Function Documentation

void SiPixelDQMRocLevelAnalyzer::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 18 of file SiPixelDQMRocLevelAnalyzer.cc.

References edm::EventSetup::get(), and edm::Event::getByLabel().

{
   using namespace edm;


#ifdef THIS_IS_AN_EVENT_EXAMPLE
   Handle<ExampleData> pIn;
   iEvent.getByLabel("example",pIn);
#endif
   
#ifdef THIS_IS_AN_EVENTSETUP_EXAMPLE
   ESHandle<SetupData> pSetup;
   iSetup.get<SetupRecord>().get(pSetup);
#endif
}

void SiPixelDQMRocLevelAnalyzer::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 37 of file SiPixelDQMRocLevelAnalyzer.cc.

References bhGainMean, bhGainMean_dist, bhGainRMS, bhGainRMS_dist, bhNoiseMean, bhNoiseMean_dist, bhNoiseRMS, bhNoiseRMS_dist, bhPedestalMean, bhPedestalMean_dist, bhPedestalRMS, bhPedestalRMS_dist, bhPixelAlive, bhPixelAlive_dist, bhThresholdMean, bhThresholdMean_dist, bhThresholdRMS, bhThresholdRMS_dist, bPixelAlive, bRS, conf_, gather_cfg::cout, dbe, ehGainMean, ehGainMean_dist, ehGainRMS, ehGainRMS_dist, ehNoiseMean, ehNoiseMean_dist, ehNoiseRMS, ehNoiseRMS_dist, ehPedestalMean, ehPedestalMean_dist, ehPedestalRMS, ehPedestalRMS_dist, ehPixelAlive, ehPixelAlive_dist, ehThresholdMean, ehThresholdMean_dist, ehThresholdRMS, ehThresholdRMS_dist, lut2db_cfg::filename, FillRocLevelHistos(), fRS, fs_, DQMStore::getAllContents(), edm::ParameterSet::getUntrackedParameter(), i, TFileDirectory::make(), mes, text2workspace::mod, DQMStore::open(), cppFunctionSkipper::operator, RocSummary(), groupFilesInBlocks::temp, vbpixCN, vbpixM, vbpixSD, vfpixCN, vfpixM, and vfpixSD.

{
  std::string filename = conf_.getUntrackedParameter<std::string>("fileName");
  bRS = conf_.getUntrackedParameter<bool>("barrelRocStud");
  fRS = conf_.getUntrackedParameter<bool>("endcapRocStud");
  bPixelAlive = conf_.getUntrackedParameter<bool>("pixelAliveStudy");
  double pixelAliveThresh = conf_.getUntrackedParameter<double>("pixelAliveThreshold");
  bool bThreshold = conf_.getUntrackedParameter<bool>("thresholdStudy");
  bool bNoise = conf_.getUntrackedParameter<bool>("noiseStudy");
  bool bGain = conf_.getUntrackedParameter<bool>("gainStudy");
  bool bPedestal = conf_.getUntrackedParameter<bool>("pedestalStudy");

  if(!bRS && !fRS) return;

  //Open file and get MEs
  dbe = edm::Service<DQMStore>().operator->();
  dbe->open(filename);
  mes = dbe->getAllContents("");
  std::cout << "found " << mes.size() << " monitoring elements!" << std::endl;


  //BARREL
  if(bRS){
    //PIXELALIVE
    if(bPixelAlive){
      bhPixelAlive = fs_->make<TH1F>("bpixAlive","PixelAliveSummary_mean_Barrel",11520, 0., 11520.);
      bhPixelAlive_dist = fs_->make<TH1F>("bpixAliveDist","Mean PixelAliveSummary_mean_Barrel Disbribution", 110, 0., 1.1);
    }
    //THRESHOLD
    if(bThreshold){
      bhThresholdMean = fs_->make<TH1F>("bthreshMean","ScurveThresholdSummary_mean_Barrel",11520,0.,11520.);
      bhThresholdMean_dist = fs_->make<TH1F>("bthreshMeanDist","Mean ScurveThresholdSummary_mean_Barrel Distribution",600,0.,150.);
      bhThresholdRMS = fs_->make<TH1F>("bthreshRMS","ScurveThresholdSummary_RMS_Barrel",11520,0.,11520.);
      bhThresholdRMS_dist = fs_->make<TH1F>("bthreshRMSDist","Mean ScurveThresholdSummary_RMS_Barrel Distribution",800,0.,80.);
    }
    //NOISE
    if(bNoise){
      bhNoiseMean = fs_->make<TH1F>("bnoiseMean","ScurveSigmasSummary_mean_Barrel",11520,0.,11520.);
      bhNoiseMean_dist = fs_->make<TH1F>("bnoiseMeanDist","Mean ScurveSigmasSummary_mean_Barrel Distribution",256,-2.,6.);
      bhNoiseRMS = fs_->make<TH1F>("bnoiseRMS","ScurveSigmasSummary_RMS_Barrel",11520,0.,11520.);
      bhNoiseRMS_dist = fs_->make<TH1F>("bnoiseRMSDist","Mean ScurveSigmasSummary_RMS_Barrel Distribution",768,0.,8.);
    }
    //GAIN
    if(bGain){
      bhGainMean = fs_->make<TH1F>("bgainMean","ScurveGainSummary_mean_Barrel",11520,0.,11520.);
      bhGainMean_dist = fs_->make<TH1F>("bgainMeanDist","Mean ScurveGainSummary_mean_Barrel Distribution",80,0.,8.);
      bhGainRMS = fs_->make<TH1F>("bgainRMS","ScurveGainSummary_RMS_Barrel",11520,0.,11520.);
      bhGainRMS_dist = fs_->make<TH1F>("bgainRMSDist","Mean ScurveGainSummary_RMS_Barrel Distribution",100,0.,10.);
    }
    //PEDESTAL
    if(bPedestal){
      bhPedestalMean = fs_->make<TH1F>("bpedestalMean","ScurvePedestalSummary_mean_Barrel",11520,0.,11520.);
      bhPedestalMean_dist = fs_->make<TH1F>("bpedestalMeanDist","Mean ScurvePedestalSummary_mean_Barrel Distribution",600,0.,300.);
      bhPedestalRMS = fs_->make<TH1F>("bpedestalRMS","ScurvePedestalSummary_RMS_Barrel",11520,0.,11520.);
      bhPedestalRMS_dist = fs_->make<TH1F>("bpedestalRMSDist","Mean ScurvePedestalSummary_RMS_Barrel Distribution",1000,0.,100.);
    }
  }

  //ENDCAP
  if(fRS){
    //PIXELALIVE
    if(bPixelAlive){
      ehPixelAlive = fs_->make<TH1F>("fpixAlive","PixelAliveSummary_mean_Endcap",4320, 0., 4320.);
      ehPixelAlive_dist = fs_->make<TH1F>("fpixAliveDist","Mean PixelAliveSummary_mean_Endcap Disbribution", 110, 0., 1.1);
    }
    //THRESHOLD
    if(bThreshold){
      ehThresholdMean = fs_->make<TH1F>("fthreshMean","ScurveThresholdSummary_mean_Endcap",4320,0.,4320.);
      ehThresholdMean_dist = fs_->make<TH1F>("fthreshMeanDist","Mean ScurveThresholdSummary_mean_Endcap Distribution",600,0.,150.);
      ehThresholdRMS = fs_->make<TH1F>("fthreshRMS","ScurveThresholdSummary_RMS_Endcap",4320,0.,4320.);
      ehThresholdRMS_dist = fs_->make<TH1F>("fthreshRMSDist","Mean ScurveThresholdSummary_RMS_Endcap Distribution",800,0.,80.);
    }
    //NOISE
    if(bNoise){
      ehNoiseMean = fs_->make<TH1F>("fnoiseMean","ScurveSigmasSummary_mean_Endcap",4320,0.,4320.);
      ehNoiseMean_dist = fs_->make<TH1F>("fnoiseMeanDist","Mean ScurveSigmasSummary_mean_Endcap Distribution",256,-2.,6.);
      ehNoiseRMS = fs_->make<TH1F>("fnoiseRMS","ScurveSigmasSummary_RMS_Endcap",4320,0.,4320.);
      ehNoiseRMS_dist = fs_->make<TH1F>("fnoiseRMSDist","Mean ScurveSigmasSummary_RMS_Endcap Distribution",384,0.,4.);
    }
    //GAIN
    if(bGain){
      ehGainMean = fs_->make<TH1F>("fgainMean","ScurveGainSummary_mean_Endcap",4320,0.,4320.);
      ehGainMean_dist = fs_->make<TH1F>("fgainMeanDist","Mean ScurveGainSummary_mean_Endcap Distribution",600,0.,150.);
      ehGainRMS = fs_->make<TH1F>("fgainRMS","ScurveGainSummary_RMS_Endcap",4320,0.,4320.);
      ehGainRMS_dist = fs_->make<TH1F>("fgainRMSDist","Mean ScurveGainSummary_RMS_Endcap Distribution",800,0.,80.);
    }
    //PEDESTAL
    if(bPedestal){
      ehPedestalMean = fs_->make<TH1F>("fpedestalMean","ScurvePedestalSummary_mean_Endcap",4320,0.,4320.);
      ehPedestalMean_dist = fs_->make<TH1F>("fpedestalMeanDist","Mean ScurvePedestalSummary_mean_Endcap Distribution",600,0.,150.);
      ehPedestalRMS = fs_->make<TH1F>("fpedestalRMS","ScurvePedestalSummary_RMS_Endcap",4320,0.,4320.);
      ehPedestalRMS_dist = fs_->make<TH1F>("fpedestalRMSDist","Mean ScurvePedestalSummary_RMS_Endcap Distribution",800,0.,80.);
    }
  }
  //PIXELALIVE
  if(bPixelAlive){
    RocSummary("pixelAlive_siPixelCalibDigis_");
    if(bRS){
      FillRocLevelHistos(bhPixelAlive, bhPixelAlive_dist, vbpixCN, vbpixM);

      //print a list of pixels with pixelAlive quantity below pixelAliveThresh
      for(unsigned int i=0; i<vbpixCN.size(); i++){
    
    if(vbpixM[i]<pixelAliveThresh){
      double temp = vbpixCN[i];
      int shell = (int)((temp-1)/2880); //0 mi, 1 mo, 2 pi, 3 po
      temp -= shell *2880;
      int lay = 1;
      if(temp>576){
        temp -= 576; lay++; 
      } 
      if(temp>960){
        temp -= 960; lay++;
      }
      int lad = 1; 
      if(temp > 32){
        temp -=32; lad++;
      }
      while(temp>64){
        temp -=64; lad++;
      }
      int mod =1; int modsize = 16;
      if(lad ==1 || (lay == 1 && lad == 10) || (lay == 2 && lad == 16) || (lay == 3 && lad == 22) ) modsize = 8; 
      while(temp>modsize){
        temp -= modsize; mod ++;
      }
      std::cout << vbpixCN[i] << " " << vbpixM[i] << ":\n";
      std::cout << "Shell "; 
      switch(shell){
      case 0: std::cout << "mI"; break;
      case 1: std::cout << "mO"; break;
      case 2: std::cout << "pI"; break;
      case 3: std::cout << "pO"; break;
      } 
      std::cout << " Lay" << lay << " Lad" << lad << " Mod" << mod << " Chip" << temp << "\n\n"; 
      
    }
    
      }
    }
    if(fRS){
      FillRocLevelHistos(ehPixelAlive, ehPixelAlive_dist, vfpixCN, vfpixM);
      //print a list of pixels with pixelAlive quantity below pixelAliveThresh
      for(unsigned int i=0; i<vfpixCN.size(); i++){
    if(vfpixM[i]<pixelAliveThresh){
      double temp = vfpixCN[i];
      int hcyl = (int)((temp-1)/1080); //0 mi, 1 mo, 2 pi, 3 po
      temp -= hcyl *1080; 
      int disk = 1; 
      if(temp > 540){
        temp -= 540; disk++;
      }
      int blade = 1;
      while(temp>45){
        temp-=45; blade++;
      }
      int panel = 1, mod = 1;
      if(temp<22){
        //panel 1
        if(temp>16){
          temp-=16; mod = 4;
        }
        else if(temp>8){
          temp-=8; mod = 3;
        }
        else if(temp>2){
          temp-=2; mod = 2;
        }
      } 
      else{
        //panel 2
        temp-=21; panel++;
        if(temp>14){
          temp-=14; mod = 3;
        }
        else if(temp>6){
          temp-=6; mod = 2;
        }
      }
      
      std::cout << vfpixCN[i] << " " << vfpixM[i] << ":\n";
      std::cout << "HalfCylinder "; 
      switch(hcyl){
      case 0: std::cout << "mI"; break;
      case 1: std::cout << "mO"; break;
      case 2: std::cout << "pI"; break;
      case 3: std::cout << "pO"; break;
      } 
      std::cout << " Disk" << disk << " Blade" << blade << " Panel" << panel << " Mod" << mod << " Chip" << temp << "\n\n"; 
      
    }
      }
    }

  }
  //THRESHOLD
  if(bThreshold) {
    vbpixCN.clear(); 
    vbpixM.clear();
    vbpixSD.clear();
    vfpixCN.clear(); 
    vfpixM.clear();
    vfpixSD.clear();

    RocSummary("ScurveThresholds_siPixelCalibDigis_");
    if(bRS){
      FillRocLevelHistos(bhThresholdMean, bhThresholdMean_dist, vbpixCN, vbpixM);
      FillRocLevelHistos(bhThresholdRMS, bhThresholdRMS_dist, vbpixCN, vbpixSD);
  }
    if(fRS){
      FillRocLevelHistos(ehThresholdMean, ehThresholdMean_dist, vfpixCN, vfpixM);
      FillRocLevelHistos(ehThresholdRMS, ehThresholdRMS_dist, vfpixCN, vfpixSD);
    }
  }
  //NOISE
  if(bNoise){
    vbpixCN.clear(); 
    vbpixM.clear();
    vbpixSD.clear();
    vfpixCN.clear(); 
    vfpixM.clear();
    vfpixSD.clear();

    RocSummary("ScurveSigmas_siPixelCalibDigis_");
    if(bRS){
      FillRocLevelHistos(bhNoiseMean, bhNoiseMean_dist, vbpixCN, vbpixM);
      FillRocLevelHistos(bhNoiseRMS, bhNoiseRMS_dist, vbpixCN, vbpixSD);
    }
    if(fRS){
      FillRocLevelHistos(ehNoiseMean, ehNoiseMean_dist, vfpixCN, vfpixM);
      FillRocLevelHistos(ehNoiseRMS, ehNoiseRMS_dist, vfpixCN, vfpixSD);
    }
  }
  //GAIN
  if(bGain){
    vbpixCN.clear(); 
    vbpixM.clear();
    vbpixSD.clear();
    vfpixCN.clear(); 
    vfpixM.clear();
    vfpixSD.clear();

    RocSummary("Gain2d_siPixelCalibDigis_");
    if(bRS){
      FillRocLevelHistos(bhGainMean, bhGainMean_dist, vbpixCN, vbpixM);
      FillRocLevelHistos(bhGainRMS, bhGainRMS_dist, vbpixCN, vbpixSD);
    }
    if(fRS){
      FillRocLevelHistos(ehGainMean, ehGainMean_dist, vfpixCN, vfpixM);
      FillRocLevelHistos(ehGainRMS, ehGainRMS_dist, vfpixCN, vfpixSD);
    }
  }
  //PEDESTAL
  if(bPedestal){
    vbpixCN.clear(); 
    vbpixM.clear();
    vbpixSD.clear();
    vfpixCN.clear(); 
    vfpixM.clear();
    vfpixSD.clear();

    RocSummary("Pedestal2d_siPixelCalibDigis_");
    if(bRS){
      FillRocLevelHistos(bhPedestalMean, bhPedestalMean_dist, vbpixCN, vbpixM);
      FillRocLevelHistos(bhPedestalRMS, bhPedestalRMS_dist, vbpixCN, vbpixSD);
    }
    if(fRS){
      FillRocLevelHistos(ehPedestalMean, ehPedestalMean_dist, vfpixCN, vfpixM);
      FillRocLevelHistos(ehPedestalRMS, ehPedestalRMS_dist, vfpixCN, vfpixSD);
    }
  }
}

void SiPixelDQMRocLevelAnalyzer::endJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 312 of file SiPixelDQMRocLevelAnalyzer.cc.

                                   {
}

void SiPixelDQMRocLevelAnalyzer::FillRocLevelHistos ( TH1F *  hrocdep, TH1F *  hdist, std::vector< double > &  vecx, std::vector< double > &  vecy  )

private

Definition at line 496 of file SiPixelDQMRocLevelAnalyzer.cc.

References i.

Referenced by beginJob().

                                                                                                                          {
  if(vecx.size() == vecy.size()){
    for(unsigned int i=0; i<vecx.size(); i++){
      hrocdep->Fill(vecx[i],vecy[i]);
      hdist->Fill(vecy[i]);
    }
  }
}

void SiPixelDQMRocLevelAnalyzer::RocSummary ( std::string  tagname )

private

Definition at line 316 of file SiPixelDQMRocLevelAnalyzer.cc.

References a, bPixelAlive, bRS, gather_cfg::cout, fRS, i, mes, text2workspace::mod, align::tec::moduleNumber(), mergeVDriftHistosByStation::name, align::tpe::panelNumber(), scaleCards::path, RocSumOneModule(), vbpixCN, vbpixM, vbpixSD, vfpixCN, vfpixM, and vfpixSD.

Referenced by beginJob().

                                                       {
    
  int maxcrow, maxccol;
  std::string name, path = "first"; 
  std::string oldPath = ""; 
  int moduleNumber = 0;
  int bchipNumber = 0; 
  int fchipNumber = 0;
  bool bbarrel = false, bforward = false, bhalfMod = false, bwasHM = false;
  bool bPFFM = false;
  bool bMNCS = false; 
  int panelNumber = -1;

  for(std::vector<MonitorElement*>::const_iterator ime = mes.begin(); ime!=mes.end(); ++ime){
    bwasHM = bhalfMod;
    //set default values
    bMNCS = false;
    bbarrel = false; bforward = false; bhalfMod = false;
    //set name, (old) path
    name = (*ime)->getName();
    oldPath = path; 
    path = (*ime)->getPathname();
 
    //determine module number if any
    if(path.find("Module_")<path.size()){
      if(path!=oldPath) {
    moduleNumber++;
    if(moduleNumber!=1) bMNCS = true;
    
      } 
    }
    else {
      if(moduleNumber>0) bMNCS = true;
      moduleNumber =0;
    }

    
    //find out location (barrel (hm), endcap)
    if(path.find("Barrel/") < path.size()) {
      bbarrel = true;
      if(path.find("H/") < path.size()) bhalfMod = true;
      }
    if(path.find("Endcap/") < path.size()) {
      bforward = true;
      panelNumber = -1;
      if(path.find("Panel_1") < path.size()) panelNumber = 1;
      if(path.find("Panel_2") < path.size()) panelNumber = 2; 

    }
    
    
    //find tagname in histoname
    if(name.find(tagname)<name.size()) bPFFM = true;
    else{
      //adjust chip number if necessary, skip ME
      if(bMNCS){
    if(!bPFFM){
      if(bbarrel && bRS){
        if(bPixelAlive){
          int a = 16; 
          if(bwasHM) a = 8; 
          for(int i=0; i<a; i++){
        bchipNumber++;
        vbpixCN.push_back(bchipNumber);
        vbpixM.push_back(0);
          }
        }
        else{
          if(bwasHM) bchipNumber += 8;
          else bchipNumber += 16; 
        }
      }
      if(bforward && fRS){   

        int maxcol = 2; 
        int mod = moduleNumber; 
        if(mod>1) mod--;
        else{
          if(panelNumber == 1) mod = 4; 
          else mod = 3;
        }
        if(panelNumber==1 && (mod==1 || mod==4)) maxcol = 1;
        if(bPixelAlive){
          for(int i=0; i<maxcol*(mod + panelNumber); i++){
        fchipNumber++;
        vfpixCN.push_back(fchipNumber);
        vfpixM.push_back(0);
          }
        }
        else{
          fchipNumber += maxcol * (mod + panelNumber);
        }
      }
    }
    else bPFFM = false;
      }     
      continue;
    }

    //BARREL ROC LEVEL PLOTS

    if(bbarrel && bRS){
      maxccol=8, maxcrow=2;
      if(bhalfMod) {
    maxcrow=1;
      }

      RocSumOneModule(maxcrow, maxccol, (*ime), vbpixCN, vbpixM, vbpixSD, bchipNumber);

    }


    //ENDCAP ROC LEVEL PLOTS
    if(bforward && fRS){
      maxccol = moduleNumber + panelNumber; 
      maxcrow = 2; 
      if(panelNumber==1 && (moduleNumber==1 || moduleNumber==4)) maxcrow = 1;

      RocSumOneModule(maxcrow, maxccol, (*ime), vfpixCN, vfpixM, vfpixSD, fchipNumber);
    }
  }

  std::cout << "Number of Chips: b" << bchipNumber << " f" << fchipNumber << " " << tagname << std::endl;

}

void SiPixelDQMRocLevelAnalyzer::RocSumOneModule ( int  maxr, int  maxc, MonitorElement *const &  me, std::vector< double > &  vecCN, std::vector< double > &  vecMean, std::vector< double > &  vecSD, int &  chipNumber  )

private

Definition at line 443 of file SiPixelDQMRocLevelAnalyzer.cc.

References bPixelAlive, trackerHits::c, MonitorElement::getBinContent(), timingPdfMaker::mean, Vispa.Plugins.EdmBrowser.EdmDataAccessor::ne(), alignCSCRings::r, mathSSE::sqrt(), and groupFilesInBlocks::temp.

Referenced by RocSummary().

                                                                                                                                                                                         {

  float temp, sum, nentries; 
  for(int cr=0; cr<maxr; cr++){
    for(int cc=0; cc<maxc; cc++){
      //compute mean of 1 ROC
      chipNumber++;
      sum = 0; 
      nentries = 0;
      //sum for 1 ROC
      for(int c=1; c<53; c++){
    for(int r=1; r<81; r++){     
      
      temp = me->getBinContent(52*cc+c, 80*cr+r);
      
      if(temp!=0.){
        sum += temp;
        nentries++;
      }
    }
      }
      if(nentries==0 && bPixelAlive){
    vecCN.push_back(chipNumber);
    vecMean.push_back(0);
      }
      if(nentries!=0){
    double mean = sum/nentries;
    double avsd = 0.; 
    int ne = 0; 
    vecCN.push_back(chipNumber);
    vecMean.push_back(mean);
    
    //computing std dev.
    for(int c=1; c<53; c++){
      for(int r=1; r<81; r++){     
        temp = me->getBinContent(52*cc+c, 80*cr+r);
        if(temp!=0){ 
          avsd += (temp-mean)*(temp-mean);
          ne++;
          
        }
      }
    }
    avsd = avsd/ne; 
    avsd = sqrt(avsd);
    vecSD.push_back(avsd);
      }
      
    }
  }
}

Member Data Documentation

TH1F* SiPixelDQMRocLevelAnalyzer::bhGainMean

private

Definition at line 91 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::bhGainMean_dist

private

Definition at line 92 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::bhGainRMS

private

Definition at line 93 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::bhGainRMS_dist

private

Definition at line 94 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::bhNoiseMean

private

Definition at line 87 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::bhNoiseMean_dist

private

Definition at line 88 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::bhNoiseRMS

private

Definition at line 89 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::bhNoiseRMS_dist

private

Definition at line 90 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::bhPedestalMean

private

Definition at line 95 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::bhPedestalMean_dist

private

Definition at line 96 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::bhPedestalRMS

private

Definition at line 97 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::bhPedestalRMS_dist

private

Definition at line 98 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::bhPixelAlive

private

Definition at line 81 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::bhPixelAlive_dist

private

Definition at line 82 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::bhThresholdMean

private

Definition at line 83 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::bhThresholdMean_dist

private

Definition at line 84 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::bhThresholdRMS

private

Definition at line 85 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::bhThresholdRMS_dist

private

Definition at line 86 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

bool SiPixelDQMRocLevelAnalyzer::bPixelAlive

private

Definition at line 71 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob(), RocSummary(), and RocSumOneModule().

bool SiPixelDQMRocLevelAnalyzer::bRS

private

Definition at line 71 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob(), and RocSummary().

edm::ParameterSet SiPixelDQMRocLevelAnalyzer::conf_

private

Definition at line 66 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

DQMStore* SiPixelDQMRocLevelAnalyzer::dbe

private

Definition at line 67 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::ehGainMean

private

Definition at line 110 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::ehGainMean_dist

private

Definition at line 111 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::ehGainRMS

private

Definition at line 112 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::ehGainRMS_dist

private

Definition at line 113 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::ehNoiseMean

private

Definition at line 106 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::ehNoiseMean_dist

private

Definition at line 107 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::ehNoiseRMS

private

Definition at line 108 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::ehNoiseRMS_dist

private

Definition at line 109 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::ehPedestalMean

private

Definition at line 114 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::ehPedestalMean_dist

private

Definition at line 115 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::ehPedestalRMS

private

Definition at line 116 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::ehPedestalRMS_dist

private

Definition at line 117 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::ehPixelAlive

private

Definition at line 100 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::ehPixelAlive_dist

private

Definition at line 101 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::ehThresholdMean

private

Definition at line 102 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::ehThresholdMean_dist

private

Definition at line 103 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::ehThresholdRMS

private

Definition at line 104 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

TH1F* SiPixelDQMRocLevelAnalyzer::ehThresholdRMS_dist

private

Definition at line 105 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

bool SiPixelDQMRocLevelAnalyzer::fRS

private

Definition at line 71 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob(), and RocSummary().

edm::Service<TFileService> SiPixelDQMRocLevelAnalyzer::fs_

private

Definition at line 68 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

std::vector<MonitorElement*> SiPixelDQMRocLevelAnalyzer::mes

private

Definition at line 70 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob(), and RocSummary().

std::vector<double> SiPixelDQMRocLevelAnalyzer::vbpixCN

private

Definition at line 73 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob(), and RocSummary().

std::vector<double> SiPixelDQMRocLevelAnalyzer::vbpixM

private

Definition at line 74 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob(), and RocSummary().

std::vector<double> SiPixelDQMRocLevelAnalyzer::vbpixSD

private

Definition at line 75 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob(), and RocSummary().

std::vector<double> SiPixelDQMRocLevelAnalyzer::vfpixCN

private

Definition at line 76 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob(), and RocSummary().

std::vector<double> SiPixelDQMRocLevelAnalyzer::vfpixM

private

Definition at line 77 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob(), and RocSummary().

std::vector<double> SiPixelDQMRocLevelAnalyzer::vfpixSD

private

Definition at line 78 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob(), and RocSummary().