/data/refman/pasoursint/CMSSW_6_1_2_SLHC2/src/RecoTracker/RoadSearchCloudMaker/plugins/RoadSearchCloudMakerAlgorithm.h

Go to the documentation of this file.

#ifndef CloudMakerAlgorithm_h
#define CloudMakerAlgorithm_h

//
// Package:         RecoTracker/RoadSearchCloudMaker
// Class:           RoadSearchCloudMakerAlgorithm
// 
// Description:     
//                  Road categories determined by outer Seed RecHit
//                      RPhi: outer Seed RecHit in the Barrel
//                      ZPhi: outer Seed RecHit in the Disks
//                  use inner and outer Seed RecHit and BeamSpot to calculate extrapolation
//                      RPhi: phi = phi0 + asin(k r)
//                      ZPhi: phi = phi0 + C z
//                  Loop over RoadSet, access Rings of Road
//                      get average radius of Ring
//                      use extrapolation to calculate phi_ref at average Ring radius
//                      determine window in phi for DetId lookup in the Ring
//                              phi_ref ± phi_window
//                                      PARAMETER: phi_window = pi/24
//                      loop over DetId's in phi window in Ring
//                              two cases (problem of phi = 0/2pi):
//                                      lower window border < upper window border
//                                      upper window border < lower window border
//                  Loop over RecHits of DetId
//                      check compatibility of RecHit with extrapolation (delta phi cut)
//                      single layer sensor treatment
//                              RPhi:
//                                      stereo and barrel single layer sensor
//                                              calculate delta-phi
//                                      disk single layer sensor (r coordinate not well defined)
//                                              calculate delta phi between global positions of maximal strip extension and reference
//                              ZPhi:
//                                      stereo sensor
//                                              calculate delta-phi
//                                      barrel single layer sensor (z coordinate not well defined)
//                                              calculate delta phi between global positions of maximal strip extension and reference
//                                      disk single layer sensor (tilted strips relative to local coordinate system of sensor
//                                              calculate delta phi between global positions of maximal strip extension and reference
//                  Check delta phi cut
//                      cut value can be calculated based on extrapolation and Seed (pT dependent delta phi cut)
//                      currently: constant delta phi cut (PARAMETER)
//                              fill RecHit into Cloud
//                                      do not fill more than 32 RecHits per DetID into cloud (PARAMETER)
//                  first stage of Cloud cleaning cuts:
//                      number of layers with hits in cloud (PARAMETER)
//                      number of layers with no hits in cloud (PARAMETER)
//                      number of consecutive layers with no hits in cloud (PARAMETER)
//
// Original Author: Oliver Gutsche, gutsche@fnal.gov
// Created:         Sat Jan 14 22:00:00 UTC 2006
//
// $Author: eulisse $
// $Date: 2012/10/24 08:10:40 $
// $Revision: 1.1 $
//

#include <string>
#include <sstream>

#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/EventSetup.h"

#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/RoadSearchSeed/interface/RoadSearchSeedCollection.h"
#include "DataFormats/TrackerRecHit2D/interface/SiStripRecHit2DCollection.h"
#include "DataFormats/TrackerRecHit2D/interface/SiPixelRecHitCollection.h"
#include "DataFormats/TrackerRecHit2D/interface/SiStripMatchedRecHit2DCollection.h"

#include "DataFormats/RoadSearchCloud/interface/RoadSearchCloudCollection.h"
#include "DataFormats/RoadSearchCloud/interface/RoadSearchCloud.h"
#include "RecoTracker/RoadMapRecord/interface/Roads.h"

#include "TrackingTools/RoadSearchHitAccess/interface/DetHitAccess.h"

#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"

#include "RecoLocalTracker/SiStripRecHitConverter/interface/SiStripRecHitMatcher.h"
 
class RoadSearchCloudMakerAlgorithm 
{
 public:
  
  RoadSearchCloudMakerAlgorithm(const edm::ParameterSet& conf);
  ~RoadSearchCloudMakerAlgorithm();

  void run(edm::Handle<RoadSearchSeedCollection> input,
           const SiStripRecHit2DCollection* rphiRecHits,
           const SiStripRecHit2DCollection* stereoRecHits,
           const SiStripMatchedRecHit2DCollection* matchedRecHits,
           const SiPixelRecHitCollection *pixRecHits,
           const edm::EventSetup& es,
           RoadSearchCloudCollection &output);

  unsigned int FillRecHitsIntoCloudGeneral(DetId id, double d0, double phi0, double k0, double phi1, double k1,
                                           Roads::type roadType, double ringPhi,
                                           const TrackerGeometry *tracker, const SiStripRecHitMatcher* theHitMatcher, RoadSearchCloud &cloud);

  unsigned int FillRecHitsIntoCloud(DetId id, const SiStripRecHit2DCollection* inputRecHits, 
                                    double d0, double phi0, double k0, Roads::type roadType, double ringPhi,
                                    const TrackerGeometry *tracker, RoadSearchCloud &cloud);

  
  unsigned int FillPixRecHitsIntoCloud(DetId id, 
                                       const SiPixelRecHitCollection *inputRecHits, 
                                       double d0, double phi0, double k0, Roads::type roadType, double ringPhi,
                                       const TrackerGeometry *tracker, RoadSearchCloud &cloud);

  bool isSingleLayer(DetId id);

  bool isBarrelSensor(DetId id);

  double phiFromExtrapolation(double d0, double phi0, double k0, double ringRadius, Roads::type roadType);

  double phiMax(Roads::type roadType, double phi0, double k0);

  double map_phi(double phi);
  double map_phi2(double phi);

  void makecircle(double x1_cs, double y1_cs,double x2_cs, double y2_cs,
                                             double x3_cs, double y3_cs);

  double CheckXYIntersection(LocalPoint& ip1, LocalPoint& op1, LocalPoint& ip2, LocalPoint& op2);

  double CheckZPhiIntersection(double iPhi1, double iZ1, double oPhi1, double oZ1,
                               double iPhi2, double iZ2, double oPhi2, double oZ2);

  double ZPhiDeltaPhi(double phi1, double phi2, double phiExpect);

  RoadSearchCloudCollection Clean(RoadSearchCloudCollection *rawColl);

  SiStripMatchedRecHit2D*  CorrectMatchedHit(const TrackingRecHit* originalRH,
                                             const GluedGeomDet* gluedDet,
                                             const TrackerGeometry *tracker,
                                             const SiStripRecHitMatcher* theHitMatcher,
                                             double k0, double phi0);

 private:

  edm::ParameterSet conf_;
  static double epsilon;
  double d0h, phi0h, omegah;
  double rphicsq;
  int rphinhits;
  const SiPixelRecHitCollection thePixRecHits;
  
  // general hit access for road search
  DetHitAccess recHitVectorClass;
  
  double theRPhiRoadSize;
  double theZPhiRoadSize;
  double theMinimumHalfRoad;
  bool UsePixels;
  bool NoFieldCosmic;
  unsigned int maxDetHitsInCloudPerDetId;
  double       minFractionOfUsedLayersPerCloud;
  double       maxFractionOfMissedLayersPerCloud;
  double       maxFractionOfConsecutiveMissedLayersPerCloud;
  unsigned int increaseMaxNumberOfConsecutiveMissedLayersPerCloud;
  unsigned int increaseMaxNumberOfMissedLayersPerCloud;

  bool doCleaning_;
  double mergingFraction_;
  unsigned int maxRecHitsInCloud_;

  std::ostringstream output_;
  double scalefactorRoadSeedWindow_;

  std::string roadsLabel_;

};

#endif