Inheritance diagram for OverlapValidation:

Public Member Functions
	OverlapValidation (const edm::ParameterSet &)

	~OverlapValidation () override

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Types
typedef TransientTrackingRecHit::ConstRecHitPointer	ConstRecHitPointer

typedef vector< Trajectory >	TrajectoryCollection

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

virtual void	analyzeTrajectory (const Trajectory &, const Propagator &, TrackerHitAssociator &, const TrackerTopology *const tTopo)

void	endJob () override

int	layerFromId (const DetId &, const TrackerTopology *const tTopo) const

Private Attributes
vector< bool >	acceptLayer

const bool	addExtraBranches_

bool	barrelOnly_

float	chi2_ [2]

const float	chi2ProbCut_

uint	clusterCharge_ [2]

float	clusterSize_ [2]

float	clusterWidthX_ [2]

float	clusterWidthY_ [2]

TrajectoryStateCombiner	combiner_

const int	compressionSettings_

edm::ParameterSet	config_

SiStripDetInfo	detInfo_

bool	doSimHit_

int	edge_ [2]

uint	event_

edm::FileInPath	FileInPath_

const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	geomToken_

unsigned short	hitCounts_ [2]

float	hitErrors_ [2]

float	hitErrorsY_ [2]

float	hitPositions_ [2]

float	hitPositionsY_ [2]

uint	layer_

float	localxdotglobalphi_ [2]

float	localxdotglobalr_ [2]

float	localxdotglobalx_ [2]

float	localxdotglobaly_ [2]

float	localxdotglobalz_ [2]

float	localydotglobalphi_ [2]

float	localydotglobalr_ [2]

float	localydotglobalx_ [2]

float	localydotglobaly_ [2]

float	localydotglobalz_ [2]

const MagneticField *	magField_

const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	magFieldToken_

const int	minHitsCut_

float	moduleX_ [2]

float	moduleY_ [2]

float	moduleZ_ [2]

float	momentum_

const Point2DBase< float, LocalTag >	onezero = Point2DBase<float, LocalTag>(1, 0)

int	overlapCounts_ [3]

unsigned int	overlapIds_ [2]

float	overlapPath_

float	predictedDeltaXError_

float	predictedDeltaYError_

float	predictedLocalErrors_ [5][2]

float	predictedLocalParameters_ [5][2]

float	predictedPositions_ [3][2]

char	relativeXSign_

char	relativeYSign_

TTree *	rootTree_

uint	run_

float	simHitPositions_ [2]

float	simHitPositionsY_ [2]

int	subdetID

const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	topoToken_

const TrackerGeometry *	trackerGeometry_

edm::InputTag	trajectoryTag_

edm::EDGetTokenT< TrajectoryCollection >	trajectoryToken_

const Point2DBase< float, LocalTag >	zeroone = Point2DBase<float, LocalTag>(0, 1)

const Point2DBase< float, LocalTag >	zerozero = Point2DBase<float, LocalTag>(0, 0)

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 82 of file OverlapValidation.cc.

Member Typedef Documentation

◆ ConstRecHitPointer

typedef TransientTrackingRecHit::ConstRecHitPointer OverlapValidation::ConstRecHitPointer

private

Definition at line 92 of file OverlapValidation.cc.

◆ TrajectoryCollection

typedef vector<Trajectory> OverlapValidation::TrajectoryCollection

private

Definition at line 89 of file OverlapValidation.cc.

Constructor & Destructor Documentation

◆ OverlapValidation()

OverlapValidation::OverlapValidation ( const edm::ParameterSet & iConfig )

explicit

Definition at line 187 of file OverlapValidation.cc.

References acceptLayer, addExtraBranches_, barrelOnly_, chi2_, clusterCharge_, clusterSize_, clusterWidthX_, clusterWidthY_, compressionSettings_, edm::ConsumesCollector::consumes(), edm::EDConsumerBase::consumesCollector(), detInfo_, doSimHit_, edge_, event_, FileInPath_, compareTotals::fs, edm::FileInPath::fullPath(), edm::ParameterSet::getParameter(), hitCounts_, hitErrors_, hitErrorsY_, hitPositions_, hitPositionsY_, TFileService::kSharedResource, layer_, localxdotglobalphi_, localxdotglobalr_, localxdotglobalx_, localxdotglobaly_, localxdotglobalz_, localydotglobalphi_, localydotglobalr_, localydotglobalx_, localydotglobaly_, localydotglobalz_, moduleX_, moduleY_, moduleZ_, momentum_, overlapCounts_, overlapIds_, overlapPath_, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, predictedDeltaXError_, predictedDeltaYError_, predictedLocalErrors_, predictedLocalParameters_, predictedPositions_, SiStripDetInfoFileReader::read(), relativeXSign_, relativeYSign_, rootTree_, run_, simHitPositions_, simHitPositionsY_, subdetID, StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, trajectoryTag_, and trajectoryToken_.

     : geomToken_(esConsumes()),
       magFieldToken_(esConsumes()),
       topoToken_(esConsumes()),
       config_(iConfig),
       rootTree_(nullptr),
       FileInPath_(SiStripDetInfoFileReader::kDefaultFile),
       compressionSettings_(iConfig.getUntrackedParameter<int>("compressionSettings", -1)),
       addExtraBranches_(false),
       minHitsCut_(6),
       chi2ProbCut_(0.001) {
   usesResource(TFileService::kSharedResource);
   edm::ConsumesCollector&& iC = consumesCollector();
   //now do what ever initialization is needed
   trajectoryTag_ = iConfig.getParameter<edm::InputTag>("trajectories");
   trajectoryToken_ = iC.consumes<TrajectoryCollection>(trajectoryTag_);
   doSimHit_ = iConfig.getParameter<bool>("associateStrip");
   detInfo_ = SiStripDetInfoFileReader::read(FileInPath_.fullPath());
 
   overlapCounts_[0] = 0;  // #trajectories
   overlapCounts_[1] = 0;  // #hits
   overlapCounts_[2] = 0;  // #overlap hits
   acceptLayer.resize(7, false);
   acceptLayer[PixelSubdetector::PixelBarrel] = iConfig.getParameter<bool>("usePXB");
   acceptLayer[PixelSubdetector::PixelEndcap] = iConfig.getParameter<bool>("usePXF");
   acceptLayer[StripSubdetector::TIB] = iConfig.getParameter<bool>("useTIB");
   acceptLayer[StripSubdetector::TOB] = iConfig.getParameter<bool>("useTOB");
   acceptLayer[StripSubdetector::TID] = iConfig.getParameter<bool>("useTID");
   acceptLayer[StripSubdetector::TEC] = iConfig.getParameter<bool>("useTEC");
   barrelOnly_ = iConfig.getParameter<bool>("barrelOnly");
 
   edm::Service<TFileService> fs;
   //
   // root output
   //
   if (compressionSettings_ > 0) {
     fs->file().SetCompressionSettings(compressionSettings_);
   }
 
   rootTree_ = fs->make<TTree>("Overlaps", "Overlaps");
   if (addExtraBranches_) {
     rootTree_->Branch("hitCounts", hitCounts_, "found/s:lost/s");
     rootTree_->Branch("chi2", chi2_, "chi2/F:ndf/F");
     rootTree_->Branch("path", &overlapPath_, "path/F");
   }
   rootTree_->Branch("layer", &layer_, "layer/i");
   rootTree_->Branch("detids", overlapIds_, "id[2]/i");
   rootTree_->Branch("predPos", predictedPositions_, "gX[2]/F:gY[2]/F:gZ[2]/F");
   rootTree_->Branch("predPar", predictedLocalParameters_, "predQP[2]/F:predDX[2]/F:predDY[2]/F:predX[2]/F:predY[2]/F");
   rootTree_->Branch("predErr", predictedLocalErrors_, "predEQP[2]/F:predEDX[2]/F:predEDY[2]/F:predEX[2]/F:predEY[2]/F");
   rootTree_->Branch("predEDeltaX", &predictedDeltaXError_, "sigDeltaX/F");
   rootTree_->Branch("predEDeltaY", &predictedDeltaYError_, "sigDeltaY/F");
   rootTree_->Branch("relSignX", &relativeXSign_, "relSignX/B");
   rootTree_->Branch("relSignY", &relativeYSign_, "relSignY/B");
   rootTree_->Branch("hitX", hitPositions_, "hitX[2]/F");
   rootTree_->Branch("hitEX", hitErrors_, "hitEX[2]/F");
   rootTree_->Branch("hitY", hitPositionsY_, "hitY[2]/F");
   rootTree_->Branch("hitEY", hitErrorsY_, "hitEY[2]/F");
   if (addExtraBranches_) {
     rootTree_->Branch("simX", simHitPositions_, "simX[2]/F");
     rootTree_->Branch("simY", simHitPositionsY_, "simY[2]/F");
     rootTree_->Branch("clusterSize", clusterSize_, "clusterSize[2]/F");
     rootTree_->Branch("clusterWidthX", clusterWidthX_, "clusterWidthX[2]/F");
     rootTree_->Branch("clusterWidthY", clusterWidthY_, "clusterWidthY[2]/F");
     rootTree_->Branch("clusterCharge", clusterCharge_, "clusterCharge[2]/i");
     rootTree_->Branch("edge", edge_, "edge[2]/I");
   }
   rootTree_->Branch("momentum", &momentum_, "momentum/F");
   rootTree_->Branch("run", &run_, "run/i");
   rootTree_->Branch("event", &event_, "event/i");
   rootTree_->Branch("subdetID", &subdetID, "subdetID/I");
   rootTree_->Branch("moduleX", moduleX_, "moduleX[2]/F");
   rootTree_->Branch("moduleY", moduleY_, "moduleY[2]/F");
   rootTree_->Branch("moduleZ", moduleZ_, "moduleZ[2]/F");
   rootTree_->Branch("localxdotglobalphi", localxdotglobalphi_, "localxdotglobalphi[2]/F");
   rootTree_->Branch("localxdotglobalr", localxdotglobalr_, "localxdotglobalr[2]/F");
   rootTree_->Branch("localxdotglobalz", localxdotglobalz_, "localxdotglobalz[2]/F");
   rootTree_->Branch("localxdotglobalx", localxdotglobalx_, "localxdotglobalx[2]/F");
   rootTree_->Branch("localxdotglobaly", localxdotglobaly_, "localxdotglobaly[2]/F");
   rootTree_->Branch("localydotglobalphi", localydotglobalphi_, "localydotglobalphi[2]/F");
   rootTree_->Branch("localydotglobalr", localydotglobalr_, "localydotglobalr[2]/F");
   rootTree_->Branch("localydotglobalz", localydotglobalz_, "localydotglobalz[2]/F");
   rootTree_->Branch("localydotglobalx", localydotglobalx_, "localydotglobalx[2]/F");
   rootTree_->Branch("localydotglobaly", localydotglobaly_, "localydotglobaly[2]/F");
 }

