GlobalDigisAnalyzer Class Reference

#include <GlobalDigisAnalyzer.h>

Inheritance diagram for GlobalDigisAnalyzer:

Public Types
typedef std::vector< double >	DoubleVector
typedef std::vector< float >	FloatVector
typedef std::vector< int >	IntVector
typedef std::map< uint32_t, float, std::less< uint32_t > >	MapType
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T...>
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T...>
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
	GlobalDigisAnalyzer (const edm::ParameterSet &)
	~GlobalDigisAnalyzer () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Protected Member Functions
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Member Functions
void	fillECal (const edm::Event &, const edm::EventSetup &)
void	fillHCal (const edm::Event &, const edm::EventSetup &)
void	fillMuon (const edm::Event &, const edm::EventSetup &)
void	fillTrk (const edm::Event &, const edm::EventSetup &)

Private Attributes
unsigned int	count
edm::EDGetTokenT < CrossingFrame< PCaloHit > >	EBHits_Token_
edm::ESGetToken < EcalADCToGeVConstant, EcalADCToGeVConstantRcd >	ecalADCtoGevToken_
double	ECalbarrelADCtoGeV_
edm::InputTag	ECalEBSrc_
edm::EDGetTokenT < EBDigiCollection >	ECalEBSrc_Token_
edm::InputTag	ECalEESrc_
edm::EDGetTokenT < EEDigiCollection >	ECalEESrc_Token_
double	ECalendcapADCtoGeV_
edm::InputTag	ECalESSrc_
edm::EDGetTokenT < ESDigiCollection >	ECalESSrc_Token_
std::map< int, double, std::less< int > >	ECalgainConv_
edm::EDGetTokenT < CrossingFrame< PCaloHit > >	EEHits_Token_
edm::EDGetTokenT < CrossingFrame< PCaloHit > >	ESHits_Token_
std::string	fName
int	frequency
bool	getAllProvenances
edm::EDGetTokenT < edm::SortedCollection < HBHEDataFrame > >	HBHEDigi_Token_
edm::ESGetToken< HcalDbService, HcalDbRecord >	hcaldbToken_
edm::InputTag	HCalDigi_
edm::InputTag	HCalSrc_
edm::EDGetTokenT < edm::PCaloHitContainer >	HCalSrc_Token_
edm::EDGetTokenT < edm::SortedCollection < HFDataFrame > >	HFDigi_Token_
std::string	hitsProducer
edm::EDGetTokenT < edm::SortedCollection < HODataFrame > >	HODigi_Token_
std::string	label
MonitorElement *	mehCSCStripADC
MonitorElement *	mehCSCStripn
MonitorElement *	mehCSCWiren
MonitorElement *	mehCSCWireTime
MonitorElement *	mehDtMuonLayer [4]
MonitorElement *	mehDtMuonn [4]
MonitorElement *	mehDtMuonTime [4]
MonitorElement *	mehDtMuonTimevLayer [4]
MonitorElement *	mehEcalAEE [2]
MonitorElement *	mehEcalMaxPos [2]
MonitorElement *	mehEcalMultvAEE [2]
MonitorElement *	mehEcaln [2]
MonitorElement *	mehEcalSHE [2]
MonitorElement *	mehEcalSHEvAEESHE [2]
MonitorElement *	mehEScalADC [3]
MonitorElement *	mehEScaln
MonitorElement *	mehHcalAEE [4]
MonitorElement *	mehHcalAEESHE [4]
MonitorElement *	mehHcaln [4]
MonitorElement *	mehHcalSHE [4]
MonitorElement *	mehHcalSHEvAEE [4]
MonitorElement *	mehRPCMuonn
MonitorElement *	mehRPCRes [5]
MonitorElement *	mehSiPixelADC [7]
MonitorElement *	mehSiPixelCol [7]
MonitorElement *	mehSiPixeln [7]
MonitorElement *	mehSiPixelRow [7]
MonitorElement *	mehSiStripADC [19]
MonitorElement *	mehSiStripn [19]
MonitorElement *	mehSiStripStrip [19]
edm::InputTag	MuCSCStripSrc_
edm::EDGetTokenT < CSCStripDigiCollection >	MuCSCStripSrc_Token_
edm::InputTag	MuCSCWireSrc_
edm::EDGetTokenT < CSCWireDigiCollection >	MuCSCWireSrc_Token_
edm::InputTag	MuDTSrc_
edm::EDGetTokenT < DTDigiCollection >	MuDTSrc_Token_
edm::InputTag	MuRPCSrc_
edm::EDGetTokenT < RPCDigiCollection >	MuRPCSrc_Token_
bool	printProvenanceInfo
edm::ESGetToken< RPCGeometry, MuonGeometryRecord >	rpcGeomToken_
edm::EDGetTokenT < edm::PSimHitContainer >	RPCSimHit_Token_
edm::InputTag	SiPxlSrc_
edm::EDGetTokenT < edm::DetSetVector< PixelDigi > >	SiPxlSrc_Token_
edm::InputTag	SiStripSrc_
edm::EDGetTokenT < edm::DetSetVector < SiStripDigi > >	SiStripSrc_Token_
int	theCSCStripPedestalCount
float	theCSCStripPedestalSum
edm::ESGetToken < TrackerTopology, TrackerTopologyRcd >	tTopoToken_
int	verbosity

Additional Inherited Members
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr < DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 106 of file GlobalDigisAnalyzer.h.

Member Typedef Documentation

typedef std::vector<double> GlobalDigisAnalyzer::DoubleVector

Definition at line 109 of file GlobalDigisAnalyzer.h.

typedef std::vector<float> GlobalDigisAnalyzer::FloatVector

Definition at line 108 of file GlobalDigisAnalyzer.h.

typedef std::vector<int> GlobalDigisAnalyzer::IntVector

Definition at line 110 of file GlobalDigisAnalyzer.h.

typedef std::map<uint32_t, float, std::less<uint32_t> > GlobalDigisAnalyzer::MapType

Definition at line 111 of file GlobalDigisAnalyzer.h.

Constructor & Destructor Documentation

GlobalDigisAnalyzer::GlobalDigisAnalyzer ( const edm::ParameterSet &  iPSet )

explicit

Definition at line 14 of file GlobalDigisAnalyzer.cc.

References EBHits_Token_, ecalADCtoGevToken_, ECalbarrelADCtoGeV_, ECalEBSrc_, ECalEBSrc_Token_, ECalEESrc_, ECalEESrc_Token_, ECalendcapADCtoGeV_, ECalESSrc_, ECalESSrc_Token_, ECalgainConv_, EEHits_Token_, DeDxTools::esConsumes(), ESHits_Token_, fName, frequency, EcalMGPAGainRatio::gain12Over6(), EcalMGPAGainRatio::gain6Over1(), getAllProvenances, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HBHEDigi_Token_, hcaldbToken_, HCalDigi_, HCalSrc_, HCalSrc_Token_, HFDigi_Token_, hitsProducer, HODigi_Token_, HLT_FULL_cff::InputTag, edm::InputTag::instance(), edm::InputTag::label(), MuCSCStripSrc_, MuCSCStripSrc_Token_, MuCSCWireSrc_, MuCSCWireSrc_Token_, MuDTSrc_, MuDTSrc_Token_, MuRPCSrc_, MuRPCSrc_Token_, printProvenanceInfo, rpcGeomToken_, RPCSimHit_Token_, SiPxlSrc_, SiPxlSrc_Token_, SiStripSrc_, SiStripSrc_Token_, AlCaHLTBitMon_QueryRunRegistry::string, tTopoToken_, and verbosity.

    : fName(""),
      verbosity(0),
      frequency(0),
      label(""),
      getAllProvenances(false),
      printProvenanceInfo(false),
      hitsProducer(""),
      theCSCStripPedestalSum(0),
      theCSCStripPedestalCount(0),
      count(0) {
  std::string MsgLoggerCat = "GlobalDigisAnalyzer_GlobalDigisAnalyzer";

  // get information from parameter set
  fName = iPSet.getUntrackedParameter<std::string>("Name");
  verbosity = iPSet.getUntrackedParameter<int>("Verbosity");
  frequency = iPSet.getUntrackedParameter<int>("Frequency");
  edm::ParameterSet m_Prov = iPSet.getParameter<edm::ParameterSet>("ProvenanceLookup");
  getAllProvenances = m_Prov.getUntrackedParameter<bool>("GetAllProvenances");
  printProvenanceInfo = m_Prov.getUntrackedParameter<bool>("PrintProvenanceInfo");
  hitsProducer = iPSet.getParameter<std::string>("hitsProducer");

  // get Labels to use to extract information
  ECalEBSrc_ = iPSet.getParameter<edm::InputTag>("ECalEBSrc");
  ECalEESrc_ = iPSet.getParameter<edm::InputTag>("ECalEESrc");
  ECalESSrc_ = iPSet.getParameter<edm::InputTag>("ECalESSrc");
  HCalSrc_ = iPSet.getParameter<edm::InputTag>("HCalSrc");
  HCalDigi_ = iPSet.getParameter<edm::InputTag>("HCalDigi");
  SiStripSrc_ = iPSet.getParameter<edm::InputTag>("SiStripSrc");
  SiPxlSrc_ = iPSet.getParameter<edm::InputTag>("SiPxlSrc");
  MuDTSrc_ = iPSet.getParameter<edm::InputTag>("MuDTSrc");
  MuCSCStripSrc_ = iPSet.getParameter<edm::InputTag>("MuCSCStripSrc");
  MuCSCWireSrc_ = iPSet.getParameter<edm::InputTag>("MuCSCWireSrc");
  MuRPCSrc_ = iPSet.getParameter<edm::InputTag>("MuRPCSrc");

  ECalEBSrc_Token_ = consumes<EBDigiCollection>(iPSet.getParameter<edm::InputTag>("ECalEBSrc"));
  ECalEESrc_Token_ = consumes<EEDigiCollection>(iPSet.getParameter<edm::InputTag>("ECalEESrc"));
  ECalESSrc_Token_ = consumes<ESDigiCollection>(iPSet.getParameter<edm::InputTag>("ECalESSrc"));
  HCalSrc_Token_ = consumes<edm::PCaloHitContainer>(iPSet.getParameter<edm::InputTag>("HCalSrc"));
  HBHEDigi_Token_ = consumes<edm::SortedCollection<HBHEDataFrame>>(iPSet.getParameter<edm::InputTag>("HCalDigi"));
  HODigi_Token_ = consumes<edm::SortedCollection<HODataFrame>>(iPSet.getParameter<edm::InputTag>("HCalDigi"));
  HFDigi_Token_ = consumes<edm::SortedCollection<HFDataFrame>>(iPSet.getParameter<edm::InputTag>("HCalDigi"));
  SiStripSrc_Token_ = consumes<edm::DetSetVector<SiStripDigi>>(iPSet.getParameter<edm::InputTag>("SiStripSrc"));
  SiPxlSrc_Token_ = consumes<edm::DetSetVector<PixelDigi>>(iPSet.getParameter<edm::InputTag>("SiPxlSrc"));
  MuDTSrc_Token_ = consumes<DTDigiCollection>(iPSet.getParameter<edm::InputTag>("MuDTSrc"));
  MuCSCStripSrc_Token_ = consumes<CSCStripDigiCollection>(iPSet.getParameter<edm::InputTag>("MuCSCStripSrc"));
  MuCSCWireSrc_Token_ = consumes<CSCWireDigiCollection>(iPSet.getParameter<edm::InputTag>("MuCSCWireSrc"));
  MuRPCSrc_Token_ = consumes<RPCDigiCollection>(iPSet.getParameter<edm::InputTag>("MuRPCSrc"));
  //
  const std::string barrelHitsName(hitsProducer + "EcalHitsEB");
  const std::string endcapHitsName(hitsProducer + "EcalHitsEE");
  const std::string preshowerHitsName(hitsProducer + "EcalHitsES");
  EBHits_Token_ = consumes<CrossingFrame<PCaloHit>>(edm::InputTag(std::string("mix"), std::string("barrelHitsName")));
  EEHits_Token_ = consumes<CrossingFrame<PCaloHit>>(edm::InputTag(std::string("mix"), std::string("endcapHitsName")));
  ESHits_Token_ =
      consumes<CrossingFrame<PCaloHit>>(edm::InputTag(std::string("mix"), std::string("preshowerHitsName")));

  RPCSimHit_Token_ =
      consumes<edm::PSimHitContainer>(edm::InputTag(std::string("g4SimHits"), std::string("MuonRPCHits")));

  ecalADCtoGevToken_ = esConsumes();
  rpcGeomToken_ = esConsumes();
  tTopoToken_ = esConsumes();
  hcaldbToken_ = esConsumes();

  // use value of first digit to determine default output level (inclusive)
  // 0 is none, 1 is basic, 2 is fill output, 3 is gather output
  verbosity %= 10;

  // print out Parameter Set information being used
  if (verbosity >= 0) {
    edm::LogInfo(MsgLoggerCat) << "\n===============================\n"
                               << "Initialized as EDAnalyzer with parameter values:\n"
                               << "    Name          = " << fName << "\n"
                               << "    Verbosity     = " << verbosity << "\n"
                               << "    Frequency     = " << frequency << "\n"
                               << "    GetProv       = " << getAllProvenances << "\n"
                               << "    PrintProv     = " << printProvenanceInfo << "\n"
                               << "    ECalEBSrc     = " << ECalEBSrc_.label() << ":" << ECalEBSrc_.instance() << "\n"
                               << "    ECalEESrc     = " << ECalEESrc_.label() << ":" << ECalEESrc_.instance() << "\n"
                               << "    ECalESSrc     = " << ECalESSrc_.label() << ":" << ECalESSrc_.instance() << "\n"
                               << "    HCalSrc       = " << HCalSrc_.label() << ":" << HCalSrc_.instance() << "\n"
                               << "    HCalDigi       = " << HCalDigi_.label() << ":" << HCalDigi_.instance() << "\n"
                               << "    SiStripSrc    = " << SiStripSrc_.label() << ":" << SiStripSrc_.instance() << "\n"
                               << "    SiPixelSrc    = " << SiPxlSrc_.label() << ":" << SiPxlSrc_.instance() << "\n"
                               << "    MuDTSrc       = " << MuDTSrc_.label() << ":" << MuDTSrc_.instance() << "\n"
                               << "    MuCSCStripSrc = " << MuCSCStripSrc_.label() << ":" << MuCSCStripSrc_.instance()
                               << "\n"
                               << "    MuCSCWireSrc  = " << MuCSCWireSrc_.label() << ":" << MuCSCWireSrc_.instance()
                               << "\n"
                               << "    MuRPCSrc      = " << MuRPCSrc_.label() << ":" << MuRPCSrc_.instance() << "\n"
                               << "===============================\n";
  }

  // set default constants
  // ECal

  ECalbarrelADCtoGeV_ = 0.035;
  ECalendcapADCtoGeV_ = 0.06;

  EcalMGPAGainRatio defaultRatios;
  ECalgainConv_[0] = 0.;
  ECalgainConv_[1] = 1.;
  ECalgainConv_[2] = defaultRatios.gain12Over6();
  ECalgainConv_[3] = ECalgainConv_[2] * (defaultRatios.gain6Over1());

  if (verbosity >= 0) {
    edm::LogInfo(MsgLoggerCat) << "Modified Calorimeter gain constants: g0 = " << ECalgainConv_[0]
                               << ", g1 = " << ECalgainConv_[1] << ", g2 = " << ECalgainConv_[2]
                               << ", g3 = " << ECalgainConv_[3];
  }
}

