CastorClusterProducer Class Reference

#include <RecoLocalCalo/Castor/src/CastorClusterProducer.cc>

Inheritance diagram for CastorClusterProducer:

Public Member Functions
	CastorClusterProducer (const edm::ParameterSet &)
	~CastorClusterProducer () override
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
	~EDProducer () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
std::vector< bool > const &	recordProvenanceList () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...
	~ProducerBase () noexcept(false) override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Types
typedef std::vector< reco::CastorCluster >	CastorClusterCollection
typedef std::vector< reco::CastorTower >	CastorTowerCollection
typedef ROOT::Math::RhoZPhiPoint	CellPoint
typedef math::XYZPointD	Point
typedef ROOT::Math::RhoEtaPhiPoint	TowerPoint

Private Member Functions
void	beginJob () override
void	endJob () override
double	phiangle (double testphi)
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
std::string	basicjets_
bool	clusteralgo_
std::string	input_
edm::EDGetTokenT< CastorTowerCollection >	tok_input_
edm::EDGetTokenT< reco::BasicJetCollection >	tok_jets_
edm::EDGetTokenT< CastorTowerCollection >	tok_tower_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex > >
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsInputProcessBlocks ()
static bool	wantsProcessBlocks ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
ProducesCollector	producesCollector ()
Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Description: CastorCluster Reconstruction Producer. Produces Clusters from Towers Implementation:

Definition at line 48 of file CastorClusterProducer.cc.

Member Typedef Documentation

CastorClusterCollection

typedef std::vector<reco::CastorCluster> CastorClusterProducer::CastorClusterCollection

private

Definition at line 64 of file CastorClusterProducer.cc.

CastorTowerCollection

typedef std::vector<reco::CastorTower> CastorClusterProducer::CastorTowerCollection

private

Definition at line 63 of file CastorClusterProducer.cc.

CellPoint

typedef ROOT::Math::RhoZPhiPoint CastorClusterProducer::CellPoint

private

Definition at line 62 of file CastorClusterProducer.cc.

Point

typedef math::XYZPointD CastorClusterProducer::Point

private

Definition at line 60 of file CastorClusterProducer.cc.

TowerPoint

typedef ROOT::Math::RhoEtaPhiPoint CastorClusterProducer::TowerPoint

private

Definition at line 61 of file CastorClusterProducer.cc.

Constructor & Destructor Documentation

CastorClusterProducer()

CastorClusterProducer::CastorClusterProducer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 84 of file CastorClusterProducer.cc.

    : input_(iConfig.getUntrackedParameter<std::string>("inputtowers", "")),
      basicjets_(iConfig.getUntrackedParameter<std::string>("basicjets", "")),
      clusteralgo_(iConfig.getUntrackedParameter<bool>("ClusterAlgo", false)) {
  // register for data access
  tok_input_ = consumes<CastorTowerCollection>(edm::InputTag(input_));
  tok_jets_ = consumes<reco::BasicJetCollection>(edm::InputTag(basicjets_));
  tok_tower_ = consumes<CastorTowerCollection>(edm::InputTag("CastorTowerReco"));
  // register your products
  produces<CastorClusterCollection>();
  // now do what ever other initialization is needed
}

References basicjets_, input_, tok_input_, tok_jets_, and tok_tower_.

~CastorClusterProducer()

CastorClusterProducer::~CastorClusterProducer ( )

override

Definition at line 97 of file CastorClusterProducer.cc.

                                              {
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)
}

Member Function Documentation

beginJob()

void CastorClusterProducer::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 252 of file CastorClusterProducer.cc.

{ LogDebug("CastorClusterProducer") << "Starting CastorClusterProducer"; }

References LogDebug.

endJob()

void CastorClusterProducer::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 255 of file CastorClusterProducer.cc.

{ LogDebug("CastorClusterProducer") << "Ending CastorClusterProducer"; }

References LogDebug.

phiangle()

double CastorClusterProducer::phiangle ( double  testphi )

private