◆ ~OverlapValidation()

OverlapValidation::~OverlapValidation ( )

override

Definition at line 291 of file OverlapValidation.cc.

References overlapCounts_, and w().

                                       {
   edm::LogWarning w("Overlaps");
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
   w << "Counters =";
   w << " Number of tracks: " << overlapCounts_[0];
   w << " Number of valid hits: " << overlapCounts_[1];
   w << " Number of overlaps: " << overlapCounts_[2];
 }

Member Function Documentation

◆ analyze()

void OverlapValidation::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 306 of file OverlapValidation.cc.

References AnalyticalPropagator_cfi::AnalyticalPropagator, analyzeTrajectory(), anyDirection, LCToCPAssociation_cfi::associator, config_, edm::EDConsumerBase::consumesCollector(), doSimHit_, event_, geomToken_, edm::EventSetup::getData(), iEvent, magField_, magFieldToken_, edm::Handle< T >::product(), TrackCandidateProducer_cfi::propagator, run_, topoToken_, trackerGeometry_, and trajectoryToken_.

                                                                                    {
   using namespace edm;
   //
   // mag field & search tracker
   //
   magField_ = &iSetup.getData(magFieldToken_);
   //
   // propagator
   //
   AnalyticalPropagator propagator(magField_, anyDirection);
   //
   // geometry
   //
   trackerGeometry_ = &iSetup.getData(geomToken_);
   //
   // make associator for SimHits
   //
   TrackerHitAssociator* associator;
   if (doSimHit_) {
     TrackerHitAssociator::Config hitassociatorconfig(config_, consumesCollector());
     associator = new TrackerHitAssociator(iEvent, hitassociatorconfig);
   } else {
     associator = nullptr;
   }
 
   //if(doSimHit_) associator = new TrackerHitAssociator(iEvent, config_); else associator = 0;
 
   //
   // trajectories (from refit)
   //
   //typedef vector<Trajectory> TrajectoryCollection;
   edm::Handle<TrajectoryCollection> trajectoryCollectionHandle;
   iEvent.getByToken(trajectoryToken_, trajectoryCollectionHandle);
   const TrajectoryCollection* const trajectoryCollection = trajectoryCollectionHandle.product();
 
   //
   // loop over trajectories from refit
   const TrackerTopology* const tTopo = &iSetup.getData(topoToken_);
   for (const auto& trajectory : *trajectoryCollection)
     analyzeTrajectory(trajectory, propagator, *associator, tTopo);
 
   run_ = iEvent.id().run();
   event_ = iEvent.id().event();
 }

◆ analyzeTrajectory()

void OverlapValidation::analyzeTrajectory	(	const Trajectory &	trajectory,
		const Propagator &	propagator,
		TrackerHitAssociator &	associator,
		const TrackerTopology *const	tTopo
	)

privatevirtual

Definition at line 351 of file OverlapValidation.cc.