GlobalDigisAnalyzer::~GlobalDigisAnalyzer ( )

override

Definition at line 127 of file GlobalDigisAnalyzer.cc.

{}

Member Function Documentation

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 427 of file GlobalDigisAnalyzer.cc.

References count, ecalADCtoGevToken_, ECalbarrelADCtoGeV_, ECalendcapADCtoGeV_, edm::EventID::event(), fillECal(), fillHCal(), fillMuon(), fillTrk(), frequency, getAllProvenances, edm::Event::getAllStableProvenance(), edm::EventSetup::getData(), EcalADCToGeVConstant::getEBValue(), EcalADCToGeVConstant::getEEValue(), mps_fire::i, edm::EventBase::id(), nevt, printProvenanceInfo, edm::EventID::run(), AlCaHLTBitMon_QueryRunRegistry::string, and verbosity.

                                                                                     {
  std::string MsgLoggerCat = "GlobalDigisAnalyzer_analyze";

  // keep track of number of events processed
  ++count;

  // THIS BLOCK MIGRATED HERE FROM beginJob:
  // setup calorimeter constants from service
  const EcalADCToGeVConstant *agc = &iSetup.getData(ecalADCtoGevToken_);
  ECalbarrelADCtoGeV_ = agc->getEBValue();
  ECalendcapADCtoGeV_ = agc->getEEValue();
  if (verbosity >= 0) {
    edm::LogInfo(MsgLoggerCat) << "Modified Calorimeter ADCtoGeV constants: barrel = " << ECalbarrelADCtoGeV_
                               << ", endcap = " << ECalendcapADCtoGeV_;
  }

  // get event id information
  edm::RunNumber_t nrun = iEvent.id().run();
  edm::EventNumber_t nevt = iEvent.id().event();

  if (verbosity > 0) {
    edm::LogInfo(MsgLoggerCat) << "Processing run " << nrun << ", event " << nevt << " (" << count << " events total)";
  } else if (verbosity == 0) {
    if (nevt % frequency == 0 || nevt == 1) {
      edm::LogInfo(MsgLoggerCat) << "Processing run " << nrun << ", event " << nevt << " (" << count
                                 << " events total)";
    }
  }

  // look at information available in the event
  if (getAllProvenances) {
    std::vector<const edm::StableProvenance *> AllProv;
    iEvent.getAllStableProvenance(AllProv);

    if (verbosity >= 0)
      edm::LogInfo(MsgLoggerCat) << "Number of Provenances = " << AllProv.size();

    if (printProvenanceInfo && (verbosity >= 0)) {
      TString eventout("\nProvenance info:\n");

      for (unsigned int i = 0; i < AllProv.size(); ++i) {
        eventout += "\n       ******************************";
        eventout += "\n       Module       : ";
        eventout += AllProv[i]->moduleLabel();
        eventout += "\n       ProductID    : ";
        eventout += AllProv[i]->productID().id();
        eventout += "\n       ClassName    : ";
        eventout += AllProv[i]->className();
        eventout += "\n       InstanceName : ";
        eventout += AllProv[i]->productInstanceName();
        eventout += "\n       BranchName   : ";
        eventout += AllProv[i]->branchName();
      }
      eventout += "\n       ******************************\n";
      edm::LogInfo(MsgLoggerCat) << eventout << "\n";
      printProvenanceInfo = false;
    }
    getAllProvenances = false;
  }

  // call fill functions
  // gather Ecal information from event
  fillECal(iEvent, iSetup);
  // gather Hcal information from event
  fillHCal(iEvent, iSetup);
  // gather Track information from event
  fillTrk(iEvent, iSetup);
  // gather Muon information from event
  fillMuon(iEvent, iSetup);

  if (verbosity > 0)
    edm::LogInfo(MsgLoggerCat) << "Done gathering data from event.";

  if (verbosity > 2)
    edm::LogInfo(MsgLoggerCat) << "Saving event contents:";

  return;
}

