PreMixingSiPixelWorker Class Reference

Inheritance diagram for PreMixingSiPixelWorker:

Public Member Functions
void	addPileups (PileUpEventPrincipal const &, edm::EventSetup const &es) override
void	addSignals (edm::Event const &e, edm::EventSetup const &es) override
void	initializeEvent (edm::Event const &e, edm::EventSetup const &c) override
	PreMixingSiPixelWorker (const edm::ParameterSet &ps, edm::ProducesCollector, edm::ConsumesCollector &&iC)
void	put (edm::Event &e, edm::EventSetup const &iSetup, std::vector< PileupSummaryInfo > const &ps, int bs) override
	~PreMixingSiPixelWorker () override=default
Public Member Functions inherited from PreMixingWorker
virtual void	beginLuminosityBlock (edm::LuminosityBlock const &iLumi, edm::EventSetup const &iSetup)
virtual void	beginRun (edm::Run const &iRun, edm::EventSetup const &iSetup)
virtual void	endRun ()
virtual void	finalizeBunchCrossing (edm::Event &iEvent, edm::EventSetup const &iSetup, int bunchCrossing)
virtual void	initializeBunchCrossing (edm::Event const &iEvent, edm::EventSetup const &iSetup, int bunchCrossing)
	PreMixingWorker ()=default
virtual	~PreMixingWorker ()=default

Private Types
typedef int	Amplitude
typedef std::multimap< int, PixelDigi >	OneDetectorMap
typedef std::multimap< int, PixelSimHitExtraInfo >	OneExtraInfoMap
typedef std::map< uint32_t, OneDetectorMap >	SiGlobalIndex
typedef std::map< int, Amplitude, std::less< int > >	signal_map_type
typedef std::map< uint32_t, signal_map_type >	signalMaps
typedef std::map< uint32_t, OneExtraInfoMap >	SiPixelExtraInfo

Private Attributes
bool	applyLateReweighting_
SiPixelDigitizerAlgorithm	digitizer_
bool	firstFinalizeEvent_ = true
bool	firstInitializeEvent_ = true
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	pDDToken_
edm::InputTag	pixeldigi_collectionPile_
edm::InputTag	pixeldigi_collectionSig_
edm::InputTag	pixeldigi_extraInfo_
std::string	PixelDigiCollectionDM_
edm::EDGetTokenT< edm::DetSetVector< PixelSimHitExtraInfo > >	PixelDigiPExtraToken_
edm::EDGetTokenT< edm::DetSetVector< PixelDigi > >	PixelDigiPToken_
edm::EDGetTokenT< edm::DetSetVector< PixelDigi > >	PixelDigiToken_
SiPixelExtraInfo	SiHitExtraStorage_
SiGlobalIndex	SiHitStorage_
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	tTopoToken_

Detailed Description

Definition at line 40 of file PreMixingSiPixelWorker.cc.

Member Typedef Documentation

Amplitude

typedef int PreMixingSiPixelWorker::Amplitude

private

Definition at line 67 of file PreMixingSiPixelWorker.cc.

OneDetectorMap

typedef std::multimap<int, PixelDigi> PreMixingSiPixelWorker::OneDetectorMap

private

Definition at line 72 of file PreMixingSiPixelWorker.cc.

OneExtraInfoMap

typedef std::multimap<int, PixelSimHitExtraInfo> PreMixingSiPixelWorker::OneExtraInfoMap

private

Definition at line 74 of file PreMixingSiPixelWorker.cc.

SiGlobalIndex

typedef std::map<uint32_t, OneDetectorMap> PreMixingSiPixelWorker::SiGlobalIndex

private

Definition at line 73 of file PreMixingSiPixelWorker.cc.

signal_map_type

typedef std::map<int, Amplitude, std::less<int> > PreMixingSiPixelWorker::signal_map_type

private

Definition at line 68 of file PreMixingSiPixelWorker.cc.

signalMaps

typedef std::map<uint32_t, signal_map_type> PreMixingSiPixelWorker::signalMaps

private

Definition at line 69 of file PreMixingSiPixelWorker.cc.

SiPixelExtraInfo

