CMS 3D CMS Logo

All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
List of all members | Public Member Functions | Private Attributes
AlignmentMonitorTemplate Class Reference
Inheritance diagram for AlignmentMonitorTemplate:
AlignmentMonitorBase

Public Member Functions

void afterAlignment () override
 
 AlignmentMonitorTemplate (const edm::ParameterSet &cfg, edm::ConsumesCollector iC)
 
void book () override
 Book or retrieve histograms; MUST be reimplemented. More...
 
void event (const edm::Event &iEvent, const edm::EventSetup &iSetup, const ConstTrajTrackPairCollection &iTrajTracks) override
 Called for each event (by "run()"): may be reimplemented. More...
 
 ~AlignmentMonitorTemplate () override
 
- Public Member Functions inherited from AlignmentMonitorBase
 AlignmentMonitorBase (const edm::ParameterSet &cfg, const edm::ConsumesCollector &iC, std::string name)
 Constructor. More...
 
 AlignmentMonitorBase (const AlignmentMonitorBase &)=delete
 
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...
 
const AlignmentMonitorBaseoperator= (const AlignmentMonitorBase &)=delete
 
void startingNewLoop ()
 Called at beginning of loop: don't reimplement. More...
 
virtual ~AlignmentMonitorBase ()
 Destructor. More...
 

Private Attributes

TH1F * m_hist
 
TH1F * m_ihist
 
TH1F * m_otherdir
 
TH1F * m_otherdir2
 
std::map< Alignable *, TH1F * > m_residuals
 

Additional Inherited Members

- Public Types inherited from AlignmentMonitorBase
typedef std::pair< const Trajectory *, const reco::Track * > ConstTrajTrackPair
 
typedef std::vector< ConstTrajTrackPairConstTrajTrackPairCollection
 
- 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")
 
TFileDirectorydirectory (std::string dir)
 
int iteration ()
 
AlignableMuonpMuon ()
 
AlignableNavigatorpNavigator ()
 
AlignmentParameterStorepStore ()
 
AlignableTrackerpTracker ()
 
- Protected Attributes inherited from AlignmentMonitorBase
const edm::InputTag m_beamSpotTag
 

Detailed Description

Definition at line 28 of file AlignmentMonitorTemplate.cc.

Constructor & Destructor Documentation

◆ AlignmentMonitorTemplate()

AlignmentMonitorTemplate::AlignmentMonitorTemplate ( const edm::ParameterSet cfg,
edm::ConsumesCollector  iC 
)
inline

Definition at line 30 of file AlignmentMonitorTemplate.cc.

31  : AlignmentMonitorBase(cfg, iC, "AlignmentMonitorTemplate"){};
AlignmentMonitorBase(const edm::ParameterSet &cfg, const edm::ConsumesCollector &iC, std::string name)
Constructor.

◆ ~AlignmentMonitorTemplate()

AlignmentMonitorTemplate::~AlignmentMonitorTemplate ( )
inlineoverride

Definition at line 32 of file AlignmentMonitorTemplate.cc.

32 {};

Member Function Documentation

◆ afterAlignment()

void AlignmentMonitorTemplate::afterAlignment ( )
overridevirtual

Called after updating AlignableTracker and AlignableMuon (by "endOfLoop()"): may be reimplemented

Reimplemented from AlignmentMonitorBase.

Definition at line 130 of file AlignmentMonitorTemplate.cc.

References AlignmentMonitorBase::iteration(), and m_otherdir.

130  {
131  m_otherdir->Fill(
132  iteration()); // this one will only get one fill per iteration, because it's called in afterAlignment()
133 }

◆ book()

void AlignmentMonitorTemplate::book ( )
overridevirtual

Book or retrieve histograms; MUST be reimplemented.

Implements AlignmentMonitorBase.

Definition at line 57 of file AlignmentMonitorTemplate.cc.