void GlobalDigisAnalyzer::bookHistograms ( DQMStore::IBooker &  iBooker, edm::Run const &  , edm::EventSetup const &    )

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 129 of file GlobalDigisAnalyzer.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::bookProfile(), mps_fire::i, dqmiolumiharvest::j, mehCSCStripADC, mehCSCStripn, mehCSCWiren, mehCSCWireTime, mehDtMuonLayer, mehDtMuonn, mehDtMuonTime, mehDtMuonTimevLayer, mehEcalAEE, mehEcalMaxPos, mehEcalMultvAEE, mehEcaln, mehEcalSHE, mehEcalSHEvAEESHE, mehEScalADC, mehEScaln, mehHcalAEE, mehHcalAEESHE, mehHcaln, mehHcalSHE, mehHcalSHEvAEE, mehRPCMuonn, mehRPCRes, mehSiPixelADC, mehSiPixelCol, mehSiPixeln, mehSiPixelRow, mehSiStripADC, mehSiStripn, mehSiStripStrip, dqm::impl::MonitorElement::setAxisTitle(), dqm::implementation::NavigatorBase::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                        {
  // Si Strip
  std::string SiStripString[19] = {"TECW1",
                                   "TECW2",
                                   "TECW3",
                                   "TECW4",
                                   "TECW5",
                                   "TECW6",
                                   "TECW7",
                                   "TECW8",
                                   "TIBL1",
                                   "TIBL2",
                                   "TIBL3",
                                   "TIBL4",
                                   "TIDW1",
                                   "TIDW2",
                                   "TIDW3",
                                   "TOBL1",
                                   "TOBL2",
                                   "TOBL3",
                                   "TOBL4"};
  for (int i = 0; i < 19; ++i) {
    mehSiStripn[i] = nullptr;
    mehSiStripADC[i] = nullptr;
    mehSiStripStrip[i] = nullptr;
  }
  std::string hcharname, hchartitle;
  iBooker.setCurrentFolder("GlobalDigisV/SiStrips");
  for (int amend = 0; amend < 19; ++amend) {
    hcharname = "hSiStripn_" + SiStripString[amend];
    hchartitle = SiStripString[amend] + "  Digis";
    mehSiStripn[amend] = iBooker.book1D(hcharname, hchartitle, 5000, 0., 10000.);
    mehSiStripn[amend]->setAxisTitle("Number of Digis", 1);
    mehSiStripn[amend]->setAxisTitle("Count", 2);

    hcharname = "hSiStripADC_" + SiStripString[amend];
    hchartitle = SiStripString[amend] + " ADC";
    mehSiStripADC[amend] = iBooker.book1D(hcharname, hchartitle, 150, 0.0, 300.);
    mehSiStripADC[amend]->setAxisTitle("ADC", 1);
    mehSiStripADC[amend]->setAxisTitle("Count", 2);

    hcharname = "hSiStripStripADC_" + SiStripString[amend];
    hchartitle = SiStripString[amend] + " Strip";
    mehSiStripStrip[amend] = iBooker.book1D(hcharname, hchartitle, 200, 0.0, 800.);
    mehSiStripStrip[amend]->setAxisTitle("Strip Number", 1);
    mehSiStripStrip[amend]->setAxisTitle("Count", 2);
  }

  // HCal
  std::string HCalString[4] = {"HB", "HE", "HO", "HF"};
  float calnUpper[4] = {30000., 30000., 30000., 20000.};
  float calnLower[4] = {0., 0., 0., 0.};
  float SHEUpper[4] = {1., 1., 1., 1.};
  float SHEvAEEUpper[4] = {5000, 5000, 5000, 5000};
  float SHEvAEELower[4] = {-5000, -5000, -5000, -5000};
  int SHEvAEEnBins[4] = {200, 200, 200, 200};
  double ProfileUpper[4] = {1., 1., 1., 1.};

  for (int i = 0; i < 4; ++i) {
    mehHcaln[i] = nullptr;
    mehHcalAEE[i] = nullptr;
    mehHcalSHE[i] = nullptr;
    mehHcalAEESHE[i] = nullptr;
    mehHcalSHEvAEE[i] = nullptr;
  }
  iBooker.setCurrentFolder("GlobalDigisV/HCals");

  for (int amend = 0; amend < 4; ++amend) {
    hcharname = "hHcaln_" + HCalString[amend];
    hchartitle = HCalString[amend] + "  digis";
    mehHcaln[amend] = iBooker.book1D(hcharname, hchartitle, 10000, calnLower[amend], calnUpper[amend]);
    mehHcaln[amend]->setAxisTitle("Number of Digis", 1);
    mehHcaln[amend]->setAxisTitle("Count", 2);

    hcharname = "hHcalAEE_" + HCalString[amend];
    hchartitle = HCalString[amend] + "Cal AEE";
    mehHcalAEE[amend] = iBooker.book1D(hcharname, hchartitle, 60, -10., 50.);
    mehHcalAEE[amend]->setAxisTitle("Analog Equivalent Energy", 1);
    mehHcalAEE[amend]->setAxisTitle("Count", 2);

    hcharname = "hHcalSHE_" + HCalString[amend];
    hchartitle = HCalString[amend] + "Cal SHE";
    mehHcalSHE[amend] = iBooker.book1D(hcharname, hchartitle, 1000, 0.0, SHEUpper[amend]);
    mehHcalSHE[amend]->setAxisTitle("Simulated Hit Energy", 1);
    mehHcalSHE[amend]->setAxisTitle("Count", 2);

    hcharname = "hHcalAEESHE_" + HCalString[amend];
    hchartitle = HCalString[amend] + "Cal AEE/SHE";
    mehHcalAEESHE[amend] =
        iBooker.book1D(hcharname, hchartitle, SHEvAEEnBins[amend], SHEvAEELower[amend], SHEvAEEUpper[amend]);
    mehHcalAEESHE[amend]->setAxisTitle("ADC / SHE", 1);
    mehHcalAEESHE[amend]->setAxisTitle("Count", 2);

    hcharname = "hHcalSHEvAEE_" + HCalString[amend];
    hchartitle = HCalString[amend] + "Cal SHE vs. AEE";
    mehHcalSHEvAEE[amend] =
        iBooker.bookProfile(hcharname, hchartitle, 60, -10., 50., 100, 0., (float)ProfileUpper[amend], "");
    mehHcalSHEvAEE[amend]->setAxisTitle("AEE / SHE", 1);
    mehHcalSHEvAEE[amend]->setAxisTitle("SHE", 2);
  }

  // Ecal
  std::string ECalString[2] = {"EB", "EE"};

  for (int i = 0; i < 2; ++i) {
    mehEcaln[i] = nullptr;
    mehEcalAEE[i] = nullptr;
    mehEcalSHE[i] = nullptr;
    mehEcalMaxPos[i] = nullptr;
    mehEcalMultvAEE[i] = nullptr;
    mehEcalSHEvAEESHE[i] = nullptr;
  }
  iBooker.setCurrentFolder("GlobalDigisV/ECals");

  for (int amend = 0; amend < 2; ++amend) {
    hcharname = "hEcaln_" + ECalString[amend];
    hchartitle = ECalString[amend] + "  digis";
    mehEcaln[amend] = iBooker.book1D(hcharname, hchartitle, 3000, 0., 40000.);
    mehEcaln[amend]->setAxisTitle("Number of Digis", 1);
    mehEcaln[amend]->setAxisTitle("Count", 2);

    hcharname = "hEcalAEE_" + ECalString[amend];
    hchartitle = ECalString[amend] + "Cal AEE";
    mehEcalAEE[amend] = iBooker.book1D(hcharname, hchartitle, 1000, 0., 100.);
    mehEcalAEE[amend]->setAxisTitle("Analog Equivalent Energy", 1);
    mehEcalAEE[amend]->setAxisTitle("Count", 2);

    hcharname = "hEcalSHE_" + ECalString[amend];
    hchartitle = ECalString[amend] + "Cal SHE";
    mehEcalSHE[amend] = iBooker.book1D(hcharname, hchartitle, 500, 0., 50.);
    mehEcalSHE[amend]->setAxisTitle("Simulated Hit Energy", 1);
    mehEcalSHE[amend]->setAxisTitle("Count", 2);

    hcharname = "hEcalMaxPos_" + ECalString[amend];
    hchartitle = ECalString[amend] + "Cal MaxPos";
    mehEcalMaxPos[amend] = iBooker.book1D(hcharname, hchartitle, 10, 0., 10.);
    mehEcalMaxPos[amend]->setAxisTitle("Maximum Position", 1);
    mehEcalMaxPos[amend]->setAxisTitle("Count", 2);

    hcharname = "hEcalSHEvAEESHE_" + ECalString[amend];
    hchartitle = ECalString[amend] + "Cal SHE vs. AEE/SHE";
    mehEcalSHEvAEESHE[amend] = iBooker.bookProfile(hcharname, hchartitle, 1000, 0., 100., 500, 0., 50., "");
    mehEcalSHEvAEESHE[amend]->setAxisTitle("AEE / SHE", 1);
    mehEcalSHEvAEESHE[amend]->setAxisTitle("SHE", 2);

    hcharname = "hEcalMultvAEE_" + ECalString[amend];
    hchartitle = ECalString[amend] + "Cal Multi vs. AEE";
    mehEcalMultvAEE[amend] = iBooker.bookProfile(hcharname, hchartitle, 1000, 0., 100., 4000, 0., 40000., "");
    mehEcalMultvAEE[amend]->setAxisTitle("Analog Equivalent Energy", 1);
    mehEcalMultvAEE[amend]->setAxisTitle("Number of Digis", 2);
  }
  mehEScaln = nullptr;

  hcharname = "hEcaln_ES";
  hchartitle = "ESCAL  digis";
  mehEScaln = iBooker.book1D(hcharname, hchartitle, 1000, 0., 5000.);
  mehEScaln->setAxisTitle("Number of Digis", 1);
  mehEScaln->setAxisTitle("Count", 2);

  std::string ADCNumber[3] = {"0", "1", "2"};
  for (int i = 0; i < 3; ++i) {
    mehEScalADC[i] = nullptr;
    hcharname = "hEcalADC" + ADCNumber[i] + "_ES";
    hchartitle = "ESCAL  ADC" + ADCNumber[i];
    mehEScalADC[i] = iBooker.book1D(hcharname, hchartitle, 1500, 0., 1500.);
    mehEScalADC[i]->setAxisTitle("ADC" + ADCNumber[i], 1);
    mehEScalADC[i]->setAxisTitle("Count", 2);
  }

  // Si Pixels ***DONE***
  std::string SiPixelString[7] = {"BRL1", "BRL2", "BRL3", "FWD1n", "FWD1p", "FWD2n", "FWD2p"};
  for (int j = 0; j < 7; ++j) {
    mehSiPixeln[j] = nullptr;
    mehSiPixelADC[j] = nullptr;
    mehSiPixelRow[j] = nullptr;
    mehSiPixelCol[j] = nullptr;
  }

  iBooker.setCurrentFolder("GlobalDigisV/SiPixels");
  for (int amend = 0; amend < 7; ++amend) {
    hcharname = "hSiPixeln_" + SiPixelString[amend];
    hchartitle = SiPixelString[amend] + " Digis";
    if (amend < 3)
      mehSiPixeln[amend] = iBooker.book1D(hcharname, hchartitle, 500, 0., 1000.);
    else
      mehSiPixeln[amend] = iBooker.book1D(hcharname, hchartitle, 500, 0., 1000.);
    mehSiPixeln[amend]->setAxisTitle("Number of Digis", 1);
    mehSiPixeln[amend]->setAxisTitle("Count", 2);

    hcharname = "hSiPixelADC_" + SiPixelString[amend];
    hchartitle = SiPixelString[amend] + " ADC";
    mehSiPixelADC[amend] = iBooker.book1D(hcharname, hchartitle, 150, 0.0, 300.);
    mehSiPixelADC[amend]->setAxisTitle("ADC", 1);
    mehSiPixelADC[amend]->setAxisTitle("Count", 2);

    hcharname = "hSiPixelRow_" + SiPixelString[amend];
    hchartitle = SiPixelString[amend] + " Row";
    mehSiPixelRow[amend] = iBooker.book1D(hcharname, hchartitle, 100, 0.0, 100.);
    mehSiPixelRow[amend]->setAxisTitle("Row Number", 1);
    mehSiPixelRow[amend]->setAxisTitle("Count", 2);

    hcharname = "hSiPixelColumn_" + SiPixelString[amend];
    hchartitle = SiPixelString[amend] + " Column";
    mehSiPixelCol[amend] = iBooker.book1D(hcharname, hchartitle, 200, 0.0, 500.);
    mehSiPixelCol[amend]->setAxisTitle("Column Number", 1);
    mehSiPixelCol[amend]->setAxisTitle("Count", 2);
  }

  // Muons
  iBooker.setCurrentFolder("GlobalDigisV/Muons");

  // DT
  std::string MuonString[4] = {"MB1", "MB2", "MB3", "MB4"};

  for (int i = 0; i < 4; ++i) {
    mehDtMuonn[i] = nullptr;
    mehDtMuonLayer[i] = nullptr;
    mehDtMuonTime[i] = nullptr;
    mehDtMuonTimevLayer[i] = nullptr;
  }

  for (int j = 0; j < 4; ++j) {
    hcharname = "hDtMuonn_" + MuonString[j];
    hchartitle = MuonString[j] + "  digis";
    mehDtMuonn[j] = iBooker.book1D(hcharname, hchartitle, 250, 0., 500.);
    mehDtMuonn[j]->setAxisTitle("Number of Digis", 1);
    mehDtMuonn[j]->setAxisTitle("Count", 2);

    hcharname = "hDtLayer_" + MuonString[j];
    hchartitle = MuonString[j] + "  Layer";
    mehDtMuonLayer[j] = iBooker.book1D(hcharname, hchartitle, 12, 1., 13.);
    mehDtMuonLayer[j]->setAxisTitle("4 * (SuperLayer - 1) + Layer", 1);
    mehDtMuonLayer[j]->setAxisTitle("Count", 2);

    hcharname = "hDtMuonTime_" + MuonString[j];
    hchartitle = MuonString[j] + "  Time";
    mehDtMuonTime[j] = iBooker.book1D(hcharname, hchartitle, 300, 400., 1000.);
    mehDtMuonTime[j]->setAxisTitle("Time", 1);
    mehDtMuonTime[j]->setAxisTitle("Count", 2);

    hcharname = "hDtMuonTimevLayer_" + MuonString[j];
    hchartitle = MuonString[j] + "  Time vs. Layer";
    mehDtMuonTimevLayer[j] = iBooker.bookProfile(hcharname, hchartitle, 12, 1., 13., 300, 400., 1000., "");
    mehDtMuonTimevLayer[j]->setAxisTitle("4 * (SuperLayer - 1) + Layer", 1);
    mehDtMuonTimevLayer[j]->setAxisTitle("Time", 2);
  }

  // CSC
  mehCSCStripn = nullptr;
  hcharname = "hCSCStripn";
  hchartitle = "CSC Strip digis";
  mehCSCStripn = iBooker.book1D(hcharname, hchartitle, 250, 0., 500.);
  mehCSCStripn->setAxisTitle("Number of Digis", 1);
  mehCSCStripn->setAxisTitle("Count", 2);

  mehCSCStripADC = nullptr;
  hcharname = "hCSCStripADC";
  hchartitle = "CSC Strip ADC";
  mehCSCStripADC = iBooker.book1D(hcharname, hchartitle, 110, 0., 1100.);
  mehCSCStripADC->setAxisTitle("ADC", 1);
  mehCSCStripADC->setAxisTitle("Count", 2);

  mehCSCWiren = nullptr;
  hcharname = "hCSCWiren";
  hchartitle = "CSC Wire digis";
  mehCSCWiren = iBooker.book1D(hcharname, hchartitle, 250, 0., 500.);
  mehCSCWiren->setAxisTitle("Number of Digis", 1);
  mehCSCWiren->setAxisTitle("Count", 2);

  mehCSCWireTime = nullptr;
  hcharname = "hCSCWireTime";
  hchartitle = "CSC Wire Time";
  mehCSCWireTime = iBooker.book1D(hcharname, hchartitle, 10, 0., 10.);
  mehCSCWireTime->setAxisTitle("Time", 1);
  mehCSCWireTime->setAxisTitle("Count", 2);

  // RPC
  mehRPCMuonn = nullptr;
  hcharname = "hRPCMuonn";
  hchartitle = "RPC digis";
  mehCSCStripn = iBooker.book1D(hcharname, hchartitle, 250, 0., 500.);
  mehCSCStripn->setAxisTitle("Number of Digis", 1);
  mehCSCStripn->setAxisTitle("Count", 2);

  std::string MuonRPCString[5] = {"Wmin2", "Wmin1", "W0", "Wpu1", "Wpu2"};
  for (int i = 0; i < 5; ++i) {
    mehRPCRes[i] = nullptr;
  }

  for (int j = 0; j < 5; ++j) {
    hcharname = "hRPCRes_" + MuonRPCString[j];
    hchartitle = MuonRPCString[j] + "  Digi - Sim";
    mehRPCRes[j] = iBooker.book1D(hcharname, hchartitle, 200, -8., 8.);
    mehRPCRes[j]->setAxisTitle("Digi - Sim center of strip x", 1);
    mehRPCRes[j]->setAxisTitle("Count", 2);
  }
}

