PFJetBenchmark Class Reference

#include <PFJetBenchmark.h>

Public Member Functions
void	gettrue (const reco::GenJet *truth, double &true_ChargedHadEnergy, double &true_NeutralHadEnergy, double &true_NeutralEmEnergy)
	PFJetBenchmark ()
void	printGenJet (const reco::GenJet *)
void	printPFJet (const reco::PFJet *)
void	process (const reco::PFJetCollection &, const reco::GenJetCollection &)
double	resChargedHadEnergyMax () const
double	resNeutralEmEnergyMax () const
double	resNeutralHadEnergyMax () const
double	resPtMax () const
void	setup (std::string Filename, bool debug, bool plotAgainstReco=0, bool onlyTwoJets=1, double deltaRMax=0.1, std::string benchmarkLabel_="ParticleFlow", double recPt=-1, double maxEta=-1, DQMStore *dbe_store=NULL)
void	write ()
virtual	~PFJetBenchmark ()

Protected Attributes
PFBenchmarkAlgo *	algo_
DQMStore *	dbe_
bool	debug_
double	deltaRMax_
unsigned int	entry_
double	maxEta_cut
bool	onlyTwoJets_
bool	plotAgainstReco_
double	recPt_cut
double	resChargedHadEnergyMax_
double	resNeutralEmEnergyMax_
double	resNeutralHadEnergyMax_
double	resPtMax_

Private Attributes
TFile *	file_
TH2F *	hBDEtavsPt
TH2F *	hBDPhivsPt
TH1F *	hBNCH
TH2F *	hBNCH0vsPt
TH2F *	hBNCH1vsPt
TH2F *	hBNCH2vsPt
TH2F *	hBNCH3vsPt
TH2F *	hBNCH4vsPt
TH2F *	hBNCH5vsPt
TH2F *	hBNCH6vsPt
TH2F *	hBNCH7vsPt
TH2F *	hBNCHvsPt
TH1F *	hBRCHE
TH2F *	hBRCHEvsPt
TH2F *	hBRHCALvsP
TH2F *	hBRHONLvsP
TH1F *	hBRNEE
TH2F *	hBRNEEvsPt
TH1F *	hBRneut
TH2F *	hBRNEUTvsP
TH2F *	hBRneutvsPt
TH1F *	hBRNHE
TH2F *	hBRNHEvsPt
TH2F *	hBRNONLvsP
TH1F *	hBRPt
TH1F *	hBRPt100_150
TH1F *	hBRPt1250_2000
TH1F *	hBRPt150_200
TH1F *	hBRPt2000_5000
TH1F *	hBRPt200_250
TH1F *	hBRPt20_40
TH1F *	hBRPt250_300
TH1F *	hBRPt300_400
TH1F *	hBRPt400_500
TH1F *	hBRPt40_60
TH1F *	hBRPt500_750
TH1F *	hBRPt60_80
TH1F *	hBRPt750_1250
TH1F *	hBRPt80_100
TH2F *	hBRPtvsPt
TH2F *	hDEtavsEta
TH2F *	hDPhivsEta
TH2F *	hEDEtavsPt
TH2F *	hEDPhivsPt
TH1F *	hENCH
TH2F *	hENCH0vsPt
TH2F *	hENCH1vsPt
TH2F *	hENCH2vsPt
TH2F *	hENCH3vsPt
TH2F *	hENCH4vsPt
TH2F *	hENCH5vsPt
TH2F *	hENCH6vsPt
TH2F *	hENCH7vsPt
TH2F *	hENCHvsPt
TH1F *	hERCHE
TH2F *	hERCHEvsPt
TH2F *	hERHCALvsP
TH2F *	hERHONLvsP
TH1F *	hERNEE
TH2F *	hERNEEvsPt
TH1F *	hERneut
TH2F *	hERNEUTvsP
TH2F *	hERneutvsPt
TH1F *	hERNHE
TH2F *	hERNHEvsPt
TH2F *	hERNONLvsP
TH1F *	hERPt
TH1F *	hERPt100_150
TH1F *	hERPt1250_2000
TH1F *	hERPt150_200
TH1F *	hERPt2000_5000
TH1F *	hERPt200_250
TH1F *	hERPt20_40
TH1F *	hERPt250_300
TH1F *	hERPt300_400
TH1F *	hERPt400_500
TH1F *	hERPt40_60
TH1F *	hERPt500_750
TH1F *	hERPt60_80
TH1F *	hERPt750_1250
TH1F *	hERPt80_100
TH2F *	hERPtvsPt
TH2F *	hFDEtavsPt
TH2F *	hFDPhivsPt
TH1F *	hFNCH
TH2F *	hFNCH0vsPt
TH2F *	hFNCH1vsPt
TH2F *	hFNCH2vsPt
TH2F *	hFNCH3vsPt
TH2F *	hFNCH4vsPt
TH2F *	hFNCH5vsPt
TH2F *	hFNCH6vsPt
TH2F *	hFNCH7vsPt
TH2F *	hFNCHvsPt
TH1F *	hFRCHE
TH2F *	hFRCHEvsPt
TH2F *	hFRHCALvsP
TH2F *	hFRHONLvsP
TH1F *	hFRNEE
TH2F *	hFRNEEvsPt
TH1F *	hFRneut
TH2F *	hFRNEUTvsP
TH2F *	hFRneutvsPt
TH1F *	hFRNHE
TH2F *	hFRNHEvsPt
TH2F *	hFRNONLvsP
TH1F *	hFRPt
TH1F *	hFRPt100_150
TH1F *	hFRPt1250_2000
TH1F *	hFRPt150_200
TH1F *	hFRPt2000_5000
TH1F *	hFRPt200_250
TH1F *	hFRPt20_40
TH1F *	hFRPt250_300
TH1F *	hFRPt300_400
TH1F *	hFRPt400_500
TH1F *	hFRPt40_60
TH1F *	hFRPt500_750
TH1F *	hFRPt60_80
TH1F *	hFRPt750_1250
TH1F *	hFRPt80_100
TH2F *	hFRPtvsPt
TH1F *	hjetsEta
TH1F *	hjetsPt
TH2F *	hNCH0vsEta
TH2F *	hNCH1vsEta
TH2F *	hNCH2vsEta
TH2F *	hNCH3vsEta
TH2F *	hNCH4vsEta
TH2F *	hNCH5vsEta
TH2F *	hNCH6vsEta
TH2F *	hNCH7vsEta
TH2F *	hNCHvsEta
TH1F *	hNjets
TH2F *	hRCHEvsEta
TH2F *	hRHCALvsEta
TH2F *	hRHONLvsEta
TH2F *	hRNeutvsEta
TH2F *	hRNEUTvsEta
TH2F *	hRNONLvsEta
TH2F *	hRPtvsEta
std::string	outputFile_

Detailed Description

Definition at line 38 of file PFJetBenchmark.h.

Constructor & Destructor Documentation

PFJetBenchmark::PFJetBenchmark

(

)

Definition at line 40 of file PFJetBenchmark.cc.

: file_(0), entry_(0) {}

PFJetBenchmark::~PFJetBenchmark ( )

virtual

Definition at line 42 of file PFJetBenchmark.cc.

References file_.

                                {
  if(file_) file_->Close();
}

Member Function Documentation

void PFJetBenchmark::gettrue	(	const reco::GenJet *	truth,
		double &	true_ChargedHadEnergy,
		double &	true_NeutralHadEnergy,
		double &	true_NeutralEmEnergy
	)

Definition at line 677 of file PFJetBenchmark.cc.

References abs, alignCSCRings::e, reco::LeafCandidate::energy(), reco::GenJet::getGenConstituents(), i, and reco::LeafCandidate::pdgId().

Referenced by process().

                                                                              {
  std::vector <const GenParticle*> mcparts = truth->getGenConstituents ();
  true_NeutralEmEnergy = 0.;
  true_ChargedHadEnergy = 0.;
  true_NeutralHadEnergy = 0.;
  // for each MC particle in turn  
  for (unsigned i = 0; i < mcparts.size (); i++) {
    const GenParticle* mcpart = mcparts[i];
    int PDG = std::abs( mcpart->pdgId());
    double e = mcpart->energy(); 
    switch(PDG){  // start PDG switch
    case 22: // photon
      true_NeutralEmEnergy += e;
      break;
    case 211: // pi
    case 321: // K
    case 2212: // p
    case 11: //electrons (until recognised)
      true_ChargedHadEnergy += e;
      break;
    case 310: // K_S0
    case 130: // K_L0
    case 3122: // Lambda0
    case 2112: // n0
      true_NeutralHadEnergy += e;
    default:
      break;
    }  // end PDG switch        
  }  // end loop on constituents.
}

