SiStripGainCalibTableProducer Class Reference

Inheritance diagram for SiStripGainCalibTableProducer:

Public Member Functions
void	fillTable (const std::vector< OnTrackCluster > &clusters, const edm::View< reco::Track > &tracks, nanoaod::FlatTable *table, const edm::EventSetup &iSetup) final
	SiStripGainCalibTableProducer (const edm::ParameterSet &params)
Public Member Functions inherited from SiStripOnTrackClusterTableProducerBase
void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) final
	SiStripOnTrackClusterTableProducerBase (const edm::ParameterSet &params)
	~SiStripOnTrackClusterTableProducerBase () override
Public Member Functions inherited from edm::stream::EDProducer<>
	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

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from SiStripOnTrackClusterTableProducerBase
template<typename VALUES >
static void	addColumn (nanoaod::FlatTable *table, const std::string &name, VALUES &&values, const std::string &doc)

Private Member Functions
double	thickness (DetId id, const TrackerGeometry *tGeom)

Private Attributes
const edm::ESGetToken< SiStripGain, SiStripGainRcd >	m_gainToken
std::map< DetId, double >	m_thicknessMap
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	m_tkGeomToken

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
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

Detailed Description

Definition at line 15 of file SiStripGainCalibTableProducer.cc.

Constructor & Destructor Documentation

SiStripGainCalibTableProducer()

SiStripGainCalibTableProducer::SiStripGainCalibTableProducer ( const edm::ParameterSet &  params )

inlineexplicit

Definition at line 17 of file SiStripGainCalibTableProducer.cc.

      : SiStripOnTrackClusterTableProducerBase(params), m_tkGeomToken{esConsumes<>()}, m_gainToken{esConsumes<>()} {}

Member Function Documentation

fillDescriptions()

static void SiStripGainCalibTableProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

inlinestatic

Definition at line 20 of file SiStripGainCalibTableProducer.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                           {
    edm::ParameterSetDescription desc;
    desc.add<std::string>("name", "cluster");
    desc.add<std::string>("doc", "");
    desc.add<bool>("extension", false);
    desc.add<edm::InputTag>("Tracks", edm::InputTag{"generalTracks"});
    descriptions.add("siStripGainCalibTable", desc);
  }

fillTable()

void SiStripGainCalibTableProducer::fillTable ( const std::vector< OnTrackCluster > &  clusters, const edm::View< reco::Track > &  tracks, nanoaod::FlatTable *  table, const edm::EventSetup &  iSetup  )

finalvirtual

Implements SiStripOnTrackClusterTableProducerBase.

Definition at line 93 of file SiStripGainCalibTableProducer.cc.

