---

HITrackingRegionProducer Class Reference

#include <RecoHI/HiTracking/plugins/HITrackingRegionProducer.h>

Inheritance diagram for HITrackingRegionProducer:

List of all members.

Public Member Functions
int	estimateMultiplicity (const edm::Event &ev, const edm::EventSetup &es) const
	HITrackingRegionProducer (const edm::ParameterSet &cfg)
virtual std::vector < TrackingRegion * >	regions (const edm::Event &ev, const edm::EventSetup &es) const
virtual	~HITrackingRegionProducer ()
Private Attributes
GlobalVector	theDirection
GlobalPoint	theOrigin
double	theOriginHalfLength
double	theOriginRadius
bool	thePrecise
double	thePtMin
std::string	theSiPixelRecHits

---

Detailed Description

Definition at line 18 of file HITrackingRegionProducer.h.

---

Constructor & Destructor Documentation

HITrackingRegionProducer::HITrackingRegionProducer

(

const edm::ParameterSet &

cfg

)

[inline]

Definition at line 22 of file HITrackingRegionProducer.h.

References edm::ParameterSet::getParameter(), theDirection, theOrigin, theOriginHalfLength, theOriginRadius, thePrecise, thePtMin, and theSiPixelRecHits.

                                                       { 

    edm::ParameterSet regionPSet = cfg.getParameter<edm::ParameterSet>("RegionPSet");

    thePtMin            = regionPSet.getParameter<double>("ptMin");
    theOriginRadius     = regionPSet.getParameter<double>("originRadius");
    theOriginHalfLength = regionPSet.getParameter<double>("originHalfLength");
    double xPos         = regionPSet.getParameter<double>("originXPos");
    double yPos         = regionPSet.getParameter<double>("originYPos");
    double zPos         = regionPSet.getParameter<double>("originZPos");
    double xDir         = regionPSet.getParameter<double>("directionXCoord");
    double yDir         = regionPSet.getParameter<double>("directionYCoord");
    double zDir         = regionPSet.getParameter<double>("directionZCoord");
    thePrecise          = regionPSet.getParameter<bool>("precise"); 
    theSiPixelRecHits   = regionPSet.getParameter<std::string>("siPixelRecHits");
    theOrigin = GlobalPoint(xPos,yPos,zPos);
    theDirection = GlobalVector(xDir, yDir, zDir);
  }   

virtual HITrackingRegionProducer::~HITrackingRegionProducer

(

)

[inline, virtual]

Definition at line 41 of file HITrackingRegionProducer.h.

{}

---

Member Function Documentation

int HITrackingRegionProducer::estimateMultiplicity	(	const edm::Event &	ev,
		const edm::EventSetup &	es
	)			const [inline]

Definition at line 44 of file HITrackingRegionProducer.h.

References DetId::det(), detId, edm::Event::getByLabel(), PXBDetId::layer(), edm::Handle< T >::product(), DetId::subdetId(), and theSiPixelRecHits.

Referenced by regions().

  {
    //rechits
    edm::Handle<SiPixelRecHitCollection> recHitColl;
    ev.getByLabel(theSiPixelRecHits, recHitColl);
 
    SiPixelRecHitCollection::id_iterator recHitIdIterator;
    SiPixelRecHitCollection::id_iterator recHitIdIteratorBegin = (recHitColl.product())->id_begin();
    SiPixelRecHitCollection::id_iterator recHitIdIteratorEnd   = (recHitColl.product())->id_end();

    int numRecHits = 0;
    for(recHitIdIterator = recHitIdIteratorBegin;
        recHitIdIterator != recHitIdIteratorEnd; recHitIdIterator++) {

         DetId detId = DetId((*recHitIdIterator).rawId()); // Get the Detid object
         unsigned int detType=detId.det();    // det type, tracker=1
         unsigned int subid=detId.subdetId(); //subdetector type, barrel=1, fpix=2
         PXBDetId pdetId = PXBDetId(detId);
         unsigned int layer=0;
         layer=pdetId.layer();
         if(detType==1 && subid==1 && layer==1) {
           SiPixelRecHitCollection::range pixelrechitRange = (recHitColl.product())->get(*recHitIdIterator);
           SiPixelRecHitCollection::const_iterator pixelrechitRangeIteratorBegin = pixelrechitRange.first;
           SiPixelRecHitCollection::const_iterator pixelrechitRangeIteratorEnd = pixelrechitRange.second;
           SiPixelRecHitCollection::const_iterator pixeliter;
           for (pixeliter = pixelrechitRangeIteratorBegin ;
              pixeliter != pixelrechitRangeIteratorEnd; ++pixeliter) {
              numRecHits++;
           }
         }
    }
    return numRecHits;
  }