void PFJetBenchmark::printGenJet

(

const reco::GenJet *

truth

)

Definition at line 730 of file PFJetBenchmark.cc.

References gather_cfg::cout, reco::LeafCandidate::eta(), reco::GenJet::getGenConstituents(), i, reco::LeafCandidate::p(), reco::LeafCandidate::pdgId(), reco::LeafCandidate::phi(), reco::LeafCandidate::pt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), and reco::LeafCandidate::pz().

Referenced by process().

                                                        {
  std::vector <const GenParticle*> mcparts = truth->getGenConstituents ();
  cout << "GenJet p/px/py/pz/pt: " << truth->p() << '/' << truth->px () 
       << '/' << truth->py() << '/' << truth->pz() << '/' << truth->pt() << endl
       << "    eta/phi: " << truth->eta () << '/' << truth->phi () << endl
       << "    # of constituents: " << mcparts.size() << endl;
  cout << "    constituents:" << endl;
  for (unsigned i = 0; i < mcparts.size (); i++) {
    const GenParticle* mcpart = mcparts[i];
    cout << "      #" << i << "  PDG code:" << mcpart->pdgId() 
     << ", p/pt/eta/phi: " << mcpart->p() << '/' << mcpart->pt() 
     << '/' << mcpart->eta() << '/' << mcpart->phi() << endl;   
  }    
}

void PFJetBenchmark::printPFJet

(

const reco::PFJet *

pfj

)

Definition at line 709 of file PFJetBenchmark.cc.

References reco::PFJet::chargedEmEnergy(), reco::PFJet::chargedHadronEnergy(), reco::PFJet::chargedMuEnergy(), reco::PFJet::chargedMultiplicity(), gather_cfg::cout, reco::LeafCandidate::eta(), reco::PFJet::neutralEmEnergy(), reco::PFJet::neutralHadronEnergy(), reco::PFJet::neutralMultiplicity(), reco::LeafCandidate::p(), reco::LeafCandidate::phi(), reco::PFJet::print(), reco::LeafCandidate::pt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), and reco::LeafCandidate::pz().

Referenced by process().

                                                   {
  cout<<setiosflags(ios::right);
  cout<<setiosflags(ios::fixed);
  cout<<setprecision(3);

  cout << "PFJet  p/px/py/pz/pt: " << pfj->p() << "/" << pfj->px () 
       << "/" << pfj->py() << "/" << pfj->pz() << "/" << pfj->pt() << endl
       << "    eta/phi: " << pfj->eta () << "/" << pfj->phi () << endl          
       << "    PFJet specific:" << std::endl
       << "      charged/neutral hadrons energy: " << pfj->chargedHadronEnergy () << '/' << pfj->neutralHadronEnergy () << endl
       << "      charged/neutral em energy: " << pfj->chargedEmEnergy () << '/' << pfj->neutralEmEnergy () << endl
       << "      charged muon energy: " << pfj->chargedMuEnergy () << '/' << endl
       << "      charged/neutral multiplicity: " << pfj->chargedMultiplicity () << '/' << pfj->neutralMultiplicity () << endl;
  
  // And print the constituents
  std::cout << pfj->print() << std::endl;

  cout<<resetiosflags(ios::right|ios::fixed);
}

void PFJetBenchmark::process	(	const reco::PFJetCollection &	pfJets,
		const reco::GenJetCollection &	genJets
	)

Definition at line 230 of file PFJetBenchmark.cc.

References abs, algo_, sipixelsummary::Barrel, reco::PFJet::chargedHadronEnergy(), reco::PFJet::chargedMultiplicity(), gather_cfg::cout, debug_, deltaR(), PFBenchmarkAlgo::deltaR(), deltaRMax_, entry_, reco::LeafCandidate::eta(), reco::PFJet::getPFConstituents(), gettrue(), hBDEtavsPt, hBDPhivsPt, hBNCH, hBNCH0vsPt, hBNCH1vsPt, hBNCH2vsPt, hBNCH3vsPt, hBNCH4vsPt, hBNCH5vsPt, hBNCH6vsPt, hBNCH7vsPt, hBNCHvsPt, hBRCHE, hBRCHEvsPt, hBRHCALvsP, hBRHONLvsP, hBRNEE, hBRNEEvsPt, hBRneut, hBRNEUTvsP, hBRneutvsPt, hBRNHE, hBRNHEvsPt, hBRNONLvsP, hBRPt, hBRPt100_150, hBRPt1250_2000, hBRPt150_200, hBRPt2000_5000, hBRPt200_250, hBRPt20_40, hBRPt250_300, hBRPt300_400, hBRPt400_500, hBRPt40_60, hBRPt500_750, hBRPt60_80, hBRPt750_1250, hBRPt80_100, hBRPtvsPt, hDEtavsEta, hDPhivsEta, hEDEtavsPt, hEDPhivsPt, hENCH, hENCH0vsPt, hENCH1vsPt, hENCH2vsPt, hENCH3vsPt, hENCH4vsPt, hENCH5vsPt, hENCH6vsPt, hENCH7vsPt, hENCHvsPt, hERCHE, hERCHEvsPt, hERHCALvsP, hERHONLvsP, hERNEE, hERNEEvsPt, hERneut, hERNEUTvsP, hERneutvsPt, hERNHE, hERNHEvsPt, hERNONLvsP, hERPt, hERPt100_150, hERPt1250_2000, hERPt150_200, hERPt2000_5000, hERPt200_250, hERPt20_40, hERPt250_300, hERPt300_400, hERPt400_500, hERPt40_60, hERPt500_750, hERPt60_80, hERPt750_1250, hERPt80_100, hERPtvsPt, hFDEtavsPt, hFDPhivsPt, hFRCHE, hFRCHEvsPt, hFRHCALvsP, hFRHONLvsP, hFRNEE, hFRNEEvsPt, hFRneut, hFRNEUTvsP, hFRneutvsPt, hFRNHE, hFRNHEvsPt, hFRNONLvsP, hFRPt, hFRPt100_150, hFRPt1250_2000, hFRPt150_200, hFRPt2000_5000, hFRPt200_250, hFRPt20_40, hFRPt250_300, hFRPt300_400, hFRPt400_500, hFRPt40_60, hFRPt500_750, hFRPt60_80, hFRPt750_1250, hFRPt80_100, hFRPtvsPt, hjetsEta, hjetsPt, hNCH0vsEta, hNCH1vsEta, hNCH2vsEta, hNCH3vsEta, hNCH4vsEta, hNCH5vsEta, hNCH6vsEta, hNCH7vsEta, hNCHvsEta, hNjets, hRCHEvsEta, hRHCALvsEta, hRHONLvsEta, hRNeutvsEta, hRNEUTvsEta, hRNONLvsEta, hRPtvsEta, i, edm::Ref< C, T, F >::isNull(), j, PFBenchmarkAlgo::matchByDeltaR(), maxEta_cut, reco::PFJet::neutralEmEnergy(), reco::PFJet::neutralHadronEnergy(), onlyTwoJets_, reco::LeafCandidate::p(), reco::LeafCandidate::phi(), plotAgainstReco_, printGenJet(), printPFJet(), reco::LeafCandidate::pt(), recPt_cut, resChargedHadEnergyMax_, resNeutralEmEnergyMax_, resNeutralHadEnergyMax_, and resPtMax_.