void GlobalDigisAnalyzer::fillECal ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

private

Definition at line 506 of file GlobalDigisAnalyzer.cc.

References ESSample::adc(), EcalMGPASample::adc(), HLT_FULL_cff::barrelHits, EBHits_Token_, ECalbarrelADCtoGeV_, ECalEBSrc_Token_, ECalEESrc_Token_, ECalESSrc_Token_, ECalgainConv_, EEHits_Token_, HLT_FULL_cff::endcapHits, ESHits_Token_, dqm::impl::MonitorElement::Fill(), EcalMGPASample::gainId(), edm::Event::getByToken(), mps_fire::i, EEDataFrame::id(), EBDataFrame::id(), PixelPluginsPhase0_cfi::isBarrel, edm::HandleBase::isValid(), LogDebug, ESDataFrame::MAXSAMPLES, EcalDataFrame::MAXSAMPLES, mehEcalAEE, mehEcalMaxPos, mehEcalMultvAEE, mehEcaln, mehEcalSHE, mehEcalSHEvAEESHE, mehEScalADC, mehEScaln, ootPhotons_cff::preshowerHits, edm::Handle< T >::product(), DetId::rawId(), makeMEIFBenchmarkPlots::sample, ESDataFrame::size(), EcalDataFrame::size(), AlCaHLTBitMon_QueryRunRegistry::string, and verbosity.

Referenced by analyze().

                                                                                      {
  std::string MsgLoggerCat = "GlobalDigisAnalyzer_fillECal";

  TString eventout;
  if (verbosity > 0)
    eventout = "\nGathering info:";

  // extract crossing frame from event
  edm::Handle<CrossingFrame<PCaloHit>> crossingFrame;

  // extract EB information
  bool isBarrel = true;
  edm::Handle<EBDigiCollection> EcalDigiEB;
  iEvent.getByToken(ECalEBSrc_Token_, EcalDigiEB);
  bool validDigiEB = true;
  if (!EcalDigiEB.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find EcalDigiEB in event!";
    validDigiEB = false;
  }
  if (validDigiEB) {
    if (EcalDigiEB->empty())
      isBarrel = false;

    if (isBarrel) {
      // loop over simhits
      MapType ebSimMap;
      iEvent.getByToken(EBHits_Token_, crossingFrame);
      bool validXFrame = true;
      if (!crossingFrame.isValid()) {
        LogDebug(MsgLoggerCat) << "Unable to find cal barrel crossingFrame in event!";
        validXFrame = false;
      }
      if (validXFrame) {
        const MixCollection<PCaloHit> barrelHits(crossingFrame.product());

        // keep track of sum of simhit energy in each crystal
        for (auto const &iHit : barrelHits) {
          EBDetId ebid = EBDetId(iHit.id());

          uint32_t crystid = ebid.rawId();
          ebSimMap[crystid] += iHit.energy();
        }
      }

      // loop over digis
      const EBDigiCollection *barrelDigi = EcalDigiEB.product();

      std::vector<double> ebAnalogSignal;
      std::vector<double> ebADCCounts;
      std::vector<double> ebADCGains;
      ebAnalogSignal.reserve(EBDataFrame::MAXSAMPLES);
      ebADCCounts.reserve(EBDataFrame::MAXSAMPLES);
      ebADCGains.reserve(EBDataFrame::MAXSAMPLES);

      int i = 0;
      for (unsigned int digis = 0; digis < EcalDigiEB->size(); ++digis) {
        ++i;

        EBDataFrame ebdf = (*barrelDigi)[digis];
        int nrSamples = ebdf.size();

        EBDetId ebid = ebdf.id();

        double Emax = 0;
        int Pmax = 0;
        double pedestalPreSample = 0.;
        double pedestalPreSampleAnalog = 0.;

        for (int sample = 0; sample < nrSamples; ++sample) {
          ebAnalogSignal[sample] = 0.;
          ebADCCounts[sample] = 0.;
          ebADCGains[sample] = -1.;
        }

        // calculate maximum energy and pedestal
        for (int sample = 0; sample < nrSamples; ++sample) {
          EcalMGPASample thisSample = ebdf[sample];
          ebADCCounts[sample] = (thisSample.adc());
          ebADCGains[sample] = (thisSample.gainId());
          ebAnalogSignal[sample] = (ebADCCounts[sample] * ECalgainConv_[(int)ebADCGains[sample]] * ECalbarrelADCtoGeV_);
          if (Emax < ebAnalogSignal[sample]) {
            Emax = ebAnalogSignal[sample];
            Pmax = sample;
          }
          if (sample < 3) {
            pedestalPreSample += ebADCCounts[sample];
            pedestalPreSampleAnalog +=
                ebADCCounts[sample] * ECalgainConv_[(int)ebADCGains[sample]] * ECalbarrelADCtoGeV_;
          }
        }
        pedestalPreSample /= 3.;
        pedestalPreSampleAnalog /= 3.;

        // calculate pedestal subtracted digi energy in the crystal
        double Erec = Emax - pedestalPreSampleAnalog * ECalgainConv_[(int)ebADCGains[Pmax]];

        // gather necessary information
        mehEcalMaxPos[0]->Fill(Pmax);
        mehEcalSHE[0]->Fill(ebSimMap[ebid.rawId()]);
        mehEcalAEE[0]->Fill(Erec);
        // Adding protection against FPE
        if (ebSimMap[ebid.rawId()] != 0) {
          mehEcalSHEvAEESHE[0]->Fill(Erec / ebSimMap[ebid.rawId()], ebSimMap[ebid.rawId()]);
        }
        // else {
        //  std::cout<<"Would have been an FPE! with
        //  ebSimMap[ebid.rawId()]==0\n";
        //}

        mehEcalMultvAEE[0]->Fill(Pmax, (float)i);
      }

      if (verbosity > 1) {
        eventout += "\n          Number of EBDigis collected:.............. ";
        eventout += i;
      }
      mehEcaln[0]->Fill((float)i);
    }
  }

  // extract EE information
  bool isEndCap = true;
  edm::Handle<EEDigiCollection> EcalDigiEE;
  iEvent.getByToken(ECalEESrc_Token_, EcalDigiEE);
  bool validDigiEE = true;
  if (!EcalDigiEE.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find EcalDigiEE in event!";
    validDigiEE = false;
  }
  if (validDigiEE) {
    if (EcalDigiEE->empty())
      isEndCap = false;

    if (isEndCap) {
      // loop over simhits
      MapType eeSimMap;
      iEvent.getByToken(EEHits_Token_, crossingFrame);
      bool validXFrame = true;
      if (!crossingFrame.isValid()) {
        LogDebug(MsgLoggerCat) << "Unable to find cal endcap crossingFrame in event!";
        validXFrame = false;
      }
      if (validXFrame) {
        const MixCollection<PCaloHit> endcapHits(crossingFrame.product());

        // keep track of sum of simhit energy in each crystal
        for (auto const &iHit : endcapHits) {
          EEDetId eeid = EEDetId(iHit.id());

          uint32_t crystid = eeid.rawId();
          eeSimMap[crystid] += iHit.energy();
        }
      }

      // loop over digis
      const EEDigiCollection *endcapDigi = EcalDigiEE.product();

      std::vector<double> eeAnalogSignal;
      std::vector<double> eeADCCounts;
      std::vector<double> eeADCGains;
      eeAnalogSignal.reserve(EEDataFrame::MAXSAMPLES);
      eeADCCounts.reserve(EEDataFrame::MAXSAMPLES);
      eeADCGains.reserve(EEDataFrame::MAXSAMPLES);

      int inc = 0;
      for (unsigned int digis = 0; digis < EcalDigiEE->size(); ++digis) {
        ++inc;

        EEDataFrame eedf = (*endcapDigi)[digis];
        int nrSamples = eedf.size();

        EEDetId eeid = eedf.id();

        double Emax = 0;
        int Pmax = 0;
        double pedestalPreSample = 0.;
        double pedestalPreSampleAnalog = 0.;

        for (int sample = 0; sample < nrSamples; ++sample) {
          eeAnalogSignal[sample] = 0.;
          eeADCCounts[sample] = 0.;
          eeADCGains[sample] = -1.;
        }

        // calculate maximum enery and pedestal
        for (int sample = 0; sample < nrSamples; ++sample) {
          EcalMGPASample thisSample = eedf[sample];

          eeADCCounts[sample] = (thisSample.adc());
          eeADCGains[sample] = (thisSample.gainId());
          eeAnalogSignal[sample] = (eeADCCounts[sample] * ECalgainConv_[(int)eeADCGains[sample]] * ECalbarrelADCtoGeV_);
          if (Emax < eeAnalogSignal[sample]) {
            Emax = eeAnalogSignal[sample];
            Pmax = sample;
          }
          if (sample < 3) {
            pedestalPreSample += eeADCCounts[sample];
            pedestalPreSampleAnalog +=
                eeADCCounts[sample] * ECalgainConv_[(int)eeADCGains[sample]] * ECalbarrelADCtoGeV_;
          }
        }
        pedestalPreSample /= 3.;
        pedestalPreSampleAnalog /= 3.;

        // calculate pedestal subtracted digi energy in the crystal
        double Erec = Emax - pedestalPreSampleAnalog * ECalgainConv_[(int)eeADCGains[Pmax]];

        // gather necessary information
        mehEcalMaxPos[1]->Fill(Pmax);
        mehEcalSHE[1]->Fill(eeSimMap[eeid.rawId()]);
        mehEcalAEE[1]->Fill(Erec);
        // Adding protection against FPE
        if (eeSimMap[eeid.rawId()] != 0) {
          mehEcalSHEvAEESHE[1]->Fill(Erec / eeSimMap[eeid.rawId()], eeSimMap[eeid.rawId()]);
        }
        // else{
        //  std::cout<<"Would have been an FPE! with
        //  eeSimMap[eeid.rawId()]==0\n";
        //}
        mehEcalMultvAEE[1]->Fill(Pmax, (float)inc);
      }

      if (verbosity > 1) {
        eventout += "\n          Number of EEDigis collected:.............. ";
        eventout += inc;
      }

      mehEcaln[1]->Fill((float)inc);
    }
  }

  // extract ES information
  bool isPreshower = true;
  edm::Handle<ESDigiCollection> EcalDigiES;
  iEvent.getByToken(ECalESSrc_Token_, EcalDigiES);
  bool validDigiES = true;
  if (!EcalDigiES.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find EcalDigiES in event!";
    validDigiES = false;
  }

  // ONLY WHILE GEOMETRY IS REMOVED
  validDigiES = false;

  if (validDigiES) {
    if (EcalDigiES->empty())
      isPreshower = false;

    if (isPreshower) {
      // loop over simhits
      iEvent.getByToken(ESHits_Token_, crossingFrame);
      bool validXFrame = true;
      if (!crossingFrame.isValid()) {
        LogDebug(MsgLoggerCat) << "Unable to find cal preshower crossingFrame in event!";
        validXFrame = false;
      }
      if (validXFrame) {
        const MixCollection<PCaloHit> preshowerHits(crossingFrame.product());

        // keep track of sum of simhit energy in each crystal
        MapType esSimMap;
        for (auto const &iHit : preshowerHits) {
          ESDetId esid = ESDetId(iHit.id());

          uint32_t crystid = esid.rawId();
          esSimMap[crystid] += iHit.energy();
        }
      }

      // loop over digis
      const ESDigiCollection *preshowerDigi = EcalDigiES.product();

      std::vector<double> esADCCounts;
      esADCCounts.reserve(ESDataFrame::MAXSAMPLES);

      int i = 0;
      for (unsigned int digis = 0; digis < EcalDigiES->size(); ++digis) {
        ++i;

        ESDataFrame esdf = (*preshowerDigi)[digis];
        int nrSamples = esdf.size();

        for (int sample = 0; sample < nrSamples; ++sample) {
          esADCCounts[sample] = 0.;
        }

        // gether ADC counts
        for (int sample = 0; sample < nrSamples; ++sample) {
          ESSample thisSample = esdf[sample];
          esADCCounts[sample] = (thisSample.adc());
        }

        mehEScalADC[0]->Fill(esADCCounts[0]);
        mehEScalADC[1]->Fill(esADCCounts[1]);
        mehEScalADC[2]->Fill(esADCCounts[2]);
      }

      if (verbosity > 1) {
        eventout += "\n          Number of ESDigis collected:.............. ";
        eventout += i;
      }

      mehEScaln->Fill((float)i);
    }
  }
  if (verbosity > 0)
    edm::LogInfo(MsgLoggerCat) << eventout << "\n";

  return;
}

