HDetIdAssociator Class Reference

#include <DetIdAssociator.h>

Inheritance diagram for HDetIdAssociator:

Public Member Functions
virtual std::set< DetId >	getCrossedDetIds (const std::set< DetId > &, const std::vector< GlobalPoint > &trajectory)
virtual std::set< DetId >	getDetIdsCloseToAPoint (const GlobalPoint &, const int idR=0)
virtual std::set< DetId >	getDetIdsCloseToAPoint (const GlobalPoint &point, const double dR=0)
virtual std::set< DetId >	getDetIdsInACone (const std::set< DetId > &, const std::vector< GlobalPoint > &trajectory, const double)
virtual std::set< DetId >	getMaxEDetId (const std::set< DetId > &, edm::Handle< CaloTowerCollection > caloTowers)
virtual std::set< DetId >	getMaxEDetId (const std::set< DetId > &, edm::Handle< HBHERecHitCollection > recHits)
virtual std::vector< GlobalPoint >	getTrajectory (const FreeTrajectoryState &, const std::vector< GlobalPoint > &)
	HDetIdAssociator ()
	HDetIdAssociator (const int nPhi, const int nEta, const double etaBinSize)
virtual int	iEta (const GlobalPoint &)
virtual int	iPhi (const GlobalPoint &)
virtual void	setPropagator (Propagator *ptr)
virtual	~HDetIdAssociator ()

Protected Member Functions
virtual void	buildMap ()
virtual void	check_setup ()
virtual std::set< DetId >	getASetOfValidDetIds ()=0
virtual std::vector< GlobalPoint >	getDetIdPoints (const DetId &)=0
virtual GlobalPoint	getPosition (const DetId &)=0
virtual bool	insideElement (const GlobalPoint &, const DetId &)=0
virtual bool	nearElement (const GlobalPoint &point, const DetId &id, const double distance)

Protected Attributes
const double	etaBinSize_
Propagator *	ivProp_
const int	nEta_
const int	nPhi_
std::vector< std::vector< std::set< DetId > > > *	theMap_

Detailed Description

Definition at line 33 of file DetIdAssociator.h.

Constructor & Destructor Documentation

HDetIdAssociator::HDetIdAssociator ( )

inline

Definition at line 35 of file DetIdAssociator.h.

:theMap_(nullptr),nPhi_(0),nEta_(0),etaBinSize_(0),ivProp_(nullptr){};

HDetIdAssociator::HDetIdAssociator ( const int  nPhi, const int  nEta, const double  etaBinSize  )

inline

Definition at line 36 of file DetIdAssociator.h.

     :theMap_(nullptr),nPhi_(nPhi),nEta_(nEta),etaBinSize_(etaBinSize),ivProp_(nullptr){};

virtual HDetIdAssociator::~HDetIdAssociator ( )

inlinevirtual

Definition at line 39 of file DetIdAssociator.h.

References getDetIdsCloseToAPoint(), and getTrajectory().

{};

Member Function Documentation

void HDetIdAssociator::buildMap ( )

protectedvirtual

Definition at line 249 of file DetIdAssociator.cc.

References check_setup(), PV3DBase< T, PVType, FrameType >::eta(), etaBinSize_, ALCARECOTkAlBeamHalo_cff::etaMax, ALCARECOTkAlBeamHalo_cff::etaMin, getASetOfValidDetIds(), getPosition(), iEta(), iPhi(), LogTrace, nEta_, nPhi_, AlignmentTrackSelector_cfi::phiMax, AlignmentTrackSelector_cfi::phiMin, point, std::swap(), and theMap_.

Referenced by check_setup(), and getDetIdsCloseToAPoint().

