JetAnalyzer_HeavyIons_matching Class Reference

#include <JetAnalyzer_HeavyIons_matching.h>

Inheritance diagram for JetAnalyzer_HeavyIons_matching:

Classes
struct	CompareMatchedJets
struct	MyJet

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	JetAnalyzer_HeavyIons_matching (const edm::ParameterSet &)
	~JetAnalyzer_HeavyIons_matching () override
Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default
	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
	~DQMEDAnalyzer () override=default

Private Types
typedef std::multiset< ABJetPair >::iterator	ABItr
typedef std::pair< MyJet, MyJet >	ABJetPair
typedef std::multiset< ABJetPair, CompareMatchedJets >	ABMatchedJets

Static Private Member Functions
static float	deltaRR (float eta1, float phi1, float eta2, float phi2)

Private Attributes
edm::EDGetTokenT< reco::BasicJetCollection >	basicJetsToken_
edm::EDGetTokenT< reco::CaloJetCollection >	caloJet1Token_
edm::EDGetTokenT< reco::CaloJetCollection >	caloJet2Token_
std::string	JetType1
std::string	JetType2
edm::EDGetTokenT< reco::JPTJetCollection >	jptJetsToken_
MonitorElement *	mChargedEmEnergy_Jet1_unmatched
MonitorElement *	mChargedEmEnergy_Jet2_unmatched
MonitorElement *	mChargedHadEnergyFraction_Jet1_unmatched
MonitorElement *	mChargedHadEnergyFraction_Jet2_unmatched
MonitorElement *	mChargedHadronEnergy_Jet1_unmatched
MonitorElement *	mChargedHadronEnergy_Jet2_unmatched
MonitorElement *	mChargedMuEnergy_Jet1_unmatched
MonitorElement *	mChargedMuEnergy_Jet2_unmatched
MonitorElement *	mElectronEnergyFraction_Jet1_unmatched
MonitorElement *	mElectronEnergyFraction_Jet2_unmatched
MonitorElement *	mEmEnergy_Jet1_unmatched
MonitorElement *	mEmEnergy_Jet2_unmatched
MonitorElement *	mHadEnergy_Jet1_unmatched
MonitorElement *	mHadEnergy_Jet2_unmatched
edm::InputTag	mInputJet1Collection
edm::InputTag	mInputJet2Collection
MonitorElement *	mMuonEnergyFraction_Jet1_unmatched
MonitorElement *	mMuonEnergyFraction_Jet2_unmatched
MonitorElement *	mNeutralEmEnergy_Jet1_unmatched
MonitorElement *	mNeutralEmEnergy_Jet2_unmatched
MonitorElement *	mNeutralHadEnergyFraction_Jet1_unmatched
MonitorElement *	mNeutralHadEnergyFraction_Jet2_unmatched
MonitorElement *	mNeutralHadronEnergy_Jet1_unmatched
MonitorElement *	mNeutralHadronEnergy_Jet2_unmatched
std::string	mOutputFile
MonitorElement *	mPhotonEnergyFraction_Jet1_unmatched
MonitorElement *	mPhotonEnergyFraction_Jet2_unmatched
MonitorElement *	mpT_Jet1_matched
MonitorElement *	mpT_Jet1_unmatched
MonitorElement *	mpT_Jet2_matched
MonitorElement *	mpT_Jet2_unmatched
MonitorElement *	mpT_ratio_Jet1Jet2
double	mRecoDelRMatch
double	mRecoJetEtaCut
double	mRecoJetPtThreshold
edm::EDGetTokenT< reco::PFJetCollection >	pfJetsToken_

Detailed Description

Definition at line 83 of file JetAnalyzer_HeavyIons_matching.h.

Member Typedef Documentation

typedef std::multiset< ABJetPair >::iterator JetAnalyzer_HeavyIons_matching::ABItr

private

Definition at line 176 of file JetAnalyzer_HeavyIons_matching.h.

typedef std::pair< MyJet, MyJet > JetAnalyzer_HeavyIons_matching::ABJetPair

private

Definition at line 157 of file JetAnalyzer_HeavyIons_matching.h.

typedef std::multiset< ABJetPair, CompareMatchedJets > JetAnalyzer_HeavyIons_matching::ABMatchedJets

private

Definition at line 175 of file JetAnalyzer_HeavyIons_matching.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 19 of file JetAnalyzer_HeavyIons_matching.cc.

