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 noexcept final
bool	wantsGlobalRuns () const noexcept final
bool	wantsInputProcessBlocks () const noexcept final
bool	wantsProcessBlocks () const noexcept 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 noexcept
bool	wantsStreamRuns () const noexcept
	~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
ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESResolverIndex > 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
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProductResolverIndices 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_
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 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 ()
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)
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 54 of file TrackerRemapper.cc.

Member Enumeration Documentation

AnalyzeData

enum TrackerRemapper::AnalyzeData

Enumerator
RECHITS
DIGIS
CLUSTERS

Definition at line 73 of file TrackerRemapper.cc.

                   {
    RECHITS = 1,
    DIGIS,
    CLUSTERS,
  };

OpMode

enum TrackerRemapper::OpMode

Enumerator
MODE_ANALYZE
MODE_REMAP

Definition at line 79 of file TrackerRemapper.cc.

{ MODE_ANALYZE = 0, MODE_REMAP = 1 };

PixelLayerEnum

enum TrackerRemapper::PixelLayerEnum

Enumerator
INVALID
PXB_L1
PXB_L2
PXB_L3
PXB_L4
PXF_R1
PXF_R2

Definition at line 61 of file TrackerRemapper.cc.

                      {
    INVALID = 0,

    PXB_L1,
    PXB_L2,
    PXB_L3,
    PXB_L4,

    PXF_R1,
    PXF_R2
  };

Constructor & Destructor Documentation

TrackerRemapper()

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()),
      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::~TrackerRemapper ( )

override

Definition at line 309 of file TrackerRemapper.cc.

{}

Member Function Documentation

analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 312 of file TrackerRemapper.cc.

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

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

analyzeGeneric() [1/2]

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 iEvent, input, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, nano_mu_digi_cff::rawId, and TrackRefitter_38T_cff::src.

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

analyzeGeneric() [2/2]

template<>

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

Definition at line 157 of file TrackerRemapper.cc.

References iEvent, and submitPVValidationJobs::t.

{
  analyzeRechits(iEvent, t);
}

analyzeRechits()

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

private

Definition at line 346 of file TrackerRemapper.cc.

References mps_fire::i, iEvent, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, nano_mu_digi_cff::rawId, rpcPointValidation_cfi::recHit, TrackRefitter_38T_cff::src, HLT_2024v14_cff::track, trackerMap, and DiMuonV_cfg::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);
    }
  }
}

beginJob()

void TrackerRemapper::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 619 of file TrackerRemapper.cc.

{}

bookBins()

void TrackerRemapper::bookBins ( )

private

Definition at line 377 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
}

endJob()

void TrackerRemapper::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 621 of file TrackerRemapper.cc.

{}

fillBarrelRemap()

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

private

Definition at line 537 of file TrackerRemapper.cc.

References TrackerGeometry::detsPXB(), GeomDet::geographicalId(), mps_fire::i, PixelBarrelName::ladderName(), nano_mu_digi_cff::layer, m_pixelHistnameMap, DetId::rawId(), PixelBarrelName::shell(), trackerMap, and heppy_batch::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);
  }
}

fillDescriptions()

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

static

Definition at line 624 of file TrackerRemapper.cc.

References edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, ProducerED_cfi::InputTag, and or.

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

fillEndcapRemap()

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

private

Definition at line 577 of file TrackerRemapper.cc.

References PixelEndcapName::bladeName(), TrackerGeometry::detsPXF(), GeomDet::geographicalId(), PixelEndcapName::halfCylinder(), mps_fire::i, nano_mu_digi_cff::layer, m_pixelHistnameMap, PixelEndcapName::pannelName(), PixelEndcapName::ringName(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::side(), trackerMap, and heppy_batch::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);
  }
}

fillPixelRemap()

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

private

Definition at line 521 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();
}

fillStripRemap()

void TrackerRemapper::fillStripRemap ( )

private

Definition at line 481 of file TrackerRemapper.cc.

References TkDetMap::getComponents(), TkDetMap::getDetFromBin(), nano_mu_digi_cff::layer, m_stripHistnameMap, m_stripRemapFile, m_tkdetmap, trackerMap, and heppy_batch::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();
}

preparePixelNames()

void TrackerRemapper::preparePixelNames ( )

private

Definition at line 291 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";
}

prepareStripNames()

void TrackerRemapper::prepareStripNames ( )