{
// modified version: take only detector central position
   check_setup();
   LogTrace("HDetIdAssociator")<<"building map" << "\n";
   if(theMap_) delete theMap_;
   theMap_ = new std::vector<std::vector<std::set<DetId> > >(nEta_,std::vector<std::set<DetId> >(nPhi_));
   int numberOfDetIdsOutsideEtaRange = 0;
   int numberOfDetIdsActive = 0;
   std::set<DetId> validIds = getASetOfValidDetIds();
   for (std::set<DetId>::const_iterator id_itr = validIds.begin(); id_itr!=validIds.end(); id_itr++) {   
//      std::vector<GlobalPoint> points = getDetIdPoints(*id_itr);
      GlobalPoint point = getPosition(*id_itr);
// reject fake DetIds (eta=0 - what are they anyway???)
      if(point.eta()==0)continue;

      int ieta = iEta(point);
      int iphi = iPhi(point);
      int etaMax(-1);
      int etaMin(nEta_);
      int phiMax(-1);
      int phiMin(nPhi_);
      if ( iphi >= nPhi_ ) iphi = iphi % nPhi_;
      assert (iphi>=0);
      if ( etaMin > ieta) etaMin = ieta;
      if ( etaMax < ieta) etaMax = ieta;
      if ( phiMin > iphi) phiMin = iphi;
      if ( phiMax < iphi) phiMax = iphi;
// for abs(eta)>1.8 one tower covers two phi segments
      if ((ieta>54||ieta<15) && iphi%2==0) phiMax++;
      if ((ieta>54||ieta<15) && iphi%2==1) phiMin--;

      if (etaMax<0||phiMax<0||etaMin>=nEta_||phiMin>=nPhi_) {
     LogTrace("HDetIdAssociator")<<"Out of range: DetId:" << id_itr->rawId() <<
       "\n\teta (min,max): " << etaMin << "," << etaMax <<
       "\n\tphi (min,max): " << phiMin << "," << phiMax <<
       "\nTower id: " << id_itr->rawId() << "\n";
     numberOfDetIdsOutsideEtaRange++;
     continue;
      }
      
      if (phiMax-phiMin > phiMin+nPhi_-phiMax){
     phiMin += nPhi_;
     std::swap(phiMin,phiMax);
      }
      for(int ieta = etaMin; ieta <= etaMax; ieta++)
    for(int iphi = phiMin; iphi <= phiMax; iphi++)
      (*theMap_)[ieta][iphi%nPhi_].insert(*id_itr);
      numberOfDetIdsActive++;
   }
   LogTrace("HDetIdAssociator") << "Number of elements outside the allowed range ( |eta|>"<<
     nEta_/2*etaBinSize_ << "): " << numberOfDetIdsOutsideEtaRange << "\n";
   LogTrace("HDetIdAssociator") << "Number of active DetId's mapped: " << 
     numberOfDetIdsActive << "\n";
}

virtual void HDetIdAssociator::check_setup ( )

inlineprotectedvirtual

Reimplemented in HCaloDetIdAssociator.

Definition at line 73 of file DetIdAssociator.h.

References buildMap(), etaBinSize_, Exception, getASetOfValidDetIds(), getDetIdPoints(), getPosition(), insideElement(), ivProp_, nEta_, and nPhi_.

Referenced by buildMap(), HCaloDetIdAssociator::check_setup(), getCrossedDetIds(), getDetIdsCloseToAPoint(), getDetIdsInACone(), getMaxEDetId(), and getTrajectory().

     {
    if (nEta_==0) throw cms::Exception("FatalError") << "Number of eta bins is not set.\n";
    if (nPhi_==0) throw cms::Exception("FatalError") << "Number of phi bins is not set.\n";
    if (ivProp_==nullptr) throw cms::Exception("FatalError") << "Track propagator is not defined\n";
    if (etaBinSize_==0) throw cms::Exception("FatalError") << "Eta bin size is not set.\n";
     }

virtual std::set<DetId> HDetIdAssociator::getASetOfValidDetIds ( )

protectedpure virtual

Implemented in HCaloDetIdAssociator, HHcalDetIdAssociator, and HEcalDetIdAssociator.

Referenced by buildMap(), and check_setup().

std::set< DetId > HDetIdAssociator::getCrossedDetIds ( const std::set< DetId > &  inset, const std::vector< GlobalPoint > &  trajectory  )