typedef std::map<uint32_t, OneExtraInfoMap> PreMixingSiPixelWorker::SiPixelExtraInfo

private

Definition at line 75 of file PreMixingSiPixelWorker.cc.

Constructor & Destructor Documentation

PreMixingSiPixelWorker()

PreMixingSiPixelWorker::PreMixingSiPixelWorker	(	const edm::ParameterSet &	ps,
		edm::ProducesCollector	producesCollector,
		edm::ConsumesCollector &&	iC
	)

Definition at line 86 of file PreMixingSiPixelWorker.cc.

References applyLateReweighting_, edm::ParameterSet::getParameter(), LogDebug, pixeldigi_collectionPile_, pixeldigi_collectionSig_, pixeldigi_extraInfo_, PixelDigiCollectionDM_, PixelDigiPExtraToken_, PixelDigiPToken_, PixelDigiToken_, edm::ProducesCollector::produces(), SiHitExtraStorage_, SiHitStorage_, and AlCaHLTBitMon_QueryRunRegistry::string.

    : tTopoToken_(iC.esConsumes()),
      pDDToken_(iC.esConsumes(edm::ESInputTag("", ps.getParameter<std::string>("PixGeometryType")))),
      digitizer_(ps, iC) {
  // declare the products to produce

  pixeldigi_collectionSig_ = ps.getParameter<edm::InputTag>("pixeldigiCollectionSig");
  pixeldigi_collectionPile_ = ps.getParameter<edm::InputTag>("pixeldigiCollectionPile");
  pixeldigi_extraInfo_ = ps.getParameter<edm::InputTag>("pixeldigiExtraCollectionPile");
  PixelDigiCollectionDM_ = ps.getParameter<std::string>("PixelDigiCollectionDM");
  applyLateReweighting_ = ps.getParameter<bool>("applyLateReweighting");
  LogDebug("PreMixingSiPixelWorker") << "applyLateReweighting_ in PreMixingSiPixelWorker  " << applyLateReweighting_;

  PixelDigiToken_ = iC.consumes<edm::DetSetVector<PixelDigi>>(pixeldigi_collectionSig_);
  PixelDigiPToken_ = iC.consumes<edm::DetSetVector<PixelDigi>>(pixeldigi_collectionPile_);
  PixelDigiPExtraToken_ = iC.consumes<edm::DetSetVector<PixelSimHitExtraInfo>>(pixeldigi_extraInfo_);

  producesCollector.produces<edm::DetSetVector<PixelDigi>>(PixelDigiCollectionDM_);
  producesCollector.produces<PixelFEDChannelCollection>(PixelDigiCollectionDM_);

  // clear local storage for this event
  SiHitStorage_.clear();
  SiHitExtraStorage_.clear();
}

~PreMixingSiPixelWorker()

PreMixingSiPixelWorker::~PreMixingSiPixelWorker ( )

overridedefault

Member Function Documentation

addPileups()

void PreMixingSiPixelWorker::addPileups ( PileUpEventPrincipal const &  pep, edm::EventSetup const &  es  )

overridevirtual

Implements PreMixingWorker.

Definition at line 154 of file PreMixingSiPixelWorker.cc.