References funct::abs(), SiStripOnTrackClusterTableProducerBase::addColumn(), ALCARECOTkAlJpsiMuMu_cff::charge, bsc_activity_cfg::clusters, SiStripGain::getApvGain(), edm::EventSetup::getHandle(), SiStripGain::getRange(), deDxTools::isFarFromBorder(), edm::ESHandleBase::isValid(), TrajectoryStateOnSurface::localDirection(), m_gainToken, m_tkGeomToken, me0TriggerPseudoDigis_cff::nStrips, EnsembleCalibrationLA_cfg::path, edm::ESHandle< T >::product(), EcalSCDynamicDPhiParametersESProducer_cfi::saturation, TableParser::table, thickness(), and PV3DBase< T, PVType, FrameType >::z().

                                                                           {
  edm::ESHandle<TrackerGeometry> tGeom = iSetup.getHandle(m_tkGeomToken);
  edm::ESHandle<SiStripGain> stripGains = iSetup.getHandle(m_gainToken);

  std::vector<double> c_localdirx;
  std::vector<double> c_localdiry;
  std::vector<double> c_localdirz;
  std::vector<uint16_t> c_firststrip;
  std::vector<uint16_t> c_nstrips;
  std::vector<bool> c_saturation;
  std::vector<bool> c_overlapping;
  std::vector<bool> c_farfromedge;
  std::vector<int> c_charge;
  std::vector<double> c_path;
  // NOTE only very few types are supported by NanoAOD, but more could be added (to discuss with XPOG / core software)
  // NOTE removed amplitude vector, I don't think it was used anywhere
  std::vector<float> c_gainused;      // NOTE was double
  std::vector<float> c_gainusedTick;  // NOTE was double
  for (const auto clus : clusters) {
    const auto& ampls = clus.cluster->amplitudes();
    const int firstStrip = clus.cluster->firstStrip();
    const int nStrips = ampls.size();
    double prevGain = -1;
    double prevGainTick = -1;
    if (stripGains.isValid()) {
      prevGain = stripGains->getApvGain(firstStrip / 128, stripGains->getRange(clus.det, 1), 1);
      prevGainTick = stripGains->getApvGain(firstStrip / 128, stripGains->getRange(clus.det, 0), 1);
    }
    const unsigned int charge = clus.cluster->charge();
    const bool saturation = std::any_of(ampls.begin(), ampls.end(), [](uint8_t amp) { return amp >= 254; });

    const bool overlapping = (((firstStrip % 128) == 0) || ((firstStrip / 128) != ((firstStrip + nStrips) / 128)) ||
                              (((firstStrip + nStrips) % 128) == 127));
    const auto& trajState = clus.measurement.updatedState();
    const auto trackDir = trajState.localDirection();
    const auto cosine = trackDir.z() / trackDir.mag();
    const auto path = (10. * thickness(clus.det, tGeom.product())) / std::abs(cosine);
    const auto farFromEdge = isFarFromBorder(trajState, clus.det, tGeom.product());
    c_localdirx.push_back(trackDir.x());
    c_localdiry.push_back(trackDir.y());
    c_localdirz.push_back(trackDir.z());
    c_firststrip.push_back(firstStrip);
    c_nstrips.push_back(nStrips);
    c_saturation.push_back(saturation);
    c_overlapping.push_back(overlapping);
    c_farfromedge.push_back(farFromEdge);
    c_charge.push_back(charge);
    c_path.push_back(path);
    c_gainused.push_back(prevGain);
    c_gainusedTick.push_back(prevGainTick);
  }
  // addColumn(table, "localdirx", c_localdirx, "<doc>");
  // addColumn(table, "localdiry", c_localdiry, "<doc>");
  // addColumn(table, "localdirz", c_localdirz, "<doc>");
  // addColumn(table, "firststrip", c_firststrip, "<doc>");
  // addColumn(table, "nstrips", c_nstrips, "<doc>");
  addColumn(table, "saturation", c_saturation, "<doc>");
  addColumn(table, "overlapping", c_overlapping, "<doc>");
  addColumn(table, "farfromedge", c_farfromedge, "<doc>");
  addColumn(table, "charge", c_charge, "<doc>");
  // addColumn(table, "path", c_path, "<doc>");
  // ExtendedCalibTree: also charge/path
  addColumn(table, "gainused", c_gainused, "<doc>");
  addColumn(table, "gainusedTick", c_gainusedTick, "<doc>");
}

thickness()

double SiStripGainCalibTableProducer::thickness ( DetId  id, const TrackerGeometry *  tGeom  )

private

Definition at line 73 of file SiStripGainCalibTableProducer.cc.

References EcalPhiSymFlatTableProducers_cfi::id, TrackerGeometry::idToDetUnit(), fastTrackerRecHitType::isPixel(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, and m_thicknessMap.

Referenced by fillTable().

                                                                                      {
  const auto it = m_thicknessMap.find(id);
  if (m_thicknessMap.end() != it) {
    return it->second;
  } else {
    double detThickness = 1.;
    const auto gdu = tGeom->idToDetUnit(id);
    const auto isPixel = (dynamic_cast<const PixelGeomDetUnit*>(gdu) != nullptr);
    const auto isStrip = (dynamic_cast<const StripGeomDetUnit*>(gdu) != nullptr);
    if (!isPixel && !isStrip) {
      edm::LogWarning("SiStripGainCalibTableProducer")
          << "DetId " << id.rawId() << " doesn't seem to belong to the Tracker";
    } else {
      detThickness = gdu->surface().bounds().thickness();
    }
    m_thicknessMap[id] = detThickness;
    return detThickness;
  }
}

Member Data Documentation

m_gainToken

const edm::ESGetToken<SiStripGain, SiStripGainRcd> SiStripGainCalibTableProducer::m_gainToken

private

Definition at line 36 of file SiStripGainCalibTableProducer.cc.

Referenced by fillTable().

m_thicknessMap

std::map<DetId, double> SiStripGainCalibTableProducer::m_thicknessMap

private

Definition at line 38 of file SiStripGainCalibTableProducer.cc.

Referenced by thickness().

m_tkGeomToken

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> SiStripGainCalibTableProducer::m_tkGeomToken

private

Definition at line 35 of file SiStripGainCalibTableProducer.cc.

Referenced by fillTable().