References funct::abs(), acceptLayer, LCToCPAssociation_cfi::associator, barrelOnly_, ALCARECOTkAlJpsiMuMu_cff::charge, chi2_, chi2ProbCut_, Trajectory::chiSquared(), ChiSquaredProbability(), SiStripRecHit1D::cluster(), SiPixelRecHit::cluster(), clusterCharge_, clusterSize_, clusterWidthX_, clusterWidthY_, TrajectoryStateCombiner::combine(), combiner_, detInfo_, change_name::diff, doSimHit_, edge_, Trajectory::foundHits(), SiStripDetInfo::getNumberOfApvsAndStripLength(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalParameters(), TrajectoryStateOnSurface::globalPosition(), SiStripDetId::glued(), runTauDisplay::gp, hitCounts_, hitErrors_, hitErrorsY_, hitPositions_, hitPositionsY_, mps_fire::i, globals_cff::id1, globals_cff::id2, PixelTopology::isItEdgePixelInX(), PixelTopology::isItEdgePixelInY(), TrajectoryStateOnSurface::isValid(), dqmiolumiharvest::j, JacobianCurvilinearToLocal::jacobian(), nano_mu_digi_cff::layer, layer_, layerFromId(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), localxdotglobalphi_, localxdotglobalr_, localxdotglobalx_, localxdotglobaly_, localxdotglobalz_, localydotglobalphi_, localydotglobalr_, localydotglobalx_, localydotglobaly_, localydotglobalz_, Trajectory::lostHits(), visualization-live-secondInstance_cfg::m, PV3DBase< T, PVType, FrameType >::mag(), magField_, LocalTrajectoryError::matrix(), Trajectory::measurements(), minHitsCut_, moduleX_, moduleY_, moduleZ_, momentum_, GluedGeomDet::monoDet(), Trajectory::ndof(), onezero, HLT_2024v12_cff::Overlap, overlapCounts_, overlapIds_, overlapPath_, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, position, predictedDeltaXError_, predictedDeltaYError_, predictedLocalErrors_, predictedLocalParameters_, predictedPositions_, TrackCandidateProducer_cfi::propagator, nano_mu_digi_cff::rawId, DetId::rawId(), relativeXSign_, relativeYSign_, rootTree_, pfClustersFromCombinedCaloHF_cfi::scale, simHitPositions_, simHitPositionsY_, findQualityFiles::size, PixelGeomDetUnit::specificTopology(), mathSSE::sqrt(), SiStripDetId::stereo(), DetId::subdetId(), subdetID, GeomDet::surface(), TrajectoryStateOnSurface::surface(), StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, Surface::toGlobal(), GloballyPositioned< T >::toLocal(), parallelization::uint, LocalTrajectoryParameters::vector(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), PV3DBase< T, PVType, FrameType >::z(), zeroone, and zerozero.

Referenced by analyze().

                                                                               {
   typedef std::pair<const TrajectoryMeasurement*, const TrajectoryMeasurement*> Overlap;
   typedef vector<Overlap> OverlapContainer;
   ++overlapCounts_[0];
 
   OverlapContainer overlapHits;
 
   // quality cuts on trajectory
   // min. # hits / matched hits
 
   if (trajectory.foundHits() < minHitsCut_)
     return;
   if (ChiSquaredProbability((double)(trajectory.chiSquared()), (double)(trajectory.ndof(false))) < chi2ProbCut_)
     return;
   //
   // loop over measurements in the trajectory and calculate residuals
   //
 
   vector<TrajectoryMeasurement> measurements(trajectory.measurements());
   for (vector<TrajectoryMeasurement>::const_iterator itm = measurements.begin(); itm != measurements.end(); ++itm) {
     //
     // skip "invalid" (i.e. missing) hits
     //
     ConstRecHitPointer hit = itm->recHit();
     DetId id = hit->geographicalId();
     int layer(layerFromId(id, tTopo));
     int subDet = id.subdetId();
 
     if (!hit->isValid()) {
       edm::LogVerbatim("OverlapValidation") << "Invalid";
       continue;
     }
     if (barrelOnly_ && (subDet == StripSubdetector::TID || subDet == StripSubdetector::TEC))
       return;
 
     //edm::LogVerbatim("OverlapValidation")  << "Check " <<subDet << ", layer = " << layer<<" stereo: "<< ((subDet > 2)?(SiStripDetId(id).stereo()):2);
     //cout << "Check SubID " <<subDet << ", layer = " << layer<<" stereo: "<< ((subDet > 2)?(SiStripDetId(id).stereo()):2) << endl;
 
     //
     // check for overlap: same subdet-id && layer number for
     // two consecutive hits
     //
     ++overlapCounts_[1];
     if ((layer != -1) && (acceptLayer[subDet])) {
       for (vector<TrajectoryMeasurement>::const_iterator itmCompare = itm - 1;
            itmCompare >= measurements.begin() && itmCompare > itm - 4;
            --itmCompare) {
         DetId compareId = itmCompare->recHit()->geographicalId();
 
         if (subDet != compareId.subdetId() || layer != layerFromId(compareId, tTopo))
           break;
         if (!itmCompare->recHit()->isValid())
           continue;
         if ((subDet == PixelSubdetector::PixelBarrel || subDet == PixelSubdetector::PixelEndcap) ||
             (SiStripDetId(id).stereo() == SiStripDetId(compareId).stereo())) {
           overlapHits.push_back(std::make_pair(&(*itmCompare), &(*itm)));
           //edm::LogVerbatim("OverlapValidation") << "adding pair "<< ((subDet > 2)?(SiStripDetId(id).stereo()) : 2)
           //                 << " from layer = " << layer;
           //cout << "adding pair "<< ((subDet > 2)?(SiStripDetId(id).stereo()) : 2) << " from subDet = " << subDet << " and layer = " << layer;
           //cout << " \t"<<run_<< "\t"<<event_<<"\t";
           //cout << min(id.rawId(),compareId.rawId())<<"\t"<<max(id.rawId(),compareId.rawId())<<endl;
           if (SiStripDetId(id).glued() == id.rawId())
             edm::LogInfo("Overlaps") << "BAD GLUED: Have glued layer with id = " << id.rawId()
                                      << " and glued id = " << SiStripDetId(id).glued()
                                      << "  and stereo = " << SiStripDetId(id).stereo() << endl;
           if (SiStripDetId(compareId).glued() == compareId.rawId())
             edm::LogInfo("Overlaps") << "BAD GLUED: Have glued layer with id = " << compareId.rawId()
                                      << " and glued id = " << SiStripDetId(compareId).glued()
                                      << "  and stereo = " << SiStripDetId(compareId).stereo() << endl;
           break;
         }
       }
     }
   }
 
   //
   // Loop over all overlap pairs.
   //
   hitCounts_[0] = trajectory.foundHits();
   hitCounts_[1] = trajectory.lostHits();
   chi2_[0] = trajectory.chiSquared();
   chi2_[1] = trajectory.ndof(false);
 
   for (const auto& overlapHit : overlapHits) {
     //
     // create reference state @ module 1 (no info from overlap hits)
     //
     ++overlapCounts_[2];
     // backward predicted state at module 1
     TrajectoryStateOnSurface bwdPred1 = overlapHit.first->backwardPredictedState();
     if (!bwdPred1.isValid())
       continue;
     //cout << "momentum from backward predicted state = " << bwdPred1.globalMomentum().mag() << endl;
     // forward predicted state at module 2
     TrajectoryStateOnSurface fwdPred2 = overlapHit.second->forwardPredictedState();
     //cout << "momentum from forward predicted state = " << fwdPred2.globalMomentum().mag() << endl;
     if (!fwdPred2.isValid())
       continue;
     // extrapolate fwdPred2 to module 1
     TrajectoryStateOnSurface fwdPred2At1 = propagator.propagate(fwdPred2, bwdPred1.surface());
     if (!fwdPred2At1.isValid())
       continue;
     // combine fwdPred2At1 with bwdPred1 (ref. state, best estimate without hits 1 and 2)
     TrajectoryStateOnSurface comb1 = combiner_.combine(bwdPred1, fwdPred2At1);
     if (!comb1.isValid())
       continue;
     //
     // propagation of reference parameters to module 2
     //
     std::pair<TrajectoryStateOnSurface, double> tsosWithS = propagator.propagateWithPath(comb1, fwdPred2.surface());
     TrajectoryStateOnSurface comb1At2 = tsosWithS.first;
     if (!comb1At2.isValid())
       continue;
     //distance of propagation from one surface to the next==could cut here
     overlapPath_ = tsosWithS.second;
     if (abs(overlapPath_) > 15)
       continue;  //cut to remove hit pairs > 15 cm apart
     // global position on module 1
     GlobalPoint position = comb1.globalPosition();
     predictedPositions_[0][0] = position.x();
     predictedPositions_[1][0] = position.y();
     predictedPositions_[2][0] = position.z();
     momentum_ = comb1.globalMomentum().mag();
     //cout << "momentum from combination = " << momentum_ << endl;
     //cout << "magnetic field from TSOS = " << comb1.magneticField()->inTesla(position).mag() << endl;
     // local parameters and errors on module 1
     AlgebraicVector5 pars = comb1.localParameters().vector();
     AlgebraicSymMatrix55 errs = comb1.localError().matrix();
     for (int i = 0; i < 5; ++i) {
       predictedLocalParameters_[i][0] = pars[i];
       predictedLocalErrors_[i][0] = sqrt(errs(i, i));
     }
     // global position on module 2
     position = comb1At2.globalPosition();
     predictedPositions_[0][1] = position.x();
     predictedPositions_[1][1] = position.y();
     predictedPositions_[2][1] = position.z();
     // local parameters and errors on module 2
     pars = comb1At2.localParameters().vector();
     errs = comb1At2.localError().matrix();
     for (int i = 0; i < 5; ++i) {
       predictedLocalParameters_[i][1] = pars[i];
       predictedLocalErrors_[i][1] = sqrt(errs(i, i));
     }
 
     //print out local errors in X to check
     //cout << "Predicted local error in X at 1 = " << predictedLocalErrors_[3][0] << "   and predicted local error in X at 2 is = " <<  predictedLocalErrors_[3][1] << endl;
     //cout << "Predicted local error in Y at 1 = " << predictedLocalErrors_[4][0] << "   and predicted local error in Y at 2 is = " <<  predictedLocalErrors_[4][1] << endl;
 
     //
     // jacobians (local-to-global@1,global 1-2,global-to-local@2)
     //
     JacobianLocalToCurvilinear jacLocToCurv(comb1.surface(), comb1.localParameters(), *magField_);
     AnalyticalCurvilinearJacobian jacCurvToCurv(
         comb1.globalParameters(), comb1At2.globalPosition(), comb1At2.globalMomentum(), tsosWithS.second);
     JacobianCurvilinearToLocal jacCurvToLoc(comb1At2.surface(), comb1At2.localParameters(), *magField_);
     // combined jacobian local-1-to-local-2
     AlgebraicMatrix55 jacobian = jacLocToCurv.jacobian() * jacCurvToCurv.jacobian() * jacCurvToLoc.jacobian();
     // covariance on module 1
     AlgebraicSymMatrix55 covComb1 = comb1.localError().matrix();
     // variance and correlations for predicted local_x on modules 1 and 2
     double c00 = covComb1(3, 3);
     double c10(0.);
     double c11(0.);
     for (int i = 1; i < 5; ++i) {
       c10 += jacobian(3, i) * covComb1(i, 3);
       for (int j = 1; j < 5; ++j)
         c11 += jacobian(3, i) * covComb1(i, j) * jacobian(3, j);
     }
     // choose relative sign in order to minimize error on difference
     double diff = c00 - 2 * fabs(c10) + c11;
     diff = diff > 0 ? sqrt(diff) : -sqrt(-diff);
     predictedDeltaXError_ = diff;
     relativeXSign_ = c10 > 0 ? -1 : 1;
     //
     // now find variance and correlations for predicted local_y
     double c00Y = covComb1(4, 4);
     double c10Y(0.);
     double c11Y(0.);
     for (int i = 1; i < 5; ++i) {
       c10Y += jacobian(4, i) * covComb1(i, 4);
       for (int j = 1; j < 5; ++j)
         c11Y += jacobian(4, i) * covComb1(i, j) * jacobian(4, j);
     }
     double diffY = c00Y - 2 * fabs(c10Y) + c11Y;
     diffY = diffY > 0 ? sqrt(diffY) : -sqrt(-diffY);
     predictedDeltaYError_ = diffY;
     relativeYSign_ = c10Y > 0 ? -1 : 1;
 
     // information on modules and hits
     overlapIds_[0] = overlapHit.first->recHit()->geographicalId().rawId();
     overlapIds_[1] = overlapHit.second->recHit()->geographicalId().rawId();
 
     //added by Heshy and Jared
     moduleX_[0] = overlapHit.first->recHit()->det()->surface().position().x();
     moduleX_[1] = overlapHit.second->recHit()->det()->surface().position().x();
     moduleY_[0] = overlapHit.first->recHit()->det()->surface().position().y();
     moduleY_[1] = overlapHit.second->recHit()->det()->surface().position().y();
     moduleZ_[0] = overlapHit.first->recHit()->det()->surface().position().z();
     moduleZ_[1] = overlapHit.second->recHit()->det()->surface().position().z();
     subdetID = overlapHit.first->recHit()->geographicalId().subdetId();
     localxdotglobalphi_[0] = overlapHit.first->recHit()->det()->surface().toGlobal(onezero).phi() -
                              overlapHit.first->recHit()->det()->surface().toGlobal(zerozero).phi();
     localxdotglobalphi_[1] = overlapHit.second->recHit()->det()->surface().toGlobal(onezero).phi() -
                              overlapHit.second->recHit()->det()->surface().toGlobal(zerozero).phi();
     //added by Jason
     localxdotglobalr_[0] = overlapHit.first->recHit()->det()->surface().toGlobal(onezero).perp() -
                            overlapHit.first->recHit()->det()->surface().toGlobal(zerozero).perp();
     localxdotglobalr_[1] = overlapHit.second->recHit()->det()->surface().toGlobal(onezero).perp() -
                            overlapHit.second->recHit()->det()->surface().toGlobal(zerozero).perp();
     localxdotglobalz_[0] = overlapHit.first->recHit()->det()->surface().toGlobal(onezero).z() -
                            overlapHit.first->recHit()->det()->surface().toGlobal(zerozero).z();
     localxdotglobalz_[1] = overlapHit.second->recHit()->det()->surface().toGlobal(onezero).z() -
                            overlapHit.second->recHit()->det()->surface().toGlobal(zerozero).z();
     localxdotglobalx_[0] = overlapHit.first->recHit()->det()->surface().toGlobal(onezero).x() -
                            overlapHit.first->recHit()->det()->surface().toGlobal(zerozero).x();
     localxdotglobalx_[1] = overlapHit.second->recHit()->det()->surface().toGlobal(onezero).x() -
                            overlapHit.second->recHit()->det()->surface().toGlobal(zerozero).x();
     localxdotglobaly_[0] = overlapHit.first->recHit()->det()->surface().toGlobal(onezero).y() -
                            overlapHit.first->recHit()->det()->surface().toGlobal(zerozero).y();
     localxdotglobaly_[1] = overlapHit.second->recHit()->det()->surface().toGlobal(onezero).y() -
                            overlapHit.second->recHit()->det()->surface().toGlobal(zerozero).y();
     localydotglobalr_[0] = overlapHit.first->recHit()->det()->surface().toGlobal(zeroone).perp() -
                            overlapHit.first->recHit()->det()->surface().toGlobal(zerozero).perp();
     localydotglobalr_[1] = overlapHit.second->recHit()->det()->surface().toGlobal(zeroone).perp() -
                            overlapHit.second->recHit()->det()->surface().toGlobal(zerozero).perp();
     localydotglobalz_[0] = overlapHit.first->recHit()->det()->surface().toGlobal(zeroone).z() -
                            overlapHit.first->recHit()->det()->surface().toGlobal(zerozero).z();
     localydotglobalz_[1] = overlapHit.second->recHit()->det()->surface().toGlobal(zeroone).z() -
                            overlapHit.second->recHit()->det()->surface().toGlobal(zerozero).z();
     localydotglobalx_[0] = overlapHit.first->recHit()->det()->surface().toGlobal(zeroone).x() -
                            overlapHit.first->recHit()->det()->surface().toGlobal(zerozero).x();
     localydotglobalx_[1] = overlapHit.second->recHit()->det()->surface().toGlobal(zeroone).x() -
                            overlapHit.second->recHit()->det()->surface().toGlobal(zerozero).x();
     localydotglobaly_[0] = overlapHit.first->recHit()->det()->surface().toGlobal(zeroone).y() -
                            overlapHit.first->recHit()->det()->surface().toGlobal(zerozero).y();
     localydotglobaly_[1] = overlapHit.second->recHit()->det()->surface().toGlobal(zeroone).y() -
                            overlapHit.second->recHit()->det()->surface().toGlobal(zerozero).y();
     localydotglobalphi_[0] = overlapHit.first->recHit()->det()->surface().toGlobal(zeroone).phi() -
                              overlapHit.first->recHit()->det()->surface().toGlobal(zerozero).phi();
     localydotglobalphi_[1] = overlapHit.second->recHit()->det()->surface().toGlobal(zeroone).phi() -
                              overlapHit.second->recHit()->det()->surface().toGlobal(zerozero).phi();
 
     if (overlapHit.first->recHit()->geographicalId().subdetId() == StripSubdetector::TIB)
       layer_ = layerFromId(overlapHit.first->recHit()->geographicalId().rawId(), tTopo);
     else if (overlapHit.first->recHit()->geographicalId().subdetId() == StripSubdetector::TOB)
       layer_ = layerFromId(overlapHit.first->recHit()->geographicalId().rawId(), tTopo) + 4;
     else if (overlapHit.first->recHit()->geographicalId().subdetId() == StripSubdetector::TID)
       layer_ = layerFromId(overlapHit.first->recHit()->geographicalId().rawId(), tTopo) + 10;
     else if (overlapHit.first->recHit()->geographicalId().subdetId() == StripSubdetector::TEC)
       layer_ = layerFromId(overlapHit.first->recHit()->geographicalId().rawId(), tTopo) + 13;
     else if (overlapHit.first->recHit()->geographicalId().subdetId() == PixelSubdetector::PixelBarrel)
       layer_ = layerFromId(overlapHit.first->recHit()->geographicalId().rawId(), tTopo) + 20;
     else if (overlapHit.first->recHit()->geographicalId().subdetId() == PixelSubdetector::PixelEndcap)
       layer_ = layerFromId(overlapHit.first->recHit()->geographicalId().rawId(), tTopo) + 30;
     else
       layer_ = 99;
 
     if (overlapIds_[0] == SiStripDetId(overlapHit.first->recHit()->geographicalId()).glued())
       edm::LogWarning("Overlaps") << "BAD GLUED: First Id = " << overlapIds_[0] << " has glued = "
                                   << SiStripDetId(overlapHit.first->recHit()->geographicalId()).glued()
                                   << "  and stereo = "
                                   << SiStripDetId(overlapHit.first->recHit()->geographicalId()).stereo() << endl;
     if (overlapIds_[1] == SiStripDetId(overlapHit.second->recHit()->geographicalId()).glued())
       edm::LogWarning("Overlaps") << "BAD GLUED: Second Id = " << overlapIds_[1] << " has glued = "
                                   << SiStripDetId(overlapHit.second->recHit()->geographicalId()).glued()
                                   << "  and stereo = "
                                   << SiStripDetId(overlapHit.second->recHit()->geographicalId()).stereo() << endl;
 
     const TransientTrackingRecHit::ConstRecHitPointer firstRecHit = overlapHit.first->recHit();
     const TransientTrackingRecHit::ConstRecHitPointer secondRecHit = overlapHit.second->recHit();
     hitPositions_[0] = firstRecHit->localPosition().x();
     hitErrors_[0] = sqrt(firstRecHit->localPositionError().xx());
     hitPositions_[1] = secondRecHit->localPosition().x();
     hitErrors_[1] = sqrt(secondRecHit->localPositionError().xx());
 
     hitPositionsY_[0] = firstRecHit->localPosition().y();
     hitErrorsY_[0] = sqrt(firstRecHit->localPositionError().yy());
     hitPositionsY_[1] = secondRecHit->localPosition().y();
     hitErrorsY_[1] = sqrt(secondRecHit->localPositionError().yy());
 
     //cout << "printing local X hit position and error for the overlap hits. Hit 1 = " << hitPositions_[0] << "+-" << hitErrors_[0] << "  and hit 2 is "  << hitPositions_[1] << "+-" << hitErrors_[1] << endl;
 
     DetId id1 = overlapHit.first->recHit()->geographicalId();
     DetId id2 = overlapHit.second->recHit()->geographicalId();
     int layer1 = layerFromId(id1, tTopo);
     int subDet1 = id1.subdetId();
     int layer2 = layerFromId(id2, tTopo);
     int subDet2 = id2.subdetId();
     if (abs(hitPositions_[0]) > 5)
       edm::LogInfo("Overlaps") << "BAD: Bad hit position: Id = " << id1.rawId()
                                << " stereo = " << SiStripDetId(id1).stereo()
                                << "  glued = " << SiStripDetId(id1).glued() << " from subdet = " << subDet1
                                << " and layer = " << layer1 << endl;
     if (abs(hitPositions_[1]) > 5)
       edm::LogInfo("Overlaps") << "BAD: Bad hit position: Id = " << id2.rawId()
                                << " stereo = " << SiStripDetId(id2).stereo()
                                << "  glued = " << SiStripDetId(id2).glued() << " from subdet = " << subDet2
                                << " and layer = " << layer2 << endl;
 
     // get track momentum
     momentum_ = comb1.globalMomentum().mag();
 
     // get cluster size
     if (!(subDet1 == PixelSubdetector::PixelBarrel || subDet1 == PixelSubdetector::PixelEndcap)) {  //strip
       const TransientTrackingRecHit::ConstRecHitPointer thit1 = overlapHit.first->recHit();
       const SiStripRecHit1D* hit1 = dynamic_cast<const SiStripRecHit1D*>((*thit1).hit());
       if (hit1) {
         //check cluster width
         const SiStripRecHit1D::ClusterRef& cluster1 = hit1->cluster();
         clusterSize_[0] = (cluster1->amplitudes()).size();
         clusterWidthX_[0] = (cluster1->amplitudes()).size();
         clusterWidthY_[0] = -1;
 
         //check if cluster at edge of sensor
         uint16_t firstStrip = cluster1->firstStrip();
         uint16_t lastStrip = firstStrip + (cluster1->amplitudes()).size() - 1;
         unsigned short Nstrips;
         Nstrips = detInfo_.getNumberOfApvsAndStripLength(id1).first * 128;
         bool atEdge = false;
         if (firstStrip == 0 || lastStrip == (Nstrips - 1))
           atEdge = true;
         if (atEdge)
           edge_[0] = 1;
         else
           edge_[0] = -1;
 
         // get cluster total charge
         const auto& stripCharges = cluster1->amplitudes();
         uint16_t charge = 0;
         for (uint i = 0; i < stripCharges.size(); i++) {
           charge += stripCharges[i];
         }
         clusterCharge_[0] = charge;
       }
 
       const TransientTrackingRecHit::ConstRecHitPointer thit2 = overlapHit.second->recHit();
       const SiStripRecHit1D* hit2 = dynamic_cast<const SiStripRecHit1D*>((*thit2).hit());
       if (hit2) {
         const SiStripRecHit1D::ClusterRef& cluster2 = hit2->cluster();
         clusterSize_[1] = (cluster2->amplitudes()).size();
         clusterWidthX_[1] = (cluster2->amplitudes()).size();
         clusterWidthY_[1] = -1;
 
         uint16_t firstStrip = cluster2->firstStrip();
         uint16_t lastStrip = firstStrip + (cluster2->amplitudes()).size() - 1;
         unsigned short Nstrips;
         Nstrips = detInfo_.getNumberOfApvsAndStripLength(id2).first * 128;
         bool atEdge = false;
         if (firstStrip == 0 || lastStrip == (Nstrips - 1))
           atEdge = true;
         if (atEdge)
           edge_[1] = 1;
         else
           edge_[1] = -1;
 
         // get cluster total charge
         const auto& stripCharges = cluster2->amplitudes();
         uint16_t charge = 0;
         for (uint i = 0; i < stripCharges.size(); i++) {
           charge += stripCharges[i];
         }
         clusterCharge_[1] = charge;
       }
       //cout << "strip cluster size2 = " << clusterWidthX_[0] << "  and size 2 = " << clusterWidthX_[1] << endl;
     }
 
     if (subDet2 == PixelSubdetector::PixelBarrel || subDet2 == PixelSubdetector::PixelEndcap) {  //pixel
 
       const TransientTrackingRecHit::ConstRecHitPointer thit1 = overlapHit.first->recHit();
       const SiPixelRecHit* recHitPix1 = dynamic_cast<const SiPixelRecHit*>((*thit1).hit());
       if (recHitPix1) {
         // check for cluster size and width
         SiPixelRecHit::ClusterRef const& cluster1 = recHitPix1->cluster();
 
         clusterSize_[0] = cluster1->size();
         clusterWidthX_[0] = cluster1->sizeX();
         clusterWidthY_[0] = cluster1->sizeY();
 
         // check for cluster at edge
         const PixelGeomDetUnit* theGeomDet = dynamic_cast<const PixelGeomDetUnit*>((*trackerGeometry_).idToDet(id1));
         const PixelTopology* topol = (&(theGeomDet->specificTopology()));
 
         int minPixelRow = cluster1->minPixelRow();  //x
         int maxPixelRow = cluster1->maxPixelRow();
         int minPixelCol = cluster1->minPixelCol();  //y
         int maxPixelCol = cluster1->maxPixelCol();
 
         bool edgeHitX = (topol->isItEdgePixelInX(minPixelRow)) || (topol->isItEdgePixelInX(maxPixelRow));
         bool edgeHitY = (topol->isItEdgePixelInY(minPixelCol)) || (topol->isItEdgePixelInY(maxPixelCol));
         if (edgeHitX || edgeHitY)
           edge_[0] = 1;
         else
           edge_[0] = -1;
 
         clusterCharge_[0] = (uint)cluster1->charge();
 
       } else {
         edm::LogWarning("Overlaps") << "didn't find pixel cluster" << endl;
         continue;
       }
 
       const TransientTrackingRecHit::ConstRecHitPointer thit2 = overlapHit.second->recHit();
       const SiPixelRecHit* recHitPix2 = dynamic_cast<const SiPixelRecHit*>((*thit2).hit());
       if (recHitPix2) {
         SiPixelRecHit::ClusterRef const& cluster2 = recHitPix2->cluster();
 
         clusterSize_[1] = cluster2->size();
         clusterWidthX_[1] = cluster2->sizeX();
         clusterWidthY_[1] = cluster2->sizeY();
         //cout << "second pixel cluster is " << clusterSize_[1] << " pixels with x width = " << clusterWidthX_[1] << " and y width = " << clusterWidthY_[1] << endl;
 
         const PixelGeomDetUnit* theGeomDet = dynamic_cast<const PixelGeomDetUnit*>((*trackerGeometry_).idToDet(id2));
         const PixelTopology* topol = (&(theGeomDet->specificTopology()));
 
         int minPixelRow = cluster2->minPixelRow();  //x
         int maxPixelRow = cluster2->maxPixelRow();
         int minPixelCol = cluster2->minPixelCol();  //y
         int maxPixelCol = cluster2->maxPixelCol();
 
         bool edgeHitX = (topol->isItEdgePixelInX(minPixelRow)) || (topol->isItEdgePixelInX(maxPixelRow));
         bool edgeHitY = (topol->isItEdgePixelInY(minPixelCol)) || (topol->isItEdgePixelInY(maxPixelCol));
         if (edgeHitX || edgeHitY)
           edge_[1] = 1;
         else
           edge_[1] = -1;
 
         clusterCharge_[1] = (uint)cluster2->charge();
 
       } else {
         edm::LogWarning("Overlaps") << "didn't find pixel cluster" << endl;
         continue;
       }
     }
 
     //also check for edge pixels
 
     //try writing out the SimHit info (for MC only)
     if (doSimHit_) {
       std::vector<PSimHit> psimHits1;
       std::vector<PSimHit> psimHits2;
       //calculate layer
       DetId id = overlapHit.first->recHit()->geographicalId();
       int layer(-1);
       layer = layerFromId(id, tTopo);
       int subDet = id.subdetId();
       edm::LogVerbatim("OverlapValidation") << "Subdet = " << subDet << " ; layer = " << layer;
 
       psimHits1 = associator.associateHit(*(firstRecHit->hit()));
       edm::LogVerbatim("OverlapValidation") << "single hit ";
       edm::LogVerbatim("OverlapValidation") << "length of psimHits1: " << psimHits1.size();
       if (!psimHits1.empty()) {
         float closest_dist = 99999.9;
         std::vector<PSimHit>::const_iterator closest_simhit = psimHits1.begin();
         for (std::vector<PSimHit>::const_iterator m = psimHits1.begin(); m < psimHits1.end(); m++) {
           //find closest simHit to the recHit
           float simX = (*m).localPosition().x();
           float dist = fabs(simX - (overlapHit.first->recHit()->localPosition().x()));
           edm::LogVerbatim("OverlapValidation")
               << "simHit1 simX = " << simX << "   hitX = " << overlapHit.first->recHit()->localPosition().x()
               << "   distX = " << dist << "   layer = " << layer;
           if (dist < closest_dist) {
             //cout << "found newest closest dist for simhit1" << endl;
             closest_dist = dist;
             closest_simhit = m;
           }
         }
         //if glued layer, convert sim hit position to matchedhit surface
         //layer index from 1-4 for TIB, 1-6 for TOB
         // Are the sim hits on the glued layers or are they split???
         if (subDet > 2 && !SiStripDetId(id).glued()) {
           const GluedGeomDet* gluedDet =
               (const GluedGeomDet*)(*trackerGeometry_).idToDet((*firstRecHit).hit()->geographicalId());
           const StripGeomDetUnit* stripDet = (StripGeomDetUnit*)gluedDet->monoDet();
           GlobalPoint gp = stripDet->surface().toGlobal((*closest_simhit).localPosition());
           LocalPoint lp = gluedDet->surface().toLocal(gp);
           LocalVector localdirection = (*closest_simhit).localDirection();
           GlobalVector globaldirection = stripDet->surface().toGlobal(localdirection);
           LocalVector direction = gluedDet->surface().toLocal(globaldirection);
           float scale = -lp.z() / direction.z();
           LocalPoint projectedPos = lp + scale * direction;
           simHitPositions_[0] = projectedPos.x();
           edm::LogVerbatim("OverlapValidation") << "simhit position from matched layer = " << simHitPositions_[0];
           simHitPositionsY_[0] = projectedPos.y();
         } else {
           simHitPositions_[0] = (*closest_simhit).localPosition().x();
           simHitPositionsY_[0] = (*closest_simhit).localPosition().y();
           edm::LogVerbatim("OverlapValidation") << "simhit position from non-matched layer = " << simHitPositions_[0];
         }
         edm::LogVerbatim("OverlapValidation") << "hit position = " << hitPositions_[0];
       } else {
         simHitPositions_[0] = -99.;
         simHitPositionsY_[0] = -99.;
         //cout << " filling simHitX: " << -99 << endl;
       }
 
       psimHits2 = associator.associateHit(*(secondRecHit->hit()));
       if (!psimHits2.empty()) {
         float closest_dist = 99999.9;
         std::vector<PSimHit>::const_iterator closest_simhit = psimHits2.begin();
         for (std::vector<PSimHit>::const_iterator m = psimHits2.begin(); m < psimHits2.end(); m++) {
           float simX = (*m).localPosition().x();
           float dist = fabs(simX - (overlapHit.second->recHit()->localPosition().x()));
           if (dist < closest_dist) {
             closest_dist = dist;
             closest_simhit = m;
           }
         }
         //if glued layer, convert sim hit position to matchedhit surface
         // if no sim hits on matched layers then this section can be removed
         if (subDet > 2 && !SiStripDetId(id).glued()) {
           const GluedGeomDet* gluedDet =
               (const GluedGeomDet*)(*trackerGeometry_).idToDet((*secondRecHit).hit()->geographicalId());
           const StripGeomDetUnit* stripDet = (StripGeomDetUnit*)gluedDet->monoDet();
           GlobalPoint gp = stripDet->surface().toGlobal((*closest_simhit).localPosition());
           LocalPoint lp = gluedDet->surface().toLocal(gp);
           LocalVector localdirection = (*closest_simhit).localDirection();
           GlobalVector globaldirection = stripDet->surface().toGlobal(localdirection);
           LocalVector direction = gluedDet->surface().toLocal(globaldirection);
           float scale = -lp.z() / direction.z();
           LocalPoint projectedPos = lp + scale * direction;
           simHitPositions_[1] = projectedPos.x();
           simHitPositionsY_[1] = projectedPos.y();
         } else {
           simHitPositions_[1] = (*closest_simhit).localPosition().x();
           simHitPositionsY_[1] = (*closest_simhit).localPosition().y();
         }
       } else {
         simHitPositions_[1] = -99.;
         simHitPositionsY_[1] = -99.;
       }
     }
     rootTree_->Fill();
   }
 }

◆ endJob()

void OverlapValidation::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 898 of file OverlapValidation.cc.

References rootTree_.

Referenced by o2olib.O2ORunMgr::executeJob().

                                {
   if (rootTree_) {
     rootTree_->GetDirectory()->cd();
     rootTree_->Write();
     delete rootTree_;
   }
 }

◆ fillDescriptions()

void OverlapValidation::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 274 of file OverlapValidation.cc.

References edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, TrackerHitAssociator::fillPSetDescription(), and ProducerED_cfi::InputTag.

                                                                                    {
   edm::ParameterSetDescription desc;
   desc.setComment("Overlap Validation analyzer plugin.");
   TrackerHitAssociator::fillPSetDescription(desc);
   //desc.add<bool>("associateStrip",false);  // NB this comes from the fillPSetDescription
   desc.addUntracked<int>("compressionSettings", -1);
   desc.add<edm::InputTag>("trajectories", edm::InputTag("FinalTrackRefitter"));
   desc.add<bool>("barrelOnly", false);
   desc.add<bool>("usePXB", true);
   desc.add<bool>("usePXF", true);
   desc.add<bool>("useTIB", true);
   desc.add<bool>("useTOB", true);
   desc.add<bool>("useTID", true);
   desc.add<bool>("useTEC", true);
   descriptions.addWithDefaultLabel(desc);
 }

◆ layerFromId()

int OverlapValidation::layerFromId	(	const DetId &	id,
		const TrackerTopology *const	tTopo
	)		const

private

Definition at line 890 of file OverlapValidation.cc.

References nano_mu_digi_cff::layer, and TrackerAlignableId::typeAndLayerFromDetId().

Referenced by analyzeTrajectory().

                                                                                             {
   TrackerAlignableId aliid;
   std::pair<int, int> subdetandlayer = aliid.typeAndLayerFromDetId(id, tTopo);
   int layer = subdetandlayer.second;
 
   return layer;
 }

Member Data Documentation

◆ acceptLayer

vector<bool> OverlapValidation::acceptLayer

private

Definition at line 149 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ addExtraBranches_

const bool OverlapValidation::addExtraBranches_

private

Definition at line 121 of file OverlapValidation.cc.

Referenced by OverlapValidation().

◆ barrelOnly_

bool OverlapValidation::barrelOnly_

private

Definition at line 152 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ chi2_

float OverlapValidation::chi2_[2]

private

Definition at line 128 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ chi2ProbCut_

const float OverlapValidation::chi2ProbCut_

private

Definition at line 123 of file OverlapValidation.cc.

Referenced by analyzeTrajectory().

◆ clusterCharge_

uint OverlapValidation::clusterCharge_[2]

private

Definition at line 146 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ clusterSize_

float OverlapValidation::clusterSize_[2]

private

Definition at line 145 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ clusterWidthX_

float OverlapValidation::clusterWidthX_[2]

private

Definition at line 143 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ clusterWidthY_

float OverlapValidation::clusterWidthY_[2]

private

Definition at line 144 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ combiner_

TrajectoryStateCombiner OverlapValidation::combiner_

private

Definition at line 115 of file OverlapValidation.cc.

Referenced by analyzeTrajectory().

◆ compressionSettings_

const int OverlapValidation::compressionSettings_

private

Definition at line 120 of file OverlapValidation.cc.

Referenced by OverlapValidation().

◆ config_

edm::ParameterSet OverlapValidation::config_

private

Definition at line 108 of file OverlapValidation.cc.

Referenced by analyze().

◆ detInfo_

SiStripDetInfo OverlapValidation::detInfo_

private

Definition at line 110 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ doSimHit_

bool OverlapValidation::doSimHit_

private

Definition at line 111 of file OverlapValidation.cc.

Referenced by analyze(), analyzeTrajectory(), and OverlapValidation().

◆ edge_

int OverlapValidation::edge_[2]

private

Definition at line 147 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ event_

uint OverlapValidation::event_

private

Definition at line 151 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

◆ FileInPath_

edm::FileInPath OverlapValidation::FileInPath_

private

Definition at line 119 of file OverlapValidation.cc.

Referenced by OverlapValidation().

◆ geomToken_

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> OverlapValidation::geomToken_

private

Definition at line 104 of file OverlapValidation.cc.

Referenced by analyze().

◆ hitCounts_

unsigned short OverlapValidation::hitCounts_[2]

private

Definition at line 127 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ hitErrors_

float OverlapValidation::hitErrors_[2]

private

Definition at line 138 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ hitErrorsY_

float OverlapValidation::hitErrorsY_[2]

private

Definition at line 140 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ hitPositions_

float OverlapValidation::hitPositions_[2]

private

Definition at line 137 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ hitPositionsY_

float OverlapValidation::hitPositionsY_[2]

private

Definition at line 139 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ layer_

uint OverlapValidation::layer_

private

Definition at line 126 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ localxdotglobalphi_

float OverlapValidation::localxdotglobalphi_[2]

private

Definition at line 163 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ localxdotglobalr_

float OverlapValidation::localxdotglobalr_[2]

private

Definition at line 164 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ localxdotglobalx_

float OverlapValidation::localxdotglobalx_[2]

private

Definition at line 166 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ localxdotglobaly_

float OverlapValidation::localxdotglobaly_[2]

private

Definition at line 167 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ localxdotglobalz_

float OverlapValidation::localxdotglobalz_[2]

private

Definition at line 165 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ localydotglobalphi_

float OverlapValidation::localydotglobalphi_[2]

private

Definition at line 168 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ localydotglobalr_

float OverlapValidation::localydotglobalr_[2]

private

Definition at line 169 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ localydotglobalx_

float OverlapValidation::localydotglobalx_[2]

private

Definition at line 171 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ localydotglobaly_

float OverlapValidation::localydotglobaly_[2]

private

Definition at line 172 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ localydotglobalz_

float OverlapValidation::localydotglobalz_[2]

private

Definition at line 170 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ magField_

const MagneticField* OverlapValidation::magField_

private

Definition at line 113 of file OverlapValidation.cc.

Referenced by analyze(), and analyzeTrajectory().

◆ magFieldToken_

const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> OverlapValidation::magFieldToken_

private

Definition at line 105 of file OverlapValidation.cc.

Referenced by analyze().

◆ minHitsCut_

const int OverlapValidation::minHitsCut_

private

Definition at line 122 of file OverlapValidation.cc.

Referenced by analyzeTrajectory().

◆ moduleX_

float OverlapValidation::moduleX_[2]

private

Definition at line 155 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ moduleY_

float OverlapValidation::moduleY_[2]

private

Definition at line 156 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ moduleZ_

float OverlapValidation::moduleZ_[2]

private

Definition at line 157 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ momentum_

float OverlapValidation::momentum_

private

Definition at line 150 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ onezero

const Point2DBase<float, LocalTag> OverlapValidation::onezero = Point2DBase<float, LocalTag>(1, 0)

private

Definition at line 160 of file OverlapValidation.cc.

Referenced by analyzeTrajectory().

◆ overlapCounts_

int OverlapValidation::overlapCounts_[3]

private

Definition at line 116 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), OverlapValidation(), and ~OverlapValidation().

