pat::PATMHTProducer Class Reference

#include <PATMHTProducer.h>

Inheritance diagram for pat::PATMHTProducer:

Classes
class	uncertaintyFunctions

Public Member Functions
	PATMHTProducer (const edm::ParameterSet &)
	~PATMHTProducer ()
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Member Functions
virtual void	beginJob ()
virtual void	endJob ()
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 &)
virtual 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::EDProducer
typedef EDProducer	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 64 of file PATMHTProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 6 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_, 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");

  produces<pat::MHTCollection>();

}

PATMHTProducer::~PATMHTProducer

(

)

Definition at line 53 of file PATMHTProducer.cc.

                                   {
}

Member Function Documentation

void PATMHTProducer::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 56 of file PATMHTProducer.cc.

                                 {
  setUncertaintyParameters();
}

void PATMHTProducer::endJob ( void  )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 59 of file PATMHTProducer.cc.

                               {
}

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

private

Definition at line 224 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<IdealGeometryRecord>().get(cttopo);
  // const CaloTowerConstituentsMap* caloTowerMap = cttopo.product();

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

  // 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 132 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 317 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  )

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 64 of file PATMHTProducer.cc.

References metsig::significanceAlgo::addObjects(), funct::cos(), gather_cfg::cout, HWWFunctions::getJets(), 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_.size()>0){
    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;


  std::auto_ptr<pat::MHTCollection>  themetsigcoll (new pat::MHTCollection);

  if(physobjvector_.size() >= 1) { // 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( themetsigcoll);


}

void PATMHTProducer::setUncertaintyParameters ( )

private

Definition at line 380 of file PATMHTProducer.cc.

                                                {

  // 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 = new 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 = new TF1("ecalEBphiFunc","[0]*x",1);
  ecalEBUncertainty.phiUncertainty->SetParameter(0,0.0174);

  //-- Ecal Endcap --//
  ecalEEUncertainty.etUncertainty = new 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 = new TF1("ecalEEphiFunc","[0]*x",1);
  ecalEEUncertainty.phiUncertainty->SetParameter(0,0.087);

  //-- Hcal Uncertainty Functions --------------------------------------//
  //-- From: FastSimulation/Calorimetry/data/HcalResponse.cfi --//
  //-- Hcal Barrel --//
  hcalHBUncertainty.etUncertainty = new 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 = new TF1("ecalHBphiFunc","[0]*x",1);
  hcalHBUncertainty.phiUncertainty->SetParameter(0,0.087);

  //-- Hcal Endcap --//
  hcalHEUncertainty.etUncertainty = new 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 = new TF1("ecalHEphiFunc","[0]*x",1);
  hcalHEUncertainty.phiUncertainty->SetParameter(0,0.087);

  //-- Hcal Outer --//
  hcalHOUncertainty.etUncertainty = new 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 = new TF1("ecalHOphiFunc","[0]*x",1);
  hcalHOUncertainty.phiUncertainty->SetParameter(0,0.087);

  //-- Hcal Forward --//
  hcalHFUncertainty.etUncertainty = new 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 = new TF1("ecalHFphiFunc","[0]*x",1);
  hcalHFUncertainty.phiUncertainty->SetParameter(0,0.174);

  //--- Jet Uncertainty Functions --------------------------------------//
  jetUncertainty.etUncertainty = new 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 = 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 = new 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 = new TF1("jetCorrEtFunc","[0]*x",1);
  jetCorrUncertainty.etUncertainty->SetParameter(0,0.0);
  jetCorrUncertainty.phiUncertainty = 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 = new TF1("eleEtFunc","[0] * x",1);
  //  eleUncertainty.etUncertainty->SetParameter(0,0.034);
  eleUncertainty.etUncertainty->SetParameter(0, eleEtUncertaintyParameter0_);


  eleUncertainty.phiUncertainty = new 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 = new TF1("muonEtFunc","[0] * x",1);
  //  muonUncertainty.etUncertainty->SetParameter(0,0.01);
  muonUncertainty.etUncertainty->SetParameter(0, muonEtUncertaintyParameter0_);
  muonUncertainty.phiUncertainty = new 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 = new TF1("muonCorrEtFunc","[0] * x",1);
  muonCorrUncertainty.etUncertainty->SetParameter(0,0.0);
  muonCorrUncertainty.phiUncertainty = 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 161 of file PATMHTProducer.h.

edm::InputTag pat::PATMHTProducer::CaloTowerTag_

private

Definition at line 166 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

bool pat::PATMHTProducer::controlledUncertainty_

private

Definition at line 95 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

edm::InputTag pat::PATMHTProducer::CorJetAlgorithmTag_

private

Definition at line 162 of file PATMHTProducer.h.

uncertaintyFunctions pat::PATMHTProducer::ecalEBUncertainty

private

Definition at line 107 of file PATMHTProducer.h.

uncertaintyFunctions pat::PATMHTProducer::ecalEEUncertainty

private

Definition at line 108 of file PATMHTProducer.h.

edm::InputTag pat::PATMHTProducer::ElectronTag_

private

Definition at line 165 of file PATMHTProducer.h.

double pat::PATMHTProducer::eleEtaMax_

private

Definition at line 135 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::eleEtUncertaintyParameter0_

private

Definition at line 155 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::elePhiUncertaintyParameter0_

private

Definition at line 156 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::elePtMin_

private

Definition at line 134 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

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

private

Definition at line 87 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

uncertaintyFunctions pat::PATMHTProducer::eleUncertainty

private

Definition at line 116 of file PATMHTProducer.h.

uncertaintyFunctions pat::PATMHTProducer::hcalHBUncertainty

private

Definition at line 109 of file PATMHTProducer.h.

uncertaintyFunctions pat::PATMHTProducer::hcalHEUncertainty

private

Definition at line 110 of file PATMHTProducer.h.

uncertaintyFunctions pat::PATMHTProducer::hcalHFUncertainty

private

Definition at line 112 of file PATMHTProducer.h.

uncertaintyFunctions pat::PATMHTProducer::hcalHOUncertainty

private

Definition at line 111 of file PATMHTProducer.h.

std::string pat::PATMHTProducer::JetCorrectionService_

private

Definition at line 163 of file PATMHTProducer.h.

uncertaintyFunctions pat::PATMHTProducer::jetCorrUncertainty

private

Definition at line 115 of file PATMHTProducer.h.

double pat::PATMHTProducer::jetEMfracMax_

private

Definition at line 132 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::jetEtaMax_

private

Definition at line 131 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::jetEtUncertaintyParameter0_

private

Definition at line 147 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::jetEtUncertaintyParameter1_

private

Definition at line 148 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::jetEtUncertaintyParameter2_

private

Definition at line 149 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::jetPhiUncertaintyParameter0_

private

Definition at line 151 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::jetPhiUncertaintyParameter1_

private

Definition at line 152 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::jetPhiUncertaintyParameter2_

private

Definition at line 153 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::jetPtMin_

private

Definition at line 130 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

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

private

Definition at line 86 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

uncertaintyFunctions pat::PATMHTProducer::jetUncertainty

private

Definition at line 114 of file PATMHTProducer.h.

std::string pat::PATMHTProducer::metCollectionLabel_

private

Definition at line 167 of file PATMHTProducer.h.

edm::InputTag pat::PATMHTProducer::mhtLabel_

private

Definition at line 85 of file PATMHTProducer.h.

double pat::PATMHTProducer::muonChiSqMax_

private

Definition at line 143 of file PATMHTProducer.h.

uncertaintyFunctions pat::PATMHTProducer::muonCorrUncertainty

private

Definition at line 118 of file PATMHTProducer.h.

double pat::PATMHTProducer::muonDPtMax_

private

Definition at line 142 of file PATMHTProducer.h.

double pat::PATMHTProducer::muonEtaMax_

private

Definition at line 138 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::muonEtUncertaintyParameter0_

private

Definition at line 158 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

int pat::PATMHTProducer::muonNHitsMin_

private

Definition at line 141 of file PATMHTProducer.h.

double pat::PATMHTProducer::muonPhiUncertaintyParameter0_

private

Definition at line 159 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::muonPtMin_

private

Definition at line 137 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

edm::InputTag pat::PATMHTProducer::MuonTag_

private

Definition at line 164 of file PATMHTProducer.h.

double pat::PATMHTProducer::muonTrackD0Max_

private

Definition at line 139 of file PATMHTProducer.h.

double pat::PATMHTProducer::muonTrackDzMax_

private

Definition at line 140 of file PATMHTProducer.h.

uncertaintyFunctions pat::PATMHTProducer::muonUncertainty

private

Definition at line 117 of file PATMHTProducer.h.

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

private

Definition at line 88 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

bool pat::PATMHTProducer::noHF_

private

Definition at line 128 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

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

private

Definition at line 90 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

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

private

Definition at line 92 of file PATMHTProducer.h.

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

private

Definition at line 126 of file PATMHTProducer.h.

std::string pat::PATMHTProducer::significanceLabel_

private

Definition at line 168 of file PATMHTProducer.h.

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

private

Definition at line 89 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::towerEtThreshold_

private

Definition at line 174 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

double pat::PATMHTProducer::uncertaintyScaleFactor_

private

Definition at line 94 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

bool pat::PATMHTProducer::useCaloTowers_

private

Definition at line 122 of file PATMHTProducer.h.

bool pat::PATMHTProducer::useElectrons_

private

Definition at line 124 of file PATMHTProducer.h.

bool pat::PATMHTProducer::useHO_

private

Definition at line 175 of file PATMHTProducer.h.

Referenced by PATMHTProducer().

bool pat::PATMHTProducer::useJets_

private

Definition at line 123 of file PATMHTProducer.h.

bool pat::PATMHTProducer::useMuons_

private

Definition at line 125 of file PATMHTProducer.h.

double pat::PATMHTProducer::verbose_

private

Definition at line 82 of file PATMHTProducer.h.

Referenced by PATMHTProducer().