OverlapValidation Class Reference

Inheritance diagram for OverlapValidation:

Public Member Functions
	OverlapValidation (const edm::ParameterSet &)
	~OverlapValidation () override
Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () 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
ESProxyIndex const *	esGetTokenIndices (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::vector< ModuleDescription const * > &modules, 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
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Types
typedef TransientTrackingRecHit::ConstRecHitPointer	ConstRecHitPointer
typedef vector< Trajectory >	TrajectoryCollection

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
virtual void	analyze (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_
edm::ParameterSet	config_
bool	doSimHit_
int	edge_ [2]
uint	event_
edm::FileInPath	FileInPath_
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 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]
SiStripDetInfoFileReader *	reader
char	relativeXSign_
char	relativeYSign_
TTree *	rootTree_
uint	run_
float	simHitPositions_ [2]
float	simHitPositionsY_ [2]
int	subdetID
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
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
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 ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
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<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)

Detailed Description

Definition at line 89 of file OverlapValidation.cc.

Member Typedef Documentation

typedef TransientTrackingRecHit::ConstRecHitPointer OverlapValidation::ConstRecHitPointer

private

Definition at line 100 of file OverlapValidation.cc.

typedef vector<Trajectory> OverlapValidation::TrajectoryCollection

private

Definition at line 95 of file OverlapValidation.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 187 of file OverlapValidation.cc.

