#include <HitResol.h>

Inheritance diagram for HitResol:

Public Types
typedef std::vector< Trajectory >	TrajectoryCollection

Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
double	checkConsistency (const StripClusterParameterEstimator::LocalValues &parameters, double xx, double xerr)

bool	getPairParameters (const MagneticField magField_, AnalyticalPropagator &propagator, const TrajectoryMeasurement traj1, const TrajectoryMeasurement *traj2, float &pairPath, float &hitDX, float &trackDX, float &trackDXE, float &trackParamX, float &trackParamY, float &trackParamDXDZ, float &trackParamDYDZ, float &trackParamXE, float &trackParamYE, float &trackParamDXDZE, float &trackParamDYDZE)

void	getSimHitRes (const GeomDetUnit det, const LocalVector &trackdirection, const TrackingRecHit &recHit, float &trackWidth, float pitch, LocalVector &drift)

double	getSimpleRes (const TrajectoryMeasurement *traj1)

	HitResol (const edm::ParameterSet &conf)

	~HitResol () override=default

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

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

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

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

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

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

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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)

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

void	beginJob () override

void	endJob () override

Private Attributes
const bool	addLumi_

float	atEdge

float	atEdge_2

unsigned int	clusterCharge

unsigned int	clusterCharge_2

unsigned int	clusterWidth

unsigned int	clusterWidth_2

const edm::EDGetTokenT< reco::TrackCollection >	combinatorialTracks_token_

const int	compSettings_

const edm::ESGetToken< StripClusterParameterEstimator, TkStripCPERecord >	cpeToken_

const bool	cutOnTracks_

const bool	DEBUG_

float	driftAlpha

float	driftAlpha_2

int	events

int	EventTrackCKF

float	expWidth

float	expWidth_2

const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	geomToken_

std::map< TString, TH2F * >	histos2d_

float	hitDX

unsigned int	iidd1

unsigned int	iidd2

const unsigned int	layers_

const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	magFieldToken_

const double	momentumCut_

float	mymom

float	mypitch1

float	mypitch2

float	N1

float	N1uProj

float	N2

float	N2uProj

int	Nstrips

int	Nstrips_2

int	NumberOf_tracks

int	numHits

float	pairPath

unsigned int	pairsOnly

float	ProbTrackChi2

TTree *	reso

const edm::EDGetTokenT< LumiScalersCollection >	scalerToken_

float	simpleRes

const edm::ESGetToken< SiStripQuality, SiStripQualityRcd >	siStripQualityToken_

float	StripCPE1_smp_pos_error

float	StripCPE2_smp_pos_error

float	StripErrorSquared1

float	StripErrorSquared2

float	thickness

float	thickness_2

const edm::EDGetTokenT< std::vector< Trajectory > >	tjToken_

const edm::EDGetTokenT< reco::TrackCollection >	tkToken_

const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	topoToken_

float	track_eta

float	track_momentum

float	track_phi

float	track_pt

float	track_trackChi2

float	track_width

float	trackDX

float	trackDXE

const unsigned int	trackMultiplicityCut_

float	trackParamDXDZ

float	trackParamDXDZE

float	trackParamDYDZ

float	trackParamDYDZE

float	trackParamX

float	trackParamXE

float	trackParamY

float	trackParamYE

float	trackWidth

float	trackWidth_2

TTree *	treso

float	uerr2

float	uerr2_2

const unsigned int	usePairsOnly_

Static Private Attributes
static constexpr float	cmToUm = 10000.f

Additional Inherited Members
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

Definition at line 33 of file HitResol.h.

Member Typedef Documentation

◆ TrajectoryCollection

typedef std::vector<Trajectory> HitResol::TrajectoryCollection

Definition at line 62 of file HitResol.h.

Constructor & Destructor Documentation

◆ HitResol()

HitResol::HitResol ( const edm::ParameterSet & conf )

explicit

Definition at line 71 of file HitResol.cc.

References TFileService::kSharedResource.

     : scalerToken_(consumes<LumiScalersCollection>(conf.getParameter<edm::InputTag>("lumiScalers"))),
       combinatorialTracks_token_(
           consumes<reco::TrackCollection>(conf.getParameter<edm::InputTag>("combinatorialTracks"))),
       tjToken_(consumes<std::vector<Trajectory> >(conf.getParameter<edm::InputTag>("trajectories"))),
       topoToken_(esConsumes()),
       geomToken_(esConsumes()),
       cpeToken_(esConsumes(edm::ESInputTag("", "StripCPEfromTrackAngle"))),
       siStripQualityToken_(esConsumes()),
       magFieldToken_(esConsumes()),
       addLumi_(conf.getUntrackedParameter<bool>("addLumi", false)),
       DEBUG_(conf.getParameter<bool>("Debug")),
       cutOnTracks_(conf.getUntrackedParameter<bool>("cutOnTracks", false)),
       momentumCut_(conf.getUntrackedParameter<double>("MomentumCut", 3.)),
       compSettings_(conf.getUntrackedParameter<int>("CompressionSettings", -1)),
       usePairsOnly_(conf.getUntrackedParameter<unsigned int>("UsePairsOnly", 1)),
       layers_(conf.getParameter<int>("Layer")),
       trackMultiplicityCut_(conf.getUntrackedParameter<unsigned int>("trackMultiplicity", 100)) {
   usesResource(TFileService::kSharedResource);
 }

◆ ~HitResol()

HitResol::~HitResol ( )

overridedefault

Member Function Documentation

◆ analyze()

void HitResol::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 159 of file HitResol.cc.

