TestTrackHits Class Reference

#include <RecoTracker/TestTrackHits/src/TestTrackHits.cc>

Inheritance diagram for TestTrackHits:

Public Member Functions
	TestTrackHits (const edm::ParameterSet &)
	~TestTrackHits () override
Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns >
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
bool	wantsGlobalLuminosityBlocks () const noexcept final
bool	wantsGlobalRuns () const noexcept final
bool	wantsInputProcessBlocks () const noexcept final
bool	wantsProcessBlocks () const noexcept 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 noexcept
bool	wantsStreamRuns () const noexcept
	~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)

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	beginRun (edm::Run const &run, const edm::EventSetup &) override
template<unsigned int D>
double	computeChi2Increment (MeasurementExtractor, TransientTrackingRecHit::ConstRecHitPointer)
double	computeChi2Increment (MeasurementExtractor me, TransientTrackingRecHit::ConstRecHitPointer hit)
void	endJob () override
void	endRun (edm::Run const &run, const edm::EventSetup &) override
std::pair< LocalPoint, LocalVector >	projectHit (const PSimHit &, const StripGeomDetUnit *, const BoundPlane &)

Private Attributes
std::string	builderName
TH1F *	energyLossRatio
TFile *	file
TH1F *	goodbadmerged
std::map< std::string, TH1F * >	hChi2BadHit
std::map< std::string, TH1F * >	hChi2DeltaHit
std::map< std::string, TH1F * >	hChi2GoodHit
std::map< std::string, TH1F * >	hChi2Increment
std::map< std::string, TH1F * >	hChi2Increment_mono
std::map< std::string, TH1F * >	hChi2Increment_stereo
std::map< std::string, TH2F * >	hChi2IncrementVsEta
std::map< std::string, TH1F * >	hChi2NSharedHit
std::map< std::string, TH1F * >	hChi2SharedHit
TH2F *	hClsize_vs_Chi2
TH1F *	hClusterSize
TH2F *	hGoodHit_vs_Chi2
TH2F *	hPixClsize_vs_Chi2
TH1F *	hPixClusterSize
TH1F *	hPixSimHitVecSize
TH2F *	hPrjClsize_vs_Chi2
TH1F *	hPrjClusterSize
TH1F *	hPrjSimHitVecSize
TH2F *	hProcess_vs_Chi2
std::map< std::string, TH1F * >	hPullGM_X_ts
std::map< std::string, TH1F * >	hPullGM_X_ts_mono
std::map< std::string, TH1F * >	hPullGM_X_ts_stereo
std::map< std::string, TH1F * >	hPullGM_Y_ts
std::map< std::string, TH1F * >	hPullGM_Y_ts_mono
std::map< std::string, TH1F * >	hPullGM_Y_ts_stereo
std::map< std::string, TH1F * >	hPullGM_Z_ts
std::map< std::string, TH1F * >	hPullGM_Z_ts_mono
std::map< std::string, TH1F * >	hPullGM_Z_ts_stereo
std::map< std::string, TH1F * >	hPullGP_X_rs
std::map< std::string, TH1F * >	hPullGP_X_rs_mono
std::map< std::string, TH1F * >	hPullGP_X_rs_stereo
std::map< std::string, TH1F * >	hPullGP_X_tr
std::map< std::string, TH1F * >	hPullGP_X_tr_mono
std::map< std::string, TH1F * >	hPullGP_X_tr_stereo
std::map< std::string, TH1F * >	hPullGP_X_ts
std::map< std::string, TH1F * >	hPullGP_X_ts_mono
std::map< std::string, TH1F * >	hPullGP_X_ts_stereo
std::map< std::string, TH1F * >	hPullGP_Y_rs
std::map< std::string, TH1F * >	hPullGP_Y_rs_mono
std::map< std::string, TH1F * >	hPullGP_Y_rs_stereo
std::map< std::string, TH1F * >	hPullGP_Y_tr
std::map< std::string, TH1F * >	hPullGP_Y_tr_mono
std::map< std::string, TH1F * >	hPullGP_Y_tr_stereo
std::map< std::string, TH1F * >	hPullGP_Y_ts
std::map< std::string, TH1F * >	hPullGP_Y_ts_mono
std::map< std::string, TH1F * >	hPullGP_Y_ts_stereo
std::map< std::string, TH1F * >	hPullGP_Z_rs
std::map< std::string, TH1F * >	hPullGP_Z_rs_mono
std::map< std::string, TH1F * >	hPullGP_Z_rs_stereo
std::map< std::string, TH1F * >	hPullGP_Z_tr
std::map< std::string, TH1F * >	hPullGP_Z_tr_mono
std::map< std::string, TH1F * >	hPullGP_Z_tr_stereo
std::map< std::string, TH1F * >	hPullGP_Z_ts
std::map< std::string, TH1F * >	hPullGP_Z_ts_mono
std::map< std::string, TH1F * >	hPullGP_Z_ts_stereo
TH1F *	hSimHitVecSize
TH2F *	hSt1Clsize_vs_Chi2
TH1F *	hSt1ClusterSize
TH1F *	hSt1SimHitVecSize
TH2F *	hSt2Clsize_vs_Chi2
TH1F *	hSt2ClusterSize
TH1F *	hSt2SimHitVecSize
TH1F *	hTotChi2BadHit
TH1F *	hTotChi2DeltaHit
TH1F *	hTotChi2GoodHit
TH1F *	hTotChi2Increment
TH1F *	hTotChi2NSharedHit
TH1F *	hTotChi2SharedHit
TH1F *	mergedPull
std::string	out
TH1F *	probXbad
TH1F *	probXdelta
TH1F *	probXgood
TH1F *	probXnoshare
TH1F *	probXshared
TH1F *	probYbad
TH1F *	probYdelta
TH1F *	probYgood
TH1F *	probYnoshare
TH1F *	probYshared
std::string	propagatorName
std::string	srcName
edm::ESHandle< TransientTrackingRecHitBuilder >	theBuilder
edm::ESGetToken< TransientTrackingRecHitBuilder, TransientRecHitRecord >	theBuilderToken
edm::ESHandle< TrackerGeometry >	theG
edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	theGToken
edm::ESHandle< MagneticField >	theMF
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	theMFToken
edm::ESHandle< Propagator >	thePropagator
edm::ESGetToken< Propagator, TrackingComponentsRecord >	thePropagatorToken
edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	theTopoToken
edm::ESHandle< TrajectoryStateUpdator >	theUpdator
edm::ESGetToken< TrajectoryStateUpdator, TrackingComponentsRecord >	theUpdatorToken
std::stringstream	title
std::string	tpName
edm::Handle< reco::TrackToTrackingParticleAssociator >	trackAssociator
edm::Handle< edm::View< reco::Track > >	trackCollectionHandle
TrackerHitAssociator::Config	trackerHitAssociatorConfig_
edm::Handle< TrackingParticleCollection >	trackingParticleCollectionHandle
edm::Handle< std::vector< Trajectory > >	trajCollectionHandle
edm::Handle< TrajTrackAssociationCollection >	trajTrackAssociationCollectionHandle
std::string	updatorName

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)
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

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 54 of file TestTrackHits.h.

Constructor & Destructor Documentation

TestTrackHits()

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

explicit

Definition at line 18 of file TestTrackHits.cc.

References builderName, edm::EDConsumerBase::esConsumes(), edm::ParameterSet::getParameter(), LogTrace, out, propagatorName, srcName, AlCaHLTBitMon_QueryRunRegistry::string, theBuilderToken, theGToken, theMFToken, thePropagatorToken, theTopoToken, theUpdatorToken, tpName, and updatorName.

                                                           : trackerHitAssociatorConfig_(consumesCollector()) {
  LogTrace("TestTrackHits") << iConfig;
  propagatorName = iConfig.getParameter<std::string>("Propagator");
  builderName = iConfig.getParameter<std::string>("TTRHBuilder");
  srcName = iConfig.getParameter<std::string>("src");
  tpName = iConfig.getParameter<std::string>("tpname");
  updatorName = iConfig.getParameter<std::string>("updator");
  out = iConfig.getParameter<std::string>("out");
  theGToken = esConsumes<edm::Transition::BeginRun>();
  theMFToken = esConsumes<edm::Transition::BeginRun>();
  thePropagatorToken = esConsumes<edm::Transition::BeginRun>(edm::ESInputTag("", propagatorName));
  theBuilderToken = esConsumes<edm::Transition::BeginRun>(edm::ESInputTag("", builderName));
  theUpdatorToken = esConsumes<edm::Transition::BeginRun>(edm::ESInputTag("", updatorName));
  theTopoToken = esConsumes();
  //   ParameterSet cuts = iConfig.getParameter<ParameterSet>("RecoTracksCuts");
  //   selectRecoTracks = RecoTrackSelector(cuts.getParameter<double>("ptMin"),
  //                       cuts.getParameter<double>("minRapidity"),
  //                       cuts.getParameter<double>("maxRapidity"),
  //                       cuts.getParameter<double>("tip"),
  //                       cuts.getParameter<double>("lip"),
  //                       cuts.getParameter<int>("minHit"),
  //                       cuts.getParameter<double>("maxChi2"));
}

~TestTrackHits()

TestTrackHits::~TestTrackHits ( )

override

Definition at line 42 of file TestTrackHits.cc.

{}

Member Function Documentation

analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 260 of file TestTrackHits.cc.

