31 int getContainedHits(
const std::vector<VertexHit> &hits,
double z0,
double &chi);
59 inputTag_ (config.getParameter<edm::InputTag>(
"inputTag")),
60 minZ_ (config.getParameter<double>(
"minZ")),
61 maxZ_ (config.getParameter<double>(
"maxZ")),
62 zStep_ (config.getParameter<double>(
"zStep")),
63 clusterPars_ (config.getParameter< std::vector<double> >(
"clusterPars")),
64 nhitsTrunc_ (config.getParameter<int>(
"nhitsTrunc")),
65 clusterTrunc_ (config.getParameter<double>(
"clusterTrunc"))
102 std::vector<VertexHit> vhits;
114 std::vector<SiPixelCluster::Pixel> pixels(
hit->cluster()->pixels());
115 bool pixelOnEdge =
false;
116 for(std::vector<SiPixelCluster::Pixel>::const_iterator pixel = pixels.begin();
117 pixel != pixels.end(); ++pixel) {
118 int pixelX = pixel->x;
119 int pixelY = pixel->y;
130 hit->localPosition().
y(),
131 hit->localPosition().
z());
136 vh.
w =
hit->cluster()->sizeY();
142 int nhits = 0, nhits_max = 0;
143 double chi = 0, chi_max = 1
e+9;
148 if(nhits > nhits_max) {
152 if(nhits >= nhits_max && chi < chi_max) {
159 int nbest=0, nminus=0, nplus=0;
164 double clusVtxQual=0.0;
165 if ((nminus+nplus)> 0)
166 clusVtxQual = (2.0*nbest)/(nminus+nplus);
168 clusVtxQual = 1000.0;
182 if (clusVtxQual < polyCut) accept =
false;
195 for(std::vector<VertexHit>::const_iterator
hit = hits.begin();
hit!= hits.end();
hit++) {
196 double p = 2 * fabs(
hit->
z - z0)/
hit->r + 0.5;
197 if(fabs(p -
hit->w) <= 1.) {
198 chi += fabs(p -
hit->w);
boost::transform_iterator< IterHelp, const_IdIter > const_iterator
#define DEFINE_FWK_MODULE(type)
GlobalPoint toGlobal(const Local2DPoint &lp) const
Conversion to the global R.F. from the R.F. of the GeomDet.
bool accept(const edm::Event &event, const edm::TriggerResults &triggerTable, const std::string &triggerPath)
virtual bool isItEdgePixelInX(int ixbin) const =0
std::vector< double > clusterPars_
data_type const * data(size_t cell) const
How EventSelector::AcceptEvent() decides whether to accept an event for output otherwise it is excluding the probing of A single or multiple positive and the trigger will pass if any such matching triggers are PASS or EXCEPTION[A criterion thatmatches no triggers at all is detected and causes a throw.] A single negative with an expectation of appropriate bit checking in the decision and the trigger will pass if any such matching triggers are FAIL or EXCEPTION A wildcarded negative criterion that matches more than one trigger in the trigger but the state exists so we define the behavior If all triggers are the negative crieriion will lead to accepting the event(this again matches the behavior of"!*"before the partial wildcard feature was incorporated).The per-event"cost"of each negative criterion with multiple relevant triggers is about the same as!*was in the past
virtual const GeomDet * idToDet(DetId) const
int getContainedHits(const std::vector< VertexHit > &hits, double z0, double &chi)
virtual const PixelTopology & specificTopology() const
Returns a reference to the pixel proxy topology.
T const * product() const
void addCollectionTag(const edm::InputTag &collectionTag)
collectionTags
T const * product() const
virtual bool isItEdgePixelInY(int iybin) const =0
virtual bool hltFilter(edm::Event &, const edm::EventSetup &, trigger::TriggerFilterObjectWithRefs &filterproduct)
HLTPixelClusterShapeFilter(const edm::ParameterSet &)
~HLTPixelClusterShapeFilter()
Power< A, B >::type pow(const A &a, const B &b)