#include <myTKAnalyses/DigiInvestigator/src/OccupancyPlots.cc>

Inheritance diagram for OccupancyPlots:

Public Member Functions
	OccupancyPlots (const edm::ParameterSet &)

	~OccupancyPlots ()

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &) override

virtual void	beginJob () override

virtual void	beginRun (const edm::Run &, const edm::EventSetup &) override

virtual void	endJob () override

virtual void	endRun (const edm::Run &, const edm::EventSetup &) override

Private Attributes
TProfile **	m_avemultiplicity

TProfile **	m_aveoccupancy

TProfile **	m_averadius

TProfile **	m_avex

TProfile **	m_avey

TProfile **	m_avez

edm::FileInPath	m_fp

std::vector< edm::EDGetTokenT < std::map< unsigned int, int > > >	m_multiplicityMapTokens

TH1F **	m_nchannels_ideal

TH1F **	m_nchannels_real

std::vector< edm::EDGetTokenT < std::map< unsigned int, int > > >	m_occupancyMapTokens

RunHistogramManager	m_rhm

std::map< unsigned int, DetIdSelector >	m_wantedsubdets

TProfile **	m_xavedr

TProfile **	m_xavedrphi

TProfile **	m_xavedz

TProfile **	m_yavedr

TProfile **	m_yavedrphi

TProfile **	m_yavedz

TProfile **	m_zavedr

TProfile **	m_zavedrphi

TProfile **	m_zavedz

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

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 ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

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

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

Definition at line 67 of file OccupancyPlots.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 123 of file OccupancyPlots.cc.

References GlobalPosition_Frontier_DevDB_cff::tag.

123 :

124 m_multiplicityMapTokens(edm::vector_transform(iConfig.getParameter<std::vector<edm::InputTag> >("multiplicityMaps"), [this](edm::InputTag const & tag){return consumes<std::map<unsigned int, int> >(tag);})),

edm::ParameterSet::getParameter

T getParameter(std::string const &) const

edm::vector_transform

auto vector_transform(std::vector< InputType > const &input, Function predicate) -> std::vector< typename std::remove_cv< typename std::remove_reference< decltype(predicate(input.front()))>::type >::type >

Definition: transform.h:11

OccupancyPlots::m_multiplicityMapTokens

std::vector< edm::EDGetTokenT< std::map< unsigned int, int > > > m_multiplicityMapTokens

Definition: OccupancyPlots.cc:82

edm::EDConsumerBase::consumes

EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)

Definition: EDConsumerBase.h:86

GlobalPosition_Frontier_DevDB_cff.tag

tuple tag

Definition: GlobalPosition_Frontier_DevDB_cff.py:11

edm::InputTag

Definition: InputTag.h:17

OccupancyPlots::~OccupancyPlots ( )

Definition at line 166 of file OccupancyPlots.cc.

 {
 
    // do anything here that needs to be done at desctruction time
    // (e.g. close files, deallocate resources etc.)
 
 }

Member Function Documentation

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

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 181 of file OccupancyPlots.cc.

References edm::Event::getByToken(), m_avemultiplicity, m_aveoccupancy, m_multiplicityMapTokens, m_occupancyMapTokens, and VarParsing::mult.

 {
   using namespace edm;
 
 
   for(std::vector<edm::EDGetTokenT<std::map<unsigned int, int> > >::const_iterator mapToken = m_multiplicityMapTokens.begin();mapToken!=m_multiplicityMapTokens.end();++mapToken) {
 
     Handle<std::map<unsigned int, int> > mults;
     iEvent.getByToken(*mapToken,mults);
 
     for(std::map<unsigned int,int>::const_iterator mult=mults->begin();mult!=mults->end();mult++) {
       if(m_avemultiplicity && *m_avemultiplicity) (*m_avemultiplicity)->Fill(mult->first,mult->second);
     }
   }
 
 
 
   for(std::vector<edm::EDGetTokenT<std::map<unsigned int, int> > >::const_iterator mapToken = m_occupancyMapTokens.begin();mapToken!=m_occupancyMapTokens.end();++mapToken) {
 
     Handle<std::map<unsigned int, int> > occus;
     iEvent.getByToken(*mapToken,occus);
 
     for(std::map<unsigned int,int>::const_iterator occu=occus->begin();occu!=occus->end();occu++) {
       if(m_aveoccupancy && *m_aveoccupancy) (*m_aveoccupancy)->Fill(occu->first,occu->second);
     }
   }
 
 }

void OccupancyPlots::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 213 of file OccupancyPlots.cc.

 {
 
 }

void OccupancyPlots::beginRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 219 of file OccupancyPlots.cc.

References RunHistogramManager::beginRun(), and m_rhm.

                                                                         {
 
   m_rhm.beginRun(iRun);
 
 }

void OccupancyPlots::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 352 of file OccupancyPlots.cc.

352 {

353 }

void OccupancyPlots::endRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 226 of file OccupancyPlots.cc.