References basicJetsToken_, caloJet1Token_, caloJet2Token_, JetType1, JetType2, edm::InputTag::label(), mChargedEmEnergy_Jet1_unmatched, mChargedEmEnergy_Jet2_unmatched, mChargedHadEnergyFraction_Jet1_unmatched, mChargedHadEnergyFraction_Jet2_unmatched, mChargedHadronEnergy_Jet1_unmatched, mChargedHadronEnergy_Jet2_unmatched, mChargedMuEnergy_Jet1_unmatched, mChargedMuEnergy_Jet2_unmatched, mElectronEnergyFraction_Jet1_unmatched, mElectronEnergyFraction_Jet2_unmatched, mEmEnergy_Jet1_unmatched, mEmEnergy_Jet2_unmatched, mHadEnergy_Jet1_unmatched, mHadEnergy_Jet2_unmatched, mInputJet1Collection, mInputJet2Collection, mMuonEnergyFraction_Jet1_unmatched, mMuonEnergyFraction_Jet2_unmatched, mNeutralEmEnergy_Jet1_unmatched, mNeutralEmEnergy_Jet2_unmatched, mNeutralHadEnergyFraction_Jet1_unmatched, mNeutralHadEnergyFraction_Jet2_unmatched, mNeutralHadronEnergy_Jet1_unmatched, mNeutralHadronEnergy_Jet2_unmatched, mPhotonEnergyFraction_Jet1_unmatched, mPhotonEnergyFraction_Jet2_unmatched, mpT_Jet1_matched, mpT_Jet1_unmatched, mpT_Jet2_matched, mpT_Jet2_unmatched, mpT_ratio_Jet1Jet2, pfJetsToken_, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                             :
  mInputJet1Collection           (iConfig.getParameter<edm::InputTag>       ("src_Jet1")),
  mInputJet2Collection           (iConfig.getParameter<edm::InputTag>       ("src_Jet2")),
  JetType1                       (iConfig.getUntrackedParameter<std::string>("Jet1")),
  JetType2                       (iConfig.getUntrackedParameter<std::string>("Jet2")),
  mRecoJetPtThreshold            (iConfig.getParameter<double>              ("recoJetPtThreshold")),
  mRecoDelRMatch                 (iConfig.getParameter<double>              ("recoDelRMatch")),
  mRecoJetEtaCut                 (iConfig.getParameter<double>              ("recoJetEtaCut"))
{
  std::string inputCollectionLabelJet1(mInputJet1Collection.label());
  std::string inputCollectionLabelJet2(mInputJet2Collection.label());
  
  //consumes
  
  if(std::string("VsCalo") == JetType1) caloJet1Token_ = consumes<reco::CaloJetCollection>(mInputJet1Collection);
  if(std::string("VsPF") == JetType1) pfJetsToken_ = consumes<reco::PFJetCollection>(mInputJet1Collection);
  if(std::string("PuCalo") == JetType1) caloJet2Token_ = consumes<reco::CaloJetCollection>(mInputJet1Collection);
  if(std::string("PuPF") == JetType1) basicJetsToken_ = consumes<reco::BasicJetCollection>(mInputJet1Collection);

  if(std::string("VsCalo") == JetType2) caloJet1Token_ = consumes<reco::CaloJetCollection>(mInputJet2Collection);
  if(std::string("VsPF") == JetType2) pfJetsToken_ = consumes<reco::PFJetCollection>(mInputJet2Collection);
  if(std::string("PuCalo") == JetType2) caloJet2Token_ = consumes<reco::CaloJetCollection>(mInputJet2Collection);
  if(std::string("PuPF") == JetType2) basicJetsToken_ = consumes<reco::BasicJetCollection>(mInputJet2Collection);

  // initialize the Jet matching histograms

  mpT_ratio_Jet1Jet2 = nullptr;
  mpT_Jet1_matched = nullptr;
  mpT_Jet2_matched = nullptr;
  mpT_Jet1_unmatched = nullptr;
  mpT_Jet2_unmatched = nullptr;

  // we need to add histograms which will hold the hadronic and electromagnetic energy content for the unmatched Jets.
  if(std::string("VsCalo") == JetType1 || std::string("PuCalo") == JetType1){
    mHadEnergy_Jet1_unmatched = nullptr;
    mEmEnergy_Jet1_unmatched = nullptr; 
  }
  if(std::string("VsCalo") == JetType2 || std::string("PuCalo") == JetType2) {
    mHadEnergy_Jet2_unmatched = nullptr;
    mEmEnergy_Jet2_unmatched = nullptr; 
  }

  if(std::string("VsPF") == JetType1){
    mChargedHadronEnergy_Jet1_unmatched = nullptr; 
    mNeutralHadronEnergy_Jet1_unmatched = nullptr;
    mChargedEmEnergy_Jet1_unmatched = nullptr;
    mNeutralEmEnergy_Jet1_unmatched = nullptr;
    mChargedMuEnergy_Jet1_unmatched = nullptr;
    mChargedHadEnergyFraction_Jet1_unmatched = nullptr;
    mNeutralHadEnergyFraction_Jet1_unmatched = nullptr;
    mPhotonEnergyFraction_Jet1_unmatched = nullptr;
    mElectronEnergyFraction_Jet1_unmatched = nullptr;
    mMuonEnergyFraction_Jet1_unmatched = nullptr;
  }
  
  if(std::string("VsPF") == JetType2){
    mChargedHadronEnergy_Jet2_unmatched = nullptr; 
    mNeutralHadronEnergy_Jet2_unmatched = nullptr;
    mChargedEmEnergy_Jet2_unmatched = nullptr;
    mNeutralEmEnergy_Jet2_unmatched = nullptr;
    mChargedMuEnergy_Jet2_unmatched = nullptr;
    mChargedHadEnergyFraction_Jet2_unmatched = nullptr;
    mNeutralHadEnergyFraction_Jet2_unmatched = nullptr;
    mPhotonEnergyFraction_Jet2_unmatched = nullptr;
    mElectronEnergyFraction_Jet2_unmatched = nullptr;
    mMuonEnergyFraction_Jet2_unmatched = nullptr;
  }
  
  
}

