SiPixelDQMRocLevelAnalyzer Class Reference

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

Inheritance diagram for SiPixelDQMRocLevelAnalyzer:

Public Types
typedef dqm::legacy::DQMStore	DQMStore
typedef dqm::legacy::MonitorElement	MonitorElement
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
	SiPixelDQMRocLevelAnalyzer (const edm::ParameterSet &)
	~SiPixelDQMRocLevelAnalyzer () override
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
	~EDAnalyzer () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	beginJob () override
void	endJob () override
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
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsInputProcessBlocks ()
static bool	wantsProcessBlocks ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

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

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

Definition at line 44 of file SiPixelDQMRocLevelAnalyzer.h.

Member Typedef Documentation

DQMStore

typedef dqm::legacy::DQMStore SiPixelDQMRocLevelAnalyzer::DQMStore

Definition at line 47 of file SiPixelDQMRocLevelAnalyzer.h.

MonitorElement

typedef dqm::legacy::MonitorElement SiPixelDQMRocLevelAnalyzer::MonitorElement

Definition at line 46 of file SiPixelDQMRocLevelAnalyzer.h.

Constructor & Destructor Documentation

SiPixelDQMRocLevelAnalyzer()

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

explicit

Definition at line 3 of file SiPixelDQMRocLevelAnalyzer.cc.

: conf_(iConfig) {}

~SiPixelDQMRocLevelAnalyzer()

SiPixelDQMRocLevelAnalyzer::~SiPixelDQMRocLevelAnalyzer ( )

override

Definition at line 5 of file SiPixelDQMRocLevelAnalyzer.cc.

{}

Member Function Documentation

analyze()

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

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 8 of file SiPixelDQMRocLevelAnalyzer.cc.

                                                                                            {
  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
}

References edm::EventSetup::get(), get, and iEvent.

beginJob()

void SiPixelDQMRocLevelAnalyzer::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 23 of file SiPixelDQMRocLevelAnalyzer.cc.

                                          {
  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);
    }
  }
}

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, corrVsCorr::filename, FillRocLevelHistos(), fRS, fs_, dqm::implementation::IGetter::getAllContents(), edm::ParameterSet::getUntrackedParameter(), mps_fire::i, createfilelist::int, TFileService::make(), mes, mod(), dqm::implementation::DQMStore::open(), RocSummary(), AlCaHLTBitMon_QueryRunRegistry::string, groupFilesInBlocks::temp, vbpixCN, vbpixM, vbpixSD, vfpixCN, vfpixM, and vfpixSD.

endJob()

void SiPixelDQMRocLevelAnalyzer::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 334 of file SiPixelDQMRocLevelAnalyzer.cc.

{}

Referenced by o2olib.O2ORunMgr::executeJob().

FillRocLevelHistos()

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

private

Definition at line 520 of file SiPixelDQMRocLevelAnalyzer.cc.

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

References mps_fire::i.

Referenced by beginJob().

RocSummary()

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

private

Definition at line 336 of file SiPixelDQMRocLevelAnalyzer.cc.

                                                             {
  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;
}

References a, bPixelAlive, bRS, gather_cfg::cout, fRS, mps_fire::i, mes, mod(), Skims_PA_cff::name, castor_dqm_sourceclient_file_cfg::path, RocSumOneModule(), AlCaHLTBitMon_QueryRunRegistry::string, createPayload::tagname, vbpixCN, vbpixM, vbpixSD, vfpixCN, vfpixM, and vfpixSD.

Referenced by beginJob().

RocSumOneModule()

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 466 of file SiPixelDQMRocLevelAnalyzer.cc.

                                                                  {
  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);
      }
    }
  }
}

References bPixelAlive, HltBtagPostValidation_cff::c, hlt_dqm_clientPB-live_cfg::me, SiStripPI::mean, alignCSCRings::r, mathSSE::sqrt(), and groupFilesInBlocks::temp.

Referenced by RocSummary().

Member Data Documentation

bhGainMean

TH1F* SiPixelDQMRocLevelAnalyzer::bhGainMean

private

Definition at line 93 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

bhGainMean_dist

TH1F* SiPixelDQMRocLevelAnalyzer::bhGainMean_dist

private

Definition at line 94 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

bhGainRMS

TH1F* SiPixelDQMRocLevelAnalyzer::bhGainRMS

private

Definition at line 95 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

bhGainRMS_dist

TH1F* SiPixelDQMRocLevelAnalyzer::bhGainRMS_dist

private

Definition at line 96 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

bhNoiseMean

TH1F* SiPixelDQMRocLevelAnalyzer::bhNoiseMean

private

Definition at line 89 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

bhNoiseMean_dist

TH1F* SiPixelDQMRocLevelAnalyzer::bhNoiseMean_dist

private

Definition at line 90 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

bhNoiseRMS

TH1F* SiPixelDQMRocLevelAnalyzer::bhNoiseRMS

private

Definition at line 91 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

bhNoiseRMS_dist

TH1F* SiPixelDQMRocLevelAnalyzer::bhNoiseRMS_dist

private

Definition at line 92 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

bhPedestalMean

TH1F* SiPixelDQMRocLevelAnalyzer::bhPedestalMean

private

Definition at line 97 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

bhPedestalMean_dist

TH1F* SiPixelDQMRocLevelAnalyzer::bhPedestalMean_dist

private

Definition at line 98 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

bhPedestalRMS

TH1F* SiPixelDQMRocLevelAnalyzer::bhPedestalRMS

private

Definition at line 99 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

bhPedestalRMS_dist

