FWLiteJetProducer.h

Go to the documentation of this file.

#ifndef Demo_PFRootEvent_FWLiteJetProducer_h
#define Demo_PFRootEvent_FWLiteJetProducer_h

#include "DataFormats/JetReco/interface/BasicJet.h"
#include "DataFormats/JetReco/interface/CaloJet.h"
#include "DataFormats/JetReco/interface/GenJet.h"
#include "DataFormats/JetReco/interface/BasicJetCollection.h"
#include "DataFormats/JetReco/interface/CaloJetCollection.h"
#include "DataFormats/JetReco/interface/GenJetCollection.h"
#include "DataFormats/RecoCandidate/interface/RecoCaloTowerCandidate.h" 
#include "RecoParticleFlow/PFRootEvent/interface/CMSIterativeConeAlgorithm.h"
#include "RecoParticleFlow/PFRootEvent/interface/CMSMidpointAlgorithm.h"
#include "RecoParticleFlow/PFRootEvent/interface/FastJetFWLiteWrapper.h"
#include "RecoParticleFlow/PFRootEvent/interface/JetRecoTypes.h"

// #include "FWCore/Framework/interface/OrphanHandle.h"
#include "DataFormats/Common/interface/OrphanHandle.h"
// #include "DataFormats/Common/interface/ProductID.h"
#include "DataFormats/Provenance/interface/ProductID.h"

class Utils;
class FWLiteJetProducer{
 public:
  FWLiteJetProducer();  
  ~FWLiteJetProducer();
  void applyCuts( const reco::CandidatePtrVector& Candidates, 
                  JetReco::InputCollection* input );         
  void makeIterativeConeJets(const JetReco::InputCollection& fInput, 
                             JetReco::OutputCollection* fOutput);
  void makeFastJets(const JetReco::InputCollection& fInput, 
                    JetReco::OutputCollection* fOutput);
  void makeMidpointJets(const JetReco::InputCollection& fInput, 
                        JetReco::OutputCollection* fOutput);
  void print();
  void updateParameter();
  CMSIterativeConeAlgorithm* algoIC_;   
  FastJetFWLiteWrapper algoF_;
  CMSMidpointAlgorithm* algoMC_;

  // Get methods --------------------------------------------
  double getmEtInputCut(){return mEtInputCut_;}
  double getmEInputCut(){return mEInputCut_;}
  double getSeedThreshold(){return seedThreshold_;}
  double  getConeRadius(){return coneRadius_;}
  double getConeAreaFraction(){return coneAreaFraction_;}
  int getMaxPairSize(){return maxPairSize_;}
  int getMaxIterations(){return maxIterations_;} 
  double getOverlapThreshold(){return overlapThreshold_;}
  double getPtMin (){return ptMin_ ;}
  double getRParam(){return  rparam_;}
  bool getDebug(){return  debug_;}

  // Set methods --------------------------------------------   
  void setmEtInputCut (double amEtInputCut){mEtInputCut_=amEtInputCut;}
  void setmEInputCut (double amEInputCut){mEInputCut_=amEInputCut;}
  void setSeedThreshold(double aSeedThreshold) {seedThreshold_=aSeedThreshold;}
  void setConeRadius(double aConeRadius) {coneRadius_=aConeRadius;}
  void setConeAreaFraction(double aConeAreaFraction) {coneAreaFraction_=aConeAreaFraction;}
  void setMaxPairSize(int aMaxPairSize) {maxPairSize_=aMaxPairSize;}
  void setMaxIterations(int amaxIteration) {maxIterations_=amaxIteration;}
  void setOverlapThreshold(double aOverlapThreshold) {overlapThreshold_=aOverlapThreshold;}
  void setPtMin (double aPtMin){ptMin_=aPtMin;}
  void setRParam (double aRparam){rparam_=aRparam;}
  void setDebug(bool aDebug){debug_=aDebug;}  

  // jets parameters ----------------------------------------
 private:       
  double mEtInputCut_;
  double mEInputCut_ ;
  double seedThreshold_;
  double coneRadius_;
  double coneAreaFraction_;
  int maxPairSize_;
  int maxIterations_;
  double overlapThreshold_;
  double ptMin_;
  double rparam_;
  bool debug_;
  
};
#endif