void GlobalDigisAnalyzer::fillHCal ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

private

Definition at line 823 of file GlobalDigisAnalyzer.cc.

References HcalCoderDb::adc2fC(), dqm::impl::MonitorElement::Fill(), edm::Event::getByToken(), edm::EventSetup::getHandle(), HBHEDigi_Token_, hcaldbToken_, HCalSrc_Token_, HFDigi_Token_, hcalSimParameters_cfi::ho, HODigi_Token_, ihf, cuy::ii, edm::HandleBase::isValid(), LogDebug, mehHcalAEE, mehHcalAEESHE, mehHcaln, mehHcalSHE, mehHcalSHEvAEE, HcalCalibrations::pedestal(), edm::Handle< T >::product(), DetId::rawId(), sdHcalBrl, sdHcalEC, sdHcalFwd, sdHcalOut, CaloSamples::size(), AlCaHLTBitMon_QueryRunRegistry::string, and verbosity.

Referenced by analyze().

                                                                                      {
  std::string MsgLoggerCat = "GlobalDigisAnalyzer_fillHCal";

  TString eventout;
  if (verbosity > 0)
    eventout = "\nGathering info:";

  // get calibration info
  const auto &HCalconditions = iSetup.getHandle(hcaldbToken_);
  if (!HCalconditions.isValid()) {
    edm::LogWarning(MsgLoggerCat) << "Unable to find HCalconditions in event!";
    return;
  }
  // HcalCalibrations calibrations;
  CaloSamples tool;

  // extract simhit info
  edm::Handle<edm::PCaloHitContainer> hcalHits;
  iEvent.getByToken(HCalSrc_Token_, hcalHits);
  bool validhcalHits = true;
  if (!hcalHits.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find hcalHits in event!";
    validhcalHits = false;
  }
  MapType fHBEnergySimHits;
  MapType fHEEnergySimHits;
  MapType fHOEnergySimHits;
  MapType fHFEnergySimHits;
  if (validhcalHits) {
    const edm::PCaloHitContainer *simhitResult = hcalHits.product();

    for (std::vector<PCaloHit>::const_iterator simhits = simhitResult->begin(); simhits != simhitResult->end();
         ++simhits) {
      HcalDetId detId(simhits->id());
      uint32_t cellid = detId.rawId();

      if (detId.subdet() == sdHcalBrl) {
        fHBEnergySimHits[cellid] += simhits->energy();
      }
      if (detId.subdet() == sdHcalEC) {
        fHEEnergySimHits[cellid] += simhits->energy();
      }
      if (detId.subdet() == sdHcalOut) {
        fHOEnergySimHits[cellid] += simhits->energy();
      }
      if (detId.subdet() == sdHcalFwd) {
        fHFEnergySimHits[cellid] += simhits->energy();
      }
    }
  }

  // get HBHE information
  edm::Handle<edm::SortedCollection<HBHEDataFrame>> hbhe;
  iEvent.getByToken(HBHEDigi_Token_, hbhe);
  bool validHBHE = true;
  if (!hbhe.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find HBHEDataFrame in event!";
    validHBHE = false;
  }

  if (validHBHE) {
    edm::SortedCollection<HBHEDataFrame>::const_iterator ihbhe;

    int iHB = 0;
    int iHE = 0;
    for (ihbhe = hbhe->begin(); ihbhe != hbhe->end(); ++ihbhe) {
      HcalDetId cell(ihbhe->id());

      if ((cell.subdet() == sdHcalBrl) || (cell.subdet() == sdHcalEC)) {
        // HCalconditions->makeHcalCalibration(cell, &calibrations);
        const HcalCalibrations &calibrations = HCalconditions->getHcalCalibrations(cell);
        const HcalQIECoder *channelCoder = HCalconditions->getHcalCoder(cell);
        const HcalQIEShape *shape = HCalconditions->getHcalShape(channelCoder);
        HcalCoderDb coder(*channelCoder, *shape);
        coder.adc2fC(*ihbhe, tool);

        // get HB info
        if (cell.subdet() == sdHcalBrl) {
          ++iHB;
          float fDigiSum = 0.0;
          for (int ii = 0; ii < tool.size(); ++ii) {
            // default ped is 4.5
            int capid = (*ihbhe)[ii].capid();
            fDigiSum += (tool[ii] - calibrations.pedestal(capid));
          }

          mehHcalSHE[0]->Fill(fHFEnergySimHits[cell.rawId()]);
          mehHcalAEE[0]->Fill(fDigiSum);
          // Adding protection against FPE
          if (fHFEnergySimHits[cell.rawId()] != 0) {
            mehHcalAEESHE[0]->Fill(fDigiSum / fHFEnergySimHits[cell.rawId()]);
          }
          // else {
          //  std::cout<<"It would have been an FPE! with
          //  fHFEnergySimHits[cell.rawId()]==0!\n";
          //}
          mehHcalSHEvAEE[0]->Fill(fDigiSum, fHFEnergySimHits[cell.rawId()]);
        }

        // get HE info
        if (cell.subdet() == sdHcalEC) {
          ++iHE;
          float fDigiSum = 0.0;
          for (int ii = 0; ii < tool.size(); ++ii) {
            int capid = (*ihbhe)[ii].capid();
            fDigiSum += (tool[ii] - calibrations.pedestal(capid));
          }

          mehHcalSHE[1]->Fill(fHFEnergySimHits[cell.rawId()]);
          mehHcalAEE[1]->Fill(fDigiSum);
          // Adding protection against FPE
          if (fHFEnergySimHits[cell.rawId()] != 0) {
            mehHcalAEESHE[1]->Fill(fDigiSum / fHFEnergySimHits[cell.rawId()]);
          }
          // else{
          //  std::cout<<"It would have been an FPE! with
          //  fHFEnergySimHits[cell.rawId()]==0!\n";
          //}
          mehHcalSHEvAEE[1]->Fill(fDigiSum, fHFEnergySimHits[cell.rawId()]);
        }
      }
    }

    if (verbosity > 1) {
      eventout += "\n          Number of HBDigis collected:.............. ";
      eventout += iHB;
    }
    mehHcaln[0]->Fill((float)iHB);

    if (verbosity > 1) {
      eventout += "\n          Number of HEDigis collected:.............. ";
      eventout += iHE;
    }
    mehHcaln[1]->Fill((float)iHE);
  }

  // get HO information
  edm::Handle<edm::SortedCollection<HODataFrame>> ho;
  iEvent.getByToken(HODigi_Token_, ho);
  bool validHO = true;
  if (!ho.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find HODataFrame in event!";
    validHO = false;
  }
  if (validHO) {
    edm::SortedCollection<HODataFrame>::const_iterator iho;

    int iHO = 0;
    for (iho = ho->begin(); iho != ho->end(); ++iho) {
      HcalDetId cell(iho->id());

      if (cell.subdet() == sdHcalOut) {
        // HCalconditions->makeHcalCalibration(cell, &calibrations);
        const HcalCalibrations &calibrations = HCalconditions->getHcalCalibrations(cell);
        const HcalQIECoder *channelCoder = HCalconditions->getHcalCoder(cell);
        const HcalQIEShape *shape = HCalconditions->getHcalShape(channelCoder);
        HcalCoderDb coder(*channelCoder, *shape);
        coder.adc2fC(*iho, tool);

        ++iHO;
        float fDigiSum = 0.0;
        for (int ii = 0; ii < tool.size(); ++ii) {
          // default ped is 4.5
          int capid = (*iho)[ii].capid();
          fDigiSum += (tool[ii] - calibrations.pedestal(capid));
        }

        mehHcalSHE[2]->Fill(fHFEnergySimHits[cell.rawId()]);
        mehHcalAEE[2]->Fill(fDigiSum);
        // Adding protection against FPE
        if (fHFEnergySimHits[cell.rawId()] != 0) {
          mehHcalAEESHE[2]->Fill(fDigiSum / fHFEnergySimHits[cell.rawId()]);
        }
        // else{
        //  std::cout<<"It would have been an FPE! with
        //  fHFEnergySimHits[cell.rawId()]==0!\n";
        //}
        mehHcalSHEvAEE[2]->Fill(fDigiSum, fHFEnergySimHits[cell.rawId()]);
      }
    }

    if (verbosity > 1) {
      eventout += "\n          Number of HODigis collected:.............. ";
      eventout += iHO;
    }
    mehHcaln[2]->Fill((float)iHO);
  }

  // get HF information
  edm::Handle<edm::SortedCollection<HFDataFrame>> hf;
  iEvent.getByToken(HFDigi_Token_, hf);
  bool validHF = true;
  if (!hf.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find HFDataFrame in event!";
    validHF = false;
  }
  if (validHF) {
    edm::SortedCollection<HFDataFrame>::const_iterator ihf;

    int iHF = 0;
    for (ihf = hf->begin(); ihf != hf->end(); ++ihf) {
      HcalDetId cell(ihf->id());

      if (cell.subdet() == sdHcalFwd) {
        // HCalconditions->makeHcalCalibration(cell, &calibrations);
        const HcalCalibrations &calibrations = HCalconditions->getHcalCalibrations(cell);
        const HcalQIECoder *channelCoder = HCalconditions->getHcalCoder(cell);
        const HcalQIEShape *shape = HCalconditions->getHcalShape(channelCoder);
        HcalCoderDb coder(*channelCoder, *shape);
        coder.adc2fC(*ihf, tool);

        ++iHF;
        float fDigiSum = 0.0;
        for (int ii = 0; ii < tool.size(); ++ii) {
          // default ped is 1.73077
          int capid = (*ihf)[ii].capid();
          fDigiSum += (tool[ii] - calibrations.pedestal(capid));
        }

        mehHcalSHE[3]->Fill(fHFEnergySimHits[cell.rawId()]);
        mehHcalAEE[3]->Fill(fDigiSum);
        // Adding protection against FPE
        if (fHFEnergySimHits[cell.rawId()] != 0) {
          mehHcalAEESHE[3]->Fill(fDigiSum / fHFEnergySimHits[cell.rawId()]);
        }
        // else{
        //  std::cout<<"It would have been an FPE! with
        //  fHFEnergySimHits[cell.rawId()]==0!\n";
        //}
        mehHcalSHEvAEE[3]->Fill(fDigiSum, fHFEnergySimHits[cell.rawId()]);
      }
    }

    if (verbosity > 1) {
      eventout += "\n          Number of HFDigis collected:.............. ";
      eventout += iHF;
    }
    mehHcaln[3]->Fill((float)iHF);
  }

  if (verbosity > 0)
    edm::LogInfo(MsgLoggerCat) << eventout << "\n";

  return;
}