Definition at line 242 of file CastorClusterProducer.cc.

                                                     {
  double phi = testphi;
  while (phi > M_PI)
    phi -= (2 * M_PI);
  while (phi < -M_PI)
    phi += (2 * M_PI);
  return phi;
}

References M_PI, and phi.

Referenced by produce().

produce()

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

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 107 of file CastorClusterProducer.cc.

                                                                                 {
  using namespace edm;
  using namespace reco;
  using namespace TMath;
 
  LogDebug("CastorClusterProducer") << "3. entering CastorClusterProducer";
 
  if (!input_.empty()) {
    // Produce CastorClusters from CastorTowers
 
    edm::Handle<CastorTowerCollection> InputTowers;
    iEvent.getByToken(tok_input_, InputTowers);
 
    auto OutputClustersfromClusterAlgo = std::make_unique<CastorClusterCollection>();
 
    // get and check input size
    int nTowers = InputTowers->size();
 
    if (nTowers == 0)
      LogDebug("CastorClusterProducer") << "Warning: You are trying to run the Cluster algorithm with 0 input towers.";
 
    CastorTowerRefVector posInputTowers, negInputTowers;
 
    for (size_t i = 0; i < InputTowers->size(); ++i) {
      reco::CastorTowerRef tower_p = reco::CastorTowerRef(InputTowers, i);
      if (tower_p->eta() > 0.)
        posInputTowers.push_back(tower_p);
      if (tower_p->eta() < 0.)
        negInputTowers.push_back(tower_p);
    }
 
    // build cluster from ClusterAlgo
    if (clusteralgo_ == true) {
      // code
      iEvent.put(std::move(OutputClustersfromClusterAlgo));
    }
  }
 
  if (!basicjets_.empty()) {
    Handle<BasicJetCollection> bjCollection;
    iEvent.getByToken(tok_jets_, bjCollection);
 
    Handle<CastorTowerCollection> ctCollection;
    iEvent.getByToken(tok_tower_, ctCollection);
 
    auto OutputClustersfromBasicJets = std::make_unique<CastorClusterCollection>();
 
    if (bjCollection->empty())
      LogDebug("CastorClusterProducer")
          << "Warning: You are trying to run the Cluster algorithm with 0 input basicjets.";
 
    for (unsigned i = 0; i < bjCollection->size(); i++) {
      const BasicJet* bj = &(*bjCollection)[i];
 
      double energy = bj->energy();
      TowerPoint temp(88.5, bj->eta(), bj->phi());
      Point position(temp);
      double emEnergy = 0.;
      double hadEnergy = 0.;
      double width = 0.;
      double depth = 0.;
      double fhot = 0.;
      double sigmaz = 0.;
      CastorTowerRefVector usedTowers;
      double zmean = 0.;
      double z2mean = 0.;
 
      std::vector<CandidatePtr> ccp = bj->getJetConstituents();
      std::vector<CandidatePtr>::const_iterator itParticle;
      for (itParticle = ccp.begin(); itParticle != ccp.end(); ++itParticle) {
        const CastorTower* castorcand = dynamic_cast<const CastorTower*>(itParticle->get());
        //cout << " castortowercandidate reference energy = " << castorcand->castorTower()->energy() << endl;
        //cout << " castortowercandidate reference eta = " << castorcand->castorTower()->eta() << endl;
        //cout << " castortowercandidate reference phi = " << castorcand->castorTower()->phi() << endl;
        //cout << " castortowercandidate reference depth = " << castorcand->castorTower()->depth() << endl;
 
        //CastorTowerCollection *ctc = new CastorTowerCollection();
        //ctc->push_back(*castorcand);
        //CastorTowerRef towerref = CastorTowerRef(ctc,0);
 
        size_t thisone = 0;
        for (size_t l = 0; l < ctCollection->size(); l++) {
          const CastorTower ct = (*ctCollection)[l];
          if (std::abs(ct.phi() - castorcand->phi()) < 0.0001) {
            thisone = l;
          }
        }
 
        CastorTowerRef towerref(ctCollection, thisone);
        usedTowers.push_back(towerref);
        emEnergy += castorcand->emEnergy();
        hadEnergy += castorcand->hadEnergy();
        depth += castorcand->depth() * castorcand->energy();
        width += pow(phiangle(castorcand->phi() - bj->phi()), 2) * castorcand->energy();
        fhot += castorcand->fhot() * castorcand->energy();
 
        // loop over rechits
        for (edm::RefVector<edm::SortedCollection<CastorRecHit> >::iterator it = castorcand->rechitsBegin();
             it != castorcand->rechitsEnd();
             it++) {
          edm::Ref<edm::SortedCollection<CastorRecHit> > rechit_p = *it;
          double Erechit = rechit_p->energy();
          HcalCastorDetId id = rechit_p->id();
          int module = id.module();
          double zrechit = 0;
          if (module < 3)
            zrechit = -14390 - 24.75 - 49.5 * (module - 1);
          if (module > 2)
            zrechit = -14390 - 99 - 49.5 - 99 * (module - 3);
          zmean += Erechit * zrechit;
          z2mean += Erechit * zrechit * zrechit;
        }  // end loop over rechits
      }
      //cout << "" << endl;
 
      depth = depth / energy;
      width = sqrt(width / energy);
      fhot = fhot / energy;
 
      zmean = zmean / energy;
      z2mean = z2mean / energy;
      double sigmaz2 = z2mean - zmean * zmean;
      if (sigmaz2 > 0)
        sigmaz = sqrt(sigmaz2);
 
      CastorCluster cc(
          energy, position, emEnergy, hadEnergy, emEnergy / energy, width, depth, fhot, sigmaz, usedTowers);
      OutputClustersfromBasicJets->push_back(cc);
    }
 
    iEvent.put(std::move(OutputClustersfromBasicJets));
  }
}

