|
|
|
#include <GenericTripletGenerator.h>
Public Member Functions | |
| void | clear () |
| GenericTripletGenerator (const edm::ParameterSet &conf) | |
| virtual const OrderedSeedingHits & | run (const TrackingRegion ®ion, const edm::Event &ev, const edm::EventSetup &es) |
| virtual | ~GenericTripletGenerator () |
Private Member Functions | |
| ctfseeding::SeedingLayerSets | init (const edm::EventSetup &es) |
| std::pair< bool, float > | qualityFilter (const OrderedHitTriplet &oht, const std::map< float, OrderedHitTriplet > &map, const ctfseeding::SeedingLayers &ls) const |
Private Attributes | |
| OrderedHitTriplets | hitTriplets |
| SeedingLayerSetsBuilder | theLsb |
Definition at line 16 of file GenericTripletGenerator.h.
| GenericTripletGenerator::GenericTripletGenerator | ( | const edm::ParameterSet & | conf | ) |
Definition at line 11 of file GenericTripletGenerator.cc.
:
//conf_(conf),
theLsb(conf.getParameter<edm::ParameterSet>("LayerPSet")){
edm::LogInfo("CtfSpecialSeedGenerator|GenericTripletGenerator") << "Constructing GenericTripletGenerator";
}
| virtual GenericTripletGenerator::~GenericTripletGenerator | ( | ) | [inline, virtual] |
Definition at line 19 of file GenericTripletGenerator.h.
{};
| void GenericTripletGenerator::clear | ( | void | ) | [inline, virtual] |
Reimplemented from OrderedHitsGenerator.
Definition at line 23 of file GenericTripletGenerator.h.
References hitTriplets.
{hitTriplets.clear();}
| SeedingLayerSets GenericTripletGenerator::init | ( | const edm::EventSetup & | es | ) | [private] |
Definition at line 18 of file GenericTripletGenerator.cc.
References SeedingLayerSetsBuilder::layers(), and theLsb.
Referenced by run().
{
//edm::ParameterSet leyerPSet = conf_.getParameter<edm::ParameterSet>("LayerPSet");
//SeedingLayerSetsBuilder lsBuilder(leyerPSet);
SeedingLayerSets lss = theLsb.layers(es);
return lss;
}
| std::pair< bool, float > GenericTripletGenerator::qualityFilter | ( | const OrderedHitTriplet & | oht, |
| const std::map< float, OrderedHitTriplet > & | map, | ||
| const ctfseeding::SeedingLayers & | ls | ||
| ) | const [private] |
Definition at line 77 of file GenericTripletGenerator.cc.
References OrderedHitTriplet::inner(), PV3DBase< T, PVType, FrameType >::mag(), mag(), OrderedHitTriplet::middle(), OrderedHitTriplet::outer(), PV3DBase< T, PVType, FrameType >::phi(), FastCircle::rho(), mathSSE::sqrt(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().
Referenced by run().
{
//first check the radius
GlobalPoint innerpos = oht.inner()->globalPosition();
GlobalPoint middlepos = oht.middle()->globalPosition();
GlobalPoint outerpos = oht.outer()->globalPosition();
std::vector<const TrackingRecHit*> ohttrh;
ohttrh.push_back(&(*(oht.inner()))); ohttrh.push_back(&(*(oht.middle()))); ohttrh.push_back(&(*(oht.outer())));
std::vector<const TrackingRecHit*>::const_iterator ioht;
//now chech that the change in phi is reasonable. the limiting angular distance is the one in case
//one of the two points is a tangent point.
float limit_phi_distance1 = sqrt((middlepos.x()-outerpos.x())*(middlepos.x()-outerpos.x()) +
(middlepos.y()-outerpos.y())*(middlepos.y()-outerpos.y()))/middlepos.mag();//actually this is the tangent of the limiting angle
float limit_phi_distance2 = sqrt((middlepos.x()-innerpos.x())*(middlepos.x()-innerpos.x()) +
(middlepos.y()-innerpos.y())*(middlepos.y()-innerpos.y()))/innerpos.mag();
//if (fabs(tan(outerpos.phi()-middlepos.phi()))>limit_phi_distance1 ||
// fabs(tan(innerpos.phi()-middlepos.phi()))>limit_phi_distance2) {
if (fabs(outerpos.phi()-middlepos.phi())>fabs(atan(limit_phi_distance1)) ||
fabs(innerpos.phi()-middlepos.phi())>fabs(atan(limit_phi_distance2)) ) {
//std::cout << "rejected because phi" << std::endl;
return std::make_pair(false, 0.);
}
//should we add a control on the r-z plane too?
/*
//now check that there is no big change in the r-z projection
float dz1 = outerpos.z()-middlepos.z();
float dr1 = sqrt(outerpos.x()*outerpos.x()+outerpos.y()*outerpos.y())-
sqrt(middlepos.x()*middlepos.x()+middlepos.y()*middlepos.y());
float dz2 = middlepos.z()-innerpos.z();
float dr2 = sqrt(middlepos.x()*middlepos.x()+middlepos.y()*middlepos.y())-
sqrt(innerpos.x()*innerpos.x()+innerpos.y()*innerpos.y());
float tan1 = dz1/dr1;
float tan2 = dz2/dr2;
//how much should we allow? should we make it configurable?
if (fabs(tan1-tan2)/tan1>0.5){
//std::cout << "rejected because z" << std::endl;
return std::make_pair(false, 0.);
}
*/
//now check the radius is not too small
FastCircle circle(innerpos, middlepos, outerpos);
if (circle.rho() < 200 && circle.rho() != 0) return std::make_pair(false, circle.rho()); //to small radius
//now check if at least 2 hits are shared with an existing triplet
//look for similar radii in the map
std::map<float, OrderedHitTriplet>::const_iterator lower_bound = map.lower_bound((1-0.01)*circle.rho());
std::map<float, OrderedHitTriplet>::const_iterator upper_bound = map.upper_bound((1+0.01)*circle.rho());
std::map<float, OrderedHitTriplet>::const_iterator iter;
for (iter = lower_bound; iter != upper_bound && iter->first <= upper_bound->first; iter++){
int shared=0;
std::vector<const TrackingRecHit*> curtrh;
curtrh.push_back(&*(iter->second.inner()));curtrh.push_back(&*(iter->second.middle()));curtrh.push_back(&*(iter->second.outer()));
std::vector<const TrackingRecHit*>::const_iterator curiter;
for (curiter = curtrh.begin(); curiter != curtrh.end(); curiter++){
for (ioht = ohttrh.begin(); ioht != ohttrh.end(); ioht++){
if ((*ioht)->geographicalId()==(*curiter)->geographicalId() &&
((*ioht)->localPosition()-(*curiter)->localPosition()).mag()<1e-5) shared++;
}
}
if (shared>1) return std::make_pair(false, circle.rho());
}
return std::make_pair(true,circle.rho());
}
| const OrderedSeedingHits & GenericTripletGenerator::run | ( | const TrackingRegion & | region, |
| const edm::Event & | ev, | ||
| const edm::EventSetup & | es | ||
| ) | [virtual] |
Implements OrderedHitsGenerator.
Definition at line 26 of file GenericTripletGenerator.cc.
References Exception, TrackingRegion::hits(), hitTriplets, init(), python::rootplot::utilities::ls(), and qualityFilter().
{
hitTriplets.clear();
hitTriplets.reserve(0);
SeedingLayerSets lss = init(es);
SeedingLayerSets::const_iterator iLss;
std::map<float, OrderedHitTriplet> radius_triplet_map;
for (iLss = lss.begin(); iLss != lss.end(); iLss++){
SeedingLayers ls = *iLss;
if (ls.size() != 3){
throw cms::Exception("CtfSpecialSeedGenerator") << "You are using " << ls.size() <<" layers in set instead of 3 ";
}
std::vector<SeedingHit> innerHits = region.hits(e, es, &ls[0]);
//std::cout << "innerHits.size()=" << innerHits.size() << std::endl;
std::vector<SeedingHit> middleHits = region.hits(e, es, &ls[1]);
//std::cout << "middleHits.size()=" << middleHits.size() << std::endl;
std::vector<SeedingHit> outerHits = region.hits(e, es, &ls[2]);
//std::cout << "outerHits.size()=" << outerHits.size() << std::endl;
//std::cout << "trying " << innerHits.size()*middleHits.size()*outerHits.size() << " combinations "<<std::endl;
std::vector<SeedingHit>::const_iterator iOuterHit;
for (iOuterHit = outerHits.begin(); iOuterHit != outerHits.end(); iOuterHit++){
std::vector<SeedingHit>::const_iterator iMiddleHit;
for (iMiddleHit = middleHits.begin(); iMiddleHit != middleHits.end(); iMiddleHit++){
std::vector<SeedingHit>::const_iterator iInnerHit;
for (iInnerHit = innerHits.begin(); iInnerHit != innerHits.end(); iInnerHit++){
//GlobalPoint innerpos = ls[0].hitBuilder()->build(&(**iInnerHit))->globalPosition();
//GlobalPoint middlepos = ls[1].hitBuilder()->build(&(**iMiddleHit))->globalPosition();
//GlobalPoint outerpos = ls[2].hitBuilder()->build(&(**iOuterHit))->globalPosition();
//FastCircle circle(innerpos,
// middlepos,
// outerpos);
//do a first check on the radius of curvature to reduce the combinatorics
OrderedHitTriplet oht(*iInnerHit,*iMiddleHit,*iOuterHit);
std::pair<bool,float> val_radius = qualityFilter(oht,radius_triplet_map,ls);
if (val_radius.first){
//if (circle.rho() > 200 || circle.rho() == 0) { //0 radius means straight line
//hitTriplets.push_back(OrderedHitTriplet(*iInnerHit,
// *iMiddleHit,
// *iOuterHit));
hitTriplets.push_back(oht);
radius_triplet_map.insert(std::make_pair(val_radius.second,oht));
}
}
}
}
}
//std::cout << "ending with " << hitTriplets.size() << " triplets" << std::endl;
return hitTriplets;
}
Definition at line 31 of file GenericTripletGenerator.h.
Definition at line 30 of file GenericTripletGenerator.h.
Referenced by init().
1.7.3