References funct::abs(), anyDirection, reco::TrackToTrackingParticleAssociator::associateRecoToSim(), pwdgSkimBPark_cfi::beamSpot, TransientTrackingRecHitBuilder::build(), TrajectoryStateOnSurface::cartesianError(), heavyFlavorValidationHarvestingSequence_cff::combiner, computeChi2Increment(), submitPVResolutionJobs::count, GlobalErrorBase< T, ErrorWeightType >::cxx(), GlobalErrorBase< T, ErrorWeightType >::cyy(), GlobalErrorBase< T, ErrorWeightType >::czz(), dumpMFGeometry_cfg::delta, DetId::det(), edm::AssociationMap< Tag >::end(), energyLossRatio, spr::find(), edm::AssociationMap< Tag >::find(), dqmdumpme::first, edm::EventSetup::getHandle(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), goodbadmerged, hChi2BadHit, hChi2DeltaHit, hChi2GoodHit, hChi2Increment, hChi2Increment_mono, hChi2Increment_stereo, hChi2IncrementVsEta, hChi2NSharedHit, hChi2SharedHit, hClsize_vs_Chi2, hClusterSize, hGoodHit_vs_Chi2, TrackQuality_cff::hitAssociator, hPixClsize_vs_Chi2, hPixClusterSize, hPixSimHitVecSize, hPrjClsize_vs_Chi2, hPrjClusterSize, hPrjSimHitVecSize, hProcess_vs_Chi2, hPullGM_X_ts, hPullGM_X_ts_mono, hPullGM_X_ts_stereo, hPullGM_Y_ts, hPullGM_Y_ts_mono, hPullGM_Y_ts_stereo, hPullGM_Z_ts, hPullGM_Z_ts_mono, hPullGM_Z_ts_stereo, hPullGP_X_rs, hPullGP_X_rs_mono, hPullGP_X_rs_stereo, hPullGP_X_tr, hPullGP_X_tr_mono, hPullGP_X_tr_stereo, hPullGP_X_ts, hPullGP_X_ts_mono, hPullGP_X_ts_stereo, hPullGP_Y_rs, hPullGP_Y_rs_mono, hPullGP_Y_rs_stereo, hPullGP_Y_tr, hPullGP_Y_tr_mono, hPullGP_Y_tr_stereo, hPullGP_Y_ts, hPullGP_Y_ts_mono, hPullGP_Y_ts_stereo, hPullGP_Z_rs, hPullGP_Z_rs_mono, hPullGP_Z_rs_stereo, hPullGP_Z_tr, hPullGP_Z_tr_mono, hPullGP_Z_tr_stereo, hPullGP_Z_ts, hPullGP_Z_ts_mono, hPullGP_Z_ts_stereo, hSimHitVecSize, hSt1Clsize_vs_Chi2, hSt1ClusterSize, hSt1SimHitVecSize, hSt2Clsize_vs_Chi2, hSt2ClusterSize, hSt2SimHitVecSize, hTotChi2BadHit, hTotChi2DeltaHit, hTotChi2GoodHit, hTotChi2Increment, hTotChi2NSharedHit, hTotChi2SharedHit, mps_fire::i, iEvent, TrajectoryStateOnSurface::isValid(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, dqmiolumiharvest::j, findQualityFiles::jj, TrackerTopology::layer(), PSimHit::localPosition(), LogDebug, LogTrace, visualization-live-secondInstance_cfg::m, CartesianTrajectoryError::matrix(), mergedPull, PSimHit::momentumAtEntry(), GloballyPositioned< T >::position(), CartesianTrajectoryError::position(), createTree::pp, SiPixelRecHit::probabilityX(), SiPixelRecHit::probabilityY(), probXbad, probXdelta, probXgood, probXnoshare, probXshared, probYbad, probYdelta, probYgood, probYnoshare, probYshared, PSimHit::processType(), edm::ESHandle< T >::product(), projectHit(), Propagator::propagate(), alignCSCRings::s, jetUpdater_cfi::sort, mathSSE::sqrt(), srcName, theBuilder, theMF, theTopoToken, title, Surface::toGlobal(), cmsswSequenceInfo::tp, tpName, HLT_2024v14_cff::track, trackAssociator, trackCollectionHandle, trackerHitAssociatorConfig_, trackingParticleCollectionHandle, trajCollectionHandle, trajTrackAssociationCollectionHandle, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), geometryCSVtoXML::yy, and PV3DBase< T, PVType, FrameType >::z().

                                                                               {
  //Retrieve tracker topology from geometry
  edm::ESHandle<TrackerTopology> tTopo = iSetup.getHandle(theTopoToken);

  LogDebug("TestTrackHits") << "new event";
  iEvent.getByLabel(srcName, trajCollectionHandle);
  iEvent.getByLabel(srcName, trackCollectionHandle);
  iEvent.getByLabel(srcName, trajTrackAssociationCollectionHandle);
  iEvent.getByLabel(tpName, trackingParticleCollectionHandle);
  edm::Handle<reco::BeamSpot> beamSpot;
  iEvent.getByLabel("offlineBeamSpot", beamSpot);

  LogTrace("TestTrackHits") << "Tr collection size=" << trackCollectionHandle->size();
  LogTrace("TestTrackHits") << "TP collection size=" << trackingParticleCollectionHandle->size();

  iEvent.getByLabel("trackAssociatorByHits", trackAssociator);

  TrackerHitAssociator hitAssociator(iEvent, trackerHitAssociatorConfig_);

  reco::RecoToSimCollection recSimColl =
      trackAssociator->associateRecoToSim(trackCollectionHandle, trackingParticleCollectionHandle);

  TrajectoryStateCombiner combiner;

  int evtHits = 0;
  int i = 0;
  int yy = 0;
  int yyy = 0;
  for (std::vector<Trajectory>::const_iterator it = trajCollectionHandle->begin(); it != trajCollectionHandle->end();
       it++) {
    LogTrace("TestTrackHits") << "\n*****************new trajectory********************";
    double tchi2 = 0;

    std::vector<TrajectoryMeasurement> tmColl = it->measurements();

    edm::Ref<std::vector<Trajectory> > traj(trajCollectionHandle, i);
    reco::TrackRef tmptrack = (*trajTrackAssociationCollectionHandle.product())[traj];
    edm::RefToBase<reco::Track> track(tmptrack);
    //     if ( !selectRecoTracks( *track,beamSpot.product() ) ) {
    //       LogTrace("TestTrackHits") << "track does not pass quality cuts: skippingtrack #" << ++yyy;
    //       i++;
    //       continue;
    //     }

    std::vector<std::pair<TrackingParticleRef, double> > tP;
    if (recSimColl.find(track) != recSimColl.end()) {
      tP = recSimColl[track];
      if (!tP.empty()) {
        edm::LogVerbatim("TestTrackHits") << "reco::Track #" << ++yyy << " with pt=" << track->pt()
                                          << " associated with quality:" << tP.begin()->second << " good track #"
                                          << ++yy << " has hits:" << track->numberOfValidHits() << "\n";
      }
    } else {
      edm::LogVerbatim("TestTrackHits") << "reco::Track #" << ++yyy << " with pt=" << track->pt()
                                        << " NOT associated to any TrackingParticle"
                                        << "\n";
      i++;
      continue;
    }
    //     if(recSimColl.find(track) != recSimColl.end()) {
    //       tP = recSimColl[track];
    //     } else {
    //       LogTrace("TestTrackHits") << "fake track: skipping track " << ++yyy;
    //       continue;//skip fake tracks
    //     }
    //     if (tP.size()==0) {
    //       LogTrace("TestTrackHits") << "fake track: skipping track " << ++yyy;
    //       continue;
    //     }
    TrackingParticleRef tp = tP.begin()->first;
    LogTrace("TestTrackHits") << "a tp is associated with fraction=" << tP.begin()->second;
    //LogTrace("TestTrackHits") << "last tp is associated with fraction=" << (tP.end()-1)->second;
    std::vector<unsigned int> tpids;
    for (TrackingParticle::g4t_iterator g4T = tp->g4Track_begin(); g4T != tp->g4Track_end(); ++g4T) {
      LogTrace("TestTrackHits") << "tp id=" << g4T->trackId();
      tpids.push_back(g4T->trackId());
    }

    //LogTrace("TestTrackHits") << "Analyzing hits of track number " << ++yyy << " good track number " << ++yy;
    int pp = 0;
    for (std::vector<TrajectoryMeasurement>::iterator tm = tmColl.begin(); tm != tmColl.end(); ++tm) {
      tchi2 += tm->estimate();

      LogTrace("TestTrackHits") << "+++++++++++++++++new hit+++++++++++++++++";
      CTTRHp rhit = tm->recHit();
      //TSOS state = tm->backwardPredictedState();
      //TSOS state = tm->forwardPredictedState();
      TSOS state = combiner(tm->backwardPredictedState(), tm->forwardPredictedState());

      if (rhit->isValid() == 0 && rhit->det() != nullptr)
        continue;
      evtHits++;
      LogTrace("TestTrackHits") << "valid hit #" << ++pp << "of hits=" << track->numberOfValidHits();

      int subdetId = rhit->det()->geographicalId().subdetId();
      DetId id = rhit->det()->geographicalId();
      int layerId = tTopo->layer(id);
      LogTrace("TestTrackHits") << "subdetId=" << subdetId << " layerId=" << layerId;

      const Surface* surf = rhit->surface();
      if (surf == nullptr)
        continue;

      double energyLoss_ = 0.;
      unsigned int monoId = 0;
      std::vector<double> energyLossM;
      std::vector<double> energyLossS;
      std::vector<PSimHit> assSimHits = hitAssociator.associateHit(*(rhit)->hit());
      unsigned int simhitvecsize = assSimHits.size();
      if (simhitvecsize == 0)
        continue;
      PSimHit shit;
      std::vector<unsigned int> trackIds;
      energyLossS.clear();
      energyLossM.clear();
      for (std::vector<PSimHit>::const_iterator m = assSimHits.begin(); m < assSimHits.end(); m++) {
        unsigned int tId = m->trackId();
        if (find(trackIds.begin(), trackIds.end(), tId) == trackIds.end())
          trackIds.push_back(tId);
        if (m->energyLoss() > energyLoss_) {
          shit = *m;
          energyLoss_ = m->energyLoss();
        }
        if (dynamic_cast<const SiStripMatchedRecHit2D*>(rhit->hit())) {
          if (monoId == 0)
            monoId = m->detUnitId();
          if (monoId == m->detUnitId()) {
            energyLossM.push_back(m->energyLoss());
          } else {
            energyLossS.push_back(m->energyLoss());
          }
          //std::cout << "detUnitId="  << m->detUnitId() << " trackId=" << m->trackId() << " energyLoss=" << m->energyLoss() << std::endl;
        } else {
          energyLossM.push_back(m->energyLoss());
        }
      }
      //double delta = 99999;
      //LocalPoint rhitLPv = rhit->localPosition();
      //vector<PSimHit> assSimHits = hitAssociator.associateHit(*(rhit)->hit());
      //unsigned int  simhitvecsize = assSimHits.size();
      //if (simhitvecsize==0) continue;
      //PSimHit shit;
      //for(std::vector<PSimHit>::const_iterator m=assSimHits.begin(); m<assSimHits.end(); m++){
      //  if ((m->localPosition()-rhitLPv).mag()<delta) {
      //    shit=*m;
      //    delta = (m->localPosition()-rhitLPv).mag();
      //  }
      //}

      //plot chi2 increment
      double chi2increment = tm->estimate();
      LogTrace("TestTrackHits") << "tm->estimate()=" << tm->estimate();
      title.str("");
      title << "Chi2Increment_" << subdetId << "-" << layerId;
      hChi2Increment[title.str()]->Fill(chi2increment);
      title.str("");
      title << "Chi2IncrementVsEta_" << subdetId << "-" << layerId;
      hChi2IncrementVsEta[title.str()]->Fill(track->eta(), chi2increment);
      hTotChi2Increment->Fill(chi2increment);
      hProcess_vs_Chi2->Fill(chi2increment, shit.processType());

      int clustersize = 0;
      bool mergedhit = false;
      if (dynamic_cast<const SiPixelRecHit*>(rhit->hit())) {
        clustersize = ((const SiPixelRecHit*)(rhit->hit()))->cluster()->size();
        hPixClsize_vs_Chi2->Fill(chi2increment, clustersize);
        hPixClusterSize->Fill(clustersize);
        hPixSimHitVecSize->Fill(simhitvecsize);
        if (simhitvecsize > 1)
          mergedhit = true;
      } else if (dynamic_cast<const SiStripRecHit2D*>(rhit->hit())) {
        clustersize = ((const SiStripRecHit2D*)(rhit->hit()))->cluster()->amplitudes().size();
        hSt1Clsize_vs_Chi2->Fill(chi2increment, clustersize);
        hSt1ClusterSize->Fill(clustersize);
        hSt1SimHitVecSize->Fill(simhitvecsize);
        if (simhitvecsize > 1)
          mergedhit = true;
      } else if (dynamic_cast<const SiStripMatchedRecHit2D*>(rhit->hit())) {
        int clsize1 = ((const SiStripMatchedRecHit2D*)(rhit->hit()))->monoCluster().amplitudes().size();
        int clsize2 = ((const SiStripMatchedRecHit2D*)(rhit->hit()))->stereoCluster().amplitudes().size();
        if (clsize1 > clsize2)
          clustersize = clsize1;
        else
          clustersize = clsize2;
        hSt2Clsize_vs_Chi2->Fill(chi2increment, clustersize);
        hSt2ClusterSize->Fill(clustersize);
        hSt2SimHitVecSize->Fill(simhitvecsize);
        if (simhitvecsize > 2)
          mergedhit = true;
      } else if (dynamic_cast<const ProjectedSiStripRecHit2D*>(rhit->hit())) {
        clustersize = ((const ProjectedSiStripRecHit2D*)(rhit->hit()))->originalHit().cluster()->amplitudes().size();
        hPrjClsize_vs_Chi2->Fill(chi2increment, clustersize);
        hPrjClusterSize->Fill(clustersize);
        hPrjSimHitVecSize->Fill(simhitvecsize);
        if (simhitvecsize > 1)
          mergedhit = true;
      }
      hClsize_vs_Chi2->Fill(chi2increment, clustersize);
      hClusterSize->Fill(clustersize);
      hSimHitVecSize->Fill(simhitvecsize);

      //       if (dynamic_cast<const SiPixelRecHit*>(rhit->hit()))
      //    hClsize_vs_Chi2->Fill( chi2increment, ((const SiPixelRecHit*)(rhit->hit()))->cluster()->size() );
      //       if (dynamic_cast<const SiStripRecHit2D*>(rhit->hit()))
      //    hClsize_vs_Chi2->Fill( chi2increment, ((const SiStripRecHit2D*)(rhit->hit()))->cluster()->amplitudes().size() );

      std::vector<SimHitIdpr> simTrackIds = hitAssociator.associateHitId(*(rhit)->hit());
      bool goodhit = false;
      for (size_t j = 0; j < simTrackIds.size(); j++) {
        LogTrace("TestTrackHits") << "hit id=" << simTrackIds[j].first;
        for (size_t jj = 0; jj < tpids.size(); jj++) {
          if (simTrackIds[j].first == tpids[jj])
            goodhit = true;
          break;
        }
      }
      bool shared = true;
      bool ioniOnly = true;
      const SiPixelRecHit* pix = dynamic_cast<const SiPixelRecHit*>(rhit->hit());
      if (goodhit) {
        if (energyLossM.size() > 1 && energyLossS.size() <= 1) {
          sort(energyLossM.begin(), energyLossM.end(), greater<double>());
          energyLossRatio->Fill(energyLossM[1] / energyLossM[0]);
        } else if (energyLossS.size() > 1 && energyLossM.size() <= 1) {
          sort(energyLossS.begin(), energyLossS.end(), greater<double>());
          energyLossRatio->Fill(energyLossS[1] / energyLossS[0]);
        } else if (energyLossS.size() > 1 && energyLossM.size() > 1) {
          sort(energyLossM.begin(), energyLossM.end(), greater<double>());
          sort(energyLossS.begin(), energyLossS.end(), greater<double>());
          energyLossM[1] / energyLossM[0] > energyLossS[1] / energyLossS[0]
              ? energyLossRatio->Fill(energyLossM[1] / energyLossM[0])
              : energyLossRatio->Fill(energyLossS[1] / energyLossS[0]);
        }
        if (mergedhit) {
          //not optimized for matched hits
          LogVerbatim("TestTrackHits") << "MERGED HIT" << std::endl;
          unsigned int idc = 0;
          for (size_t jj = 0; jj < tpids.size(); jj++) {
            idc += std::count(trackIds.begin(), trackIds.end(), tpids[jj]);
          }
          if (idc == trackIds.size()) {
            shared = false;
          }
          for (std::vector<PSimHit>::const_iterator m = assSimHits.begin() + 1; m < assSimHits.end(); m++) {
            if ((m->processType() != 7 && m->processType() != 8 && m->processType() != 9) &&
                abs(m->particleType()) != 11) {
              ioniOnly = false;
              break;
            }
          }
          if (ioniOnly && !shared) {
            title.str("");
            title << "Chi2DeltaHit_" << subdetId << "-" << layerId;
            hChi2DeltaHit[title.str()]->Fill(chi2increment);
            hTotChi2DeltaHit->Fill(chi2increment);
            if (pix) {
              probXdelta->Fill(pix->probabilityX());
              probYdelta->Fill(pix->probabilityY());
            }
          } else if (!ioniOnly && !shared) {
            title.str("");
            title << "Chi2NSharedHit_" << subdetId << "-" << layerId;
            hChi2NSharedHit[title.str()]->Fill(chi2increment);
            hTotChi2NSharedHit->Fill(chi2increment);
            if (pix) {
              probXnoshare->Fill(pix->probabilityX());
              probYnoshare->Fill(pix->probabilityY());
            }
          } else {
            title.str("");
            title << "Chi2SharedHit_" << subdetId << "-" << layerId;
            hChi2SharedHit[title.str()]->Fill(chi2increment);
            hTotChi2SharedHit->Fill(chi2increment);
            if (pix) {
              probXshared->Fill(pix->probabilityX());
              probYshared->Fill(pix->probabilityY());
            }
          }

          for (std::vector<PSimHit>::const_iterator m = assSimHits.begin(); m < assSimHits.end(); m++) {
            unsigned int tId = m->trackId();
            LogVerbatim("TestTrackHits") << "component with id=" << tId << " eLoss=" << m->energyLoss()
                                         << " pType=" << m->processType();
            if (find(tpids.begin(), tpids.end(), tId) == tpids.end())
              continue;
            if (m->processType() == 2) {
              GlobalPoint gpr = rhit->globalPosition();
              AlgebraicSymMatrix33 ger = rhit->globalPositionError().matrix();
              GlobalPoint gps = surf->toGlobal(m->localPosition());
              LogVerbatim("TestTrackHits") << gpr << " " << gps << " " << ger;
              ROOT::Math::SVector<double, 3> delta;
              delta[0] = gpr.x() - gps.x();
              delta[1] = gpr.y() - gps.y();
              delta[2] = gpr.z() - gps.z();
              LogVerbatim("TestTrackHits") << delta << " " << ger;
              double mpull = sqrt(delta[0] * delta[0] / ger[0][0] + delta[1] * delta[1] / ger[1][1] +
                                  delta[2] * delta[2] / ger[2][2]);
              LogVerbatim("TestTrackHits") << "hit pull=" << mpull;  //ger.similarity(delta);
              mergedPull->Fill(mpull);
              break;
            }
          }
        } else {
          LogVerbatim("TestTrackHits") << "good hit";
          title.str("");
          title << "Chi2GoodHit_" << subdetId << "-" << layerId;
          hChi2GoodHit[title.str()]->Fill(chi2increment);
          hTotChi2GoodHit->Fill(chi2increment);
          if (pix) {
            probXgood->Fill(pix->probabilityX());
            probYgood->Fill(pix->probabilityY());
          }
        }
      } else {
        LogVerbatim("TestTrackHits") << "bad hit";
        title.str("");
        title << "Chi2BadHit_" << subdetId << "-" << layerId;
        hChi2BadHit[title.str()]->Fill(chi2increment);
        hTotChi2BadHit->Fill(chi2increment);
        goodbadmerged->Fill(2);
        if (pix) {
          probXbad->Fill(pix->probabilityX());
          probYbad->Fill(pix->probabilityY());
        }
      }
      hGoodHit_vs_Chi2->Fill(chi2increment, goodhit);

      LocalVector shitLMom;
      LocalPoint shitLPos;
      if (dynamic_cast<const SiStripMatchedRecHit2D*>(rhit->hit())) {
        if (simhitvecsize > 2 && goodhit) {
          if (ioniOnly && !shared)
            goodbadmerged->Fill(3);
          else if (!ioniOnly && !shared)
            goodbadmerged->Fill(4);
          else
            goodbadmerged->Fill(5);
        } else if (goodhit)
          goodbadmerged->Fill(1);
        double rechitmatchedx = rhit->localPosition().x();
        double rechitmatchedy = rhit->localPosition().y();
        double mindist = 999999;
        float distx, disty;
        std::pair<LocalPoint, LocalVector> closestPair;
        const StripGeomDetUnit* stripDet = (StripGeomDetUnit*)((const GluedGeomDet*)(rhit)->det())->stereoDet();
        const BoundPlane& plane = (rhit)->det()->surface();
        for (std::vector<PSimHit>::const_iterator m = assSimHits.begin(); m < assSimHits.end(); m++) {
          //project simhit;
          std::pair<LocalPoint, LocalVector> hitPair = projectHit((*m), stripDet, plane);
          distx = fabs(rechitmatchedx - hitPair.first.x());
          disty = fabs(rechitmatchedy - hitPair.first.y());
          double dist = distx * distx + disty * disty;
          if (sqrt(dist) < mindist) {
            mindist = dist;
            closestPair = hitPair;
          }
        }
        shitLPos = closestPair.first;
        shitLMom = closestPair.second;
      } else {
        if (simhitvecsize > 1 && goodhit) {
          if (ioniOnly && !shared)
            goodbadmerged->Fill(3);
          else if (!ioniOnly && !shared)
            goodbadmerged->Fill(4);
          else
            goodbadmerged->Fill(5);
        } else if (goodhit)
          goodbadmerged->Fill(1);
        shitLPos = shit.localPosition();
        shitLMom = shit.momentumAtEntry();
      }
      GlobalVector shitGMom = surf->toGlobal(shitLMom);
      GlobalPoint shitGPos = surf->toGlobal(shitLPos);

      GlobalVector tsosGMom = state.globalMomentum();
      GlobalError tsosGMEr(state.cartesianError().matrix().Sub<AlgebraicSymMatrix33>(3, 3));
      GlobalPoint tsosGPos = state.globalPosition();
      GlobalError tsosGPEr = state.cartesianError().position();

      GlobalPoint rhitGPos = (rhit)->globalPosition();
      GlobalError rhitGPEr = (rhit)->globalPositionError();

      LogVerbatim("TestTrackHits") << "assSimHits.size()=" << assSimHits.size();
      LogVerbatim("TestTrackHits") << "tsos globalPos   =" << tsosGPos;
      LogVerbatim("TestTrackHits") << "sim hit globalPos=" << shitGPos;
      LogVerbatim("TestTrackHits") << "rec hit globalPos=" << rhitGPos;
      LogVerbatim("TestTrackHits") << "geographicalId   =" << rhit->det()->geographicalId().rawId();
      LogVerbatim("TestTrackHits") << "surface position =" << surf->position();

#if 0
      if (rhit->detUnit()) LogTrace("TestTrackHits") << "rhit->detUnit()->geographicalId()=" 
                             << rhit->detUnit()->geographicalId().rawId() ;
      LogTrace("TestTrackHits") << "rhit->det()->surface().position()=" 
                << rhit->det()->surface().position() ;
      if (rhit->detUnit()) LogTrace("TestTrackHits") << "rhit->detUnit()->surface().position()="  
                             << rhit->detUnit()->surface().position() ;
      LogTrace("TestTrackHits") << "rhit->det()->position()=" << rhit->det()->position() ;
      if (rhit->detUnit()) LogTrace("TestTrackHits") << "rhit->detUnit()->position()="  
                             << rhit->detUnit()->position() ;
      LogTrace("TestTrackHits") << "rhit->det()->surface().bounds().length()=" 
                << rhit->det()->surface().bounds().length() ;
      if (rhit->detUnit()) LogTrace("TestTrackHits") << "rhit->detUnit()->surface().bounds().length()="  
                             << rhit->detUnit()->surface().bounds().length() ;
      LogTrace("TestTrackHits") << "rhit->det()->surface().bounds().width()=" 
                << rhit->det()->surface().bounds().width() ;
      if (rhit->detUnit()) LogTrace("TestTrackHits") << "rhit->detUnit()->surface().bounds().width()="  
                             << rhit->detUnit()->surface().bounds().width() ;
      LogTrace("TestTrackHits") << "rhit->det()->surface().bounds().thickness()=" 
                << rhit->det()->surface().bounds().thickness() ;
      if (rhit->detUnit()) LogTrace("TestTrackHits") << "rhit->detUnit()->surface().bounds().thickness()="  
                             << rhit->detUnit()->surface().bounds().thickness() ;
#endif

      double pullGPX_rs = (rhitGPos.x() - shitGPos.x()) / sqrt(rhitGPEr.cxx());
      double pullGPY_rs = (rhitGPos.y() - shitGPos.y()) / sqrt(rhitGPEr.cyy());
      double pullGPZ_rs = (rhitGPos.z() - shitGPos.z()) / sqrt(rhitGPEr.czz());
      //double pullGPX_rs = (rhitGPos.x()-shitGPos.x());
      //double pullGPY_rs = (rhitGPos.y()-shitGPos.y());
      //double pullGPZ_rs = (rhitGPos.z()-shitGPos.z());

      LogTrace("TestTrackHits") << "rs";

      title.str("");
      title << "PullGP_X_" << subdetId << "-" << layerId << "_rs";
      hPullGP_X_rs[title.str()]->Fill(pullGPX_rs);
      title.str("");
      title << "PullGP_Y_" << subdetId << "-" << layerId << "_rs";
      hPullGP_Y_rs[title.str()]->Fill(pullGPY_rs);
      title.str("");
      title << "PullGP_Z_" << subdetId << "-" << layerId << "_rs";
      hPullGP_Z_rs[title.str()]->Fill(pullGPZ_rs);

      double pullGPX_tr = (tsosGPos.x() - rhitGPos.x()) / sqrt(tsosGPEr.cxx() + rhitGPEr.cxx());
      double pullGPY_tr = (tsosGPos.y() - rhitGPos.y()) / sqrt(tsosGPEr.cyy() + rhitGPEr.cyy());
      double pullGPZ_tr = (tsosGPos.z() - rhitGPos.z()) / sqrt(tsosGPEr.czz() + rhitGPEr.czz());
      //double pullGPX_tr = (tsosGPos.x()-rhitGPos.x());
      //double pullGPY_tr = (tsosGPos.y()-rhitGPos.y());
      //double pullGPZ_tr = (tsosGPos.z()-rhitGPos.z());

      LogTrace("TestTrackHits") << "tr";

      title.str("");
      title << "PullGP_X_" << subdetId << "-" << layerId << "_tr";
      hPullGP_X_tr[title.str()]->Fill(pullGPX_tr);
      title.str("");
      title << "PullGP_Y_" << subdetId << "-" << layerId << "_tr";
      hPullGP_Y_tr[title.str()]->Fill(pullGPY_tr);
      title.str("");
      title << "PullGP_Z_" << subdetId << "-" << layerId << "_tr";
      hPullGP_Z_tr[title.str()]->Fill(pullGPZ_tr);

      double pullGPX_ts = (tsosGPos.x() - shitGPos.x()) / sqrt(tsosGPEr.cxx());
      double pullGPY_ts = (tsosGPos.y() - shitGPos.y()) / sqrt(tsosGPEr.cyy());
      double pullGPZ_ts = (tsosGPos.z() - shitGPos.z()) / sqrt(tsosGPEr.czz());
      //double pullGPX_ts = (tsosGPos.x()-shitGPos.x());
      //double pullGPY_ts = (tsosGPos.y()-shitGPos.y());
      //double pullGPZ_ts = (tsosGPos.z()-shitGPos.z());

      LogTrace("TestTrackHits") << "ts1";

      title.str("");
      title << "PullGP_X_" << subdetId << "-" << layerId << "_ts";
      hPullGP_X_ts[title.str()]->Fill(pullGPX_ts);
      title.str("");
      title << "PullGP_Y_" << subdetId << "-" << layerId << "_ts";
      hPullGP_Y_ts[title.str()]->Fill(pullGPY_ts);
      title.str("");
      title << "PullGP_Z_" << subdetId << "-" << layerId << "_ts";
      hPullGP_Z_ts[title.str()]->Fill(pullGPZ_ts);

      double pullGMX_ts = (tsosGMom.x() - shitGMom.x()) / sqrt(tsosGMEr.cxx());
      double pullGMY_ts = (tsosGMom.y() - shitGMom.y()) / sqrt(tsosGMEr.cyy());
      double pullGMZ_ts = (tsosGMom.z() - shitGMom.z()) / sqrt(tsosGMEr.czz());
      //double pullGMX_ts = (tsosGMom.x()-shitGMom.x());
      //double pullGMY_ts = (tsosGMom.y()-shitGMom.y());
      //double pullGMZ_ts = (tsosGMom.z()-shitGMom.z());

      LogTrace("TestTrackHits") << "ts2";

      title.str("");
      title << "PullGM_X_" << subdetId << "-" << layerId << "_ts";
      hPullGM_X_ts[title.str()]->Fill(pullGMX_ts);
      title.str("");
      title << "PullGM_Y_" << subdetId << "-" << layerId << "_ts";
      hPullGM_Y_ts[title.str()]->Fill(pullGMY_ts);
      title.str("");
      title << "PullGM_Z_" << subdetId << "-" << layerId << "_ts";
      hPullGM_Z_ts[title.str()]->Fill(pullGMZ_ts);
      if (dynamic_cast<const SiStripMatchedRecHit2D*>(rhit->hit())) {
        Propagator* thePropagatorAnyDir = new PropagatorWithMaterial(anyDirection, 0.105, theMF.product(), 1.6);
        //mono
        LogTrace("TestTrackHits") << "MONO HIT";
        auto m = dynamic_cast<const SiStripMatchedRecHit2D*>(rhit->hit())->monoHit();
        CTTRHp tMonoHit = theBuilder->build(&m);
        if (tMonoHit == nullptr)
          continue;
        vector<PSimHit> assMonoSimHits = hitAssociator.associateHit(*tMonoHit->hit());
        if (assMonoSimHits.empty())
          continue;
        const PSimHit sMonoHit = *(assMonoSimHits.begin());
        const Surface* monoSurf = &(tMonoHit->det()->surface());
        if (monoSurf == nullptr)
          continue;
        TSOS monoState = thePropagatorAnyDir->propagate(state, *monoSurf);
        if (monoState.isValid() == 0)
          continue;

        LocalVector monoShitLMom = sMonoHit.momentumAtEntry();
        GlobalVector monoShitGMom = monoSurf->toGlobal(monoShitLMom);
        LocalPoint monoShitLPos = sMonoHit.localPosition();
        GlobalPoint monoShitGPos = monoSurf->toGlobal(monoShitLPos);

        //LogTrace("TestTrackHits") << "assMonoSimHits.size()=" << assMonoSimHits.size() ;
        //LogTrace("TestTrackHits") << "mono shit=" << monoShitGPos ;

        GlobalVector monoTsosGMom = monoState.globalMomentum();
        GlobalError monoTsosGMEr(monoState.cartesianError().matrix().Sub<AlgebraicSymMatrix33>(3, 3));
        GlobalPoint monoTsosGPos = monoState.globalPosition();
        GlobalError monoTsosGPEr = monoState.cartesianError().position();

        GlobalPoint monoRhitGPos = tMonoHit->globalPosition();
        GlobalError monoRhitGPEr = tMonoHit->globalPositionError();

        double pullGPX_rs_mono = (monoRhitGPos.x() - monoShitGPos.x()) / sqrt(monoRhitGPEr.cxx());
        double pullGPY_rs_mono = (monoRhitGPos.y() - monoShitGPos.y()) / sqrt(monoRhitGPEr.cyy());
        double pullGPZ_rs_mono = (monoRhitGPos.z() - monoShitGPos.z()) / sqrt(monoRhitGPEr.czz());
        //double pullGPX_rs_mono = (monoRhitGPos.x()-monoShitGPos.x());
        //double pullGPY_rs_mono = (monoRhitGPos.y()-monoShitGPos.y());
        //double pullGPZ_rs_mono = (monoRhitGPos.z()-monoShitGPos.z());

        MeasurementExtractor meMo(monoState);
        double chi2mono = computeChi2Increment(meMo, tMonoHit);

        title.str("");
        title << "Chi2Increment_mono_" << subdetId << "-" << layerId;
        hChi2Increment_mono[title.str()]->Fill(chi2mono);

        title.str("");
        title << "PullGP_X_" << subdetId << "-" << layerId << "_rs_mono";
        hPullGP_X_rs_mono[title.str()]->Fill(pullGPX_rs_mono);
        title.str("");
        title << "PullGP_Y_" << subdetId << "-" << layerId << "_rs_mono";
        hPullGP_Y_rs_mono[title.str()]->Fill(pullGPY_rs_mono);
        title.str("");
        title << "PullGP_Z_" << subdetId << "-" << layerId << "_rs_mono";
        hPullGP_Z_rs_mono[title.str()]->Fill(pullGPZ_rs_mono);

        double pullGPX_tr_mono = (monoTsosGPos.x() - monoRhitGPos.x()) / sqrt(monoTsosGPEr.cxx() + monoRhitGPEr.cxx());
        double pullGPY_tr_mono = (monoTsosGPos.y() - monoRhitGPos.y()) / sqrt(monoTsosGPEr.cyy() + monoRhitGPEr.cyy());
        double pullGPZ_tr_mono = (monoTsosGPos.z() - monoRhitGPos.z()) / sqrt(monoTsosGPEr.czz() + monoRhitGPEr.czz());
        //double pullGPX_tr_mono = (monoTsosGPos.x()-monoRhitGPos.x());
        //double pullGPY_tr_mono = (monoTsosGPos.y()-monoRhitGPos.y());
        //double pullGPZ_tr_mono = (monoTsosGPos.z()-monoRhitGPos.z());

        title.str("");
        title << "PullGP_X_" << subdetId << "-" << layerId << "_tr_mono";
        hPullGP_X_tr_mono[title.str()]->Fill(pullGPX_tr_mono);
        title.str("");
        title << "PullGP_Y_" << subdetId << "-" << layerId << "_tr_mono";
        hPullGP_Y_tr_mono[title.str()]->Fill(pullGPY_tr_mono);
        title.str("");
        title << "PullGP_Z_" << subdetId << "-" << layerId << "_tr_mono";
        hPullGP_Z_tr_mono[title.str()]->Fill(pullGPZ_tr_mono);

        double pullGPX_ts_mono = (monoTsosGPos.x() - monoShitGPos.x()) / sqrt(monoTsosGPEr.cxx());
        double pullGPY_ts_mono = (monoTsosGPos.y() - monoShitGPos.y()) / sqrt(monoTsosGPEr.cyy());
        double pullGPZ_ts_mono = (monoTsosGPos.z() - monoShitGPos.z()) / sqrt(monoTsosGPEr.czz());
        //double pullGPX_ts_mono = (monoTsosGPos.x()-monoShitGPos.x());
        //double pullGPY_ts_mono = (monoTsosGPos.y()-monoShitGPos.y());
        //double pullGPZ_ts_mono = (monoTsosGPos.z()-monoShitGPos.z());

        title.str("");
        title << "PullGP_X_" << subdetId << "-" << layerId << "_ts_mono";
        hPullGP_X_ts_mono[title.str()]->Fill(pullGPX_ts_mono);
        title.str("");
        title << "PullGP_Y_" << subdetId << "-" << layerId << "_ts_mono";
        hPullGP_Y_ts_mono[title.str()]->Fill(pullGPY_ts_mono);
        title.str("");
        title << "PullGP_Z_" << subdetId << "-" << layerId << "_ts_mono";
        hPullGP_Z_ts_mono[title.str()]->Fill(pullGPZ_ts_mono);

        double pullGMX_ts_mono = (monoTsosGMom.x() - monoShitGMom.x()) / sqrt(monoTsosGMEr.cxx());
        double pullGMY_ts_mono = (monoTsosGMom.y() - monoShitGMom.y()) / sqrt(monoTsosGMEr.cyy());
        double pullGMZ_ts_mono = (monoTsosGMom.z() - monoShitGMom.z()) / sqrt(monoTsosGMEr.czz());
        //double pullGMX_ts_mono = (monoTsosGMom.x()-monoShitGMom.x());
        //double pullGMY_ts_mono = (monoTsosGMom.y()-monoShitGMom.y());
        //double pullGMZ_ts_mono = (monoTsosGMom.z()-monoShitGMom.z());

        title.str("");
        title << "PullGM_X_" << subdetId << "-" << layerId << "_ts_mono";
        hPullGM_X_ts_mono[title.str()]->Fill(pullGMX_ts_mono);
        title.str("");
        title << "PullGM_Y_" << subdetId << "-" << layerId << "_ts_mono";
        hPullGM_Y_ts_mono[title.str()]->Fill(pullGMY_ts_mono);
        title.str("");
        title << "PullGM_Z_" << subdetId << "-" << layerId << "_ts_mono";
        hPullGM_Z_ts_mono[title.str()]->Fill(pullGMZ_ts_mono);

        //stereo
        LogTrace("TestTrackHits") << "STEREO HIT";
        auto s = dynamic_cast<const SiStripMatchedRecHit2D*>(rhit->hit())->stereoHit();
        CTTRHp tStereoHit = theBuilder->build(&s);
        if (tStereoHit == nullptr)
          continue;
        vector<PSimHit> assStereoSimHits = hitAssociator.associateHit(*tStereoHit->hit());
        if (assStereoSimHits.empty())
          continue;
        const PSimHit sStereoHit = *(assStereoSimHits.begin());
        const Surface* stereoSurf = &(tStereoHit->det()->surface());
        if (stereoSurf == nullptr)
          continue;
        TSOS stereoState = thePropagatorAnyDir->propagate(state, *stereoSurf);
        if (stereoState.isValid() == 0)
          continue;

        LocalVector stereoShitLMom = sStereoHit.momentumAtEntry();
        GlobalVector stereoShitGMom = stereoSurf->toGlobal(stereoShitLMom);
        LocalPoint stereoShitLPos = sStereoHit.localPosition();
        GlobalPoint stereoShitGPos = stereoSurf->toGlobal(stereoShitLPos);

        //LogTrace("TestTrackHits") << "assStereoSimHits.size()=" << assStereoSimHits.size() ;
        //LogTrace("TestTrackHits") << "stereo shit=" << stereoShitGPos ;

        GlobalVector stereoTsosGMom = stereoState.globalMomentum();
        GlobalError stereoTsosGMEr(stereoState.cartesianError().matrix().Sub<AlgebraicSymMatrix33>(3, 3));
        GlobalPoint stereoTsosGPos = stereoState.globalPosition();
        GlobalError stereoTsosGPEr = stereoState.cartesianError().position();

        GlobalPoint stereoRhitGPos = tStereoHit->globalPosition();
        GlobalError stereoRhitGPEr = tStereoHit->globalPositionError();

        MeasurementExtractor meSt(stereoState);
        double chi2stereo = computeChi2Increment(meSt, tStereoHit);

        title.str("");
        title << "Chi2Increment_stereo_" << subdetId << "-" << layerId;
        hChi2Increment_stereo[title.str()]->Fill(chi2stereo);

        double pullGPX_rs_stereo = (stereoRhitGPos.x() - stereoShitGPos.x()) / sqrt(stereoRhitGPEr.cxx());
        double pullGPY_rs_stereo = (stereoRhitGPos.y() - stereoShitGPos.y()) / sqrt(stereoRhitGPEr.cyy());
        double pullGPZ_rs_stereo = (stereoRhitGPos.z() - stereoShitGPos.z()) / sqrt(stereoRhitGPEr.czz());
        //  double pullGPX_rs_stereo = (stereoRhitGPos.x()-stereoShitGPos.x())/*/sqrt(stereoRhitGPEr.cxx())*/;
        //  double pullGPY_rs_stereo = (stereoRhitGPos.y()-stereoShitGPos.y())/*/sqrt(stereoRhitGPEr.cyy())*/;
        //  double pullGPZ_rs_stereo = (stereoRhitGPos.z()-stereoShitGPos.z())/*/sqrt(stereoRhitGPEr.czz())*/;

        title.str("");
        title << "PullGP_X_" << subdetId << "-" << layerId << "_rs_stereo";
        hPullGP_X_rs_stereo[title.str()]->Fill(pullGPX_rs_stereo);
        title.str("");
        title << "PullGP_Y_" << subdetId << "-" << layerId << "_rs_stereo";
        hPullGP_Y_rs_stereo[title.str()]->Fill(pullGPY_rs_stereo);
        title.str("");
        title << "PullGP_Z_" << subdetId << "-" << layerId << "_rs_stereo";
        hPullGP_Z_rs_stereo[title.str()]->Fill(pullGPZ_rs_stereo);

        double pullGPX_tr_stereo =
            (stereoTsosGPos.x() - stereoRhitGPos.x()) / sqrt(stereoTsosGPEr.cxx() + stereoRhitGPEr.cxx());
        double pullGPY_tr_stereo =
            (stereoTsosGPos.y() - stereoRhitGPos.y()) / sqrt(stereoTsosGPEr.cyy() + stereoRhitGPEr.cyy());
        double pullGPZ_tr_stereo =
            (stereoTsosGPos.z() - stereoRhitGPos.z()) / sqrt(stereoTsosGPEr.czz() + stereoRhitGPEr.czz());
        //double pullGPX_tr_stereo = (stereoTsosGPos.x()-stereoRhitGPos.x());
        //double pullGPY_tr_stereo = (stereoTsosGPos.y()-stereoRhitGPos.y());
        //double pullGPZ_tr_stereo = (stereoTsosGPos.z()-stereoRhitGPos.z());

        title.str("");
        title << "PullGP_X_" << subdetId << "-" << layerId << "_tr_stereo";
        hPullGP_X_tr_stereo[title.str()]->Fill(pullGPX_tr_stereo);
        title.str("");
        title << "PullGP_Y_" << subdetId << "-" << layerId << "_tr_stereo";
        hPullGP_Y_tr_stereo[title.str()]->Fill(pullGPY_tr_stereo);
        title.str("");
        title << "PullGP_Z_" << subdetId << "-" << layerId << "_tr_stereo";
        hPullGP_Z_tr_stereo[title.str()]->Fill(pullGPZ_tr_stereo);

        double pullGPX_ts_stereo = (stereoTsosGPos.x() - stereoShitGPos.x()) / sqrt(stereoTsosGPEr.cxx());
        double pullGPY_ts_stereo = (stereoTsosGPos.y() - stereoShitGPos.y()) / sqrt(stereoTsosGPEr.cyy());
        double pullGPZ_ts_stereo = (stereoTsosGPos.z() - stereoShitGPos.z()) / sqrt(stereoTsosGPEr.czz());
        //double pullGPX_ts_stereo = (stereoTsosGPos.x()-stereoShitGPos.x());
        //double pullGPY_ts_stereo = (stereoTsosGPos.y()-stereoShitGPos.y());
        //double pullGPZ_ts_stereo = (stereoTsosGPos.z()-stereoShitGPos.z());

        title.str("");
        title << "PullGP_X_" << subdetId << "-" << layerId << "_ts_stereo";
        hPullGP_X_ts_stereo[title.str()]->Fill(pullGPX_ts_stereo);
        title.str("");
        title << "PullGP_Y_" << subdetId << "-" << layerId << "_ts_stereo";
        hPullGP_Y_ts_stereo[title.str()]->Fill(pullGPY_ts_stereo);
        title.str("");
        title << "PullGP_Z_" << subdetId << "-" << layerId << "_ts_stereo";
        hPullGP_Z_ts_stereo[title.str()]->Fill(pullGPZ_ts_stereo);

        double pullGMX_ts_stereo = (stereoTsosGMom.x() - stereoShitGMom.x()) / sqrt(stereoTsosGMEr.cxx());
        double pullGMY_ts_stereo = (stereoTsosGMom.y() - stereoShitGMom.y()) / sqrt(stereoTsosGMEr.cyy());
        double pullGMZ_ts_stereo = (stereoTsosGMom.z() - stereoShitGMom.z()) / sqrt(stereoTsosGMEr.czz());
        //double pullGMX_ts_stereo = (stereoTsosGMom.x()-stereoShitGMom.x());
        //double pullGMY_ts_stereo = (stereoTsosGMom.y()-stereoShitGMom.y());
        //double pullGMZ_ts_stereo = (stereoTsosGMom.z()-stereoShitGMom.z());

        title.str("");
        title << "PullGM_X_" << subdetId << "-" << layerId << "_ts_stereo";
        hPullGM_X_ts_stereo[title.str()]->Fill(pullGMX_ts_stereo);
        title.str("");
        title << "PullGM_Y_" << subdetId << "-" << layerId << "_ts_stereo";
        hPullGM_Y_ts_stereo[title.str()]->Fill(pullGMY_ts_stereo);
        title.str("");
        title << "PullGM_Z_" << subdetId << "-" << layerId << "_ts_stereo";
        hPullGM_Z_ts_stereo[title.str()]->Fill(pullGMZ_ts_stereo);
      }
    }
    LogTrace("TestTrackHits") << "traj chi2=" << tchi2;
    LogTrace("TestTrackHits") << "track chi2=" << track->chi2();
    i++;
  }
  LogTrace("TestTrackHits") << "end of event: processd hits=" << evtHits;
}