Referenced by PFJetBenchmarkAnalyzer::analyze(), ConfigBuilder.ConfigBuilder.PrintAllModules::leave(), and PFRootEventManager::processEntry().

                                                                                                   {
  // loop over reco  pf  jets
  resPtMax_ = 0.;
  resChargedHadEnergyMax_ = 0.;
  resNeutralHadEnergyMax_ = 0.;
  resNeutralEmEnergyMax_ = 0.; 
  int NPFJets = 0;
    
  for(unsigned i=0; i<pfJets.size(); i++) {   

    // Count the number of jets with a larger energy
    unsigned highJets = 0;
    for(unsigned j=0; j<pfJets.size(); j++) { 
      if ( j != i && pfJets[j].pt() > pfJets[i].pt() ) highJets++;
    }
    if ( onlyTwoJets_ && highJets > 1 ) continue;
        
        
    const reco::PFJet& pfj = pfJets[i];
    double rec_pt = pfj.pt();
    double rec_eta = pfj.eta();
    double rec_phi = pfj.phi();

    // skip PFjets with pt < recPt_cut GeV
    if (rec_pt<recPt_cut and recPt_cut != -1.) continue;
    // skip PFjets with eta > maxEta_cut
    if (fabs(rec_eta)>maxEta_cut and maxEta_cut != -1.) continue;

    NPFJets++;
        
    // fill inclusive PFjet distribution pt > 20 GeV
    hNjets->Fill(NPFJets);
    hjetsPt->Fill(rec_pt);
    hjetsEta->Fill(rec_eta);

    // separate Barrel PFJets from Endcap PFJets
    bool Barrel = false;
    bool Endcap = false;
    bool Forward = false;
    if (std::abs(rec_eta) < 1.4 ) Barrel = true;
    if (std::abs (rec_eta) > 1.6 && std::abs (rec_eta) < 2.4 ) Endcap = true;
    if (std::abs (rec_eta) > 2.5 && std::abs (rec_eta) < 2.9 ) Forward = true;
    if (std::abs (rec_eta) > 3.1 && std::abs (rec_eta) < 4.7 ) Forward = true;

    // do only barrel for now
    //  if(!Barrel) continue;

    // look for the closets gen Jet : truth
    const GenJet *truth = algo_->matchByDeltaR(&pfj,&genJets);
    if(!truth) continue;   
    double deltaR = algo_->deltaR(&pfj, truth);
    // check deltaR is small enough
    if(deltaR < deltaRMax_ || (abs(rec_eta)>2.5 && deltaR < 0.2) || deltaRMax_ == -1.0 ) {//start case deltaR < deltaRMax

      // generate histograms comparing the reco and truth candidate (truth = closest in delta-R) 
      // get the quantities to place on the denominator and/or divide by
      double pt_denom;
      double true_E = truth->p();
      double true_pt = truth->pt();
      double true_eta = truth->eta();
      double true_phi = truth->phi();

      if (plotAgainstReco_) {pt_denom = rec_pt;}
      else {pt_denom = true_pt;}
      // get true specific quantities
      double true_ChargedHadEnergy;
      double true_NeutralHadEnergy;
      double true_NeutralEmEnergy;
      gettrue (truth, true_ChargedHadEnergy, true_NeutralHadEnergy, true_NeutralEmEnergy);
      double true_NeutralEnergy = true_NeutralHadEnergy + true_NeutralEmEnergy;
      double rec_ChargedHadEnergy = pfj.chargedHadronEnergy();
      double rec_NeutralHadEnergy = pfj.neutralHadronEnergy();
      double rec_NeutralEmEnergy = pfj.neutralEmEnergy();
      double rec_NeutralEnergy = rec_NeutralHadEnergy + rec_NeutralEmEnergy;
      double rec_ChargedMultiplicity = pfj.chargedMultiplicity();
      std::vector <PFCandidatePtr> constituents = pfj.getPFConstituents ();
      std::vector <unsigned int> chMult(9, static_cast<unsigned int>(0)); 
      for (unsigned ic = 0; ic < constituents.size (); ++ic) {
    if ( constituents[ic]->particleId() > 3 ) continue;
    reco::TrackRef trackRef = constituents[ic]->trackRef();
    if ( trackRef.isNull() ) {
      //std::cout << "Warning in entry " << entry_ 
      //        << " : a track with Id " << constituents[ic]->particleId() 
      //        << " has no track ref.." << std::endl;
      continue;
    }
    unsigned int iter = 0; 
    switch (trackRef->algo()) {
    case TrackBase::ctf:
    case TrackBase::iter0:
      iter = 0;
      break;
    case TrackBase::iter1:
      iter = 1;
      break;
    case TrackBase::iter2:
      iter = 2;
      break;
    case TrackBase::iter3:
      iter = 3;
      break;
    case TrackBase::iter4:
      iter = 4;
      break;
    case TrackBase::iter5:
      iter = 5;
      break;
    case TrackBase::iter6:
      iter = 6;
      break;
    case TrackBase::iter8:
      iter = 7;
      //std::cout << "Warning in entry " << entry_ << " : iter = " << trackRef->algo() << std::endl;
      //std::cout << ic << " " << *(constituents[ic]) << std::endl;
      break;
    default:
      iter = 8;
      std::cout << "Warning in entry " << entry_ << " : iter = " << trackRef->algo() << std::endl;
      std::cout << ic << " " << *(constituents[ic]) << std::endl;
      break;
    }
    ++(chMult[iter]);
      }

      bool plot1 = false;
      bool plot2 = false;
      bool plot3 = false;
      bool plot4 = false;
      bool plot5 = false;
      bool plot6 = false;
      bool plot7 = false;
      double cut1 = 0.0001;
      double cut2 = 0.0001;
      double cut3 = 0.0001;
      double cut4 = 0.0001;
      double cut5 = 0.0001;
      double cut6 = 0.0001;
      double cut7 = 0.0001;
      double resPt =0.;
      double resChargedHadEnergy= 0.;
      double resNeutralHadEnergy= 0.;
      double resNeutralEmEnergy= 0.;
      double resNeutralEnergy= 0.;

      double resHCALEnergy = 0.;
      double resNEUTEnergy = 0.;
      if ( rec_NeutralHadEnergy > cut6 && rec_ChargedHadEnergy < cut1 ) { 
    double true_NEUTEnergy = true_NeutralHadEnergy + true_NeutralEmEnergy;
    double true_HCALEnergy = true_NEUTEnergy - rec_NeutralEmEnergy;
    double rec_NEUTEnergy = rec_NeutralHadEnergy+rec_NeutralEmEnergy; 
    double rec_HCALEnergy = rec_NeutralHadEnergy; 
    resHCALEnergy = (rec_HCALEnergy-true_HCALEnergy)/rec_HCALEnergy;
    resNEUTEnergy = (rec_NEUTEnergy-true_NEUTEnergy)/rec_NEUTEnergy;
    if ( rec_NeutralEmEnergy > cut7 ) {
      plot6 = true;
    } else {
      plot7 = true;
    }
      }

      // get relative delta quantities (protect against division by zero!)
      if (true_pt > 0.0001){
    resPt = (rec_pt -true_pt)/true_pt ; 
    plot1 = true;}
      if (true_ChargedHadEnergy > cut1){
    resChargedHadEnergy = (rec_ChargedHadEnergy- true_ChargedHadEnergy)/true_ChargedHadEnergy;
    plot2 = true;}
      if (true_NeutralHadEnergy > cut2){
    resNeutralHadEnergy = (rec_NeutralHadEnergy- true_NeutralHadEnergy)/true_NeutralHadEnergy;
    plot3 = true;}
      else 
    if (rec_NeutralHadEnergy > cut3){
      resNeutralHadEnergy = (rec_NeutralHadEnergy- true_NeutralHadEnergy)/rec_NeutralHadEnergy;
      plot3 = true;}
      if (true_NeutralEmEnergy > cut4){
    resNeutralEmEnergy = (rec_NeutralEmEnergy- true_NeutralEmEnergy)/true_NeutralEmEnergy;
    plot4 = true;}
      if (true_NeutralEnergy > cut5){
    resNeutralEnergy = (rec_NeutralEnergy- true_NeutralEnergy)/true_NeutralEnergy;
    plot5 = true;}
      
      //double deltaEta = algo_->deltaEta(&pfj, truth);
      //double deltaPhi = algo_->deltaPhi(&pfj, truth);

      // Print outliers for further debugging
      if ( ( resPt > 0.2 && true_pt > 100. ) || 
       ( resPt < -0.5 && true_pt > 100. ) ) {
    //if ( ( true_pt > 50. && 
    //     ( ( truth->eta()>3.0 && rec_eta-truth->eta() < -0.1 ) || 
    //       ( truth->eta()<-3.0 && rec_eta-truth->eta() > 0.1 ) ))) {
    std::cout << "Entry " << entry_ 
          << " resPt = " << resPt
          <<" resCharged  " << resChargedHadEnergy
          <<" resNeutralHad  " << resNeutralHadEnergy
          << " resNeutralEm  " << resNeutralEmEnergy
          << " pT (T/R) " << true_pt << "/" << rec_pt 
          << " Eta (T/R) " << truth->eta() << "/" << rec_eta 
          << " Phi (T/R) " << truth->phi() << "/" << rec_phi 
          << std::endl;

    // check overlapping PF jets
    const reco::PFJet* pfoj = 0; 
    double dRo = 1E9;
    for(unsigned j=0; j<pfJets.size(); j++) { 
      const reco::PFJet& pfo = pfJets[j];
      if ( j != i &&  algo_->deltaR(&pfj,&pfo) < dRo && pfo.pt() > 0.25*pfj.pt()) { 
        dRo = algo_->deltaR(&pfj,&pfo); 
        pfoj = &pfo;
      }
    }
    
    // Check overlapping Gen Jet 
    math::XYZTLorentzVector overlappinGenJet(0.,0.,0.,0.);
    const reco::GenJet* genoj = 0;
    double dRgo = 1E9;
    for(unsigned j=0; j<genJets.size(); j++) { 
      const reco::GenJet* gjo = &(genJets[j]);
      if ( gjo != truth && algo_->deltaR(truth,gjo) < dRgo && gjo->pt() > 0.25*truth->pt() ) { 
        dRgo = algo_->deltaR(truth,gjo);
        genoj = gjo;
      }
    }
    
    if ( dRo < 0.8 && dRgo < 0.8 && algo_->deltaR(genoj,pfoj) < 2.*deltaRMax_ ) 
      std::cout << "Excess probably due to overlapping jets (DR = " <<   algo_->deltaR(genoj,pfoj) << "),"
            << " at DeltaR(T/R) = " << dRgo << "/" << dRo  
            << " with pT(T/R) " << genoj->pt() << "/" << pfoj->pt()
            << " and Eta (T/R) " << genoj->eta() << "/" << pfoj->eta()
            << " and Phi (T/R) " << genoj->phi() << "/" << pfoj->phi()
            << std::endl;
      }

      if(std::abs(resPt) > std::abs(resPtMax_)) resPtMax_ = resPt;
      if(std::abs(resChargedHadEnergy) > std::abs(resChargedHadEnergyMax_) ) resChargedHadEnergyMax_ = resChargedHadEnergy;
      if(std::abs(resNeutralHadEnergy) > std::abs(resNeutralHadEnergyMax_) ) resNeutralHadEnergyMax_ = resNeutralHadEnergy;
      if(std::abs(resNeutralEmEnergy) > std::abs(resNeutralEmEnergyMax_) ) resNeutralEmEnergyMax_ = resNeutralEmEnergy;
      if (debug_) {
    cout << i <<"  =========PFJet Pt "<< rec_pt
         << " eta " << rec_eta
         << " phi " << rec_phi
         << " Charged Had Energy " << rec_ChargedHadEnergy
         << " Neutral Had Energy " << rec_NeutralHadEnergy
         << " Neutral elm Energy " << rec_NeutralEmEnergy << endl;
    cout << " matching Gen Jet Pt " << true_pt
         << " eta " << truth->eta()
         << " phi " << truth->phi()
         << " Charged Had Energy " << true_ChargedHadEnergy
         << " Neutral Had Energy " << true_NeutralHadEnergy
         << " Neutral elm Energy " << true_NeutralEmEnergy << endl;
    printPFJet(&pfj);
    //      cout<<pfj.print()<<endl;
    printGenJet(truth);
    //cout <<truth->print()<<endl;
                
    cout << "==============deltaR " << deltaR << "  resPt " << resPt
         << " resChargedHadEnergy " << resChargedHadEnergy
         << " resNeutralHadEnergy " << resNeutralHadEnergy
         << " resNeutralEmEnergy " << resNeutralEmEnergy
         << endl;
      }
            

      if(plot1) {
    if ( rec_eta > 0. ) 
      hDEtavsEta->Fill(true_eta,rec_eta-true_eta);
    else
      hDEtavsEta->Fill(true_eta,-rec_eta+true_eta);
    hDPhivsEta->Fill(true_eta,rec_phi-true_phi);

    hRPtvsEta->Fill(true_eta, resPt);
    hNCHvsEta->Fill(true_eta, rec_ChargedMultiplicity);
    hNCH0vsEta->Fill(true_eta,chMult[0]);
    hNCH1vsEta->Fill(true_eta,chMult[1]);
    hNCH2vsEta->Fill(true_eta,chMult[2]);
    hNCH3vsEta->Fill(true_eta,chMult[3]);
    hNCH4vsEta->Fill(true_eta,chMult[4]);
    hNCH5vsEta->Fill(true_eta,chMult[5]);
    hNCH6vsEta->Fill(true_eta,chMult[6]);
    hNCH7vsEta->Fill(true_eta,chMult[7]);
      }
      if(plot2)hRCHEvsEta->Fill(true_eta, resChargedHadEnergy);
      if(plot5)hRNeutvsEta->Fill(true_eta, resNeutralEnergy);
      if(plot6) { 
    hRHCALvsEta->Fill(true_eta, resHCALEnergy);
    hRNEUTvsEta->Fill(true_eta, resNEUTEnergy);
      }
      if(plot7) {  
    hRHONLvsEta->Fill(true_eta, resHCALEnergy);
    hRNONLvsEta->Fill(true_eta, resNEUTEnergy);
      }

      // fill histograms for relative delta quantitites of matched jets
      // delta Pt or E quantities for Barrel
      if (Barrel){
    if(plot1) { 
      hBRPt->Fill (resPt);
      if ( pt_denom >  20. && pt_denom <  40. ) hBRPt20_40->Fill (resPt);
      if ( pt_denom >  40. && pt_denom <  60. ) hBRPt40_60->Fill (resPt);
      if ( pt_denom >  60. && pt_denom <  80. ) hBRPt60_80->Fill (resPt);
      if ( pt_denom >  80. && pt_denom < 100. ) hBRPt80_100->Fill (resPt);
      if ( pt_denom > 100. && pt_denom < 150. ) hBRPt100_150->Fill (resPt);
      if ( pt_denom > 150. && pt_denom < 200. ) hBRPt150_200->Fill (resPt);
      if ( pt_denom > 200. && pt_denom < 250. ) hBRPt200_250->Fill (resPt);
      if ( pt_denom > 250. && pt_denom < 300. ) hBRPt250_300->Fill (resPt);
      if ( pt_denom > 300. && pt_denom < 400. ) hBRPt300_400->Fill (resPt);
      if ( pt_denom > 400. && pt_denom < 500. ) hBRPt400_500->Fill (resPt);
      if ( pt_denom > 500. && pt_denom < 750. ) hBRPt500_750->Fill (resPt);
      if ( pt_denom > 750. && pt_denom < 1250. ) hBRPt750_1250->Fill (resPt);
      if ( pt_denom > 1250. && pt_denom < 2000. ) hBRPt1250_2000->Fill (resPt);
      if ( pt_denom > 2000. && pt_denom < 5000. ) hBRPt2000_5000->Fill (resPt);
      hBNCH->Fill(rec_ChargedMultiplicity);
      hBNCH0vsPt->Fill(pt_denom,chMult[0]);
      hBNCH1vsPt->Fill(pt_denom,chMult[1]);
      hBNCH2vsPt->Fill(pt_denom,chMult[2]);
      hBNCH3vsPt->Fill(pt_denom,chMult[3]);
      hBNCH4vsPt->Fill(pt_denom,chMult[4]);
      hBNCH5vsPt->Fill(pt_denom,chMult[5]);
      hBNCH6vsPt->Fill(pt_denom,chMult[6]);
      hBNCH7vsPt->Fill(pt_denom,chMult[7]);
      hBNCHvsPt->Fill(pt_denom,rec_ChargedMultiplicity);
      if ( rec_eta > 0. ) 
        hBDEtavsPt->Fill(pt_denom,rec_eta-true_eta);
      else
        hBDEtavsPt->Fill(pt_denom,-rec_eta+true_eta);
      hBDPhivsPt->Fill(pt_denom,rec_phi-true_phi);
    }
    if(plot2)hBRCHE->Fill(resChargedHadEnergy);
    if(plot3)hBRNHE->Fill(resNeutralHadEnergy);
    if(plot4)hBRNEE->Fill(resNeutralEmEnergy);
    if(plot5)hBRneut->Fill(resNeutralEnergy);
    if(plot1)hBRPtvsPt->Fill(pt_denom, resPt);
    if(plot2)hBRCHEvsPt->Fill(pt_denom, resChargedHadEnergy);
    if(plot3)hBRNHEvsPt->Fill(pt_denom, resNeutralHadEnergy);
    if(plot4)hBRNEEvsPt->Fill(pt_denom, resNeutralEmEnergy);
    if(plot5)hBRneutvsPt->Fill(pt_denom, resNeutralEnergy);
    if(plot6) { 
      hBRHCALvsP->Fill(true_E, resHCALEnergy);
      hBRNEUTvsP->Fill(true_E, resNEUTEnergy);
    }
    if(plot7) { 
      hBRHONLvsP->Fill(true_E, resHCALEnergy);
      hBRNONLvsP->Fill(true_E, resNEUTEnergy);
    }

      }
      // delta Pt or E quantities for Endcap
      if (Endcap){
    if(plot1) {
      hERPt->Fill (resPt);
      if ( pt_denom >  20. && pt_denom <  40. ) hERPt20_40->Fill (resPt);
      if ( pt_denom >  40. && pt_denom <  60. ) hERPt40_60->Fill (resPt);
      if ( pt_denom >  60. && pt_denom <  80. ) hERPt60_80->Fill (resPt);
      if ( pt_denom >  80. && pt_denom < 100. ) hERPt80_100->Fill (resPt);
      if ( pt_denom > 100. && pt_denom < 150. ) hERPt100_150->Fill (resPt);
      if ( pt_denom > 150. && pt_denom < 200. ) hERPt150_200->Fill (resPt);
      if ( pt_denom > 200. && pt_denom < 250. ) hERPt200_250->Fill (resPt);
      if ( pt_denom > 250. && pt_denom < 300. ) hERPt250_300->Fill (resPt);
      if ( pt_denom > 300. && pt_denom < 400. ) hERPt300_400->Fill (resPt);
      if ( pt_denom > 400. && pt_denom < 500. ) hERPt400_500->Fill (resPt);
      if ( pt_denom > 500. && pt_denom < 750. ) hERPt500_750->Fill (resPt);
      if ( pt_denom > 750. && pt_denom < 1250. ) hERPt750_1250->Fill (resPt);
      if ( pt_denom > 1250. && pt_denom < 2000. ) hERPt1250_2000->Fill (resPt);
      if ( pt_denom > 2000. && pt_denom < 5000. ) hERPt2000_5000->Fill (resPt);
      hENCH->Fill(rec_ChargedMultiplicity);
      hENCHvsPt->Fill(pt_denom,rec_ChargedMultiplicity);
      hENCH0vsPt->Fill(pt_denom,chMult[0]);
      hENCH1vsPt->Fill(pt_denom,chMult[1]);
      hENCH2vsPt->Fill(pt_denom,chMult[2]);
      hENCH3vsPt->Fill(pt_denom,chMult[3]);
      hENCH4vsPt->Fill(pt_denom,chMult[4]);
      hENCH5vsPt->Fill(pt_denom,chMult[5]);
      hENCH6vsPt->Fill(pt_denom,chMult[6]);
      hENCH7vsPt->Fill(pt_denom,chMult[7]);
      if ( rec_eta > 0. ) 
        hEDEtavsPt->Fill(pt_denom,rec_eta-true_eta);
      else
        hEDEtavsPt->Fill(pt_denom,-rec_eta+true_eta);
      hEDPhivsPt->Fill(pt_denom,rec_phi-true_phi);
    }
    if(plot2)hERCHE->Fill(resChargedHadEnergy);
    if(plot3)hERNHE->Fill(resNeutralHadEnergy);
    if(plot4)hERNEE->Fill(resNeutralEmEnergy);
    if(plot5)hERneut->Fill(resNeutralEnergy);
    if(plot1)hERPtvsPt->Fill(pt_denom, resPt);
    if(plot2)hERCHEvsPt->Fill(pt_denom, resChargedHadEnergy);
    if(plot3)hERNHEvsPt->Fill(pt_denom, resNeutralHadEnergy);
    if(plot4)hERNEEvsPt->Fill(pt_denom, resNeutralEmEnergy);
    if(plot5)hERneutvsPt->Fill(pt_denom, resNeutralEnergy);
    if(plot6) {
      hERHCALvsP->Fill(true_E, resHCALEnergy);
      hERNEUTvsP->Fill(true_E, resNEUTEnergy);
    }
    if(plot7) {
      hERHONLvsP->Fill(true_E, resHCALEnergy);
      hERNONLvsP->Fill(true_E, resNEUTEnergy);
    }
      }                     
      // delta Pt or E quantities for Forward
      if (Forward){
    if(plot1) {
      hFRPt->Fill (resPt);
      if ( pt_denom >  20. && pt_denom <  40. ) hFRPt20_40->Fill (resPt);
      if ( pt_denom >  40. && pt_denom <  60. ) hFRPt40_60->Fill (resPt);
      if ( pt_denom >  60. && pt_denom <  80. ) hFRPt60_80->Fill (resPt);
      if ( pt_denom >  80. && pt_denom < 100. ) hFRPt80_100->Fill (resPt);
      if ( pt_denom > 100. && pt_denom < 150. ) hFRPt100_150->Fill (resPt);
      if ( pt_denom > 150. && pt_denom < 200. ) hFRPt150_200->Fill (resPt);
      if ( pt_denom > 200. && pt_denom < 250. ) hFRPt200_250->Fill (resPt);
      if ( pt_denom > 250. && pt_denom < 300. ) hFRPt250_300->Fill (resPt);
      if ( pt_denom > 300. && pt_denom < 400. ) hFRPt300_400->Fill (resPt);
      if ( pt_denom > 400. && pt_denom < 500. ) hFRPt400_500->Fill (resPt);
      if ( pt_denom > 500. && pt_denom < 750. ) hFRPt500_750->Fill (resPt);
      if ( pt_denom > 750. && pt_denom < 1250. ) hFRPt750_1250->Fill (resPt);
      if ( pt_denom > 1250. && pt_denom < 2000. ) hFRPt1250_2000->Fill (resPt);
      if ( pt_denom > 2000. && pt_denom < 5000. ) hFRPt2000_5000->Fill (resPt);
      if ( rec_eta > 0. ) 
        hFDEtavsPt->Fill(pt_denom,rec_eta-true_eta);
      else
        hFDEtavsPt->Fill(pt_denom,-rec_eta+true_eta);
      hFDPhivsPt->Fill(pt_denom,rec_phi-true_phi);
    }
    if(plot2)hFRCHE->Fill(resChargedHadEnergy);
    if(plot3)hFRNHE->Fill(resNeutralHadEnergy);
    if(plot4)hFRNEE->Fill(resNeutralEmEnergy);
    if(plot5)hFRneut->Fill(resNeutralEnergy);
    if(plot1)hFRPtvsPt->Fill(pt_denom, resPt);
    if(plot2)hFRCHEvsPt->Fill(pt_denom, resChargedHadEnergy);
    if(plot3)hFRNHEvsPt->Fill(pt_denom, resNeutralHadEnergy);
    if(plot4)hFRNEEvsPt->Fill(pt_denom, resNeutralEmEnergy);
    if(plot5)hFRneutvsPt->Fill(pt_denom, resNeutralEnergy);
    if(plot6) {
      hFRHCALvsP->Fill(true_E, resHCALEnergy);
      hFRNEUTvsP->Fill(true_E, resNEUTEnergy);
    }
    if(plot7) {
      hFRHONLvsP->Fill(true_E, resHCALEnergy);
      hFRNONLvsP->Fill(true_E, resNEUTEnergy);
    }
      }                     
    } // end case deltaR < deltaRMax
        
  } // i loop on pf Jets    

  // Increment counter
  entry_++;
}