void GlobalDigisAnalyzer::fillMuon ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

private

Definition at line 1322 of file GlobalDigisAnalyzer.cc.

References RPCRoll::centreOfStrip(), dqm::impl::MonitorElement::Fill(), HcalObjRepresent::Fill(), edm::Event::getByToken(), edm::EventSetup::getHandle(), edm::HandleBase::isValid(), dqmiolumiharvest::j, dqmdumpme::last, phase1PixelTopology::layer, LogDebug, mehCSCStripADC, mehCSCStripn, mehCSCWiren, mehCSCWireTime, mehDtMuonLayer, mehDtMuonn, mehDtMuonTime, mehDtMuonTimevLayer, mehRPCMuonn, mehRPCRes, MuCSCStripSrc_Token_, MuCSCWireSrc_Token_, MuDTSrc_Token_, MuRPCSrc_Token_, me0TriggerPseudoDigis_cff::nStrips, EcalCondDBWriter_cfi::pedestal, sistrip::SpyUtilities::range(), RPCDetId::region(), RPCDetId::ring(), RPCDetId, rpcGeomToken_, RPCSimHit_Token_, AlCaHLTBitMon_QueryRunRegistry::string, DigiDM_cff::strips, theCSCStripPedestalCount, theCSCStripPedestalSum, verbosity, DigiDM_cff::wires, and x.

Referenced by analyze().

                                                                                      {
  std::string MsgLoggerCat = "GlobalDigisAnalyzer_fillMuon";

  TString eventout;
  if (verbosity > 0)
    eventout = "\nGathering info:";

  // get DT information
  edm::Handle<DTDigiCollection> dtDigis;
  iEvent.getByToken(MuDTSrc_Token_, dtDigis);
  bool validdtDigis = true;
  if (!dtDigis.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find dtDigis in event!";
    validdtDigis = false;
  }
  if (validdtDigis) {
    int nDt0 = 0;
    int nDt1 = 0;
    int nDt2 = 0;
    int nDt3 = 0;
    int nDt = 0;
    DTDigiCollection::DigiRangeIterator detUnitIt;
    for (detUnitIt = dtDigis->begin(); detUnitIt != dtDigis->end(); ++detUnitIt) {
      const DTLayerId &id = (*detUnitIt).first;
      const DTDigiCollection::Range &range = (*detUnitIt).second;
      for (DTDigiCollection::const_iterator digiIt = range.first; digiIt != range.second; ++digiIt) {
        ++nDt;

        DTWireId wireId(id, (*digiIt).wire());
        if (wireId.station() == 1) {
          mehDtMuonLayer[0]->Fill(id.layer());
          mehDtMuonTime[0]->Fill((*digiIt).time());
          mehDtMuonTimevLayer[0]->Fill(id.layer(), (*digiIt).time());
          ++nDt0;
        }
        if (wireId.station() == 2) {
          mehDtMuonLayer[1]->Fill(id.layer());
          mehDtMuonTime[1]->Fill((*digiIt).time());
          mehDtMuonTimevLayer[1]->Fill(id.layer(), (*digiIt).time());
          ++nDt1;
        }
        if (wireId.station() == 3) {
          mehDtMuonLayer[2]->Fill(id.layer());
          mehDtMuonTime[2]->Fill((*digiIt).time());
          mehDtMuonTimevLayer[2]->Fill(id.layer(), (*digiIt).time());
          ++nDt2;
        }
        if (wireId.station() == 4) {
          mehDtMuonLayer[3]->Fill(id.layer());
          mehDtMuonTime[3]->Fill((*digiIt).time());
          mehDtMuonTimevLayer[3]->Fill(id.layer(), (*digiIt).time());
          ++nDt3;
        }
      }
    }
    mehDtMuonn[0]->Fill((float)nDt0);
    mehDtMuonn[1]->Fill((float)nDt1);
    mehDtMuonn[2]->Fill((float)nDt2);
    mehDtMuonn[3]->Fill((float)nDt3);

    if (verbosity > 1) {
      eventout += "\n          Number of DtMuonDigis collected:.......... ";
      eventout += nDt;
    }
  }

  // get CSC Strip information
  edm::Handle<CSCStripDigiCollection> strips;
  iEvent.getByToken(MuCSCStripSrc_Token_, strips);
  bool validstrips = true;
  if (!strips.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find muon strips in event!";
    validstrips = false;
  }

  if (validstrips) {
    int nStrips = 0;
    for (CSCStripDigiCollection::DigiRangeIterator j = strips->begin(); j != strips->end(); ++j) {
      std::vector<CSCStripDigi>::const_iterator digiItr = (*j).second.first;
      std::vector<CSCStripDigi>::const_iterator last = (*j).second.second;

      for (; digiItr != last; ++digiItr) {
        ++nStrips;

        // average pedestals
        std::vector<int> adcCounts = digiItr->getADCCounts();
        theCSCStripPedestalSum += adcCounts[0];
        theCSCStripPedestalSum += adcCounts[1];
        theCSCStripPedestalCount += 2;

        // if there are enough pedestal statistics
        if (theCSCStripPedestalCount > 100) {
          float pedestal = theCSCStripPedestalSum / theCSCStripPedestalCount;
          if (adcCounts[5] > (pedestal + 100))
            mehCSCStripADC->Fill(adcCounts[4] - pedestal);
        }
      }
    }

    if (verbosity > 1) {
      eventout += "\n          Number of CSCStripDigis collected:........ ";
      eventout += nStrips;
    }
    mehCSCStripn->Fill((float)nStrips);
  }

  // get CSC Wire information
  edm::Handle<CSCWireDigiCollection> wires;
  iEvent.getByToken(MuCSCWireSrc_Token_, wires);
  bool validwires = true;
  if (!wires.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find muon wires in event!";
    validwires = false;
  }

  if (validwires) {
    int nWires = 0;
    for (CSCWireDigiCollection::DigiRangeIterator j = wires->begin(); j != wires->end(); ++j) {
      std::vector<CSCWireDigi>::const_iterator digiItr = (*j).second.first;
      std::vector<CSCWireDigi>::const_iterator endDigi = (*j).second.second;

      for (; digiItr != endDigi; ++digiItr) {
        ++nWires;
        mehCSCWireTime->Fill(digiItr->getTimeBin());
      }
    }

    if (verbosity > 1) {
      eventout += "\n          Number of CSCWireDigis collected:......... ";
      eventout += nWires;
    }
    mehCSCWiren->Fill((float)nWires);
  }

  // get RPC information
  // Get the RPC Geometry
  const auto &rpcGeom = iSetup.getHandle(rpcGeomToken_);
  if (!rpcGeom.isValid()) {
    edm::LogWarning(MsgLoggerCat) << "Unable to find RPCGeometryRecord in event!";
    return;
  }

  edm::Handle<edm::PSimHitContainer> rpcsimHit;
  iEvent.getByToken(RPCSimHit_Token_, rpcsimHit);
  bool validrpcsim = true;
  if (!rpcsimHit.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find rpcsimHit in event!";
    validrpcsim = false;
  }

  edm::Handle<RPCDigiCollection> rpcDigis;
  iEvent.getByToken(MuRPCSrc_Token_, rpcDigis);
  bool validrpcdigi = true;
  if (!rpcDigis.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find rpcDigis in event!";
    validrpcdigi = false;
  }

  // ONLY UNTIL PROBLEM WITH RPC DIGIS IS FIGURED OUT
  validrpcdigi = false;

  // Loop on simhits
  edm::PSimHitContainer::const_iterator rpcsimIt;
  std::map<RPCDetId, std::vector<double>> allsims;

  if (validrpcsim) {
    for (rpcsimIt = rpcsimHit->begin(); rpcsimIt != rpcsimHit->end(); rpcsimIt++) {
      RPCDetId Rsid = (RPCDetId)(*rpcsimIt).detUnitId();
      int ptype = rpcsimIt->particleType();

      if (ptype == 13 || ptype == -13) {
        std::vector<double> buff;
        if (allsims.find(Rsid) != allsims.end()) {
          buff = allsims[Rsid];
        }
        buff.push_back(rpcsimIt->localPosition().x());
        allsims[Rsid] = buff;
      }
    }
  }

  // CRASH HAPPENS SOMEWHERE HERE IN FOR DECLARATION
  // WHAT IS WRONG WITH rpcDigis?????
  if (validrpcdigi) {
    int nRPC = 0;
    RPCDigiCollection::DigiRangeIterator rpcdetUnitIt;
    for (rpcdetUnitIt = rpcDigis->begin(); rpcdetUnitIt != rpcDigis->end(); ++rpcdetUnitIt) {
      const RPCDetId Rsid = (*rpcdetUnitIt).first;
      const RPCRoll *roll = dynamic_cast<const RPCRoll *>(rpcGeom->roll(Rsid));
      const RPCDigiCollection::Range &range = (*rpcdetUnitIt).second;

      std::vector<double> sims;
      if (allsims.find(Rsid) != allsims.end()) {
        sims = allsims[Rsid];
      }

      int ndigi = 0;
      for (RPCDigiCollection::const_iterator rpcdigiIt = range.first; rpcdigiIt != range.second; ++rpcdigiIt) {
        ++ndigi;
        ++nRPC;
      }

      if (sims.size() == 1 && ndigi == 1) {
        double dis = roll->centreOfStrip(range.first->strip()).x() - sims[0];

        if (Rsid.region() == 0) {
          if (Rsid.ring() == -2)
            mehRPCRes[0]->Fill((float)dis);
          else if (Rsid.ring() == -1)
            mehRPCRes[1]->Fill((float)dis);
          else if (Rsid.ring() == 0)
            mehRPCRes[2]->Fill((float)dis);
          else if (Rsid.ring() == 1)
            mehRPCRes[3]->Fill((float)dis);
          else if (Rsid.ring() == 2)
            mehRPCRes[4]->Fill((float)dis);
        }
      }
    }

    if (verbosity > 1) {
      eventout += "\n          Number of RPCDigis collected:.............. ";
      eventout += nRPC;
    }
    mehRPCMuonn->Fill(float(nRPC));
  }

  if (verbosity > 0)
    edm::LogInfo(MsgLoggerCat) << eventout << "\n";

  return;
}

void GlobalDigisAnalyzer::fillTrk ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

private

Definition at line 1077 of file GlobalDigisAnalyzer.cc.

References SplitLinear::begin, edm::DetSetVector< T >::begin(), edm::DetSetVector< T >::end(), dataset::end, dqm::impl::MonitorElement::Fill(), edm::Event::getByToken(), edm::EventSetup::getData(), mps_fire::i, edm::HandleBase::isValid(), LogDebug, mehSiPixelADC, mehSiPixelCol, mehSiPixeln, mehSiPixelRow, mehSiStripADC, mehSiStripn, mehSiStripStrip, TrackerTopology::pxbLayer(), TrackerTopology::pxfDisk(), TrackerTopology::pxfSide(), sdPxlBrl, sdPxlFwd, sdSiTEC, sdSiTIB, sdSiTID, sdSiTOB, SiPxlSrc_Token_, SiStripSrc_Token_, AlCaHLTBitMon_QueryRunRegistry::string, DetId::subdetId(), TrackerTopology::tecWheel(), TrackerTopology::tibLayer(), TrackerTopology::tidWheel(), TrackerTopology::tobLayer(), tTopoToken_, and verbosity.

