CTPPSPixelDigiProducer Class Reference

#include <SimPPS/PPSPixelDigiProducer/plugins/PPSPixelDigiProducer.cc>

Inheritance diagram for CTPPSPixelDigiProducer:

Public Member Functions
	CTPPSPixelDigiProducer (const edm::ParameterSet &)
	~CTPPSPixelDigiProducer () 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)

Private Types
typedef std::map< unsigned int, std::vector< PSimHit > >	simhit_map
typedef simhit_map::iterator	simhit_map_iterator

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

Private Attributes
edm::ParameterSet	conf_
edm::ESGetToken < CTPPSPixelGainCalibrations, CTPPSPixelGainCalibrationsRcd >	gainCalibESToken_
CLHEP::HepRandomEngine *	rndEngine_ = nullptr
std::vector< std::string >	RPix_hit_containers_
std::map< uint32_t, std::unique_ptr < RPixDetDigitizer > >	theAlgoMap
edm::ESGetToken < PPSPixelTopology, PPSPixelTopologyRcd >	thePixelTopologyToken_
edm::EDGetTokenT < CrossingFrame< PSimHit > >	tokenCrossingFramePPSPixel
int	verbosity_

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

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

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

Definition at line 79 of file PPSPixelDigiProducer.cc.

Member Typedef Documentation

typedef std::map<unsigned int, std::vector<PSimHit> > CTPPSPixelDigiProducer::simhit_map

private

Definition at line 90 of file PPSPixelDigiProducer.cc.

typedef simhit_map::iterator CTPPSPixelDigiProducer::simhit_map_iterator

private

Definition at line 91 of file PPSPixelDigiProducer.cc.

Constructor & Destructor Documentation

CTPPSPixelDigiProducer::CTPPSPixelDigiProducer ( const edm::ParameterSet &  conf )

explicit

Definition at line 105 of file PPSPixelDigiProducer.cc.

References edm::ParameterSet::getParameter(), HLT_FULL_cff::InputTag, RPix_hit_containers_, tokenCrossingFramePPSPixel, and verbosity_.

    : conf_(conf), gainCalibESToken_(esConsumes()), thePixelTopologyToken_(esConsumes()) {
  produces<edm::DetSetVector<CTPPSPixelDigi>>();

  // register data to consume
  tokenCrossingFramePPSPixel = consumes<CrossingFrame<PSimHit>>(edm::InputTag("mix", "g4SimHitsCTPPSPixelHits"));

  RPix_hit_containers_.clear();
  RPix_hit_containers_ = conf.getParameter<std::vector<std::string>>("ROUList");
  verbosity_ = conf.getParameter<int>("RPixVerbosity");
}

CTPPSPixelDigiProducer::~CTPPSPixelDigiProducer ( )

override

Definition at line 117 of file PPSPixelDigiProducer.cc.

{}

Member Function Documentation

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

static

