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