TrackerRemapper Class Reference

Inheritance diagram for TrackerRemapper:

Public Types
enum	AnalyzeData { RECHITS = 1, DIGIS, CLUSTERS }
enum	OpMode { MODE_ANALYZE = 0, MODE_REMAP = 1 }
enum	PixelLayerEnum {   INVALID = 0, PXB_L1, PXB_L2, PXB_L3,   PXB_L4, PXF_R1, PXF_R2 }
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
template<>
void	analyzeGeneric (const edm::Event &iEvent, const edm::EDGetTokenT< reco::TrackCollection > &t)
	TrackerRemapper (const edm::ParameterSet &)
	~TrackerRemapper () override
Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources >
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
bool	wantsInputProcessBlocks () const final
bool	wantsProcessBlocks () const final
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDAnalyzerBase () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector < ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
template<class T >
void	analyzeGeneric (const edm::Event &iEvent, const edm::EDGetTokenT< T > &src)
void	analyzeRechits (const edm::Event &iEvent, const edm::EDGetTokenT< reco::TrackCollection > &src)
void	beginJob () override
void	bookBins ()
void	endJob () override
void	fillBarrelRemap (TFile *rootFileHandle, const TrackerGeometry *theTrackerGeometry, const TrackerTopology *tt)
void	fillEndcapRemap (TFile *rootFileHandle, const TrackerGeometry *theTrackerGeometry, const TrackerTopology *tt)
void	fillPixelRemap (const TrackerGeometry *theTrackerGeometry, const TrackerTopology *tt)
void	fillStripRemap ()
void	preparePixelNames ()
void	prepareStripNames ()
void	readVertices (double &minx, double &maxx, double &miny, double &maxy)

Private Attributes
edm::EDGetTokenT < edmNew::DetSetVector < SiStripCluster > >	clusterSrcToken
edm::EDGetTokenT < edmNew::DetSetVector < SiStripDigi > >	digiSrcToken
edm::Service< TFileService >	fs
const edm::ESGetToken < TrackerGeometry, TrackerDigiGeometryRecord >	geomToken_
const edm::ParameterSet &	iConfig
int	m_analyzeMode
std::map< unsigned, std::string >	m_analyzeModeNameMap
std::map< long, TGraph * >	m_bins
std::vector< unsigned >	m_detIdVector
int	m_opMode
std::string	m_pixelBaseDir
std::string	m_pixelDesiredHistogramBarrel
std::string	m_pixelDesiredHistogramDisk
std::map< unsigned, std::string >	m_pixelHistnameMap
std::string	m_pixelRemapFile
std::string	m_stripBaseDir
std::string	m_stripDesiredHistogram
std::map< unsigned, std::string >	m_stripHistnameMap
std::string	m_stripRemapFile
const TkDetMap *	m_tkdetmap
edm::EDGetTokenT < reco::TrackCollection >	rechitSrcToken
std::string	runString
const edm::ESGetToken < TkDetMap, TrackerTopologyRcd >	tkDetMapToken_
const edm::ESGetToken < TrackerTopology, TrackerTopologyRcd >	topoToken_
TH2Poly *	trackerMap {nullptr}

Additional Inherited Members
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
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)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 53 of file TrackerRemapper.cc.

Member Enumeration Documentation

enum TrackerRemapper::AnalyzeData

Enumerator
RECHITS
DIGIS
CLUSTERS

Definition at line 72 of file TrackerRemapper.cc.

                   {
    RECHITS = 1,
    DIGIS,
    CLUSTERS,
  };

enum TrackerRemapper::OpMode

Enumerator
MODE_ANALYZE
MODE_REMAP

Definition at line 78 of file TrackerRemapper.cc.

{ MODE_ANALYZE = 0, MODE_REMAP = 1 };

enum TrackerRemapper::PixelLayerEnum

Enumerator
INVALID
PXB_L1
PXB_L2
PXB_L3
PXB_L4
PXF_R1
PXF_R2

Definition at line 60 of file TrackerRemapper.cc.

                      {
    INVALID = 0,

    PXB_L1,
    PXB_L2,
    PXB_L3,
    PXB_L4,

    PXF_R1,
    PXF_R2
  };

Constructor & Destructor Documentation

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

explicit

Definition at line 164 of file TrackerRemapper.cc.