◆ overlapIds_

unsigned int OverlapValidation::overlapIds_[2]

private

Definition at line 129 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ overlapPath_

float OverlapValidation::overlapPath_

private

Definition at line 125 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ predictedDeltaXError_

float OverlapValidation::predictedDeltaXError_

private

Definition at line 133 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ predictedDeltaYError_

float OverlapValidation::predictedDeltaYError_

private

Definition at line 134 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ predictedLocalErrors_

float OverlapValidation::predictedLocalErrors_[5][2]

private

Definition at line 132 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ predictedLocalParameters_

float OverlapValidation::predictedLocalParameters_[5][2]

private

Definition at line 131 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ predictedPositions_

float OverlapValidation::predictedPositions_[3][2]

private

Definition at line 130 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ relativeXSign_

char OverlapValidation::relativeXSign_

private

Definition at line 135 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ relativeYSign_

char OverlapValidation::relativeYSign_

private

Definition at line 136 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ rootTree_

TTree* OverlapValidation::rootTree_

private

Definition at line 118 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), endJob(), and OverlapValidation().

◆ run_

uint OverlapValidation::run_

private

Definition at line 151 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

◆ simHitPositions_

float OverlapValidation::simHitPositions_[2]

private

Definition at line 141 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ simHitPositionsY_