Definition at line 119 of file PPSPixelDigiProducer.cc.

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

                                                                                        {
  edm::ParameterSetDescription desc;
  // all distances in [mm]
  // RPDigiProducer
  desc.add<std::vector<std::string>>("ROUList", {"CTPPSPixelHits"});
  desc.add<int>("RPixVerbosity", 0);
  desc.add<bool>("CTPPSPixelDigiSimHitRelationsPersistence",
                 false);  // save links betweend digi, clusters and OSCAR/Geant4 hits

  // RPDetDigitizer
  desc.add<double>("RPixEquivalentNoiseCharge", 1000.0);
  desc.add<bool>("RPixNoNoise", false);

  // RPDisplacementGenerator
  desc.add<double>("RPixGeVPerElectron", 3.61e-09);
  desc.add<std::vector<double>>("RPixInterSmearing", {0.011});
  desc.add<bool>("RPixLandauFluctuations", true);
  desc.add<int>("RPixChargeDivisions", 20);
  desc.add<double>("RPixDeltaProductionCut", 0.120425);  // [MeV]

  // RPixChargeShare
  desc.add<std::string>("ChargeMapFile2E", "SimPPS/PPSPixelDigiProducer/data/PixelChargeMap.txt");
  desc.add<std::string>("ChargeMapFile2E_2X", "SimPPS/PPSPixelDigiProducer/data/PixelChargeMap_2X.txt");
  desc.add<std::string>("ChargeMapFile2E_2Y", "SimPPS/PPSPixelDigiProducer/data/PixelChargeMap_2Y.txt");
  desc.add<std::string>("ChargeMapFile2E_2X2Y", "SimPPS/PPSPixelDigiProducer/data/PixelChargeMap_2X2Y.txt");
  desc.add<double>(
      "RPixCoupling",
      0.250);  // fraction of the remaining charge going to the closer neighbour pixel. Value = 0.135, Value = 0.0 bypass the charge map and the charge sharing approach

  // RPixDummyROCSimulator
  desc.add<double>("RPixDummyROCThreshold", 1900.0);
  desc.add<double>("RPixDummyROCElectronPerADC", 135.0);  // 210.0 to be verified
  desc.add<int>("VCaltoElectronGain", 50);                // same values as in RPixDetClusterizer
  desc.add<int>("VCaltoElectronOffset", -411);            //
  desc.add<bool>("doSingleCalibration", false);           //
  desc.add<double>("RPixDeadPixelProbability", 0.001);
  desc.add<bool>("RPixDeadPixelSimulationOn", true);

  desc.add<std::string>("mixLabel", "mix");
  desc.add<std::string>("InputCollection", "g4SimHitsCTPPSPixelHits");
  descriptions.add("RPixDetDigitizer", desc);
}

void CTPPSPixelDigiProducer::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 166 of file PPSPixelDigiProducer.cc.

References MixCollection< T >::MixItr::bunch(), cuy::col, conf_, submitPVResolutionJobs::count, edm::DetSet< T >::data, Exception, gainCalibESToken_, edm::Event::getByToken(), edm::EventSetup::getData(), edm::RandomNumberGenerator::getEngine(), MixCollection< T >::MixItr::getTrigger(), edm::Service< T >::isAvailable(), eostools::move(), edm::Handle< T >::product(), edm::Event::put(), rndEngine_, edm::Event::streamID(), theAlgoMap, thePixelTopologyToken_, tokenCrossingFramePPSPixel, and verbosity_.

                                                                                  {
  if (!rndEngine_) {
    edm::Service<edm::RandomNumberGenerator> rng;
    if (!rng.isAvailable()) {
      throw cms::Exception("Configuration")
          << "This class requires the RandomNumberGeneratorService\n"
             "which is not present in the configuration file.  You must add the service\n"
             "in the configuration file or remove the modules that require it.";
    }
    rndEngine_ = &(rng->getEngine(iEvent.streamID()));
  }

  // get calibration DB
  const auto& gainCalibration = iSetup.getData(gainCalibESToken_);
  const auto& thePixelTopology = iSetup.getData(thePixelTopologyToken_);

  // Step A: Get Inputs
  edm::Handle<CrossingFrame<PSimHit>> cf;

  iEvent.getByToken(tokenCrossingFramePPSPixel, cf);

  if (verbosity_) {
    edm::LogInfo("PPS") << "PixelDigiProducer \n\n=================== Starting SimHit access"
                        << "  ===================";

    MixCollection<PSimHit> col{cf.product(), std::pair(-0, 0)};
    MixCollection<PSimHit>::iterator cfi;
    int count = 0;
    for (cfi = col.begin(); cfi != col.end(); cfi++) {
      edm::LogInfo("PPS") << "PixelDigiProducer"
                          << " Hit " << count << " has tof " << cfi->timeOfFlight() << " trackid " << cfi->trackId()
                          << " bunchcr " << cfi.bunch() << " trigger " << cfi.getTrigger()
                          << ", from EncodedEventId: " << cfi->eventId().bunchCrossing() << " "
                          << cfi->eventId().event() << " bcr from MixCol " << cfi.bunch();
      edm::LogInfo("PPS") << " PixelDigiProducer Hit: " << (*cfi) << "  " << cfi->exitPoint();
      count++;
    }
  }

  MixCollection<PSimHit> allRPixHits{cf.product(), std::pair(0, 0)};

  if (verbosity_)
    edm::LogInfo("PPS") << "PixelDigiProducer Input MixCollection size = " << allRPixHits.size();

  //Loop on PSimHit
  simhit_map SimHitMap;
  SimHitMap.clear();

  MixCollection<PSimHit>::iterator isim;
  for (isim = allRPixHits.begin(); isim != allRPixHits.end(); ++isim) {
    SimHitMap[(*isim).detUnitId()].push_back((*isim));
  }

  // Step B: LOOP on hits in event
  std::vector<edm::DetSet<CTPPSPixelDigi>> theDigiVector;
  theDigiVector.reserve(400);
  theDigiVector.clear();

  for (simhit_map_iterator it = SimHitMap.begin(); it != SimHitMap.end(); ++it) {
    edm::DetSet<CTPPSPixelDigi> digi_collector(it->first);

    if (theAlgoMap.find(it->first) == theAlgoMap.end()) {
      theAlgoMap[it->first] = std::make_unique<RPixDetDigitizer>(
          conf_, *rndEngine_, it->first, thePixelTopology);  //a digitizer for any detector
    }

    std::vector<int> input_links;
    std::vector<std::vector<std::pair<int, double>>> output_digi_links;  // links to simhits

    theAlgoMap.at(it->first)->run(
        SimHitMap[it->first], input_links, digi_collector.data, output_digi_links, &gainCalibration, &thePixelTopology);

    if (!digi_collector.data.empty()) {
      theDigiVector.push_back(digi_collector);
    }
  }

  std::unique_ptr<edm::DetSetVector<CTPPSPixelDigi>> digi_output(new edm::DetSetVector<CTPPSPixelDigi>(theDigiVector));

  if (verbosity_) {
    edm::LogInfo("PPS") << "PixelDigiProducer digi_output->size()=" << digi_output->size();
  }

  iEvent.put(std::move(digi_output));
}