References AlignmentParameterStore::alignables(), AlignmentMonitorBase::book1D(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, m_hist, m_ihist, m_otherdir, m_otherdir2, m_residuals, Skims_PA_cff::name, AlignmentMonitorBase::pStore(), and runGCPTkAlMap::title.

57  {
58  m_hist = book1D("/", "hist", "hist", 10, 0.5, 10.5); // there's only one of these per job
59  m_ihist = book1D("/iterN/", "ihist", "ihist", 10, 0, 1); // there's a new one of these for each iteration
60  // "/iterN/" is a special directory name, in which the "N" gets replaced by the current iteration number.
61 
62  m_otherdir = book1D("/otherdir/", "hist", "this is a histogram in another directory", 10, 0.5, 10.5);
63  m_otherdir2 =
64  book1D("/iterN/otherdir/", "hist", "here's one in another directory inside the iterN directories", 10, 0.5, 10.5);
65 
66  // This is a procedure that makes one histogram for each selected alignable, and puts them in the iterN directory.
67  // This is not a constant-time lookup. If you need something faster, see AlignmentMonitorMuonHIP, which has a
68  // dynamically-allocated array of TH1F*s.
69  const auto& alignables = pStore()->alignables();
70  for (const auto& it : alignables) {
71  char name[256], title[256];
72  snprintf(name, sizeof(name), "xresid%d", it->geomDetId().rawId());
73  snprintf(title, sizeof(title), "x track-hit for DetId %d", it->geomDetId().rawId());
74 
75  m_residuals[it] = book1D("/iterN/", name, title, 100, -5., 5.);
76  }
77 
78  // Important: you create TObject pointers with the "new" operator, but YOU don't delete them. They're deleted by the
79  // base class, which knows when they are no longer needed (based on whether they are in the /iterN/ directory or not).
80 
81  // For more detail, see the twiki page:
82  // https://twiki.cern.ch/twiki/bin/view/CMS/SWGuideAlignmentMonitors for creating a plugin like this one
83  // https://twiki.cern.ch/twiki/bin/view/CMS/SWGuideAlignmentAlgorithms#Monitoring for configuring it
84 }
AlignmentParameterStore * pStore()
std::map< Alignable *, TH1F * > m_residuals
const align::Alignables & alignables(void) const
get all alignables
TH1F * book1D(std::string dir, std::string name, std::string title, int nchX, double lowX, double highX)

◆ event()

void AlignmentMonitorTemplate::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 86 of file AlignmentMonitorTemplate.cc.

References AlignableNavigator::alignableFromDetId(), TrajectoryMeasurement::backwardPredictedState(), TrajectoryStateCombiner::combine(), AlignableNavigator::detAndSubdetInMap(), TrajectoryMeasurement::forwardPredictedState(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, AlignmentMonitorBase::iteration(), TrajectoryStateOnSurface::localPosition(), m_hist, m_ihist, m_residuals, Trajectory::measurements(), Alignable::mother(), AlignmentMonitorBase::pNavigator(), TrajectoryMeasurement::recHit(), cond::persistency::search(), DiMuonV_cfg::tracks, PV3DBase< T, PVType, FrameType >::x(), and hit::x.

Referenced by Types.EventID::cppID().

88  {
89  m_hist->Fill(iteration()); // get the number of events per iteration
90  m_ihist->Fill(0.5); // get the number of events in this iteration in the central bin
91 
92  TrajectoryStateCombiner tsoscomb;
93 
94  // This is a procedure that loops over tracks/hits, calculates residuals, and fills the appropriate histogram.
95  for (ConstTrajTrackPairCollection::const_iterator it = tracks.begin(); it != tracks.end(); ++it) {
96  const Trajectory* traj = it->first;
97  // const reco::Track* track = it->second;
98  // If your tracks are refit using the producer in RecoTracker, you'll get updated reco::Track objects with
99  // each iteration, and it makes sense to make plots using these.
100  // If your tracks are refit using TrackingTools/TrackRefitter, only the Trajectories will be updated.
101  // We're working on that. I'll try to remember to change this comment when the update is ready.
102 
103  std::vector<TrajectoryMeasurement> measurements = traj->measurements();
104  for (std::vector<TrajectoryMeasurement>::const_iterator im = measurements.begin(); im != measurements.end(); ++im) {
105  const TrajectoryMeasurement& meas = *im;
106  const TransientTrackingRecHit* hit = &(*meas.recHit());
107  const DetId id = hit->geographicalId();
108 
109  if (hit->isValid() && pNavigator()->detAndSubdetInMap(id)) {
110  // Combine the forward-propagated state with the backward-propagated state to get a minimally-biased residual.
111  // This is the same procedure that is used to calculate residuals in the HIP algorithm
113 
114  // Search for our histogram using the Alignable* -> TH1F* map
115  // The "alignable = alignable->mother()" part ascends the alignable tree, because hits are on the lowest-level
116  // while our histograms may be associated with higher-level Alignables.
117  Alignable* alignable = pNavigator()->alignableFromDetId(id);
118  std::map<Alignable*, TH1F*>::const_iterator search = m_residuals.find(alignable);
119  while (search == m_residuals.end() && (alignable = alignable->mother()))
120  search = m_residuals.find(alignable);
121 
122  if (search != m_residuals.end()) {
123  search->second->Fill(tsosc.localPosition().x() - hit->localPosition().x());
124  }
125  } // end if hit is valid
126  } // end loop over hits
127  } // end loop over tracks/trajectories
128 }
AlignableNavigator * pNavigator()
Alignable * mother() const
Return pointer to container alignable (if any)
Definition: Alignable.h:91
std::vector< T >::const_iterator search(const cond::Time_t &val, const std::vector< T > &container)
Definition: IOVProxy.cc:21
std::map< Alignable *, TH1F * > m_residuals
TrajectoryStateOnSurface const & forwardPredictedState() const
Access to forward predicted state (from fitter or builder)
bool detAndSubdetInMap(const DetId &detid) const
Given a DetId, returns true if DetIds with this detector and subdetector id are in the map (not neces...
DataContainer const & measurements() const
Definition: Trajectory.h:178
T x() const
Definition: PV3DBase.h:59
TrajectoryStateOnSurface const & backwardPredictedState() const
Access to backward predicted state (from smoother)
Definition: DetId.h:17
AlignableDetOrUnitPtr alignableFromDetId(const DetId &detid)
Returns AlignableDetOrUnitPtr corresponding to given DetId.
TSOS combine(const TSOS &pTsos1, const TSOS &pTsos2) const
ConstRecHitPointer const & recHit() const

Member Data Documentation

◆ m_hist

TH1F* AlignmentMonitorTemplate::m_hist
private

Definition at line 41 of file AlignmentMonitorTemplate.cc.

Referenced by book(), and event().

◆ m_ihist

TH1F * AlignmentMonitorTemplate::m_ihist
private

Definition at line 41 of file AlignmentMonitorTemplate.cc.

Referenced by book(), and event().

◆ m_otherdir

TH1F * AlignmentMonitorTemplate::m_otherdir
private

Definition at line 41 of file AlignmentMonitorTemplate.cc.

Referenced by afterAlignment(), and book().

◆ m_otherdir2

TH1F * AlignmentMonitorTemplate::m_otherdir2
private

Definition at line 41 of file AlignmentMonitorTemplate.cc.

Referenced by book().

◆ m_residuals

std::map<Alignable*, TH1F*> AlignmentMonitorTemplate::m_residuals
private

Definition at line 42 of file AlignmentMonitorTemplate.cc.

Referenced by book(), and event().