HiSuperClusterProducer Class Reference

Inheritance diagram for HiSuperClusterProducer:

Public Member Functions
virtual void	endJob ()
	HiSuperClusterProducer (const edm::ParameterSet &ps)
void	produce (edm::Event &, const edm::EventSetup &) 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

Private Member Functions
void	getClusterPtrVector (edm::Event &evt, const edm::EDGetTokenT< reco::BasicClusterCollection > &clustersToken, reco::CaloClusterPtrVector *)
void	outputValidationInfo (reco::SuperClusterCollection &superclusterCollection)
void	produceSuperclustersForECALPart (edm::Event &evt, const edm::EDGetTokenT< reco::BasicClusterCollection > &clustersToken, edm::EDPutTokenT< reco::SuperClusterCollection > const &putToken)

Private Attributes
const float	barrelBCEnergyThreshold_
const float	barrelEtaSearchRoad_
const float	barrelPhiSearchRoad_
const edm::EDPutTokenT< reco::SuperClusterCollection >	barrelSuperclusterPut_
std::unique_ptr< HiBremRecoveryClusterAlgo >	bremAlgo_p
const bool	doBarrel_
const bool	doEndcaps_
const edm::EDGetTokenT< reco::BasicClusterCollection >	ebClustersToken_
const edm::EDGetTokenT< reco::BasicClusterCollection >	eeClustersToken_
const float	endcapBCEnergyThreshold_
const float	endcapEtaSearchRoad_
const float	endcapPhiSearchRoad_
const edm::EDPutTokenT< reco::SuperClusterCollection >	endcapSuperclusterPut_
int	noSuperClusters = 0
const float	seedTransverseEnergyThreshold_
double	totalE = 0.0

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 20 of file HiSuperClusterProducer.cc.

Constructor & Destructor Documentation

HiSuperClusterProducer()

HiSuperClusterProducer::HiSuperClusterProducer

(

const edm::ParameterSet &

ps

)

Definition at line 65 of file HiSuperClusterProducer.cc.

References edm::ParameterSet::getParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

    : endcapSuperclusterPut_{produces<reco::SuperClusterCollection>(
          ps.getParameter<std::string>("endcapSuperclusterCollection"))},
      barrelSuperclusterPut_{
          produces<reco::SuperClusterCollection>(ps.getParameter<std::string>("barrelSuperclusterCollection"))},
      eeClustersToken_{consumes<reco::BasicClusterCollection>(
          edm::InputTag(ps.getParameter<std::string>("endcapClusterProducer"),
                        ps.getParameter<std::string>("endcapClusterCollection")))},
      ebClustersToken_{consumes<reco::BasicClusterCollection>(
          edm::InputTag(ps.getParameter<std::string>("barrelClusterProducer"),
                        ps.getParameter<std::string>("barrelClusterCollection")))},
      barrelEtaSearchRoad_{float(ps.getParameter<double>("barrelEtaSearchRoad"))},
      barrelPhiSearchRoad_{float(ps.getParameter<double>("barrelPhiSearchRoad"))},
      endcapEtaSearchRoad_{float(ps.getParameter<double>("endcapEtaSearchRoad"))},
      endcapPhiSearchRoad_{float(ps.getParameter<double>("endcapPhiSearchRoad"))},
      seedTransverseEnergyThreshold_{float(ps.getParameter<double>("seedTransverseEnergyThreshold"))},
      barrelBCEnergyThreshold_{float(ps.getParameter<double>("barrelBCEnergyThreshold"))},
      endcapBCEnergyThreshold_{float(ps.getParameter<double>("endcapBCEnergyThreshold"))},
      doBarrel_{ps.getParameter<bool>("doBarrel")},
      doEndcaps_{ps.getParameter<bool>("doEndcaps")}

