#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 final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

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

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

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

std::vector< ESProxyIndex > 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

bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices 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 ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

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

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

template<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(), 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(), 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, L1EGammaClusterEmuProducer_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 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