Referenced by analyze().

                                                                                     {
  // Retrieve tracker topology from geometry
  const TrackerTopology *const tTopo = &iSetup.getData(tTopoToken_);
  std::string MsgLoggerCat = "GlobalDigisAnalyzer_fillTrk";
  TString eventout;
  if (verbosity > 0)
    eventout = "\nGathering info:";

  // get strip information
  edm::Handle<edm::DetSetVector<SiStripDigi>> stripDigis;
  iEvent.getByToken(SiStripSrc_Token_, stripDigis);
  bool validstripDigis = true;
  if (!stripDigis.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find stripDigis in event!";
    validstripDigis = false;
  }

  if (validstripDigis) {
    int nStripBrl = 0, nStripFwd = 0;
    edm::DetSetVector<SiStripDigi>::const_iterator DSViter;
    for (DSViter = stripDigis->begin(); DSViter != stripDigis->end(); ++DSViter) {
      unsigned int id = DSViter->id;
      DetId detId(id);
      edm::DetSet<SiStripDigi>::const_iterator begin = DSViter->data.begin();
      edm::DetSet<SiStripDigi>::const_iterator end = DSViter->data.end();
      edm::DetSet<SiStripDigi>::const_iterator iter;

      // get TIB
      if (detId.subdetId() == sdSiTIB) {
        for (iter = begin; iter != end; ++iter) {
          ++nStripBrl;
          if (tTopo->tibLayer(id) == 1) {
            mehSiStripADC[0]->Fill((*iter).adc());
            mehSiStripStrip[0]->Fill((*iter).strip());
          }
          if (tTopo->tibLayer(id) == 2) {
            mehSiStripADC[1]->Fill((*iter).adc());
            mehSiStripStrip[1]->Fill((*iter).strip());
          }
          if (tTopo->tibLayer(id) == 3) {
            mehSiStripADC[2]->Fill((*iter).adc());
            mehSiStripStrip[2]->Fill((*iter).strip());
          }
          if (tTopo->tibLayer(id) == 4) {
            mehSiStripADC[3]->Fill((*iter).adc());
            mehSiStripStrip[3]->Fill((*iter).strip());
          }
        }
      }

      // get TOB
      if (detId.subdetId() == sdSiTOB) {
        for (iter = begin; iter != end; ++iter) {
          ++nStripBrl;
          if (tTopo->tobLayer(id) == 1) {
            mehSiStripADC[4]->Fill((*iter).adc());
            mehSiStripStrip[4]->Fill((*iter).strip());
          }
          if (tTopo->tobLayer(id) == 2) {
            mehSiStripADC[5]->Fill((*iter).adc());
            mehSiStripStrip[5]->Fill((*iter).strip());
          }
          if (tTopo->tobLayer(id) == 3) {
            mehSiStripADC[6]->Fill((*iter).adc());
            mehSiStripStrip[6]->Fill((*iter).strip());
          }
          if (tTopo->tobLayer(id) == 4) {
            mehSiStripADC[7]->Fill((*iter).adc());
            mehSiStripStrip[7]->Fill((*iter).strip());
          }
        }
      }

      // get TID
      if (detId.subdetId() == sdSiTID) {
        for (iter = begin; iter != end; ++iter) {
          ++nStripFwd;
          if (tTopo->tidWheel(id) == 1) {
            mehSiStripADC[8]->Fill((*iter).adc());
            mehSiStripStrip[8]->Fill((*iter).strip());
          }
          if (tTopo->tidWheel(id) == 2) {
            mehSiStripADC[9]->Fill((*iter).adc());
            mehSiStripStrip[9]->Fill((*iter).strip());
          }
          if (tTopo->tidWheel(id) == 3) {
            mehSiStripADC[10]->Fill((*iter).adc());
            mehSiStripStrip[10]->Fill((*iter).strip());
          }
        }
      }

      // get TEC
      if (detId.subdetId() == sdSiTEC) {
        for (iter = begin; iter != end; ++iter) {
          ++nStripFwd;
          if (tTopo->tecWheel(id) == 1) {
            mehSiStripADC[11]->Fill((*iter).adc());
            mehSiStripStrip[11]->Fill((*iter).strip());
          }
          if (tTopo->tecWheel(id) == 2) {
            mehSiStripADC[12]->Fill((*iter).adc());
            mehSiStripStrip[12]->Fill((*iter).strip());
          }
          if (tTopo->tecWheel(id) == 3) {
            mehSiStripADC[13]->Fill((*iter).adc());
            mehSiStripStrip[13]->Fill((*iter).strip());
          }
          if (tTopo->tecWheel(id) == 4) {
            mehSiStripADC[14]->Fill((*iter).adc());
            mehSiStripStrip[14]->Fill((*iter).strip());
          }
          if (tTopo->tecWheel(id) == 5) {
            mehSiStripADC[15]->Fill((*iter).adc());
            mehSiStripStrip[15]->Fill((*iter).strip());
          }
          if (tTopo->tecWheel(id) == 6) {
            mehSiStripADC[16]->Fill((*iter).adc());
            mehSiStripStrip[16]->Fill((*iter).strip());
          }
          if (tTopo->tecWheel(id) == 7) {
            mehSiStripADC[17]->Fill((*iter).adc());
            mehSiStripStrip[17]->Fill((*iter).strip());
          }
          if (tTopo->tecWheel(id) == 8) {
            mehSiStripADC[18]->Fill((*iter).adc());
            mehSiStripStrip[18]->Fill((*iter).strip());
          }
        }
      }
    }  // end loop over DataSetVector

    if (verbosity > 1) {
      eventout += "\n          Number of BrlStripDigis collected:........ ";
      eventout += nStripBrl;
    }
    for (int i = 0; i < 8; ++i) {
      mehSiStripn[i]->Fill((float)nStripBrl);
    }

    if (verbosity > 1) {
      eventout += "\n          Number of FrwdStripDigis collected:....... ";
      eventout += nStripFwd;
    }
    for (int i = 8; i < 19; ++i) {
      mehSiStripn[i]->Fill((float)nStripFwd);
    }
  }

  // get pixel information
  edm::Handle<edm::DetSetVector<PixelDigi>> pixelDigis;
  iEvent.getByToken(SiPxlSrc_Token_, pixelDigis);
  bool validpixelDigis = true;
  if (!pixelDigis.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find pixelDigis in event!";
    validpixelDigis = false;
  }
  if (validpixelDigis) {
    int nPxlBrl = 0, nPxlFwd = 0;
    edm::DetSetVector<PixelDigi>::const_iterator DPViter;
    for (DPViter = pixelDigis->begin(); DPViter != pixelDigis->end(); ++DPViter) {
      unsigned int id = DPViter->id;
      DetId detId(id);
      edm::DetSet<PixelDigi>::const_iterator begin = DPViter->data.begin();
      edm::DetSet<PixelDigi>::const_iterator end = DPViter->data.end();
      edm::DetSet<PixelDigi>::const_iterator iter;

      // get Barrel pixels
      if (detId.subdetId() == sdPxlBrl) {
        for (iter = begin; iter != end; ++iter) {
          ++nPxlBrl;
          if (tTopo->pxbLayer(id) == 1) {
            mehSiPixelADC[0]->Fill((*iter).adc());
            mehSiPixelRow[0]->Fill((*iter).row());
            mehSiPixelCol[0]->Fill((*iter).column());
          }
          if (tTopo->pxbLayer(id) == 2) {
            mehSiPixelADC[1]->Fill((*iter).adc());
            mehSiPixelRow[1]->Fill((*iter).row());
            mehSiPixelCol[1]->Fill((*iter).column());
          }
          if (tTopo->pxbLayer(id) == 3) {
            mehSiPixelADC[2]->Fill((*iter).adc());
            mehSiPixelRow[2]->Fill((*iter).row());
            mehSiPixelCol[2]->Fill((*iter).column());
          }
        }
      }

      // get Forward pixels
      if (detId.subdetId() == sdPxlFwd) {
        for (iter = begin; iter != end; ++iter) {
          ++nPxlFwd;
          if (tTopo->pxfDisk(id) == 1) {
            if (tTopo->pxfSide(id) == 1) {
              mehSiPixelADC[4]->Fill((*iter).adc());
              mehSiPixelRow[4]->Fill((*iter).row());
              mehSiPixelCol[4]->Fill((*iter).column());
            }
            if (tTopo->pxfSide(id) == 2) {
              mehSiPixelADC[3]->Fill((*iter).adc());
              mehSiPixelRow[3]->Fill((*iter).row());
              mehSiPixelCol[3]->Fill((*iter).column());
            }
          }
          if (tTopo->pxfDisk(id) == 2) {
            if (tTopo->pxfSide(id) == 1) {
              mehSiPixelADC[6]->Fill((*iter).adc());
              mehSiPixelRow[6]->Fill((*iter).row());
              mehSiPixelCol[6]->Fill((*iter).column());
            }
            if (tTopo->pxfSide(id) == 2) {
              mehSiPixelADC[5]->Fill((*iter).adc());
              mehSiPixelRow[5]->Fill((*iter).row());
              mehSiPixelCol[5]->Fill((*iter).column());
            }
          }
        }
      }
    }

    if (verbosity > 1) {
      eventout += "\n          Number of BrlPixelDigis collected:........ ";
      eventout += nPxlBrl;
    }
    for (int i = 0; i < 3; ++i) {
      mehSiPixeln[i]->Fill((float)nPxlBrl);
    }

    if (verbosity > 1) {
      eventout += "\n          Number of FrwdPixelDigis collected:....... ";
      eventout += nPxlFwd;
    }

    for (int i = 3; i < 7; ++i) {
      mehSiPixeln[i]->Fill((float)nPxlFwd);
    }
  }

  if (verbosity > 0)
    edm::LogInfo(MsgLoggerCat) << eventout << "\n";

  return;
}

Member Data Documentation

unsigned int GlobalDigisAnalyzer::count

private

Definition at line 245 of file GlobalDigisAnalyzer.h.

Referenced by analyze().

edm::EDGetTokenT<CrossingFrame<PCaloHit> > GlobalDigisAnalyzer::EBHits_Token_

private

Definition at line 239 of file GlobalDigisAnalyzer.h.

Referenced by fillECal(), and GlobalDigisAnalyzer().

edm::ESGetToken<EcalADCToGeVConstant, EcalADCToGeVConstantRcd> GlobalDigisAnalyzer::ecalADCtoGevToken_

private

Definition at line 158 of file GlobalDigisAnalyzer.h.

Referenced by analyze(), and GlobalDigisAnalyzer().

double GlobalDigisAnalyzer::ECalbarrelADCtoGeV_

private

Definition at line 164 of file GlobalDigisAnalyzer.h.

Referenced by analyze(), fillECal(), and GlobalDigisAnalyzer().

edm::InputTag GlobalDigisAnalyzer::ECalEBSrc_

private

Definition at line 155 of file GlobalDigisAnalyzer.h.

Referenced by GlobalDigisAnalyzer().

edm::EDGetTokenT<EBDigiCollection> GlobalDigisAnalyzer::ECalEBSrc_Token_

private

Definition at line 226 of file GlobalDigisAnalyzer.h.

Referenced by fillECal(), and GlobalDigisAnalyzer().

edm::InputTag GlobalDigisAnalyzer::ECalEESrc_

private

Definition at line 156 of file GlobalDigisAnalyzer.h.

Referenced by GlobalDigisAnalyzer().

edm::EDGetTokenT<EEDigiCollection> GlobalDigisAnalyzer::ECalEESrc_Token_

private

Definition at line 227 of file GlobalDigisAnalyzer.h.

Referenced by fillECal(), and GlobalDigisAnalyzer().

double GlobalDigisAnalyzer::ECalendcapADCtoGeV_