{
  // The verbosity level
  HiBremRecoveryClusterAlgo::VerbosityLevel verbosity;
  std::string verbosityString = ps.getParameter<std::string>("VerbosityLevel");
  if (verbosityString == "DEBUG")
    verbosity = HiBremRecoveryClusterAlgo::pDEBUG;
  else if (verbosityString == "WARNING")
    verbosity = HiBremRecoveryClusterAlgo::pWARNING;
  else if (verbosityString == "INFO")
    verbosity = HiBremRecoveryClusterAlgo::pINFO;
  else
    verbosity = HiBremRecoveryClusterAlgo::pERROR;

  if (verbosityString == "INFO") {
    std::cout << "Barrel BC Energy threshold = " << barrelBCEnergyThreshold_ << std::endl;
    std::cout << "Endcap BC Energy threshold = " << endcapBCEnergyThreshold_ << std::endl;
  }

  bremAlgo_p = std::make_unique<HiBremRecoveryClusterAlgo>(barrelEtaSearchRoad_,
                                                           barrelPhiSearchRoad_,
                                                           endcapEtaSearchRoad_,
                                                           endcapPhiSearchRoad_,
                                                           seedTransverseEnergyThreshold_,
                                                           barrelBCEnergyThreshold_,
                                                           endcapBCEnergyThreshold_,
                                                           verbosity);
}

Member Function Documentation

endJob()

void HiSuperClusterProducer::endJob ( void  )

virtual

Definition at line 114 of file HiSuperClusterProducer.cc.

References noSuperClusters, str, and totalE.

                                    {
  double averEnergy = 0.;
  std::ostringstream str;
  str << "HiSuperClusterProducer::endJob()\n"
      << "  total # reconstructed super clusters: " << noSuperClusters << "\n"
      << "  total energy of all clusters: " << totalE << "\n";
  if (noSuperClusters > 0) {
    averEnergy = totalE / noSuperClusters;
    str << "  average SuperCluster energy = " << averEnergy << "\n";
  }
  edm::LogInfo("HiSuperClusterProducerInfo") << str.str() << "\n";
}

getClusterPtrVector()

void HiSuperClusterProducer::getClusterPtrVector ( edm::Event &  evt, const edm::EDGetTokenT< reco::BasicClusterCollection > &  clustersToken, reco::CaloClusterPtrVector *  clusterPtrVector_p  )

private

Definition at line 156 of file HiSuperClusterProducer.cc.

References edm::Event::getHandle(), mps_fire::i, edm::HandleBase::isValid(), edm::Handle< T >::product(), and edm::PtrVector< T >::push_back().

Referenced by produceSuperclustersForECALPart().

                                                                                               {
  edm::Handle<reco::BasicClusterCollection> bccHandle = evt.getHandle(clustersToken);

  if (!(bccHandle.isValid())) {
    edm::LogError("HiSuperClusterProducerError") << "could not get a handle on the BasicCluster Collection!";
    clusterPtrVector_p = nullptr;
  }

  const reco::BasicClusterCollection* clusterCollection_p = bccHandle.product();
  for (unsigned int i = 0; i < clusterCollection_p->size(); i++) {
    clusterPtrVector_p->push_back(reco::CaloClusterPtr(bccHandle, i));
  }
}

outputValidationInfo()

void HiSuperClusterProducer::outputValidationInfo ( reco::SuperClusterCollection &  superclusterCollection )

private

produce()

void HiSuperClusterProducer::produce ( edm::Event &  evt, const edm::EventSetup &  es  )

override

Definition at line 127 of file HiSuperClusterProducer.cc.

References barrelSuperclusterPut_, doBarrel_, doEndcaps_, ebClustersToken_, eeClustersToken_, endcapSuperclusterPut_, and produceSuperclustersForECALPart().

                                                                           {
  if (doEndcaps_)
    produceSuperclustersForECALPart(evt, eeClustersToken_, endcapSuperclusterPut_);

  if (doBarrel_)
    produceSuperclustersForECALPart(evt, ebClustersToken_, barrelSuperclusterPut_);
}

produceSuperclustersForECALPart()

void HiSuperClusterProducer::produceSuperclustersForECALPart ( edm::Event &  evt, const edm::EDGetTokenT< reco::BasicClusterCollection > &  clustersToken, edm::EDPutTokenT< reco::SuperClusterCollection > const &  putToken  )

private

Definition at line 135 of file HiSuperClusterProducer.cc.

