AlignmentMonitorGeneric Class Reference

#include <AlignmentMonitorGeneric.h>

Inheritance diagram for AlignmentMonitorGeneric:

Public Member Functions
	AlignmentMonitorGeneric (const edm::ParameterSet &)
void	book () override
	Book or retrieve histograms; MUST be reimplemented. More...
void	event (const edm::Event &, const edm::EventSetup &, const ConstTrajTrackPairCollection &) override
	Called for each event (by "run()"): may be reimplemented. More...
Public Member Functions inherited from AlignmentMonitorBase
virtual void	afterAlignment ()
	AlignmentMonitorBase (const edm::ParameterSet &cfg, std::string name)
	Constructor. More...
void	beginOfJob (AlignableTracker *pTracker, AlignableMuon *pMuon, AlignmentParameterStore *pStore)
	Called at beginning of job: don't reimplement. More...
void	duringLoop (const edm::Event &iEvent, const edm::EventSetup &iSetup, const ConstTrajTrackPairCollection &iTrajTracks)
	Called for each event: don't reimplement. More...
void	endOfJob ()
	Called at end of processing: don't implement. More...
void	endOfLoop ()
	Called at end of loop: don't reimplement. More...
void	startingNewLoop ()
	Called at beginning of loop: don't reimplement. More...
virtual	~AlignmentMonitorBase ()
	Destructor. More...

Private Types
typedef std::vector< TH1F * >	Hist1Ds

Private Attributes
std::map< const Alignable *, Hist1Ds >	m_resHists
Hist1Ds	m_trkHists

Static Private Attributes
static const unsigned int	nBin_ = 50

Additional Inherited Members
Public Types inherited from AlignmentMonitorBase
typedef std::pair< const Trajectory *, const reco::Track * >	ConstTrajTrackPair
typedef std::vector< ConstTrajTrackPair >	ConstTrajTrackPairCollection
Protected Member Functions inherited from AlignmentMonitorBase
TH1F *	book1D (std::string dir, std::string name, std::string title, int nchX, double lowX, double highX)
TH2F *	book2D (std::string dir, std::string name, std::string title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
TProfile *	bookProfile (std::string dir, std::string name, std::string title, int nchX, double lowX, double highX, int nchY=1, double lowY=0., double highY=0., const char *option="s")
TFileDirectory *	directory (std::string dir)
int	iteration ()
AlignableMuon *	pMuon ()
AlignableNavigator *	pNavigator ()
AlignmentParameterStore *	pStore ()
AlignableTracker *	pTracker ()
Protected Attributes inherited from AlignmentMonitorBase
const edm::InputTag	m_beamSpotTag

Detailed Description

Definition at line 28 of file AlignmentMonitorGeneric.h.

Member Typedef Documentation

typedef std::vector<TH1F*> AlignmentMonitorGeneric::Hist1Ds

private

Definition at line 29 of file AlignmentMonitorGeneric.h.

Constructor & Destructor Documentation

AlignmentMonitorGeneric::AlignmentMonitorGeneric

(

const edm::ParameterSet &

cfg

)

Definition at line 11 of file AlignmentMonitorGeneric.cc.

    : AlignmentMonitorBase(cfg, "AlignmentMonitorGeneric") {}

Member Function Documentation

void AlignmentMonitorGeneric::book ( )

overridevirtual

Book or retrieve histograms; MUST be reimplemented.

Implements AlignmentMonitorBase.

Definition at line 14 of file AlignmentMonitorGeneric.cc.

References Alignable::alignableObjectId(), AlignmentParameterStore::alignables(), AlignmentMonitorBase::book1D(), AlignableObjectId::commonObjectIdProvider(), compare::hists, B2GTnPMonitor_cfi::histTitle, mps_fire::i, triggerObjects_cff::id, Alignable::id(), M_PI, m_resHists, m_trkHists, dqmiodumpmetadata::n, Skims_PA_cff::name, nBin_, AlignmentMonitorBase::pMuon(), AlignmentMonitorBase::pStore(), AlignmentMonitorBase::pTracker(), AlCaHLTBitMon_QueryRunRegistry::string, and DetId::subdetId().

                                   {
  std::vector<std::string> residNames;  // names of residual histograms

  residNames.push_back("x hit residuals pos track");
  residNames.push_back("x hit residuals neg track");
  residNames.push_back("y hit residuals pos track");
  residNames.push_back("y hit residuals neg track");

  auto alignableObjectId = AlignableObjectId::commonObjectIdProvider(pTracker(), pMuon());

  const auto& alignables = pStore()->alignables();

  unsigned int nAlignable = alignables.size();
  unsigned int nResidName = residNames.size();

  for (unsigned int i = 0; i < nAlignable; ++i) {
    const Alignable* ali = alignables[i];

    Hist1Ds& hists = m_resHists[ali];

    hists.resize(nResidName, nullptr);

    align::ID id = ali->id();
    align::StructureType type = ali->alignableObjectId();

    for (unsigned int n = 0; n < nResidName; ++n) {
      const std::string& name = residNames[n];

      TString histName(name.c_str());
      histName += Form("_%s_%d", alignableObjectId.idToString(type), id);
      histName.ReplaceAll(" ", "");

      TString histTitle(name.c_str());
      histTitle += Form(" for %s with ID %d (subdet %d)", alignableObjectId.idToString(type), id, DetId(id).subdetId());

      hists[n] = book1D(std::string("/iterN/") + std::string(name) + std::string("/"),
                        std::string(histName.Data()),
                        std::string(histTitle.Data()),
                        nBin_,
                        -5.,
                        5.);
    }
  }

  m_trkHists.resize(6, nullptr);

  m_trkHists[0] = book1D("/iterN/", "pt", "track p_{t} (GeV)", nBin_, 0.0, 100.0);
  m_trkHists[1] = book1D("/iterN/", "eta", "track #eta", nBin_, -3.0, 3.0);
  m_trkHists[2] = book1D("/iterN/", "phi", "track #phi", nBin_, -M_PI, M_PI);
  m_trkHists[3] = book1D("/iterN/", "d0", "track d0 (cm)", nBin_, -0.02, 0.02);
  m_trkHists[4] = book1D("/iterN/", "dz", "track dz (cm)", nBin_, -20.0, 20.0);
  m_trkHists[5] = book1D("/iterN/", "chi2", "track #chi^{2}/dof", nBin_, 0.0, 20.0);
}

void AlignmentMonitorGeneric::event ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup, const ConstTrajTrackPairCollection &  iTrajTracks  )