References CLUSTERS, clusterSrcToken, DIGIS, digiSrcToken, Exception, edm::ParameterSet::getParameter(), m_analyzeMode, m_analyzeModeNameMap, m_opMode, m_pixelBaseDir, m_pixelDesiredHistogramBarrel, m_pixelDesiredHistogramDisk, m_pixelRemapFile, m_stripBaseDir, m_stripDesiredHistogram, m_stripRemapFile, MODE_ANALYZE, MODE_REMAP, preparePixelNames(), prepareStripNames(), RECHITS, rechitSrcToken, runString, and AlCaHLTBitMon_QueryRunRegistry::string.

    : geomToken_(esConsumes()),
      topoToken_(esConsumes()),
      tkDetMapToken_(esConsumes()),
      iConfig(iConfig),
      m_opMode(iConfig.getParameter<int>("opMode")),
      m_analyzeMode(iConfig.getParameter<int>("analyzeMode")) {
  usesResource("TFileService");

  if (m_opMode == MODE_REMAP) {
    m_stripRemapFile = iConfig.getParameter<std::string>("stripRemapFile");
    m_stripDesiredHistogram = iConfig.getParameter<std::string>("stripHistogram");
    runString = iConfig.getParameter<std::string>("runString");

    m_pixelRemapFile = std::string("DQM_V0001_PixelPhase1_R000305516.root");

    m_stripBaseDir = std::string("DQMData/Run " + runString + "/SiStrip/Run summary/MechanicalView/");
    m_pixelBaseDir = std::string("DQMData/Run " + runString + "/PixelPhase1/Run summary/Phase1_MechanicalView/");

    m_pixelDesiredHistogramBarrel = std::string("adc_per_SignedModule_per_SignedLadder");
    m_pixelDesiredHistogramDisk = std::string("adc_per_PXDisk_per_SignedBladePanel");

    prepareStripNames();
    preparePixelNames();
  } else if (m_opMode == MODE_ANALYZE) {
    m_analyzeModeNameMap[RECHITS] = "# Rechits";
    m_analyzeModeNameMap[DIGIS] = "# Digis";
    m_analyzeModeNameMap[CLUSTERS] = "# Clusters";

    switch (m_analyzeMode) {
      case RECHITS:
        rechitSrcToken = consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("src"));
        break;
      case DIGIS:
        digiSrcToken = consumes<edmNew::DetSetVector<SiStripDigi>>(iConfig.getParameter<edm::InputTag>("src"));
        break;
      case CLUSTERS:
        clusterSrcToken = consumes<edmNew::DetSetVector<SiStripCluster>>(iConfig.getParameter<edm::InputTag>("src"));
        break;
      default:
        edm::LogError("LogicError") << "Unrecognized analyze mode!" << std::endl;
    }
  } else {
    throw cms::Exception("TrackerRemapper") << "Unrecognized operations mode!" << std::endl;
  }

  // TColor::SetPalette(1);
}

TrackerRemapper::~TrackerRemapper ( )

override

Definition at line 310 of file TrackerRemapper.cc.

{}

Member Function Documentation

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 313 of file TrackerRemapper.cc.

References analyzeGeneric(), bookBins(), clusterSrcToken, digiSrcToken, fillPixelRemap(), fillStripRemap(), geomToken_, edm::EventSetup::getData(), m_analyzeMode, m_opMode, m_tkdetmap, MODE_ANALYZE, MODE_REMAP, rechitSrcToken, tkDetMapToken_, topoToken_, trackerMap, and groupFilesInBlocks::tt.

{
  // get the ES products

  const TrackerGeometry* theTrackerGeometry = &iSetup.getData(geomToken_);
  const TrackerTopology* tt = &iSetup.getData(topoToken_);
  m_tkdetmap = &iSetup.getData(tkDetMapToken_);

  if (!trackerMap)
    bookBins();

  if (m_opMode == MODE_ANALYZE) {
    switch (m_analyzeMode) {
      case AnalyzeData::RECHITS:
        analyzeGeneric(iEvent, rechitSrcToken);
        break;
      case AnalyzeData::DIGIS:
        analyzeGeneric(iEvent, digiSrcToken);
        break;
      case AnalyzeData::CLUSTERS:
        analyzeGeneric(iEvent, clusterSrcToken);
        break;
      default:
        edm::LogError("LogicError") << "Unrecognized Analyze mode!" << std::endl;
        return;
    }
  } else if (m_opMode == MODE_REMAP) {
    fillStripRemap();
    fillPixelRemap(theTrackerGeometry, tt);
  }
}

template<class T >

void TrackerRemapper::analyzeGeneric ( const edm::Event &  iEvent, const edm::EDGetTokenT< T > &  src  )