References applyLateReweighting_, edm::DetSetVector< T >::begin(), PileUpEventPrincipal::bunchCrossing(), digitizer_, edm::DetSetVector< T >::end(), mps_fire::end, Exception, PileUpEventPrincipal::getByLabel(), edm::EventSetup::getData(), edm::RandomNumberGenerator::getEngine(), edm::EventPrincipal::id(), input, edm::DetSetVector< T >::insert(), edm::HandleBase::isValid(), SiPixelDigitizerAlgorithm::lateSignalReweight(), LogDebug, pDDToken_, pixeldigi_collectionPile_, pixeldigi_extraInfo_, PileUpEventPrincipal::principal(), SiHitExtraStorage_, SiHitStorage_, edm::EventPrincipal::streamID(), and tTopoToken_.

                                                                                                {
  LogDebug("PreMixingSiPixelWorker") << "\n===============> adding pileups from event  " << pep.principal().id()
                                     << " for bunchcrossing " << pep.bunchCrossing();

  // fill in maps of hits; same code as addSignals, except now applied to the pileup events

  edm::Handle<edm::DetSetVector<PixelDigi>> inputHandle;
  pep.getByLabel(pixeldigi_collectionPile_, inputHandle);

  // added for the Late CR
  edm::Handle<edm::DetSetVector<PixelSimHitExtraInfo>> pixelAddInfo;
  pep.getByLabel(pixeldigi_extraInfo_, pixelAddInfo);
  const TrackerTopology* tTopo = &es.getData(tTopoToken_);
  auto const& pDD = es.getData(pDDToken_);
  edm::Service<edm::RandomNumberGenerator> rng;
  CLHEP::HepRandomEngine* engine = &rng->getEngine(pep.principal().streamID());

  if (inputHandle.isValid()) {
    const auto& input = *inputHandle;

    bool loadExtraInformation = false;

    if (pixelAddInfo.isValid() && applyLateReweighting_) {
      // access the extra information
      loadExtraInformation = true;
      // Iterate on detector units
      edm::DetSetVector<PixelSimHitExtraInfo>::const_iterator detIdIter;
      for (detIdIter = pixelAddInfo->begin(); detIdIter != pixelAddInfo->end(); detIdIter++) {
        uint32_t detid = detIdIter->id;  // = rawid
        OneExtraInfoMap LocalExtraMap;
        edm::DetSet<PixelSimHitExtraInfo>::const_iterator di;
        for (di = detIdIter->data.begin(); di != detIdIter->data.end(); di++) {
          LocalExtraMap.insert(OneExtraInfoMap::value_type((di->hitIndex()), *di));
        }
        SiHitExtraStorage_.insert(SiPixelExtraInfo::value_type(detid, LocalExtraMap));
      }  // end loop on detIdIter
    }    // end if applyLateReweighting_
    else if (!pixelAddInfo.isValid() && applyLateReweighting_) {
      edm::LogError("PreMixingSiPixelWorker") << " Problem in accessing the Extra Pixel SimHit Collection  !!!! ";
      edm::LogError("PreMixingSiPixelWorker") << " The Late Charge Reweighting can not be applied ";
      throw cms::Exception("PreMixingSiPixelWorker")
          << " Problem in accessing the Extra Pixel SimHit Collection for Late Charge Reweighting \n";
    }

    //loop on all detsets (detectorIDs) inside the input collection
    edm::DetSetVector<PixelDigi>::const_iterator DSViter = input.begin();
    for (; DSViter != input.end(); DSViter++) {
#ifdef DEBUG
      LogDebug("PreMixingSiPixelWorker") << "Pileups: Processing DetID " << DSViter->id;
#endif

      uint32_t detID = DSViter->id;
      edm::DetSet<PixelDigi>::const_iterator begin = (DSViter->data).begin();
      edm::DetSet<PixelDigi>::const_iterator end = (DSViter->data).end();
      edm::DetSet<PixelDigi>::const_iterator icopy;

      // find correct local map (or new one) for this detector ID

      SiGlobalIndex::const_iterator itest;

      itest = SiHitStorage_.find(detID);

      std::vector<PixelDigi> TempDigis;
      for (icopy = begin; icopy != end; icopy++) {
        TempDigis.push_back(*icopy);
      }
      if (loadExtraInformation) {
        // apply the Late Charge Reweighthing on Pile-up digi
        SiPixelExtraInfo::const_iterator jtest;
        jtest = SiHitExtraStorage_.find(detID);
        OneExtraInfoMap LocalSimHitExtraMap = jtest->second;
        std::vector<PixelSimHitExtraInfo> TempSimExtra;
        for (auto& iLocal : LocalSimHitExtraMap) {
          TempSimExtra.push_back(iLocal.second);
        }

        for (const auto& iu : pDD.detUnits()) {
          if (iu->type().isTrackerPixel()) {
            uint32_t detIDinLoop = iu->geographicalId().rawId();
            if (detIDinLoop == detID) {
              digitizer_.lateSignalReweight(
                  dynamic_cast<const PixelGeomDetUnit*>(iu), TempDigis, TempSimExtra, tTopo, engine);
              break;
            }
          }
        }
      }

      if (itest != SiHitStorage_.end()) {  // this detID already has hits, add to existing map

        OneDetectorMap LocalMap = itest->second;

        // fill in local map with extra channels
        for (unsigned int ij = 0; ij < TempDigis.size(); ij++) {
          LocalMap.insert(OneDetectorMap::value_type((TempDigis[ij].channel()), TempDigis[ij]));
        }

        SiHitStorage_[detID] = LocalMap;

      } else {  // fill local storage with this information, put in global collection

        OneDetectorMap LocalMap;

        for (unsigned int ij = 0; ij < TempDigis.size(); ij++) {
          LocalMap.insert(OneDetectorMap::value_type((TempDigis[ij].channel()), TempDigis[ij]));
        }

        SiHitStorage_.insert(SiGlobalIndex::value_type(detID, LocalMap));
      }
    }
  }
}