double PFJetBenchmark::resChargedHadEnergyMax ( ) const

inline

Definition at line 62 of file PFJetBenchmark.h.

References resChargedHadEnergyMax_.

Referenced by PFRootEventManager::processEntry().

{return resChargedHadEnergyMax_;};

double PFJetBenchmark::resNeutralEmEnergyMax ( ) const

inline

Definition at line 64 of file PFJetBenchmark.h.

References resNeutralEmEnergyMax_.

Referenced by PFRootEventManager::processEntry().

{return resNeutralEmEnergyMax_;};

double PFJetBenchmark::resNeutralHadEnergyMax ( ) const

inline

Definition at line 63 of file PFJetBenchmark.h.

References resNeutralHadEnergyMax_.

Referenced by PFRootEventManager::processEntry().

{return resNeutralHadEnergyMax_;};

double PFJetBenchmark::resPtMax ( ) const

inline

Definition at line 61 of file PFJetBenchmark.h.

References resPtMax_.

Referenced by PFRootEventManager::processEntry().

{return resPtMax_;};

void PFJetBenchmark::setup	(	std::string	Filename,
		bool	debug,
		bool	plotAgainstReco = 0,
		bool	onlyTwoJets = 1,
		double	deltaRMax = 0.1,
		std::string	benchmarkLabel_ = "ParticleFlow",
		double	recPt = -1,
		double	maxEta = -1,
		DQMStore *	dbe_store = NULL
	)

