if it is a cluster to project
Definition at line 222 of file FastPrimaryVertexWithWeightsProducer.cc.
224 using namespace reco;
227 const float barrel_lenght = 30;
263 std::vector<float> zProjections;
264 std::vector<float> zWeights;
268 float px = (*jit)->px();
269 float py = (*jit)->py();
270 float pz = (*jit)->pz();
271 float pt = (*jit)->pt();
272 float eta = (*jit)->eta();
273 float jetZOverRho = (*jit)->momentum().Z() / (*jit)->momentum().Rho();
278 DetId id = it->detId();
281 float zmodule = modulepos.
z() -
285 for (
size_t j = 0;
j < detset.
size();
j++) {
298 pp->localParametersV(aCluster, (*trackerGeometry->
idToDetUnit(
id)))[0].first);
309 zProjections.push_back(
z);
312 if (
std::abs(modulepos.
z()) < barrel_lenght) {
318 float weight_sizeY_up = (sizeY_up -
sizeY) / (sizeY_up - sizeY_peak);
319 float weight_sizeY_down = (
sizeY - sizeY_down) / (sizeY_peak - sizeY_down);
320 weight_sizeY_down = weight_sizeY_down * (weight_sizeY_down > 0) * (weight_sizeY_down < 1);
321 weight_sizeY_up = weight_sizeY_up * (weight_sizeY_up > 0) * (weight_sizeY_up < 1);
322 float weight_sizeY = weight_sizeY_up + weight_sizeY_down;
325 float rho =
sqrt(v_bs.x() * v_bs.x() + v_bs.y() * v_bs.y());
337 float weightCluster_down =
339 weightCluster_down = weightCluster_down * (weightCluster_down > 0) * (weightCluster_down < 1);
340 weightCluster_up = weightCluster_up * (weightCluster_up > 0) * (weightCluster_up < 1);
341 float weight_charge = weightCluster_up + weightCluster_down;
344 weight =
weight_dPhi * weight_sizeY * weight_rho * weight_sizeX1 * weight_charge;
345 }
else if (
std::abs(modulepos.
z()) > barrel_lenght)
354 zWeights.push_back(
weight);
365 std::multimap<float, float> zWithW;
367 for (
i = 0;
i < zProjections.size();
i++)
368 zWithW.insert(std::pair<float, float>(zProjections[
i], zWeights[
i]));
370 for (std::multimap<float, float>::iterator it = zWithW.begin(); it != zWithW.end(); it++,
i++) {
371 zProjections[
i] = it->first;
372 zWeights[
i] = it->second;
376 std::vector<float> zWeightsSquared;
377 for (std::vector<float>::iterator it = zWeights.begin(); it != zWeights.end(); it++) {
378 zWeightsSquared.push_back((*it) * (*it));
392 float centerWMax = res_step3;
397 if (zProjections.size() > 2) {
400 e(0, 0) = 0.0015 * 0.0015;
401 e(1, 1) = 0.0015 * 0.0015;
405 auto pOut = std::make_unique<reco::VertexCollection>();
406 pOut->push_back(thePV);
410 e(0, 0) = 0.0015 * 0.0015;
411 e(1, 1) = 0.0015 * 0.0015;
415 auto pOut = std::make_unique<reco::VertexCollection>();
416 pOut->push_back(thePV);
422 const float half_width_peak = 1;
423 float nWeightedTot = 0;
424 float nWeightedTotPeak = 0;
425 for (std::vector<float>::iterator it = zProjections.begin(); it != zProjections.end(); it++) {
426 nWeightedTot += zWeights[it - zProjections.begin()];
427 if ((
res - half_width_peak) <= (*it) && (*it) <= (
res + half_width_peak)) {
428 nWeightedTotPeak += zWeights[it - zProjections.begin()];
432 auto zClusterQuality = std::make_unique<float>();
433 *zClusterQuality = -1;
434 if (nWeightedTot != 0) {
436 *zClusterQuality = nWeightedTotPeak /
sqrt(nWeightedTot / (2 * half_width_peak));
References funct::abs(), pwdgSkimBPark_cfi::beamSpot, beamSpotToken, ALCARECOTkAlJpsiMuMu_cff::charge, SiPixelCluster::charge(), clustersToken, SiPixelRawToDigiRegional_cfi::deltaPhi, MillePedeFileConverter_cfg::e, PVValHelper::eta, JetChargeProducer_cfi::exp, FindPeakFastPV(), edm::EventSetup::getData(), mps_fire::i, TrackerGeometry::idToDet(), TrackerGeometry::idToDetUnit(), iEvent, dqmiolumiharvest::j, singleTopDQM_cfi::jets, jetsToken, m_barrel, m_EC_weight, m_endCap, m_geomToken, m_maxDeltaPhi, m_maxDeltaPhi_EC, m_maxJetEta, m_maxJetEta_EC, m_maxSizeX, m_maxSizeY_q, m_maxZ, m_minJetEta_EC, m_minJetPt, m_minSizeY_q, m_njets, m_peakSizeY_q, m_PixelCellHeightOverWidth, m_pixelCPEToken, m_ptWeighting, m_ptWeighting_offset, m_ptWeighting_slope, m_weight_charge_down, m_weight_charge_peak, m_weight_charge_up, m_weight_dPhi, m_weight_dPhi_EC, m_weight_rho_up, m_weight_SizeX1, m_weightCut_step2, m_weightCut_step3, m_zClusterSearchArea_step2, m_zClusterSearchArea_step3, m_zClusterWidth_step1, m_zClusterWidth_step2, m_zClusterWidth_step3, eostools::move(), AlCaHLTBitMon_ParallelJobs::p, PV3DBase< T, PVType, FrameType >::phi(), LumiMonitor_cff::pixelClusters, GeomDet::position(), L1TrackObjectNtupleMaker_cfg::pOut, createTree::pp, edm::Handle< T >::product(), DiDispStaMuonMonitor_cfi::pt, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, rho, objectSelection_cff::selectedJets, edmNew::DetSet< T >::size(), SiPixelCluster::sizeX(), SiPixelCluster::sizeY(), RecoTauValidation_cfi::sizeY, mathSSE::sqrt(), GeomDet::surface(), Surface::toGlobal(), findQualityFiles::v, HLT_FULL_cff::weight_dPhi, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), z, and PV3DBase< T, PVType, FrameType >::z().
Public Member Functions inherited from edm::stream::EDProducer<>
1.8.16