Member Data Documentation

edm::ParameterSet CTPPSPixelDigiProducer::conf_

private

Definition at line 93 of file PPSPixelDigiProducer.cc.

Referenced by produce().

edm::ESGetToken<CTPPSPixelGainCalibrations, CTPPSPixelGainCalibrationsRcd> CTPPSPixelDigiProducer::gainCalibESToken_

private

Definition at line 101 of file PPSPixelDigiProducer.cc.

Referenced by produce().

CLHEP::HepRandomEngine* CTPPSPixelDigiProducer::rndEngine_ = nullptr

private

Definition at line 97 of file PPSPixelDigiProducer.cc.

Referenced by produce().

std::vector<std::string> CTPPSPixelDigiProducer::RPix_hit_containers_

private

Definition at line 89 of file PPSPixelDigiProducer.cc.

Referenced by CTPPSPixelDigiProducer().

std::map<uint32_t, std::unique_ptr<RPixDetDigitizer> > CTPPSPixelDigiProducer::theAlgoMap

private

Definition at line 95 of file PPSPixelDigiProducer.cc.

Referenced by produce().

edm::ESGetToken<PPSPixelTopology, PPSPixelTopologyRcd> CTPPSPixelDigiProducer::thePixelTopologyToken_

private

Definition at line 102 of file PPSPixelDigiProducer.cc.

Referenced by produce().

edm::EDGetTokenT<CrossingFrame<PSimHit> > CTPPSPixelDigiProducer::tokenCrossingFramePPSPixel

private

Definition at line 100 of file PPSPixelDigiProducer.cc.

Referenced by CTPPSPixelDigiProducer(), and produce().

int CTPPSPixelDigiProducer::verbosity_

private

Definition at line 98 of file PPSPixelDigiProducer.cc.

Referenced by CTPPSPixelDigiProducer(), and produce().