TH1F* SiPixelDQMRocLevelAnalyzer::bhPedestalRMS_dist

private

Definition at line 100 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

bhPixelAlive

TH1F* SiPixelDQMRocLevelAnalyzer::bhPixelAlive

private

Definition at line 83 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

bhPixelAlive_dist

TH1F* SiPixelDQMRocLevelAnalyzer::bhPixelAlive_dist

private

Definition at line 84 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

bhThresholdMean

TH1F* SiPixelDQMRocLevelAnalyzer::bhThresholdMean

private

Definition at line 85 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

bhThresholdMean_dist

TH1F* SiPixelDQMRocLevelAnalyzer::bhThresholdMean_dist

private

Definition at line 86 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

bhThresholdRMS

TH1F* SiPixelDQMRocLevelAnalyzer::bhThresholdRMS

private

Definition at line 87 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

bhThresholdRMS_dist

TH1F* SiPixelDQMRocLevelAnalyzer::bhThresholdRMS_dist

private

Definition at line 88 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

bPixelAlive

bool SiPixelDQMRocLevelAnalyzer::bPixelAlive

private

Definition at line 73 of file SiPixelDQMRocLevelAnalyzer.h.

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

bRS

bool SiPixelDQMRocLevelAnalyzer::bRS

private

Definition at line 73 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob(), and RocSummary().

conf_

edm::ParameterSet SiPixelDQMRocLevelAnalyzer::conf_

private

Definition at line 68 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

dbe

DQMStore* SiPixelDQMRocLevelAnalyzer::dbe

private

Definition at line 69 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

ehGainMean

TH1F* SiPixelDQMRocLevelAnalyzer::ehGainMean

private

Definition at line 112 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

ehGainMean_dist

TH1F* SiPixelDQMRocLevelAnalyzer::ehGainMean_dist

private

Definition at line 113 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

ehGainRMS

TH1F* SiPixelDQMRocLevelAnalyzer::ehGainRMS

private

Definition at line 114 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

ehGainRMS_dist

TH1F* SiPixelDQMRocLevelAnalyzer::ehGainRMS_dist

private

Definition at line 115 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

ehNoiseMean

TH1F* SiPixelDQMRocLevelAnalyzer::ehNoiseMean

private

Definition at line 108 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

ehNoiseMean_dist

TH1F* SiPixelDQMRocLevelAnalyzer::ehNoiseMean_dist

private

Definition at line 109 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

ehNoiseRMS

TH1F* SiPixelDQMRocLevelAnalyzer::ehNoiseRMS

private

Definition at line 110 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

ehNoiseRMS_dist

TH1F* SiPixelDQMRocLevelAnalyzer::ehNoiseRMS_dist

private

Definition at line 111 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

ehPedestalMean

TH1F* SiPixelDQMRocLevelAnalyzer::ehPedestalMean

private

Definition at line 116 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

ehPedestalMean_dist

TH1F* SiPixelDQMRocLevelAnalyzer::ehPedestalMean_dist

private

Definition at line 117 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

ehPedestalRMS

TH1F* SiPixelDQMRocLevelAnalyzer::ehPedestalRMS

private

Definition at line 118 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

ehPedestalRMS_dist

TH1F* SiPixelDQMRocLevelAnalyzer::ehPedestalRMS_dist

private

Definition at line 119 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

ehPixelAlive

TH1F* SiPixelDQMRocLevelAnalyzer::ehPixelAlive

private

Definition at line 102 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

ehPixelAlive_dist

TH1F* SiPixelDQMRocLevelAnalyzer::ehPixelAlive_dist

private

Definition at line 103 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

ehThresholdMean

TH1F* SiPixelDQMRocLevelAnalyzer::ehThresholdMean

private

Definition at line 104 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

ehThresholdMean_dist

TH1F* SiPixelDQMRocLevelAnalyzer::ehThresholdMean_dist

private

Definition at line 105 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

ehThresholdRMS

TH1F* SiPixelDQMRocLevelAnalyzer::ehThresholdRMS

private

Definition at line 106 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

ehThresholdRMS_dist

TH1F* SiPixelDQMRocLevelAnalyzer::ehThresholdRMS_dist

private

Definition at line 107 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

fRS

bool SiPixelDQMRocLevelAnalyzer::fRS

private

Definition at line 73 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob(), and RocSummary().

fs_

edm::Service<TFileService> SiPixelDQMRocLevelAnalyzer::fs_

private

Definition at line 70 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob().

mes

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

private

Definition at line 72 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob(), and RocSummary().

vbpixCN

std::vector<double> SiPixelDQMRocLevelAnalyzer::vbpixCN

private

Definition at line 75 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob(), and RocSummary().

vbpixM

std::vector<double> SiPixelDQMRocLevelAnalyzer::vbpixM

private

Definition at line 76 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob(), and RocSummary().

vbpixSD

std::vector<double> SiPixelDQMRocLevelAnalyzer::vbpixSD

private

Definition at line 77 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob(), and RocSummary().

vfpixCN

std::vector<double> SiPixelDQMRocLevelAnalyzer::vfpixCN

private

Definition at line 78 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob(), and RocSummary().

vfpixM

std::vector<double> SiPixelDQMRocLevelAnalyzer::vfpixM

private

Definition at line 79 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob(), and RocSummary().

vfpixSD

std::vector<double> SiPixelDQMRocLevelAnalyzer::vfpixSD

private

Definition at line 80 of file SiPixelDQMRocLevelAnalyzer.h.

Referenced by beginJob(), and RocSummary().