EgammaRandomSeeds.cc

Go to the documentation of this file.

#include "RecoEgamma/EgammaTools/interface/EgammaRandomSeeds.h"

#include "DataFormats/EgammaReco/interface/SuperCluster.h"

uint32_t egamma::getRandomSeedFromSC(const edm::Event& iEvent, const reco::SuperClusterRef scRef) {
  const int offset = 0;  //for future expansion
  std::seed_seq seeder = {int(iEvent.id().event()),
                          int(iEvent.id().luminosityBlock()),
                          int(iEvent.id().run()),
                          int(scRef->seed()->seed().rawId()),
                          int(scRef->seed()->hitsAndFractions().size()),
                          offset};
  uint32_t seed = 0, tries = 10;
  do {
    seeder.generate(&seed, &seed + 1);
    tries++;
  } while (seed == 0 && tries < 10);
  return seed ? seed : iEvent.id().event() + 10000 * scRef.key();
}