References AnalyticalPropagator_cfi::AnalyticalPropagator, anyDirection, atEdge, atEdge_2, ChiSquaredProbability(), clusterWidth, clusterWidth_2, combinatorialTracks_token_, cpeToken_, shallow::drift(), StripClusterParameterEstimator::driftDirection(), MillePedeFileConverter_cfg::e, events, expWidth, expWidth_2, geomToken_, edm::EventSetup::getData(), edm::EventSetup::getHandle(), getPairParameters(), getSimHitRes(), getSimpleRes(), TrackerTopology::glued(), histos2d_, hitDX, globals_cff::id1, globals_cff::id2, TrackerGeometry::idToDetUnit(), iidd1, iidd2, TrackerTopology::isStereo(), TrajectoryStateOnSurface::localDirection(), LogDebug, magFieldToken_, momentumCut_, mymom, mypitch1, StripTopology::nstrips(), Nstrips, Nstrips_2, numHits, pairPath, pairsOnly, ProbTrackChi2, edm::Handle< T >::product(), reso, simpleRes, siStripQualityToken_, findQualityFiles::size, StripSubdetector::TEC, StripSubdetector::TIB, tjToken_, topoToken_, track_eta, track_momentum, track_phi, track_pt, track_trackChi2, trackDX, trackDXE, PbPb_ZMuSkimMuonDPG_cff::tracker, trackParamDXDZ, trackParamDXDZE, trackParamDYDZ, trackParamDYDZE, trackParamX, trackParamXE, trackParamY, trackParamYE, treso, and usePairsOnly_.

                                                                  {
   //Retrieve tracker topology from geometry
   const TrackerTopology* tTopo = &es.getData(topoToken_);
 
   LogDebug("SiStripHitResolution:HitResol") << "beginning analyze from HitResol" << endl;
 
   using namespace edm;
   using namespace reco;
 
   // Step A: Get Inputs
 
   int run_nr = e.id().run();
   int ev_nr = e.id().event();
 
   // get the tracks
   edm::Handle<reco::TrackCollection> trackCollectionCKF;
   e.getByToken(combinatorialTracks_token_, trackCollectionCKF);
   const reco::TrackCollection* tracksCKF = trackCollectionCKF.product();
 
   // get the trajectory collection
   edm::Handle<std::vector<Trajectory> > trajectoryCollectionHandle;
   e.getByToken(tjToken_, trajectoryCollectionHandle);
   const TrajectoryCollection* trajectoryCollection = trajectoryCollectionHandle.product();
 
   //get tracker geometry
   edm::ESHandle<TrackerGeometry> tracker = es.getHandle(geomToken_);
   const TrackerGeometry* tkgeom = &(*tracker);
 
   //get Cluster Parameter Estimator
   edm::ESHandle<StripClusterParameterEstimator> parameterestimator = es.getHandle(cpeToken_);
   const StripClusterParameterEstimator& stripcpe(*parameterestimator);
 
   // get the SiStripQuality records
   edm::ESHandle<SiStripQuality> SiStripQuality_ = es.getHandle(siStripQualityToken_);
 
   // get the magnetic field
   const MagneticField* magField_ = &es.getData(magFieldToken_);
 
   events++;
 
   // List of variables for SiStripHitResolution ntuple
   mymom = 0;
   numHits = 0;
   ProbTrackChi2 = 0;
   iidd1 = 0;
   mypitch1 = 0;
   clusterWidth = 0;
   expWidth = 0;
   atEdge = 0;
   simpleRes = 0;
   iidd2 = 0;
   clusterWidth_2 = 0;
   expWidth_2 = 0;
   atEdge_2 = 0;
   pairPath = 0;
   hitDX = 0;
   trackDX = 0;
   trackDXE = 0;
   trackParamX = 0;
   trackParamY = 0;
   trackParamDXDZ = 0;
   trackParamDYDZ = 0;
   trackParamXE = 0;
   trackParamYE = 0;
   trackParamDXDZE = 0;
   trackParamDYDZE = 0;
   pairsOnly = 0;
 
   LogDebug("HitResol") << "Starting analysis, nrun nevent, tracksCKF->size(): " << run_nr << " " << ev_nr << " "
                        << tracksCKF->size() << std::endl;
 
   for (unsigned int iT = 0; iT < tracksCKF->size(); ++iT) {
     track_momentum = tracksCKF->at(iT).p();
     track_pt = tracksCKF->at(iT).p();
     track_eta = tracksCKF->at(iT).eta();
     track_phi = tracksCKF->at(iT).phi();
     track_trackChi2 = ChiSquaredProbability((double)(tracksCKF->at(iT).chi2()), (double)(tracksCKF->at(iT).ndof()));
     treso->Fill();
   }
 
   histos2d_["track_phi_vs_eta"]->Fill(track_phi, track_eta);
 
   // loop over trajectories from refit
   for (const auto& traj : *trajectoryCollection) {
     const auto& TMeas = traj.measurements();
     // Loop on each measurement and take it into consideration
     //--------------------------------------------------------
     for (auto itm = TMeas.cbegin(); itm != TMeas.cend(); ++itm) {
       if (!itm->updatedState().isValid()) {
         LogDebug("HitResol") << "trajectory measurement not valid" << std::endl;
         continue;
       }
 
       const TransientTrackingRecHit::ConstRecHitPointer mypointhit = itm->recHit();
       const TrackingRecHit* myhit = (*itm->recHit()).hit();
 
       ProbTrackChi2 = 0;
       numHits = 0;
 
       LogDebug("HitResol") << "TrackChi2 =  "
                            << ChiSquaredProbability((double)(traj.chiSquared()), (double)(traj.ndof(false))) << "\n"
                            << "itm->updatedState().globalMomentum().perp(): "
                            << itm->updatedState().globalMomentum().perp() << "\n"
                            << "numhits " << traj.foundHits() << std::endl;
 
       numHits = traj.foundHits();
       ProbTrackChi2 = ChiSquaredProbability((double)(traj.chiSquared()), (double)(traj.ndof(false)));
 
       mymom = itm->updatedState().globalMomentum().perp();
 
       //Now for the first hit
       TrajectoryStateOnSurface mytsos = itm->updatedState();
       const auto hit1 = itm->recHit();
       DetId id1 = hit1->geographicalId();
       if (id1.subdetId() < StripSubdetector::TIB || id1.subdetId() > StripSubdetector::TEC)
         continue;
 
       if (hit1->isValid() && mymom > momentumCut_ &&
           (id1.subdetId() >= StripSubdetector::TIB && id1.subdetId() <= StripSubdetector::TEC)) {
         const auto stripdet = dynamic_cast<const StripGeomDetUnit*>(tkgeom->idToDetUnit(hit1->geographicalId()));
         const StripTopology& Topo = stripdet->specificTopology();
         int Nstrips = Topo.nstrips();
         mypitch1 = stripdet->surface().bounds().width() / Topo.nstrips();
 
         const auto det = dynamic_cast<const StripGeomDetUnit*>(tkgeom->idToDetUnit(mypointhit->geographicalId()));
 
         TrajectoryStateOnSurface mytsos = itm->updatedState();
         LocalVector trackDirection = mytsos.localDirection();
         LocalVector drift = stripcpe.driftDirection(stripdet);
 
         const auto hit1d = dynamic_cast<const SiStripRecHit1D*>(myhit);
 
         if (hit1d) {
           getSimHitRes(det, trackDirection, *hit1d, expWidth, &mypitch1, drift);
           clusterWidth = hit1d->cluster()->amplitudes().size();
           uint16_t firstStrip = hit1d->cluster()->firstStrip();
           uint16_t lastStrip = firstStrip + (hit1d->cluster()->amplitudes()).size() - 1;
           atEdge = (firstStrip == 0 || lastStrip == (Nstrips - 1));
         }
 
         const auto hit2d = dynamic_cast<const SiStripRecHit2D*>(myhit);
 
         if (hit2d) {
           getSimHitRes(det, trackDirection, *hit2d, expWidth, &mypitch1, drift);
           clusterWidth = hit2d->cluster()->amplitudes().size();
           uint16_t firstStrip = hit2d->cluster()->firstStrip();
           uint16_t lastStrip = firstStrip + (hit2d->cluster()->amplitudes()).size() - 1;
           atEdge = (firstStrip == 0 || lastStrip == (Nstrips - 1));
         }
 
         simpleRes =
             getSimpleRes(&(*itm));  // simple resolution by using the track re-fit forward and backward predicted state
 
         histos2d_["residual_vs_trackMomentum"]->Fill(itm->updatedState().globalMomentum().mag(),
                                                      simpleRes * 10000);   // reso in cm *10000 == micro-meter
         histos2d_["residual_vs_trackPt"]->Fill(mymom, simpleRes * 10000);  // reso in cm *10000 == micro-meter
         histos2d_["residual_vs_trackEta"]->Fill(itm->updatedState().globalMomentum().eta(), simpleRes * 10000);
         histos2d_["residual_vs_trackPhi"]->Fill(itm->updatedState().globalMomentum().phi(), simpleRes * 10000);
         histos2d_["residual_vs_expectedWidth"]->Fill(expWidth, simpleRes * 10000);
         histos2d_["numHits_vs_residual"]->Fill(simpleRes * 10000, numHits);
 
         // Now to see if there is a match - pair method - hit in overlapping sensors
         vector<TrajectoryMeasurement>::const_iterator itTraj2 = TMeas.end();  // last hit along the fitted track
 
         for (auto itmCompare = itm - 1;
              // start to compare from the 5th hit
              itmCompare >= TMeas.cbegin() && itmCompare > itm - 4;
              --itmCompare) {
           const auto hit2 = itmCompare->recHit();
           if (!hit2->isValid())
             continue;
           DetId id2 = hit2->geographicalId();
 
           //must be from the same detector and layer
           iidd1 = hit1->geographicalId().rawId();
           iidd2 = hit2->geographicalId().rawId();
           if (id1.subdetId() != id2.subdetId() || ::checkLayer(iidd1, tTopo) != ::checkLayer(iidd2, tTopo))
             break;
           //must both be stereo if one is
           if (tTopo->isStereo(id1) != tTopo->isStereo(id2))
             continue;
           //A check i dont completely understand but might as well keep there
           if (tTopo->glued(id1) == id1.rawId())
             LogDebug("HitResol") << "BAD GLUED: Have glued layer with id = " << id1.rawId()
                                  << " and glued id = " << tTopo->glued(id1) << "  and stereo = " << tTopo->isStereo(id1)
                                  << endl;
           if (tTopo->glued(id2) == id2.rawId())
             LogDebug("HitResol") << "BAD GLUED: Have glued layer with id = " << id2.rawId()
                                  << " and glued id = " << tTopo->glued(id2) << "  and stereo = " << tTopo->isStereo(id2)
                                  << endl;
 
           itTraj2 = itmCompare;
           break;
         }
 
         if (itTraj2 == TMeas.cend()) {
         } else {
           LogDebug("HitResol") << "Found overlapping sensors " << std::endl;
           pairsOnly = usePairsOnly_;
 
           //We found one....let's fill in the truth info!
           TrajectoryStateOnSurface tsos_2 = itTraj2->updatedState();
           LocalVector trackDirection_2 = tsos_2.localDirection();
           const auto myhit2 = itTraj2->recHit();
           const auto myhit_2 = (*itTraj2->recHit()).hit();
           const auto stripdet_2 = dynamic_cast<const StripGeomDetUnit*>(tkgeom->idToDetUnit(myhit2->geographicalId()));
           const StripTopology& Topo_2 = stripdet_2->specificTopology();
           int Nstrips_2 = Topo_2.nstrips();
           float mypitch_2 = stripdet_2->surface().bounds().width() / Topo_2.nstrips();
 
           if (mypitch1 != mypitch_2)
             return;  // for PairsOnly
 
           const auto det_2 = dynamic_cast<const StripGeomDetUnit*>(tkgeom->idToDetUnit(myhit2->geographicalId()));
 
           LocalVector drift_2 = stripcpe.driftDirection(stripdet_2);
 
           const auto hit1d_2 = dynamic_cast<const SiStripRecHit1D*>(myhit_2);
           if (hit1d_2) {
             getSimHitRes(det_2, trackDirection_2, *hit1d_2, expWidth_2, &mypitch_2, drift_2);
             clusterWidth_2 = hit1d_2->cluster()->amplitudes().size();
             uint16_t firstStrip_2 = hit1d_2->cluster()->firstStrip();
             uint16_t lastStrip_2 = firstStrip_2 + (hit1d_2->cluster()->amplitudes()).size() - 1;
             atEdge_2 = (firstStrip_2 == 0 || lastStrip_2 == (Nstrips_2 - 1));
           }
 
           const auto hit2d_2 = dynamic_cast<const SiStripRecHit2D*>(myhit_2);
           if (hit2d_2) {
             getSimHitRes(det_2, trackDirection_2, *hit2d_2, expWidth_2, &mypitch_2, drift_2);
             clusterWidth_2 = hit2d_2->cluster()->amplitudes().size();
             uint16_t firstStrip_2 = hit2d_2->cluster()->firstStrip();
             uint16_t lastStrip_2 = firstStrip_2 + (hit2d_2->cluster()->amplitudes()).size() - 1;
             atEdge_2 = (firstStrip_2 == 0 || lastStrip_2 == (Nstrips_2 - 1));
           }
 
           // if(pairsOnly && (pitch != pitch2) ) fill = false;
 
           // Make AnalyticalPropagator to use in getPairParameters
           AnalyticalPropagator mypropagator(magField_, anyDirection);
 
           if (!getPairParameters(&(*magField_),
                                  mypropagator,
                                  &(*itTraj2),
                                  &(*itm),
                                  pairPath,
                                  hitDX,
                                  trackDX,
                                  trackDXE,
                                  trackParamX,
                                  trackParamY,
                                  trackParamDXDZ,
                                  trackParamDYDZ,
                                  trackParamXE,
                                  trackParamYE,
                                  trackParamDXDZE,
                                  trackParamDYDZE)) {
           } else {
             LogDebug("HitResol") << " \n\n\n"
                                  << " momentum       " << mymom << "\n"
                                  << " numHits        " << numHits << "\n"
                                  << " trackChi2      " << ProbTrackChi2 << "\n"
                                  << " detID1         " << iidd1 << "\n"
                                  << " pitch1         " << mypitch1 << "\n"
                                  << " clusterW1      " << clusterWidth << "\n"
                                  << " expectedW1     " << expWidth << "\n"
                                  << " atEdge1        " << atEdge << "\n"
                                  << " simpleRes      " << simpleRes << "\n"
                                  << " detID2         " << iidd2 << "\n"
                                  << " clusterW2      " << clusterWidth_2 << "\n"
                                  << " expectedW2     " << expWidth_2 << "\n"
                                  << " atEdge2        " << atEdge_2 << "\n"
                                  << " pairPath       " << pairPath << "\n"
                                  << " hitDX          " << hitDX << "\n"
                                  << " trackDX        " << trackDX << "\n"
                                  << " trackDXE       " << trackDXE << "\n"
                                  << " trackParamX      " << trackParamX << "\n"
                                  << " trackParamY      " << trackParamY << "\n"
                                  << " trackParamDXDZ " << trackParamDXDZ << "\n"
                                  << " trackParamDYDZ " << trackParamDYDZ << "\n"
                                  << " trackParamXE   " << trackParamXE << "\n"
                                  << " trackParamYE   " << trackParamYE << "\n"
                                  << " trackParamDXDZE" << trackParamDXDZE << "\n"
                                  << " trackParamDYDZE" << trackParamDYDZE << std::endl;
             reso->Fill();
           }
         }  //itTraj2 != TMeas.end()
       }    //hit1->isValid()....
     }      // itm
   }        // it
 }

