---

GlobalGsfElectronAlgo.h

Go to the documentation of this file.

#ifndef GlobalGsfElectronAlgo_H
#define GlobalGsfElectronAlgo_H

#include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
#include "DataFormats/EgammaReco/interface/SuperClusterFwd.h"
#include "DataFormats/GsfTrackReco/interface/GsfTrackFwd.h"
#include "DataFormats/TrajectorySeed/interface/TrajectorySeedCollection.h"

#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/Event.h"

#include "MagneticField/Engine/interface/MagneticField.h"
#include "TrackingTools/MaterialEffects/interface/PropagatorWithMaterial.h"
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
#include "RecoCaloTools/MetaCollections/interface/CaloRecHitMetaCollections.h"

#include "RecoTracker/MeasurementDet/interface/MeasurementTracker.h"
#include "Geometry/CaloGeometry/interface/CaloGeometry.h"

class MultiTrajectoryStateTransform;
class GsfPropagatorAdapter;
 
class GlobalGsfElectronAlgo {

public:

  GlobalGsfElectronAlgo(const edm::ParameterSet& conf,
                         double maxEOverPBarrel, double maxEOverPEndcaps, 
                         double minEOverPBarrel, double minEOverPEndcaps,
                         double maxDeltaEta, double maxDeltaPhi, 
                         bool highPtPresel, double highPtMin,
                         bool applyEtaCorrection, bool applyAmbResolution);
  ~GlobalGsfElectronAlgo();

  void setupES(const edm::EventSetup& setup);
  void run(edm::Event&, reco::GsfElectronCollection&);

 private:

  // create electrons from superclusters, tracks and Hcal rechits
  void process(edm::Handle<reco::GsfTrackCollection> tracksH,
               edm::Handle<reco::SuperClusterCollection> superClustersBarrelH,
               edm::Handle<reco::SuperClusterCollection> superClustersEndcapH,
               const math::XYZPoint &bs,
               reco::GsfElectronCollection & outEle);
  
  
  // preselection method
  bool preSelection(const reco::SuperCluster& clus);

  // interface to be improved...
  void createElectron(const reco::SuperClusterRef & scRef,
                      const reco::GsfTrackRef &trackRef,
                      reco::GsfElectronCollection & outEle);  

  void resolveElectrons(std::vector<reco::GsfElectron> &, reco::GsfElectronCollection & outEle);
  
  //Gsf mode calculations
  GlobalVector computeMode(const TrajectoryStateOnSurface &tsos);

  const reco::GsfTrackRef
    superClusterMatching(reco::SuperClusterRef sc, edm::Handle<reco::GsfTrackCollection> tracks);

  // intermediate calculations
  bool calculateTSOS(const reco::GsfTrack &t,const reco::SuperCluster & theClus, const math::XYZPoint & bs);

  // preselection parameters
  // maximum E/p where E is the supercluster corrected energy and p the track momentum at innermost state  
  double maxEOverPBarrel_;   
  double maxEOverPEndcaps_;   
  // minimum E/p where E is the supercluster corrected energy and p the track momentum at innermost state  
  double minEOverPBarrel_;   
  double minEOverPEndcaps_;     
  // maximum eta difference between the supercluster position and the track position at the closest impact to the supercluster 
  double maxDeltaEta_;
  // maximum phi difference between the supercluster position and the track position at the closest impact to the supercluster
  // position to the supercluster
  double maxDeltaPhi_;

  // high pt preselection parameters
  bool highPtPreselection_;
  double highPtMin_;
  
  // if this parameter is false, only SC level Escale correctoins are applied
  bool applyEtaCorrection_;
  
  // if this parameter is true, "double" electrons are resolved
  bool applyAmbResolution_;

  // input configuration
  edm::InputTag barrelSuperClusters_;
  edm::InputTag endcapSuperClusters_;
  edm::InputTag tracks_;
  edm::InputTag hcalRecHits_;

  edm::ESHandle<MagneticField>                theMagField;
  edm::ESHandle<CaloGeometry>                 theCaloGeom;
  edm::ESHandle<TrackerGeometry>              trackerHandle_;

  const MultiTrajectoryStateTransform *mtsTransform_;
  const GsfPropagatorAdapter *geomPropBw_;
  const GsfPropagatorAdapter *geomPropFw_;

  // internal variables 
  int subdet_; //subdetector for this cluster
  GlobalPoint sclPos_;
  GlobalVector vtxMom_;
  TrajectoryStateOnSurface innTSOS_;
  TrajectoryStateOnSurface outTSOS_;
  TrajectoryStateOnSurface vtxTSOS_;
  TrajectoryStateOnSurface sclTSOS_;
  TrajectoryStateOnSurface seedTSOS_;

  HBHERecHitMetaCollection *mhbhe_;

  unsigned long long cacheIDGeom_;
  unsigned long long cacheIDTDGeom_;
  unsigned long long cacheIDMagField_;
};

#endif // GlobalGsfElectronAlgo_H

---