#include <RecoTracker/DebugTools/src/TestSmoothHits.cc>

Inheritance diagram for TestSmoothHits:

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

	~TestSmoothHits () 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

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

TFile *	file

edm::ESHandle< TrajectoryFitter >	fit

edm::ESGetToken< TrajectoryFitter, TrajectoryFitter::Record >	fitToken

std::string	fname

TH2F *	hChi2_vs_clsize

TH2F *	hChi2_vs_Process

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

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 *	hTotChi2Increment

double	maxeta

double	mineta

std::string	propagatorName

edm::ESHandle< TrajectorySmoother >	smooth

edm::ESGetToken< TrajectorySmoother, TrajectoryFitter::Record >	smoothToken

std::string	sname

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::Handle< TrackCandidateCollection >	theTCCollection

edm::EDGetTokenT< TrackCandidateCollection >	theTCCollectionToken

std::stringstream	title

TrackerHitAssociator::Config	trackerHitAssociatorConfig_

edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	tTopoToken

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 50 of file TestSmoothHits.h.

Constructor & Destructor Documentation

◆ TestSmoothHits()

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

explicit

Definition at line 23 of file TestSmoothHits.cc.

References builderName, edm::EDConsumerBase::consumes(), edm::EDConsumerBase::esConsumes(), fitToken, fname, edm::ParameterSet::getParameter(), LogTrace, maxeta, mineta, propagatorName, smoothToken, sname, srcName, AlCaHLTBitMon_QueryRunRegistry::string, theBuilderToken, theGToken, theMFToken, thePropagatorToken, theTCCollectionToken, and tTopoToken.

                                                              : trackerHitAssociatorConfig_(consumesCollector()) {
   LogTrace("TestSmoothHits") << iConfig << std::endl;
   propagatorName = iConfig.getParameter<std::string>("Propagator");
   builderName = iConfig.getParameter<std::string>("TTRHBuilder");
   srcName = iConfig.getParameter<std::string>("src");
   fname = iConfig.getParameter<std::string>("Fitter");
   sname = iConfig.getParameter<std::string>("Smoother");
   mineta = iConfig.getParameter<double>("mineta");
   maxeta = iConfig.getParameter<double>("maxeta");
 
   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));
   fitToken = esConsumes<edm::Transition::BeginRun>(edm::ESInputTag("", fname));
   smoothToken = esConsumes<edm::Transition::BeginRun>(edm::ESInputTag("", sname));
   tTopoToken = esConsumes();
 
   theTCCollectionToken = consumes(edm::InputTag(srcName));
 }

◆ ~TestSmoothHits()

TestSmoothHits::~TestSmoothHits ( )

override

Definition at line 44 of file TestSmoothHits.cc.

44 {}

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 201 of file TestSmoothHits.cc.