virtual

Definition at line 351 of file DetIdAssociator.cc.

References check_setup(), and insideElement().

Referenced by getDetIdsCloseToAPoint().

{
   check_setup();
   std::set<DetId> outset;
   for(std::set<DetId>::const_iterator id_iter = inset.begin(); id_iter != inset.end(); id_iter++)
     for(std::vector<GlobalPoint>::const_iterator point_iter = trajectory.begin(); point_iter != trajectory.end(); point_iter++)
       if (insideElement(*point_iter, *id_iter))  outset.insert(*id_iter);
   return outset;
}

virtual std::vector<GlobalPoint> HDetIdAssociator::getDetIdPoints ( const DetId &  )

protectedpure virtual

Implemented in HCaloDetIdAssociator.

Referenced by check_setup().

std::set< DetId > HDetIdAssociator::getDetIdsCloseToAPoint ( const GlobalPoint &  direction, const int  idR = 0  )

virtual

Definition at line 143 of file DetIdAssociator.cc.

References begin, buildMap(), check_setup(), end, PV3DBase< T, PVType, FrameType >::eta(), mps_fire::i, iEta(), iPhi(), LogTrace, nEta_, nPhi_, PV3DBase< T, PVType, FrameType >::phi(), and theMap_.

Referenced by getDetIdsCloseToAPoint(), and ~HDetIdAssociator().

{
   std::set<DetId> set;
   check_setup();
   int nDets=0;
   if (! theMap_) buildMap();
   LogTrace("MatchPoint") << "point (eta,phi): " << direction.eta() << "," << direction.phi() << "\n";
   int ieta = iEta(direction);
   int iphi = iPhi(direction);

   LogTrace("MatchPoint") << "(ieta,iphi): " << ieta << "," << iphi << "\n";
   
   if (ieta>=0 && ieta<nEta_ && iphi>=0 && iphi<nPhi_){
      set = (*theMap_)[ieta][iphi];
      nDets++;
      if (idR>0){
      LogTrace("MatchPoint") << "Add neighbors (ieta,iphi): " << ieta << "," << iphi << "\n";
     //add neighbors
     int maxIEta = ieta+idR;
     int minIEta = ieta-idR;
     if(maxIEta>=nEta_) maxIEta = nEta_-1;
     if(minIEta<0) minIEta = 0;
     int maxIPhi = iphi+idR;
     int minIPhi = iphi-idR;
     if(minIPhi<0) {
        minIPhi+=nPhi_;
        maxIPhi+=nPhi_;
     }
     LogTrace("MatchPoint") << "\tieta (min,max): " << minIEta << "," << maxIEta<< "\n";
     LogTrace("MatchPoint") << "\tiphi (min,max): " << minIPhi << "," << maxIPhi<< "\n";
     for (int i=minIEta;i<=maxIEta;i++)
       for (int j=minIPhi;j<=maxIPhi;j++) {
          if( i==ieta && j==iphi) continue; // already in the set
          set.insert((*theMap_)[i][j%nPhi_].begin(),(*theMap_)[i][j%nPhi_].end());
              nDets++;
       }
      }
   }
//     if(set.size() > 0) {
//   if (ieta+idR<55 && ieta-idR>14 && set.size() != (2*idR+1)*(2*idR+1)){
//       std::cout<<" RRRA: "<<set.size()<<" DetIds in region "<<ieta<<" "<<iphi<<std::endl;
//       for( std::set<DetId>::const_iterator itr=set.begin(); itr!=set.end(); itr++) {
//         GlobalPoint point = getPosition(*itr);
//         std::cout << "DetId: " <<itr->rawId() <<" (eta,phi): " << point.eta() << "," << point.phi()<<" "<<iEta(point)<<" "<<iPhi(point)<<std::endl;
//       }
//     }
//     else {
//       std::cout <<" HDetIdAssociator::getDetIdsCloseToAPoint::There are strange days "<<std::endl; 
//     }  
   return set;
}