JetAnalyzer_HeavyIons_matching::~JetAnalyzer_HeavyIons_matching ( )

override

Definition at line 154 of file JetAnalyzer_HeavyIons_matching.cc.

{}

Member Function Documentation

void JetAnalyzer_HeavyIons_matching::analyze ( const edm::Event &  mEvent, const edm::EventSetup &  mSetup  )

override

Definition at line 180 of file JetAnalyzer_HeavyIons_matching.cc.

References basicJetsToken_, caloJet1Token_, caloJet2Token_, deltaRR(), PVValHelper::eta, JetAnalyzer_HeavyIons_matching::MyJet::eta, MonitorElement::Fill(), edm::Event::getByToken(), JetAnalyzer_HeavyIons_matching::MyJet::id, JetType1, JetType2, mChargedEmEnergy_Jet1_unmatched, mChargedEmEnergy_Jet2_unmatched, mChargedHadEnergyFraction_Jet1_unmatched, mChargedHadEnergyFraction_Jet2_unmatched, mChargedHadronEnergy_Jet1_unmatched, mChargedHadronEnergy_Jet2_unmatched, mChargedMuEnergy_Jet1_unmatched, mChargedMuEnergy_Jet2_unmatched, mElectronEnergyFraction_Jet1_unmatched, mElectronEnergyFraction_Jet2_unmatched, mEmEnergy_Jet1_unmatched, mEmEnergy_Jet2_unmatched, mHadEnergy_Jet1_unmatched, mHadEnergy_Jet2_unmatched, mMuonEnergyFraction_Jet1_unmatched, mMuonEnergyFraction_Jet2_unmatched, mNeutralEmEnergy_Jet1_unmatched, mNeutralEmEnergy_Jet2_unmatched, mNeutralHadEnergyFraction_Jet1_unmatched, mNeutralHadEnergyFraction_Jet2_unmatched, mNeutralHadronEnergy_Jet1_unmatched, mNeutralHadronEnergy_Jet2_unmatched, mPhotonEnergyFraction_Jet1_unmatched, mPhotonEnergyFraction_Jet2_unmatched, mpT_Jet1_matched, mpT_Jet1_unmatched, mpT_Jet2_matched, mpT_Jet2_unmatched, mpT_ratio_Jet1Jet2, mRecoDelRMatch, mRecoJetEtaCut, mRecoJetPtThreshold, pfJetBenchmark_cfi::pfJets, pfJetsToken_, JetAnalyzer_HeavyIons_matching::MyJet::phi, EnergyCorrector::pt, JetAnalyzer_HeavyIons_matching::MyJet::pt, and AlCaHLTBitMon_QueryRunRegistry::string.