beginRun()

void TestTrackHits::beginRun ( edm::Run const &  run, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 44 of file TestTrackHits.cc.

References energyLossRatio, file, edm::EventSetup::getHandle(), goodbadmerged, hChi2BadHit, hChi2DeltaHit, hChi2GoodHit, hChi2Increment, hChi2Increment_mono, hChi2Increment_stereo, hChi2IncrementVsEta, hChi2NSharedHit, hChi2SharedHit, hClsize_vs_Chi2, hClusterSize, hGoodHit_vs_Chi2, hPixClsize_vs_Chi2, hPixClusterSize, hPixSimHitVecSize, hPrjClsize_vs_Chi2, hPrjClusterSize, hPrjSimHitVecSize, hProcess_vs_Chi2, hPullGM_X_ts, hPullGM_X_ts_mono, hPullGM_X_ts_stereo, hPullGM_Y_ts, hPullGM_Y_ts_mono, hPullGM_Y_ts_stereo, hPullGM_Z_ts, hPullGM_Z_ts_mono, hPullGM_Z_ts_stereo, hPullGP_X_rs, hPullGP_X_rs_mono, hPullGP_X_rs_stereo, hPullGP_X_tr, hPullGP_X_tr_mono, hPullGP_X_tr_stereo, hPullGP_X_ts, hPullGP_X_ts_mono, hPullGP_X_ts_stereo, hPullGP_Y_rs, hPullGP_Y_rs_mono, hPullGP_Y_rs_stereo, hPullGP_Y_tr, hPullGP_Y_tr_mono, hPullGP_Y_tr_stereo, hPullGP_Y_ts, hPullGP_Y_ts_mono, hPullGP_Y_ts_stereo, hPullGP_Z_rs, hPullGP_Z_rs_mono, hPullGP_Z_rs_stereo, hPullGP_Z_tr, hPullGP_Z_tr_mono, hPullGP_Z_tr_stereo, hPullGP_Z_ts, hPullGP_Z_ts_mono, hPullGP_Z_ts_stereo, hSimHitVecSize, hSt1Clsize_vs_Chi2, hSt1ClusterSize, hSt1SimHitVecSize, hSt2Clsize_vs_Chi2, hSt2ClusterSize, hSt2SimHitVecSize, hTotChi2BadHit, hTotChi2DeltaHit, hTotChi2GoodHit, hTotChi2Increment, hTotChi2NSharedHit, hTotChi2SharedHit, mps_fire::i, dqmiolumiharvest::j, mergedPull, out, probXbad, probXdelta, probXgood, probXnoshare, probXshared, probYbad, probYdelta, probYgood, probYnoshare, probYshared, theBuilder, theBuilderToken, theG, theGToken, theMF, theMFToken, thePropagator, thePropagatorToken, theUpdator, theUpdatorToken, and title.

                                                                           {
  theG = iSetup.getHandle(theGToken);
  theMF = iSetup.getHandle(theMFToken);
  thePropagator = iSetup.getHandle(thePropagatorToken);
  theBuilder = iSetup.getHandle(theBuilderToken);
  theUpdator = iSetup.getHandle(theUpdatorToken);

  file = new TFile(out.c_str(), "recreate");
  for (int i = 0; i != 6; i++)
    for (int j = 0; j != 9; j++) {
      if (i == 0 && j > 2)
        break;
      if (i == 1 && j > 1)
        break;
      if (i == 2 && j > 3)
        break;
      if (i == 3 && j > 2)
        break;
      if (i == 4 && j > 5)
        break;
      if (i == 5 && j > 8)
        break;

      title.str("");
      title << "Chi2Increment_" << i + 1 << "-" << j + 1;
      hChi2Increment[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, 0, 100);
      title.str("");
      title << "Chi2IncrementVsEta_" << i + 1 << "-" << j + 1;
      hChi2IncrementVsEta[title.str()] =
          new TH2F(title.str().c_str(), title.str().c_str(), 50, -2.5, 2.5, 1000, 0, 100);
      title.str("");
      title << "Chi2GoodHit_" << i + 1 << "-" << j + 1;
      hChi2GoodHit[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, 0, 100);
      title.str("");
      title << "Chi2BadHit_" << i + 1 << "-" << j + 1;
      hChi2BadHit[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, 0, 100);
      title.str("");
      title << "Chi2DeltaHit_" << i + 1 << "-" << j + 1;
      hChi2DeltaHit[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, 0, 100);
      title.str("");
      title << "Chi2NSharedHit_" << i + 1 << "-" << j + 1;
      hChi2NSharedHit[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, 0, 100);
      title.str("");
      title << "Chi2SharedHit_" << i + 1 << "-" << j + 1;
      hChi2SharedHit[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, 0, 100);

      title.str("");
      title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_ts";
      hPullGP_X_ts[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
      title.str("");
      title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_ts";
      hPullGP_Y_ts[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
      title.str("");
      title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_ts";
      hPullGP_Z_ts[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);

      title.str("");
      title << "PullGM_X_" << i + 1 << "-" << j + 1 << "_ts";
      hPullGM_X_ts[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
      title.str("");
      title << "PullGM_Y_" << i + 1 << "-" << j + 1 << "_ts";
      hPullGM_Y_ts[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
      title.str("");
      title << "PullGM_Z_" << i + 1 << "-" << j + 1 << "_ts";
      hPullGM_Z_ts[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);

      title.str("");
      title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_tr";
      hPullGP_X_tr[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
      title.str("");
      title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_tr";
      hPullGP_Y_tr[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
      title.str("");
      title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_tr";
      hPullGP_Z_tr[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);

      title.str("");
      title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_rs";
      hPullGP_X_rs[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
      title.str("");
      title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_rs";
      hPullGP_Y_rs[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
      title.str("");
      title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_rs";
      hPullGP_Z_rs[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);

      if (((i == 2 || i == 4) && (j == 0 || j == 1)) || (i == 3 || i == 5)) {
        //mono
        title.str("");
        title << "Chi2Increment_mono_" << i + 1 << "-" << j + 1;
        hChi2Increment_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, 0, 100);

        title.str("");
        title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_ts_mono";
        hPullGP_X_ts_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
        title.str("");
        title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_ts_mono";
        hPullGP_Y_ts_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
        title.str("");
        title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_ts_mono";
        hPullGP_Z_ts_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);

        title.str("");
        title << "PullGM_X_" << i + 1 << "-" << j + 1 << "_ts_mono";
        hPullGM_X_ts_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
        title.str("");
        title << "PullGM_Y_" << i + 1 << "-" << j + 1 << "_ts_mono";
        hPullGM_Y_ts_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
        title.str("");
        title << "PullGM_Z_" << i + 1 << "-" << j + 1 << "_ts_mono";
        hPullGM_Z_ts_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);

        title.str("");
        title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_tr_mono";
        hPullGP_X_tr_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
        title.str("");
        title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_tr_mono";
        hPullGP_Y_tr_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
        title.str("");
        title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_tr_mono";
        hPullGP_Z_tr_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);

        title.str("");
        title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_rs_mono";
        hPullGP_X_rs_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
        title.str("");
        title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_rs_mono";
        hPullGP_Y_rs_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
        title.str("");
        title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_rs_mono";
        hPullGP_Z_rs_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);

        //stereo
        title.str("");
        title << "Chi2Increment_stereo_" << i + 1 << "-" << j + 1;
        hChi2Increment_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, 0, 100);

        title.str("");
        title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_ts_stereo";
        hPullGP_X_ts_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
        title.str("");
        title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_ts_stereo";
        hPullGP_Y_ts_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
        title.str("");
        title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_ts_stereo";
        hPullGP_Z_ts_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);

        title.str("");
        title << "PullGM_X_" << i + 1 << "-" << j + 1 << "_ts_stereo";
        hPullGM_X_ts_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
        title.str("");
        title << "PullGM_Y_" << i + 1 << "-" << j + 1 << "_ts_stereo";
        hPullGM_Y_ts_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
        title.str("");
        title << "PullGM_Z_" << i + 1 << "-" << j + 1 << "_ts_stereo";
        hPullGM_Z_ts_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);

        title.str("");
        title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_tr_stereo";
        hPullGP_X_tr_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
        title.str("");
        title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_tr_stereo";
        hPullGP_Y_tr_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
        title.str("");
        title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_tr_stereo";
        hPullGP_Z_tr_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);

        title.str("");
        title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_rs_stereo";
        hPullGP_X_rs_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
        title.str("");
        title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_rs_stereo";
        hPullGP_Y_rs_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
        title.str("");
        title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_rs_stereo";
        hPullGP_Z_rs_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
      }
    }
  hTotChi2Increment = new TH1F("TotChi2Increment", "TotChi2Increment", 1000, 0, 100);
  hTotChi2GoodHit = new TH1F("TotChi2GoodHit", "TotChi2GoodHit", 1000, 0, 100);
  hTotChi2BadHit = new TH1F("TotChi2BadHit", "TotChi2BadHit", 1000, 0, 100);
  hTotChi2DeltaHit = new TH1F("TotChi2DeltaHit", "TotChi2DeltaHit", 1000, 0, 100);
  hTotChi2NSharedHit = new TH1F("TotChi2NSharedHit", "TotChi2NSharedHit", 1000, 0, 100);
  hTotChi2SharedHit = new TH1F("TotChi2SharedHit", "TotChi2SharedHit", 1000, 0, 100);
  hProcess_vs_Chi2 = new TH2F("Process_vs_Chi2", "Process_vs_Chi2", 1000, 0, 100, 17, -0.5, 16.5);
  hClsize_vs_Chi2 = new TH2F("Clsize_vs_Chi2", "Clsize_vs_Chi2", 1000, 0, 100, 17, -0.5, 16.5);
  hPixClsize_vs_Chi2 = new TH2F("PixClsize_vs_Chi2", "PixClsize_vs_Chi2", 1000, 0, 100, 17, -0.5, 16.5);
  hPrjClsize_vs_Chi2 = new TH2F("PrjClsize_vs_Chi2", "PrjClsize_vs_Chi2", 1000, 0, 100, 17, -0.5, 16.5);
  hSt1Clsize_vs_Chi2 = new TH2F("St1Clsize_vs_Chi2", "St1Clsize_vs_Chi2", 1000, 0, 100, 17, -0.5, 16.5);
  hSt2Clsize_vs_Chi2 = new TH2F("St2Clsize_vs_Chi2", "St2Clsize_vs_Chi2", 1000, 0, 100, 17, -0.5, 16.5);
  hGoodHit_vs_Chi2 = new TH2F("GoodHit_vs_Chi2", "GoodHit_vs_Chi2", 10000, 0, 1000, 2, -0.5, 1.5);
  hClusterSize = new TH1F("ClusterSize", "ClusterSize", 40, -0.5, 39.5);
  hPixClusterSize = new TH1F("PixClusterSize", "PixClusterSize", 40, -0.5, 39.5);
  hPrjClusterSize = new TH1F("PrjClusterSize", "PrjClusterSize", 40, -0.5, 39.5);
  hSt1ClusterSize = new TH1F("St1ClusterSize", "St1ClusterSize", 40, -0.5, 39.5);
  hSt2ClusterSize = new TH1F("St2ClusterSize", "St2ClusterSize", 40, -0.5, 39.5);
  hSimHitVecSize = new TH1F("hSimHitVecSize", "hSimHitVecSize", 40, -0.5, 39.5);
  hPixSimHitVecSize = new TH1F("PixSimHitVecSize", "PixSimHitVecSize", 40, -0.5, 39.5);
  hPrjSimHitVecSize = new TH1F("PrjSimHitVecSize", "PrjSimHitVecSize", 40, -0.5, 39.5);
  hSt1SimHitVecSize = new TH1F("St1SimHitVecSize", "St1SimHitVecSize", 40, -0.5, 39.5);
  hSt2SimHitVecSize = new TH1F("St2SimHitVecSize", "St2SimHitVecSize", 40, -0.5, 39.5);
  goodbadmerged = new TH1F("goodbadmerged", "goodbadmerged", 5, 0.5, 5.5);
  energyLossRatio = new TH1F("energyLossRatio", "energyLossRatio", 100, 0, 1);
  mergedPull = new TH1F("mergedPull", "mergedPull", 200, 0, 20);
  probXgood = new TH1F("probXgood", "probXgood", 110, 0, 1.1);
  probXbad = new TH1F("probXbad", "probXbad", 110, 0, 1.1);
  probXdelta = new TH1F("probXdelta", "probXdelta", 110, 0, 1.1);
  probXshared = new TH1F("probXshared", "probXshared", 110, 0, 1.1);
  probXnoshare = new TH1F("probXnoshare", "probXnoshare", 110, 0, 1.1);
  probYgood = new TH1F("probYgood", "probYgood", 110, 0, 1.1);
  probYbad = new TH1F("probYbad", "probYbad", 110, 0, 1.1);
  probYdelta = new TH1F("probYdelta", "probYdelta", 110, 0, 1.1);
  probYshared = new TH1F("probYshared", "probYshared", 110, 0, 1.1);
  probYnoshare = new TH1F("probYnoshare", "probYnoshare", 110, 0, 1.1);
}