References cond::rpcobgas::detid, edm::FileInPath::fullPath(), edm::EventSetup::get(), SiPixelDetInfoFileReader::getAllDetIds(), SiStripDetInfoFileReader::getAllDetIds(), SiStripDetInfoFileReader::getNumberOfApvsAndStripLength(), m_averadius, m_avex, m_avey, m_avez, m_fp, m_nchannels_ideal, m_nchannels_real, m_wantedsubdets, m_xavedr, m_xavedrphi, m_xavedz, m_yavedr, m_yavedrphi, m_yavedz, m_zavedr, m_zavedrphi, m_zavedz, cppFunctionSkipper::operator, PV3DBase< T, PVType, FrameType >::perp(), position, matplotRender::reader, EgammaValidation_Wenu_cff::sel, DetId::Tracker, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                       {
 
 
   //  edm::ESHandle<GlobalTrackingGeometry> trkgeo;
   //  iSetup.get<GlobalTrackingGeometryRecord>().get("",trkgeo);
   edm::ESHandle<TrackerGeometry> trkgeo;
   iSetup.get<TrackerDigiGeometryRecord>().get("",trkgeo);
 
   const Local2DPoint center(0.,0.);
   const Local3DPoint locz(0.,0.,1.);
   const Local3DPoint locx(1.,0.,0.);
   const Local3DPoint locy(0.,1.,0.);
   const GlobalPoint origin(0.,0.,0.);
 
   TrackingGeometry::DetIdContainer detunits = trkgeo->detUnitIds();
 
   for(TrackingGeometry::DetIdContainer::const_iterator det = detunits.begin(); det!=detunits.end(); ++det) {
 
     if(det->det()!=DetId::Tracker) continue;
 
     edm::LogInfo("DetIdFromGeometry") << det->rawId();
 
     GlobalPoint position = trkgeo->idToDet(*det)->toGlobal(center);
     GlobalPoint zpos = trkgeo->idToDet(*det)->toGlobal(locz);
     GlobalPoint xpos = trkgeo->idToDet(*det)->toGlobal(locx);
     GlobalPoint ypos = trkgeo->idToDet(*det)->toGlobal(locy);
     GlobalVector posvect = position - origin;
     GlobalVector dz = zpos - position;
     GlobalVector dx = xpos - position;
     GlobalVector dy = ypos - position;
 
     double dzdr = posvect.perp()>0 ? (dz.x()*posvect.x()+dz.y()*posvect.y())/posvect.perp() : 0. ;
     double dxdr = posvect.perp()>0 ? (dx.x()*posvect.x()+dx.y()*posvect.y())/posvect.perp() : 0. ;
     double dydr = posvect.perp()>0 ? (dy.x()*posvect.x()+dy.y()*posvect.y())/posvect.perp() : 0. ;
 
     double dzdrphi = posvect.perp()>0 ? (dz.y()*posvect.x()-dz.x()*posvect.y())/posvect.perp() : 0. ;
     double dxdrphi = posvect.perp()>0 ? (dx.y()*posvect.x()-dx.x()*posvect.y())/posvect.perp() : 0. ;
     double dydrphi = posvect.perp()>0 ? (dy.y()*posvect.x()-dy.x()*posvect.y())/posvect.perp() : 0. ;
 
      for(std::map<unsigned int,DetIdSelector>::const_iterator sel=m_wantedsubdets.begin();sel!=m_wantedsubdets.end();++sel) {
 
        if(sel->second.isSelected(*det)) {
      edm::LogInfo("SelectedDetId") << sel->first;
      // average positions
      if(m_averadius && *m_averadius) (*m_averadius)->Fill(sel->first,position.perp());
      if(m_avez && *m_avez) (*m_avez)->Fill(sel->first,position.z());
      if(m_avex && *m_avex) (*m_avex)->Fill(sel->first,position.x());
      if(m_avey && *m_avey) (*m_avey)->Fill(sel->first,position.y());
      if(m_zavedr && *m_zavedr) (*m_zavedr)->Fill(sel->first,dzdr);
      if(m_zavedz && *m_zavedz) (*m_zavedz)->Fill(sel->first,dz.z());
      if(m_zavedrphi && *m_zavedrphi) (*m_zavedrphi)->Fill(sel->first,dzdrphi);
      if(m_xavedr && *m_xavedr) (*m_xavedr)->Fill(sel->first,dxdr);
      if(m_xavedz && *m_xavedz) (*m_xavedz)->Fill(sel->first,dx.z());
      if(m_xavedrphi && *m_xavedrphi) (*m_xavedrphi)->Fill(sel->first,dxdrphi);
      if(m_yavedr && *m_yavedr) (*m_yavedr)->Fill(sel->first,dydr);
      if(m_yavedz && *m_yavedz) (*m_yavedz)->Fill(sel->first,dy.z());
      if(m_yavedrphi && *m_yavedrphi) (*m_yavedrphi)->Fill(sel->first,dydrphi);
        }
      }
   }
 
 
   edm::ESHandle<SiStripQuality> quality;
   iSetup.get<SiStripQualityRcd>().get("",quality);
 
 
    SiStripDetInfoFileReader * reader=edm::Service<SiStripDetInfoFileReader>().operator->();
 
    const std::vector<uint32_t>& detids = reader->getAllDetIds();
 
    for(std::vector<uint32_t>::const_iterator detid=detids.begin();detid!=detids.end();++detid) {
 
      int nchannideal = reader->getNumberOfApvsAndStripLength(*detid).first*128;
      //     int nchannreal = reader->getNumberOfApvsAndStripLength(*detid).first*128;
      int nchannreal = 0;
      for(int strip = 0; strip < nchannideal; ++strip) {
        if(!quality->IsStripBad(*detid,strip)) ++nchannreal;
      }
 
 
      for(std::map<unsigned int,DetIdSelector>::const_iterator sel=m_wantedsubdets.begin();sel!=m_wantedsubdets.end();++sel) {
 
        if(sel->second.isSelected(*detid)) {
      if(m_nchannels_ideal && *m_nchannels_ideal) (*m_nchannels_ideal)->Fill(sel->first,nchannideal);
      if(m_nchannels_real && *m_nchannels_real) (*m_nchannels_real)->Fill(sel->first,nchannreal);
        }
      }
 
    }
 
 
   edm::ESHandle<SiPixelQuality> pxlquality;
   iSetup.get<SiPixelQualityRcd>().get("",pxlquality);
 
 
    SiPixelDetInfoFileReader pxlreader(m_fp.fullPath());
 
    const std::vector<uint32_t>& pxldetids = pxlreader.getAllDetIds();
 
    for(std::vector<uint32_t>::const_iterator detid=pxldetids.begin();detid!=pxldetids.end();++detid) {
 
      int nchannideal = pxlreader.getDetUnitDimensions(*detid).first*pxlreader.getDetUnitDimensions(*detid).second;
      int nchannreal = 0;
      if(!pxlquality->IsModuleBad(*detid)) {
        nchannreal = pxlreader.getDetUnitDimensions(*detid).first*pxlreader.getDetUnitDimensions(*detid).second;
      }
      /*
      int nchannreal = 0;
      for(int strip = 0; strip < nchannideal; ++strip) {
        if(!quality->IsStripBad(*detid,strip)) ++nchannreal;
      }
      */
 
      for(std::map<unsigned int,DetIdSelector>::const_iterator sel=m_wantedsubdets.begin();sel!=m_wantedsubdets.end();++sel) {
 
        if(sel->second.isSelected(*detid)) {
      if(m_nchannels_ideal && *m_nchannels_ideal) (*m_nchannels_ideal)->Fill(sel->first,nchannideal);
      if(m_nchannels_real && *m_nchannels_real) (*m_nchannels_real)->Fill(sel->first,nchannreal);
        }
      }
 
    }
 
 }