{

  // Get the Jet collection
  //----------------------------------------------------------------------------
  
  std::vector<const Jet*> recoJet1;
  recoJet1.clear();
  std::vector<const Jet*> recoJet2;
  recoJet2.clear();
  
  edm::Handle<CaloJetCollection>  caloJet1;
  edm::Handle<CaloJetCollection>  caloJet2;
  edm::Handle<JPTJetCollection>   jptJets;
  edm::Handle<PFJetCollection>    pfJets;
  edm::Handle<BasicJetCollection> basicJets;

  if(std::string("VsCalo") == JetType1) {
    mEvent.getByToken(caloJet1Token_, caloJet1);
    for (unsigned ijet=0; ijet<caloJet1->size(); ++ijet) recoJet1.push_back(&(*caloJet1)[ijet]);
  }
  if(std::string("PuCalo") == JetType1) {
    mEvent.getByToken(caloJet2Token_, caloJet1);
    for (unsigned ijet=0; ijet<caloJet1->size(); ++ijet) recoJet1.push_back(&(*caloJet1)[ijet]);
  }
  if(std::string("VsPF") == JetType1) {
    mEvent.getByToken(pfJetsToken_, pfJets);
    for (unsigned ijet=0; ijet<pfJets->size(); ++ijet) recoJet1.push_back(&(*pfJets)[ijet]);
  }
  if(std::string("PuPF") == JetType1) {
    mEvent.getByToken(basicJetsToken_, basicJets);
    for (unsigned ijet=0; ijet<basicJets->size(); ++ijet) recoJet1.push_back(&(*basicJets)[ijet]);
  }

  if(std::string("VsCalo") == JetType2) {
    mEvent.getByToken(caloJet1Token_, caloJet2);
    for (unsigned ijet=0; ijet<caloJet2->size(); ++ijet) recoJet2.push_back(&(*caloJet2)[ijet]);
  }
  if(std::string("PuCalo") == JetType2) {
    mEvent.getByToken(caloJet2Token_, caloJet2);
    for (unsigned ijet=0; ijet<caloJet2->size(); ++ijet) recoJet2.push_back(&(*caloJet2)[ijet]);
  }
  if(std::string("VsPF") == JetType2) {
    mEvent.getByToken(pfJetsToken_, pfJets);
    for (unsigned ijet=0; ijet<pfJets->size(); ++ijet) recoJet2.push_back(&(*pfJets)[ijet]);
  }
  if(std::string("PuPF") == JetType2) {
    mEvent.getByToken(basicJetsToken_, basicJets);
    for (unsigned ijet=0; ijet<basicJets->size(); ++ijet) recoJet2.push_back(&(*basicJets)[ijet]);
  }

  // start to perform the matching - between recoJet1 and recoJet2.

  Int_t Jet1_nref = recoJet1.size();
  Int_t Jet2_nref = recoJet2.size();

  int jet1 = 0;
  int jet2 = 0; 
  
  std::vector <MyJet> vJet1, vJet2;
  std::vector <int> Jet1_ID(Jet1_nref), Jet2_ID(Jet2_nref);
  
  if(Jet1_nref == 0 || Jet2_nref == 0) return;
  
  for(unsigned ijet1 = 0; ijet1 < recoJet1.size(); ++ijet1){


    if(recoJet1[ijet1]->pt() < mRecoJetPtThreshold) continue;
    if(fabs(recoJet1[ijet1]->eta()) < mRecoJetEtaCut) continue;

    MyJet JET1;
    JET1.eta = recoJet1[ijet1]->eta();
    JET1.phi = recoJet1[ijet1]->phi();
    JET1.pt  = recoJet1[ijet1]->pt();
    JET1.id  = ijet1; 

    vJet1.push_back(JET1);
    jet1++;

  }// first jet loop

  for(unsigned ijet2 = 0; ijet2 < recoJet2.size(); ++ijet2){

    if(recoJet2[ijet2]->pt() < mRecoJetPtThreshold) continue;
    if(fabs(recoJet2[ijet2]->eta()) < mRecoJetEtaCut) continue;

    MyJet JET2;
    JET2.eta = recoJet2[ijet2]->eta();
    JET2.phi = recoJet2[ijet2]->phi();
    JET2.pt  = recoJet2[ijet2]->pt();
    JET2.id  = ijet2; 

    vJet2.push_back(JET2);
    jet2++;

  }// second jet loop

  bool onlyJet1     = (jet1>0  && jet2==0)  ? true : false;
  bool onlyJet2     = (jet1==0 && jet2 >0)  ? true : false;
  bool bothJet1Jet2 = (jet1>0  && jet2 >0)  ? true : false;
  
  int matchedJets   = 0;
  int unmatchedJet1 = 0;
  int unmatchedJet2 = 0;  

  std::vector < MyJet >::const_iterator iJet, jJet;

  if(onlyJet1) {

    for(iJet = vJet1.begin(); iJet != vJet1.end(); ++iJet){

      int pj = (*iJet).id;
      
      mpT_Jet1_unmatched->Fill(recoJet1[pj]->pt());

    }

  }else if(onlyJet2) {

    for(iJet = vJet2.begin(); iJet != vJet2.end(); ++iJet){

      int cj = (*iJet).id;

      mpT_Jet2_unmatched->Fill(recoJet2[cj]->pt());
    }
    
  }else if (bothJet1Jet2){

    ABMatchedJets mABMatchedJets; 

    for(iJet = vJet1.begin(); iJet != vJet1.end(); ++iJet){
      for(jJet = vJet2.begin(); jJet != vJet2.end(); ++jJet){
    mABMatchedJets.insert(std::make_pair(*iJet, *jJet));    
      }
    }

    
    ABItr itr;
    // matched Jets matching Jet 1 to Jet 2 
    for(itr = mABMatchedJets.begin(); itr != mABMatchedJets.end(); ++itr){

      ABJetPair jetpair = (*itr);
      MyJet Aj = jetpair.first;
      MyJet Bj = jetpair.second;

      float delr = JetAnalyzer_HeavyIons_matching::deltaRR(Bj.eta, Bj.phi, Aj.eta, Aj.phi);

      if( delr < mRecoDelRMatch && Jet1_ID[Aj.id] == 0){

    mpT_ratio_Jet1Jet2->Fill((Float_t) recoJet2[Bj.id]->pt()/recoJet1[Aj.id]->pt());

    mpT_Jet1_matched->Fill(recoJet1[Aj.id]->pt());
    mpT_Jet2_matched->Fill(recoJet2[Bj.id]->pt());
    
    Jet1_ID[Aj.id] = 1;
    Jet2_ID[Bj.id] = 1;

    matchedJets++; 


      }

    }

    // for unmatched Jets 
    for(itr = mABMatchedJets.begin(); itr != mABMatchedJets.end(); ++itr){

      ABJetPair jetpair = (*itr);

      MyJet Aj = jetpair.first;
      MyJet Bj = jetpair.second;

      if(Jet1_ID[Aj.id] == 0) {

    mpT_Jet1_unmatched->Fill(recoJet1[Aj.id]->pt());
    unmatchedJet1++;
    Jet1_ID[Aj.id] = 1;

    if(std::string("VsCalo") == JetType1 || std::string("PuCalo") == JetType1){
      mHadEnergy_Jet1_unmatched->Fill((*caloJet1)[Aj.id].hadEnergyInHO()
                      + (*caloJet1)[Aj.id].hadEnergyInHB()
                      + (*caloJet1)[Aj.id].hadEnergyInHF()
                      + (*caloJet1)[Aj.id].hadEnergyInHE()
                      );
      mEmEnergy_Jet1_unmatched->Fill((*caloJet1)[Aj.id].emEnergyInEB()
                     + (*caloJet1)[Aj.id].emEnergyInEE()
                     + (*caloJet1)[Aj.id].emEnergyInHF()
                     );
    }

    if(std::string("VsPF") == JetType1){
      mChargedHadronEnergy_Jet1_unmatched->Fill((*pfJets)[Aj.id].chargedHadronEnergy());
      mNeutralHadronEnergy_Jet1_unmatched->Fill((*pfJets)[Aj.id].neutralHadronEnergy());
      mChargedEmEnergy_Jet1_unmatched->Fill((*pfJets)[Aj.id].chargedEmEnergy());
      mNeutralEmEnergy_Jet1_unmatched->Fill((*pfJets)[Aj.id].neutralEmEnergy());
      mChargedMuEnergy_Jet1_unmatched->Fill((*pfJets)[Aj.id].chargedMuEnergy());

      mChargedHadEnergyFraction_Jet1_unmatched->Fill((*pfJets)[Aj.id].chargedHadronEnergyFraction());
      mNeutralHadEnergyFraction_Jet1_unmatched->Fill((*pfJets)[Aj.id].neutralHadronEnergyFraction());
      mPhotonEnergyFraction_Jet1_unmatched->Fill((*pfJets)[Aj.id].photonEnergyFraction());
      mElectronEnergyFraction_Jet1_unmatched->Fill((*pfJets)[Aj.id].electronEnergyFraction());
      mMuonEnergyFraction_Jet1_unmatched->Fill((*pfJets)[Aj.id].muonEnergyFraction());
    }
    
      }

      if(Jet2_ID[Bj.id] == 0) {

    mpT_Jet2_unmatched->Fill(recoJet2[Bj.id]->pt());
    unmatchedJet2++;
    Jet2_ID[Bj.id] = 2;
    if(std::string("VsCalo") == JetType2 || std::string("PuCalo") == JetType2){
      mHadEnergy_Jet2_unmatched->Fill((*caloJet2)[Bj.id].hadEnergyInHO()
                      + (*caloJet2)[Bj.id].hadEnergyInHB()
                      + (*caloJet2)[Bj.id].hadEnergyInHF()
                      + (*caloJet2)[Bj.id].hadEnergyInHE()
                      );
      mEmEnergy_Jet2_unmatched->Fill((*caloJet2)[Bj.id].emEnergyInEB()
                     + (*caloJet2)[Bj.id].emEnergyInEE()
                     + (*caloJet2)[Bj.id].emEnergyInHF()
                     );
    }
    
    
    if(std::string("VsPF") == JetType2){
      mChargedHadronEnergy_Jet2_unmatched->Fill((*pfJets)[Bj.id].chargedHadronEnergy());
      mNeutralHadronEnergy_Jet2_unmatched->Fill((*pfJets)[Bj.id].neutralHadronEnergy());
      mChargedEmEnergy_Jet2_unmatched->Fill((*pfJets)[Bj.id].chargedEmEnergy());
      mNeutralEmEnergy_Jet2_unmatched->Fill((*pfJets)[Bj.id].neutralEmEnergy());
      mChargedMuEnergy_Jet2_unmatched->Fill((*pfJets)[Bj.id].chargedMuEnergy());
      
      mChargedHadEnergyFraction_Jet2_unmatched->Fill((*pfJets)[Bj.id].chargedHadronEnergyFraction());
      mNeutralHadEnergyFraction_Jet2_unmatched->Fill((*pfJets)[Bj.id].neutralHadronEnergyFraction());
      mPhotonEnergyFraction_Jet2_unmatched->Fill((*pfJets)[Bj.id].photonEnergyFraction());
      mElectronEnergyFraction_Jet2_unmatched->Fill((*pfJets)[Bj.id].electronEnergyFraction());
      mMuonEnergyFraction_Jet2_unmatched->Fill((*pfJets)[Bj.id].muonEnergyFraction());
    }
       
      }
      
    }
    
  }// both Jet1 and Jet2 in the event. 


}

