NuclearTester.h

Go to the documentation of this file.

//----------------------------------------------------------------------------
//-----------------------------------------------------------------------------
#ifndef CD_NuclearTester_H_
#define CD_NuclearTester_H_
#include "TrackingTools/PatternTools/interface/TrajectoryMeasurement.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
#include "TrackingTools/DetLayers/interface/MeasurementEstimator.h"
 
class NuclearTester {
private:
  typedef TrajectoryMeasurement TM;
  typedef std::vector<TM> TMContainer;
  typedef TrajectoryMeasurement::ConstRecHitPointer ConstRecHitPointer;
  typedef std::pair<TrajectoryMeasurement, TMContainer> TMPair;
  typedef std::vector<TMPair> TMPairVector;
 
public:
  NuclearTester(unsigned int max_hits, const MeasurementEstimator* est, const TrackerGeometry* track_geom);
 
  bool isNuclearInteraction();
 
  double meanHitDistance(const std::vector<TrajectoryMeasurement>& vecTM) const;
 
  std::vector<GlobalPoint> HitPositions(const std::vector<TrajectoryMeasurement>& vecTM) const;
 
  double fwdEstimate(const std::vector<TrajectoryMeasurement>& vecTM) const;
 
  std::vector<TM>::const_iterator lastValidTM(const std::vector<TM>& vecTM) const;
 
  void push_back(const TM& init_tm, const TMContainer& vecTM) {
    allTM.push_back(std::make_pair(init_tm, vecTM));
    compatible_hits.push_back(vecTM.size());
  }
 
  const TMContainer& back() const { return allTM.back().second; }
 
  double meanHitDistance() const { return meanHitDistance(back()); }
 
  double fwdEstimate() const { return fwdEstimate(back()); }
 
  void reset(unsigned int nMeasurements) {
    allTM.clear();
    compatible_hits.clear();
    maxHits = (nMeasurements < maxHits) ? nMeasurements : maxHits;
  }
 
  int nuclearIndex() const { return NuclearIndex; }
 
  const TMPair& goodTMPair() const { return *(allTM.begin() + nuclearIndex() - 1); }
 
  unsigned int nHitsChecked() const { return compatible_hits.size(); }
 
  std::vector<int> compatibleHits() const { return compatible_hits; }
 
private:
  // data members
  TMPairVector allTM;
  std::vector<int> compatible_hits;
  int NuclearIndex;
 
  // input parameters
  unsigned int maxHits;
  const MeasurementEstimator* theEstimator;
  const TrackerGeometry* trackerGeom;
 
  bool checkWithMultiplicity();
};
#endif