private

Definition at line 137 of file TrackerRemapper.cc.

References edm::Event::getByToken(), input, and edm::HandleBase::isValid().

Referenced by analyze().

{
  edm::Handle<T> input;
  iEvent.getByToken(src, input);

  if (!input.isValid()) {
    edm::LogError("TrackerRemapper") << "<GENERIC> not found... Aborting...\n";
    return;
  }

  typename T::const_iterator it;
  for (it = input->begin(); it != input->end(); ++it) {
    auto id = DetId(it->detId());
    trackerMap->Fill(TString::Format("%ld", (long)id.rawId()), it->size());
  }
}

template<>

void TrackerRemapper::analyzeGeneric	(	const edm::Event &	iEvent,
		const edm::EDGetTokenT< reco::TrackCollection > &	t
	)

Definition at line 157 of file TrackerRemapper.cc.

{
  analyzeRechits(iEvent, t);
}

void TrackerRemapper::analyzeRechits ( const edm::Event &  iEvent, const edm::EDGetTokenT< reco::TrackCollection > &  src  )

private

Definition at line 347 of file TrackerRemapper.cc.

References edm::Event::getByToken(), mps_fire::i, edm::HandleBase::isValid(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, DetId::subdetId(), HLT_FULL_cff::track, trackerMap, and tracks.

{
  edm::Handle<reco::TrackCollection> tracks;
  iEvent.getByToken(src, tracks);
  if (!tracks.isValid()) {
    LogInfo("Analyzer") << "reco::TrackCollection not found... Aborting...\n";
    return;
  }

  for (auto const& track : *tracks) {
    auto recHitsBegin = track.recHitsBegin();
    for (unsigned i = 0; i < track.recHitsSize(); ++i) {
      auto recHit = *(recHitsBegin + i);
      if (!recHit->isValid())
        continue;

      DetId id = recHit->geographicalId();
      unsigned subdetId = id.subdetId();

      //reject Pixel
      if (subdetId == PixelSubdetector::PixelBarrel || subdetId == PixelSubdetector::PixelEndcap)
        continue;

      // NRECHITS
      trackerMap->Fill(TString::Format("%ld", (long)id.rawId()), 1);
    }
  }
}

void TrackerRemapper::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 623 of file TrackerRemapper.cc.

{}

void TrackerRemapper::bookBins ( )

private

Definition at line 378 of file TrackerRemapper.cc.

References fs, TFileService::getBareDirectory(), m_analyzeMode, m_analyzeModeNameMap, m_bins, m_opMode, m_stripDesiredHistogram, TFileService::make(), margin, MODE_ANALYZE, MODE_REMAP, readVertices(), runString, AlCaHLTBitMon_QueryRunRegistry::string, and trackerMap.

Referenced by analyze().

{
  // Read vertices from file
  double minx = 0xFFFFFF, maxx = -0xFFFFFF, miny = 0xFFFFFF, maxy = -0xFFFFFFF;
  readVertices(minx, maxx, miny, maxy);

  TObject* ghostObj = fs->make<TH2Poly>("ghost", "ghost", -1, 1, -1, 1);

  TDirectory* topDir = fs->getBareDirectory();
  topDir->cd();

  int margin = 50;
  std::string mapTitle;
  switch (m_opMode) {
    case MODE_ANALYZE:
      mapTitle = std::string(m_analyzeModeNameMap[m_analyzeMode] + " - " + runString);
      break;
    case MODE_REMAP:
      mapTitle = std::string(m_stripDesiredHistogram + " - " + runString);
      break;
  }

  trackerMap = new TH2Poly("Tracker Map", mapTitle.c_str(), minx - margin, maxx + margin, miny - margin, maxy + margin);
  trackerMap->SetFloat();
  trackerMap->SetOption("COLZ");
  trackerMap->SetStats(false);

  for (auto pair : m_bins) {
    trackerMap->AddBin(pair.second->Clone());
  }

  topDir->Add(trackerMap);

  ghostObj->Delete();  //not needed any more
}

void TrackerRemapper::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 625 of file TrackerRemapper.cc.

{}

void TrackerRemapper::fillBarrelRemap ( TFile *  rootFileHandle, const TrackerGeometry *  theTrackerGeometry, const TrackerTopology *  tt  )

private

Definition at line 541 of file TrackerRemapper.cc.

References TrackerGeometry::detsPXB(), GeomDet::geographicalId(), mps_fire::i, PixelBarrelName::ladderName(), phase1PixelTopology::layer, m_pixelHistnameMap, callgraph::module, TrackerTopology::pxbLayer(), TrackerTopology::pxbModule(), DetId::rawId(), PixelBarrelName::shell(), trackerMap, and hgcalPerformanceValidation::val.

Referenced by fillPixelRemap().

{
  TrackingGeometry::DetContainer pxb = theTrackerGeometry->detsPXB();

  for (auto& i : pxb) {
    const GeomDet* det = i;

    PXBDetId id = det->geographicalId();
    long rawid = id.rawId();

    int module = tt->pxbModule(id);
    int layer = tt->pxbLayer(id);

    int signedOnlineModule = module - 4;
    if (signedOnlineModule <= 0)
      --signedOnlineModule;

    PixelBarrelName pixelBarrelName = PixelBarrelName(id, tt, true);
    int onlineShell = pixelBarrelName.shell();

    int signedOnlineLadder = ((onlineShell & 1) ? -pixelBarrelName.ladderName() : pixelBarrelName.ladderName());

    const TProfile2D* histHandle = (TProfile2D*)rootFileHandle->Get(m_pixelHistnameMap[layer].c_str());

    unsigned nx = histHandle->GetNbinsX();
    unsigned ny = histHandle->GetNbinsY();

    unsigned binX = signedOnlineModule + ((nx + 1) >> 1);
    unsigned binY = (signedOnlineLadder) + ((ny + 1) >> 1);

    double val = histHandle->GetBinContent(binX, binY);

    trackerMap->Fill(TString::Format("%ld", rawid), val);
  }
}

void TrackerRemapper::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 628 of file TrackerRemapper.cc.

References edm::ParameterSetDescription::add(), edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, edm::ParameterSetDescription::ifValue(), HLT_FULL_cff::InputTag, or, edm::ParameterSetDescription::setComment(), and edm::ParameterDescriptionNode::setComment().

{
  edm::ParameterSetDescription desc;
  desc.setComment(
      "Creates TH2Poly Strip Tracker maps by either analyzing the event or remapping exising DQM historams");
  desc.add<edm::InputTag>("src", edm::InputTag("generalTracks"));
  desc.ifValue(edm::ParameterDescription<int>("opMode", 0, true),
               0 >> edm::EmptyGroupDescription() or
                   1 >> (edm::ParameterDescription<std::string>("stripRemapFile", "", true) and
                         edm::ParameterDescription<std::string>("stripHistogram", "", true) and
                         edm::ParameterDescription<std::string>("runString", "", true)))
      ->setComment("0 for Analyze, 1 for Remap");

  desc.ifValue(edm::ParameterDescription<int>("analyzeMode", 1, true), edm::allowedValues<int>(1, 2, 3))
      ->setComment("1=Rechits, 2=Digis, 3=Clusters");

  //desc.add<unsigned int>("analyzeMode", 1)->setComment("1=Rechits, 2=Digis, 3=Clusters");
  //desc.add<unsigned int>("opMode", 0)->setComment("0 for Analyze, 1 for Remap");
  //desc.addOptional<std::string>("stripRemapFile","" )->setComment("file name to analyze, will come from the config file");
  //desc.addOptional<std::string>("stripHistogram","TkHMap_NumberValidHits" )->setComment("histogram to use to remap");
  //desc.addOptional<std::string>("runString", "")->setComment("run number, will come form config file");
  descriptions.addWithDefaultLabel(desc);
}

void TrackerRemapper::fillEndcapRemap ( TFile *  rootFileHandle, const TrackerGeometry *  theTrackerGeometry, const TrackerTopology *  tt  )

private

Definition at line 581 of file TrackerRemapper.cc.

References PixelEndcapName::bladeName(), TrackerGeometry::detsPXF(), GeomDet::geographicalId(), PixelEndcapName::halfCylinder(), mps_fire::i, phase1PixelTopology::layer, TrackerTopology::layer(), m_pixelHistnameMap, PixelEndcapName::pannelName(), PixelEndcapName::ringName(), TrackerTopology::side(), trackerMap, and hgcalPerformanceValidation::val.

Referenced by fillPixelRemap().

{
  TrackingGeometry::DetContainer pxf = theTrackerGeometry->detsPXF();

  for (auto& i : pxf) {
    const GeomDet* det = i;

    PXFDetId id = det->geographicalId();

    int side = tt->side(id);
    int disk = tt->layer(id);

    long rawid = id.rawId();

    PixelEndcapName pixelEndcapName = PixelEndcapName(PXFDetId(rawid), tt, true);

    unsigned layer = pixelEndcapName.ringName() - 1 + PixelLayerEnum::PXF_R1;
    const TProfile2D* histHandle = (TProfile2D*)rootFileHandle->Get(m_pixelHistnameMap[layer].c_str());

    // ---- REMAP (Online -> Offline)
    unsigned nx = histHandle->GetNbinsX();
    unsigned ny = histHandle->GetNbinsY();

    int onlineBlade = pixelEndcapName.bladeName();
    bool isInnerOnlineBlade = !(pixelEndcapName.halfCylinder() & 1);  // inner -> blade > 0 (?)

    int signedOnlineBlade = (isInnerOnlineBlade) ? onlineBlade : -onlineBlade;
    int signedDisk = (side == 2) ? disk : -disk;
    int pannel = pixelEndcapName.pannelName() - 1;

    unsigned binX = signedDisk + ((nx + 1) >> 1);
    unsigned binY = (signedOnlineBlade * 2) + (ny >> 1);

    double val = histHandle->GetBinContent(binX, binY + pannel);

    trackerMap->Fill(TString::Format("%ld", rawid), val);
  }
}

void TrackerRemapper::fillPixelRemap ( const TrackerGeometry *  theTrackerGeometry, const TrackerTopology *  tt  )

private

Definition at line 525 of file TrackerRemapper.cc.

References fillBarrelRemap(), fillEndcapRemap(), and m_pixelRemapFile.

Referenced by analyze().

{
  TFile* rootFileHandle = new TFile(m_pixelRemapFile.c_str());

  if (!rootFileHandle) {
    edm::LogError("TrackerRemapper") << "Could not find file:\n\t" << m_pixelRemapFile << std::endl;
    return;
  }
  fillBarrelRemap(rootFileHandle, theTrackerGeometry, tt);
  fillEndcapRemap(rootFileHandle, theTrackerGeometry, tt);

  rootFileHandle->Close();
}

void TrackerRemapper::fillStripRemap ( )

private

Definition at line 485 of file TrackerRemapper.cc.

References TkDetMap::getComponents(), TkDetMap::getDetFromBin(), phase1PixelTopology::layer, m_stripHistnameMap, m_stripRemapFile, m_tkdetmap, trackerMap, and hgcalPerformanceValidation::val.

Referenced by analyze().

{
  int nchX;
  int nchY;
  double lowX, highX;
  double lowY, highY;

  TFile* rootFileHandle = new TFile(m_stripRemapFile.c_str());

  for (int layer = TkLayerMap::TkLayerEnum::TIB_L1; layer <= TkLayerMap::TkLayerEnum::TECP_W9; ++layer) {
    m_tkdetmap->getComponents(layer, nchX, lowX, highX, nchY, lowY, highY);

    const TProfile2D* histHandle = (TProfile2D*)rootFileHandle->Get(m_stripHistnameMap[layer].c_str());

    if (!histHandle) {
      edm::LogError("TrackerRemapper") << "Could not find histogram:\n\t" << m_stripHistnameMap[layer] << std::endl;
      return;
    }

    for (unsigned binx = 1; binx <= (unsigned)nchX; ++binx) {
      for (unsigned biny = 1; biny <= (unsigned)nchY; ++biny) {
        long rawid = m_tkdetmap->getDetFromBin(layer, binx, biny);

        if (rawid)  //bin represents real module -> go to file
        {
          double val = histHandle->GetBinContent(binx, biny);

          // edm::LogInfo("TrackerRemapper") << rawid << " " << val << "\n";

          trackerMap->Fill(TString::Format("%ld", rawid), val);
        }
      }
    }
  }

  rootFileHandle->Close();
}

void TrackerRemapper::preparePixelNames ( )

private

Definition at line 292 of file TrackerRemapper.cc.

References m_pixelBaseDir, m_pixelDesiredHistogramBarrel, m_pixelDesiredHistogramDisk, and m_pixelHistnameMap.

Referenced by TrackerRemapper().

{
  m_pixelHistnameMap[PixelLayerEnum::PXB_L1] =
      m_pixelBaseDir + "PXBarrel/" + m_pixelDesiredHistogramBarrel + "_PXLayer_1;1";
  m_pixelHistnameMap[PixelLayerEnum::PXB_L2] =
      m_pixelBaseDir + "PXBarrel/" + m_pixelDesiredHistogramBarrel + "_PXLayer_2;1";
  m_pixelHistnameMap[PixelLayerEnum::PXB_L3] =
      m_pixelBaseDir + "PXBarrel/" + m_pixelDesiredHistogramBarrel + "_PXLayer_3;1";
  m_pixelHistnameMap[PixelLayerEnum::PXB_L4] =
      m_pixelBaseDir + "PXBarrel/" + m_pixelDesiredHistogramBarrel + "_PXLayer_4;1";

  m_pixelHistnameMap[PixelLayerEnum::PXF_R1] =
      m_pixelBaseDir + "PXForward/" + m_pixelDesiredHistogramDisk + "_PXRing_1;1";
  m_pixelHistnameMap[PixelLayerEnum::PXF_R2] =
      m_pixelBaseDir + "PXForward/" + m_pixelDesiredHistogramDisk + "_PXRing_2;1";
}

void TrackerRemapper::prepareStripNames ( )

private

Definition at line 214 of file TrackerRemapper.cc.

References m_stripBaseDir, m_stripDesiredHistogram, and m_stripHistnameMap.

Referenced by TrackerRemapper().

{
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TIB_L1] =
      m_stripBaseDir + "TIB/layer_1/" + m_stripDesiredHistogram + "_TIB_L1;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TIB_L2] =
      m_stripBaseDir + "TIB/layer_2/" + m_stripDesiredHistogram + "_TIB_L2;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TIB_L3] =
      m_stripBaseDir + "TIB/layer_3/" + m_stripDesiredHistogram + "_TIB_L3;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TIB_L4] =
      m_stripBaseDir + "TIB/layer_4/" + m_stripDesiredHistogram + "_TIB_L4;1";

  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TIDM_D1] =
      m_stripBaseDir + "TID/MINUS/wheel_1/" + m_stripDesiredHistogram + "_TIDM_D1;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TIDM_D2] =
      m_stripBaseDir + "TID/MINUS/wheel_2/" + m_stripDesiredHistogram + "_TIDM_D2;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TIDM_D3] =
      m_stripBaseDir + "TID/MINUS/wheel_3/" + m_stripDesiredHistogram + "_TIDM_D3;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TIDP_D1] =
      m_stripBaseDir + "TID/PLUS/wheel_1/" + m_stripDesiredHistogram + "_TIDP_D1;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TIDP_D2] =
      m_stripBaseDir + "TID/PLUS/wheel_2/" + m_stripDesiredHistogram + "_TIDP_D2;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TIDP_D3] =
      m_stripBaseDir + "TID/PLUS/wheel_3/" + m_stripDesiredHistogram + "_TIDP_D3;1";

  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TOB_L1] =
      m_stripBaseDir + "TOB/layer_1/" + m_stripDesiredHistogram + "_TOB_L1;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TOB_L2] =
      m_stripBaseDir + "TOB/layer_2/" + m_stripDesiredHistogram + "_TOB_L2;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TOB_L3] =
      m_stripBaseDir + "TOB/layer_3/" + m_stripDesiredHistogram + "_TOB_L3;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TOB_L4] =
      m_stripBaseDir + "TOB/layer_4/" + m_stripDesiredHistogram + "_TOB_L4;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TOB_L5] =
      m_stripBaseDir + "TOB/layer_5/" + m_stripDesiredHistogram + "_TOB_L5;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TOB_L6] =
      m_stripBaseDir + "TOB/layer_6/" + m_stripDesiredHistogram + "_TOB_L6;1";

  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TECM_W1] =
      m_stripBaseDir + "TEC/MINUS/wheel_1/" + m_stripDesiredHistogram + "_TECM_W1;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TECM_W2] =
      m_stripBaseDir + "TEC/MINUS/wheel_2/" + m_stripDesiredHistogram + "_TECM_W2;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TECM_W3] =
      m_stripBaseDir + "TEC/MINUS/wheel_3/" + m_stripDesiredHistogram + "_TECM_W3;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TECM_W4] =
      m_stripBaseDir + "TEC/MINUS/wheel_4/" + m_stripDesiredHistogram + "_TECM_W4;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TECM_W5] =
      m_stripBaseDir + "TEC/MINUS/wheel_5/" + m_stripDesiredHistogram + "_TECM_W5;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TECM_W6] =
      m_stripBaseDir + "TEC/MINUS/wheel_6/" + m_stripDesiredHistogram + "_TECM_W6;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TECM_W7] =
      m_stripBaseDir + "TEC/MINUS/wheel_7/" + m_stripDesiredHistogram + "_TECM_W7;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TECM_W8] =
      m_stripBaseDir + "TEC/MINUS/wheel_8/" + m_stripDesiredHistogram + "_TECM_W8;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TECM_W9] =
      m_stripBaseDir + "TEC/MINUS/wheel_9/" + m_stripDesiredHistogram + "_TECM_W9;1";

  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TECP_W1] =
      m_stripBaseDir + "TEC/PLUS/wheel_1/" + m_stripDesiredHistogram + "_TECP_W1;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TECP_W2] =
      m_stripBaseDir + "TEC/PLUS/wheel_2/" + m_stripDesiredHistogram + "_TECP_W2;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TECP_W3] =
      m_stripBaseDir + "TEC/PLUS/wheel_3/" + m_stripDesiredHistogram + "_TECP_W3;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TECP_W4] =
      m_stripBaseDir + "TEC/PLUS/wheel_4/" + m_stripDesiredHistogram + "_TECP_W4;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TECP_W5] =
      m_stripBaseDir + "TEC/PLUS/wheel_5/" + m_stripDesiredHistogram + "_TECP_W5;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TECP_W6] =
      m_stripBaseDir + "TEC/PLUS/wheel_6/" + m_stripDesiredHistogram + "_TECP_W6;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TECP_W7] =
      m_stripBaseDir + "TEC/PLUS/wheel_7/" + m_stripDesiredHistogram + "_TECP_W7;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TECP_W8] =
      m_stripBaseDir + "TEC/PLUS/wheel_8/" + m_stripDesiredHistogram + "_TECP_W8;1";
  m_stripHistnameMap[TkLayerMap::TkLayerEnum::TECP_W9] =
      m_stripBaseDir + "TEC/PLUS/wheel_9/" + m_stripDesiredHistogram + "_TECP_W9;1";
}