void JetAnalyzer_HeavyIons_matching::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  iRun, edm::EventSetup const &    )

override

Definition at line 90 of file JetAnalyzer_HeavyIons_matching.cc.

References DQMStore::IBooker::book1D(), JetType1, JetType2, edm::InputTag::label(), mChargedEmEnergy_Jet1_unmatched, mChargedEmEnergy_Jet2_unmatched, mChargedHadEnergyFraction_Jet1_unmatched, mChargedHadEnergyFraction_Jet2_unmatched, mChargedHadronEnergy_Jet1_unmatched, mChargedHadronEnergy_Jet2_unmatched, mChargedMuEnergy_Jet1_unmatched, mChargedMuEnergy_Jet2_unmatched, mElectronEnergyFraction_Jet1_unmatched, mElectronEnergyFraction_Jet2_unmatched, mEmEnergy_Jet1_unmatched, mEmEnergy_Jet2_unmatched, mHadEnergy_Jet1_unmatched, mHadEnergy_Jet2_unmatched, mInputJet1Collection, mInputJet2Collection, mMuonEnergyFraction_Jet1_unmatched, mMuonEnergyFraction_Jet2_unmatched, mNeutralEmEnergy_Jet1_unmatched, mNeutralEmEnergy_Jet2_unmatched, mNeutralHadEnergyFraction_Jet1_unmatched, mNeutralHadEnergyFraction_Jet2_unmatched, mNeutralHadronEnergy_Jet1_unmatched, mNeutralHadronEnergy_Jet2_unmatched, mOutputFile, mPhotonEnergyFraction_Jet1_unmatched, mPhotonEnergyFraction_Jet2_unmatched, mpT_Jet1_matched, mpT_Jet1_unmatched, mpT_Jet2_matched, mpT_Jet2_unmatched, mpT_ratio_Jet1Jet2, DQMStore::IBooker::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

  {

    ibooker.setCurrentFolder("JetMET/HIJetValidation/"+mInputJet1Collection.label()+"_DeltaRMatched_"+mInputJet2Collection.label());
    
    mpT_ratio_Jet1Jet2 = ibooker.book1D("Ratio_Jet1pT_vs_Jet2pT",Form(";Matched %s Jet pT / %s Jet pT;Counts", "JetType1", "JetType2"),100, 0, 10);
    mpT_Jet1_matched = ibooker.book1D("Jet1_matched_jet_Spectra",Form(";Matched %s Spectra; counts", "JetType1"),100, 0, 1000);
    mpT_Jet2_matched = ibooker.book1D("Jet2_matched_jet_Spectra",Form(";Matched %s Spectra; counts", "JetType2"),100, 0, 1000);
    mpT_Jet1_unmatched = ibooker.book1D("Jet1_unmatched_jet_Spectra",Form(";Unmatched %s spectra;counts","JetType1"),100, 0, 1000);
    mpT_Jet2_unmatched = ibooker.book1D("Jet2_unmatched_jet_Spectra",Form(";Unmatched %s spectra;counts","JetType2"),100, 0, 1000);
    
    if(std::string("VsCalo") == JetType1 || std::string("PuCalo") == JetType1){
      mHadEnergy_Jet1_unmatched = ibooker.book1D("HadEnergy_Jet1_unmatched", Form("HadEnergy_Jet1_unmatched;HadEnergy unmatched %s;counts", "JetType1"), 50, 0, 200);
      mEmEnergy_Jet1_unmatched = ibooker.book1D("EmEnergy_Jet1_unmatched", Form("EmEnergy_Jet1_unmatched;EMEnergy unmatched %s;counts", "JetType1"), 50, 0, 200);
    }
    
    if(std::string("VsCalo") == JetType2 || std::string("PuCalo") == JetType2){
      mHadEnergy_Jet2_unmatched = ibooker.book1D("HadEnergy_Jet2_unmatched", Form("HadEnergy_Jet2_unmatched;HadEnergy unmatched %s;counts", "JetType2"), 50, 0, 200);
      mEmEnergy_Jet2_unmatched = ibooker.book1D("EmEnergy_Jet2_unmatched", Form("EmEnergy_Jet2_unmatched;EMEnergy unmatched %s;counts", "JetType2"), 50, 0, 200);
    }

    if(std::string("VsPF") == JetType1){
      mChargedHadronEnergy_Jet1_unmatched = ibooker.book1D("ChargedHadronEnergy_Jet1_unmatched", Form(";charged HAD energy unmatched %s;counts","JetType1"),    100, 0, 300);
      mNeutralHadronEnergy_Jet1_unmatched = ibooker.book1D("neutralHadronEnergy_Jet1_unmatched", Form(";neutral HAD energy unmatched %s;counts", "JetType1"),    100, 0, 300);
      mChargedEmEnergy_Jet1_unmatched = ibooker.book1D("ChargedEmEnergy_Jet1_unmatched", Form(";charged EM energy unmatched %s;counts", "JetType1"),    100, 0, 300);
      mNeutralEmEnergy_Jet1_unmatched = ibooker.book1D("neutralEmEnergy_Jet1_unmatched", Form(";neutral EM energy unmatched %s;counts", "JetType1"),    100, 0, 300);
      mChargedMuEnergy_Jet1_unmatched = ibooker.book1D("ChargedMuEnergy_Jet1_unmatched", Form(";charged Mu energy unmatched %s;counts", "JetType1"),    100, 0, 300);

      mChargedHadEnergyFraction_Jet1_unmatched = ibooker.book1D("ChargedHadEnergyFraction_Jet1_unmatched",Form(";h^{+/-} Energy Fraction %s;counts", "JetType1"),50, 0, 1);
      mNeutralHadEnergyFraction_Jet1_unmatched = ibooker.book1D("NeutralHadEnergyFraction_Jet1_unmatched",Form(";h^{0} Energy Fraction %s;counts", "JetType1"),50, 0, 1);
      mPhotonEnergyFraction_Jet1_unmatched = ibooker.book1D("PhotonEnergyFraction_Jet1_unmatched",Form(";#gamma Energy Fraction %s;counts", "JetType1"),50, 0, 1);
      mElectronEnergyFraction_Jet1_unmatched = ibooker.book1D("ElectronEnergyFraction_Jet1_unmatched",Form(";e Energy Fraction %s;counts", "JetType1"),50, 0, 1);
      mMuonEnergyFraction_Jet1_unmatched = ibooker.book1D("MuonoEnergyFraction_Jet1_unmatched",Form(";#mu Energy Fraction %s;counts", "JetType1"),50, 0, 1);
      
    }

    if(std::string("VsPF") == JetType2){
      mChargedHadronEnergy_Jet2_unmatched = ibooker.book1D("ChargedHadronEnergy_Jet2_unmatched", Form(";charged HAD energy unmatched %s;counts","JetType2"),    100, 0, 300);
      mNeutralHadronEnergy_Jet2_unmatched = ibooker.book1D("neutralHadronEnergy_Jet2_unmatched", Form(";neutral HAD energy unmatched %s;counts", "JetType2"),    100, 0, 300);
      mChargedEmEnergy_Jet2_unmatched = ibooker.book1D("ChargedEmEnergy_Jet2_unmatched", Form(";charged EM energy unmatched %s;counts", "JetType2"),    100, 0, 300);
      mNeutralEmEnergy_Jet2_unmatched = ibooker.book1D("neutralEmEnergy_Jet2_unmatched", Form(";neutral EM energy unmatched %s;counts", "JetType2"),    100, 0, 300);
      mChargedMuEnergy_Jet2_unmatched = ibooker.book1D("ChargedMuEnergy_Jet2_unmatched", Form(";charged Mu energy unmatched %s;counts", "JetType2"),    100, 0, 300);

      mChargedHadEnergyFraction_Jet2_unmatched = ibooker.book1D("ChargedHadEnergyFraction_Jet2_unmatched",Form(";h^{+/-} Energy Fraction %s;counts", "JetType2"),50, 0, 1);
      mNeutralHadEnergyFraction_Jet2_unmatched = ibooker.book1D("NeutralHadEnergyFraction_Jet2_unmatched",Form(";h^{0} Energy Fraction %s;counts", "JetType2"),50, 0, 1);
      mPhotonEnergyFraction_Jet2_unmatched = ibooker.book1D("PhotonEnergyFraction_Jet2_unmatched",Form(";#gamma Energy Fraction %s;counts", "JetType2"),50, 0, 1);
      mElectronEnergyFraction_Jet2_unmatched = ibooker.book1D("ElectronEnergyFraction_Jet2_unmatched",Form(";e Energy Fraction %s;counts", "JetType2"),50, 0, 1);
      mMuonEnergyFraction_Jet2_unmatched = ibooker.book1D("MuonoEnergyFraction_Jet2_unmatched",Form(";#mu Energy Fraction %s;counts", "JetType2"),50, 0, 1);

    }
    
    if (mOutputFile.empty ()) 
      LogInfo("OutputInfo") << " Histograms will NOT be saved";
    else 
      LogInfo("OutputInfo") << " Histograms will be saved to file:" << mOutputFile;


  }