computeChi2Increment() [1/2]

template<unsigned int D>

double TestTrackHits::computeChi2Increment ( MeasurementExtractor  me, TransientTrackingRecHit::ConstRecHitPointer  rhit  )

private

Definition at line 1311 of file TestTrackHits.cc.

References hlt_dqm_clientPB-live_cfg::me, dttmaxenums::R, and alignCSCRings::r.

Referenced by analyze().

                                                                                                                  {
  typedef typename AlgebraicROOTObject<D>::Vector VecD;
  typedef typename AlgebraicROOTObject<D, D>::SymMatrix SMatDD;
  VecD r = asSVector<D>(rhit->parameters()) - me.measuredParameters<D>(*rhit);

  SMatDD R = asSMatrix<D>(rhit->parametersError()) + me.measuredError<D>(*rhit);
  R.Invert();
  return ROOT::Math::Similarity(r, R);
}

computeChi2Increment() [2/2]

double TestTrackHits::computeChi2Increment ( MeasurementExtractor  me, TransientTrackingRecHit::ConstRecHitPointer  hit  )

inlineprivate

Definition at line 69 of file TestTrackHits.h.

References Exception, and hlt_dqm_clientPB-live_cfg::me.

                                                                                                      {
    switch (hit->dimension()) {
      case 1:
        return computeChi2Increment<1>(me, hit);
      case 2:
        return computeChi2Increment<2>(me, hit);
      case 3:
        return computeChi2Increment<3>(me, hit);
      case 4:
        return computeChi2Increment<4>(me, hit);
      case 5:
        return computeChi2Increment<5>(me, hit);
    }
    throw cms::Exception("CkfDebugger error: rechit of dimension not 1,2,3,4,5");
  }