float OverlapValidation::simHitPositionsY_[2]

private

Definition at line 142 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ subdetID

int OverlapValidation::subdetID

private

Definition at line 158 of file OverlapValidation.cc.

Referenced by analyzeTrajectory(), and OverlapValidation().

◆ topoToken_

const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> OverlapValidation::topoToken_

private

Definition at line 106 of file OverlapValidation.cc.

Referenced by analyze().

◆ trackerGeometry_

const TrackerGeometry* OverlapValidation::trackerGeometry_

private

Definition at line 112 of file OverlapValidation.cc.

Referenced by analyze().

◆ trajectoryTag_

edm::InputTag OverlapValidation::trajectoryTag_

private

Definition at line 109 of file OverlapValidation.cc.

Referenced by OverlapValidation().

◆ trajectoryToken_

edm::EDGetTokenT<TrajectoryCollection> OverlapValidation::trajectoryToken_

private

Definition at line 91 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

◆ zeroone

const Point2DBase<float, LocalTag> OverlapValidation::zeroone = Point2DBase<float, LocalTag>(0, 1)

private

Definition at line 161 of file OverlapValidation.cc.

Referenced by analyzeTrajectory().

◆ zerozero

const Point2DBase<float, LocalTag> OverlapValidation::zerozero = Point2DBase<float, LocalTag>(0, 0)