static float JetAnalyzer_HeavyIons_matching::deltaRR ( float  eta1, float  phi1, float  eta2, float  phi2  )

inlinestaticprivate

Definition at line 178 of file JetAnalyzer_HeavyIons_matching.h.

References runTauDisplay::dr, M_PI, and mathSSE::sqrt().

Referenced by analyze(), and JetAnalyzer_HeavyIons_matching::CompareMatchedJets::operator()().

  {
    float deta = eta1 - eta2;
    float dphi = fabs(phi1 - phi2);
    if(dphi > M_PI )dphi -= 2*M_PI;
    float dr = sqrt(deta*deta + dphi*dphi);
    return dr;
  }

Member Data Documentation

edm::EDGetTokenT<reco::BasicJetCollection> JetAnalyzer_HeavyIons_matching::basicJetsToken_

private

Definition at line 111 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), and JetAnalyzer_HeavyIons_matching().

edm::EDGetTokenT<reco::CaloJetCollection> JetAnalyzer_HeavyIons_matching::caloJet1Token_

private

Definition at line 108 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), and JetAnalyzer_HeavyIons_matching().

edm::EDGetTokenT<reco::CaloJetCollection> JetAnalyzer_HeavyIons_matching::caloJet2Token_

private

Definition at line 109 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), and JetAnalyzer_HeavyIons_matching().