endJob()

void TestTrackHits::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 978 of file TestTrackHits.cc.

References energyLossRatio, file, goodbadmerged, hChi2BadHit, hChi2DeltaHit, hChi2GoodHit, hChi2Increment, hChi2Increment_mono, hChi2Increment_stereo, hChi2IncrementVsEta, hChi2NSharedHit, hChi2SharedHit, hClsize_vs_Chi2, hClusterSize, hGoodHit_vs_Chi2, hPixClsize_vs_Chi2, hPixClusterSize, hPixSimHitVecSize, hPrjClsize_vs_Chi2, hPrjClusterSize, hPrjSimHitVecSize, hProcess_vs_Chi2, hPullGM_X_ts, hPullGM_X_ts_mono, hPullGM_X_ts_stereo, hPullGM_Y_ts, hPullGM_Y_ts_mono, hPullGM_Y_ts_stereo, hPullGM_Z_ts, hPullGM_Z_ts_mono, hPullGM_Z_ts_stereo, hPullGP_X_rs, hPullGP_X_rs_mono, hPullGP_X_rs_stereo, hPullGP_X_tr, hPullGP_X_tr_mono, hPullGP_X_tr_stereo, hPullGP_X_ts, hPullGP_X_ts_mono, hPullGP_X_ts_stereo, hPullGP_Y_rs, hPullGP_Y_rs_mono, hPullGP_Y_rs_stereo, hPullGP_Y_tr, hPullGP_Y_tr_mono, hPullGP_Y_tr_stereo, hPullGP_Y_ts, hPullGP_Y_ts_mono, hPullGP_Y_ts_stereo, hPullGP_Z_rs, hPullGP_Z_rs_mono, hPullGP_Z_rs_stereo, hPullGP_Z_tr, hPullGP_Z_tr_mono, hPullGP_Z_tr_stereo, hPullGP_Z_ts, hPullGP_Z_ts_mono, hPullGP_Z_ts_stereo, hSimHitVecSize, hSt1Clsize_vs_Chi2, hSt1ClusterSize, hSt1SimHitVecSize, hSt2Clsize_vs_Chi2, hSt2ClusterSize, hSt2SimHitVecSize, hTotChi2BadHit, hTotChi2DeltaHit, hTotChi2GoodHit, hTotChi2Increment, hTotChi2NSharedHit, hTotChi2SharedHit, mps_fire::i, dqmiolumiharvest::j, mergedPull, probXbad, probXdelta, probXgood, probXnoshare, probXshared, probYbad, probYdelta, probYgood, probYnoshare, probYshared, and title.

                           {
  //file->Write();
  TDirectory* chi2d = file->mkdir("Chi2Increment");

  TDirectory* gp_ts = file->mkdir("GP_TSOS-SimHit");
  TDirectory* gm_ts = file->mkdir("GM_TSOS-SimHit");
  TDirectory* gp_tr = file->mkdir("GP_TSOS-RecHit");
  TDirectory* gp_rs = file->mkdir("GP_RecHit-SimHit");

  TDirectory* gp_tsx = gp_ts->mkdir("X");
  TDirectory* gp_tsy = gp_ts->mkdir("Y");
  TDirectory* gp_tsz = gp_ts->mkdir("Z");
  TDirectory* gm_tsx = gm_ts->mkdir("X");
  TDirectory* gm_tsy = gm_ts->mkdir("Y");
  TDirectory* gm_tsz = gm_ts->mkdir("Z");
  TDirectory* gp_trx = gp_tr->mkdir("X");
  TDirectory* gp_try = gp_tr->mkdir("Y");
  TDirectory* gp_trz = gp_tr->mkdir("Z");
  TDirectory* gp_rsx = gp_rs->mkdir("X");
  TDirectory* gp_rsy = gp_rs->mkdir("Y");
  TDirectory* gp_rsz = gp_rs->mkdir("Z");

  TDirectory* gp_tsx_mono = gp_ts->mkdir("MONOX");
  TDirectory* gp_tsy_mono = gp_ts->mkdir("MONOY");
  TDirectory* gp_tsz_mono = gp_ts->mkdir("MONOZ");
  TDirectory* gm_tsx_mono = gm_ts->mkdir("MONOX");
  TDirectory* gm_tsy_mono = gm_ts->mkdir("MONOY");
  TDirectory* gm_tsz_mono = gm_ts->mkdir("MONOZ");
  TDirectory* gp_trx_mono = gp_tr->mkdir("MONOX");
  TDirectory* gp_try_mono = gp_tr->mkdir("MONOY");
  TDirectory* gp_trz_mono = gp_tr->mkdir("MONOZ");
  TDirectory* gp_rsx_mono = gp_rs->mkdir("MONOX");
  TDirectory* gp_rsy_mono = gp_rs->mkdir("MONOY");
  TDirectory* gp_rsz_mono = gp_rs->mkdir("MONOZ");

  TDirectory* gp_tsx_stereo = gp_ts->mkdir("STEREOX");
  TDirectory* gp_tsy_stereo = gp_ts->mkdir("STEREOY");
  TDirectory* gp_tsz_stereo = gp_ts->mkdir("STEREOZ");
  TDirectory* gm_tsx_stereo = gm_ts->mkdir("STEREOX");
  TDirectory* gm_tsy_stereo = gm_ts->mkdir("STEREOY");
  TDirectory* gm_tsz_stereo = gm_ts->mkdir("STEREOZ");
  TDirectory* gp_trx_stereo = gp_tr->mkdir("STEREOX");
  TDirectory* gp_try_stereo = gp_tr->mkdir("STEREOY");
  TDirectory* gp_trz_stereo = gp_tr->mkdir("STEREOZ");
  TDirectory* gp_rsx_stereo = gp_rs->mkdir("STEREOX");
  TDirectory* gp_rsy_stereo = gp_rs->mkdir("STEREOY");
  TDirectory* gp_rsz_stereo = gp_rs->mkdir("STEREOZ");

  chi2d->cd();
  hTotChi2Increment->Write();
  hTotChi2GoodHit->Write();
  hTotChi2BadHit->Write();
  hTotChi2DeltaHit->Write();
  hTotChi2NSharedHit->Write();
  hTotChi2SharedHit->Write();
  hProcess_vs_Chi2->Write();
  hClsize_vs_Chi2->Write();
  hPixClsize_vs_Chi2->Write();
  hPrjClsize_vs_Chi2->Write();
  hSt1Clsize_vs_Chi2->Write();
  hSt2Clsize_vs_Chi2->Write();
  hGoodHit_vs_Chi2->Write();
  hClusterSize->Write();
  hPixClusterSize->Write();
  hPrjClusterSize->Write();
  hSt1ClusterSize->Write();
  hSt2ClusterSize->Write();
  hSimHitVecSize->Write();
  hPixSimHitVecSize->Write();
  hPrjSimHitVecSize->Write();
  hSt1SimHitVecSize->Write();
  hSt2SimHitVecSize->Write();
  goodbadmerged->Write();
  energyLossRatio->Write();
  mergedPull->Write();
  probXgood->Write();
  probXbad->Write();
  probXdelta->Write();
  probXshared->Write();
  probXnoshare->Write();
  probYgood->Write();
  probYbad->Write();
  probYdelta->Write();
  probYshared->Write();
  probYnoshare->Write();
  for (int i = 0; i != 6; i++)
    for (int j = 0; j != 9; j++) {
      if (i == 0 && j > 2)
        break;
      if (i == 1 && j > 1)
        break;
      if (i == 2 && j > 3)
        break;
      if (i == 3 && j > 2)
        break;
      if (i == 4 && j > 5)
        break;
      if (i == 5 && j > 8)
        break;
      chi2d->cd();
      title.str("");
      title << "Chi2Increment_" << i + 1 << "-" << j + 1;
      hChi2Increment[title.str()]->Write();
      title.str("");
      title << "Chi2IncrementVsEta_" << i + 1 << "-" << j + 1;
      hChi2IncrementVsEta[title.str()]->Write();
      title.str("");
      title << "Chi2GoodHit_" << i + 1 << "-" << j + 1;
      hChi2GoodHit[title.str()]->Write();
      title.str("");
      title << "Chi2BadHit_" << i + 1 << "-" << j + 1;
      hChi2BadHit[title.str()]->Write();
      title.str("");
      title << "Chi2DeltaHit_" << i + 1 << "-" << j + 1;
      hChi2DeltaHit[title.str()]->Write();
      title.str("");
      title << "Chi2NSharedHit_" << i + 1 << "-" << j + 1;
      hChi2NSharedHit[title.str()]->Write();
      title.str("");
      title << "Chi2SharedHit_" << i + 1 << "-" << j + 1;
      hChi2SharedHit[title.str()]->Write();

      gp_ts->cd();
      gp_tsx->cd();
      title.str("");
      title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_ts";
      hPullGP_X_ts[title.str()]->Write();
      gp_tsy->cd();
      title.str("");
      title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_ts";
      hPullGP_Y_ts[title.str()]->Write();
      gp_tsz->cd();
      title.str("");
      title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_ts";
      hPullGP_Z_ts[title.str()]->Write();

      gm_ts->cd();
      gm_tsx->cd();
      title.str("");
      title << "PullGM_X_" << i + 1 << "-" << j + 1 << "_ts";
      hPullGM_X_ts[title.str()]->Write();
      gm_tsy->cd();
      title.str("");
      title << "PullGM_Y_" << i + 1 << "-" << j + 1 << "_ts";
      hPullGM_Y_ts[title.str()]->Write();
      gm_tsz->cd();
      title.str("");
      title << "PullGM_Z_" << i + 1 << "-" << j + 1 << "_ts";
      hPullGM_Z_ts[title.str()]->Write();

      gp_tr->cd();
      gp_trx->cd();
      title.str("");
      title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_tr";
      hPullGP_X_tr[title.str()]->Write();
      gp_try->cd();
      title.str("");
      title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_tr";
      hPullGP_Y_tr[title.str()]->Write();
      gp_trz->cd();
      title.str("");
      title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_tr";
      hPullGP_Z_tr[title.str()]->Write();

      gp_rs->cd();
      gp_rsx->cd();
      title.str("");
      title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_rs";
      hPullGP_X_rs[title.str()]->Write();
      gp_rsy->cd();
      title.str("");
      title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_rs";
      hPullGP_Y_rs[title.str()]->Write();
      gp_rsz->cd();
      title.str("");
      title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_rs";
      hPullGP_Z_rs[title.str()]->Write();

      if (((i == 2 || i == 4) && (j == 0 || j == 1)) || (i == 3 || i == 5)) {
        chi2d->cd();
        title.str("");
        title << "Chi2Increment_mono_" << i + 1 << "-" << j + 1;
        hChi2Increment_mono[title.str()]->Write();
        title.str("");
        title << "Chi2Increment_stereo_" << i + 1 << "-" << j + 1;
        hChi2Increment_stereo[title.str()]->Write();
        //mono
        gp_ts->cd();
        gp_tsx_mono->cd();
        title.str("");
        title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_ts_mono";
        hPullGP_X_ts_mono[title.str()]->Write();
        gp_tsy_mono->cd();
        title.str("");
        title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_ts_mono";
        hPullGP_Y_ts_mono[title.str()]->Write();
        gp_tsz_mono->cd();
        title.str("");
        title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_ts_mono";
        hPullGP_Z_ts_mono[title.str()]->Write();

        gm_ts->cd();
        gm_tsx_mono->cd();
        title.str("");
        title << "PullGM_X_" << i + 1 << "-" << j + 1 << "_ts_mono";
        hPullGM_X_ts_mono[title.str()]->Write();
        gm_tsy_mono->cd();
        title.str("");
        title << "PullGM_Y_" << i + 1 << "-" << j + 1 << "_ts_mono";
        hPullGM_Y_ts_mono[title.str()]->Write();
        gm_tsz_mono->cd();
        title.str("");
        title << "PullGM_Z_" << i + 1 << "-" << j + 1 << "_ts_mono";
        hPullGM_Z_ts_mono[title.str()]->Write();

        gp_tr->cd();
        gp_trx_mono->cd();
        title.str("");
        title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_tr_mono";
        hPullGP_X_tr_mono[title.str()]->Write();
        gp_try_mono->cd();
        title.str("");
        title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_tr_mono";
        hPullGP_Y_tr_mono[title.str()]->Write();
        gp_trz_mono->cd();
        title.str("");
        title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_tr_mono";
        hPullGP_Z_tr_mono[title.str()]->Write();

        gp_rs->cd();
        gp_rsx_mono->cd();
        title.str("");
        title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_rs_mono";
        hPullGP_X_rs_mono[title.str()]->Write();
        gp_rsy_mono->cd();
        title.str("");
        title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_rs_mono";
        hPullGP_Y_rs_mono[title.str()]->Write();
        gp_rsz_mono->cd();
        title.str("");
        title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_rs_mono";
        hPullGP_Z_rs_mono[title.str()]->Write();

        //stereo
        gp_ts->cd();
        gp_tsx_stereo->cd();
        title.str("");
        title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_ts_stereo";
        hPullGP_X_ts_stereo[title.str()]->Write();
        gp_tsy_stereo->cd();
        title.str("");
        title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_ts_stereo";
        hPullGP_Y_ts_stereo[title.str()]->Write();
        gp_tsz_stereo->cd();
        title.str("");
        title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_ts_stereo";
        hPullGP_Z_ts_stereo[title.str()]->Write();

        gm_ts->cd();
        gm_tsx_stereo->cd();
        title.str("");
        title << "PullGM_X_" << i + 1 << "-" << j + 1 << "_ts_stereo";
        hPullGM_X_ts_stereo[title.str()]->Write();
        gm_tsy_stereo->cd();
        title.str("");
        title << "PullGM_Y_" << i + 1 << "-" << j + 1 << "_ts_stereo";
        hPullGM_Y_ts_stereo[title.str()]->Write();
        gm_tsz_stereo->cd();
        title.str("");
        title << "PullGM_Z_" << i + 1 << "-" << j + 1 << "_ts_stereo";
        hPullGM_Z_ts_stereo[title.str()]->Write();

        gp_tr->cd();
        gp_trx_stereo->cd();
        title.str("");
        title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_tr_stereo";
        hPullGP_X_tr_stereo[title.str()]->Write();
        gp_try_stereo->cd();
        title.str("");
        title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_tr_stereo";
        hPullGP_Y_tr_stereo[title.str()]->Write();
        gp_trz_stereo->cd();
        title.str("");
        title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_tr_stereo";
        hPullGP_Z_tr_stereo[title.str()]->Write();

        gp_rs->cd();
        gp_rsx_stereo->cd();
        title.str("");
        title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_rs_stereo";
        hPullGP_X_rs_stereo[title.str()]->Write();
        gp_rsy_stereo->cd();
        title.str("");
        title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_rs_stereo";
        hPullGP_Y_rs_stereo[title.str()]->Write();
        gp_rsz_stereo->cd();
        title.str("");
        title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_rs_stereo";
        hPullGP_Z_rs_stereo[title.str()]->Write();
      }
    }

  file->Close();
}