virtual std::set<DetId> HDetIdAssociator::getDetIdsCloseToAPoint ( const GlobalPoint &  point, const double  dR = 0  )

inlinevirtual

Definition at line 47 of file DetIdAssociator.h.

References eleHcalExtractorBlocks_cff::caloTowers, PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, etaBinSize_, getCrossedDetIds(), getDetIdsCloseToAPoint(), getDetIdsInACone(), getMaxEDetId(), iEta(), createfilelist::int, iPhi(), and nPhi_.

     {
    int etaIdR = int(dR/etaBinSize_); 
    int phiIdR = int(dR/(2*3.1416)*nPhi_);
    if (etaIdR>phiIdR)
      return getDetIdsCloseToAPoint(point, 1+etaIdR);
    else
      return getDetIdsCloseToAPoint(point, 1+phiIdR);
     }

std::set< DetId > HDetIdAssociator::getDetIdsInACone ( const std::set< DetId > &  inset, const std::vector< GlobalPoint > &  trajectory, const double  dR  )

virtual

Definition at line 306 of file DetIdAssociator.cc.

References begin, check_setup(), relativeConstraints::empty, end, getPosition(), mps_fire::i, iEta(), iPhi(), nearElement(), nEta_, nPhi_, point, and theMap_.

Referenced by getDetIdsCloseToAPoint().

{
// modified version: if dR<0, returns 3x3 towers around the input one (Michal)
   check_setup();
   std::set<DetId> outset;
  
   if(dR>=0) {
     for(std::set<DetId>::const_iterator id_iter = inset.begin(); id_iter != inset.end(); id_iter++)
       for(std::vector<GlobalPoint>::const_iterator point_iter = trajectory.begin(); point_iter != trajectory.end(); point_iter++)
         if (nearElement(*point_iter,*id_iter,dR)) outset.insert(*id_iter);
   }
   else {
     if (inset.size()!=1) return outset;
     std::set<DetId>::const_iterator id_inp = inset.begin();
     int ieta;
     int iphi;
     GlobalPoint point = getPosition(*id_inp);
     ieta = iEta(point);
     iphi = iPhi(point);
     for (int i=ieta-1;i<=ieta+1;i++) {
       for (int j=iphi-1;j<=iphi+1;j++) {
//         if( i==ieta && j==iphi) continue;
         if( i<0 || i>=nEta_) continue;
         int j2fill = j%nPhi_;
         if(j2fill<0) j2fill+=nPhi_;
         if((*theMap_)[i][j2fill].empty())continue;
         outset.insert((*theMap_)[i][j2fill].begin(),(*theMap_)[i][j2fill].end());
       }
     }
   }

//   if (outset.size() > 0) {
//     std::cout<<" RRRA: DetIds in cone:"<<std::endl;
//     for( std::set<DetId>::const_iterator itr=outset.begin(); itr!=outset.end(); itr++) {
//       GlobalPoint point = getPosition(*itr);
//       std::cout << "DetId: " <<itr->rawId() <<" (eta,phi): " << point.eta() << "," << point.phi()<<std::endl;
//     }
//   }

   return outset;
}

std::set< DetId > HDetIdAssociator::getMaxEDetId ( const std::set< DetId > &  inset, edm::Handle< CaloTowerCollection >  caloTowers  )

virtual

Definition at line 363 of file DetIdAssociator.cc.

References check_setup(), and triggerObjects_cff::id.

Referenced by getDetIdsCloseToAPoint().

