---

SubSeedGenerator Class Reference

#include <RecoEgamma/EgammaElectronAlgos/interface/SubSeedGenerator.h>

List of all members.

Public Member Functions
void	run (edm::Event &, const edm::EventSetup &setup, const edm::Handle< reco::SuperClusterCollection > &, reco::ElectronPixelSeedCollection &)
void	setupES (const edm::EventSetup &setup)
	SubSeedGenerator (const edm::ParameterSet &conf)
	~SubSeedGenerator ()
Private Attributes
double	deta_
double	dphi_
double	dr_
edm::InputTag	initialSeeds_
double	pt_

---

Detailed Description

Definition at line 19 of file SubSeedGenerator.h.

---

Constructor & Destructor Documentation

SubSeedGenerator::SubSeedGenerator

(

const edm::ParameterSet &

conf

)

Definition at line 25 of file SubSeedGenerator.cc.

References deta_, dphi_, dr_, edm::ParameterSet::getParameter(), initialSeeds_, and pt_.

                                                              {
  initialSeeds_               = conf.getParameter<edm::InputTag>("initialSeeds");
  dr_                  = conf.getParameter<double>("seedDr");
  dphi_                = conf.getParameter<double>("seedDPhi");
  deta_                = conf.getParameter<double>("seedDEta");
  pt_                  = conf.getParameter<double>("seedPt");
}

SubSeedGenerator::~SubSeedGenerator

(

)

Definition at line 33 of file SubSeedGenerator.cc.

                                    {
  
}

---

Member Function Documentation

void SubSeedGenerator::run	(	edm::Event &	e,
		const edm::EventSetup &	setup,
		const edm::Handle< reco::SuperClusterCollection > &	superClusters,
		reco::ElectronPixelSeedCollection &	out
	)

Definition at line 38 of file SubSeedGenerator.cc.

References funct::cos(), deta_, dPhi(), dphi_, reco::CaloCluster::energy(), edm::EventSetup::get(), edm::Event::getByLabel(), i, initialSeeds_, PV3DBase< T, PVType, FrameType >::perp(), reco::CaloCluster::position(), pt_, funct::sin(), PV3DBase< T, PVType, FrameType >::x(), reco::CaloCluster::x(), PV3DBase< T, PVType, FrameType >::y(), reco::CaloCluster::y(), PV3DBase< T, PVType, FrameType >::z(), and reco::CaloCluster::z().

Referenced by GlobalSeedProducer::produce().

                                                                                                                                                                     {
  
  edm::ESHandle<TrackerGeometry> tracker;
  setup.get<TrackerDigiGeometryRecord>().get(tracker);
  
  // get initial TrajectorySeeds
  edm::Handle<TrajectorySeedCollection> theInitialSeedColl;
  e.getByLabel(initialSeeds_, theInitialSeedColl);
  
  //seeds selection
  for(unsigned int i=0; i< superClusters->size(); ++i) {
    reco::SuperCluster theClus = (*superClusters)[i];   
    
    std::vector<TrajectorySeed>::const_iterator seed_iter;
    for(seed_iter = theInitialSeedColl->begin(); seed_iter != theInitialSeedColl->end(); ++seed_iter) {
      
      GlobalPoint  gp = tracker->idToDet( DetId(seed_iter->startingState().detId()))->surface().toGlobal( seed_iter->startingState().parameters().position());
      GlobalVector gv = tracker->idToDet( DetId(seed_iter->startingState().detId()))->surface().toGlobal( seed_iter->startingState().parameters().momentum());
      
      math::XYZVector seedGlobalDir(gv.x(),gv.y(),gv.z());   
      math::XYZVector clusterGlobalDir(theClus.x() - gp.x(), theClus.y() - gp.y(), theClus.z() - gp.z());
      
      double clusEt = theClus.energy()*sin(clusterGlobalDir.theta());
      double clusEstimatedCurvature = clusEt/0.3/4*100;  //4 tesla (temporary solution)
      double DphiBending = theClus.position().rho()/2./clusEstimatedCurvature; //ecal radius
      
      //cout << "=== et,curvature, phiBending: " 
      //   << clusEt << " , "
      //   << clusEstimatedCurvature << " , "
      //   << DphiBending << endl;
      

      double tmpDr = ROOT::Math::VectorUtil::DeltaR(clusterGlobalDir, seedGlobalDir);
      float dEta = fabs(clusterGlobalDir.Eta() - seedGlobalDir.Eta());
      float dPhi = fabs(acos(cos(clusterGlobalDir.Phi() - seedGlobalDir.Phi())));     

      float dPhi1 = fabs(dPhi - DphiBending);
      float dPhi2 = fabs(dPhi + DphiBending);
      
      if (dEta <= deta_) {
        if (dPhi1 <= dphi_|| dPhi2 <= dphi_ ) {
          if (gv.perp() > pt_) {
            //if(tmpDr <= dr_) {  
              edm::Ref<reco::SuperClusterCollection> sclRef=edm::Ref<reco::SuperClusterCollection> (superClusters,i);
              out.push_back(reco::ElectronPixelSeed(sclRef,*seed_iter)); 
              //}
          }
        }
      }
    }
    
  }//end loop over cluster
  
  
  edm::LogVerbatim("myElectronProd") << "========== SubSeedsCollectionProducer Info ==========";
  edm::LogVerbatim("myElectronProd") << "number of initial seeds: " << theInitialSeedColl->size();
  edm::LogVerbatim("myElectronProd") << "number of filtered seeds: " << out.size();
  edm::LogVerbatim("myElectronProd") << "=================================================";
  
}

void SubSeedGenerator::setupES

(

const edm::EventSetup &

setup

)

[inline]

Definition at line 25 of file SubSeedGenerator.h.

{;} //FIXME: temporary

---

Member Data Documentation

double SubSeedGenerator::deta_ [private]

Definition at line 31 of file SubSeedGenerator.h.

Referenced by run(), and SubSeedGenerator().

double SubSeedGenerator::dphi_ [private]

Definition at line 31 of file SubSeedGenerator.h.

Referenced by run(), and SubSeedGenerator().

double SubSeedGenerator::dr_ [private]

Definition at line 31 of file SubSeedGenerator.h.

Referenced by SubSeedGenerator().

edm::InputTag SubSeedGenerator::initialSeeds_ [private]

Definition at line 29 of file SubSeedGenerator.h.

Referenced by run(), and SubSeedGenerator().

double SubSeedGenerator::pt_ [private]

Definition at line 31 of file SubSeedGenerator.h.

Referenced by run(), and SubSeedGenerator().

---

The documentation for this class was generated from the following files:

• RecoEgamma/EgammaElectronAlgos/interface/SubSeedGenerator.h
• RecoEgamma/EgammaElectronAlgos/src/SubSeedGenerator.cc

---