std::string JetAnalyzer_HeavyIons_matching::JetType1

private

Definition at line 101 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

std::string JetAnalyzer_HeavyIons_matching::JetType2

private

Definition at line 102 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

edm::EDGetTokenT<reco::JPTJetCollection> JetAnalyzer_HeavyIons_matching::jptJetsToken_

private

Definition at line 112 of file JetAnalyzer_HeavyIons_matching.h.

MonitorElement* JetAnalyzer_HeavyIons_matching::mChargedEmEnergy_Jet1_unmatched

private

Definition at line 124 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mChargedEmEnergy_Jet2_unmatched

private

Definition at line 138 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mChargedHadEnergyFraction_Jet1_unmatched

private

Definition at line 128 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mChargedHadEnergyFraction_Jet2_unmatched

private

Definition at line 142 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mChargedHadronEnergy_Jet1_unmatched

private

Definition at line 122 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mChargedHadronEnergy_Jet2_unmatched

private

Definition at line 136 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mChargedMuEnergy_Jet1_unmatched

private

Definition at line 126 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mChargedMuEnergy_Jet2_unmatched

private

Definition at line 140 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mElectronEnergyFraction_Jet1_unmatched

private

Definition at line 131 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mElectronEnergyFraction_Jet2_unmatched

private

Definition at line 145 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mEmEnergy_Jet1_unmatched

