|
|
|
#include <PixelTripletNoTipGenerator.h>
Public Member Functions | |
| virtual void | hitTriplets (const TrackingRegion ®ion, OrderedHitTriplets &trs, const edm::Event &ev, const edm::EventSetup &es) |
| virtual void | init (const HitPairGenerator &pairs, const std::vector< ctfseeding::SeedingLayer > &layers, LayerCacheType *layerCache) |
| const HitPairGenerator & | pairGenerator () const |
| PixelTripletNoTipGenerator (const edm::ParameterSet &cfg) | |
| const std::vector < ctfseeding::SeedingLayer > & | thirdLayers () const |
| virtual | ~PixelTripletNoTipGenerator () |
Private Types | |
| typedef CombinedHitTripletGenerator::LayerCacheType | LayerCacheType |
Private Attributes | |
| float | extraHitPhiToleranceForPreFiltering |
| float | extraHitRPhitolerance |
| float | extraHitRZtolerance |
| double | theChi2Cut |
| LayerCacheType * | theLayerCache |
| std::vector < ctfseeding::SeedingLayer > | theLayers |
| double | theNSigma |
| HitPairGenerator * | thePairGenerator |
Definition at line 15 of file PixelTripletNoTipGenerator.h.
typedef CombinedHitTripletGenerator::LayerCacheType PixelTripletNoTipGenerator::LayerCacheType [private] |
Reimplemented from HitTripletGeneratorFromPairAndLayers.
Definition at line 16 of file PixelTripletNoTipGenerator.h.
| PixelTripletNoTipGenerator::PixelTripletNoTipGenerator | ( | const edm::ParameterSet & | cfg | ) |
Definition at line 25 of file PixelTripletNoTipGenerator.cc.
: thePairGenerator(0), theLayerCache(0), extraHitRZtolerance(cfg.getParameter<double>("extraHitRZtolerance")), extraHitRPhitolerance(cfg.getParameter<double>("extraHitRPhitolerance")), extraHitPhiToleranceForPreFiltering(cfg.getParameter<double>("extraHitPhiToleranceForPreFiltering")), theNSigma(cfg.getParameter<double>("nSigma")), theChi2Cut(cfg.getParameter<double>("chi2Cut")) { }
| virtual PixelTripletNoTipGenerator::~PixelTripletNoTipGenerator | ( | ) | [inline, virtual] |
Definition at line 20 of file PixelTripletNoTipGenerator.h.
References thePairGenerator.
{ delete thePairGenerator; }
| void PixelTripletNoTipGenerator::hitTriplets | ( | const TrackingRegion & | region, |
| OrderedHitTriplets & | trs, | ||
| const edm::Event & | ev, | ||
| const edm::EventSetup & | es | ||
| ) | [virtual] |
Implements HitTripletGenerator.
Definition at line 44 of file PixelTripletNoTipGenerator.cc.
References Reference_intrackfit_cff::barrel, PixelRecoUtilities::curvature(), extraHitPhiToleranceForPreFiltering, extraHitRZtolerance, HitPairGenerator::hitPairs(), RecHitsSortedInPhi::hits(), HitPairGeneratorFromLayerPair::innerLayer(), PixelRecoUtilities::inversePt(), geometryCSVtoXML::line, DetLayer::location(), M_PI, TrackingRegion::origin(), TrackingRegion::originRBound(), HitPairGeneratorFromLayerPair::outerLayer(), p1, p2, p3, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), GeometricSearchDet::position(), TrackingRegion::ptMin(), CosmicsPD_Skims::radius, edm::shift, findQualityFiles::size, OrderedHitTriplets::size(), OrderedHitPairs::size(), PixelRecoRange< T >::sort(), BarrelDetLayer::specificSurface(), funct::sqr(), mathSSE::sqrt(), BarrelDetLayer::surface(), ForwardDetLayer::surface(), swap(), theChi2Cut, theLayers, OrderedHitsGenerator::theMaxElement, theNSigma, thePairGenerator, and PV3DBase< T, PVType, FrameType >::z().
{
//
//edm::Handle<reco::BeamSpot> bsHandle;
//ev.getByLabel( theBeamSpotTag, bsHandle);
//if(!bsHandle.isValid()) return;
//const reco::BeamSpot & bs = *bsHandle;
//double errorXY = sqrt( sqr(bs.BeamWidthX()) + sqr(bs.BeamWidthY()) );
//
GlobalPoint bsPoint = region.origin();
double errorXY = region.originRBound();
GlobalVector shift = bsPoint - GlobalPoint(0.,0.,0.);
OrderedHitPairs pairs; pairs.reserve(30000);
OrderedHitPairs::const_iterator ip;
thePairGenerator->hitPairs(region,pairs,ev,es);
if (pairs.size() ==0) return;
int size = theLayers.size();
const RecHitsSortedInPhi **thirdHitMap = new const RecHitsSortedInPhi*[size];
for (int il=0; il <=size-1; il++) {
thirdHitMap[il] = &(*theLayerCache)(&theLayers[il], region, ev, es);
}
const HitPairGeneratorFromLayerPair * pairGen = dynamic_cast<const HitPairGeneratorFromLayerPair *>(thePairGenerator);
const DetLayer * firstLayer = pairGen->innerLayer().detLayer();
const DetLayer * secondLayer = pairGen->outerLayer().detLayer();
if (!firstLayer || !secondLayer) return;
MultipleScatteringParametrisation sigma1RPhi( firstLayer, es);
MultipleScatteringParametrisation sigma2RPhi( secondLayer, es);
typedef RecHitsSortedInPhi::Hit Hit;
for (ip = pairs.begin(); ip != pairs.end(); ip++) {
GlobalPoint p1((*ip).inner()->globalPosition()-shift);
GlobalPoint p2((*ip).outer()->globalPosition()-shift);
ThirdHitPredictionFromInvLine predictionRPhiTMP(p1, p2 );
double pt_p1p2 = 1./PixelRecoUtilities::inversePt(predictionRPhiTMP.curvature(),es);
PixelRecoPointRZ point1(p1.perp(), p1.z());
PixelRecoPointRZ point2(p2.perp(), p2.z());
PixelRecoLineRZ line(point1, point2);
double msRPhi1 = sigma1RPhi(pt_p1p2, line.cotLine(), 0.f);
double msRPhi2 = sigma2RPhi(pt_p1p2, line.cotLine(),point1);
double sinTheta = 1/sqrt(1+sqr(line.cotLine()));
double cosTheta = fabs(line.cotLine())/sqrt(1+sqr(line.cotLine()));
double p1_errorRPhi = sqrt(sqr((*ip).inner()->errorGlobalRPhi())+sqr(msRPhi1) +sqr(errorXY));
double p2_errorRPhi = sqrt(sqr((*ip).outer()->errorGlobalRPhi())+sqr(msRPhi2) +sqr(errorXY));
ThirdHitPredictionFromInvLine predictionRPhi(p1, p2, p1_errorRPhi, p2_errorRPhi );
for (int il=0; il <=size-1; il++) {
const DetLayer * layer = theLayers[il].detLayer();
bool barrelLayer = (layer->location() == GeomDetEnumerators::barrel);
MultipleScatteringParametrisation sigma3RPhi( layer, es);
double msRPhi3 = sigma3RPhi(pt_p1p2, line.cotLine(),point2);
Range rRange;
if (barrelLayer) {
const BarrelDetLayer& bl = dynamic_cast<const BarrelDetLayer&>(*layer);
float halfThickness = bl.surface().bounds().thickness()/2;
float radius = bl.specificSurface().radius();
rRange = Range(radius-halfThickness, radius+halfThickness);
} else {
const ForwardDetLayer& fl = dynamic_cast<const ForwardDetLayer&>(*layer);
float halfThickness = fl.surface().bounds().thickness()/2;
float zLayer = fl.position().z() ;
float zMin = zLayer-halfThickness;
float zMax = zLayer+halfThickness;
GlobalVector dr = p2-p1;
GlobalPoint p3_a = p2+dr*(zMin-p2.z())/dr.z();
GlobalPoint p3_b = p2+dr*(zMax-p2.z())/dr.z();
if (zLayer * p3_a.z() < 0) continue;
rRange = Range(p3_a.perp(), p3_b.perp());
rRange.sort();
}
double displacment = shift.perp();
GlobalPoint crossing1 = predictionRPhi.crossing(rRange.min()-displacment)+shift;
GlobalPoint crossing2 = predictionRPhi.crossing(rRange.max()+displacment)+shift;
float c1_phi= crossing1.phi();
float c2_phi= crossing2.phi();
if (c2_phi < c1_phi) swap(c1_phi,c2_phi);
if (c2_phi-c1_phi > M_PI) { c2_phi -= 2*M_PI; swap(c1_phi,c2_phi); }
double extraAngle = (displacment+theNSigma*msRPhi3)/rRange.min()+extraHitPhiToleranceForPreFiltering;
c1_phi -= extraAngle;
c2_phi += extraAngle;
vector<Hit> thirdHits = thirdHitMap[il]->hits(c1_phi, c2_phi) ;
typedef vector<Hit>::const_iterator IH;
for (IH th=thirdHits.begin(), eh=thirdHits.end(); th < eh; ++th) {
GlobalPoint p3((*th)->globalPosition()-shift);
double p3_errorRPhi = sqrt(sqr((*th)->errorGlobalRPhi()) +sqr(msRPhi3) + sqr(errorXY));
predictionRPhi.add(p3,p3_errorRPhi);
double curvature = predictionRPhi.curvature();
ThirdHitZPrediction zPrediction( (*ip).inner()->globalPosition(), sqrt(sqr((*ip).inner()->errorGlobalR())+sqr(msRPhi1/cosTheta)), sqrt( sqr((*ip).inner()->errorGlobalZ())+ sqr(msRPhi1/sinTheta)),
(*ip).outer()->globalPosition(), sqrt(sqr((*ip).outer()->errorGlobalR())+sqr(msRPhi2/cosTheta)), sqrt( sqr((*ip).outer()->errorGlobalZ())+sqr(msRPhi2/sinTheta)),
curvature, theNSigma);
ThirdHitZPrediction::Range zRange = zPrediction((*th)->globalPosition(), sqrt(sqr((*th)->errorGlobalR()))+sqr(msRPhi3/cosTheta));
double z3Hit = (*th)->globalPosition().z();
double z3HitError = theNSigma*(sqrt(sqr((*th)->errorGlobalZ()) + sqr(msRPhi3/sinTheta) ))+extraHitRZtolerance;
Range hitZRange(z3Hit-z3HitError, z3Hit+z3HitError);
bool inside = hitZRange.hasIntersection(zRange);
double curvatureMS = PixelRecoUtilities::curvature(1./region.ptMin(),es);
bool ptCut = (predictionRPhi.curvature()-theNSigma*predictionRPhi.errorCurvature() < curvatureMS);
bool chi2Cut = (predictionRPhi.chi2() < theChi2Cut);
if (inside && ptCut && chi2Cut) {
if (theMaxElement!=0 && result.size() >= theMaxElement){
result.clear();
edm::LogError("TooManyTriplets")<<" number of triples exceed maximum. no triplets produced.";
delete [] thirdHitMap;
return;
}
result.push_back( OrderedHitTriplet( (*ip).inner(), (*ip).outer(), *th));
}
predictionRPhi.remove(p3,p3_errorRPhi);
}
}
}
delete [] thirdHitMap;
}
| void PixelTripletNoTipGenerator::init | ( | const HitPairGenerator & | pairs, |
| const std::vector< ctfseeding::SeedingLayer > & | layers, | ||
| LayerCacheType * | layerCache | ||
| ) | [virtual] |
Implements HitTripletGeneratorFromPairAndLayers.
Definition at line 35 of file PixelTripletNoTipGenerator.cc.
References HitPairGenerator::clone(), theLayerCache, theLayers, and thePairGenerator.
{
thePairGenerator = pairs.clone();
theLayers = layers;
theLayerCache = layerCache;
}
| const HitPairGenerator& PixelTripletNoTipGenerator::pairGenerator | ( | ) | const [inline] |
Definition at line 28 of file PixelTripletNoTipGenerator.h.
References thePairGenerator.
{ return *thePairGenerator; }
| const std::vector<ctfseeding::SeedingLayer>& PixelTripletNoTipGenerator::thirdLayers | ( | ) | const [inline] |
Definition at line 29 of file PixelTripletNoTipGenerator.h.
References theLayers.
{ return theLayers; }
float PixelTripletNoTipGenerator::extraHitPhiToleranceForPreFiltering [private] |
Definition at line 37 of file PixelTripletNoTipGenerator.h.
Referenced by hitTriplets().
float PixelTripletNoTipGenerator::extraHitRPhitolerance [private] |
Definition at line 36 of file PixelTripletNoTipGenerator.h.
float PixelTripletNoTipGenerator::extraHitRZtolerance [private] |
Definition at line 35 of file PixelTripletNoTipGenerator.h.
Referenced by hitTriplets().
double PixelTripletNoTipGenerator::theChi2Cut [private] |
Definition at line 39 of file PixelTripletNoTipGenerator.h.
Referenced by hitTriplets().
Definition at line 34 of file PixelTripletNoTipGenerator.h.
Referenced by init().
std::vector<ctfseeding::SeedingLayer> PixelTripletNoTipGenerator::theLayers [private] |
Definition at line 33 of file PixelTripletNoTipGenerator.h.
Referenced by hitTriplets(), init(), and thirdLayers().
double PixelTripletNoTipGenerator::theNSigma [private] |
Definition at line 38 of file PixelTripletNoTipGenerator.h.
Referenced by hitTriplets().
Definition at line 32 of file PixelTripletNoTipGenerator.h.
Referenced by hitTriplets(), init(), pairGenerator(), and ~PixelTripletNoTipGenerator().
1.7.3