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