virtual std::vector<TrackingRegion* > HITrackingRegionProducer::regions	(	const edm::Event &	ev,
		const edm::EventSetup &	es
	)			const [inline, virtual]

Implements TrackingRegionProducer.

Definition at line 78 of file HITrackingRegionProducer.h.

References aa, estimateMultiplicity(), LogTrace, Pi, HLT_VtxMuL3::result, theDirection, theOrigin, theOriginHalfLength, theOriginRadius, thePrecise, and thePtMin.

                                                                                                 {

    int estMult = estimateMultiplicity(ev, es);

    // fit from MC information
    float aa = 1.90935e-04;
    float bb = -2.90167e-01;
    float cc = 3.86125e+02;

    float estTracks = aa*estMult*estMult+bb*estMult+cc;

    LogTrace("heavyIonHLTVertexing")<<"[HIVertexing]";
    LogTrace("heavyIonHLTVertexing")<<" [HIVertexing: hits in the 1. layer:" << estMult << "]";
    LogTrace("heavyIonHLTVertexing")<<" [HIVertexing: estimated number of tracks:" << estTracks << "]";

    float regTracking = 400.;  //if we have more tracks -> regional tracking
    float etaB = 10.;
    float phiB = TMath::Pi()/2.;

    float decEta = estTracks/600.;
    etaB = 2.5/decEta;

    if(estTracks>regTracking) {
      LogTrace("heavyIonHLTVertexing")<<" [HIVertexing: Regional Tracking]";
      LogTrace("heavyIonHLTVertexing")<<"  [Regional Tracking: eta range: -" << etaB << ", "<< etaB <<"]";
      LogTrace("heavyIonHLTVertexing")<<"  [Regional Tracking: phi range: -" << phiB << ", "<< phiB <<"]";
      LogTrace("heavyIonHLTVertexing")<<"  [Regional Tracking: factor of decrease: " << decEta*2. << "]";  // 2:from phi
    }

    // tracking region selection
    std::vector<TrackingRegion* > result;
    if(estTracks>regTracking) {  // regional tracking
      result.push_back( 
          new RectangularEtaPhiTrackingRegion(theDirection, theOrigin, thePtMin, theOriginRadius, theOriginHalfLength, etaB, phiB, thePrecise) );
    }
    else {                       // global tracking
      LogTrace("heavyIonHLTVertexing")<<" [HIVertexing: Global Tracking]";
      result.push_back( 
          new GlobalTrackingRegion(thePtMin, theOrigin, theOriginRadius, theOriginHalfLength, thePrecise) );
    }
    return 
result;
  }

---

Member Data Documentation

GlobalVector HITrackingRegionProducer::theDirection [private]

Definition at line 129 of file HITrackingRegionProducer.h.

Referenced by HITrackingRegionProducer(), and regions().

GlobalPoint HITrackingRegionProducer::theOrigin [private]

Definition at line 125 of file HITrackingRegionProducer.h.

Referenced by HITrackingRegionProducer(), and regions().

double HITrackingRegionProducer::theOriginHalfLength [private]

Definition at line 127 of file HITrackingRegionProducer.h.

Referenced by HITrackingRegionProducer(), and regions().

double HITrackingRegionProducer::theOriginRadius [private]

Definition at line 126 of file HITrackingRegionProducer.h.

Referenced by HITrackingRegionProducer(), and regions().

bool HITrackingRegionProducer::thePrecise [private]

Definition at line 128 of file HITrackingRegionProducer.h.

Referenced by HITrackingRegionProducer(), and regions().

double HITrackingRegionProducer::thePtMin [private]

Definition at line 124 of file HITrackingRegionProducer.h.

Referenced by HITrackingRegionProducer(), and regions().

std::string HITrackingRegionProducer::theSiPixelRecHits [private]

Definition at line 123 of file HITrackingRegionProducer.h.

Referenced by estimateMultiplicity(), and HITrackingRegionProducer().

---

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

• RecoHI/HiTracking/plugins/HITrackingRegionProducer.h

---