References TransientTrackingRecHitBuilder::build(), TrajectoryStateOnSurface::cartesianError(), heavyFlavorValidationHarvestingSequence_cff::combiner, GlobalErrorBase< T, ErrorWeightType >::cxx(), GlobalErrorBase< T, ErrorWeightType >::cyy(), GlobalErrorBase< T, ErrorWeightType >::czz(), dumpMFGeometry_cfg::delta, DetId::det(), hcalRecHitTable_cff::detId, edm::EventSetup::getHandle(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), hChi2_vs_clsize, hChi2_vs_Process, hChi2Increment, TrackQuality_cff::hitAssociator, hfClusterShapes_cfi::hits, 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, hTotChi2Increment, mps_fire::i, TrackerGeometry::idToDet(), iEvent, sistrip::SpyUtilities::isValid(), TrajectoryStateOnSurface::isValid(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, TrackerTopology::layer(), PSimHit::localPosition(), LogTrace, visualization-live-secondInstance_cfg::m, mag(), CartesianTrajectoryError::matrix(), maxeta, mineta, PSimHit::momentumAtEntry(), CartesianTrajectoryError::position(), PSimHit::processType(), edm::ESHandle< T >::product(), projectHit(), Propagator::propagate(), rpcPointValidation_cfi::recHit, TrackCandidate::recHits(), mps_fire::result, alignCSCRings::s, TrackCandidate::seed(), smooth, mathSSE::sqrt(), GeomDet::surface(), theBuilder, theG, theMF, thePropagator, theTCCollection, theTCCollectionToken, title, Surface::toGlobal(), trackerHitAssociatorConfig_, TrajectorySmoother::trajectories(), TrackCandidate::trajectoryStateOnDet(), trajectoryStateTransform::transientState(), tTopoToken, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                 {
   //Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopo = iSetup.getHandle(tTopoToken);
 
   LogTrace("TestSmoothHits") << "new event" << std::endl;
 
   theTCCollection = iEvent.getHandle(theTCCollectionToken);
   TrackerHitAssociator hitAssociator(iEvent, trackerHitAssociatorConfig_);
 
   TrajectoryStateCombiner combiner;
 
   for (TrackCandidateCollection::const_iterator i = theTCCollection->begin(); i != theTCCollection->end(); i++) {
     LogTrace("TestSmoothHits") << "new candidate" << std::endl;
 
     const TrackCandidate* theTC = &(*i);
     PTrajectoryStateOnDet state = theTC->trajectoryStateOnDet();
 
     //convert PTrajectoryStateOnDet to TrajectoryStateOnSurface
 
     DetId detId(state.detId());
     TrajectoryStateOnSurface theTSOS =
         trajectoryStateTransform::transientState(state, &(theG->idToDet(detId)->surface()), theMF.product());
 
     if (theTSOS.globalMomentum().eta() > maxeta || theTSOS.globalMomentum().eta() < mineta)
       continue;
 
     //convert the TrackingRecHit vector to a TransientTrackingRecHit vector
     TransientTrackingRecHit::RecHitContainer hits;
 
     for (auto const& recHit : theTC->recHits()) {
       hits.push_back(theBuilder->build(&recHit));
     }
 
     //call the fitter
     std::vector<Trajectory> fitted = fit->fit(theTC->seed(), hits, theTSOS);
     //call the smoother
     std::vector<Trajectory> result;
     for (std::vector<Trajectory>::iterator it = fitted.begin(); it != fitted.end(); it++) {
       std::vector<Trajectory> smoothed = smooth->trajectories(*it);
       result.insert(result.end(), smoothed.begin(), smoothed.end());
     }
     if (result.empty())
       continue;
     std::vector<TrajectoryMeasurement> vtm = result[0].measurements();
 
     TSOS lastState = theTSOS;
     for (std::vector<TrajectoryMeasurement>::iterator tm = vtm.begin(); tm != vtm.end(); tm++) {
       TransientTrackingRecHit::ConstRecHitPointer rhit = tm->recHit();
       if ((rhit)->isValid() == 0 && rhit->det() != nullptr)
         continue;
       LogTrace("TestSmoothHits") << "new hit";
 
       int subdetId = rhit->det()->geographicalId().subdetId();
       DetId id = rhit->det()->geographicalId();
       int layerId = tTopo->layer(id);
       LogTrace("TestSmoothHits") << "subdetId=" << subdetId << " layerId=" << layerId;
 
       double delta = 99999;
       LocalPoint rhitLPv = rhit->localPosition();
 
       std::vector<PSimHit> assSimHits = hitAssociator.associateHit(*(rhit->hit()));
       if (assSimHits.empty())
         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();
         }
       }
 
       TSOS currentState = combiner(tm->backwardPredictedState(), tm->forwardPredictedState());
       TSOS updatedState = tm->updatedState();
 
       //plot chi2 increment
       double chi2increment = tm->estimate();
       LogTrace("TestSmoothHits") << "tm->estimate()=" << tm->estimate();
       title.str("");
       title << "Chi2Increment_" << subdetId << "-" << layerId;
       hChi2Increment[title.str()]->Fill(chi2increment);
       hTotChi2Increment->Fill(chi2increment);
       hChi2_vs_Process->Fill(chi2increment, shit.processType());
       if (dynamic_cast<const SiPixelRecHit*>(rhit->hit()))
         hChi2_vs_clsize->Fill(chi2increment, ((const SiPixelRecHit*)(rhit->hit()))->cluster()->size());
       if (dynamic_cast<const SiStripRecHit2D*>(rhit->hit()))
         hChi2_vs_clsize->Fill(chi2increment, ((const SiStripRecHit2D*)(rhit->hit()))->cluster()->amplitudes().size());
 
       //test hits
       const Surface* surf = &((rhit)->det()->surface());
       LocalVector shitLMom;
       LocalPoint shitLPos;
       if (dynamic_cast<const SiStripMatchedRecHit2D*>(rhit->hit())) {
         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 {
         shitLPos = shit.localPosition();
         shitLMom = shit.momentumAtEntry();
       }
       GlobalVector shitGMom = surf->toGlobal(shitLMom);
       GlobalPoint shitGPos = surf->toGlobal(shitLPos);
       //      if (dynamic_cast<const SiStripMatchedRecHit2D*>((rhit)->hit())) {
       //    double rechitmatchedx = rhit->localPosition().x();
       //    double rechitmatchedy = rhit->localPosition().y();
       //    double mindist = 999999;
       //    double 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++){
       //      const PSimHit& hit = *m;
       //      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 = LocalVector( plane.toLocal(stripDet->surface().toGlobal( stripDir)));
       //      std::pair<LocalPoint,LocalVector> hitPair( projectedPos, localStripDir);
       //      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 {
       //    shitLPos = shit.localPosition();
       //    shitLMom = shit.momentumAtEntry();
       //       }
       //       GlobalVector shitGMom = surf->toGlobal(shitLMom);
       //       GlobalPoint shitGPos = surf->toGlobal(shitLPos);
 
       GlobalVector tsosGMom = currentState.globalMomentum();
       GlobalError tsosGMEr(currentState.cartesianError().matrix().Sub<AlgebraicSymMatrix33>(3, 3));
       GlobalPoint tsosGPos = currentState.globalPosition();
       GlobalError tsosGPEr = currentState.cartesianError().position();
 
       GlobalPoint rhitGPos = (rhit)->globalPosition();
       GlobalError rhitGPEr = (rhit)->globalPositionError();
 
       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("TestSmoothHits") << "rs" << std::endl;
 
       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("TestSmoothHits") << "tr" << std::endl;
 
       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("TestSmoothHits") << "ts1" << std::endl;
 
       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("TestSmoothHits") << "ts2" << std::endl;
 
       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())) {
         //mono
         LogTrace("TestSmoothHits") << "MONO HIT" << std::endl;
         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 = *(assSimHits.begin());
         const Surface* monoSurf = &(tMonoHit->det()->surface());
         if (monoSurf == nullptr)
           continue;
         TSOS monoState = thePropagator->propagate(lastState, *monoSurf);
         if (monoState.isValid() == 0)
           continue;
 
         LocalVector monoShitLMom = sMonoHit.momentumAtEntry();
         GlobalVector monoShitGMom = monoSurf->toGlobal(monoShitLMom);
         LocalPoint monoShitLPos = sMonoHit.localPosition();
         GlobalPoint monoShitGPos = monoSurf->toGlobal(monoShitLPos);
 
         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());
 
         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("TestSmoothHits") << "STEREO HIT" << std::endl;
         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 = *(assSimHits.begin());
         const Surface* stereoSurf = &(tStereoHit->det()->surface());
         if (stereoSurf == nullptr)
           continue;
         TSOS stereoState = thePropagator->propagate(lastState, *stereoSurf);
         if (stereoState.isValid() == 0)
           continue;
 
         LocalVector stereoShitLMom = sStereoHit.momentumAtEntry();
         GlobalVector stereoShitGMom = stereoSurf->toGlobal(stereoShitLMom);
         LocalPoint stereoShitLPos = sStereoHit.localPosition();
         GlobalPoint stereoShitGPos = stereoSurf->toGlobal(stereoShitLPos);
 
         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();
 
         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());
         //double pullGPY_rs_stereo = (stereoRhitGPos.y()-stereoShitGPos.y());
         //double pullGPZ_rs_stereo = (stereoRhitGPos.z()-stereoShitGPos.z());
 
         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);
       }
       lastState = updatedState;
       //#endif
     }
   }
   LogTrace("TestSmoothHits") << "end of event" << std::endl;
 }