private

Definition at line 165 of file GlobalDigisAnalyzer.h.

Referenced by analyze(), and GlobalDigisAnalyzer().

edm::InputTag GlobalDigisAnalyzer::ECalESSrc_

private

Definition at line 157 of file GlobalDigisAnalyzer.h.

Referenced by GlobalDigisAnalyzer().

edm::EDGetTokenT<ESDigiCollection> GlobalDigisAnalyzer::ECalESSrc_Token_

private

Definition at line 228 of file GlobalDigisAnalyzer.h.

Referenced by fillECal(), and GlobalDigisAnalyzer().

std::map<int, double, std::less<int> > GlobalDigisAnalyzer::ECalgainConv_

private

Definition at line 163 of file GlobalDigisAnalyzer.h.

Referenced by fillECal(), and GlobalDigisAnalyzer().

edm::EDGetTokenT<CrossingFrame<PCaloHit> > GlobalDigisAnalyzer::EEHits_Token_

private

Definition at line 240 of file GlobalDigisAnalyzer.h.

Referenced by fillECal(), and GlobalDigisAnalyzer().

edm::EDGetTokenT<CrossingFrame<PCaloHit> > GlobalDigisAnalyzer::ESHits_Token_

private

Definition at line 241 of file GlobalDigisAnalyzer.h.

Referenced by fillECal(), and GlobalDigisAnalyzer().

std::string GlobalDigisAnalyzer::fName

private

Definition at line 135 of file GlobalDigisAnalyzer.h.

Referenced by GlobalDigisAnalyzer().

int GlobalDigisAnalyzer::frequency

private

Definition at line 137 of file GlobalDigisAnalyzer.h.

Referenced by analyze(), and GlobalDigisAnalyzer().

bool GlobalDigisAnalyzer::getAllProvenances

private

Definition at line 139 of file GlobalDigisAnalyzer.h.

Referenced by analyze(), and GlobalDigisAnalyzer().

edm::EDGetTokenT<edm::SortedCollection<HBHEDataFrame> > GlobalDigisAnalyzer::HBHEDigi_Token_

private

Definition at line 230 of file GlobalDigisAnalyzer.h.

Referenced by fillHCal(), and GlobalDigisAnalyzer().

edm::ESGetToken<HcalDbService, HcalDbRecord> GlobalDigisAnalyzer::hcaldbToken_

private

Definition at line 161 of file GlobalDigisAnalyzer.h.

Referenced by fillHCal(), and GlobalDigisAnalyzer().

edm::InputTag GlobalDigisAnalyzer::HCalDigi_

private

Definition at line 176 of file GlobalDigisAnalyzer.h.

Referenced by GlobalDigisAnalyzer().

edm::InputTag GlobalDigisAnalyzer::HCalSrc_

private

Definition at line 175 of file GlobalDigisAnalyzer.h.

Referenced by GlobalDigisAnalyzer().

edm::EDGetTokenT<edm::PCaloHitContainer> GlobalDigisAnalyzer::HCalSrc_Token_

private

Definition at line 229 of file GlobalDigisAnalyzer.h.

Referenced by fillHCal(), and GlobalDigisAnalyzer().

edm::EDGetTokenT<edm::SortedCollection<HFDataFrame> > GlobalDigisAnalyzer::HFDigi_Token_

private

Definition at line 232 of file GlobalDigisAnalyzer.h.

Referenced by fillHCal(), and GlobalDigisAnalyzer().

std::string GlobalDigisAnalyzer::hitsProducer

private

Definition at line 141 of file GlobalDigisAnalyzer.h.

Referenced by GlobalDigisAnalyzer().

edm::EDGetTokenT<edm::SortedCollection<HODataFrame> > GlobalDigisAnalyzer::HODigi_Token_

private

Definition at line 231 of file GlobalDigisAnalyzer.h.

Referenced by fillHCal(), and GlobalDigisAnalyzer().

std::string GlobalDigisAnalyzer::label

private

Definition at line 138 of file GlobalDigisAnalyzer.h.

MonitorElement* GlobalDigisAnalyzer::mehCSCStripADC

private

Definition at line 209 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillMuon().

MonitorElement* GlobalDigisAnalyzer::mehCSCStripn

private

Definition at line 208 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillMuon().

MonitorElement* GlobalDigisAnalyzer::mehCSCWiren

private

Definition at line 210 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillMuon().

MonitorElement* GlobalDigisAnalyzer::mehCSCWireTime

private

Definition at line 211 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillMuon().

MonitorElement* GlobalDigisAnalyzer::mehDtMuonLayer[4]

private

Definition at line 200 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillMuon().

MonitorElement* GlobalDigisAnalyzer::mehDtMuonn[4]

private

Definition at line 199 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillMuon().

MonitorElement* GlobalDigisAnalyzer::mehDtMuonTime[4]

private

Definition at line 201 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillMuon().

MonitorElement* GlobalDigisAnalyzer::mehDtMuonTimevLayer[4]

private

Definition at line 202 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillMuon().

MonitorElement* GlobalDigisAnalyzer::mehEcalAEE[2]

private

Definition at line 148 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillECal().

MonitorElement* GlobalDigisAnalyzer::mehEcalMaxPos[2]

private

Definition at line 150 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillECal().

MonitorElement* GlobalDigisAnalyzer::mehEcalMultvAEE[2]

private

Definition at line 151 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillECal().

MonitorElement* GlobalDigisAnalyzer::mehEcaln[2]

private

Definition at line 146 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillECal().

MonitorElement* GlobalDigisAnalyzer::mehEcalSHE[2]

private

Definition at line 149 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillECal().

MonitorElement* GlobalDigisAnalyzer::mehEcalSHEvAEESHE[2]

private

Definition at line 152 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillECal().

MonitorElement* GlobalDigisAnalyzer::mehEScalADC[3]

private

Definition at line 153 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillECal().

MonitorElement* GlobalDigisAnalyzer::mehEScaln

private

Definition at line 147 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillECal().

MonitorElement* GlobalDigisAnalyzer::mehHcalAEE[4]

private

Definition at line 170 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillHCal().

MonitorElement* GlobalDigisAnalyzer::mehHcalAEESHE[4]

private

Definition at line 172 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillHCal().

MonitorElement* GlobalDigisAnalyzer::mehHcaln[4]

private

Definition at line 169 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillHCal().

MonitorElement* GlobalDigisAnalyzer::mehHcalSHE[4]

private

Definition at line 171 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillHCal().

MonitorElement* GlobalDigisAnalyzer::mehHcalSHEvAEE[4]

private

Definition at line 173 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillHCal().

MonitorElement* GlobalDigisAnalyzer::mehRPCMuonn

private

Definition at line 220 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillMuon().

MonitorElement* GlobalDigisAnalyzer::mehRPCRes[5]

private

Definition at line 221 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillMuon().

MonitorElement* GlobalDigisAnalyzer::mehSiPixelADC[7]

private

Definition at line 190 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

MonitorElement* GlobalDigisAnalyzer::mehSiPixelCol[7]

private

Definition at line 192 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

MonitorElement* GlobalDigisAnalyzer::mehSiPixeln[7]

private

Definition at line 189 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

MonitorElement* GlobalDigisAnalyzer::mehSiPixelRow[7]

private

Definition at line 191 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

MonitorElement* GlobalDigisAnalyzer::mehSiStripADC[19]

private

Definition at line 182 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

MonitorElement* GlobalDigisAnalyzer::mehSiStripn[19]

private

Definition at line 181 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

MonitorElement* GlobalDigisAnalyzer::mehSiStripStrip[19]

private

Definition at line 183 of file GlobalDigisAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

edm::InputTag GlobalDigisAnalyzer::MuCSCStripSrc_

private

Definition at line 213 of file GlobalDigisAnalyzer.h.

Referenced by GlobalDigisAnalyzer().

edm::EDGetTokenT<CSCStripDigiCollection> GlobalDigisAnalyzer::MuCSCStripSrc_Token_

private

Definition at line 236 of file GlobalDigisAnalyzer.h.

Referenced by fillMuon(), and GlobalDigisAnalyzer().

edm::InputTag GlobalDigisAnalyzer::MuCSCWireSrc_

private

Definition at line 217 of file GlobalDigisAnalyzer.h.

Referenced by GlobalDigisAnalyzer().

edm::EDGetTokenT<CSCWireDigiCollection> GlobalDigisAnalyzer::MuCSCWireSrc_Token_

private

Definition at line 237 of file GlobalDigisAnalyzer.h.

Referenced by fillMuon(), and GlobalDigisAnalyzer().

edm::InputTag GlobalDigisAnalyzer::MuDTSrc_

private

Definition at line 204 of file GlobalDigisAnalyzer.h.

Referenced by GlobalDigisAnalyzer().

edm::EDGetTokenT<DTDigiCollection> GlobalDigisAnalyzer::MuDTSrc_Token_

private

Definition at line 235 of file GlobalDigisAnalyzer.h.

Referenced by fillMuon(), and GlobalDigisAnalyzer().

edm::InputTag GlobalDigisAnalyzer::MuRPCSrc_

private

Definition at line 223 of file GlobalDigisAnalyzer.h.

Referenced by GlobalDigisAnalyzer().

edm::EDGetTokenT<RPCDigiCollection> GlobalDigisAnalyzer::MuRPCSrc_Token_

private

Definition at line 238 of file GlobalDigisAnalyzer.h.

Referenced by fillMuon(), and GlobalDigisAnalyzer().

bool GlobalDigisAnalyzer::printProvenanceInfo

private

Definition at line 140 of file GlobalDigisAnalyzer.h.

Referenced by analyze(), and GlobalDigisAnalyzer().

edm::ESGetToken<RPCGeometry, MuonGeometryRecord> GlobalDigisAnalyzer::rpcGeomToken_

private

Definition at line 159 of file GlobalDigisAnalyzer.h.

Referenced by fillMuon(), and GlobalDigisAnalyzer().

edm::EDGetTokenT<edm::PSimHitContainer> GlobalDigisAnalyzer::RPCSimHit_Token_

private

Definition at line 242 of file GlobalDigisAnalyzer.h.

Referenced by fillMuon(), and GlobalDigisAnalyzer().

edm::InputTag GlobalDigisAnalyzer::SiPxlSrc_

private

Definition at line 194 of file GlobalDigisAnalyzer.h.

Referenced by GlobalDigisAnalyzer().

edm::EDGetTokenT<edm::DetSetVector<PixelDigi> > GlobalDigisAnalyzer::SiPxlSrc_Token_

private

Definition at line 234 of file GlobalDigisAnalyzer.h.

Referenced by fillTrk(), and GlobalDigisAnalyzer().

edm::InputTag GlobalDigisAnalyzer::SiStripSrc_

private

Definition at line 185 of file GlobalDigisAnalyzer.h.

Referenced by GlobalDigisAnalyzer().

edm::EDGetTokenT<edm::DetSetVector<SiStripDigi> > GlobalDigisAnalyzer::SiStripSrc_Token_

private

Definition at line 233 of file GlobalDigisAnalyzer.h.

Referenced by fillTrk(), and GlobalDigisAnalyzer().

int GlobalDigisAnalyzer::theCSCStripPedestalCount

private

Definition at line 215 of file GlobalDigisAnalyzer.h.

Referenced by fillMuon().

float GlobalDigisAnalyzer::theCSCStripPedestalSum

private

Definition at line 214 of file GlobalDigisAnalyzer.h.

Referenced by fillMuon().

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> GlobalDigisAnalyzer::tTopoToken_

private

Definition at line 160 of file GlobalDigisAnalyzer.h.

Referenced by fillTrk(), and GlobalDigisAnalyzer().

int GlobalDigisAnalyzer::verbosity

private

Definition at line 136 of file GlobalDigisAnalyzer.h.

Referenced by analyze(), fillECal(), fillHCal(), fillMuon(), fillTrk(), and GlobalDigisAnalyzer().