#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

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::vector< ModuleDescription const * > &modules, 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

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	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>
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<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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.

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

References LogDebug.

◆ endJob()

void CastorClusterProducer::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 255 of file CastorClusterProducer.cc.

255 { 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));
   }
 }