{
// returns the most energetic tower in the NxN box (Michal)
   check_setup();
   std::set<DetId> outset;
   std::set<DetId>::const_iterator id_max = inset.begin();
   double Ehadmax=0;

   for(std::set<DetId>::const_iterator id_iter = inset.begin(); id_iter != inset.end(); id_iter++) {
     DetId id(*id_iter);
//     GlobalPoint point = getPosition(*id_iter);
//     int ieta = iEta(point);
//     int iphi = iPhi(point);
     CaloTowerCollection::const_iterator tower = (*caloTowers).find(id);
     if(tower != (*caloTowers).end() && tower->hadEnergy()>Ehadmax) {
       id_max = id_iter;
       Ehadmax = tower->hadEnergy();
     }
   }

   if (Ehadmax > 0) outset.insert(*id_max);

//   if (outset.size() > 0) {
//     std::cout<<" RRRA: Most energetic DetId:"<<std::endl;
//     for( std::set<DetId>::const_iterator itr=outset.begin(); itr!=outset.end(); itr++) {
//       GlobalPoint point = getPosition(*itr);
//       std::cout << "DetId: " <<itr->rawId() <<" (eta,phi): " << point.eta() << "," << point.phi()<<std::endl;
//     }
//   }

   return outset;
}

std::set< DetId > HDetIdAssociator::getMaxEDetId ( const std::set< DetId > &  inset, edm::Handle< HBHERecHitCollection >  recHits  )

virtual

Definition at line 398 of file DetIdAssociator.cc.

References check_setup(), and triggerObjects_cff::id.

{
// returns the most energetic tower in the NxN box - from RecHits (Michal)
   check_setup();
   std::set<DetId> outset;
   std::set<DetId>::const_iterator id_max = inset.begin();
   double Ehadmax=0;

   for(std::set<DetId>::const_iterator id_iter = inset.begin(); id_iter != inset.end(); id_iter++) {
     DetId id(*id_iter);
//     GlobalPoint point = getPosition(*id_iter);
//     int ieta = iEta(point);
//     int iphi = iPhi(point);
     HBHERecHitCollection::const_iterator hit = (*recHits).find(id);
     if(hit != (*recHits).end() && hit->energy()>Ehadmax) {
       id_max = id_iter;
       Ehadmax = hit->energy();
     }
   }

   if (Ehadmax > 0) outset.insert(*id_max);

//   if (outset.size() > 0) {
//     std::cout<<" RRRA: Most energetic DetId:"<<std::endl;
//     for( std::set<DetId>::const_iterator itr=outset.begin(); itr!=outset.end(); itr++) {
//       GlobalPoint point = getPosition(*itr);
//       std::cout << "DetId: " <<itr->rawId() <<" (eta,phi): " << point.eta() << "," << point.phi()<<std::endl;
//     }
//   }

   return outset;
}

virtual GlobalPoint HDetIdAssociator::getPosition ( const DetId &  )

protectedpure virtual

Implemented in HCaloDetIdAssociator.

Referenced by buildMap(), check_setup(), getDetIdsInACone(), and nearElement().

std::vector< GlobalPoint > HDetIdAssociator::getTrajectory ( const FreeTrajectoryState &  ftsStart, const std::vector< GlobalPoint > &  surfaces  )

virtual

Definition at line 25 of file DetIdAssociator.cc.