private

Definition at line 213 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";
}

readVertices()

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

private

Definition at line 415 of file TrackerRemapper.cc.

References mps_check::array, hcalRecHitTable_cff::detId, Exception, contentValuesFiles::fullPath, mps_fire::i, recoMuon::in, fastTrackerRecHitType::isPixel(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::iy(), mps_splice::line, m_bins, m_detIdVector, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, AlCaHLTBitMon_QueryRunRegistry::string, 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());

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

    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

clusterSrcToken

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

private

Definition at line 132 of file TrackerRemapper.cc.

Referenced by analyze(), and TrackerRemapper().

digiSrcToken

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

private

Definition at line 131 of file TrackerRemapper.cc.

Referenced by analyze(), and TrackerRemapper().

fs

edm::Service<TFileService> TrackerRemapper::fs

private

Definition at line 106 of file TrackerRemapper.cc.

Referenced by bookBins().

geomToken_

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

private

Definition at line 102 of file TrackerRemapper.cc.

Referenced by analyze().

m_analyzeMode

int TrackerRemapper::m_analyzeMode

private

Definition at line 109 of file TrackerRemapper.cc.

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

m_analyzeModeNameMap

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

private

Definition at line 118 of file TrackerRemapper.cc.

Referenced by bookBins(), and TrackerRemapper().

m_bins

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

private

Definition at line 111 of file TrackerRemapper.cc.

Referenced by bookBins(), and readVertices().

m_detIdVector

std::vector<unsigned> TrackerRemapper::m_detIdVector

private

Definition at line 112 of file TrackerRemapper.cc.

Referenced by readVertices().

m_opMode

int TrackerRemapper::m_opMode

private

Definition at line 108 of file TrackerRemapper.cc.

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

m_pixelBaseDir

std::string TrackerRemapper::m_pixelBaseDir

private

Definition at line 124 of file TrackerRemapper.cc.

Referenced by preparePixelNames(), and TrackerRemapper().

m_pixelDesiredHistogramBarrel

std::string TrackerRemapper::m_pixelDesiredHistogramBarrel

private

Definition at line 124 of file TrackerRemapper.cc.

Referenced by preparePixelNames(), and TrackerRemapper().

m_pixelDesiredHistogramDisk

std::string TrackerRemapper::m_pixelDesiredHistogramDisk

private

Definition at line 124 of file TrackerRemapper.cc.

Referenced by preparePixelNames(), and TrackerRemapper().

m_pixelHistnameMap

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

private

Definition at line 117 of file TrackerRemapper.cc.

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

m_pixelRemapFile

std::string TrackerRemapper::m_pixelRemapFile

private

Definition at line 121 of file TrackerRemapper.cc.

Referenced by fillPixelRemap(), and TrackerRemapper().

m_stripBaseDir

std::string TrackerRemapper::m_stripBaseDir

private

Definition at line 123 of file TrackerRemapper.cc.

Referenced by prepareStripNames(), and TrackerRemapper().

m_stripDesiredHistogram

std::string TrackerRemapper::m_stripDesiredHistogram

private

Definition at line 123 of file TrackerRemapper.cc.

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

m_stripHistnameMap

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

private

Definition at line 116 of file TrackerRemapper.cc.

Referenced by fillStripRemap(), and prepareStripNames().

m_stripRemapFile

std::string TrackerRemapper::m_stripRemapFile

private

Definition at line 120 of file TrackerRemapper.cc.

Referenced by fillStripRemap(), and TrackerRemapper().

m_tkdetmap

const TkDetMap* TrackerRemapper::m_tkdetmap

private

Definition at line 114 of file TrackerRemapper.cc.

Referenced by analyze(), and fillStripRemap().

rechitSrcToken

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

private

Definition at line 130 of file TrackerRemapper.cc.

Referenced by analyze(), and TrackerRemapper().

runString

std::string TrackerRemapper::runString

private

Definition at line 126 of file TrackerRemapper.cc.

Referenced by bookBins(), and TrackerRemapper().

tkDetMapToken_

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

private

Definition at line 104 of file TrackerRemapper.cc.

Referenced by analyze().

topoToken_

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

private

Definition at line 103 of file TrackerRemapper.cc.

Referenced by analyze().

trackerMap

TH2Poly* TrackerRemapper::trackerMap {nullptr}

private

Definition at line 128 of file TrackerRemapper.cc.

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