References acceptLayer, addExtraBranches_, barrelOnly_, chi2_, clusterCharge_, clusterSize_, clusterWidthX_, clusterWidthY_, edm::ConsumesCollector::consumes(), edm::EDConsumerBase::consumesCollector(), doSimHit_, edge_, event_, FileInPath_, edm::FileInPath::fullPath(), edm::ParameterSet::getParameter(), hitCounts_, hitErrors_, hitErrorsY_, hitPositions_, hitPositionsY_, 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_, reader, relativeXSign_, relativeYSign_, rootTree_, run_, simHitPositions_, simHitPositionsY_, subdetID, StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, trajectoryTag_, and trajectoryToken_.

                                                                   :
  config_(iConfig), rootTree_(nullptr),
  FileInPath_("CalibTracker/SiStripCommon/data/SiStripDetInfo.dat"),
  addExtraBranches_(false),
  minHitsCut_(6),
  chi2ProbCut_(0.001)
{
  
  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");
  reader=new SiStripDetInfoFileReader(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
  //
  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 ( )

override

Definition at line 271 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

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

overrideprivate

Definition at line 291 of file OverlapValidation.cc.

References AnalyticalPropagator_cfi::AnalyticalPropagator, anyDirection, ctfWithMaterialTrackMCMatch_cfi::associator, config_, edm::EDConsumerBase::consumesCollector(), doSimHit_, edm::EventID::event(), event_, edm::EventSetup::get(), edm::Event::getByToken(), edm::EventBase::id(), magField_, edm::Handle< T >::product(), edm::ESHandle< T >::product(), PhotonConversionTrajectorySeedProducerFromQuadruplets_cfi::propagator, edm::EventID::run(), run_, trackerGeometry_, and trajectoryToken_.

{
  using namespace edm;
  //
  // mag field & search tracker
  //
  edm::ESHandle<MagneticField> magFieldHandle;
  iSetup.get<IdealMagneticFieldRecord>().get(magFieldHandle);
  magField_ = magFieldHandle.product();
  //
  // propagator
  //
  AnalyticalPropagator propagator(magField_,anyDirection);
  //
  // geometry
  

  //
  edm::ESHandle<TrackerGeometry> geometryHandle;
  iSetup.get<TrackerDigiGeometryRecord>().get(geometryHandle);
  trackerGeometry_ = geometryHandle.product();
  //
  // 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
 edm::ESHandle<TrackerTopology> tTopoHandle;
 iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
 const TrackerTopology* const tTopo = tTopoHandle.product(); 
 for (const auto& trajectory : *trajectoryCollection)  analyze(trajectory,propagator,*associator, tTopo);
  
  run_ = iEvent.id().run();
  event_ = iEvent.id().event();
  

  
}

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

privatevirtual

Definition at line 358 of file OverlapValidation.cc.

References funct::abs(), acceptLayer, TrackerHitAssociator::associateHit(), barrelOnly_, ALCARECOTkAlJpsiMuMu_cff::charge, chi2_, chi2ProbCut_, Trajectory::chiSquared(), ChiSquaredProbability(), SiStripRecHit1D::cluster(), SiPixelRecHit::cluster(), clusterCharge_, clusterSize_, clusterWidthX_, clusterWidthY_, TrajectoryStateCombiner::combine(), combiner_, diffTreeTool::diff, doSimHit_, edge_, Trajectory::foundHits(), SiStripDetInfoFileReader::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(), JacobianCurvilinearToLocal::jacobian(), layer_, layerFromId(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), localxdotglobalphi_, localxdotglobalr_, localxdotglobalx_, localxdotglobaly_, localxdotglobalz_, localydotglobalphi_, localydotglobalr_, localydotglobalx_, localydotglobaly_, localydotglobalz_, Trajectory::lostHits(), funct::m, PV3DBase< T, PVType, FrameType >::mag(), magField_, LocalTrajectoryError::matrix(), Trajectory::measurements(), minHitsCut_, moduleX_, moduleY_, moduleZ_, momentum_, GluedGeomDet::monoDet(), Trajectory::ndof(), onezero, muons1stStep_cfi::Overlap, overlapCounts_, overlapIds_, overlapPath_, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, position, predictedDeltaXError_, predictedDeltaYError_, predictedLocalErrors_, predictedLocalParameters_, predictedPositions_, Propagator::propagate(), Propagator::propagateWithPath(), DetId::rawId(), reader, relativeXSign_, relativeYSign_, rootTree_, Scenarios_cff::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.

{
  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=reader->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 std::vector<uint8_t>& stripCharges = cluster1->amplitudes();
    uint16_t charge = 0;
    for (uint i = 0; i < stripCharges.size(); i++) {
      charge += stripCharges.at(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=reader->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 std::vector<uint8_t>& stripCharges = cluster2->amplitudes();
    uint16_t charge = 0;
    for (uint i = 0; i < stripCharges.size(); i++) {
      charge += stripCharges.at(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();
  }
}

void OverlapValidation::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 865 of file OverlapValidation.cc.

References DEFINE_FWK_MODULE, and rootTree_.

Referenced by o2olib.O2ORunMgr::executeJob().

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

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

private

Definition at line 853 of file OverlapValidation.cc.

References TrackerAlignableId::typeAndLayerFromDetId().

Referenced by analyze().

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


  return layer;
  
}

Member Data Documentation

vector<bool> OverlapValidation::acceptLayer

private

Definition at line 148 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

const bool OverlapValidation::addExtraBranches_

private

Definition at line 120 of file OverlapValidation.cc.

Referenced by OverlapValidation().

bool OverlapValidation::barrelOnly_

private

Definition at line 151 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::chi2_[2]

private

Definition at line 127 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

const float OverlapValidation::chi2ProbCut_

private

Definition at line 122 of file OverlapValidation.cc.

Referenced by analyze().

uint OverlapValidation::clusterCharge_[2]

private

Definition at line 145 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::clusterSize_[2]

private

Definition at line 144 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::clusterWidthX_[2]

private

Definition at line 142 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::clusterWidthY_[2]

private

Definition at line 143 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

TrajectoryStateCombiner OverlapValidation::combiner_

private

Definition at line 115 of file OverlapValidation.cc.

Referenced by analyze().

edm::ParameterSet OverlapValidation::config_

private

Definition at line 108 of file OverlapValidation.cc.

Referenced by analyze().

bool OverlapValidation::doSimHit_

private

Definition at line 111 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

int OverlapValidation::edge_[2]

private

Definition at line 146 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

uint OverlapValidation::event_

private

Definition at line 150 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

edm::FileInPath OverlapValidation::FileInPath_

private

Definition at line 119 of file OverlapValidation.cc.

Referenced by OverlapValidation().

unsigned short OverlapValidation::hitCounts_[2]

private

Definition at line 126 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::hitErrors_[2]

private

Definition at line 137 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::hitErrorsY_[2]

private

Definition at line 139 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::hitPositions_[2]

private

Definition at line 136 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::hitPositionsY_[2]

private

Definition at line 138 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

uint OverlapValidation::layer_

private

Definition at line 125 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::localxdotglobalphi_[2]

private

Definition at line 162 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::localxdotglobalr_[2]

private

Definition at line 163 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::localxdotglobalx_[2]

private

Definition at line 165 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::localxdotglobaly_[2]

private

Definition at line 166 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::localxdotglobalz_[2]

private

Definition at line 164 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::localydotglobalphi_[2]

private

Definition at line 167 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::localydotglobalr_[2]

private

Definition at line 168 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::localydotglobalx_[2]

private

Definition at line 170 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::localydotglobaly_[2]

private

Definition at line 171 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::localydotglobalz_[2]

private

Definition at line 169 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

const MagneticField* OverlapValidation::magField_

private

Definition at line 113 of file OverlapValidation.cc.

Referenced by analyze().

const int OverlapValidation::minHitsCut_

private

Definition at line 121 of file OverlapValidation.cc.

Referenced by analyze().

float OverlapValidation::moduleX_[2]

private

Definition at line 154 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::moduleY_[2]

private

Definition at line 155 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::moduleZ_[2]

private

Definition at line 156 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::momentum_

private

Definition at line 149 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

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

private

Definition at line 159 of file OverlapValidation.cc.

Referenced by analyze().

int OverlapValidation::overlapCounts_[3]

private

Definition at line 116 of file OverlapValidation.cc.

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

unsigned int OverlapValidation::overlapIds_[2]

private

Definition at line 128 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::overlapPath_

private

Definition at line 124 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::predictedDeltaXError_

private

Definition at line 132 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::predictedDeltaYError_

private

Definition at line 133 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::predictedLocalErrors_[5][2]

private

Definition at line 131 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::predictedLocalParameters_[5][2]

private

Definition at line 130 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::predictedPositions_[3][2]

private

Definition at line 129 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

SiStripDetInfoFileReader* OverlapValidation::reader

private

Definition at line 110 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

char OverlapValidation::relativeXSign_

private

Definition at line 134 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

char OverlapValidation::relativeYSign_

private

Definition at line 135 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

TTree* OverlapValidation::rootTree_

private

Definition at line 118 of file OverlapValidation.cc.

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

uint OverlapValidation::run_

private

Definition at line 150 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::simHitPositions_[2]

private

Definition at line 140 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

float OverlapValidation::simHitPositionsY_[2]

private

Definition at line 141 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

int OverlapValidation::subdetID

private

Definition at line 157 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

const TrackerGeometry* OverlapValidation::trackerGeometry_

private

Definition at line 112 of file OverlapValidation.cc.

Referenced by analyze().

edm::InputTag OverlapValidation::trajectoryTag_

private

Definition at line 109 of file OverlapValidation.cc.

Referenced by OverlapValidation().

edm::EDGetTokenT<TrajectoryCollection> OverlapValidation::trajectoryToken_

private

Definition at line 99 of file OverlapValidation.cc.

Referenced by analyze(), and OverlapValidation().

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

private

Definition at line 160 of file OverlapValidation.cc.

Referenced by analyze().

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

private

Definition at line 158 of file OverlapValidation.cc.

Referenced by analyze().