void TrackerRemapper::readVertices ( double &  minx, double &  maxx, double &  miny, double &  maxy  )

private

Definition at line 416 of file TrackerRemapper.cc.

References mps_check::array, submitPVResolutionJobs::count, Exception, contentValuesFiles::fullPath, mps_fire::i, recoMuon::in, fastTrackerRecHitType::isPixel(), geometryCSVtoXML::line, m_bins, m_detIdVector, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, AlCaHLTBitMon_QueryRunRegistry::string, DetId::subdetId(), x, and y.

Referenced by bookBins().

{
  std::ifstream in;

  // TPolyline vertices stored at https://github.com/cms-data/DQM-SiStripMonitorClient
  in.open(edm::FileInPath("DQM/SiStripMonitorClient/data/Geometry/tracker_map_bare").fullPath().c_str());

  unsigned count = 0;

  if (!in.good()) {
    throw cms::Exception("TrackerRemapper") << "Error Reading File" << std::endl;
  }
  while (in.good()) {
    long detid = 0;
    double x[5], y[5];

    std::string line;
    std::getline(in, line);
    ++count;

    TString string(line);
    TObjArray* array = string.Tokenize(" ");
    int ix{0}, iy{0};
    bool isPixel{false};
    for (int i = 0; i < array->GetEntries(); ++i) {
      if (i == 0) {
        detid = static_cast<TObjString*>(array->At(i))->String().Atoll();

        // Drop Pixel Data
        DetId detId(detid);
        if (detId.subdetId() == PixelSubdetector::PixelBarrel || detId.subdetId() == PixelSubdetector::PixelEndcap) {
          isPixel = true;
          break;
        }
      } else {
        if (i % 2 == 0) {
          x[ix] = static_cast<TObjString*>(array->At(i))->String().Atof();

          if (x[ix] < minx)
            minx = x[ix];
          if (x[ix] > maxx)
            maxx = x[ix];

          ++ix;
        } else {
          y[iy] = static_cast<TObjString*>(array->At(i))->String().Atof();

          if (y[iy] < miny)
            miny = y[iy];
          if (y[iy] > maxy)
            maxy = y[iy];

          ++iy;
        }
      }
    }

    if (isPixel)
      continue;

    m_detIdVector.push_back(detid);
    m_bins[detid] = new TGraph(ix, x, y);
    m_bins[detid]->SetName(TString::Format("%ld", detid));
    m_bins[detid]->SetTitle(TString::Format("Module ID=%ld", detid));
  }
}

