EmbeddingKineReweightNtupleProducer Class Reference

#include <EmbeddingKineReweightNtupleProducer.h>

Inheritance diagram for EmbeddingKineReweightNtupleProducer:

Classes
struct	branchEntryType

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &)
void	beginJob ()
	EmbeddingKineReweightNtupleProducer (const edm::ParameterSet &)
	~EmbeddingKineReweightNtupleProducer ()
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	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 (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Types
typedef std::map< std::string, branchEntryType >	branchMap
typedef std::vector < edm::InputTag >	vInputTag

Private Member Functions
void	addBranch_EnPxPyPz (const std::string &)
void	addBranch_MEtResProjections (const std::string &)
void	addBranchF (const std::string &)
void	addBranchI (const std::string &)
void	printBranches (std::ostream &)
void	setValue_EnPxPyPz (const std::string &, const reco::Candidate::LorentzVector &)
void	setValue_MEtResProjections (const std::string &, const reco::Candidate::LorentzVector &, const reco::Candidate::LorentzVector &, const reco::Candidate::LorentzVector &)
void	setValueF (const std::string &, double)
void	setValueI (const std::string &, int)

Private Attributes
branchMap	branches_
std::string	moduleLabel_
TTree *	ntuple_
edm::InputTag	srcGenCaloMEt_
edm::InputTag	srcGenDiTaus_
edm::InputTag	srcGenFilterInfo_
edm::InputTag	srcGenParticles_
edm::InputTag	srcGenPFMEt_
edm::InputTag	srcRecCaloMEt_
edm::InputTag	srcRecJets_
edm::InputTag	srcRecLeg1_
edm::InputTag	srcRecLeg2_
edm::InputTag	srcRecPFMEt_
edm::InputTag	srcSelectedMuons_
vInputTag	srcWeights_

Static Private Attributes
static int	verbosity_ = 0

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
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

Produce Ntuple for computing embeddingKineReweight LUTs

Author

Christian Veelken, LLR

Version

Revision:

1.2

Id:

EmbeddingKineReweightNtupleProducer.h,v 1.2 2013/04/17 07:11:37 veelken Exp

Definition at line 31 of file EmbeddingKineReweightNtupleProducer.h.

Member Typedef Documentation

typedef std::map<std::string, branchEntryType> EmbeddingKineReweightNtupleProducer::branchMap

private

Definition at line 82 of file EmbeddingKineReweightNtupleProducer.h.

typedef std::vector<edm::InputTag> EmbeddingKineReweightNtupleProducer::vInputTag

private

Definition at line 72 of file EmbeddingKineReweightNtupleProducer.h.

Constructor & Destructor Documentation

EmbeddingKineReweightNtupleProducer::EmbeddingKineReweightNtupleProducer

(

const edm::ParameterSet &

cfg

)

Definition at line 28 of file EmbeddingKineReweightNtupleProducer.cc.

References edm::ParameterSet::getParameter(), srcGenCaloMEt_, srcGenDiTaus_, srcGenFilterInfo_, srcGenParticles_, srcGenPFMEt_, srcRecCaloMEt_, srcRecJets_, srcRecLeg1_, srcRecLeg2_, srcRecPFMEt_, srcSelectedMuons_, and srcWeights_.

  : moduleLabel_(cfg.getParameter<std::string>("@module_label")),
    ntuple_(0)
{
  srcGenDiTaus_ = cfg.getParameter<edm::InputTag>("srcGenDiTaus");

  srcGenParticles_ = cfg.getParameter<edm::InputTag>("srcGenParticles"); 
  srcSelectedMuons_ = cfg.getParameter<edm::InputTag>("srcSelectedMuons");

  srcRecLeg1_ = cfg.getParameter<edm::InputTag>("srcRecLeg1");
  srcRecLeg2_ = cfg.getParameter<edm::InputTag>("srcRecLeg2");

  srcRecJets_ = cfg.getParameter<edm::InputTag>("srcRecJets");

  srcGenCaloMEt_ = cfg.getParameter<edm::InputTag>("srcGenCaloMEt");
  srcRecCaloMEt_ = cfg.getParameter<edm::InputTag>("srcRecCaloMEt");
  srcGenPFMEt_ = cfg.getParameter<edm::InputTag>("srcGenPFMEt");
  srcRecPFMEt_ = cfg.getParameter<edm::InputTag>("srcRecPFMEt");

  srcWeights_ = cfg.getParameter<vInputTag>("srcWeights");
  srcGenFilterInfo_ = cfg.getParameter<edm::InputTag>("srcGenFilterInfo");
}

EmbeddingKineReweightNtupleProducer::~EmbeddingKineReweightNtupleProducer

(

)

Definition at line 51 of file EmbeddingKineReweightNtupleProducer.cc.

{
// nothing to be done yet...
}

Member Function Documentation

void EmbeddingKineReweightNtupleProducer::addBranch_EnPxPyPz ( const std::string &  name )

private

Definition at line 361 of file EmbeddingKineReweightNtupleProducer.cc.

References addBranchF(), addBranchI(), bitset_utilities::append(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by beginJob().

{
  addBranchF(std::string(name).append("En"));
  addBranchF(std::string(name).append("Px"));
  addBranchF(std::string(name).append("Py"));
  addBranchF(std::string(name).append("Pz"));
  addBranchF(std::string(name).append("Pt"));
  addBranchF(std::string(name).append("Eta"));
  addBranchF(std::string(name).append("Phi"));
  addBranchF(std::string(name).append("M"));
  addBranchI(std::string(name).append("IsValid"));
}

void EmbeddingKineReweightNtupleProducer::addBranch_MEtResProjections ( const std::string &  name )

private

Definition at line 374 of file EmbeddingKineReweightNtupleProducer.cc.

References addBranchF(), bitset_utilities::append(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by beginJob().

{
  addBranchF(std::string(name).append("ParlZ"));
  addBranchF(std::string(name).append("PerpZ"));
}

void EmbeddingKineReweightNtupleProducer::addBranchF ( const std::string &  name )

private

Definition at line 302 of file EmbeddingKineReweightNtupleProducer.cc.

References assert(), branches_, mergeVDriftHistosByStation::name, ntuple_, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by addBranch_EnPxPyPz(), addBranch_MEtResProjections(), and beginJob().

{
  //std::cout << "<EmbeddingKineReweightNtupleProducer::addBranchF>:" << std::endl;
  //std::cout << " moduleLabel = " << moduleLabel_ << std::endl;
  //std::cout << " branchName = " << name << std::endl;
  assert(branches_.count(name) == 0);
  std::string name_and_format = name + "/F";
  ntuple_->Branch(name.c_str(), &branches_[name].valueF_, name_and_format.c_str());
}

void EmbeddingKineReweightNtupleProducer::addBranchI ( const std::string &  name )

private

Definition at line 312 of file EmbeddingKineReweightNtupleProducer.cc.

References assert(), branches_, mergeVDriftHistosByStation::name, ntuple_, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by addBranch_EnPxPyPz(), and beginJob().

{
  //std::cout << "<EmbeddingKineReweightNtupleProducer::addBranchI>:" << std::endl;
  //std::cout << " moduleLabel = " << moduleLabel_ << std::endl;
  //std::cout << " branchName = " << name << std::endl;
  assert(branches_.count(name) == 0);
  std::string name_and_format = name + "/I";
  ntuple_->Branch(name.c_str(), &branches_[name].valueI_, name_and_format.c_str());
}

void EmbeddingKineReweightNtupleProducer::analyze ( const edm::Event &  evt, const edm::EventSetup &  es  )

virtual

Implements edm::EDAnalyzer.

Definition at line 116 of file EmbeddingKineReweightNtupleProducer.cc.

References Abs(), assert(), reco::Candidate::begin(), branches_, ecal_dqm_sourceclient-live_cfg::cerr, reco::Candidate::charge(), reco::CompositeCandidate::daughter(), reco::Candidate::end(), genParticleCandidates2GenParticles_cfi::genParticles, edm::Event::getByLabel(), getSelMuons(), getTheMuMinus(), getTheMuPlus(), edm::RefToBase< T >::isNonnull(), edm::InputTag::label(), ntuple_, reco::CompositeCandidate::numberOfDaughters(), reco::Candidate::p4(), reco::LeafCandidate::p4(), setValue_EnPxPyPz(), setValue_MEtResProjections(), setValueF(), setValueI(), srcGenCaloMEt_, srcGenDiTaus_, srcGenFilterInfo_, srcGenParticles_, srcGenPFMEt_, srcRecCaloMEt_, srcRecJets_, srcRecLeg1_, srcRecLeg2_, srcRecPFMEt_, srcSelectedMuons_, srcWeights_, AlCaHLTBitMon_QueryRunRegistry::string, and puppiForMET_cff::weight.

{
//--- reset all branches
  for ( branchMap::iterator branch = branches_.begin();
    branch != branches_.end(); ++branch ) {
    branch->second.valueF_ = 0.;
    branch->second.valueI_ = 0;
  }

//--- set branches containing generator level Z->tautau information
  const reco::Candidate* genDiTau_ref = 0;
  if ( srcGenDiTaus_.label() != "" ) {
    edm::Handle<CandidateView> genDiTaus;
    evt.getByLabel(srcGenDiTaus_, genDiTaus);

    for ( CandidateView::const_iterator genDiTau = genDiTaus->begin();
      genDiTau != genDiTaus->end(); ++genDiTau ) {
      const reco::CompositeCandidate* genDiTau_composite = dynamic_cast<const reco::CompositeCandidate*>(&(*genDiTau));
      if ( !(genDiTau_composite && genDiTau_composite->numberOfDaughters() == 2) ) {
    std::cerr << "genDiTau has " << genDiTau_composite->numberOfDaughters() << " daughters --> skipping !!" << std::endl;
    continue;   
      }
      
      const reco::Candidate* genTau1 = genDiTau_composite->daughter(0);
      const reco::Candidate* genTau2 = genDiTau_composite->daughter(1);
      if ( !(genTau1 && genTau2) ) {
    std::cerr << "Failed to find daughters of genDiTau --> skipping !!" << std::endl;
    continue;
      }
      const reco::Candidate* genTauPlus  = 0;
      const reco::Candidate* genTauMinus = 0;
      if ( genTau1->charge() > +0.5 && genTau2->charge() < -0.5 ) {
    genTauPlus  = genTau1;
    genTauMinus = genTau2;
      } else if ( genTau1->charge() < -0.5 && genTau2->charge() > +0.5 ) {
    genTauPlus  = genTau2;
    genTauMinus = genTau1;
      } else {
    std::cerr << "Failed to find daughters of genDiTau --> skipping !!" << std::endl;
    continue;
      }
      
      setValue_EnPxPyPz("genDiTau", genDiTau->p4());
      setValue_EnPxPyPz("genTauPlus", genTauPlus->p4());
      setValue_EnPxPyPz("genTauMinus", genTauMinus->p4());
      genDiTau_ref = &(*genDiTau);
    }
  }

  if ( !genDiTau_ref ) {
    std::cerr << "Failed to find genDiTau --> skipping !!" << std::endl;
    return;
  }

//--- set branches containing generator level Z->mumu information
  if ( srcGenParticles_.label() != "" ) {
    edm::Handle<reco::GenParticleCollection> genParticles;
    evt.getByLabel(srcGenParticles_, genParticles);
    
    for ( reco::GenParticleCollection::const_iterator genParticle = genParticles->begin();
      genParticle != genParticles->end(); ++genParticle ) {
      if ( TMath::Abs(genParticle->pdgId()) == 23 ) {
    std::vector<const reco::GenParticle*> daughters;
    findDaughters(&(*genParticle), daughters, -1);
    const reco::GenParticle* genMuPlus   = 0;
    const reco::GenParticle* genMuMinus  = 0;
    for ( std::vector<const reco::GenParticle*>::const_iterator daughter = daughters.begin();
          daughter != daughters.end(); ++daughter ) {
      if ( (*daughter)->pdgId() == -13 ) genMuPlus   = (*daughter);
      if ( (*daughter)->pdgId() == +13 ) genMuMinus  = (*daughter);
    }
    if ( genMuPlus && genMuMinus ) {
      setValue_EnPxPyPz("genDiMuon", genParticle->p4());
      setValue_EnPxPyPz("genMuonPlus", genMuPlus->p4());
      setValue_EnPxPyPz("genMuonMinus", genMuMinus->p4());
      
//--- set branches containing information on reconstructed Z->mumu decays
      if ( srcSelectedMuons_.label() != "" ) {
        std::vector<reco::CandidateBaseRef> selMuons = getSelMuons(evt, srcSelectedMuons_);
        const reco::CandidateBaseRef recMuPlus  = getTheMuPlus(selMuons);
        const reco::CandidateBaseRef recMuMinus = getTheMuMinus(selMuons);
        if ( recMuPlus.isNonnull() && recMuMinus.isNonnull() ) {
          setValue_EnPxPyPz("recDiMuon", recMuPlus->p4() + recMuMinus->p4());
          setValue_EnPxPyPz("recMuonPlus", recMuPlus->p4());
          setValue_EnPxPyPz("recMuonMinus", recMuMinus->p4());
        }
      }
    }
      }
    }
  }

  edm::Handle<CandidateView> recLeg1;
  evt.getByLabel(srcRecLeg1_, recLeg1);
  if ( recLeg1->size() >= 1 ) {
    setValue_EnPxPyPz("recLeg1", recLeg1->at(0).p4());
  }
  edm::Handle<CandidateView> recLeg2;
  evt.getByLabel(srcRecLeg2_, recLeg2);
  if ( srcRecLeg2_ == srcRecLeg1_ && recLeg2->size() >= 2 ) {
    setValue_EnPxPyPz("recLeg2", recLeg1->at(1).p4());
  } else if ( recLeg2->size() >= 1 ) {
    setValue_EnPxPyPz("recLeg2", recLeg2->at(0).p4());
  }

  edm::Handle<pat::JetCollection> recJets;
  evt.getByLabel(srcRecJets_, recJets);
  int numJetsRawPtGt20  = 0;
  int numJetsCorrPtGt20 = 0;
  int numJetsRawPtGt30  = 0;
  int numJetsCorrPtGt30 = 0;
  for ( pat::JetCollection::const_iterator recJet = recJets->begin();
    recJet != recJets->end(); ++recJet ) {
    
    reco::Candidate::LorentzVector rawJetP4 = recJet->correctedP4("Uncorrected");
    double rawJetPt = rawJetP4.pt();
    double rawJetAbsEta = TMath::Abs(rawJetP4.eta());
    if ( rawJetAbsEta < 4.5 ) {
      if ( rawJetPt > 20. ) ++numJetsRawPtGt20;
      if ( rawJetPt > 30. ) ++numJetsRawPtGt30;
    }
    
    reco::Candidate::LorentzVector corrJetP4 = recJet->p4();
    double corrJetPt = corrJetP4.pt();
    double corrJetAbsEta = TMath::Abs(corrJetP4.eta());
    if ( corrJetAbsEta < 4.5 ) {
      if ( corrJetPt > 20. ) ++numJetsCorrPtGt20;
      if ( corrJetPt > 30. ) ++numJetsCorrPtGt30;
    }
  }
  setValueI("numJetsRawPtGt20", numJetsRawPtGt20);
  setValueI("numJetsCorrPtGt20", numJetsCorrPtGt20);
  setValueI("numJetsRawPtGt30", numJetsRawPtGt30);
  setValueI("numJetsCorrPtGt30", numJetsCorrPtGt30);

  if ( genDiTau_ref ) {
    edm::Handle<METView> genCaloMET;
    evt.getByLabel(srcGenCaloMEt_, genCaloMET);
    const reco::Candidate::LorentzVector& genCaloMEtP4 = genCaloMET->front().p4();
    setValue_EnPxPyPz("genCaloMEt", genCaloMEtP4);
    edm::Handle<METView> recCaloMET;
    evt.getByLabel(srcRecCaloMEt_, recCaloMET);
    const reco::Candidate::LorentzVector& recCaloMEtP4 = recCaloMET->front().p4();
    //std::cout << "<EmbeddingKineReweightNtupleProducer>:" << std::endl;
    //std::cout << " recCaloMEt(" << srcRecCaloMEt_.label() << "): Pt = " << recCaloMEtP4.pt() << ", phi = " << recCaloMEtP4.phi() 
    //          << " (Px = " << recCaloMEtP4.px() << ", Py = " << recCaloMEtP4.py() << ")" << std::endl;
    setValue_EnPxPyPz("recCaloMEt", recCaloMEtP4);
    setValue_MEtResProjections("recMinusGenCaloMEt", genCaloMEtP4, recCaloMEtP4, genDiTau_ref->p4());

    edm::Handle<METView> genPFMET;
    evt.getByLabel(srcGenPFMEt_, genPFMET);
    const reco::Candidate::LorentzVector& genPFMEtP4 = genPFMET->front().p4();
    setValue_EnPxPyPz("genPFMEt", genPFMEtP4);
    edm::Handle<METView> recPFMET;
    evt.getByLabel(srcRecPFMEt_, recPFMET);
    const reco::Candidate::LorentzVector& recPFMEtP4 = recPFMET->front().p4();
    setValue_EnPxPyPz("recPFMEt", recPFMEtP4);
    setValue_MEtResProjections("recMinusGenPFMEt", genPFMEtP4, recPFMEtP4, genDiTau_ref->p4());
  }

//--- set branches containing event weight information 
  for ( vInputTag::const_iterator srcWeight = srcWeights_.begin();
    srcWeight != srcWeights_.end(); ++srcWeight ) {
    edm::Handle<double> weight;
    evt.getByLabel(*srcWeight, weight);
    std::string branchName = Form("%s_%s", srcWeight->label().data(), srcWeight->instance().data());
    setValueF(branchName.data(), *weight);
  }
  if ( srcGenFilterInfo_.label() != "" ) {
    edm::Handle<GenFilterInfo> genFilterInfo;
    evt.getByLabel(srcGenFilterInfo_, genFilterInfo);
    double genFilterInfo_weight = 1.;
    bool genFilterInfo_isValid = false;
    if ( genFilterInfo->numEventsTried() > 0 ) {
      genFilterInfo_weight = genFilterInfo->filterEfficiency();
      genFilterInfo_isValid = true;
    }
    setValueF("genFilterInfo", genFilterInfo_weight);
    setValueI("genFilterInfoIsValid", genFilterInfo_isValid);
  }

//--- finally, fill all computed quantities into TTree
  assert(ntuple_);
  ntuple_->Fill();
}

void EmbeddingKineReweightNtupleProducer::beginJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 56 of file EmbeddingKineReweightNtupleProducer.cc.

References addBranch_EnPxPyPz(), addBranch_MEtResProjections(), addBranchF(), addBranchI(), edm::InputTag::label(), TFileService::make(), ntuple_, srcGenFilterInfo_, srcWeights_, and AlCaHLTBitMon_QueryRunRegistry::string.

{
//--- create TTree
  edm::Service<TFileService> fs;
  ntuple_ = fs->make<TTree>("embeddingKineReweightNtuple", "embeddingKineReweightNtuple");

//--- add branches storing quantities for computing embeddingKineReweight LUTs
//     1.) for embeddingKineReweightGENtoEmbedded:
  addBranch_EnPxPyPz("genDiTau");
  addBranch_EnPxPyPz("genTau1");
  addBranch_EnPxPyPz("genTau2");  
//     2.) for embeddingKineReweightGENtoREC:
  addBranch_EnPxPyPz("genDiMuon");
  addBranch_EnPxPyPz("genMuonPlus");
  addBranch_EnPxPyPz("genMuonMinus");
  addBranch_EnPxPyPz("recDiMuon");
  addBranch_EnPxPyPz("recMuonPlus");
  addBranch_EnPxPyPz("recMuonMinus");

  addBranch_EnPxPyPz("recLeg1");
  addBranch_EnPxPyPz("recLeg2");

  addBranchI("numJetsRawPtGt20");
  addBranchI("numJetsCorrPtGt20");
  addBranchI("numJetsRawPtGt30");
  addBranchI("numJetsCorrPtGt30");

  addBranch_EnPxPyPz("genCaloMEt");
  addBranch_EnPxPyPz("recCaloMEt");
  addBranch_MEtResProjections("recMinusGenCaloMEt");
  addBranch_EnPxPyPz("genPFMEt");
  addBranch_EnPxPyPz("recPFMEt");
  addBranch_MEtResProjections("recMinusGenPFMEt");

  for ( vInputTag::const_iterator srcWeight = srcWeights_.begin();
    srcWeight != srcWeights_.end(); ++srcWeight ) {
    std::string branchName = Form("%s_%s", srcWeight->label().data(), srcWeight->instance().data());
    addBranchF(branchName.data());
  }
  if ( srcGenFilterInfo_.label() != "" ) {
    addBranchF("genFilterInfo");
    addBranchI("genFilterInfoIsValid");
  }
}

void EmbeddingKineReweightNtupleProducer::printBranches ( std::ostream &  stream )

private

Definition at line 322 of file EmbeddingKineReweightNtupleProducer.cc.

References branches_, and moduleLabel_.

{
  stream << "<EmbeddingKineReweightNtupleProducer::printBranches>:" << std::endl;
  stream << " registered branches for module = " << moduleLabel_ << std::endl;
  for ( branchMap::const_iterator branch = branches_.begin();
        branch != branches_.end(); ++branch ) {
    stream << " " << branch->first << std::endl;
  }
  stream << std::endl;
}

void EmbeddingKineReweightNtupleProducer::setValue_EnPxPyPz ( const std::string &  name, const reco::Candidate::LorentzVector &  p4  )

private

Definition at line 384 of file EmbeddingKineReweightNtupleProducer.cc.

References bitset_utilities::append(), setValueF(), setValueI(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by analyze().

{
  setValueF(std::string(name).append("En"), p4.E());
  setValueF(std::string(name).append("Px"), p4.px());
  setValueF(std::string(name).append("Py"), p4.py());
  setValueF(std::string(name).append("Pz"), p4.pz());
  setValueF(std::string(name).append("Pt"), p4.pt());
  setValueF(std::string(name).append("Eta"), p4.eta());
  setValueF(std::string(name).append("Phi"), p4.phi());
  setValueF(std::string(name).append("M"), p4.M());
  setValueI(std::string(name).append("IsValid"), true);
}

void EmbeddingKineReweightNtupleProducer::setValue_MEtResProjections ( const std::string &  name, const reco::Candidate::LorentzVector &  genMEtP4, const reco::Candidate::LorentzVector &  recMEtP4, const reco::Candidate::LorentzVector &  zP4  )

private

Definition at line 397 of file EmbeddingKineReweightNtupleProducer.cc.

References bitset_utilities::append(), setValueF(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by analyze().

{
  if ( zP4.pt() > 0. ) {
    double qX = zP4.px();
    double qY = zP4.py();
    double qT = TMath::Sqrt(qX*qX + qY*qY);
    double dX = recMEtP4.px() - genMEtP4.px();
    double dY = recMEtP4.py() - genMEtP4.py();
    double dParl = (dX*qX + dY*qY)/qT;
    double dPerp = (dX*qY - dY*qX)/qT;
    setValueF(std::string(name).append("ParlZ"), dParl);
    setValueF(std::string(name).append("PerpZ"), dPerp);
  }
}

void EmbeddingKineReweightNtupleProducer::setValueF ( const std::string &  name, double  value  )

private

Definition at line 333 of file EmbeddingKineReweightNtupleProducer.cc.

References branches_, gather_cfg::cout, Exception, relativeConstraints::value, and verbosity_.

Referenced by analyze(), setValue_EnPxPyPz(), and setValue_MEtResProjections().

{
  if ( verbosity_ ) std::cout << "branch = " << name << ": value = " << value << std::endl;
  branchMap::iterator branch = branches_.find(name);
  if ( branch != branches_.end() ) {
    branch->second.valueF_ = value;
  } else {
    throw cms::Exception("InvalidParameter") 
      << "No branch with name = " << name << " defined !!\n";
  }
}

void EmbeddingKineReweightNtupleProducer::setValueI ( const std::string &  name, int  value  )

private

Definition at line 345 of file EmbeddingKineReweightNtupleProducer.cc.

References branches_, gather_cfg::cout, Exception, relativeConstraints::value, and verbosity_.

Referenced by analyze(), and setValue_EnPxPyPz().

{
  if ( verbosity_ ) std::cout << "branch = " << name << ": value = " << value << std::endl;
  branchMap::iterator branch = branches_.find(name);
  if ( branch != branches_.end() ) {
    branch->second.valueI_ = value;
  } else {
    throw cms::Exception("InvalidParameter") 
      << "No branch with name = " << name << " defined !!\n";
  }
}

Member Data Documentation

branchMap EmbeddingKineReweightNtupleProducer::branches_

private

Definition at line 83 of file EmbeddingKineReweightNtupleProducer.h.

Referenced by addBranchF(), addBranchI(), analyze(), printBranches(), setValueF(), and setValueI().

std::string EmbeddingKineReweightNtupleProducer::moduleLabel_

private

Definition at line 55 of file EmbeddingKineReweightNtupleProducer.h.

Referenced by Mixins._TypedParameterizable::insertInto(), Modules.ESSource::nameInProcessDesc_(), Modules.ESProducer::nameInProcessDesc_(), and printBranches().

TTree* EmbeddingKineReweightNtupleProducer::ntuple_

private

Definition at line 85 of file EmbeddingKineReweightNtupleProducer.h.

Referenced by addBranchF(), addBranchI(), analyze(), and beginJob().

edm::InputTag EmbeddingKineReweightNtupleProducer::srcGenCaloMEt_

private

Definition at line 67 of file EmbeddingKineReweightNtupleProducer.h.

Referenced by analyze(), and EmbeddingKineReweightNtupleProducer().

edm::InputTag EmbeddingKineReweightNtupleProducer::srcGenDiTaus_

private

Definition at line 57 of file EmbeddingKineReweightNtupleProducer.h.

Referenced by analyze(), and EmbeddingKineReweightNtupleProducer().

edm::InputTag EmbeddingKineReweightNtupleProducer::srcGenFilterInfo_

private

Definition at line 74 of file EmbeddingKineReweightNtupleProducer.h.

Referenced by analyze(), beginJob(), and EmbeddingKineReweightNtupleProducer().

edm::InputTag EmbeddingKineReweightNtupleProducer::srcGenParticles_

private

Definition at line 59 of file EmbeddingKineReweightNtupleProducer.h.

Referenced by analyze(), and EmbeddingKineReweightNtupleProducer().

edm::InputTag EmbeddingKineReweightNtupleProducer::srcGenPFMEt_

private

Definition at line 69 of file EmbeddingKineReweightNtupleProducer.h.

Referenced by analyze(), and EmbeddingKineReweightNtupleProducer().

edm::InputTag EmbeddingKineReweightNtupleProducer::srcRecCaloMEt_

private

Definition at line 68 of file EmbeddingKineReweightNtupleProducer.h.

Referenced by analyze(), and EmbeddingKineReweightNtupleProducer().

edm::InputTag EmbeddingKineReweightNtupleProducer::srcRecJets_

private

Definition at line 65 of file EmbeddingKineReweightNtupleProducer.h.

Referenced by analyze(), and EmbeddingKineReweightNtupleProducer().

edm::InputTag EmbeddingKineReweightNtupleProducer::srcRecLeg1_

private

Definition at line 62 of file EmbeddingKineReweightNtupleProducer.h.

Referenced by analyze(), and EmbeddingKineReweightNtupleProducer().

edm::InputTag EmbeddingKineReweightNtupleProducer::srcRecLeg2_

private

Definition at line 63 of file EmbeddingKineReweightNtupleProducer.h.

Referenced by analyze(), and EmbeddingKineReweightNtupleProducer().

edm::InputTag EmbeddingKineReweightNtupleProducer::srcRecPFMEt_

private

Definition at line 70 of file EmbeddingKineReweightNtupleProducer.h.

Referenced by analyze(), and EmbeddingKineReweightNtupleProducer().

edm::InputTag EmbeddingKineReweightNtupleProducer::srcSelectedMuons_

private

Definition at line 60 of file EmbeddingKineReweightNtupleProducer.h.

Referenced by analyze(), and EmbeddingKineReweightNtupleProducer().

vInputTag EmbeddingKineReweightNtupleProducer::srcWeights_

private

Definition at line 73 of file EmbeddingKineReweightNtupleProducer.h.

Referenced by analyze(), beginJob(), and EmbeddingKineReweightNtupleProducer().

int EmbeddingKineReweightNtupleProducer::verbosity_ = 0

staticprivate

Definition at line 87 of file EmbeddingKineReweightNtupleProducer.h.

Referenced by setValueF(), and setValueI().