addSignals()

void PreMixingSiPixelWorker::addSignals ( edm::Event const &  e, edm::EventSetup const &  es  )

overridevirtual

Implements PreMixingWorker.

Definition at line 123 of file PreMixingSiPixelWorker.cc.

References MillePedeFileConverter_cfg::e, mps_fire::end, input, LogDebug, PixelDigiToken_, and SiHitStorage_.

                                                                                  {
  // fill in maps of hits

  LogDebug("PreMixingSiPixelWorker") << "===============> adding MC signals for " << e.id();

  edm::Handle<edm::DetSetVector<PixelDigi>> input;

  if (e.getByToken(PixelDigiToken_, input)) {
    //loop on all detsets (detectorIDs) inside the input collection
    edm::DetSetVector<PixelDigi>::const_iterator DSViter = input->begin();
    for (; DSViter != input->end(); DSViter++) {
#ifdef DEBUG
      LogDebug("PreMixingSiPixelWorker") << "Processing DetID " << DSViter->id;
#endif

      uint32_t detID = DSViter->id;
      edm::DetSet<PixelDigi>::const_iterator begin = (DSViter->data).begin();
      edm::DetSet<PixelDigi>::const_iterator end = (DSViter->data).end();
      edm::DetSet<PixelDigi>::const_iterator icopy;

      OneDetectorMap LocalMap;

      for (icopy = begin; icopy != end; icopy++) {
        LocalMap.insert(OneDetectorMap::value_type((icopy->channel()), *icopy));
      }

      SiHitStorage_.insert(SiGlobalIndex::value_type(detID, LocalMap));
    }
  }
}  // end of addSiPixelSignals

initializeEvent()

void PreMixingSiPixelWorker::initializeEvent ( edm::Event const &  e, edm::EventSetup const &  c  )

overridevirtual

Implements PreMixingWorker.

Definition at line 115 of file PreMixingSiPixelWorker.cc.

References digitizer_, firstInitializeEvent_, SiPixelDigitizerAlgorithm::init(), and SiPixelDigitizerAlgorithm::initializeEvent().

                                                                                           {
  if (firstInitializeEvent_) {
    digitizer_.init(iSetup);
    firstInitializeEvent_ = false;
  }
  digitizer_.initializeEvent();
}

put()

void PreMixingSiPixelWorker::put ( edm::Event &  e, edm::EventSetup const &  iSetup, std::vector< PileupSummaryInfo > const &  ps, int  bs  )

overridevirtual

Implements PreMixingWorker.

Definition at line 267 of file PreMixingSiPixelWorker.cc.