◆ beginRun()

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

overrideprivate

Definition at line 46 of file TestSmoothHits.cc.

References file, fitToken, edm::EventSetup::getHandle(), hChi2_vs_clsize, hChi2_vs_Process, hChi2Increment, 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, hTotChi2Increment, mps_fire::i, dqmiolumiharvest::j, smooth, smoothToken, theBuilder, theBuilderToken, theG, theGToken, theMF, theMFToken, thePropagator, thePropagatorToken, and title.

                                                                             {
   theG = iSetup.getHandle(theGToken);
   theMF = iSetup.getHandle(theMFToken);
   thePropagator = iSetup.getHandle(thePropagatorToken);
   theBuilder = iSetup.getHandle(theBuilderToken);
   fit = iSetup.getHandle(fitToken);
   smooth = iSetup.getHandle(smoothToken);
 
   file = new TFile("testSmoothHits.root", "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 << "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 << "Chi2Increment_" << i + 1 << "-" << j + 1;
       hChi2Increment[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, 0, 100);
 
       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 << "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 << "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);
   hChi2_vs_Process = new TH2F("Chi2_vs_Process", "Chi2_vs_Process", 1000, 0, 100, 17, -0.5, 16.5);
   hChi2_vs_clsize = new TH2F("Chi2_vs_clsize", "Chi2_vs_clsize", 1000, 0, 100, 17, -0.5, 16.5);
 }