private

Definition at line 121 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mEmEnergy_Jet2_unmatched

private

Definition at line 135 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mHadEnergy_Jet1_unmatched

private

Definition at line 120 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mHadEnergy_Jet2_unmatched

private

Definition at line 134 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

edm::InputTag JetAnalyzer_HeavyIons_matching::mInputJet1Collection

private

Definition at line 97 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by bookHistograms(), and JetAnalyzer_HeavyIons_matching().

edm::InputTag JetAnalyzer_HeavyIons_matching::mInputJet2Collection

private

Definition at line 98 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mMuonEnergyFraction_Jet1_unmatched

private

Definition at line 132 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mMuonEnergyFraction_Jet2_unmatched

private

Definition at line 146 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mNeutralEmEnergy_Jet1_unmatched

private

Definition at line 125 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mNeutralEmEnergy_Jet2_unmatched

private

Definition at line 139 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mNeutralHadEnergyFraction_Jet1_unmatched

private

Definition at line 129 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mNeutralHadEnergyFraction_Jet2_unmatched

private

Definition at line 143 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mNeutralHadronEnergy_Jet1_unmatched

private

Definition at line 123 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mNeutralHadronEnergy_Jet2_unmatched

private

Definition at line 137 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

std::string JetAnalyzer_HeavyIons_matching::mOutputFile

private

Definition at line 100 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by bookHistograms().

MonitorElement* JetAnalyzer_HeavyIons_matching::mPhotonEnergyFraction_Jet1_unmatched

private

Definition at line 130 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mPhotonEnergyFraction_Jet2_unmatched

private

Definition at line 144 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mpT_Jet1_matched

private

Definition at line 115 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mpT_Jet1_unmatched

private

Definition at line 117 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mpT_Jet2_matched

private

Definition at line 116 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mpT_Jet2_unmatched

private

Definition at line 118 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

MonitorElement* JetAnalyzer_HeavyIons_matching::mpT_ratio_Jet1Jet2

private

Definition at line 114 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), bookHistograms(), and JetAnalyzer_HeavyIons_matching().

double JetAnalyzer_HeavyIons_matching::mRecoDelRMatch

private

Definition at line 104 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze().

double JetAnalyzer_HeavyIons_matching::mRecoJetEtaCut

private

Definition at line 105 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze().

double JetAnalyzer_HeavyIons_matching::mRecoJetPtThreshold

private

Definition at line 103 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze().

edm::EDGetTokenT<reco::PFJetCollection> JetAnalyzer_HeavyIons_matching::pfJetsToken_

private

Definition at line 110 of file JetAnalyzer_HeavyIons_matching.h.

Referenced by analyze(), and JetAnalyzer_HeavyIons_matching().