References gpuClustering::adc, cms::cuda::bs, SiPixelDigitizerAlgorithm::calculateInstlumiFactor(), SiPixelDigitizerAlgorithm::chooseScenario(), SiPixelDigitizerAlgorithm::digitize(), digitizer_, MillePedeFileConverter_cfg::e, Exception, firstFinalizeEvent_, edm::EventSetup::getData(), edm::RandomNumberGenerator::getEngine(), SiPixelDigitizerAlgorithm::init_DynIneffDB(), SiPixelDigitizerAlgorithm::killBadFEDChannels(), eostools::move(), pDDToken_, PixelDigiCollectionDM_, SiPixelDigitizerAlgorithm::setSimAccumulator(), SiHitExtraStorage_, SiHitStorage_, and tTopoToken_.

                                         {
  // collection of Digis to put in the event

  std::vector<edm::DetSet<PixelDigi>> vPixelDigi;

  // loop through our collection of detectors, merging hits and putting new ones in the output
  signalMaps signal;

  // big loop over Detector IDs:

  for (SiGlobalIndex::const_iterator IDet = SiHitStorage_.begin(); IDet != SiHitStorage_.end(); IDet++) {
    uint32_t detID = IDet->first;

    OneDetectorMap LocalMap = IDet->second;

    signal_map_type Signals;
    Signals.clear();

    //counter variables
    int formerPixel = -1;
    int currentPixel;
    int ADCSum = 0;

    OneDetectorMap::const_iterator iLocalchk;

    for (OneDetectorMap::const_iterator iLocal = LocalMap.begin(); iLocal != LocalMap.end(); ++iLocal) {
      currentPixel = iLocal->first;

      if (currentPixel == formerPixel) {  // we have to add these digis together
        ADCSum += (iLocal->second).adc();
      } else {
        if (formerPixel != -1) {  // ADC info stolen from SiStrips...
          if (ADCSum > 511)
            ADCSum = 255;
          else if (ADCSum > 253 && ADCSum < 512)
            ADCSum = 254;

          Signals.insert(std::make_pair(formerPixel, ADCSum));
        }
        // save pointers for next iteration
        formerPixel = currentPixel;
        ADCSum = (iLocal->second).adc();
      }

      iLocalchk = iLocal;
      if ((++iLocalchk) == LocalMap.end()) {  //make sure not to lose the last one
        if (ADCSum > 511)
          ADCSum = 255;
        else if (ADCSum > 253 && ADCSum < 512)
          ADCSum = 254;
        Signals.insert(std::make_pair(formerPixel, ADCSum));
      }

    }  // end of loop over one detector

    // stick this into the global vector of detector info
    signal.insert(std::make_pair(detID, Signals));

  }  // end of big loop over all detector IDs

  // put the collection of digis in the event
  edm::LogInfo("PreMixingSiPixelWorker") << "total # Merged Pixels: " << signal.size();

  std::vector<edm::DetSet<PixelDigi>> theDigiVector;

  // Load inefficiency constants (1st pass), set pileup information.
  if (firstFinalizeEvent_) {
    digitizer_.init_DynIneffDB(iSetup);
    firstFinalizeEvent_ = false;
  }

  digitizer_.calculateInstlumiFactor(ps, bs);
  digitizer_.setSimAccumulator(signal);

  edm::Service<edm::RandomNumberGenerator> rng;
  CLHEP::HepRandomEngine* engine = &rng->getEngine(e.streamID());

  auto const& pDD = iSetup.getData(pDDToken_);
  const TrackerTopology* tTopo = &iSetup.getData(tTopoToken_);

  if (digitizer_.killBadFEDChannels()) {
    std::unique_ptr<PixelFEDChannelCollection> PixelFEDChannelCollection_ = digitizer_.chooseScenario(ps, engine);
    if (PixelFEDChannelCollection_ == nullptr) {
      throw cms::Exception("NullPointerError") << "PixelFEDChannelCollection not set in chooseScenario function.\n";
    }
    e.put(std::move(PixelFEDChannelCollection_), PixelDigiCollectionDM_);
  }

  for (const auto& iu : pDD.detUnits()) {
    if (iu->type().isTrackerPixel()) {
      edm::DetSet<PixelDigi> collector(iu->geographicalId().rawId());
      edm::DetSet<PixelDigiSimLink> linkcollector(
          iu->geographicalId().rawId());                // ignored as DigiSimLinks are combined separately
      std::vector<PixelDigiAddTempInfo> tempcollector;  // to be ignored ?

      digitizer_.digitize(
          dynamic_cast<const PixelGeomDetUnit*>(iu), collector.data, linkcollector.data, tempcollector, tTopo, engine);
      if (!collector.data.empty()) {
        theDigiVector.push_back(std::move(collector));
      }
    }
  }

  e.put(std::make_unique<edm::DetSetVector<PixelDigi>>(theDigiVector), PixelDigiCollectionDM_);

  // clear local storage for this event
  SiHitStorage_.clear();
  SiHitExtraStorage_.clear();
}