◆ endJob()

void TestSmoothHits::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 651 of file TestSmoothHits.cc.

References file, hChi2_vs_clsize, hChi2_vs_Process, hChi2Increment, 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, hTotChi2Increment, mps_fire::i, dqmiolumiharvest::j, and title.

                             {
   //file->Write();
   TDirectory* chi2i = file->mkdir("Chi2_Increment");
 
   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");
 
   chi2i->cd();
   hTotChi2Increment->Write();
   hChi2_vs_Process->Write();
   hChi2_vs_clsize->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;
       chi2i->cd();
       title.str("");
       title << "Chi2Increment_" << i + 1 << "-" << j + 1;
       hChi2Increment[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)) {
         //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 TestSmoothHits::endRun	(	edm::Run const &	run,
		const edm::EventSetup &
	)

inlineoverrideprivate

Definition at line 58 of file TestSmoothHits.h.

58 {}

◆ projectHit()

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

private

Definition at line 900 of file TestSmoothHits.cc.

References funct::cos(), DeadROC_duringRun::dir, isotrackNtupler::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 TestSmoothHits::builderName

private

Definition at line 68 of file TestSmoothHits.h.

Referenced by TestSmoothHits().

◆ file

TFile* TestSmoothHits::file

private

Definition at line 90 of file TestSmoothHits.h.

Referenced by beginRun(), and endJob().

◆ fit

edm::ESHandle<TrajectoryFitter> TestSmoothHits::fit

private

Definition at line 78 of file TestSmoothHits.h.

Referenced by trackingPlots.Iteration::modules().

◆ fitToken

edm::ESGetToken<TrajectoryFitter, TrajectoryFitter::Record> TestSmoothHits::fitToken

private

Definition at line 86 of file TestSmoothHits.h.

Referenced by beginRun(), and TestSmoothHits().

◆ fname

std::string TestSmoothHits::fname

private

Definition at line 70 of file TestSmoothHits.h.

Referenced by TestSmoothHits().

◆ hChi2_vs_clsize

TH2F * TestSmoothHits::hChi2_vs_clsize

private

Definition at line 106 of file TestSmoothHits.h.

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

◆ hChi2_vs_Process

TH2F* TestSmoothHits::hChi2_vs_Process

private

Definition at line 106 of file TestSmoothHits.h.

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

◆ hChi2Increment

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

private

Definition at line 104 of file TestSmoothHits.h.

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

◆ hPullGM_X_ts

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

private

Definition at line 95 of file TestSmoothHits.h.

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

◆ hPullGM_X_ts_mono

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

private

Definition at line 111 of file TestSmoothHits.h.

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

◆ hPullGM_X_ts_stereo

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

private

Definition at line 124 of file TestSmoothHits.h.

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

◆ hPullGM_Y_ts

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

private

Definition at line 96 of file TestSmoothHits.h.

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

