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

Inheritance diagram for CastorClusterProducer:

Public Member Functions
	CastorClusterProducer (const edm::ParameterSet &)

	~CastorClusterProducer ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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
virtual void	beginJob () override

virtual void	endJob () override

double	phiangle (double testphi)

virtual 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
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)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, 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 ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

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

Definition at line 50 of file CastorClusterProducer.cc.

Member Typedef Documentation

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

private

Definition at line 66 of file CastorClusterProducer.cc.

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

private

Definition at line 65 of file CastorClusterProducer.cc.

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

private

Definition at line 64 of file CastorClusterProducer.cc.

typedef math::XYZPointD CastorClusterProducer::Point

private

Definition at line 62 of file CastorClusterProducer.cc.

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

private

Definition at line 63 of file CastorClusterProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 87 of file CastorClusterProducer.cc.

References basicjets_, input_, HLT_25ns14e33_v1_cff::InputTag, tok_input_, tok_jets_, and tok_tower_.

                                                                            :
   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
 }

CastorClusterProducer::~CastorClusterProducer ( )

Definition at line 102 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

void CastorClusterProducer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 251 of file CastorClusterProducer.cc.

References LogDebug.

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

void CastorClusterProducer::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 257 of file CastorClusterProducer.cc.

References LogDebug.

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

double CastorClusterProducer::phiangle ( double testphi )

private

Definition at line 243 of file CastorClusterProducer.cc.

References M_PI, and phi.

Referenced by produce().

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

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

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 114 of file CastorClusterProducer.cc.

References funct::abs(), basicjets_, clusteralgo_, reco::CastorTower::depth(), HLT_25ns14e33_v1_cff::depth, reco::CastorTower::emEnergy(), relval_parameters_module::energy, reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), reco::CastorTower::fhot(), edm::Event::getByToken(), reco::Jet::getJetConstituents(), reco::CastorTower::hadEnergy(), i, edm::Ref< C, T, F >::id(), input_, prof2calltree::l, LogDebug, module::module(), ecaldqm::nTowers, reco::LeafCandidate::phi(), phiangle(), position, funct::pow(), edm::RefVector< C, T, F >::push_back(), edm::Event::put(), reco::CastorTower::rechitsBegin(), reco::CastorTower::rechitsEnd(), dt_dqm_sourceclient_common_cff::reco, mathSSE::sqrt(), groupFilesInBlocks::temp, tok_input_, tok_jets_, tok_tower_, and create_public_lumi_plots::width.

                                                                                  {
   
   using namespace edm;
   using namespace reco;
   using namespace TMath;
   
   LogDebug("CastorClusterProducer")
     <<"3. entering CastorClusterProducer";
   
   if ( input_ != "") {
   
   // Produce CastorClusters from CastorTowers
   
   edm::Handle<CastorTowerCollection> InputTowers;
   iEvent.getByToken(tok_input_, InputTowers);
 
   std::auto_ptr<CastorClusterCollection> OutputClustersfromClusterAlgo (new 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(OutputClustersfromClusterAlgo);
   }
   
   }
   
   if ( basicjets_ != "") {
   
     Handle<BasicJetCollection> bjCollection;
     iEvent.getByToken(tok_jets_,bjCollection);
     
     Handle<CastorTowerCollection> ctCollection;
     iEvent.getByToken(tok_tower_,ctCollection);
     
     std::auto_ptr<CastorClusterCollection> OutputClustersfromBasicJets (new CastorClusterCollection);
     
     if (bjCollection->size()==0) 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(OutputClustersfromBasicJets);
   
   }
  
 }