Member Data Documentation

applyLateReweighting_

bool PreMixingSiPixelWorker::applyLateReweighting_

private

Definition at line 82 of file PreMixingSiPixelWorker.cc.

Referenced by addPileups(), and PreMixingSiPixelWorker().

digitizer_

SiPixelDigitizerAlgorithm PreMixingSiPixelWorker::digitizer_

private

Definition at line 62 of file PreMixingSiPixelWorker.cc.

Referenced by addPileups(), initializeEvent(), and put().

firstFinalizeEvent_

bool PreMixingSiPixelWorker::firstFinalizeEvent_ = true

private

Definition at line 81 of file PreMixingSiPixelWorker.cc.

Referenced by put().

firstInitializeEvent_

bool PreMixingSiPixelWorker::firstInitializeEvent_ = true

private

Definition at line 80 of file PreMixingSiPixelWorker.cc.

Referenced by initializeEvent().

pDDToken_

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> PreMixingSiPixelWorker::pDDToken_

private

Definition at line 60 of file PreMixingSiPixelWorker.cc.

Referenced by addPileups(), and put().

pixeldigi_collectionPile_

edm::InputTag PreMixingSiPixelWorker::pixeldigi_collectionPile_

private

Definition at line 52 of file PreMixingSiPixelWorker.cc.

Referenced by addPileups(), and PreMixingSiPixelWorker().

pixeldigi_collectionSig_

edm::InputTag PreMixingSiPixelWorker::pixeldigi_collectionSig_

private

Definition at line 51 of file PreMixingSiPixelWorker.cc.

Referenced by PreMixingSiPixelWorker().

pixeldigi_extraInfo_

edm::InputTag PreMixingSiPixelWorker::pixeldigi_extraInfo_

private

Definition at line 53 of file PreMixingSiPixelWorker.cc.

Referenced by addPileups(), and PreMixingSiPixelWorker().

PixelDigiCollectionDM_

std::string PreMixingSiPixelWorker::PixelDigiCollectionDM_

private

Definition at line 54 of file PreMixingSiPixelWorker.cc.

Referenced by PreMixingSiPixelWorker(), and put().

PixelDigiPExtraToken_

edm::EDGetTokenT<edm::DetSetVector<PixelSimHitExtraInfo> > PreMixingSiPixelWorker::PixelDigiPExtraToken_

private

Definition at line 58 of file PreMixingSiPixelWorker.cc.

Referenced by PreMixingSiPixelWorker().

PixelDigiPToken_

edm::EDGetTokenT<edm::DetSetVector<PixelDigi> > PreMixingSiPixelWorker::PixelDigiPToken_

private

Definition at line 57 of file PreMixingSiPixelWorker.cc.

Referenced by PreMixingSiPixelWorker().

PixelDigiToken_

edm::EDGetTokenT<edm::DetSetVector<PixelDigi> > PreMixingSiPixelWorker::PixelDigiToken_

private

Definition at line 56 of file PreMixingSiPixelWorker.cc.

Referenced by addSignals(), and PreMixingSiPixelWorker().

SiHitExtraStorage_

SiPixelExtraInfo PreMixingSiPixelWorker::SiHitExtraStorage_

private

Definition at line 78 of file PreMixingSiPixelWorker.cc.

Referenced by addPileups(), PreMixingSiPixelWorker(), and put().

SiHitStorage_

SiGlobalIndex PreMixingSiPixelWorker::SiHitStorage_

private

Definition at line 77 of file PreMixingSiPixelWorker.cc.

Referenced by addPileups(), addSignals(), PreMixingSiPixelWorker(), and put().

tTopoToken_

const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> PreMixingSiPixelWorker::tTopoToken_

private

Definition at line 59 of file PreMixingSiPixelWorker.cc.

Referenced by addPileups(), and put().