overridevirtual

Called for each event (by "run()"): may be reimplemented.

Reimplemented from AlignmentMonitorBase.

Definition at line 68 of file AlignmentMonitorGeneric.cc.

References AlignableNavigator::alignableFromDetId(), TrajectoryMeasurement::backwardPredictedState(), ALCARECOTkAlJpsiMuMu_cff::charge, reco::TrackBase::d0(), DEFINE_EDM_PLUGIN, TrackingRecHit::det(), reco::TrackBase::dz(), fftjetvertexadder_cfi::errX, fftjetvertexadder_cfi::errY, reco::TrackBase::eta(), TrajectoryMeasurement::forwardPredictedState(), TrackingRecHit::geographicalId(), h, TrackingRecHit::isValid(), TrackerGeomDet::localAlignmentError(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localPosition(), TrackingRecHit::localPosition(), TrackingRecHit::localPositionError(), visualization-live-secondInstance_cfg::m, m_resHists, m_trkHists, Alignable::mother(), reco::TrackBase::normalizedChi2(), reco::TrackBase::phi(), AlignmentMonitorBase::pNavigator(), LocalTrajectoryError::positionError(), reco::TrackBase::pt(), TrajectoryMeasurement::recHit(), mathSSE::sqrt(), OrderedSet::t, HLT_2018_cff::track, LocalError::xx(), LocalError::xy(), and LocalError::yy().

Referenced by Types.EventID::cppID().

                                                                                {
  TrajectoryStateCombiner tsoscomb;

  for (unsigned int t = 0; t < tracks.size(); ++t) {
    const reco::Track* track = tracks[t].second;

    float charge = tracks[t].second->charge();

    const std::vector<TrajectoryMeasurement>& meass = tracks[t].first->measurements();

    for (unsigned int m = 0; m < meass.size(); ++m) {
      const TrajectoryMeasurement& meas = meass[m];
      const TransientTrackingRecHit& hit = *meas.recHit();

      if (hit.isValid()) {
        const Alignable* ali = pNavigator()->alignableFromDetId(hit.geographicalId());

        while (ali) {
          std::map<const Alignable*, Hist1Ds>::iterator h = m_resHists.find(ali);
          if (h != m_resHists.end()) {
            TrajectoryStateOnSurface tsos = tsoscomb(meas.forwardPredictedState(), meas.backwardPredictedState());

            align::LocalVector res = tsos.localPosition() - hit.localPosition();
            LocalError err1 = tsos.localError().positionError();
            LocalError err2 = hit.localPositionError();  // CPE+APE

            // subtract APEs from err2 (if existing) from covariance matrix
            auto det = static_cast<const TrackerGeomDet*>(hit.det());
            const auto localAPE = det->localAlignmentError();
            if (localAPE.valid()) {
              err2 = LocalError(err2.xx() - localAPE.xx(), err2.xy() - localAPE.xy(), err2.yy() - localAPE.yy());
            }

            float errX = std::sqrt(err1.xx() + err2.xx());
            float errY = std::sqrt(err1.yy() + err2.yy());

            h->second[charge > 0 ? 0 : 1]->Fill(res.x() / errX);
            h->second[charge > 0 ? 2 : 3]->Fill(res.y() / errY);
          }
          ali = ali->mother();
        }
      }
    }

    m_trkHists[0]->Fill(track->pt());
    m_trkHists[1]->Fill(track->eta());
    m_trkHists[2]->Fill(track->phi());
    m_trkHists[3]->Fill(track->d0());
    m_trkHists[4]->Fill(track->dz());
    m_trkHists[5]->Fill(track->normalizedChi2());
  }
}

Member Data Documentation

std::map<const Alignable*, Hist1Ds> AlignmentMonitorGeneric::m_resHists

private

Definition at line 43 of file AlignmentMonitorGeneric.h.

Referenced by book(), and event().

Hist1Ds AlignmentMonitorGeneric::m_trkHists

private

Definition at line 41 of file AlignmentMonitorGeneric.h.

Referenced by book(), and event().

const unsigned int AlignmentMonitorGeneric::nBin_ = 50

staticprivate

Definition at line 39 of file AlignmentMonitorGeneric.h.

Referenced by book().