/data/refman/pasoursint/CMSSW_4_2_9_HLT1_bphpatch4/src/RecoTracker/TrackProducer/src/ClusterRemovalRefSetter.cc

Go to the documentation of this file.

#include "RecoTracker/TrackProducer/interface/ClusterRemovalRefSetter.h"
#include "DataFormats/SiPixelDetId/interface/PixelSubdetector.h"

ClusterRemovalRefSetter::ClusterRemovalRefSetter(const edm::Event &iEvent, const edm::InputTag tag) {
    edm::Handle<reco::ClusterRemovalInfo> hCRI;
    iEvent.getByLabel(tag, hCRI);
    cri_ = &*hCRI; 

    //std::cout << "Rekeying PIXEL ProdID " << cri_->pixelNewRefProd().id() << " => " << cri_->pixelRefProd().id() << std::endl;
    //std::cout << "Rekeying STRIP ProdID " << cri_->stripNewRefProd().id() << " => " << cri_->stripRefProd().id() << std::endl;
}

void ClusterRemovalRefSetter::reKey(TrackingRecHit *hit) const {
    if (!hit->isValid()) return;
    DetId detid = hit->geographicalId(); 
    uint32_t subdet = detid.subdetId();
    if ((subdet == PixelSubdetector::PixelBarrel) || (subdet == PixelSubdetector::PixelEndcap)) {
        if (!cri_->hasPixel()) return;
        reKey(reinterpret_cast<SiPixelRecHit *>(hit), detid.rawId());
    } else {
        if (!cri_->hasStrip()) return;
        const std::type_info &type = typeid(*hit);
        if (type == typeid(SiStripRecHit2D)) {
            reKey(reinterpret_cast<SiStripRecHit2D *>(hit), detid.rawId());
        } else if (type == typeid(SiStripRecHit1D)) {
            reKey(reinterpret_cast<SiStripRecHit1D *>(hit), detid.rawId());
        } else if (type == typeid(SiStripMatchedRecHit2D)) {
            SiStripMatchedRecHit2D *mhit = reinterpret_cast<SiStripMatchedRecHit2D *>(hit);
            // const_cast is needed: monoHit() e stereoHit() are const only - at least for now
            reKey(mhit->monoHit(), mhit->monoHit()->geographicalId().rawId());
            reKey(mhit->stereoHit(), mhit->stereoHit()->geographicalId().rawId());
        } else if (type == typeid(ProjectedSiStripRecHit2D)) {
            ProjectedSiStripRecHit2D *phit = reinterpret_cast<ProjectedSiStripRecHit2D *>(hit);
            reKey(&phit->originalHit(), phit->originalHit().geographicalId().rawId());
        } else throw cms::Exception("Unknown RecHit Type") << "RecHit of type " << type.name() << " not supported. (use c++filt to demangle the name)";
    }
}

void ClusterRemovalRefSetter::reKey(SiStripRecHit2D *hit, uint32_t detid) const {
    using reco::ClusterRemovalInfo;
    const ClusterRemovalInfo::Indices &indices = cri_->stripIndices();
    SiStripRecHit2D::ClusterRef newRef = hit->cluster();
    // "newRef" as it refs to the "new"(=cleaned) collection, instead of the old one

    if (newRef.id() != cri_->stripNewRefProd().id()) {   // this is a cfg error in the tracking configuration, much more likely
        throw cms::Exception("Inconsistent Data") << "ClusterRemovalRefSetter: " << 
            "Existing strip cluster refers to product ID " << newRef.id() << 
            " while the ClusterRemovalInfo expects as *new* cluster collection the ID " << cri_->stripNewRefProd().id() << "\n";
    }

    size_t newIndex = newRef.key();
    assert(newIndex < indices.size());
    size_t oldIndex = indices[newIndex];
    SiStripRecHit2D::ClusterRef oldRef(cri_->stripRefProd(), oldIndex);
    hit->setClusterRef(oldRef);
}

void ClusterRemovalRefSetter::reKey(SiStripRecHit1D *hit, uint32_t detid) const {
    using reco::ClusterRemovalInfo;
    const ClusterRemovalInfo::Indices &indices = cri_->stripIndices();
    SiStripRecHit1D::ClusterRef newRef = hit->cluster();
    // "newRef" as it refs to the "new"(=cleaned) collection, instead of the old one

    if (newRef.id() != cri_->stripNewRefProd().id()) {   // this is a cfg error in the tracking configuration, much more likely
        throw cms::Exception("Inconsistent Data") << "ClusterRemovalRefSetter: " << 
            "Existing strip cluster refers to product ID " << newRef.id() << 
            " while the ClusterRemovalInfo expects as *new* cluster collection the ID " << cri_->stripNewRefProd().id() << "\n";
    }

    size_t newIndex = newRef.key();
    assert(newIndex < indices.size());
    size_t oldIndex = indices[newIndex];
    SiStripRecHit1D::ClusterRef oldRef(cri_->stripRefProd(), oldIndex);
    hit->setClusterRef(oldRef);
}

void ClusterRemovalRefSetter::reKey(SiPixelRecHit *hit, uint32_t detid) const {
    using reco::ClusterRemovalInfo;
    const ClusterRemovalInfo::Indices &indices = cri_->pixelIndices();
    SiPixelRecHit::ClusterRef newRef = hit->cluster();
    // "newRef" as it refs to the "new"(=cleaned) collection, instead of the old one

    if (newRef.id()  != cri_->pixelNewRefProd().id()) {
        throw cms::Exception("Inconsistent Data") << "ClusterRemovalRefSetter: " << 
            "Existing pixel cluster refers to product ID " << newRef.id() << 
            " while the ClusterRemovalInfo expects as *new* cluster collection the ID " << cri_->pixelNewRefProd().id() << "\n";
    }
    size_t newIndex = newRef.key();
    assert(newIndex < indices.size());
    size_t oldIndex = indices[newIndex];
    SiPixelRecHit::ClusterRef oldRef(cri_->pixelRefProd(), oldIndex);
    hit->setClusterRef(oldRef);
}