References bremAlgo_p, edm::Event::emplace(), getClusterPtrVector(), eostools::move(), noSuperClusters, and totalE.

Referenced by produce().

                                                                {
  // get the cluster collection out and turn it to a BasicClusterRefVector:
  reco::CaloClusterPtrVector clusterPtrVector_p{};
  getClusterPtrVector(evt, clustersToken, &clusterPtrVector_p);

  // run the brem recovery and get the SC collection
  reco::SuperClusterCollection superclusters_ap = bremAlgo_p->makeSuperClusters(clusterPtrVector_p);

  // count the total energy and the number of superclusters
  for (auto const& cluster : superclusters_ap) {
    totalE += cluster.energy();
    noSuperClusters++;
  }

  // put the SC collection in the event
  evt.emplace(putToken, std::move(superclusters_ap));
}

Member Data Documentation

barrelBCEnergyThreshold_

const float HiSuperClusterProducer::barrelBCEnergyThreshold_

private

Definition at line 39 of file HiSuperClusterProducer.cc.

barrelEtaSearchRoad_

const float HiSuperClusterProducer::barrelEtaSearchRoad_

private

Definition at line 34 of file HiSuperClusterProducer.cc.

barrelPhiSearchRoad_

const float HiSuperClusterProducer::barrelPhiSearchRoad_

private

Definition at line 35 of file HiSuperClusterProducer.cc.

barrelSuperclusterPut_

const edm::EDPutTokenT<reco::SuperClusterCollection> HiSuperClusterProducer::barrelSuperclusterPut_

private

Definition at line 29 of file HiSuperClusterProducer.cc.

Referenced by produce().

bremAlgo_p

std::unique_ptr<HiBremRecoveryClusterAlgo> HiSuperClusterProducer::bremAlgo_p

private

Definition at line 45 of file HiSuperClusterProducer.cc.

Referenced by produceSuperclustersForECALPart().

doBarrel_

const bool HiSuperClusterProducer::doBarrel_

private

Definition at line 42 of file HiSuperClusterProducer.cc.

Referenced by produce().

doEndcaps_

const bool HiSuperClusterProducer::doEndcaps_

private

Definition at line 43 of file HiSuperClusterProducer.cc.

Referenced by produce().

ebClustersToken_

const edm::EDGetTokenT<reco::BasicClusterCollection> HiSuperClusterProducer::ebClustersToken_

private

Definition at line 32 of file HiSuperClusterProducer.cc.

Referenced by produce().

eeClustersToken_

const edm::EDGetTokenT<reco::BasicClusterCollection> HiSuperClusterProducer::eeClustersToken_

private

Definition at line 31 of file HiSuperClusterProducer.cc.

Referenced by produce().

endcapBCEnergyThreshold_

const float HiSuperClusterProducer::endcapBCEnergyThreshold_

private

Definition at line 40 of file HiSuperClusterProducer.cc.

endcapEtaSearchRoad_

const float HiSuperClusterProducer::endcapEtaSearchRoad_

private

Definition at line 36 of file HiSuperClusterProducer.cc.

endcapPhiSearchRoad_

const float HiSuperClusterProducer::endcapPhiSearchRoad_

private

Definition at line 37 of file HiSuperClusterProducer.cc.

endcapSuperclusterPut_

const edm::EDPutTokenT<reco::SuperClusterCollection> HiSuperClusterProducer::endcapSuperclusterPut_

private

Definition at line 28 of file HiSuperClusterProducer.cc.

Referenced by produce().

noSuperClusters

int HiSuperClusterProducer::noSuperClusters = 0

private

Definition at line 48 of file HiSuperClusterProducer.cc.

Referenced by endJob(), and produceSuperclustersForECALPart().

seedTransverseEnergyThreshold_

const float HiSuperClusterProducer::seedTransverseEnergyThreshold_

private

Definition at line 38 of file HiSuperClusterProducer.cc.

totalE

double HiSuperClusterProducer::totalE = 0.0

private

Definition at line 47 of file HiSuperClusterProducer.cc.

Referenced by endJob(), and produceSuperclustersForECALPart().