References funct::abs(), check_setup(), TCMET_cfi::corner, PV3DBase< T, PVType, FrameType >::eta(), TrajectoryStateOnSurface::freeState(), TrajectoryStateOnSurface::globalPosition(), TrajectoryStateOnSurface::isValid(), ivProp_, LogTrace, PV3DBase< T, PVType, FrameType >::mag(), FreeTrajectoryState::momentum(), PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), point, FreeTrajectoryState::position(), Propagator::propagate(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by ~HDetIdAssociator().

{
   check_setup();
   std::vector<GlobalPoint> trajectory;
   TrajectoryStateOnSurface tSOSDest;
   FreeTrajectoryState ftsCurrent = ftsStart;

   for(std::vector<GlobalPoint>::const_iterator surface_iter = surfaces.begin(); 
       surface_iter != surfaces.end(); surface_iter++) {
      // this stuff is some weird pointer, which destroy itself
      Cylinder *cylinder = new Cylinder(surface_iter->perp(), 
                                        Surface::PositionType(0,0,0),
                    Surface::RotationType() ); 
      Plane *forwardEndcap = new Plane(Surface::PositionType(0,0,surface_iter->z()),
                       Surface::RotationType());
      Plane *backwardEndcap = new Plane(Surface::PositionType(0,0,-surface_iter->z()),
                    Surface::RotationType());

      
      LogTrace("StartingPoint")<< "Propagate from "<< "\n"
    << "\tx: " << ftsStart.position().x()<< "\n"
    << "\ty: " << ftsStart.position().y()<< "\n"
    << "\tz: " << ftsStart.position().z()<< "\n"
    << "\tmomentum eta: " << ftsStart.momentum().eta()<< "\n"
    << "\tmomentum phi: " << ftsStart.momentum().phi()<< "\n"
    << "\tmomentum: " << ftsStart.momentum().mag()<< "\n";
     
         float tanTheta = ftsCurrent.momentum().perp()/ftsCurrent.momentum().z();
         float corner = surface_iter->perp()/surface_iter->z();
/*
      std::cout<<"Propagate from "<< "\n"
        << "\tx: " << ftsCurrent.position().x()<< "\n"
        << "\ty: " << ftsCurrent.position().y()<< "\n"
        << "\tz: " << ftsCurrent.position().z()<< "\n"
        << "\tz: " << ftsCurrent.position().perp()<< "\n"
        << "\tz: " << tanTheta<<" "<< corner <<"\n"
        << "\tmomentum eta: " << ftsCurrent.momentum().eta()<< "\n"
        << "\tmomentum phi: " << ftsCurrent.momentum().phi()<< "\n"
        << "\tmomentum: " << ftsCurrent.momentum().mag()<<std::endl;
*/ 
      // First propage the track to the cylinder if |eta|<1, othewise to the encap
      // and correct depending on the result
      int ibar = 0;
      if (fabs(tanTheta) > corner)
    {
                   tSOSDest = ivProp_->propagate(ftsCurrent, *cylinder); 
 //                  std::cout<<" Propagate to cylinder "<<std::endl;
        }
      else if(tanTheta > 0.)
    {tSOSDest = ivProp_->propagate(ftsCurrent, *forwardEndcap); ibar=1; }
      else
    {tSOSDest = ivProp_->propagate(ftsCurrent, *backwardEndcap); ibar=-1; }

//       std::cout<<" Trajectory valid? "<<tSOSDest.isValid()<<" First propagation in "<<ibar<<std::endl;

      if(! tSOSDest.isValid() )
      {
// barrel
       if(ibar == 0){ 
           if (tanTheta < 0 ) tSOSDest = ivProp_->propagate( ftsCurrent,*forwardEndcap);
           if (tanTheta >= 0 ) tSOSDest = ivProp_->propagate( ftsCurrent,*backwardEndcap);
       }
         else
         {
                tSOSDest = ivProp_->propagate(ftsCurrent, *cylinder);
         }
      } else
          {
// missed target
            if(abs(ibar) > 0)
            {
              if(tSOSDest.globalPosition().perp() > surface_iter->perp())
              {
                tSOSDest = ivProp_->propagate(ftsCurrent, *cylinder);
              }           
            }
              else
              {
                if (tanTheta < 0 ) tSOSDest = ivProp_->propagate( ftsCurrent,*forwardEndcap);
                if (tanTheta >= 0 ) tSOSDest = ivProp_->propagate( ftsCurrent,*backwardEndcap);
              }
          }


      // If missed the target, propagate to again
//      if ((!tSOSDest.isValid()) && point.perp() > surface_iter->perp())
//  {tSOSDest = ivProp_->propagate(ftsCurrent, *cylinder);std::cout<<" Propagate again 1 "<<std::endl;}
//         std::cout<<" Track is ok after repropagation to cylinder or not? "<<tSOSDest.isValid()<<std::endl;
//      if ((!tSOSDest.isValid()) && ftsStart.momentum().eta()>0. && fabs(ftsStart.momentum().eta())>1.) 
//  {tSOSDest = ivProp_->propagate(ftsStart, *forwardEndcap);std::cout<<" Propagate again 2 "<<std::endl;}
//       std::cout<<" Track is ok after repropagation forward or not? "<<tSOSDest.isValid()<<std::endl;
//      if ((!tSOSDest.isValid()) && ftsStart.momentum().eta()<0.&&fabs(ftsStart.momentum().eta())>1.) 
//  {tSOSDest = ivProp_->propagate(ftsStart, *backwardEndcap);std::cout<<" Propagate again 3 "<<std::endl;}
//       std::cout<<" Track is after repropagation backward ok or not? "<<tSOSDest.isValid()<<std::endl; 


      if (! tSOSDest.isValid()) return trajectory;
      
//      std::cout<<" Propagate reach something"<<std::endl; 
      LogTrace("SuccessfullPropagation") << "Great, I reached something." << "\n"
    << "\tx: " << tSOSDest.freeState()->position().x() << "\n"
    << "\ty: " << tSOSDest.freeState()->position().y() << "\n"
    << "\tz: " << tSOSDest.freeState()->position().z() << "\n"
    << "\teta: " << tSOSDest.freeState()->position().eta() << "\n"
    << "\tphi: " << tSOSDest.freeState()->position().phi() << "\n";

//      std::cout<<" The position of trajectory "<<tSOSDest.freeState()->position().perp()<<" "<<tSOSDest.freeState()->position().z()<<std::endl;  

      GlobalPoint point = tSOSDest.freeState()->position(); 
      point = tSOSDest.freeState()->position();
      ftsCurrent = *tSOSDest.freeState();
      trajectory.push_back(point);
   }
   return trajectory;
}

int HDetIdAssociator::iEta ( const GlobalPoint &  point )

virtual

Definition at line 197 of file DetIdAssociator.cc.

References PV3DBase< T, PVType, FrameType >::eta(), etaBinSize_, createfilelist::int, and nEta_.

Referenced by buildMap(), getDetIdsCloseToAPoint(), and getDetIdsInACone().

{
// unequal bin sizes for endcap, following HCAL geometry
   int iEta1 = int(point.eta()/etaBinSize_ + nEta_/2);
   if (point.eta()>1.827 && point.eta()<=1.830) return iEta1-1;
   else if (point.eta()>1.914 && point.eta()<=1.930) return iEta1-1;
   else if (point.eta()>2.001 && point.eta()<=2.043) return iEta1-1;
   else if (point.eta()>2.088 && point.eta()<=2.172) return iEta1-1;
   else if (point.eta()>2.175 && point.eta()<=2.262) return iEta1-1;
   else if (point.eta()>2.262 && point.eta()<=2.332) return iEta1-2;
   else if (point.eta()>2.332 && point.eta()<=2.349) return iEta1-1;
   else if (point.eta()>2.349 && point.eta()<=2.436) return iEta1-2;
   else if (point.eta()>2.436 && point.eta()<=2.500) return iEta1-3;
   else if (point.eta()>2.500 && point.eta()<=2.523) return iEta1-2;
   else if (point.eta()>2.523 && point.eta()<=2.610) return iEta1-3;
   else if (point.eta()>2.610 && point.eta()<=2.650) return iEta1-4;
   else if (point.eta()>2.650 && point.eta()<=2.697) return iEta1-3;
   else if (point.eta()>2.697 && point.eta()<=2.784) return iEta1-4;
   else if (point.eta()>2.784 && point.eta()<=2.868) return iEta1-5;
   else if (point.eta()>2.868 && point.eta()<=2.871) return iEta1-4;
   else if (point.eta()>2.871 && point.eta()<=2.958) return iEta1-5;
   else if (point.eta()>2.958) return iEta1-6;
   else if (point.eta()<-1.827 && point.eta()>=-1.830) return iEta1+1;
   else if (point.eta()<-1.914 && point.eta()>=-1.930) return iEta1+1;
   else if (point.eta()<-2.001 && point.eta()>=-2.043) return iEta1+1;
   else if (point.eta()<-2.088 && point.eta()>=-2.172) return iEta1+1;
   else if (point.eta()<-2.175 && point.eta()>=-2.262) return iEta1+1;
   else if (point.eta()<-2.262 && point.eta()>=-2.332) return iEta1+2;
   else if (point.eta()<-2.332 && point.eta()>=-2.349) return iEta1+1;
   else if (point.eta()<-2.349 && point.eta()>=-2.436) return iEta1+2;
   else if (point.eta()<-2.436 && point.eta()>=-2.500) return iEta1+3;
   else if (point.eta()<-2.500 && point.eta()>=-2.523) return iEta1+2;
   else if (point.eta()<-2.523 && point.eta()>=-2.610) return iEta1+3;
   else if (point.eta()<-2.610 && point.eta()>=-2.650) return iEta1+4;
   else if (point.eta()<-2.650 && point.eta()>=-2.697) return iEta1+3;
   else if (point.eta()<-2.697 && point.eta()>=-2.784) return iEta1+4;
   else if (point.eta()<-2.784 && point.eta()>=-2.868) return iEta1+5;
   else if (point.eta()<-2.868 && point.eta()>=-2.871) return iEta1+4;
   else if (point.eta()<-2.871 && point.eta()>=-2.958) return iEta1+5;
   else if (point.eta()<-2.349) return iEta1+6;
   else return iEta1;
}

virtual bool HDetIdAssociator::insideElement ( const GlobalPoint &  , const DetId &    )

protectedpure virtual

Implemented in HCaloDetIdAssociator.

Referenced by check_setup(), and getCrossedDetIds().

int HDetIdAssociator::iPhi ( const GlobalPoint &  point )

virtual

Definition at line 241 of file DetIdAssociator.cc.

References createfilelist::int, nPhi_, PV3DBase< T, PVType, FrameType >::phi(), and pi.

Referenced by buildMap(), getDetIdsCloseToAPoint(), and getDetIdsInACone().

{
   double pi=4*atan(1.);
   int iPhi1 = int((double(point.phi())+pi)/(2*pi)*nPhi_);
   return iPhi1;
}

virtual bool HDetIdAssociator::nearElement ( const GlobalPoint &  point, const DetId &  id, const double  distance  )

inlineprotectedvirtual

Definition at line 87 of file DetIdAssociator.h.

References SoftLeptonByDistance_cfi::distance, PV3DBase< T, PVType, FrameType >::eta(), getPosition(), PV3DBase< T, PVType, FrameType >::phi(), funct::pow(), and mathSSE::sqrt().

Referenced by getDetIdsInACone().

     {
    GlobalPoint center = getPosition(id);
    return sqrt(pow(point.eta()-center.eta(),2)+pow(point.phi()-center.phi(),2)) < distance;
     };

virtual void HDetIdAssociator::setPropagator ( Propagator *  ptr )

inlinevirtual

Definition at line 70 of file DetIdAssociator.h.

References ivProp_.

{   ivProp_ = ptr; };

Member Data Documentation

const double HDetIdAssociator::etaBinSize_

protected

Definition at line 96 of file DetIdAssociator.h.

Referenced by buildMap(), check_setup(), getDetIdsCloseToAPoint(), and iEta().

Propagator* HDetIdAssociator::ivProp_

protected

Definition at line 97 of file DetIdAssociator.h.

Referenced by check_setup(), getTrajectory(), and setPropagator().

const int HDetIdAssociator::nEta_

protected

Definition at line 95 of file DetIdAssociator.h.

Referenced by buildMap(), check_setup(), getDetIdsCloseToAPoint(), getDetIdsInACone(), and iEta().

const int HDetIdAssociator::nPhi_

protected

Definition at line 94 of file DetIdAssociator.h.

Referenced by buildMap(), check_setup(), getDetIdsCloseToAPoint(), getDetIdsInACone(), and iPhi().

std::vector<std::vector<std::set<DetId> > >* HDetIdAssociator::theMap_

protected

Definition at line 91 of file DetIdAssociator.h.

Referenced by buildMap(), getDetIdsCloseToAPoint(), and getDetIdsInACone().