endRun()

void TestTrackHits::endRun ( edm::Run const &  run, const edm::EventSetup &    )

inlineoverrideprivate

Definition at line 62 of file TestTrackHits.h.

{}

projectHit()

std::pair< LocalPoint, LocalVector > TestTrackHits::projectHit ( const PSimHit &  hit, const StripGeomDetUnit *  stripDet, const BoundPlane &  plane  )

private

Definition at line 1285 of file TestTrackHits.cc.

References funct::cos(), DeadROC_duringRun::dir, l1tEGammaCrystalsEmulatorProducer_cfi::scale, funct::sin(), StripGeomDetUnit::specificTopology(), StripTopology::strip(), StripTopology::stripAngle(), GeomDet::surface(), Surface::toGlobal(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by analyze().

                                                                                      {
  const StripTopology& topol = stripDet->specificTopology();
  GlobalPoint globalpos = stripDet->surface().toGlobal(hit.localPosition());
  LocalPoint localHit = plane.toLocal(globalpos);
  //track direction
  LocalVector locdir = hit.localDirection();
  //rotate track in new frame

  GlobalVector globaldir = stripDet->surface().toGlobal(locdir);
  LocalVector dir = plane.toLocal(globaldir);
  float scale = -localHit.z() / dir.z();

  LocalPoint projectedPos = localHit + scale * dir;

  float selfAngle = topol.stripAngle(topol.strip(hit.localPosition()));

  LocalVector stripDir(sin(selfAngle), cos(selfAngle), 0);  // vector along strip in hit frame

  LocalVector localStripDir(plane.toLocal(stripDet->surface().toGlobal(stripDir)));

  return std::pair<LocalPoint, LocalVector>(projectedPos, localStripDir);
}

Member Data Documentation

builderName

std::string TestTrackHits::builderName

private

Definition at line 88 of file TestTrackHits.h.

Referenced by TestTrackHits().

energyLossRatio

TH1F * TestTrackHits::energyLossRatio

private

Definition at line 140 of file TestTrackHits.h.

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

file

TFile* TestTrackHits::file

private

Definition at line 111 of file TestTrackHits.h.

Referenced by beginRun(), and endJob().

goodbadmerged

TH1F* TestTrackHits::goodbadmerged

private

Definition at line 140 of file TestTrackHits.h.

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

hChi2BadHit

std::map<std::string, TH1F *> TestTrackHits::hChi2BadHit

private

Definition at line 130 of file TestTrackHits.h.

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

hChi2DeltaHit

std::map<std::string, TH1F *> TestTrackHits::hChi2DeltaHit

private

Definition at line 131 of file TestTrackHits.h.

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

hChi2GoodHit

std::map<std::string, TH1F *> TestTrackHits::hChi2GoodHit

private

Definition at line 129 of file TestTrackHits.h.

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

hChi2Increment

std::map<std::string, TH1F *> TestTrackHits::hChi2Increment

private

Definition at line 125 of file TestTrackHits.h.

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

hChi2Increment_mono

std::map<std::string, TH1F *> TestTrackHits::hChi2Increment_mono

private

Definition at line 126 of file TestTrackHits.h.

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

hChi2Increment_stereo

std::map<std::string, TH1F *> TestTrackHits::hChi2Increment_stereo

private

Definition at line 127 of file TestTrackHits.h.

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

hChi2IncrementVsEta

std::map<std::string, TH2F *> TestTrackHits::hChi2IncrementVsEta

private

Definition at line 128 of file TestTrackHits.h.

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

hChi2NSharedHit

std::map<std::string, TH1F *> TestTrackHits::hChi2NSharedHit

private

Definition at line 132 of file TestTrackHits.h.

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

hChi2SharedHit

std::map<std::string, TH1F *> TestTrackHits::hChi2SharedHit

private

Definition at line 133 of file TestTrackHits.h.

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

hClsize_vs_Chi2

TH2F * TestTrackHits::hClsize_vs_Chi2

private

Definition at line 136 of file TestTrackHits.h.

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

hClusterSize

TH1F* TestTrackHits::hClusterSize

private

Definition at line 138 of file TestTrackHits.h.

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

hGoodHit_vs_Chi2

TH2F * TestTrackHits::hGoodHit_vs_Chi2

private

Definition at line 136 of file TestTrackHits.h.

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

hPixClsize_vs_Chi2

TH2F* TestTrackHits::hPixClsize_vs_Chi2

private

Definition at line 137 of file TestTrackHits.h.

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

hPixClusterSize

TH1F * TestTrackHits::hPixClusterSize

private

Definition at line 138 of file TestTrackHits.h.

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

hPixSimHitVecSize

TH1F * TestTrackHits::hPixSimHitVecSize

private

Definition at line 139 of file TestTrackHits.h.

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

hPrjClsize_vs_Chi2

TH2F * TestTrackHits::hPrjClsize_vs_Chi2

private

Definition at line 137 of file TestTrackHits.h.

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

hPrjClusterSize

TH1F * TestTrackHits::hPrjClusterSize

private

Definition at line 138 of file TestTrackHits.h.

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

hPrjSimHitVecSize

TH1F * TestTrackHits::hPrjSimHitVecSize

private

Definition at line 139 of file TestTrackHits.h.

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

hProcess_vs_Chi2

TH2F* TestTrackHits::hProcess_vs_Chi2

private

Definition at line 136 of file TestTrackHits.h.

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

hPullGM_X_ts

std::map<std::string, TH1F *> TestTrackHits::hPullGM_X_ts

private

Definition at line 116 of file TestTrackHits.h.

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

hPullGM_X_ts_mono

std::map<std::string, TH1F *> TestTrackHits::hPullGM_X_ts_mono

private

Definition at line 147 of file TestTrackHits.h.

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

hPullGM_X_ts_stereo

std::map<std::string, TH1F *> TestTrackHits::hPullGM_X_ts_stereo

private

Definition at line 160 of file TestTrackHits.h.

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

hPullGM_Y_ts

std::map<std::string, TH1F *> TestTrackHits::hPullGM_Y_ts

private

Definition at line 117 of file TestTrackHits.h.

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

hPullGM_Y_ts_mono

std::map<std::string, TH1F *> TestTrackHits::hPullGM_Y_ts_mono

private

Definition at line 148 of file TestTrackHits.h.

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

hPullGM_Y_ts_stereo

std::map<std::string, TH1F *> TestTrackHits::hPullGM_Y_ts_stereo

private

Definition at line 161 of file TestTrackHits.h.

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

hPullGM_Z_ts

std::map<std::string, TH1F *> TestTrackHits::hPullGM_Z_ts

private

Definition at line 118 of file TestTrackHits.h.

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

hPullGM_Z_ts_mono

std::map<std::string, TH1F *> TestTrackHits::hPullGM_Z_ts_mono

private

Definition at line 149 of file TestTrackHits.h.

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

hPullGM_Z_ts_stereo

std::map<std::string, TH1F *> TestTrackHits::hPullGM_Z_ts_stereo

private

Definition at line 162 of file TestTrackHits.h.

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

hPullGP_X_rs

std::map<std::string, TH1F *> TestTrackHits::hPullGP_X_rs

private

Definition at line 119 of file TestTrackHits.h.

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

hPullGP_X_rs_mono

std::map<std::string, TH1F *> TestTrackHits::hPullGP_X_rs_mono

private

Definition at line 150 of file TestTrackHits.h.

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

hPullGP_X_rs_stereo

std::map<std::string, TH1F *> TestTrackHits::hPullGP_X_rs_stereo

private

Definition at line 163 of file TestTrackHits.h.

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

hPullGP_X_tr

std::map<std::string, TH1F *> TestTrackHits::hPullGP_X_tr

private

Definition at line 122 of file TestTrackHits.h.

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

hPullGP_X_tr_mono

std::map<std::string, TH1F *> TestTrackHits::hPullGP_X_tr_mono

private

Definition at line 153 of file TestTrackHits.h.

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

hPullGP_X_tr_stereo

std::map<std::string, TH1F *> TestTrackHits::hPullGP_X_tr_stereo

private

Definition at line 166 of file TestTrackHits.h.

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

hPullGP_X_ts

std::map<std::string, TH1F *> TestTrackHits::hPullGP_X_ts

private

Definition at line 113 of file TestTrackHits.h.

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

hPullGP_X_ts_mono

std::map<std::string, TH1F *> TestTrackHits::hPullGP_X_ts_mono

private

Definition at line 144 of file TestTrackHits.h.

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

hPullGP_X_ts_stereo

std::map<std::string, TH1F *> TestTrackHits::hPullGP_X_ts_stereo

private

Definition at line 157 of file TestTrackHits.h.

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

hPullGP_Y_rs

std::map<std::string, TH1F *> TestTrackHits::hPullGP_Y_rs

private

Definition at line 120 of file TestTrackHits.h.

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

hPullGP_Y_rs_mono

std::map<std::string, TH1F *> TestTrackHits::hPullGP_Y_rs_mono

private

Definition at line 151 of file TestTrackHits.h.

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

hPullGP_Y_rs_stereo

std::map<std::string, TH1F *> TestTrackHits::hPullGP_Y_rs_stereo

private

Definition at line 164 of file TestTrackHits.h.

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

hPullGP_Y_tr

std::map<std::string, TH1F *> TestTrackHits::hPullGP_Y_tr

private

Definition at line 123 of file TestTrackHits.h.

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

hPullGP_Y_tr_mono

std::map<std::string, TH1F *> TestTrackHits::hPullGP_Y_tr_mono

private

Definition at line 154 of file TestTrackHits.h.

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

hPullGP_Y_tr_stereo

std::map<std::string, TH1F *> TestTrackHits::hPullGP_Y_tr_stereo

private

Definition at line 167 of file TestTrackHits.h.

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

hPullGP_Y_ts

std::map<std::string, TH1F *> TestTrackHits::hPullGP_Y_ts

private

Definition at line 114 of file TestTrackHits.h.

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

hPullGP_Y_ts_mono

std::map<std::string, TH1F *> TestTrackHits::hPullGP_Y_ts_mono

private

Definition at line 145 of file TestTrackHits.h.

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

hPullGP_Y_ts_stereo

std::map<std::string, TH1F *> TestTrackHits::hPullGP_Y_ts_stereo

private

Definition at line 158 of file TestTrackHits.h.

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

hPullGP_Z_rs

std::map<std::string, TH1F *> TestTrackHits::hPullGP_Z_rs

private

Definition at line 121 of file TestTrackHits.h.

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

hPullGP_Z_rs_mono

std::map<std::string, TH1F *> TestTrackHits::hPullGP_Z_rs_mono

private

Definition at line 152 of file TestTrackHits.h.

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

hPullGP_Z_rs_stereo

std::map<std::string, TH1F *> TestTrackHits::hPullGP_Z_rs_stereo

private

Definition at line 165 of file TestTrackHits.h.

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

hPullGP_Z_tr

std::map<std::string, TH1F *> TestTrackHits::hPullGP_Z_tr

private

Definition at line 124 of file TestTrackHits.h.

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

hPullGP_Z_tr_mono

std::map<std::string, TH1F *> TestTrackHits::hPullGP_Z_tr_mono

private

Definition at line 155 of file TestTrackHits.h.

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

hPullGP_Z_tr_stereo

std::map<std::string, TH1F *> TestTrackHits::hPullGP_Z_tr_stereo

private

Definition at line 168 of file TestTrackHits.h.

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

hPullGP_Z_ts

std::map<std::string, TH1F *> TestTrackHits::hPullGP_Z_ts

private

Definition at line 115 of file TestTrackHits.h.

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

hPullGP_Z_ts_mono

std::map<std::string, TH1F *> TestTrackHits::hPullGP_Z_ts_mono

private

Definition at line 146 of file TestTrackHits.h.

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

hPullGP_Z_ts_stereo

std::map<std::string, TH1F *> TestTrackHits::hPullGP_Z_ts_stereo

private

Definition at line 159 of file TestTrackHits.h.

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

hSimHitVecSize

TH1F* TestTrackHits::hSimHitVecSize

private

Definition at line 139 of file TestTrackHits.h.

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

hSt1Clsize_vs_Chi2

TH2F * TestTrackHits::hSt1Clsize_vs_Chi2

private

Definition at line 137 of file TestTrackHits.h.

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

hSt1ClusterSize

TH1F * TestTrackHits::hSt1ClusterSize

private

Definition at line 138 of file TestTrackHits.h.

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

hSt1SimHitVecSize

TH1F * TestTrackHits::hSt1SimHitVecSize

private

Definition at line 139 of file TestTrackHits.h.

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

hSt2Clsize_vs_Chi2

TH2F * TestTrackHits::hSt2Clsize_vs_Chi2

private

Definition at line 137 of file TestTrackHits.h.

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

hSt2ClusterSize

TH1F * TestTrackHits::hSt2ClusterSize

private

Definition at line 138 of file TestTrackHits.h.

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

hSt2SimHitVecSize

TH1F * TestTrackHits::hSt2SimHitVecSize

private

Definition at line 139 of file TestTrackHits.h.

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

hTotChi2BadHit

TH1F * TestTrackHits::hTotChi2BadHit

private

Definition at line 134 of file TestTrackHits.h.

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

hTotChi2DeltaHit

TH1F * TestTrackHits::hTotChi2DeltaHit

private

Definition at line 134 of file TestTrackHits.h.

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

hTotChi2GoodHit

TH1F * TestTrackHits::hTotChi2GoodHit

private

Definition at line 134 of file TestTrackHits.h.

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

hTotChi2Increment

TH1F* TestTrackHits::hTotChi2Increment

private

Definition at line 134 of file TestTrackHits.h.

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

hTotChi2NSharedHit

TH1F * TestTrackHits::hTotChi2NSharedHit

private

Definition at line 134 of file TestTrackHits.h.

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

hTotChi2SharedHit

TH1F * TestTrackHits::hTotChi2SharedHit

private

Definition at line 134 of file TestTrackHits.h.

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

mergedPull

TH1F * TestTrackHits::mergedPull

private

Definition at line 140 of file TestTrackHits.h.

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

out

std::string TestTrackHits::out

private

Definition at line 92 of file TestTrackHits.h.

Referenced by beginRun(), and TestTrackHits().

probXbad

TH1F * TestTrackHits::probXbad

private

Definition at line 141 of file TestTrackHits.h.

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

probXdelta

TH1F * TestTrackHits::probXdelta

private

Definition at line 141 of file TestTrackHits.h.

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

probXgood

TH1F* TestTrackHits::probXgood

private

Definition at line 141 of file TestTrackHits.h.

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

probXnoshare

TH1F * TestTrackHits::probXnoshare

private

Definition at line 141 of file TestTrackHits.h.

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

probXshared

TH1F * TestTrackHits::probXshared

private

Definition at line 141 of file TestTrackHits.h.

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

probYbad

TH1F * TestTrackHits::probYbad

private

Definition at line 142 of file TestTrackHits.h.

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

probYdelta

TH1F * TestTrackHits::probYdelta

private

Definition at line 142 of file TestTrackHits.h.

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

probYgood

TH1F* TestTrackHits::probYgood

private

Definition at line 142 of file TestTrackHits.h.

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

probYnoshare

TH1F * TestTrackHits::probYnoshare

private

Definition at line 142 of file TestTrackHits.h.

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

probYshared

TH1F * TestTrackHits::probYshared

private

Definition at line 142 of file TestTrackHits.h.

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

propagatorName

std::string TestTrackHits::propagatorName

private

Definition at line 87 of file TestTrackHits.h.

Referenced by TestTrackHits().

srcName

std::string TestTrackHits::srcName

private

Definition at line 89 of file TestTrackHits.h.

Referenced by analyze(), and TestTrackHits().

theBuilder

edm::ESHandle<TransientTrackingRecHitBuilder> TestTrackHits::theBuilder

private

Definition at line 96 of file TestTrackHits.h.

Referenced by analyze(), and beginRun().

theBuilderToken

edm::ESGetToken<TransientTrackingRecHitBuilder, TransientRecHitRecord> TestTrackHits::theBuilderToken

private

Definition at line 101 of file TestTrackHits.h.

Referenced by beginRun(), and TestTrackHits().

theG

edm::ESHandle<TrackerGeometry> TestTrackHits::theG

private

Definition at line 93 of file TestTrackHits.h.

Referenced by beginRun().

theGToken

edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> TestTrackHits::theGToken

private

Definition at line 98 of file TestTrackHits.h.

Referenced by beginRun(), and TestTrackHits().

theMF

edm::ESHandle<MagneticField> TestTrackHits::theMF

private

Definition at line 94 of file TestTrackHits.h.

Referenced by analyze(), and beginRun().

theMFToken

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> TestTrackHits::theMFToken

private

Definition at line 99 of file TestTrackHits.h.

Referenced by beginRun(), and TestTrackHits().

thePropagator

edm::ESHandle<Propagator> TestTrackHits::thePropagator

private

Definition at line 95 of file TestTrackHits.h.

Referenced by beginRun().

thePropagatorToken

edm::ESGetToken<Propagator, TrackingComponentsRecord> TestTrackHits::thePropagatorToken

private

Definition at line 100 of file TestTrackHits.h.

Referenced by beginRun(), and TestTrackHits().

theTopoToken

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> TestTrackHits::theTopoToken

private

Definition at line 103 of file TestTrackHits.h.

Referenced by analyze(), and TestTrackHits().

theUpdator

edm::ESHandle<TrajectoryStateUpdator> TestTrackHits::theUpdator

private

Definition at line 97 of file TestTrackHits.h.

Referenced by beginRun().

theUpdatorToken

edm::ESGetToken<TrajectoryStateUpdator, TrackingComponentsRecord> TestTrackHits::theUpdatorToken

private

Definition at line 102 of file TestTrackHits.h.

Referenced by beginRun(), and TestTrackHits().

title

std::stringstream TestTrackHits::title

private

Definition at line 112 of file TestTrackHits.h.

Referenced by analyze(), beginRun(), endJob(), and Formatter.SimpleHTMLFormatter::headers().

tpName

std::string TestTrackHits::tpName

private

Definition at line 90 of file TestTrackHits.h.

Referenced by analyze(), and TestTrackHits().

trackAssociator

edm::Handle<reco::TrackToTrackingParticleAssociator> TestTrackHits::trackAssociator

private

Definition at line 104 of file TestTrackHits.h.

Referenced by analyze().

trackCollectionHandle

edm::Handle<edm::View<reco::Track> > TestTrackHits::trackCollectionHandle

private

Definition at line 106 of file TestTrackHits.h.

Referenced by analyze().

trackerHitAssociatorConfig_

TrackerHitAssociator::Config TestTrackHits::trackerHitAssociatorConfig_

private

Definition at line 85 of file TestTrackHits.h.

Referenced by analyze().

trackingParticleCollectionHandle

edm::Handle<TrackingParticleCollection> TestTrackHits::trackingParticleCollectionHandle

private

Definition at line 108 of file TestTrackHits.h.

Referenced by analyze().

trajCollectionHandle

edm::Handle<std::vector<Trajectory> > TestTrackHits::trajCollectionHandle

private

Definition at line 105 of file TestTrackHits.h.

Referenced by analyze().

trajTrackAssociationCollectionHandle

edm::Handle<TrajTrackAssociationCollection> TestTrackHits::trajTrackAssociationCollectionHandle

private

Definition at line 107 of file TestTrackHits.h.

Referenced by analyze().

updatorName

std::string TestTrackHits::updatorName

private

Definition at line 91 of file TestTrackHits.h.

Referenced by TestTrackHits().