MuonSeedCreator.h

Go to the documentation of this file.

#ifndef RecoMuon_SeedGenerator_MuonSeedCreator_H
#define RecoMuon_SeedGenerator_MuonSeedCreator_H
 
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
#include "DataFormats/TrajectorySeed/interface/TrajectorySeed.h"
#include "RecoMuon/TransientTrackingRecHit/interface/MuonTransientTrackingRecHit.h"
 
#include "FWCore/ParameterSet/interface/ParameterSet.h"
 
#include <vector>
 
class RecHit;
class Plane;
class GeomDet;
class MagneticField;
class MuonTransientTrackingRecHit;
 
class MuonSeedCreator {
public:
  typedef MuonTransientTrackingRecHit::MuonRecHitContainer SegmentContainer;
 
  explicit MuonSeedCreator(const edm::ParameterSet& pset);
 
  ~MuonSeedCreator();
 
  // Operations
 
  void setBField(const MagneticField* theField) { BField = theField; };
 
  TrajectorySeed createSeed(
      int type, const SegmentContainer& seg, const std::vector<int>& layers, int NShower, int NShowerSeg);
 
private:
  void estimatePtCSC(const SegmentContainer& seg, const std::vector<int>& layers, double& pt, double& spt);
 
  void estimatePtOverlap(const SegmentContainer& seg, const std::vector<int>& layers, double& pt, double& spt);
 
  void estimatePtDT(const SegmentContainer& seg, const std::vector<int>& layers, double& pt, double& spt);
 
  void estimatePtSingle(const SegmentContainer& seg, const std::vector<int>& layers, double& pt, double& spt);
 
  void estimatePtShowering(int& NShowers, int& NShowerSeg, double& pt, double& spt);
 
  void weightedPt(const std::vector<double>& ptEstimate,
                  const std::vector<double>& sptEstimate,
                  double& ptAvg,
                  double& sptAvg);
 
  std::vector<double> getPt(const std::vector<double>& vParameters, double eta, double dPhi);
 
  double scaledPhi(double dphi, double t1);
 
  // Miminum and maximum pt momentum of a track
  float theMinMomentum;
  float theMaxMomentum;
  float defaultMomentum;
 
  // Error on pt estimate which prevents weighted average from blowing up ( spt --> 0 )
  double sysError;
 
  // Flag for internal debugging
  bool debug;
 
  // Cache Magnetic Field for current event
  const MagneticField* BField;
 
  // seed parameters vectors
  std::vector<double> DT12;
  std::vector<double> DT13;
  std::vector<double> DT14;
  std::vector<double> DT23;
  std::vector<double> DT24;
  std::vector<double> DT34;
 
  std::vector<double> CSC01;
  std::vector<double> CSC12;
  std::vector<double> CSC02;
  std::vector<double> CSC13;
  std::vector<double> CSC03;
  std::vector<double> CSC14;
  std::vector<double> CSC23;
  std::vector<double> CSC24;
  std::vector<double> CSC34;
 
  std::vector<double> OL1213;
  std::vector<double> OL1222;
  std::vector<double> OL1232;
  std::vector<double> OL2213;
  std::vector<double> OL2222;
 
  std::vector<double> SME11;
  std::vector<double> SME12;
  std::vector<double> SME13;
  std::vector<double> SME21;
  std::vector<double> SME22;
  std::vector<double> SME31;
  std::vector<double> SME32;
  std::vector<double> SME41;
 
  std::vector<double> SMB10;
  std::vector<double> SMB11;
  std::vector<double> SMB12;
  std::vector<double> SMB20;
  std::vector<double> SMB21;
  std::vector<double> SMB22;
  std::vector<double> SMB30;
  std::vector<double> SMB31;
  std::vector<double> SMB32;
 
  // dphi scaling factors
 
  std::vector<double> CSC01_1;
  std::vector<double> CSC12_1;
  std::vector<double> CSC12_2;
  std::vector<double> CSC12_3;
  std::vector<double> CSC13_2;
  std::vector<double> CSC13_3;
  std::vector<double> CSC14_3;
  std::vector<double> CSC23_1;
  std::vector<double> CSC23_2;
  std::vector<double> CSC24_1;
  std::vector<double> CSC34_1;
 
  std::vector<double> DT12_1;
  std::vector<double> DT12_2;
  std::vector<double> DT13_1;
  std::vector<double> DT13_2;
  std::vector<double> DT14_1;
  std::vector<double> DT14_2;
  std::vector<double> DT23_1;
  std::vector<double> DT23_2;
  std::vector<double> DT24_1;
  std::vector<double> DT24_2;
  std::vector<double> DT34_1;
  std::vector<double> DT34_2;
 
  std::vector<double> OL_1213;
  std::vector<double> OL_1222;
  std::vector<double> OL_1232;
  std::vector<double> OL_2213;
  std::vector<double> OL_2222;
 
  std::vector<double> SMB_10S;
  std::vector<double> SMB_11S;
  std::vector<double> SMB_12S;
  std::vector<double> SMB_20S;
  std::vector<double> SMB_21S;
  std::vector<double> SMB_22S;
  std::vector<double> SMB_30S;
  std::vector<double> SMB_31S;
  std::vector<double> SMB_32S;
 
  std::vector<double> SME_11S;
  std::vector<double> SME_12S;
  std::vector<double> SME_13S;
  std::vector<double> SME_21S;
  std::vector<double> SME_22S;
};
#endif