/data/refman/pasoursint/CMSSW_5_3_1/src/RecoTracker/ConversionSeedGenerators/interface/SeedForPhotonConversionFromQuadruplets.h

Go to the documentation of this file.

#ifndef SeedForPhotonConversionFromQuadruplets_H
#define SeedForPhotonConversionFromQuadruplets_H

#include "RecoTracker/TkTrackingRegions/interface/GlobalTrackingRegion.h"
#include "RecoTracker/TkSeedGenerator/interface/SeedCreator.h"
#include "RecoTracker/TkSeedingLayers/interface/SeedingHitSet.h"
#include "RecoTracker/ConversionSeedGenerators/interface/PrintRecoObjects.h"
#include "RecoTracker/ConversionSeedGenerators/interface/Quad.h"
class FreeTrajectoryState;

class SeedForPhotonConversionFromQuadruplets : public SeedCreator {
public:
  static const int cotTheta_Max=99999;
  
  SeedForPhotonConversionFromQuadruplets( const edm::ParameterSet & cfg):
    thePropagatorLabel(cfg.getParameter<std::string>("propagator")),
    theBOFFMomentum(cfg.existsAs<double>("SeedMomentumForBOFF") ? cfg.getParameter<double>("SeedMomentumForBOFF") : 5.0)
      {}

  SeedForPhotonConversionFromQuadruplets( 
      const std::string & propagator = "PropagatorWithMaterial", double seedMomentumForBOFF = -5.0) 
   : thePropagatorLabel(propagator), theBOFFMomentum(seedMomentumForBOFF) { }

  //dtor
  virtual ~SeedForPhotonConversionFromQuadruplets(){}

  virtual const TrajectorySeed * trajectorySeed( TrajectorySeedCollection & seedCollection,
                                                 const SeedingHitSet & phits,
                                                 const SeedingHitSet & mhits,
                                                 const TrackingRegion & region,
                                                 const edm::EventSetup& es,
                                                 std::stringstream& ss, std::vector<Quad> & quadV);

  virtual const TrajectorySeed *trajectorySeed(
                                               TrajectorySeedCollection & seedCollection,
                                               const SeedingHitSet & hits,
                                               const TrackingRegion & region,
                                               const edm::EventSetup& es,
                                               const SeedComparitor *filter){ return 0;}
  
  double simpleGetSlope(const TransientTrackingRecHit::ConstRecHitPointer &ohit, const TransientTrackingRecHit::ConstRecHitPointer &nohit, const TransientTrackingRecHit::ConstRecHitPointer &ihit, const TransientTrackingRecHit::ConstRecHitPointer &nihit, const TrackingRegion & region, double & cotTheta, double & z0);
  double verySimpleFit(int size, double* ax, double* ay, double* e2y, double& p0, double& e2p0, double& p1);
  double getSqrEffectiveErrorOnZ(const TransientTrackingRecHit::ConstRecHitPointer &hit, const TrackingRegion & region);

  //
  // Some utility methods added by sguazz
  void stupidPrint(std::string s,float* d);
  void stupidPrint(std::string s,double* d);
  void stupidPrint(const char* s,GlobalPoint* d);
  void stupidPrint(const char* s,GlobalPoint* d, int n);
  void bubbleSortVsPhi(GlobalPoint arr[], int n, GlobalPoint vtx);
  void bubbleReverseSortVsPhi(GlobalPoint arr[], int n, GlobalPoint vtx);
  //
  //




 protected:

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

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

  virtual TransientTrackingRecHit::RecHitPointer refitHit(
                                                          const TransientTrackingRecHit::ConstRecHitPointer &hit, 
                                                          const TrajectoryStateOnSurface &state) const;

  bool similarQuadExist(Quad & thisQuad, std::vector<Quad>& quadV);

  
protected:
  std::string thePropagatorLabel;
  double theBOFFMomentum;

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