private

Definition at line 159 of file OverlapValidation.cc.

Referenced by analyzeTrajectory().

Public Member Functions

Static Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ ConstRecHitPointer

◆ TrajectoryCollection

Constructor & Destructor Documentation

◆ OverlapValidation()

◆ ~OverlapValidation()

Member Function Documentation

◆ analyze()

◆ analyzeTrajectory()

◆ endJob()

◆ fillDescriptions()

◆ layerFromId()

Member Data Documentation

◆ acceptLayer

◆ addExtraBranches_

◆ barrelOnly_

◆ chi2_

◆ chi2ProbCut_

◆ clusterCharge_

◆ clusterSize_

◆ clusterWidthX_

◆ clusterWidthY_

◆ combiner_

◆ compressionSettings_

◆ config_

◆ detInfo_

◆ doSimHit_

◆ edge_

◆ event_

◆ FileInPath_

◆ geomToken_

◆ hitCounts_

◆ hitErrors_

◆ hitErrorsY_

◆ hitPositions_

◆ hitPositionsY_

◆ layer_

◆ localxdotglobalphi_

◆ localxdotglobalr_

◆ localxdotglobalx_

◆ localxdotglobaly_

◆ localxdotglobalz_

◆ localydotglobalphi_

◆ localydotglobalr_

◆ localydotglobalx_

◆ localydotglobaly_

◆ localydotglobalz_

◆ magField_

◆ magFieldToken_

◆ minHitsCut_

◆ moduleX_

◆ moduleY_

◆ moduleZ_

◆ momentum_

◆ onezero

◆ overlapCounts_

◆ overlapIds_

◆ overlapPath_

◆ predictedDeltaXError_

◆ predictedDeltaYError_

◆ predictedLocalErrors_

◆ predictedLocalParameters_

◆ predictedPositions_

◆ relativeXSign_

◆ relativeYSign_

◆ rootTree_

◆ run_

◆ simHitPositions_

◆ simHitPositionsY_

◆ subdetID

◆ topoToken_

◆ trackerGeometry_

◆ trajectoryTag_