◆ beginJob()

void HitResol::beginJob ( )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 92 of file HitResol.cc.

References atEdge, atEdge_2, clusterWidth, clusterWidth_2, compSettings_, events, EventTrackCKF, expWidth, expWidth_2, compareTotals::fs, histos2d_, hitDX, iidd1, iidd2, mymom, mypitch1, NumberOf_tracks, numHits, pairPath, pairsOnly, ProbTrackChi2, reso, simpleRes, track_eta, track_momentum, track_phi, track_pt, track_trackChi2, trackDX, trackDXE, trackParamDXDZ, trackParamDXDZE, trackParamDYDZ, trackParamDYDZE, trackParamX, trackParamXE, trackParamY, trackParamYE, and treso.

                         {
   edm::Service<TFileService> fs;
   if (compSettings_ > 0) {
     edm::LogInfo("SiStripHitResolution:HitResol") << "the compressions settings are:" << compSettings_ << std::endl;
     fs->file().SetCompressionSettings(compSettings_);
   }
 
   reso = fs->make<TTree>("reso", "tree hit pairs for resolution studies");
   reso->Branch("momentum", &mymom, "momentum/F");
   reso->Branch("numHits", &numHits, "numHits/I");
   reso->Branch("trackChi2", &ProbTrackChi2, "trackChi2/F");
   reso->Branch("detID1", &iidd1, "detID1/I");
   reso->Branch("pitch1", &mypitch1, "pitch1/F");
   reso->Branch("clusterW1", &clusterWidth, "clusterW1/I");
   reso->Branch("expectedW1", &expWidth, "expectedW1/F");
   reso->Branch("atEdge1", &atEdge, "atEdge1/F");
   reso->Branch("simpleRes", &simpleRes, "simpleRes/F");
   reso->Branch("detID2", &iidd2, "detID2/I");
   reso->Branch("clusterW2", &clusterWidth_2, "clusterW2/I");
   reso->Branch("expectedW2", &expWidth_2, "expectedW2/F");
   reso->Branch("atEdge2", &atEdge_2, "atEdge2/F");
   reso->Branch("pairPath", &pairPath, "pairPath/F");
   reso->Branch("hitDX", &hitDX, "hitDX/F");
   reso->Branch("trackDX", &trackDX, "trackDX/F");
   reso->Branch("trackDXE", &trackDXE, "trackDXE/F");
   reso->Branch("trackParamX", &trackParamX, "trackParamX/F");
   reso->Branch("trackParamY", &trackParamY, "trackParamY/F");
   reso->Branch("trackParamDXDZ", &trackParamDXDZ, "trackParamDXDZ/F");
   reso->Branch("trackParamDYDZ", &trackParamDYDZ, "trackParamDYDZ/F");
   reso->Branch("trackParamXE", &trackParamXE, "trackParamXE/F");
   reso->Branch("trackParamYE", &trackParamYE, "trackParamYE/F");
   reso->Branch("trackParamDXDZE", &trackParamDXDZE, "trackParamDXDZE/F");
   reso->Branch("trackParamDYDZE", &trackParamDYDZE, "trackParamDYDZE/F");
   reso->Branch("pairsOnly", &pairsOnly, "pairsOnly/I");
   treso = fs->make<TTree>("treso", "tree tracks  for resolution studies");
   treso->Branch("track_momentum", &track_momentum, "track_momentum/F");
   treso->Branch("track_pt", &track_pt, "track_pt/F");
   treso->Branch("track_eta", &track_eta, "track_eta/F");
   treso->Branch("track_phi", &track_phi, "track_phi/F");
   treso->Branch("track_trackChi2", &track_trackChi2, "track_trackChi2/F");
   treso->Branch("track_width", &expWidth, "track_width/F");  // from 1D HIT
   treso->Branch("NumberOf_tracks", &NumberOf_tracks, "NumberOf_tracks/I");
 
   events = 0;
   EventTrackCKF = 0;
 
   histos2d_["track_phi_vs_eta"] = new TH2F("track_phi_vs_eta", ";track phi;track eta", 60, -3.5, 3.5, 60, -3., 3.);
   histos2d_["residual_vs_trackMomentum"] = new TH2F("residual_vs_trackMomentum",
                                                     ";track momentum [GeV]; x_{pred_track} - x_{reco_hit} [#mum]",
                                                     60,
                                                     0.,
                                                     10.,
                                                     60,
                                                     0.,
                                                     200.);
   histos2d_["residual_vs_trackPt"] = new TH2F(
       "residual_vs_trackPt", ";track p_{T}[GeV];x_{pred_track} - x_{reco_hit} [#mum]", 60, 0., 10., 60, 0., 200.);
   histos2d_["residual_vs_trackEta"] =
       new TH2F("residual_vs_trackEta", ";track #eta;x_{pred_track} - x_{reco_hit} [#mum]", 60, 0., 3., 60, 0., 200.);
   histos2d_["residual_vs_trackPhi"] =
       new TH2F("residual_vs_trackPhi", ";track #phi;x_{pred_track} - x_{reco_hit} [#mum]", 60, 0., 3.5, 60, 0., 200.);
   histos2d_["residual_vs_expectedWidth"] = new TH2F(
       "residual_vs_expectedWidth", ";track Width;x_{pred_track} - x_{reco_hit} [#mum]", 3, 0., 3., 60, 0., 200.);
   histos2d_["numHits_vs_residual"] =
       new TH2F("numHits_vs_residual", ";x_{pred_track} - x_{reco_hit} [#mum];N Hits", 60, 0., 200., 15, 0., 15.);
 }