Member Data Documentation

edm::EDGetTokenT<edmNew::DetSetVector<SiStripCluster> > TrackerRemapper::clusterSrcToken

private

Definition at line 132 of file TrackerRemapper.cc.

Referenced by analyze(), and TrackerRemapper().

edm::EDGetTokenT<edmNew::DetSetVector<SiStripDigi> > TrackerRemapper::digiSrcToken

private

Definition at line 131 of file TrackerRemapper.cc.

Referenced by analyze(), and TrackerRemapper().

edm::Service<TFileService> TrackerRemapper::fs

private

Definition at line 105 of file TrackerRemapper.cc.

Referenced by bookBins().

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> TrackerRemapper::geomToken_

private

Definition at line 101 of file TrackerRemapper.cc.

Referenced by analyze().

const edm::ParameterSet& TrackerRemapper::iConfig

private

Definition at line 106 of file TrackerRemapper.cc.

int TrackerRemapper::m_analyzeMode

private

Definition at line 109 of file TrackerRemapper.cc.

Referenced by analyze(), bookBins(), and TrackerRemapper().

std::map<unsigned, std::string> TrackerRemapper::m_analyzeModeNameMap

private

Definition at line 118 of file TrackerRemapper.cc.

Referenced by bookBins(), and TrackerRemapper().