References funct::abs(), basicjets_, clusteralgo_, LEDCalibrationChannels::depth, reco::CastorTower::depth(), reco::CastorTower::emEnergy(), HCALHighEnergyHPDFilter_cfi::energy, reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), reco::CastorTower::fhot(), reco::Jet::getJetConstituents(), reco::CastorTower::hadEnergy(), mps_fire::i, edm::Ref< C, T, F >::id(), iEvent, input_, cmsLHEtoEOSManager::l, LogDebug, callgraph::module, eostools::move(), ecaldqm::nTowers, reco::LeafCandidate::phi(), phiangle(), position, funct::pow(), edm::RefVector< C, T, F >::push_back(), reco::CastorTower::rechitsBegin(), reco::CastorTower::rechitsEnd(), mathSSE::sqrt(), groupFilesInBlocks::temp, tok_input_, tok_jets_, tok_tower_, and ApeEstimator_cff::width.

Member Data Documentation

basicjets_

std::string CastorClusterProducer::basicjets_

private

Definition at line 65 of file CastorClusterProducer.cc.

Referenced by CastorClusterProducer(), and produce().

clusteralgo_

bool CastorClusterProducer::clusteralgo_

private

Definition at line 69 of file CastorClusterProducer.cc.

Referenced by produce().

input_

std::string CastorClusterProducer::input_

private

Definition at line 65 of file CastorClusterProducer.cc.

Referenced by CastorClusterProducer(), and produce().

tok_input_

edm::EDGetTokenT<CastorTowerCollection> CastorClusterProducer::tok_input_

private

Definition at line 66 of file CastorClusterProducer.cc.

Referenced by CastorClusterProducer(), and produce().

tok_jets_

edm::EDGetTokenT<reco::BasicJetCollection> CastorClusterProducer::tok_jets_

private

Definition at line 67 of file CastorClusterProducer.cc.

Referenced by CastorClusterProducer(), and produce().

tok_tower_

edm::EDGetTokenT<CastorTowerCollection> CastorClusterProducer::tok_tower_

private

Definition at line 68 of file CastorClusterProducer.cc.

Referenced by CastorClusterProducer(), and produce().