◆ checkConsistency()

double HitResol::checkConsistency	(	const StripClusterParameterEstimator::LocalValues &	parameters,
		double	xx,
		double	xerr
	)

Definition at line 459 of file HitResol.cc.

References funct::abs(), relativeConstraints::error, mathSSE::sqrt(), and geometryCSVtoXML::xx.

                                                {
   double error = sqrt(parameters.second.xx() + xerr * xerr);
   double separation = abs(parameters.first.x() - xx);
   double consistency = separation / error;
   return consistency;
 }

◆ endJob()

void HitResol::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 450 of file HitResol.cc.

References events, EventTrackCKF, LogDebug, reso, and treso.

Referenced by o2olib.O2ORunMgr::executeJob().

                       {
   LogDebug("SiStripHitResolution:HitResol") << " Events Analysed             " << events << endl;
   LogDebug("SiStripHitResolution:HitResol") << " Number Of Tracked events    " << EventTrackCKF << endl;
 
   reso->GetDirectory()->cd();
   reso->Write();
   treso->Write();
 }

◆ fillDescriptions()

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

static

Definition at line 630 of file HitResol.cc.

References edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, and HLT_2022v15_cff::InputTag.

                                                                           {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("lumiScalers", edm::InputTag("scalersRawToDigi"));
   desc.add<edm::InputTag>("combinatorialTracks", edm::InputTag("generalTracks"));
   desc.add<edm::InputTag>("trajectories", edm::InputTag("generalTracks"));
   desc.addUntracked<int>("CompressionSettings", -1);
   desc.add<int>("Layer", 0);
   desc.add<bool>("Debug", false);
   desc.addUntracked<bool>("addLumi", false);
   desc.addUntracked<bool>("cutOnTracks", false);
   desc.addUntracked<unsigned int>("trackMultiplicity", 100);
   desc.addUntracked<double>("MomentumCut", 3.);
   desc.addUntracked<unsigned int>("UsePairsOnly", 1);
   descriptions.addWithDefaultLabel(desc);
 }