std::map<long, TGraph*> TrackerRemapper::m_bins

private

Definition at line 111 of file TrackerRemapper.cc.

Referenced by bookBins(), and readVertices().

std::vector<unsigned> TrackerRemapper::m_detIdVector

private

Definition at line 112 of file TrackerRemapper.cc.

Referenced by readVertices().

int TrackerRemapper::m_opMode

private

Definition at line 108 of file TrackerRemapper.cc.

Referenced by analyze(), bookBins(), and TrackerRemapper().

std::string TrackerRemapper::m_pixelBaseDir

private

Definition at line 124 of file TrackerRemapper.cc.

Referenced by preparePixelNames(), and TrackerRemapper().

std::string TrackerRemapper::m_pixelDesiredHistogramBarrel

private

Definition at line 124 of file TrackerRemapper.cc.

Referenced by preparePixelNames(), and TrackerRemapper().

std::string TrackerRemapper::m_pixelDesiredHistogramDisk

private

Definition at line 124 of file TrackerRemapper.cc.

Referenced by preparePixelNames(), and TrackerRemapper().

std::map<unsigned, std::string> TrackerRemapper::m_pixelHistnameMap

private

Definition at line 117 of file TrackerRemapper.cc.

Referenced by fillBarrelRemap(), fillEndcapRemap(), and preparePixelNames().