◆ hPullGM_Y_ts_mono

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

private

Definition at line 112 of file TestSmoothHits.h.

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

◆ hPullGM_Y_ts_stereo

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

private

Definition at line 125 of file TestSmoothHits.h.

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

◆ hPullGM_Z_ts

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

private

Definition at line 97 of file TestSmoothHits.h.

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

◆ hPullGM_Z_ts_mono

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

private

Definition at line 113 of file TestSmoothHits.h.

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

◆ hPullGM_Z_ts_stereo

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

private

Definition at line 126 of file TestSmoothHits.h.

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

◆ hPullGP_X_rs

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

private

Definition at line 98 of file TestSmoothHits.h.

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

◆ hPullGP_X_rs_mono

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

private

Definition at line 114 of file TestSmoothHits.h.

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

◆ hPullGP_X_rs_stereo

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

private

Definition at line 127 of file TestSmoothHits.h.

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

◆ hPullGP_X_tr

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

private

Definition at line 101 of file TestSmoothHits.h.

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

◆ hPullGP_X_tr_mono

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

private

Definition at line 117 of file TestSmoothHits.h.

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

◆ hPullGP_X_tr_stereo

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

private

Definition at line 130 of file TestSmoothHits.h.

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

◆ hPullGP_X_ts

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

private

Definition at line 92 of file TestSmoothHits.h.

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

◆ hPullGP_X_ts_mono

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

private

Definition at line 108 of file TestSmoothHits.h.

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

◆ hPullGP_X_ts_stereo

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

private

Definition at line 121 of file TestSmoothHits.h.

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

◆ hPullGP_Y_rs

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

private

Definition at line 99 of file TestSmoothHits.h.

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

◆ hPullGP_Y_rs_mono

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

private

Definition at line 115 of file TestSmoothHits.h.

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

◆ hPullGP_Y_rs_stereo

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

private

Definition at line 128 of file TestSmoothHits.h.

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

◆ hPullGP_Y_tr

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

private

Definition at line 102 of file TestSmoothHits.h.

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

◆ hPullGP_Y_tr_mono

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

private

Definition at line 118 of file TestSmoothHits.h.

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

◆ hPullGP_Y_tr_stereo

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

private

Definition at line 131 of file TestSmoothHits.h.

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

◆ hPullGP_Y_ts

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

private

Definition at line 93 of file TestSmoothHits.h.

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

◆ hPullGP_Y_ts_mono

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

private

Definition at line 109 of file TestSmoothHits.h.

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

◆ hPullGP_Y_ts_stereo

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

private

Definition at line 122 of file TestSmoothHits.h.

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

◆ hPullGP_Z_rs

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

private

Definition at line 100 of file TestSmoothHits.h.

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

◆ hPullGP_Z_rs_mono

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

private

Definition at line 116 of file TestSmoothHits.h.

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

◆ hPullGP_Z_rs_stereo

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

private

Definition at line 129 of file TestSmoothHits.h.

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

◆ hPullGP_Z_tr

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

private

Definition at line 103 of file TestSmoothHits.h.

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

◆ hPullGP_Z_tr_mono

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

private

Definition at line 119 of file TestSmoothHits.h.

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

◆ hPullGP_Z_tr_stereo

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

private

Definition at line 132 of file TestSmoothHits.h.

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

◆ hPullGP_Z_ts

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

private

Definition at line 94 of file TestSmoothHits.h.

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

◆ hPullGP_Z_ts_mono

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

private

Definition at line 110 of file TestSmoothHits.h.

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

◆ hPullGP_Z_ts_stereo

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

private

Definition at line 123 of file TestSmoothHits.h.

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

◆ hTotChi2Increment

TH1F* TestSmoothHits::hTotChi2Increment

private

Definition at line 105 of file TestSmoothHits.h.

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

◆ maxeta

double TestSmoothHits::maxeta

private

Definition at line 65 of file TestSmoothHits.h.

Referenced by analyze(), and TestSmoothHits().

◆ mineta

double TestSmoothHits::mineta

private

Definition at line 65 of file TestSmoothHits.h.

Referenced by analyze(), and TestSmoothHits().

◆ propagatorName

