CastorClusterProducer Class Reference

#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 ()
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
	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
	EDConsumerBase ()
ProductHolderIndex	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) 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 ()
virtual void	endJob ()
double	phiangle (double testphi)
virtual void	produce (edm::Event &, const edm::EventSetup &)

Private Attributes
std::string	basicjets_
bool	clusteralgo_
std::string	input_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
typedef WorkerT< EDProducer >	WorkerType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
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::EDProducer
CurrentProcessingContext const *	currentContext () const
Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
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 51 of file CastorClusterProducer.cc.

Member Typedef Documentation

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

private

Definition at line 67 of file CastorClusterProducer.cc.

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

private

Definition at line 66 of file CastorClusterProducer.cc.

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

private

Definition at line 65 of file CastorClusterProducer.cc.

typedef math::XYZPointD CastorClusterProducer::Point

private

Definition at line 63 of file CastorClusterProducer.cc.

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

private

Definition at line 64 of file CastorClusterProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 86 of file CastorClusterProducer.cc.

                                                                           :
  input_(iConfig.getUntrackedParameter<std::string>("inputtowers","")),
  basicjets_(iConfig.getUntrackedParameter<std::string>("basicjets","")),
  clusteralgo_(iConfig.getUntrackedParameter<bool>("ClusterAlgo",false))
{
  // register your products
  produces<CastorClusterCollection>();
  // now do what ever other initialization is needed
}

CastorClusterProducer::~CastorClusterProducer

(

)

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

void CastorClusterProducer::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 246 of file CastorClusterProducer.cc.

References LogDebug.

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

void CastorClusterProducer::endJob ( void  )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 252 of file CastorClusterProducer.cc.

References LogDebug.

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

double CastorClusterProducer::phiangle ( double  testphi )

private

Definition at line 238 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  )

privatevirtual

Implements edm::EDProducer.

Definition at line 109 of file CastorClusterProducer.cc.

References abs, basicjets_, clusteralgo_, reco::CastorTower::depth(), reco::CastorTower::emEnergy(), relval_parameters_module::energy, reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), reco::CastorTower::fhot(), edm::Event::getByLabel(), reco::Jet::getJetConstituents(), reco::CastorTower::hadEnergy(), i, edm::Ref< C, T, F >::id(), input_, prof2calltree::l, LogDebug, module::module(), reco::LeafCandidate::phi(), phiangle(), position, funct::pow(), edm::RefVector< C, T, F >::push_back(), edm::Event::put(), reco::CastorTower::rechitsBegin(), reco::CastorTower::rechitsEnd(), L1Trigger_dataformats::reco, mathSSE::sqrt(), groupFilesInBlocks::temp, 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.getByLabel(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.getByLabel(basicjets_,bjCollection);
    
    Handle<CastorTowerCollection> ctCollection;
    iEvent.getByLabel("CastorTowerReco",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);
  
  }
 
}

Member Data Documentation

std::string CastorClusterProducer::basicjets_

private

Definition at line 68 of file CastorClusterProducer.cc.

Referenced by produce().

bool CastorClusterProducer::clusteralgo_

private

Definition at line 69 of file CastorClusterProducer.cc.

Referenced by produce().

std::string CastorClusterProducer::input_

private

Definition at line 68 of file CastorClusterProducer.cc.

Referenced by produce().