---

MuonStandaloneAlgorithm.h

Go to the documentation of this file.

#ifndef Alignment_MuonStandaloneAlgorithm_MuonStandaloneAlgorithm_H
#define Alignment_MuonStandaloneAlgorithm_MuonStandaloneAlgorithm_H


// Base Class Headers
#include "FWCore/Framework/interface/EDProducer.h"

#include "DataFormats/DTRecHit/interface/DTRecSegment4DCollection.h"
#include "DataFormats/CSCRecHit/interface/CSCSegmentCollection.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "Alignment/MuonStandaloneAlgorithm/interface/TrackForAlignment.h"
#include "RecoMuon/TrackingTools/interface/MuonServiceProxy.h"
#include "RecoMuon/TrackingTools/interface/MuonUpdatorAtVertex.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "TMatrixD.h"
#include "DataFormats/CLHEP/interface/AlgebraicObjects.h"


namespace edm {
  class ParameterSet;
  class Event;
  class EventSetup;
  class InputTag;
}


using namespace edm;


class Propagator;
class TrajectoryStateOnSurface;
class FreeTrajectoryState;

typedef std::vector< std::vector<int> > intDVector;
typedef std::vector<TrackingRecHit *> RecHitVector;

class MuonStandaloneAlgorithm: public edm::EDProducer {
public:
  MuonStandaloneAlgorithm(const edm::ParameterSet& pset);

  virtual ~MuonStandaloneAlgorithm();

  // Operations

  virtual void produce(edm::Event & event, const edm::EventSetup& eventSetup);

protected:

private:
  
  bool Derivative(TrajectoryStateOnSurface , const GeomDet *, int, std::vector<double>, TMatrixD *);

  bool calculateJacobian(TrajectoryStateOnSurface, const GeomDet *, std::vector<double>, TMatrixD *);

  void calculateError(AlgebraicSymMatrix55, TMatrixD *);
  
  RecHitVector doMatching(const reco::Track &, edm::Handle<DTRecSegment4DCollection> &, edm::Handle<CSCSegmentCollection> &, intDVector *, intDVector *, edm::ESHandle<GlobalTrackingGeometry> &); 

  FreeTrajectoryState produceFreeTrajectoryState(TrajectoryStateOnSurface, int, double, std::vector<double> );
  
  edm::InputTag theSTAMuonTag;
  
  std::string theSeedCollectionLabel;
  std::string theRecHits4DLabelDT;
  std::string theRecHits2DLabelCSC;
  int theSelectorOfFirstPoint;
  int verbosity;
  double curvStep;
  double thetaStep;
  double phiStep;
  double transStep;
  double longiStep;

  Propagator *thePropagator1;

  const static int NDOFTrack = 5;
  const static int NDOFAlign = 6;
  const static int NDOFChamber = 4;
  const static int NDOFCoor = 4;

      
};
#endif

---