std::string TestSmoothHits::propagatorName

private

Definition at line 67 of file TestSmoothHits.h.

Referenced by TestSmoothHits().

◆ smooth

edm::ESHandle<TrajectorySmoother> TestSmoothHits::smooth

private

Definition at line 79 of file TestSmoothHits.h.

Referenced by analyze(), and beginRun().

◆ smoothToken

edm::ESGetToken<TrajectorySmoother, TrajectoryFitter::Record> TestSmoothHits::smoothToken

private

Definition at line 87 of file TestSmoothHits.h.

Referenced by beginRun(), and TestSmoothHits().

◆ sname

std::string TestSmoothHits::sname

private

Definition at line 71 of file TestSmoothHits.h.

Referenced by TestSmoothHits().

◆ srcName

std::string TestSmoothHits::srcName

private

Definition at line 69 of file TestSmoothHits.h.

Referenced by TestSmoothHits().

◆ theBuilder

edm::ESHandle<TransientTrackingRecHitBuilder> TestSmoothHits::theBuilder

private

Definition at line 76 of file TestSmoothHits.h.

Referenced by analyze(), and beginRun().

◆ theBuilderToken

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

private

Definition at line 84 of file TestSmoothHits.h.

Referenced by beginRun(), and TestSmoothHits().

◆ theG

edm::ESHandle<TrackerGeometry> TestSmoothHits::theG

private

Definition at line 73 of file TestSmoothHits.h.

Referenced by analyze(), and beginRun().

◆ theGToken

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

private

Definition at line 81 of file TestSmoothHits.h.

Referenced by beginRun(), and TestSmoothHits().

◆ theMF

edm::ESHandle<MagneticField> TestSmoothHits::theMF

private

Definition at line 74 of file TestSmoothHits.h.

Referenced by analyze(), and beginRun().

◆ theMFToken

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

private

Definition at line 82 of file TestSmoothHits.h.

Referenced by beginRun(), and TestSmoothHits().

◆ thePropagator

edm::ESHandle<Propagator> TestSmoothHits::thePropagator

private

Definition at line 75 of file TestSmoothHits.h.

Referenced by analyze(), and beginRun().

◆ thePropagatorToken

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

private

Definition at line 83 of file TestSmoothHits.h.

Referenced by beginRun(), and TestSmoothHits().

◆ theTCCollection

edm::Handle<TrackCandidateCollection> TestSmoothHits::theTCCollection

private

Definition at line 77 of file TestSmoothHits.h.

Referenced by analyze().

◆ theTCCollectionToken

edm::EDGetTokenT<TrackCandidateCollection> TestSmoothHits::theTCCollectionToken

private

Definition at line 85 of file TestSmoothHits.h.

Referenced by analyze(), and TestSmoothHits().

◆ title

std::stringstream TestSmoothHits::title

private

Definition at line 91 of file TestSmoothHits.h.

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

◆ trackerHitAssociatorConfig_

TrackerHitAssociator::Config TestSmoothHits::trackerHitAssociatorConfig_

private

Definition at line 63 of file TestSmoothHits.h.

Referenced by analyze().

◆ tTopoToken

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> TestSmoothHits::tTopoToken

private

Definition at line 88 of file TestSmoothHits.h.

Referenced by analyze(), and TestSmoothHits().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ TestSmoothHits()

◆ ~TestSmoothHits()

Member Function Documentation

◆ analyze()

◆ beginRun()

◆ endJob()

◆ endRun()

◆ projectHit()

Member Data Documentation

◆ builderName

◆ file

◆ fit

◆ fitToken

◆ fname

◆ hChi2_vs_clsize

◆ hChi2_vs_Process

◆ hChi2Increment

◆ 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

◆ hTotChi2Increment

◆ maxeta

◆ mineta

◆ propagatorName

◆ smooth

◆ smoothToken

◆ sname

◆ srcName

◆ theBuilder

◆ theBuilderToken

◆ theG

◆ theGToken

◆ theMF

◆ theMFToken

◆ thePropagator

◆ thePropagatorToken

◆ theTCCollection

◆ theTCCollectionToken

◆ title

◆ trackerHitAssociatorConfig_

◆ tTopoToken