161 std::vector<edm::Handle<reco::TrackToTrackingParticleAssociator>>
associators;
172 std::unordered_set<unsigned> PFTrackToGeneralTrack;
175 PFTrackToGeneralTrack.insert(ptr->trackRef().key());
185 std::unordered_set<unsigned> MuonTrackToGeneralTrack;
186 for (
auto const&
mu : *
muons.product()) {
189 MuonTrackToGeneralTrack.insert(muTrkRef.
key());
223 const std::vector<reco::PFCluster>& SimClusters = *SimClustersH;
225 std::unordered_map<uint64_t, size_t> hashToSimCluster;
227 for (
unsigned i = 0;
i < SimClustersTruth.size(); ++
i) {
228 const auto& simTruth = SimClustersTruth[
i];
229 hashToSimCluster[hashSimInfo(simTruth)] =
i;
233 std::vector<reco::RecoToSimCollection> associatedTracks, associatedTracksGsf;
235 associatedTracks.emplace_back(
associator->associateRecoToSim(TrackCollectionH, TPCollectionH));
241 auto superclusters = std::make_unique<reco::SuperClusterCollection>();
242 auto blocks = std::make_unique<reco::PFBlockCollection>();
243 std::unordered_map<size_t, size_t> simCluster2Block;
244 std::unordered_map<size_t, size_t> simCluster2BlockIndex;
245 std::unordered_multimap<size_t, size_t> caloParticle2SimCluster;
246 std::vector<int> caloParticle2SuperCluster;
247 for (
unsigned icp = 0; icp < CaloParticles.size(); ++icp) {
250 const auto&
simclusters = CaloParticles[icp].simClusters();
253 std::vector<size_t> good_simclusters;
254 for (
unsigned isc = 0; isc <
simclusters.size(); ++isc) {
260 good_simclusters.push_back(isc);
261 etot += clusterRef->energy();
262 pttot += clusterRef->pt();
264 block.addElement(bec.get());
265 simCluster2Block[simc.key()] = icp;
266 simCluster2BlockIndex[simc.key()] = bec->index();
267 caloParticle2SimCluster.emplace(CaloParticles[icp].g4Tracks()[0].trackId(), simc.key());
273 caloParticle2SuperCluster.push_back(-1);
275 caloParticle2SuperCluster[icp] = superclusters->size();
280 for (
auto idx : good_simclusters) {
283 if (
seed.isNull() ||
seed->energy() < ptr->energy()) {
285 seedpos = ptr->position();
288 superclusters->emplace_back(etot, seedpos,
seed,
clusters);
293 auto superClustersHandle = evt.
put(
std::move(superclusters),
"perfect");
298 usedSimCluster(SimClusters.size(),
false);
300 auto candidates = std::make_unique<reco::PFCandidateCollection>();
305 if (PFTrackToGeneralTrack.count(itk) == 0)
313 if (assoc_tps == associatedTracks.back().end())
319 const auto charge = tkRef->charge();
320 const auto three_mom = tkRef->momentum();
347 candidate.setTime((*trackTimeH)[tkRef], (*trackTimeErrH)[tkRef]);
349 candidate.setTime(0., -1.);
356 if (hashToSimCluster.count(
hash)) {
357 auto simcHash = hashToSimCluster[
hash];
359 if (!usedSimCluster[simcHash]) {
360 if (simCluster2Block.count(simcHash) && simCluster2BlockIndex.count(simcHash)) {
361 size_t block = simCluster2Block.find(simcHash)->second;
362 size_t blockIdx = simCluster2BlockIndex.find(simcHash)->second;
364 candidate.addElementInBlock(blockRef,
blockIdx);
365 usedSimCluster[simcHash] =
true;
368 if (absPdgId == 11) {
369 if (simCluster2Block.count(simcHash)) {
370 auto block_index = simCluster2Block.find(simcHash)->second;
371 auto supercluster_index = caloParticle2SuperCluster[block_index];
372 if (supercluster_index != -1) {
374 for (
const auto& elem : blockRef->elements()) {
375 const auto& ref = elem.clusterRef();
376 if (!usedSimCluster[ref.key()]) {
377 candidate.addElementInBlock(blockRef, elem.index());
378 usedSimCluster[ref.key()] =
true;
392 if (caloParticle2SimCluster.count(
match.first->g4Tracks()[0].trackId())) {
393 auto range = caloParticle2SimCluster.equal_range(
match.first->g4Tracks()[0].trackId());
395 if (!usedSimCluster[
it->second]) {
396 usedSimCluster[
it->second] =
true;
397 if (simCluster2Block.find(
it->second) != simCluster2Block.end()) {
398 size_t block = simCluster2Block.find(
it->second)->second;
399 size_t blockIdx = simCluster2BlockIndex.find(
it->second)->second;
401 candidate.addElementInBlock(blockRef,
blockIdx);
407 usedTrack[tkRef.key()] =
true;
409 if (MuonTrackToGeneralTrack.count(itk) || absPdgId == 13)
415 const auto& theblocks = *blocksHandle;
416 for (
unsigned ibl = 0; ibl < theblocks.size(); ++ibl) {
418 const auto&
elements = theblocks[ibl].elements();
420 const auto& ref = elem.clusterRef();
421 const auto& simtruth = SimClustersTruth[ref.key()];
423 if (!usedSimCluster[ref.key()]) {
424 auto absPdgId =
std::abs(simtruth.pdgId());
433 const auto three_mom = (ref->position() -
math::XYZPoint(0, 0, 0)).
unit() * ref->correctedEnergy();
435 candidates->emplace_back(0, clu_p4, part_type);
437 candidate.addElementInBlock(blref, elem.index());
const edm::EDGetTokenT< edm::ValueMap< float > > srcTrackTimeError_
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
ParticleType
particle types
std::vector< Track > TrackCollection
collection of Tracks
bool isNonnull() const
Checks for non-null.
bool getByToken(EDGetToken token, Handle< PROD > &result) const
muons
the two sets of parameters below are mutually exclusive, depending if RECO or ALCARECO is used the us...
key_type key() const
Accessor for product key.
const double neutralHADThreshold_
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
const edm::EDGetTokenT< edm::ValueMap< float > > srcTrackTimeQuality_
const double timingQualityThreshold_
const edm::EDGetTokenT< CaloParticleCollection > caloParticles_
std::tuple< layerClusterToCaloParticle, caloParticleToLayerCluster > association
const std::vector< edm::EDGetTokenT< reco::TrackToTrackingParticleAssociator > > associators_
Abs< T >::type abs(const T &t)
key
prepare the HTCondor submission files and eventually submit them
const edm::EDGetTokenT< edm::View< reco::Track > > tracks_
const edm::EDGetTokenT< edm::ValueMap< float > > srcGsfTrackTime_
Basic3DVector unit() const
const edm::EDGetTokenT< edm::View< reco::PFRecTrack > > pfRecTracks_
const edm::EDGetTokenT< SimClusterCollection > simClustersTruth_
const edm::EDGetTokenT< TrackingParticleCollection > trackingParticles_
unsigned long long uint64_t
const bool useTimingQuality_
XYZPointD XYZPoint
point in space with cartesian internal representation
const double superClusterThreshold_
const edm::EDGetTokenT< reco::MuonCollection > muons_
const edm::EDGetTokenT< edm::ValueMap< float > > srcGsfTrackTimeQuality_
std::vector< l1t::PFTrack > PFTrackCollection
const edm::EDGetTokenT< edm::ValueMap< float > > srcTrackTime_
std::vector< CaloParticle > CaloParticleCollection
const edm::EDGetTokenT< edm::ValueMap< float > > srcGsfTrackTimeError_
std::pair< typename Association::data_type::first_type, double > match(Reference key, Association association, bool bestMatchByMaxValue)
Generic matching function.
const double neutralEMThreshold_
std::vector< SimCluster > SimClusterCollection
const edm::EDGetTokenT< std::vector< reco::PFCluster > > simClusters_
Static Public Member Functions inherited from edm::global::EDProducerBase
1.8.14