std::string TrackerRemapper::m_pixelRemapFile

private

Definition at line 121 of file TrackerRemapper.cc.

Referenced by fillPixelRemap(), and TrackerRemapper().

std::string TrackerRemapper::m_stripBaseDir

private

Definition at line 123 of file TrackerRemapper.cc.

Referenced by prepareStripNames(), and TrackerRemapper().

std::string TrackerRemapper::m_stripDesiredHistogram

private

Definition at line 123 of file TrackerRemapper.cc.

Referenced by bookBins(), prepareStripNames(), and TrackerRemapper().

std::map<unsigned, std::string> TrackerRemapper::m_stripHistnameMap

private

Definition at line 116 of file TrackerRemapper.cc.

Referenced by fillStripRemap(), and prepareStripNames().

std::string TrackerRemapper::m_stripRemapFile

private

Definition at line 120 of file TrackerRemapper.cc.

Referenced by fillStripRemap(), and TrackerRemapper().

const TkDetMap* TrackerRemapper::m_tkdetmap

private

Definition at line 114 of file TrackerRemapper.cc.

Referenced by analyze(), and fillStripRemap().

edm::EDGetTokenT<reco::TrackCollection> TrackerRemapper::rechitSrcToken

private

Definition at line 130 of file TrackerRemapper.cc.

Referenced by analyze(), and TrackerRemapper().

std::string TrackerRemapper::runString

private

Definition at line 126 of file TrackerRemapper.cc.

Referenced by bookBins(), and TrackerRemapper().

const edm::ESGetToken<TkDetMap, TrackerTopologyRcd> TrackerRemapper::tkDetMapToken_

private

Definition at line 103 of file TrackerRemapper.cc.

Referenced by analyze().

const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> TrackerRemapper::topoToken_

private

Definition at line 102 of file TrackerRemapper.cc.

Referenced by analyze().

TH2Poly* TrackerRemapper::trackerMap {nullptr}

private

Definition at line 128 of file TrackerRemapper.cc.

Referenced by analyze(), analyzeRechits(), bookBins(), fillBarrelRemap(), fillEndcapRemap(), and fillStripRemap().