SeedForPhotonConversion1Leg.h

Go to the documentation of this file.

#ifndef SeedForPhotonConversion1Leg_H
#define SeedForPhotonConversion1Leg_H

#include "RecoTracker/TkSeedingLayers/interface/SeedingHitSet.h"
#include "PrintRecoObjects.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/TrajectorySeed/interface/TrajectorySeedCollection.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "MagneticField/Engine/interface/MagneticField.h"

#include "RecoTracker/TransientTrackingRecHit/interface/TkTransientTrackingRecHitBuilder.h"


#include "FWCore/Utilities/interface/GCC11Compatibility.h"

class FreeTrajectoryState;

//
// this class need to be cleaned and optimized as those in RecoTracker/TkSeedGenerator
//
class dso_hidden SeedForPhotonConversion1Leg  {
public:
  static const int cotTheta_Max=99999;
  
  SeedForPhotonConversion1Leg( const edm::ParameterSet & cfg):
    thePropagatorLabel (cfg.getParameter<std::string>("propagator"))
    ,theBOFFMomentum   (cfg.getParameter<double>("SeedMomentumForBOFF"))
    ,TTRHBuilder       (cfg.getParameter<std::string>("TTRHBuilder"))
      {}

  //dtor
  ~SeedForPhotonConversion1Leg(){}

  const TrajectorySeed * trajectorySeed( TrajectorySeedCollection & seedCollection,
                         const SeedingHitSet & hits,
                         const GlobalPoint & vertex,
                         const GlobalVector & vertexBounds,
                         float ptmin,
                         const edm::EventSetup& es,
                         float cotTheta,
                         std::stringstream& ss);

  
 protected:

  bool checkHit(
            const TrajectoryStateOnSurface &,
            const SeedingHitSet::ConstRecHitPointer &hit,
            const edm::EventSetup& es) const { return true; }

  GlobalTrajectoryParameters initialKinematic(
                              const SeedingHitSet & hits, 
                              const GlobalPoint & vertexPos, 
                              const edm::EventSetup& es,
                              const float cotTheta) const;
  
  CurvilinearTrajectoryError initialError(
                          const GlobalVector& vertexBounds, 
                          float ptMin,  
                          float sinTheta) const;
  
  const TrajectorySeed * buildSeed(
                       TrajectorySeedCollection & seedCollection,
                       const SeedingHitSet & hits,
                       const FreeTrajectoryState & fts,
                       const edm::EventSetup& es) const;

  SeedingHitSet::RecHitPointer refitHit( SeedingHitSet::ConstRecHitPointer hit, 
                     const TrajectoryStateOnSurface &state, const TkClonerImpl& cloner) const;
  
protected:
  std::string thePropagatorLabel;
  double theBOFFMomentum;
  std::string TTRHBuilder;

  std::stringstream * pss;
  PrintRecoObjects po;
};
#endif