diff --git a/Reconstruction/SiliconTracking/src/ForwardTrackingAlg.cpp b/Reconstruction/SiliconTracking/src/ForwardTrackingAlg.cpp index 9badc084e5c1e1be9c660286b63329c79689f914..6833e4041883b5bfc42ec21a9ac053bf4b5b8d70 100644 --- a/Reconstruction/SiliconTracking/src/ForwardTrackingAlg.cpp +++ b/Reconstruction/SiliconTracking/src/ForwardTrackingAlg.cpp @@ -216,9 +216,10 @@ StatusCode ForwardTrackingAlg::execute(){ debug() << "\t\t---Reading in Collections---" << endmsg; Navigation::Instance()->Initialize(); std::vector<const edm4hep::TrackerHitCollection*> hitFTDCollections; - + int pixelCollectionID = -1; try { auto hitFTDPixelCol = _inFTDPixelColHdl.get(); + pixelCollectionID = hitFTDPixelCol->getID(); hitFTDCollections.push_back(hitFTDPixelCol); Navigation::Instance()->AddTrackerHitCollection(hitFTDPixelCol); } @@ -262,6 +263,9 @@ StatusCode ForwardTrackingAlg::execute(){ debug() << "Number of hits in collection " << hitFTDCollections[iCol]->getID() << ": " << nHits << endmsg; for(auto trackerHit : *hitFTDCollections[iCol]){ + if(pixelCollectionID==hitFTDCollections[iCol]->getID()){ + if ( UTIL::BitSet32( trackerHit.getType() )[ UTIL::ILDTrkHitTypeBit::ONE_DIMENSIONAL ] ) continue; + } edm4hep::ConstTrackerHit hit = trackerHit; debug() << "hit " << trackerHit.id() << " " << KiTrackMarlin::getCellID0Info( trackerHit.getCellID() ) << " " << KiTrackMarlin::getPositionInfo( hit )<< endmsg; diff --git a/Reconstruction/SiliconTracking/src/SiliconTrackingAlg.cpp b/Reconstruction/SiliconTracking/src/SiliconTrackingAlg.cpp index 2e7ae95c7ce6dc0752f646d9feb8b086d839ee4a..9dc2964076fb4020e1a99220e9d16839450b7e11 100644 --- a/Reconstruction/SiliconTracking/src/SiliconTrackingAlg.cpp +++ b/Reconstruction/SiliconTracking/src/SiliconTrackingAlg.cpp @@ -435,7 +435,8 @@ int SiliconTrackingAlg::InitialiseFTD() { _nTotalFTDHits = nelem; //for (int ielem=0; ielem<nelem; ++ielem) { - for(auto hit : *hitFTDPixelCol){ + for(auto hit : *hitFTDPixelCol){ + if ( UTIL::BitSet32( hit.getType() )[ UTIL::ILDTrkHitTypeBit::ONE_DIMENSIONAL ] ) continue; //dm4hep::ConstTrackerHit hit = hitFTDPixelCol->at(ielem); TrackerHitExtended * hitExt = new TrackerHitExtended( hit ); //gear::Vector3D U(1.0,hit->getU()[1],hit->getU()[0],gear::Vector3D::spherical);