Definition at line 63 of file PFJetBenchmark.cc.

References BOOK1D, BOOK2D, gather_cfg::cout, dbe_, DBOOK1D, DBOOK2D, debug, debug_, deltaRMax, deltaRMax_, DSETAXES, ECAL, eta(), file_, HCAL, maxEta, maxEta_cut, NULL, onlyTwoJets, onlyTwoJets_, outputFile_, P, scaleCards::path, phi, plotAgainstReco, plotAgainstReco_, PT, recPt, recPt_cut, SETAXES, and DQMStore::setCurrentFolder().

Referenced by PFJetBenchmarkAnalyzer::PFJetBenchmarkAnalyzer(), and PFRootEventManager::readOptions().

                                     {
  debug_ = debug; 
  plotAgainstReco_ = plotAgainstReco;
  onlyTwoJets_ = onlyTwoJets;
  deltaRMax_ = deltaRMax;
  outputFile_=Filename;
  recPt_cut = recPt;
  maxEta_cut= maxEta;
  file_ = NULL;
  dbe_ = dbe_store;
  // print parameters
  cout<< "PFJetBenchmark Setup parameters =============================================="<<endl;
  cout << "Filename to write histograms " << Filename<<endl;
  cout << "PFJetBenchmark debug " << debug_<< endl;
  cout << "plotAgainstReco " << plotAgainstReco_ << endl;
  cout << "onlyTwoJets " << onlyTwoJets_ << endl;
  cout << "deltaRMax " << deltaRMax << endl;
  cout << "benchmarkLabel " << benchmarkLabel_ << endl;
  cout << "recPt_cut " << recPt_cut << endl;
  cout << "maxEta_cut " << maxEta_cut << endl;
  
  // Book histogram

  // Establish DQM Store
  string path = "PFTask/Benchmarks/"+ benchmarkLabel_ + "/";
  if (plotAgainstReco) path += "Reco"; else path += "Gen";
  if (dbe_) {
    dbe_->setCurrentFolder(path.c_str());
  }
  else {
    file_ = new TFile(outputFile_.c_str(), "recreate");
//    TTree * tr = new TTree("PFTast");
//    tr->Branch("Benchmarks/ParticleFlow")
    cout << "Info: DQM is not available to provide data storage service. Using TFile to save histograms. "<<endl;
  }
  // Jets inclusive  distributions  (Pt > 20 or specified recPt GeV)
  char cutString[35];
  sprintf(cutString,"Jet multiplicity P_{T}>%4.1f GeV", recPt_cut);
  BOOK1D(Njets,cutString,50, 0, 50);

  BOOK1D(jetsPt,"Jets P_{T} Distribution",100, 0, 500);

  sprintf(cutString,"Jets #eta Distribution |#eta|<%4.1f", maxEta_cut);
  BOOK1D(jetsEta,cutString,100, -5, 5);
    
  BOOK2D(RPtvsEta,"#DeltaP_{T}/P_{T} vs #eta",200, -5., 5., 200,-1,1); 
  BOOK2D(RNeutvsEta,"R_{Neutral} vs #eta",200, -5., 5., 200,-1,1); 
  BOOK2D(RNEUTvsEta,"R_{HCAL+ECAL} vs #eta",200, -5., 5., 200,-1,1); 
  BOOK2D(RNONLvsEta,"R_{HCAL+ECAL - Hcal Only} vs #eta",200, -5., 5., 200,-1,1); 
  BOOK2D(RHCALvsEta,"R_{HCAL} vs #eta",200, -5., 5., 200,-1,1); 
  BOOK2D(RHONLvsEta,"R_{HCAL only} vs #eta",200, -5., 5., 200,-1,1); 
  BOOK2D(RCHEvsEta,"R_{Charged} vs #eta",200, -5., 5., 200,-1,1); 
  BOOK2D(NCHvsEta,"N_{Charged} vs #eta",200, -5., 5., 200,0.,200.);
  BOOK2D(NCH0vsEta,"N_{Charged} vs #eta, iter 0",200, -5., 5., 200,0.,200.);
  BOOK2D(NCH1vsEta,"N_{Charged} vs #eta, iter 1",200, -5., 5., 200,0.,200.);
  BOOK2D(NCH2vsEta,"N_{Charged} vs #eta, iter 2",200, -5., 5., 200,0.,200.);
  BOOK2D(NCH3vsEta,"N_{Charged} vs #eta, iter 3",200, -5., 5., 200,0.,200.);
  BOOK2D(NCH4vsEta,"N_{Charged} vs #eta, iter 4",200, -5., 5., 200,0.,200.);
  BOOK2D(NCH5vsEta,"N_{Charged} vs #eta, iter 5",200, -5., 5., 200,0.,200.);
  BOOK2D(NCH6vsEta,"N_{Charged} vs #eta, iter 6",200, -5., 5., 200,0.,200.);
  BOOK2D(NCH7vsEta,"N_{Charged} vs #eta, iter 7",200, -5., 5., 200,0.,200.);
  // delta Pt or E quantities for Barrel
  DBOOK1D(RPt,#DeltaP_{T}/P_{T},80,-1,1);
  DBOOK1D(RCHE,#DeltaE/E (charged had),80,-2,2);
  DBOOK1D(RNHE,#DeltaE/E (neutral had),80,-2,2);
  DBOOK1D(RNEE,#DeltaE/E (neutral em),80,-2,2);
  DBOOK1D(Rneut,#DeltaE/E (neutral),80,-2,2);
  DBOOK1D(NCH, #N_{charged},200,0.,200.);
  DBOOK2D(RPtvsPt,#DeltaP_{T}/P_{T} vs P_{T},250, 0, 500, 200,-1,1);       //used to be 50 bin for each in x-direction
  DBOOK2D(RCHEvsPt,#DeltaE/E (charged had) vs P_{T},250, 0, 500, 120,-1,2);
  DBOOK2D(RNHEvsPt,#DeltaE/E (neutral had) vs P_{T},250, 0, 500, 120,-1,2);
  DBOOK2D(RNEEvsPt,#DeltaE/E (neutral em) vs P_{T},250, 0, 500, 120,-1,2);
  DBOOK2D(RneutvsPt,#DeltaE/E (neutral) vs P_{T},250, 0, 500, 120,-1,2);
  DBOOK2D(NCHvsPt,N_{charged} vs P_{T},250,0,500,200,0.,200.);
  DBOOK2D(NCH0vsPt, N_{charged} vs P_{T} iter 0,250,0,500,200,0.,200.);
  DBOOK2D(NCH1vsPt, N_{charged} vs P_{T} iter 1,250,0,500,200,0.,200.);
  DBOOK2D(NCH2vsPt, N_{charged} vs P_{T} iter 2,250,0,500,200,0.,200.);
  DBOOK2D(NCH3vsPt, N_{charged} vs P_{T} iter 3,250,0,500,200,0.,200.);
  DBOOK2D(NCH4vsPt, N_{charged} vs P_{T} iter 4,250,0,500,200,0.,200.);
  DBOOK2D(NCH5vsPt, N_{charged} vs P_{T} iter 5,250,0,500,200,0.,200.);
  DBOOK2D(NCH6vsPt, N_{charged} vs P_{T} iter 6,250,0,500,200,0.,200.);
  DBOOK2D(NCH7vsPt, N_{charged} vs P_{T} iter 7,250,0,500,200,0.,200.);
  

  DBOOK2D(RNEUTvsP,#DeltaE/E (ECAL+HCAL) vs P,250, 0, 1000, 150,-1.5,1.5);
  DBOOK2D(RNONLvsP,#DeltaE/E (ECAL+HCAL-only) vs P,250, 0, 1000, 150,-1.5,1.5);
  DBOOK2D(RHCALvsP,#DeltaE/E (HCAL) vs P,250, 0, 1000, 150,-1.5,1.5);
  DBOOK2D(RHONLvsP,#DeltaE/E (HCAL only) vs P,250, 0, 1000, 150,-1.5,1.5);
  DBOOK1D(RPt20_40,#DeltaP_{T}/P_{T},80,-1,1);
  DBOOK1D(RPt40_60,#DeltaP_{T}/P_{T},80,-1,1);
  DBOOK1D(RPt60_80,#DeltaP_{T}/P_{T},80,-1,1);
  DBOOK1D(RPt80_100,#DeltaP_{T}/P_{T},80,-1,1);
  DBOOK1D(RPt100_150,#DeltaP_{T}/P_{T},80,-1,1);
  DBOOK1D(RPt150_200,#DeltaP_{T}/P_{T},80,-1,1);
  DBOOK1D(RPt200_250,#DeltaP_{T}/P_{T},80,-1,1);
  DBOOK1D(RPt250_300,#DeltaP_{T}/P_{T},80,-1,1);
  DBOOK1D(RPt300_400,#DeltaP_{T}/P_{T},160,-1,1);
  DBOOK1D(RPt400_500,#DeltaP_{T}/P_{T},160,-1,1);
  DBOOK1D(RPt500_750,#DeltaP_{T}/P_{T},160,-1,1);
  DBOOK1D(RPt750_1250,#DeltaP_{T}/P_{T},160,-1,1);
  DBOOK1D(RPt1250_2000,#DeltaP_{T}/P_{T},160,-1,1);
  DBOOK1D(RPt2000_5000,#DeltaP_{T}/P_{T},160,-1,1);

  DBOOK2D(DEtavsPt,#Delta#eta vs P_{T},1000,0,2000,500,-0.5,0.5);
  DBOOK2D(DPhivsPt,#Delta#phi vs P_{T},1000,0,2000,500,-0.5,0.5);
  BOOK2D(DEtavsEta,"#Delta#eta vs P_{T}",1000,-5.,+5.,500,-0.5,0.5);
  BOOK2D(DPhivsEta,"#Delta#phi vs P_{T}",1000,-5.,+5.,500,-0.5,0.5);
    
 // Set Axis Titles
 
 // Jets inclusive  distributions  (Pt > 20 GeV)
  SETAXES(Njets,"","Multiplicity");
  SETAXES(jetsPt, PT, "Number of Events");
  SETAXES(jetsEta, "#eta", "Number of Events");
  SETAXES(RNeutvsEta, "#eta", "#DeltaE/E (Neutral)");
  SETAXES(RNEUTvsEta, "#eta", "#DeltaE/E (ECAL+HCAL)");
  SETAXES(RNONLvsEta, "#eta", "#DeltaE/E (ECAL+HCAL-only)");
  SETAXES(RHCALvsEta, "#eta", "#DeltaE/E (HCAL)");
  SETAXES(RHONLvsEta, "#eta", "#DeltaE/E (HCAL Only)");
  SETAXES(RCHEvsEta, "#eta", "#DeltaE/E (Charged)");
  SETAXES(RPtvsEta, "#eta", "#DeltaP_{T}/P_{T}");
  SETAXES(DEtavsEta, "#eta", "#Delta#eta");
  SETAXES(DPhivsEta,"#eta", "#Delta#phi");
  // delta Pt or E quantities for Barrel and Endcap
  DSETAXES(RPt, "#DeltaP_{T}/P_{T}", "Events");
  DSETAXES(RPt20_40, "#DeltaP_{T}/P_{T}", "Events");
  DSETAXES(RPt40_60, "#DeltaP_{T}/P_{T}", "Events");
  DSETAXES(RPt60_80, "#DeltaP_{T}/P_{T}", "Events");
  DSETAXES(RPt80_100, "#DeltaP_{T}/P_{T}", "Events");
  DSETAXES(RPt100_150, "#DeltaP_{T}/P_{T}", "Events");
  DSETAXES(RPt150_200, "#DeltaP_{T}/P_{T}", "Events");
  DSETAXES(RPt200_250, "#DeltaP_{T}/P_{T}", "Events");
  DSETAXES(RPt250_300, "#DeltaP_{T}/P_{T}", "Events");
  DSETAXES(RPt300_400, "#DeltaP_{T}/P_{T}", "Events");
  DSETAXES(RPt400_500, "#DeltaP_{T}/P_{T}", "Events");
  DSETAXES(RPt500_750, "#DeltaP_{T}/P_{T}", "Events");
  DSETAXES(RPt750_1250, "#DeltaP_{T}/P_{T}", "Events");
  DSETAXES(RPt1250_2000, "#DeltaP_{T}/P_{T}", "Events");
  DSETAXES(RPt2000_5000, "#DeltaP_{T}/P_{T}", "Events");
  DSETAXES(RCHE, "#DeltaE/E(charged had)", "Events");
  DSETAXES(RNHE, "#DeltaE/E(neutral had)", "Events");
  DSETAXES(RNEE, "#DeltaE/E(neutral em)", "Events");
  DSETAXES(Rneut, "#DeltaE/E(neutral)", "Events");
  DSETAXES(RPtvsPt, PT, "#DeltaP_{T}/P_{T}");
  DSETAXES(RCHEvsPt, PT, "#DeltaE/E(charged had)");
  DSETAXES(RNHEvsPt, PT, "#DeltaE/E(neutral had)");
  DSETAXES(RNEEvsPt, PT, "#DeltaE/E(neutral em)");
  DSETAXES(RneutvsPt, PT, "#DeltaE/E(neutral)");
  DSETAXES(RHCALvsP, P, "#DeltaE/E(HCAL)");
  DSETAXES(RHONLvsP, P, "#DeltaE/E(HCAL-only)");
  DSETAXES(RNEUTvsP, P, "#DeltaE/E(ECAL+HCAL)");
  DSETAXES(RNONLvsP, P, "#DeltaE/E(ECAL+HCAL-only)");
  DSETAXES(DEtavsPt, PT, "#Delta#eta");
  DSETAXES(DPhivsPt, PT, "#Delta#phi");

}

void PFJetBenchmark::write

(

void

)

Definition at line 46 of file PFJetBenchmark.cc.

References gather_cfg::cout, dbe_, file_, outputFile_, and DQMStore::save().

Referenced by PFJetBenchmarkAnalyzer::endJob(), and PFRootEventManager::write().

                           {
   // Store the DAQ Histograms 
  if (outputFile_.size() != 0) {
    if (dbe_)
          dbe_->save(outputFile_.c_str());
    // use bare Root if no DQM (FWLite applications)
    else if (file_) {
       file_->Write(outputFile_.c_str());
       cout << "Benchmark output written to file " << outputFile_.c_str() << endl;
       file_->Close();
       }
  }
  else 
    cout << "No output file specified ("<<outputFile_<<"). Results will not be saved!" << endl;
    
} 

Member Data Documentation

PFBenchmarkAlgo* PFJetBenchmark::algo_

protected

Definition at line 225 of file PFJetBenchmark.h.

Referenced by process().

DQMStore* PFJetBenchmark::dbe_

protected

Definition at line 237 of file PFJetBenchmark.h.

Referenced by setup(), and write().

bool PFJetBenchmark::debug_

protected

Definition at line 226 of file PFJetBenchmark.h.

Referenced by process(), and setup().

double PFJetBenchmark::deltaRMax_

protected

Definition at line 229 of file PFJetBenchmark.h.

Referenced by process(), and setup().

unsigned int PFJetBenchmark::entry_

protected

Definition at line 236 of file PFJetBenchmark.h.

Referenced by process().

TFile* PFJetBenchmark::file_

private

Definition at line 70 of file PFJetBenchmark.h.

Referenced by setup(), write(), and ~PFJetBenchmark().

TH2F* PFJetBenchmark::hBDEtavsPt

private

Definition at line 125 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hBDPhivsPt

private

Definition at line 126 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hBNCH

private

Definition at line 128 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hBNCH0vsPt

private

Definition at line 129 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hBNCH1vsPt

private

Definition at line 130 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hBNCH2vsPt

private

Definition at line 131 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hBNCH3vsPt

private

Definition at line 132 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hBNCH4vsPt

private

Definition at line 133 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hBNCH5vsPt

private

Definition at line 134 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hBNCH6vsPt

private

Definition at line 135 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hBNCH7vsPt

private

Definition at line 136 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hBNCHvsPt

private

Definition at line 127 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hBRCHE

private

Definition at line 112 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hBRCHEvsPt

private

Definition at line 117 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hBRHCALvsP

private

Definition at line 123 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hBRHONLvsP

private

Definition at line 124 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hBRNEE

private

Definition at line 114 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hBRNEEvsPt

private

Definition at line 119 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hBRneut

private

Definition at line 115 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hBRNEUTvsP

private

Definition at line 121 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hBRneutvsPt

private

Definition at line 120 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hBRNHE

private

Definition at line 113 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hBRNHEvsPt

private

Definition at line 118 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hBRNONLvsP

private

Definition at line 122 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hBRPt

private

Definition at line 97 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hBRPt100_150

private

Definition at line 102 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hBRPt1250_2000

private

Definition at line 110 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hBRPt150_200

private

Definition at line 103 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hBRPt2000_5000

private

Definition at line 111 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hBRPt200_250

private

Definition at line 104 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hBRPt20_40

private

Definition at line 98 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hBRPt250_300

private

Definition at line 105 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hBRPt300_400

private

Definition at line 106 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hBRPt400_500

private

Definition at line 107 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hBRPt40_60

private

Definition at line 99 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hBRPt500_750

private

Definition at line 108 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hBRPt60_80

private

Definition at line 100 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hBRPt750_1250

private

Definition at line 109 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hBRPt80_100

private

Definition at line 101 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hBRPtvsPt

private

Definition at line 116 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hDEtavsEta

private

Definition at line 78 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hDPhivsEta

private

Definition at line 79 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hEDEtavsPt

private

Definition at line 167 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hEDPhivsPt

private

Definition at line 168 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hENCH

private

Definition at line 170 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hENCH0vsPt

private

Definition at line 171 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hENCH1vsPt

private

Definition at line 172 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hENCH2vsPt

private

Definition at line 173 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hENCH3vsPt

private

Definition at line 174 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hENCH4vsPt

private

Definition at line 175 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hENCH5vsPt

private

Definition at line 176 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hENCH6vsPt

private

Definition at line 177 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hENCH7vsPt

private

Definition at line 178 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hENCHvsPt

private

Definition at line 169 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hERCHE

private

Definition at line 154 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hERCHEvsPt

private

Definition at line 159 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hERHCALvsP

private

Definition at line 165 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hERHONLvsP

private

Definition at line 166 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hERNEE

private

Definition at line 156 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hERNEEvsPt

private

Definition at line 161 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hERneut

private

Definition at line 157 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hERNEUTvsP

private

Definition at line 163 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hERneutvsPt

private

Definition at line 162 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hERNHE

private

Definition at line 155 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hERNHEvsPt

private

Definition at line 160 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hERNONLvsP

private

Definition at line 164 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hERPt

private

Definition at line 139 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hERPt100_150

private

Definition at line 144 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hERPt1250_2000

private

Definition at line 152 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hERPt150_200

private

Definition at line 145 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hERPt2000_5000

private

Definition at line 153 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hERPt200_250

private

Definition at line 146 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hERPt20_40

private

Definition at line 140 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hERPt250_300

private

Definition at line 147 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hERPt300_400

private

Definition at line 148 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hERPt400_500

private

Definition at line 149 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hERPt40_60

private

Definition at line 141 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hERPt500_750

private

Definition at line 150 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hERPt60_80

private

Definition at line 142 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hERPt750_1250

private

Definition at line 151 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hERPt80_100

private

Definition at line 143 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hERPtvsPt

private

Definition at line 158 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hFDEtavsPt

private

Definition at line 209 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hFDPhivsPt

private

Definition at line 210 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hFNCH

private

Definition at line 212 of file PFJetBenchmark.h.

TH2F* PFJetBenchmark::hFNCH0vsPt

private

Definition at line 213 of file PFJetBenchmark.h.

TH2F* PFJetBenchmark::hFNCH1vsPt

private

Definition at line 214 of file PFJetBenchmark.h.

TH2F* PFJetBenchmark::hFNCH2vsPt

private

Definition at line 215 of file PFJetBenchmark.h.

TH2F* PFJetBenchmark::hFNCH3vsPt

private

Definition at line 216 of file PFJetBenchmark.h.

TH2F* PFJetBenchmark::hFNCH4vsPt

private

Definition at line 217 of file PFJetBenchmark.h.

TH2F* PFJetBenchmark::hFNCH5vsPt

private

Definition at line 218 of file PFJetBenchmark.h.

TH2F* PFJetBenchmark::hFNCH6vsPt

private

Definition at line 219 of file PFJetBenchmark.h.

TH2F* PFJetBenchmark::hFNCH7vsPt

private

Definition at line 220 of file PFJetBenchmark.h.

TH2F* PFJetBenchmark::hFNCHvsPt

private

Definition at line 211 of file PFJetBenchmark.h.

TH1F* PFJetBenchmark::hFRCHE

private

Definition at line 196 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hFRCHEvsPt

private

Definition at line 201 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hFRHCALvsP

private

Definition at line 207 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hFRHONLvsP

private

Definition at line 208 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hFRNEE

private

Definition at line 198 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hFRNEEvsPt

private

Definition at line 203 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hFRneut

private

Definition at line 199 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hFRNEUTvsP

private

Definition at line 205 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hFRneutvsPt

private

Definition at line 204 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hFRNHE

private

Definition at line 197 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hFRNHEvsPt

private

Definition at line 202 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hFRNONLvsP

private

Definition at line 206 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hFRPt

private

Definition at line 181 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hFRPt100_150

private

Definition at line 186 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hFRPt1250_2000

private

Definition at line 194 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hFRPt150_200

private

Definition at line 187 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hFRPt2000_5000

private

Definition at line 195 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hFRPt200_250

private

Definition at line 188 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hFRPt20_40

private

Definition at line 182 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hFRPt250_300

private

Definition at line 189 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hFRPt300_400

private

Definition at line 190 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hFRPt400_500

private

Definition at line 191 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hFRPt40_60

private

Definition at line 183 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hFRPt500_750

private

Definition at line 192 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hFRPt60_80

private

Definition at line 184 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hFRPt750_1250

private

Definition at line 193 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hFRPt80_100

private

Definition at line 185 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hFRPtvsPt

private

Definition at line 200 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hjetsEta

private

Definition at line 76 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hjetsPt

private

Definition at line 75 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hNCH0vsEta

private

Definition at line 87 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hNCH1vsEta

private

Definition at line 88 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hNCH2vsEta

private

Definition at line 89 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hNCH3vsEta

private

Definition at line 90 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hNCH4vsEta

private

Definition at line 91 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hNCH5vsEta

private

Definition at line 92 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hNCH6vsEta

private

Definition at line 93 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hNCH7vsEta

private

Definition at line 94 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hNCHvsEta

private

Definition at line 86 of file PFJetBenchmark.h.

Referenced by process().

TH1F* PFJetBenchmark::hNjets

private

Definition at line 74 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hRCHEvsEta

private

Definition at line 85 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hRHCALvsEta

private

Definition at line 83 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hRHONLvsEta

private

Definition at line 84 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hRNeutvsEta

private

Definition at line 80 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hRNEUTvsEta

private

Definition at line 81 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hRNONLvsEta

private

Definition at line 82 of file PFJetBenchmark.h.

Referenced by process().

TH2F* PFJetBenchmark::hRPtvsEta

private

Definition at line 77 of file PFJetBenchmark.h.

Referenced by process().

double PFJetBenchmark::maxEta_cut

protected

Definition at line 235 of file PFJetBenchmark.h.

Referenced by process(), and setup().

bool PFJetBenchmark::onlyTwoJets_

protected

Definition at line 228 of file PFJetBenchmark.h.

Referenced by process(), and setup().

std::string PFJetBenchmark::outputFile_

private

Definition at line 222 of file PFJetBenchmark.h.

Referenced by setup(), and write().

bool PFJetBenchmark::plotAgainstReco_

protected

Definition at line 227 of file PFJetBenchmark.h.

Referenced by process(), and setup().

double PFJetBenchmark::recPt_cut

protected

Definition at line 234 of file PFJetBenchmark.h.

Referenced by process(), and setup().

double PFJetBenchmark::resChargedHadEnergyMax_

protected

Definition at line 231 of file PFJetBenchmark.h.

Referenced by process(), and resChargedHadEnergyMax().

double PFJetBenchmark::resNeutralEmEnergyMax_

protected

Definition at line 233 of file PFJetBenchmark.h.

Referenced by process(), and resNeutralEmEnergyMax().

double PFJetBenchmark::resNeutralHadEnergyMax_

protected

Definition at line 232 of file PFJetBenchmark.h.

Referenced by process(), and resNeutralHadEnergyMax().

double PFJetBenchmark::resPtMax_

protected

Definition at line 230 of file PFJetBenchmark.h.

Referenced by process(), and resPtMax().