Member Data Documentation

TProfile** OccupancyPlots::m_avemultiplicity

private

Definition at line 89 of file OccupancyPlots.cc.

Referenced by analyze().

TProfile** OccupancyPlots::m_aveoccupancy

private

Definition at line 90 of file OccupancyPlots.cc.

Referenced by analyze().

TProfile** OccupancyPlots::m_averadius

private

Definition at line 95 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_avex

private

Definition at line 97 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_avey

private

Definition at line 98 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_avez

private

Definition at line 96 of file OccupancyPlots.cc.

Referenced by endRun().

edm::FileInPath OccupancyPlots::m_fp

private

Definition at line 84 of file OccupancyPlots.cc.

Referenced by endRun().

std::vector<edm::EDGetTokenT<std::map<unsigned int, int> > > OccupancyPlots::m_multiplicityMapTokens

private

Definition at line 82 of file OccupancyPlots.cc.

Referenced by analyze().

TH1F** OccupancyPlots::m_nchannels_ideal

private

Definition at line 92 of file OccupancyPlots.cc.

Referenced by endRun().

TH1F** OccupancyPlots::m_nchannels_real

private

Definition at line 93 of file OccupancyPlots.cc.

Referenced by endRun().

std::vector<edm::EDGetTokenT<std::map<unsigned int, int> > > OccupancyPlots::m_occupancyMapTokens

private

Definition at line 83 of file OccupancyPlots.cc.

Referenced by analyze().

RunHistogramManager OccupancyPlots::m_rhm

private

Definition at line 86 of file OccupancyPlots.cc.

Referenced by beginRun().

std::map<unsigned int,DetIdSelector> OccupancyPlots::m_wantedsubdets

private

Definition at line 87 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_xavedr

private

Definition at line 105 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_xavedrphi

private

Definition at line 107 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_xavedz

private

Definition at line 106 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_yavedr

private

Definition at line 102 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_yavedrphi

private

Definition at line 104 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_yavedz

private

Definition at line 103 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_zavedr

private

Definition at line 99 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_zavedrphi

private

Definition at line 101 of file OccupancyPlots.cc.

Referenced by endRun().

TProfile** OccupancyPlots::m_zavedz

private

Definition at line 100 of file OccupancyPlots.cc.

Referenced by endRun().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation