pat::PATMHTProducer Class Reference

Inheritance diagram for pat::PATMHTProducer:

Classes
class	uncertaintyFunctions

Public Member Functions
	PATMHTProducer (const edm::ParameterSet &)
	~PATMHTProducer () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Private Member Functions
double	getElectrons (edm::Event &, const edm::EventSetup &)
double	getJets (edm::Event &, const edm::EventSetup &)
double	getMuons (edm::Event &, const edm::EventSetup &)
void	getTowers (edm::Event &, const edm::EventSetup &)
void	produce (edm::Event &, const edm::EventSetup &) override
void	setUncertaintyParameters ()

Private Attributes
edm::InputTag	CaloJetAlgorithmTag_
edm::InputTag	CaloTowerTag_
bool	controlledUncertainty_
edm::InputTag	CorJetAlgorithmTag_
uncertaintyFunctions	ecalEBUncertainty
uncertaintyFunctions	ecalEEUncertainty
edm::InputTag	ElectronTag_
double	eleEtaMax_
double	eleEtUncertaintyParameter0_
double	elePhiUncertaintyParameter0_
double	elePtMin_
edm::EDGetTokenT< edm::View < pat::Electron > >	eleToken_
uncertaintyFunctions	eleUncertainty
uncertaintyFunctions	hcalHBUncertainty
uncertaintyFunctions	hcalHEUncertainty
uncertaintyFunctions	hcalHFUncertainty
uncertaintyFunctions	hcalHOUncertainty
std::string	JetCorrectionService_
uncertaintyFunctions	jetCorrUncertainty
double	jetEMfracMax_
double	jetEtaMax_
double	jetEtUncertaintyParameter0_
double	jetEtUncertaintyParameter1_
double	jetEtUncertaintyParameter2_
double	jetPhiUncertaintyParameter0_
double	jetPhiUncertaintyParameter1_
double	jetPhiUncertaintyParameter2_
double	jetPtMin_
edm::EDGetTokenT< edm::View < pat::Jet > >	jetToken_
uncertaintyFunctions	jetUncertainty
std::string	metCollectionLabel_
edm::InputTag	mhtLabel_
double	muonChiSqMax_
uncertaintyFunctions	muonCorrUncertainty
double	muonDPtMax_
double	muonEtaMax_
double	muonEtUncertaintyParameter0_
int	muonNHitsMin_
double	muonPhiUncertaintyParameter0_
double	muonPtMin_
edm::InputTag	MuonTag_
double	muonTrackD0Max_
double	muonTrackDzMax_
uncertaintyFunctions	muonUncertainty
edm::EDGetTokenT< edm::View < pat::Muon > >	muoToken_
bool	noHF_
edm::EDGetTokenT< edm::View < pat::Photon > >	phoToken_
std::vector< metsig::SigInputObj >	physobjvector_
std::set< CaloTowerDetId >	s_clusteredTowers
std::string	significanceLabel_
edm::EDGetTokenT< edm::View < pat::Tau > >	tauToken_
double	towerEtThreshold_
double	uncertaintyScaleFactor_
bool	useCaloTowers_
bool	useElectrons_
bool	useHO_
bool	useJets_
bool	useMuons_
double	verbose_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T...>
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T...>
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 39 of file PATMHTProducer.cc.

Constructor & Destructor Documentation

PATMHTProducer::PATMHTProducer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 154 of file PATMHTProducer.cc.

References CaloTowerTag_, controlledUncertainty_, eleEtaMax_, eleEtUncertaintyParameter0_, elePhiUncertaintyParameter0_, elePtMin_, eleToken_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), jetEMfracMax_, jetEtaMax_, jetEtUncertaintyParameter0_, jetEtUncertaintyParameter1_, jetEtUncertaintyParameter2_, jetPhiUncertaintyParameter0_, jetPhiUncertaintyParameter1_, jetPhiUncertaintyParameter2_, jetPtMin_, jetToken_, muonEtaMax_, muonEtUncertaintyParameter0_, muonPhiUncertaintyParameter0_, muonPtMin_, muoToken_, noHF_, phoToken_, setUncertaintyParameters(), tauToken_, towerEtThreshold_, uncertaintyScaleFactor_, useHO_, and verbose_.

                                                                {
  // Initialize the configurables
  verbose_ = iConfig.getParameter<double>("verbose");

  jetToken_ = consumes<edm::View<pat::Jet> >(iConfig.getUntrackedParameter<edm::InputTag>("jetTag"));
  eleToken_ = consumes<edm::View<pat::Electron> >(iConfig.getUntrackedParameter<edm::InputTag>("electronTag"));
  muoToken_ = consumes<edm::View<pat::Muon> >(iConfig.getUntrackedParameter<edm::InputTag>("muonTag"));
  tauToken_ = consumes<edm::View<pat::Tau> >(iConfig.getUntrackedParameter<edm::InputTag>("tauTag"));
  phoToken_ = consumes<edm::View<pat::Photon> >(iConfig.getUntrackedParameter<edm::InputTag>("photonTag"));

  uncertaintyScaleFactor_ = iConfig.getParameter<double>("uncertaintyScaleFactor");
  controlledUncertainty_ = iConfig.getParameter<bool>("controlledUncertainty");

  jetPtMin_ = iConfig.getParameter<double>("jetPtMin");
  jetEtaMax_ = iConfig.getParameter<double>("jetEtaMax");
  jetEMfracMax_ = iConfig.getParameter<double>("jetEMfracMax");
  elePtMin_ = iConfig.getParameter<double>("elePtMin");
  eleEtaMax_ = iConfig.getParameter<double>("eleEtaMax");
  muonPtMin_ = iConfig.getParameter<double>("muonPtMin");
  muonEtaMax_ = iConfig.getParameter<double>("muonEtaMax");

  jetEtUncertaintyParameter0_ = iConfig.getParameter<double>("jetEtUncertaintyParameter0");
  jetEtUncertaintyParameter1_ = iConfig.getParameter<double>("jetEtUncertaintyParameter1");
  jetEtUncertaintyParameter2_ = iConfig.getParameter<double>("jetEtUncertaintyParameter2");
  jetPhiUncertaintyParameter0_ = iConfig.getParameter<double>("jetPhiUncertaintyParameter0");
  jetPhiUncertaintyParameter1_ = iConfig.getParameter<double>("jetPhiUncertaintyParameter1");
  jetPhiUncertaintyParameter2_ = iConfig.getParameter<double>("jetPhiUncertaintyParameter2");

  eleEtUncertaintyParameter0_ = iConfig.getParameter<double>("eleEtUncertaintyParameter0");
  elePhiUncertaintyParameter0_ = iConfig.getParameter<double>("elePhiUncertaintyParameter0");

  muonEtUncertaintyParameter0_ = iConfig.getParameter<double>("muonEtUncertaintyParameter0");
  muonPhiUncertaintyParameter0_ = iConfig.getParameter<double>("muonPhiUncertaintyParameter0");

  CaloTowerTag_ = iConfig.getParameter<edm::InputTag>("CaloTowerTag");
  noHF_ = iConfig.getParameter<bool>("noHF");

  //  muonCalo_ = iConfig.getParameter<bool>("muonCalo");
  towerEtThreshold_ = iConfig.getParameter<double>("towerEtThreshold");
  useHO_ = iConfig.getParameter<bool>("useHO");

  setUncertaintyParameters();

  produces<pat::MHTCollection>();
}

PATMHTProducer::~PATMHTProducer ( )

override

Definition at line 200 of file PATMHTProducer.cc.

{}

Member Function Documentation

double PATMHTProducer::getElectrons ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

private

Definition at line 343 of file PATMHTProducer.cc.

References edm::View< T >::begin(), gather_cfg::cout, HI_PhotonSkim_cff::electrons, edm::View< T >::end(), edm::Event::getByToken(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                    {
  std::string objectname = "electron";

  double number_of_electrons_ = 0.0;

  // edm::ESHandle<CaloTowerConstituentsMap> cttopo;
  // iSetup.get<HcalRecNumberingRecord>().get(cttopo);
  // const CaloTowerConstituentsMap* caloTowerMap = cttopo.product();

  edm::Handle<edm::View<pat::Electron> > electronHandle;
  iEvent.getByToken(eleToken_, electronHandle);
  edm::View<pat::Electron> electrons = *electronHandle;

  // Fill Input Vector with Electrons
  for (edm::View<pat::Electron>::const_iterator electron_iter = electrons.begin(); electron_iter != electrons.end();
       ++electron_iter) {
    // Select electrons
    if (electron_iter->et() < elePtMin_ || TMath::Abs(electron_iter->eta()) > eleEtaMax_)
      continue;

    if (verbose_ == 3.) {
      std::cout << "electron pt = " << electron_iter->pt() << " eta : " << electron_iter->eta() << std::endl;
    }

    double electron_et = electron_iter->et();
    double electron_phi = electron_iter->phi();

    double sigma_et, sigma_phi;

    if (controlledUncertainty_) {
      sigma_et = eleUncertainty.etUncertainty->Eval(electron_et);
      sigma_phi = eleUncertainty.phiUncertainty->Eval(electron_et);
    } else {
      sigma_et = 0.0;   //electron_iter->resolutionEt();
      sigma_phi = 0.0;  // electron_iter->resolutionPhi();
    }

    if (verbose_ == 3.) {
      std::cout << "electron sigma_et : " << sigma_et << ", electron sigma_phi : " << sigma_phi << std::endl;
    }

    if (sigma_et < 0 || sigma_phi < 0)
      edm::LogWarning("PATMHTProducer") << " uncertainties for " << objectname << " are (et, phi): " << sigma_et << ","
                                        << sigma_phi << " (et,phi): " << electron_et << "," << electron_phi;

    if (uncertaintyScaleFactor_ != 1.0) {
      sigma_et = sigma_et * uncertaintyScaleFactor_;
      sigma_phi = sigma_phi * uncertaintyScaleFactor_;
    }

    metsig::SigInputObj tmp_electron(objectname, electron_et, electron_phi, sigma_et, sigma_phi);
    physobjvector_.push_back(tmp_electron);
    number_of_electrons_++;

    //-- Store tower DetId's to be removed from Calo Tower sum later --//
    /*
    const reco::SuperCluster& eleSC = *( electron_iter->superCluster() );

    std::vector<DetId> v_eleDetIds = eleSC.getHitsByDetId();

    //-- Convert cells to calo towers and add to set --//
    for( std::vector<DetId>::iterator cellId = v_eleDetIds.begin();
         cellId != v_eleDetIds.end();
         cellId++) {

      CaloTowerDetId towerId = caloTowerMap->towerOf(*cellId);
      if (towerId != nullDetId) {
    //std::cout << ">>> electron towerId: " << towerId << std::endl;
    std::pair<std::_Rb_tree_const_iterator<CaloTowerDetId>,bool> p1 = s_clusteredTowers.insert(towerId);
      }
      else
    std::cerr<<"No matching tower found for electron cell!\n";
    }

    */
  }

  return number_of_electrons_;
}

double PATMHTProducer::getJets ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

private

Definition at line 263 of file PATMHTProducer.cc.

References edm::View< T >::begin(), gather_cfg::cout, edm::View< T >::end(), edm::Event::getByToken(), cuy::ii, fwrapper::jets, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                               {
  std::string objectname = "jet";

  double number_of_jets_ = 0.0;

  edm::Handle<edm::View<pat::Jet> > jetHandle;
  iEvent.getByToken(jetToken_, jetHandle);
  edm::View<pat::Jet> jets = *jetHandle;

  // Fill Input Vector with Jets
  for (edm::View<pat::Jet>::const_iterator jet_iter = jets.begin(); jet_iter != jets.end(); ++jet_iter) {
    if ((jet_iter->pt() < jetPtMin_) || (TMath::Abs(jet_iter->eta()) > jetEtaMax_) ||
        (jet_iter->emEnergyFraction() > jetEMfracMax_))
      continue;

    double jet_et = jet_iter->et();
    double jet_phi = jet_iter->phi();

    if (verbose_ == 3.) {
      std::cout << "jet pt : " << jet_iter->pt() << " eta : " << jet_iter->eta()
                << " EMF: " << jet_iter->emEnergyFraction() << std::endl;
    }

    double sigma_et, sigma_phi;

    if (controlledUncertainty_) {
      sigma_et = jetUncertainty.etUncertainty->Eval(jet_et);
      sigma_phi = jetUncertainty.phiUncertainty->Eval(jet_et);
    } else {
      sigma_et = 0.0;   // jet_iter->resolutionEt();
      sigma_phi = 0.0;  //jet_iter->resolutionPhi();
    }

    if (verbose_ == 3.) {
      std::cout << "jet sigma_et : " << sigma_et << ", jet sigma_phi : " << sigma_phi << std::endl;
    }

    if (sigma_et <= 0 || sigma_phi <= 0)
      edm::LogWarning("PATMHTProducer") << " uncertainties for " << objectname << " are (et, phi): " << sigma_et << ","
                                        << sigma_phi << " (et,phi): " << jet_et << "," << jet_phi;
    // try to read out the jet resolution from the root file at PatUtils
    //-- Store jet for Significance Calculation --//

    if (uncertaintyScaleFactor_ != 1.0) {
      sigma_et = sigma_et * uncertaintyScaleFactor_;
      sigma_phi = sigma_phi * uncertaintyScaleFactor_;
      // edm::LogWarning("PATMHTProducer") << " using uncertainty scale factor: " << uncertaintyScaleFactor_ <<
      //" , uncertainties for " << objectname <<" changed to (et, phi): " << sigma_et << "," << sigma_phi;
    }

    if (verbose_ == 101.) {  // Study the Jets behavior

      std::cout << "v101> " << number_of_jets_ << "  " << jet_et << "  " << sigma_et << "  " << jet_phi << "  "
                << sigma_phi << std::endl;
    }

    metsig::SigInputObj tmp_jet(objectname, jet_et, jet_phi, sigma_et, sigma_phi);
    physobjvector_.push_back(tmp_jet);
    number_of_jets_++;

    //-- Store tower DetId's to be removed from Calo Tower sum later --//
    std::vector<CaloTowerPtr> v_towers = jet_iter->getCaloConstituents();
    //std::cout << "tower size = " << v_towers.size() << std::endl;

    for (unsigned int ii = 0; ii < v_towers.size(); ii++) {
      s_clusteredTowers.insert((*v_towers.at(ii)).id());
      //std::cout << "tower id = " << (*v_towers.at(ii)).id() << std::endl;
    }
  }

  if (verbose_ == 101.) {  // Study the Jets behavior - seperate events
    std::cout << "v101> --------------------------------------------" << std::endl;
  }

  return number_of_jets_;
}

double PATMHTProducer::getMuons ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

private

Definition at line 427 of file PATMHTProducer.cc.

References edm::View< T >::begin(), gather_cfg::cout, edm::View< T >::end(), edm::Event::getByToken(), edm::HandleBase::isValid(), patZpeak::muons, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                {
  std::string objectname = "muon";
  edm::Handle<edm::View<pat::Muon> > muonHandle;
  iEvent.getByToken(muoToken_, muonHandle);
  edm::View<pat::Muon> muons = *muonHandle;

  if (!muonHandle.isValid()) {
    std::cout << ">>> PATMHTSelector not valid muon Handle!" << std::endl;
    return 0.0;
  }

  double number_of_muons_ = 0.0;

  for (edm::View<pat::Muon>::const_iterator muon_iter = muons.begin(); muon_iter != muons.end(); ++muon_iter) {
    if (muon_iter->pt() < muonPtMin_ || TMath::Abs(muon_iter->eta()) > muonEtaMax_)
      continue;

    if (verbose_ == 3.) {
      std::cout << "muon pt = " << muon_iter->pt() << " eta : " << muon_iter->eta() << std::endl;
    }

    double muon_pt = muon_iter->pt();
    double muon_phi = muon_iter->phi();

    double sigma_et, sigma_phi;

    if (controlledUncertainty_) {
      sigma_et = muonUncertainty.etUncertainty->Eval(muon_pt);
      sigma_phi = muonUncertainty.phiUncertainty->Eval(muon_pt);
    } else {
      sigma_et = 0.0;   //muon_iter->resolutionEt();
      sigma_phi = 0.0;  // muon_iter->resolutionPhi();
    }

    if (verbose_ == 3.) {
      std::cout << "muon sigma_et : " << sigma_et << ", muon sigma_phi : " << sigma_phi << std::endl;
    }

    if (sigma_et < 0 || sigma_phi < 0)
      edm::LogWarning("PATMHTProducer") << " uncertainties for " << objectname << " are (et, phi): " << sigma_et << ","
                                        << sigma_phi << " (pt,phi): " << muon_pt << "," << muon_phi;

    if (uncertaintyScaleFactor_ != 1.0) {
      sigma_et = sigma_et * uncertaintyScaleFactor_;
      sigma_phi = sigma_phi * uncertaintyScaleFactor_;
    }

    metsig::SigInputObj tmp_muon(objectname, muon_pt, muon_phi, sigma_et, sigma_phi);
    physobjvector_.push_back(tmp_muon);
    number_of_muons_++;

  }  // end Muon loop

  return number_of_muons_;
}

void pat::PATMHTProducer::getTowers ( edm::Event &  , const edm::EventSetup &    )

private

void PATMHTProducer::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 202 of file PATMHTProducer.cc.

References metsig::significanceAlgo::addObjects(), funct::cos(), gather_cfg::cout, eostools::move(), edm::Event::put(), pat::MHT::setNumberOfElectrons(), pat::MHT::setNumberOfJets(), pat::MHT::setNumberOfMuons(), metsig::significanceAlgo::significance(), and funct::sin().

                                                                             {
  // make sure the SigInputObj container is empty
  while (!physobjvector_.empty()) {
    physobjvector_.erase(physobjvector_.begin(), physobjvector_.end());
  }

  // Clean the clustered towers
  s_clusteredTowers.clear();

  double number_of_jets = getJets(iEvent, iSetup);

  double number_of_electrons = getElectrons(iEvent, iSetup);

  double number_of_muons = getMuons(iEvent, iSetup);

  if (verbose_ == 1.) {
    std::cout << ">>>---> Number of jets: " << number_of_jets << std::endl;
    std::cout << ">>>---> Number of electrons: " << number_of_jets << std::endl;
    std::cout << ">>>---> Number of muons: " << number_of_muons << std::endl;
  }

  double met_x = 0;
  double met_y = 0;
  double met_et = 0;
  double met_phi = 0;
  double met_set = 0;

  auto themetsigcoll = std::make_unique<pat::MHTCollection>();

  if (!physobjvector_.empty()) {  // Only when the vector is not empty

    // calculate the MHT significance

    metsig::significanceAlgo signifAlgo;
    signifAlgo.addObjects(physobjvector_);
    double significance = signifAlgo.significance(met_et, met_phi, met_set);

    met_x = met_et * cos(met_phi);
    met_y = met_et * sin(met_phi);

    if (verbose_ == 1.) {
      std::cout << ">>>----> MHT Sgificance = " << significance << std::endl;
    }

    pat::MHT themetsigobj(reco::Particle::LorentzVector(met_x, met_y, 0, met_et), met_set, significance);

    // Store the number of jets, electrons, muons
    themetsigobj.setNumberOfJets(number_of_jets);
    themetsigobj.setNumberOfElectrons(number_of_electrons);
    themetsigobj.setNumberOfMuons(number_of_muons);

    themetsigcoll->push_back(themetsigobj);

  }  // If the vector is empty, just put empty product.

  iEvent.put(std::move(themetsigcoll));
}

void PATMHTProducer::setUncertaintyParameters ( )

private

Definition at line 484 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

                                                 {
  // set the various functions here:

  //-- For Et functions, [0]= par_n, [1]=par_s, [2]= par_c ---//
  //-- Ecal Uncertainty Functions ------------------------------------//
  //-- From: FastSimulation/Calorimetry/data/HcalResponse.cfi --//
  //-- Ecal Barrel --//
  ecalEBUncertainty.etUncertainty =
      std::make_unique<TF1>("ecalEBEtFunc", "x*sqrt(([0]*[0]/(x*x))+([1]*[1]/x)+([2]*[2]))", 3);
  ecalEBUncertainty.etUncertainty->SetParameter(0, 0.2);
  ecalEBUncertainty.etUncertainty->SetParameter(1, 0.03);
  ecalEBUncertainty.etUncertainty->SetParameter(2, 0.005);

  ecalEBUncertainty.phiUncertainty = std::make_unique<TF1>("ecalEBphiFunc", "[0]*x", 1);
  ecalEBUncertainty.phiUncertainty->SetParameter(0, 0.0174);

  //-- Ecal Endcap --//
  ecalEEUncertainty.etUncertainty =
      std::make_unique<TF1>("ecalEEEtFunc", "x*sqrt(([0]*[0]/(x*x))+([1]*[1]/x)+([2]*[2]))", 3);
  ecalEEUncertainty.etUncertainty->SetParameter(0, 0.2);
  ecalEEUncertainty.etUncertainty->SetParameter(1, 0.03);
  ecalEEUncertainty.etUncertainty->SetParameter(2, 0.005);

  ecalEEUncertainty.phiUncertainty = std::make_unique<TF1>("ecalEEphiFunc", "[0]*x", 1);
  ecalEEUncertainty.phiUncertainty->SetParameter(0, 0.087);

  //-- Hcal Uncertainty Functions --------------------------------------//
  //-- From: FastSimulation/Calorimetry/data/HcalResponse.cfi --//
  //-- Hcal Barrel --//
  hcalHBUncertainty.etUncertainty =
      std::make_unique<TF1>("hcalHBEtFunc", "x*sqrt(([0]*[0]/(x*x))+([1]*[1]/x)+([2]*[2]))", 3);
  hcalHBUncertainty.etUncertainty->SetParameter(0, 0.);
  hcalHBUncertainty.etUncertainty->SetParameter(1, 1.22);
  hcalHBUncertainty.etUncertainty->SetParameter(2, 0.05);

  hcalHBUncertainty.phiUncertainty = std::make_unique<TF1>("ecalHBphiFunc", "[0]*x", 1);
  hcalHBUncertainty.phiUncertainty->SetParameter(0, 0.087);

  //-- Hcal Endcap --//
  hcalHEUncertainty.etUncertainty =
      std::make_unique<TF1>("hcalHEEtFunc", "x*sqrt(([0]*[0]/(x*x))+([1]*[1]/x)+([2]*[2]))", 3);
  hcalHEUncertainty.etUncertainty->SetParameter(0, 0.);
  hcalHEUncertainty.etUncertainty->SetParameter(1, 1.3);
  hcalHEUncertainty.etUncertainty->SetParameter(2, 0.05);

  hcalHEUncertainty.phiUncertainty = std::make_unique<TF1>("ecalHEphiFunc", "[0]*x", 1);
  hcalHEUncertainty.phiUncertainty->SetParameter(0, 0.087);

  //-- Hcal Outer --//
  hcalHOUncertainty.etUncertainty =
      std::make_unique<TF1>("hcalHOEtFunc", "x*sqrt(([0]*[0]/(x*x))+([1]*[1]/x)+([2]*[2]))", 3);
  hcalHOUncertainty.etUncertainty->SetParameter(0, 0.);
  hcalHOUncertainty.etUncertainty->SetParameter(1, 1.82);
  hcalHOUncertainty.etUncertainty->SetParameter(2, 0.09);

  hcalHOUncertainty.phiUncertainty = std::make_unique<TF1>("ecalHOphiFunc", "[0]*x", 1);
  hcalHOUncertainty.phiUncertainty->SetParameter(0, 0.087);

  //-- Hcal Forward --//
  hcalHFUncertainty.etUncertainty =
      std::make_unique<TF1>("hcalHFEtFunc", "x*sqrt(([0]*[0]/(x*x))+([1]*[1]/x)+([2]*[2]))", 3);
  hcalHFUncertainty.etUncertainty->SetParameter(0, 0.);
  hcalHFUncertainty.etUncertainty->SetParameter(1, 1.82);
  hcalHFUncertainty.etUncertainty->SetParameter(2, 0.09);

  hcalHFUncertainty.phiUncertainty = std::make_unique<TF1>("ecalHFphiFunc", "[0]*x", 1);
  hcalHFUncertainty.phiUncertainty->SetParameter(0, 0.174);

  //--- Jet Uncertainty Functions --------------------------------------//
  jetUncertainty.etUncertainty = std::make_unique<TF1>("jetEtFunc", "x*sqrt(([0]*[0]/(x*x))+([1]*[1]/x)+([2]*[2]))", 3);
  //-- values from PTDR 1, ch 11.4 --//
  jetUncertainty.etUncertainty->SetParameter(0, jetEtUncertaintyParameter0_);
  jetUncertainty.etUncertainty->SetParameter(1, jetEtUncertaintyParameter1_);
  jetUncertainty.etUncertainty->SetParameter(2, jetEtUncertaintyParameter2_);

  //-- phi value from our own fits --//
  //jetUncertainty.phiUncertainty.reset( new TF1("jetPhiFunc","[0]*x",1) );
  //jetUncertainty.phiUncertainty->SetParameter(0, jetPhiUncertaintyParameter0_);

  //-- phi Functions and values from
  // http://indico.cern.ch/getFile.py/access?contribId=9&sessionId=0&resId=0&materialId=slides&confId=46394
  jetUncertainty.phiUncertainty =
      std::make_unique<TF1>("jetPhiFunc", "x*sqrt(([0]*[0]/(x*x))+([1]*[1]/x)+([2]*[2]))", 3);
  jetUncertainty.phiUncertainty->SetParameter(0, jetPhiUncertaintyParameter0_);
  jetUncertainty.phiUncertainty->SetParameter(1, jetPhiUncertaintyParameter1_);
  jetUncertainty.phiUncertainty->SetParameter(2, jetPhiUncertaintyParameter2_);

  //-- Jet corrections are assumed not to have an error --//
  /*jetCorrUncertainty.etUncertainty.reset( new TF1("jetCorrEtFunc","[0]*x",1) );
  jetCorrUncertainty.etUncertainty->SetParameter(0,0.0);
  jetCorrUncertainty.phiUncertainty.reset( new TF1("jetCorrPhiFunc","[0]*x",1) );
  jetCorrUncertainty.phiUncertainty->SetParameter(0,0.0*(3.14159/180.));*/

  //--- Electron Uncertainty Functions ---------------------------------//
  // completely ambiguious values for electron-like jets...
  // the egamma group keeps track of these here:
  // https://twiki.cern.ch/twiki/bin/view/CMS/EgammaCMSSWVal
  // electron resolution in energy is around 3.4%, measured for 10 < pT < 50 at realistic events with pile-up.

  eleUncertainty.etUncertainty = std::make_unique<TF1>("eleEtFunc", "[0] * x", 1);
  //  eleUncertainty.etUncertainty->SetParameter(0,0.034);
  eleUncertainty.etUncertainty->SetParameter(0, eleEtUncertaintyParameter0_);

  eleUncertainty.phiUncertainty = std::make_unique<TF1>("elePhiFunc", "[0] * x", 1);
  //  eleUncertainty.phiUncertainty->SetParameter(0,1*(3.14159/180.));
  eleUncertainty.phiUncertainty->SetParameter(0, elePhiUncertaintyParameter0_);

  //--- Muon Uncertainty Functions ------------------------------------//
  // and ambiguious values for the muons...

  muonUncertainty.etUncertainty = std::make_unique<TF1>("muonEtFunc", "[0] * x", 1);
  //  muonUncertainty.etUncertainty->SetParameter(0,0.01);
  muonUncertainty.etUncertainty->SetParameter(0, muonEtUncertaintyParameter0_);
  muonUncertainty.phiUncertainty = std::make_unique<TF1>("muonPhiFunc", "[0] * x", 1);
  //  muonUncertainty.phiUncertainty->SetParameter(0,1*(3.14159/180.));
  muonUncertainty.phiUncertainty->SetParameter(0, muonPhiUncertaintyParameter0_);

  //-- Muon calo deposites are assumed not to have an error --//
  /*muonCorrUncertainty.etUncertainty.reset( new TF1("muonCorrEtFunc","[0] * x",1) );
  muonCorrUncertainty.etUncertainty->SetParameter(0,0.0);
  muonCorrUncertainty.phiUncertainty.reset( new TF1("muonCorrPhiFunc","[0] * x",1) );
  muonCorrUncertainty.phiUncertainty->SetParameter(0,0.0*(3.14159/180.)); */
}

Member Data Documentation

edm::InputTag pat::PATMHTProducer::CaloJetAlgorithmTag_

private

Definition at line 132 of file PATMHTProducer.cc.

edm::InputTag pat::PATMHTProducer::CaloTowerTag_

private

Definition at line 137 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

bool pat::PATMHTProducer::controlledUncertainty_

private

Definition at line 67 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

edm::InputTag pat::PATMHTProducer::CorJetAlgorithmTag_

private

Definition at line 133 of file PATMHTProducer.cc.

uncertaintyFunctions pat::PATMHTProducer::ecalEBUncertainty

private

Definition at line 78 of file PATMHTProducer.cc.

uncertaintyFunctions pat::PATMHTProducer::ecalEEUncertainty

private

Definition at line 79 of file PATMHTProducer.cc.

edm::InputTag pat::PATMHTProducer::ElectronTag_

private

Definition at line 136 of file PATMHTProducer.cc.

double pat::PATMHTProducer::eleEtaMax_

private

Definition at line 106 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::eleEtUncertaintyParameter0_

private

Definition at line 126 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::elePhiUncertaintyParameter0_

private

Definition at line 127 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::elePtMin_

private

Definition at line 105 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

edm::EDGetTokenT<edm::View<pat::Electron> > pat::PATMHTProducer::eleToken_

private

Definition at line 59 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

uncertaintyFunctions pat::PATMHTProducer::eleUncertainty

private

Definition at line 87 of file PATMHTProducer.cc.

uncertaintyFunctions pat::PATMHTProducer::hcalHBUncertainty

private

Definition at line 80 of file PATMHTProducer.cc.

uncertaintyFunctions pat::PATMHTProducer::hcalHEUncertainty

private

Definition at line 81 of file PATMHTProducer.cc.

uncertaintyFunctions pat::PATMHTProducer::hcalHFUncertainty

private

Definition at line 83 of file PATMHTProducer.cc.

uncertaintyFunctions pat::PATMHTProducer::hcalHOUncertainty

private

Definition at line 82 of file PATMHTProducer.cc.

std::string pat::PATMHTProducer::JetCorrectionService_

private

Definition at line 134 of file PATMHTProducer.cc.

uncertaintyFunctions pat::PATMHTProducer::jetCorrUncertainty

private

Definition at line 86 of file PATMHTProducer.cc.

double pat::PATMHTProducer::jetEMfracMax_

private

Definition at line 103 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::jetEtaMax_

private

Definition at line 102 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::jetEtUncertaintyParameter0_

private

Definition at line 118 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::jetEtUncertaintyParameter1_

private

Definition at line 119 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::jetEtUncertaintyParameter2_

private

Definition at line 120 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::jetPhiUncertaintyParameter0_

private

Definition at line 122 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::jetPhiUncertaintyParameter1_

private

Definition at line 123 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::jetPhiUncertaintyParameter2_

private

Definition at line 124 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::jetPtMin_

private

Definition at line 101 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

edm::EDGetTokenT<edm::View<pat::Jet> > pat::PATMHTProducer::jetToken_

private

Definition at line 58 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

uncertaintyFunctions pat::PATMHTProducer::jetUncertainty

private

Definition at line 85 of file PATMHTProducer.cc.

std::string pat::PATMHTProducer::metCollectionLabel_

private

Definition at line 138 of file PATMHTProducer.cc.

edm::InputTag pat::PATMHTProducer::mhtLabel_

private

Definition at line 57 of file PATMHTProducer.cc.

double pat::PATMHTProducer::muonChiSqMax_

private

Definition at line 114 of file PATMHTProducer.cc.

uncertaintyFunctions pat::PATMHTProducer::muonCorrUncertainty

private

Definition at line 89 of file PATMHTProducer.cc.

double pat::PATMHTProducer::muonDPtMax_

private

Definition at line 113 of file PATMHTProducer.cc.

double pat::PATMHTProducer::muonEtaMax_

private

Definition at line 109 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::muonEtUncertaintyParameter0_

private

Definition at line 129 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

int pat::PATMHTProducer::muonNHitsMin_

private

Definition at line 112 of file PATMHTProducer.cc.

double pat::PATMHTProducer::muonPhiUncertaintyParameter0_

private

Definition at line 130 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::muonPtMin_

private

Definition at line 108 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

edm::InputTag pat::PATMHTProducer::MuonTag_

private

Definition at line 135 of file PATMHTProducer.cc.

double pat::PATMHTProducer::muonTrackD0Max_

private

Definition at line 110 of file PATMHTProducer.cc.

double pat::PATMHTProducer::muonTrackDzMax_

private

Definition at line 111 of file PATMHTProducer.cc.

uncertaintyFunctions pat::PATMHTProducer::muonUncertainty

private

Definition at line 88 of file PATMHTProducer.cc.

edm::EDGetTokenT<edm::View<pat::Muon> > pat::PATMHTProducer::muoToken_

private

Definition at line 60 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

bool pat::PATMHTProducer::noHF_

private

Definition at line 99 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

edm::EDGetTokenT<edm::View<pat::Photon> > pat::PATMHTProducer::phoToken_

private

Definition at line 62 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

std::vector<metsig::SigInputObj> pat::PATMHTProducer::physobjvector_

private

Definition at line 64 of file PATMHTProducer.cc.

std::set<CaloTowerDetId> pat::PATMHTProducer::s_clusteredTowers

private

Definition at line 97 of file PATMHTProducer.cc.

std::string pat::PATMHTProducer::significanceLabel_

private

Definition at line 139 of file PATMHTProducer.cc.

edm::EDGetTokenT<edm::View<pat::Tau> > pat::PATMHTProducer::tauToken_

private

Definition at line 61 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::towerEtThreshold_

private

Definition at line 145 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::uncertaintyScaleFactor_

private

Definition at line 66 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

bool pat::PATMHTProducer::useCaloTowers_

private

Definition at line 93 of file PATMHTProducer.cc.

bool pat::PATMHTProducer::useElectrons_

private

Definition at line 95 of file PATMHTProducer.cc.

bool pat::PATMHTProducer::useHO_

private

Definition at line 146 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().

bool pat::PATMHTProducer::useJets_

private

Definition at line 94 of file PATMHTProducer.cc.

bool pat::PATMHTProducer::useMuons_

private

Definition at line 96 of file PATMHTProducer.cc.

double pat::PATMHTProducer::verbose_

private

Definition at line 54 of file PATMHTProducer.cc.

Referenced by PATMHTProducer().