◆ getPairParameters()

bool HitResol::getPairParameters	(	const MagneticField *	magField_,
		AnalyticalPropagator &	propagator,
		const TrajectoryMeasurement *	traj1,
		const TrajectoryMeasurement *	traj2,
		float &	pairPath,
		float &	hitDX,
		float &	trackDX,
		float &	trackDXE,
		float &	trackParamX,
		float &	trackParamY,
		float &	trackParamDXDZ,
		float &	trackParamDYDZ,
		float &	trackParamXE,
		float &	trackParamYE,
		float &	trackParamDXDZE,
		float &	trackParamDYDZE
	)

Definition at line 511 of file HitResol.cc.

References TrajectoryMeasurement::backwardPredictedState(), TrajectoryStateCombiner::combine(), change_name::diff, TrajectoryMeasurement::forwardPredictedState(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalParameters(), TrajectoryStateOnSurface::globalPosition(), hitDX, mps_fire::i, TrajectoryStateOnSurface::isValid(), dqmiolumiharvest::j, JacobianCurvilinearToLocal::jacobian(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), LogDebug, PV3DBase< T, PVType, FrameType >::mag(), LocalTrajectoryError::matrix(), pairPath, TrackCandidateProducer_cfi::propagator, TrajectoryMeasurement::recHit(), mathSSE::sqrt(), TrajectoryStateOnSurface::surface(), trackDX, trackDXE, trackParamDXDZ, trackParamDXDZE, trackParamDYDZ, trackParamDYDZE, trackParamX, trackParamXE, trackParamY, trackParamYE, and LocalTrajectoryParameters::vector().

Referenced by analyze().

                                                          {
   pairPath = 0;
   hitDX = 0;
   trackDX = 0;
   trackDXE = 0;
 
   trackParamX = 0;
   trackParamY = 0;
   trackParamDXDZ = 0;
   trackParamDYDZ = 0;
   trackParamXE = 0;
   trackParamYE = 0;
   trackParamDXDZE = 0;
   trackParamDYDZE = 0;
 
   TrajectoryStateCombiner combiner_;
 
   // backward predicted state at module 1
   const TrajectoryStateOnSurface& bwdPred1 = traj1->backwardPredictedState();
   if (!bwdPred1.isValid())
     return false;
   LogDebug("HitResol") << "momentum from backward predicted state = " << bwdPred1.globalMomentum().mag() << endl;
   // forward predicted state at module 2
   const TrajectoryStateOnSurface& fwdPred2 = traj2->forwardPredictedState();
   LogDebug("HitResol") << "momentum from forward predicted state = " << fwdPred2.globalMomentum().mag() << endl;
   if (!fwdPred2.isValid())
     return false;
   // extrapolate fwdPred2 to module 1
   TrajectoryStateOnSurface fwdPred2At1 = propagator.propagate(fwdPred2, bwdPred1.surface());
   if (!fwdPred2At1.isValid())
     return false;
   // combine fwdPred2At1 with bwdPred1 (ref. state, best estimate without hits 1 and 2)
   TrajectoryStateOnSurface comb1 = combiner_.combine(bwdPred1, fwdPred2At1);
   if (!comb1.isValid())
     return false;
 
   //
   // propagation of reference parameters to module 2
   //
   std::pair<TrajectoryStateOnSurface, double> tsosWithS = propagator.propagateWithPath(comb1, fwdPred2.surface());
   TrajectoryStateOnSurface comb1At2 = tsosWithS.first;
   if (!comb1At2.isValid())
     return false;
   //distance of propagation from one surface to the next==could cut here
   pairPath = tsosWithS.second;
   if (TMath::Abs(pairPath) > 15)
     return false;  //cut to remove hit pairs > 15 cm apart
 
   // local parameters and errors on module 1
   AlgebraicVector5 pars = comb1.localParameters().vector();
   AlgebraicSymMatrix55 errs = comb1.localError().matrix();
   //number 3 is predX
   double predX1 = pars[3];
   //track fitted parameters in local coordinates for position 0
   (trackParamX) = pars[3];
   (trackParamY) = pars[4];
   (trackParamDXDZ) = pars[1];
   (trackParamDYDZ) = pars[2];
   (trackParamXE) = TMath::Sqrt(errs(3, 3));
   (trackParamYE) = TMath::Sqrt(errs(4, 4));
   (trackParamDXDZE) = TMath::Sqrt(errs(1, 1));
   (trackParamDYDZE) = TMath::Sqrt(errs(2, 2));
 
   // local parameters and errors on module 2
   pars = comb1At2.localParameters().vector();
   errs = comb1At2.localError().matrix();
   double predX2 = pars[3];
 
   JacobianLocalToCurvilinear jacLocToCurv(comb1.surface(), comb1.localParameters(), *magField_);
   AnalyticalCurvilinearJacobian jacCurvToCurv(
       comb1.globalParameters(), comb1At2.globalPosition(), comb1At2.globalMomentum(), tsosWithS.second);
   JacobianCurvilinearToLocal jacCurvToLoc(comb1At2.surface(), comb1At2.localParameters(), *magField_);
   // combined jacobian local-1-to-local-2
   AlgebraicMatrix55 jacobian = jacLocToCurv.jacobian() * jacCurvToCurv.jacobian() * jacCurvToLoc.jacobian();
   // covariance on module 1
   AlgebraicSymMatrix55 covComb1 = comb1.localError().matrix();
   // variance and correlations for predicted local_x on modules 1 and 2
   double c00 = covComb1(3, 3);
   double c10(0.);
   double c11(0.);
   for (int i = 1; i < 5; ++i) {
     c10 += jacobian(3, i) * covComb1(i, 3);
     for (int j = 1; j < 5; ++j)
       c11 += jacobian(3, i) * covComb1(i, j) * jacobian(3, j);
   }
   // choose relative sign in order to minimize error on difference
   double diff = c00 - 2 * fabs(c10) + c11;
   diff = diff > 0 ? sqrt(diff) : -sqrt(-diff);
   (trackDXE) = diff;
   double relativeXSign_ = c10 > 0 ? -1 : 1;
 
   (trackDX) = predX1 + relativeXSign_ * predX2;
 
   double recHitX_1 = traj1->recHit()->localPosition().x();
   double recHitX_2 = traj2->recHit()->localPosition().x();
 
   (hitDX) = recHitX_1 + relativeXSign_ * recHitX_2;
 
   return true;
 }

◆ getSimHitRes()

void HitResol::getSimHitRes	(	const GeomDetUnit *	det,
		const LocalVector &	trackdirection,
		const TrackingRecHit &	recHit,
		float &	trackWidth,
		float *	pitch,
		LocalVector &	drift
	)

Definition at line 468 of file HitResol.cc.

References shallow::drift(), position, rpcPointValidation_cfi::recHit, funct::tan(), thickness, trackWidth, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by analyze().

                                                 {
   const auto stripdet = dynamic_cast<const StripGeomDetUnit*>(det);
   const auto& topol = dynamic_cast<const StripTopology&>(stripdet->topology());
 
   LocalPoint position = recHit.localPosition();
   (*pitch) = topol.localPitch(position);
 
   float anglealpha = 0;
   if (trackdirection.z() != 0) {
     anglealpha = atan(trackdirection.x() / trackdirection.z()) * TMath::RadToDeg();
   }
 
   //  LocalVector drift = stripcpe.driftDirection(stripdet);
   float thickness = stripdet->surface().bounds().thickness();
   float tanalpha = tan(anglealpha * TMath::DegToRad());
   float tanalphaL = drift.x() / drift.z();
   (trackWidth) = fabs((thickness / (*pitch)) * tanalpha - (thickness / (*pitch)) * tanalphaL);
 }

◆ getSimpleRes()

double HitResol::getSimpleRes ( const TrajectoryMeasurement * traj1 )

Definition at line 492 of file HitResol.cc.

References TrajectoryMeasurement::backwardPredictedState(), TrajectoryStateCombiner::combine(), TrajectoryMeasurement::forwardPredictedState(), TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localPosition(), TrajectoryMeasurement::recHit(), and PV3DBase< T, PVType, FrameType >::x().

Referenced by analyze().

                                                                 {
   TrajectoryStateOnSurface theCombinedPredictedState;
 
   if (traj1->backwardPredictedState().isValid())
     theCombinedPredictedState =
         TrajectoryStateCombiner().combine(traj1->forwardPredictedState(), traj1->backwardPredictedState());
   else
     theCombinedPredictedState = traj1->forwardPredictedState();
 
   if (!theCombinedPredictedState.isValid()) {
     return -100;
   }
 
   const TransientTrackingRecHit::ConstRecHitPointer& firstRecHit = traj1->recHit();
   double recHitX_1 = firstRecHit->localPosition().x();
   return (theCombinedPredictedState.localPosition().x() - recHitX_1);
 }

Member Data Documentation

◆ addLumi_

const bool HitResol::addLumi_

private

Definition at line 85 of file HitResol.h.

◆ atEdge

float HitResol::atEdge

private

Definition at line 135 of file HitResol.h.

Referenced by analyze(), and beginJob().

◆ atEdge_2

float HitResol::atEdge_2

private

Definition at line 136 of file HitResol.h.

Referenced by analyze(), and beginJob().

◆ clusterCharge

unsigned int HitResol::clusterCharge

private

Definition at line 119 of file HitResol.h.

◆ clusterCharge_2

unsigned int HitResol::clusterCharge_2

private

Definition at line 120 of file HitResol.h.

◆ clusterWidth

unsigned int HitResol::clusterWidth

private

Definition at line 117 of file HitResol.h.

Referenced by analyze(), and beginJob().

◆ clusterWidth_2

unsigned int HitResol::clusterWidth_2

private

Definition at line 118 of file HitResol.h.

Referenced by analyze(), and beginJob().

◆ cmToUm

constexpr float HitResol::cmToUm = 10000.f

staticprivate

Definition at line 100 of file HitResol.h.

◆ combinatorialTracks_token_

const edm::EDGetTokenT<reco::TrackCollection> HitResol::combinatorialTracks_token_

private

Definition at line 73 of file HitResol.h.

Referenced by analyze().

◆ compSettings_

const int HitResol::compSettings_

private

Definition at line 89 of file HitResol.h.

Referenced by beginJob().

◆ cpeToken_

const edm::ESGetToken<StripClusterParameterEstimator, TkStripCPERecord> HitResol::cpeToken_

private

Definition at line 80 of file HitResol.h.

Referenced by analyze().

◆ cutOnTracks_

const bool HitResol::cutOnTracks_

private

Definition at line 87 of file HitResol.h.

◆ DEBUG_

const bool HitResol::DEBUG_

private

Definition at line 86 of file HitResol.h.

◆ driftAlpha

float HitResol::driftAlpha

private

Definition at line 129 of file HitResol.h.

◆ driftAlpha_2

float HitResol::driftAlpha_2

private

Definition at line 130 of file HitResol.h.

◆ events

int HitResol::events

private

Definition at line 103 of file HitResol.h.

Referenced by eventsfwlite.Events::__getattr__(), eventsfwlite.Events::__getitem__(), eventsfwlite.Events::__len__(), analyze(), beginJob(), endJob(), looper.Looper::loop(), and looper.Looper::process().

◆ EventTrackCKF

int HitResol::EventTrackCKF

private

Definition at line 103 of file HitResol.h.

Referenced by beginJob(), and endJob().

◆ expWidth

float HitResol::expWidth

private

Definition at line 127 of file HitResol.h.

Referenced by analyze(), and beginJob().

◆ expWidth_2

float HitResol::expWidth_2

private

Definition at line 128 of file HitResol.h.

Referenced by analyze(), and beginJob().

◆ geomToken_

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

private

Definition at line 79 of file HitResol.h.

Referenced by analyze().

◆ histos2d_

std::map<TString, TH2F*> HitResol::histos2d_

private

Definition at line 97 of file HitResol.h.

Referenced by analyze(), and beginJob().

◆ hitDX

float HitResol::hitDX

private

Definition at line 138 of file HitResol.h.

Referenced by analyze(), beginJob(), and getPairParameters().

◆ iidd1

unsigned int HitResol::iidd1

private

Definition at line 121 of file HitResol.h.

Referenced by analyze(), and beginJob().

◆ iidd2

unsigned int HitResol::iidd2

private

Definition at line 122 of file HitResol.h.

Referenced by analyze(), and beginJob().

◆ layers_

const unsigned int HitResol::layers_

private

Definition at line 91 of file HitResol.h.

◆ magFieldToken_

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

private

Definition at line 82 of file HitResol.h.

Referenced by analyze().

◆ momentumCut_

const double HitResol::momentumCut_

private

Definition at line 88 of file HitResol.h.

Referenced by analyze().

◆ mymom

float HitResol::mymom

private

Definition at line 107 of file HitResol.h.

Referenced by analyze(), and beginJob().

◆ mypitch1

float HitResol::mypitch1

private

Definition at line 125 of file HitResol.h.

Referenced by analyze(), and beginJob().

◆ mypitch2

float HitResol::mypitch2

private

Definition at line 126 of file HitResol.h.

◆ N1

float HitResol::N1

private

Definition at line 155 of file HitResol.h.

◆ N1uProj

float HitResol::N1uProj

private

Definition at line 157 of file HitResol.h.

◆ N2

float HitResol::N2

private

Definition at line 156 of file HitResol.h.

◆ N2uProj

float HitResol::N2uProj

private

Definition at line 158 of file HitResol.h.

◆ Nstrips

int HitResol::Nstrips

private

Definition at line 159 of file HitResol.h.

Referenced by analyze().

◆ Nstrips_2

int HitResol::Nstrips_2

private

Definition at line 160 of file HitResol.h.

Referenced by analyze().

◆ NumberOf_tracks

int HitResol::NumberOf_tracks

private

Definition at line 109 of file HitResol.h.

Referenced by beginJob().

◆ numHits

int HitResol::numHits

private

Definition at line 108 of file HitResol.h.

Referenced by analyze(), and beginJob().

◆ pairPath

float HitResol::pairPath

private

Definition at line 124 of file HitResol.h.

Referenced by analyze(), beginJob(), and getPairParameters().

◆ pairsOnly

unsigned int HitResol::pairsOnly

private

Definition at line 123 of file HitResol.h.

Referenced by analyze(), and beginJob().

◆ ProbTrackChi2

float HitResol::ProbTrackChi2

private

Definition at line 110 of file HitResol.h.

Referenced by analyze(), and beginJob().

◆ reso

TTree* HitResol::reso

private

Definition at line 95 of file HitResol.h.

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

◆ scalerToken_

const edm::EDGetTokenT<LumiScalersCollection> HitResol::scalerToken_

private

Definition at line 72 of file HitResol.h.

◆ simpleRes

float HitResol::simpleRes

private

Definition at line 137 of file HitResol.h.

Referenced by analyze(), and beginJob().

◆ siStripQualityToken_

const edm::ESGetToken<SiStripQuality, SiStripQualityRcd> HitResol::siStripQualityToken_

private

Definition at line 81 of file HitResol.h.

Referenced by analyze().

◆ StripCPE1_smp_pos_error

float HitResol::StripCPE1_smp_pos_error

private

Definition at line 111 of file HitResol.h.

◆ StripCPE2_smp_pos_error

float HitResol::StripCPE2_smp_pos_error

private

Definition at line 112 of file HitResol.h.

◆ StripErrorSquared1

float HitResol::StripErrorSquared1

private

Definition at line 113 of file HitResol.h.

◆ StripErrorSquared2

float HitResol::StripErrorSquared2

private

Definition at line 114 of file HitResol.h.

◆ thickness

float HitResol::thickness

private

Definition at line 131 of file HitResol.h.

Referenced by getSimHitRes().

◆ thickness_2

float HitResol::thickness_2

private

Definition at line 132 of file HitResol.h.

◆ tjToken_

const edm::EDGetTokenT<std::vector<Trajectory> > HitResol::tjToken_

private

Definition at line 74 of file HitResol.h.

Referenced by analyze().

◆ tkToken_

const edm::EDGetTokenT<reco::TrackCollection> HitResol::tkToken_

private

Definition at line 75 of file HitResol.h.

◆ topoToken_

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

private

Definition at line 78 of file HitResol.h.

Referenced by analyze().

◆ track_eta

float HitResol::track_eta

private

Definition at line 151 of file HitResol.h.

Referenced by analyze(), and beginJob().

◆ track_momentum

float HitResol::track_momentum

private

Definition at line 149 of file HitResol.h.

Referenced by analyze(), and beginJob().

◆ track_phi

float HitResol::track_phi

private

Definition at line 153 of file HitResol.h.

Referenced by analyze(), and beginJob().

◆ track_pt

float HitResol::track_pt

private

Definition at line 150 of file HitResol.h.

Referenced by analyze(), and beginJob().

◆ track_trackChi2

float HitResol::track_trackChi2

private

Definition at line 154 of file HitResol.h.

Referenced by analyze(), and beginJob().

◆ track_width

float HitResol::track_width

private

Definition at line 152 of file HitResol.h.

◆ trackDX

float HitResol::trackDX

private

Definition at line 139 of file HitResol.h.

Referenced by analyze(), beginJob(), and getPairParameters().

◆ trackDXE

float HitResol::trackDXE

private

Definition at line 140 of file HitResol.h.

Referenced by analyze(), beginJob(), and getPairParameters().

◆ trackMultiplicityCut_

const unsigned int HitResol::trackMultiplicityCut_

private

Definition at line 92 of file HitResol.h.

◆ trackParamDXDZ

float HitResol::trackParamDXDZ

private

Definition at line 145 of file HitResol.h.

Referenced by analyze(), beginJob(), and getPairParameters().

◆ trackParamDXDZE

float HitResol::trackParamDXDZE

private

Definition at line 147 of file HitResol.h.

Referenced by analyze(), beginJob(), and getPairParameters().

◆ trackParamDYDZ

float HitResol::trackParamDYDZ

private

Definition at line 146 of file HitResol.h.

Referenced by analyze(), beginJob(), and getPairParameters().

◆ trackParamDYDZE

float HitResol::trackParamDYDZE

private

Definition at line 148 of file HitResol.h.

Referenced by analyze(), beginJob(), and getPairParameters().

◆ trackParamX

float HitResol::trackParamX

private

Definition at line 141 of file HitResol.h.

Referenced by analyze(), beginJob(), and getPairParameters().

◆ trackParamXE

float HitResol::trackParamXE

private

Definition at line 143 of file HitResol.h.

Referenced by analyze(), beginJob(), and getPairParameters().

◆ trackParamY

float HitResol::trackParamY

private

Definition at line 142 of file HitResol.h.

Referenced by analyze(), beginJob(), and getPairParameters().

◆ trackParamYE

float HitResol::trackParamYE

private

Definition at line 144 of file HitResol.h.

Referenced by analyze(), beginJob(), and getPairParameters().

◆ trackWidth

float HitResol::trackWidth

private

Definition at line 133 of file HitResol.h.

Referenced by getSimHitRes().

◆ trackWidth_2

float HitResol::trackWidth_2

private

Definition at line 134 of file HitResol.h.

◆ treso

TTree* HitResol::treso

private

Definition at line 96 of file HitResol.h.

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

◆ uerr2

float HitResol::uerr2

private

Definition at line 115 of file HitResol.h.

◆ uerr2_2

float HitResol::uerr2_2

private

Definition at line 116 of file HitResol.h.

◆ usePairsOnly_

const unsigned int HitResol::usePairsOnly_

private

Definition at line 90 of file HitResol.h.

Referenced by analyze().

Public Types

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ TrajectoryCollection

Constructor & Destructor Documentation

◆ HitResol()

◆ ~HitResol()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ checkConsistency()

◆ endJob()

◆ fillDescriptions()

◆ getPairParameters()

◆ getSimHitRes()

◆ getSimpleRes()

Member Data Documentation

◆ addLumi_

◆ atEdge

◆ atEdge_2

◆ clusterCharge

◆ clusterCharge_2

◆ clusterWidth

◆ clusterWidth_2

◆ cmToUm

◆ combinatorialTracks_token_

◆ compSettings_

◆ cpeToken_

◆ cutOnTracks_

◆ DEBUG_

◆ driftAlpha

◆ driftAlpha_2

◆ events

◆ EventTrackCKF

◆ expWidth

◆ expWidth_2

◆ geomToken_

◆ histos2d_

◆ hitDX

◆ iidd1

◆ iidd2

◆ layers_

◆ magFieldToken_

◆ momentumCut_

◆ mymom

◆ mypitch1

◆ mypitch2

◆ N1

◆ N1uProj

◆ N2

◆ N2uProj

◆ Nstrips

◆ Nstrips_2

◆ NumberOf_tracks

◆ numHits

◆ pairPath

◆ pairsOnly

◆ ProbTrackChi2

◆ reso

◆ scalerToken_

◆ simpleRes

◆ siStripQualityToken_

◆ StripCPE1_smp_pos_error

◆ StripCPE2_smp_pos_error

◆ StripErrorSquared1

◆ StripErrorSquared2

◆ thickness

◆ thickness_2

◆ tjToken_

◆ tkToken_

◆ topoToken_

◆ track_eta

◆ track_momentum

◆ track_phi