HLTBJet Class Reference

#include <HLTBJet.h>

Public Member Functions
void	analyze (const edm::Handle< edm::View< reco::Jet > > &rawBJets, const edm::Handle< edm::View< reco::Jet > > &correctedBJets, const edm::Handle< reco::JetTagCollection > &lifetimeBJetsL25, const edm::Handle< reco::JetTagCollection > &lifetimeBJetsL3, const edm::Handle< reco::JetTagCollection > &softmuonBJetsL25, const edm::Handle< reco::JetTagCollection > &softmuonBJetsL3, const edm::Handle< reco::JetTagCollection > &performanceBJetsL25, const edm::Handle< reco::JetTagCollection > &performanceBJetsL3, TTree *tree)
void	clear (void)
	HLTBJet ()
void	setup (const edm::ParameterSet &config, TTree *tree)
	~HLTBJet ()

Private Member Functions
void	analyseCorrectedJets (const edm::View< reco::Jet > &jets)
void	analyseJets (const edm::View< reco::Jet > &jets)
void	analyseLifetime (const edm::View< reco::Jet > &jets, const reco::JetTagCollection &tagsL25, const reco::JetTagCollection &tagsL3)
void	analysePerformance (const edm::View< reco::Jet > &jets, const reco::JetTagCollection &tagsL25, const reco::JetTagCollection &tagsL3)
void	analyseSoftmuon (const edm::View< reco::Jet > &jets, const reco::JetTagCollection &tagsL25, const reco::JetTagCollection &tagsL3)

Private Attributes
int	NohBJetL2
int	NohBJetL2Corrected
float *	ohBJetIPL25Tag
float *	ohBJetIPL3Tag
float *	ohBJetL2CorrectedEnergy
float *	ohBJetL2CorrectedEt
float *	ohBJetL2CorrectedEta
float *	ohBJetL2CorrectedPhi
float *	ohBJetL2CorrectedPt
float *	ohBJetL2Energy
float *	ohBJetL2Et
float *	ohBJetL2Eta
float *	ohBJetL2Phi
float *	ohBJetL2Pt
int *	ohBJetMuL25Tag
float *	ohBJetMuL3Tag
int *	ohBJetPerfL25Tag
int *	ohBJetPerfL3Tag

Detailed Description

Definition at line 14 of file HLTBJet.h.

Constructor & Destructor Documentation

HLTBJet::HLTBJet

(

)

Definition at line 24 of file HLTBJet.cc.

References kMaxBJets, NohBJetL2, NohBJetL2Corrected, ohBJetIPL25Tag, ohBJetIPL3Tag, ohBJetL2CorrectedEnergy, ohBJetL2CorrectedEt, ohBJetL2CorrectedEta, ohBJetL2CorrectedPhi, ohBJetL2CorrectedPt, ohBJetL2Energy, ohBJetL2Et, ohBJetL2Eta, ohBJetL2Phi, ohBJetL2Pt, ohBJetMuL25Tag, ohBJetMuL3Tag, ohBJetPerfL25Tag, and ohBJetPerfL3Tag.

{
  // set of variables for uncorrected L2 jets
  NohBJetL2                 = 0;
  ohBJetL2Energy            = new float[kMaxBJets];
  ohBJetL2Et                = new float[kMaxBJets];
  ohBJetL2Pt                = new float[kMaxBJets];
  ohBJetL2Eta               = new float[kMaxBJets];
  ohBJetL2Phi               = new float[kMaxBJets];
              
  // set of variables for corrected L2 jets
  NohBJetL2Corrected        = 0;
  ohBJetL2CorrectedEnergy   = new float[kMaxBJets];
  ohBJetL2CorrectedEt       = new float[kMaxBJets];
  ohBJetL2CorrectedPt       = new float[kMaxBJets];
  ohBJetL2CorrectedEta      = new float[kMaxBJets];
  ohBJetL2CorrectedPhi      = new float[kMaxBJets];
  
  // set of variables for lifetime-based b-tag
  ohBJetIPL25Tag            = new float[kMaxBJets];
  ohBJetIPL3Tag             = new float[kMaxBJets];
  
  // set of variables for soft-muon-based b-tag
  ohBJetMuL25Tag            = new int[kMaxBJets];          // do not optimize
  ohBJetMuL3Tag             = new float[kMaxBJets];        // SoftMuonbyPt
  
  // set of variables for b-tagging performance measurements
  // SoftMuonbyDR
  ohBJetPerfL25Tag          = new int[kMaxBJets];          // do not optimize 
  ohBJetPerfL3Tag           = new int[kMaxBJets];          // do not optimize
}

HLTBJet::~HLTBJet

(

)

Definition at line 81 of file HLTBJet.cc.

{ }

Member Function Documentation

void HLTBJet::analyseCorrectedJets ( const edm::View< reco::Jet > &  jets )

private

Definition at line 159 of file HLTBJet.cc.

References i, kMaxBJets, min, NohBJetL2Corrected, ohBJetL2CorrectedEnergy, ohBJetL2CorrectedEt, ohBJetL2CorrectedEta, ohBJetL2CorrectedPhi, ohBJetL2CorrectedPt, edm::View< T >::size(), and findQualityFiles::size.

Referenced by analyze().

{
  // the jets need to be persistable, so .size() returns an 'unsigned int' to be stable across the architectures
  // so, for the comparison, we cast back to size_t
  size_t size = std::min(kMaxBJets, size_t(jets.size()) );
  NohBJetL2Corrected = size;
  for (size_t i = 0; i < size; ++i) {
    ohBJetL2CorrectedEnergy[i] = jets[i].energy();
    ohBJetL2CorrectedEt[i]     = jets[i].et();
    ohBJetL2CorrectedPt[i]     = jets[i].pt();
    ohBJetL2CorrectedEta[i]    = jets[i].eta();
    ohBJetL2CorrectedPhi[i]    = jets[i].phi();
  }
}

void HLTBJet::analyseJets ( const edm::View< reco::Jet > &  jets )

private

Definition at line 144 of file HLTBJet.cc.

References i, kMaxBJets, min, NohBJetL2, ohBJetL2Energy, ohBJetL2Et, ohBJetL2Eta, ohBJetL2Phi, ohBJetL2Pt, edm::View< T >::size(), and findQualityFiles::size.

Referenced by analyze().

{
  // the jets need to be persistable, so .size() returns an 'unsigned int' to be stable across the architectures
  // so, for the comparison, we cast back to size_t
  size_t size = std::min(kMaxBJets, size_t(jets.size()) ); 
  NohBJetL2 = size;
  for (size_t i = 0; i < size; ++i) {
    ohBJetL2Energy[i] = jets[i].energy();
    ohBJetL2Et[i]     = jets[i].et();
    ohBJetL2Pt[i]     = jets[i].pt();
    ohBJetL2Eta[i]    = jets[i].eta();
    ohBJetL2Phi[i]    = jets[i].phi();
  }
}

void HLTBJet::analyseLifetime ( const edm::View< reco::Jet > &  jets, const reco::JetTagCollection &  tagsL25, const reco::JetTagCollection &  tagsL3  )

private

Definition at line 174 of file HLTBJet.cc.

References i, kBTagLifetimeBJetsL25, kBTagLifetimeBJetsL3, kMaxBJets, min, ohBJetIPL25Tag, ohBJetIPL3Tag, edm::AssociationVector< KeyRefProd, CVal, KeyRef, SizeType, KeyReferenceHelper >::size(), edm::View< T >::size(), and findQualityFiles::size.

Referenced by analyze().

{
  if (tagsL25.size() != jets.size()) {
    edm::LogWarning("OpenHLT") << kBTagLifetimeBJetsL25 << " collection has " << tagsL25.size() << " elements, but " << jets.size() << " where expected from L2" << std::endl;
    return;
  }
  if (tagsL3.size() != jets.size()) {
    edm::LogWarning("OpenHLT") << kBTagLifetimeBJetsL3 << " collection has " << tagsL3.size() << " elements, but " << jets.size() << " where expected from L2" << std::endl;
    return;
  }
  // the jets need to be persistable, so .size() returns an 'unsigned int' to be stable across the architectures
  // so, for the comparison, we cast back to size_t
  size_t size = std::min(kMaxBJets, size_t(jets.size()) );
  for (size_t i = 0; i < size; i++) {
    ohBJetIPL25Tag[i] = tagsL25[i].second;
    ohBJetIPL3Tag[i]  = tagsL3[i].second;
  }
}

void HLTBJet::analysePerformance ( const edm::View< reco::Jet > &  jets, const reco::JetTagCollection &  tagsL25, const reco::JetTagCollection &  tagsL3  )

private

Definition at line 218 of file HLTBJet.cc.

References i, kBTagPerformanceBJetsL25, kBTagPerformanceBJetsL3, kMaxBJets, min, ohBJetPerfL25Tag, ohBJetPerfL3Tag, edm::AssociationVector< KeyRefProd, CVal, KeyRef, SizeType, KeyReferenceHelper >::size(), edm::View< T >::size(), and findQualityFiles::size.

Referenced by analyze().

{
  if (tagsL25.size() != jets.size()) {
    edm::LogWarning("OpenHLT") << kBTagPerformanceBJetsL25 << " collection has " << tagsL25.size() << " elements, but " << jets.size() << " where expected from L2" << std::endl;
    return;
  }
  if (tagsL3.size() != jets.size()) {
    edm::LogWarning("OpenHLT") << kBTagPerformanceBJetsL3 << " collection has " << tagsL3.size() << " elements, but " << jets.size() << " where expected from L2" << std::endl;
    return;
  }
  // the jets need to be persistable, so .size() returns an 'unsigned int' to be stable across the architectures
  // so, for the comparison, we cast back to size_t
  size_t size = std::min(kMaxBJets, size_t(jets.size()) );
  for (size_t i = 0; i < size; i++) {
    ohBJetPerfL25Tag[i] = (tagsL25[i].second > 0.) ? 1 : 0;
    ohBJetPerfL3Tag[i]  = (tagsL3[i].second  > 0.) ? 1 : 0;
  }
}

void HLTBJet::analyseSoftmuon ( const edm::View< reco::Jet > &  jets, const reco::JetTagCollection &  tagsL25, const reco::JetTagCollection &  tagsL3  )

private

Definition at line 196 of file HLTBJet.cc.

References i, kBTagSoftmuonBJetsL25, kBTagSoftmuonBJetsL3, kMaxBJets, min, ohBJetMuL25Tag, ohBJetMuL3Tag, edm::AssociationVector< KeyRefProd, CVal, KeyRef, SizeType, KeyReferenceHelper >::size(), edm::View< T >::size(), and findQualityFiles::size.

Referenced by analyze().

{
  if (tagsL25.size() != jets.size()) {
    edm::LogWarning("OpenHLT") << kBTagSoftmuonBJetsL25 << " collection has " << tagsL25.size() << " elements, but " << jets.size() << " where expected from L2" << std::endl;
    return;
  }
  if (tagsL3.size() != jets.size()) {
    edm::LogWarning("OpenHLT") << kBTagSoftmuonBJetsL3 << " collection has " << tagsL3.size() << " elements, but " << jets.size() << " where expected from L2" << std::endl;
    return;
  }
  // the jets need to be persistable, so .size() returns an 'unsigned int' to be stable across the architectures
  // so, for the comparison, we cast back to size_t
  size_t size = std::min(kMaxBJets, size_t(jets.size()) );
  for (size_t i = 0; i < size; i++) {
    ohBJetMuL25Tag[i] = (tagsL25[i].second > 0.) ? 1 : 0;
    ohBJetMuL3Tag[i]  = tagsL3[i].second;
  }
}

void HLTBJet::analyze	(	const edm::Handle< edm::View< reco::Jet > > &	rawBJets,
		const edm::Handle< edm::View< reco::Jet > > &	correctedBJets,
		const edm::Handle< reco::JetTagCollection > &	lifetimeBJetsL25,
		const edm::Handle< reco::JetTagCollection > &	lifetimeBJetsL3,
		const edm::Handle< reco::JetTagCollection > &	softmuonBJetsL25,
		const edm::Handle< reco::JetTagCollection > &	softmuonBJetsL3,
		const edm::Handle< reco::JetTagCollection > &	performanceBJetsL25,
		const edm::Handle< reco::JetTagCollection > &	performanceBJetsL3,
		TTree *	tree
	)

Definition at line 112 of file HLTBJet.cc.

References analyseCorrectedJets(), analyseJets(), analyseLifetime(), analysePerformance(), analyseSoftmuon(), clear(), and edm::HandleBase::isValid().

Referenced by HLTAnalyzer::analyze().

{
  // reset the tree variables
  clear();
  
  // if the required collections are available, fill the corresponding tree branches
  if (rawBJets.isValid())
    analyseJets(* rawBJets);

  if (correctedBJets.isValid())
    analyseCorrectedJets(* correctedBJets);
 
  if (correctedBJets.isValid() and lifetimeBJetsL25.isValid() and lifetimeBJetsL3.isValid())
    analyseLifetime(* correctedBJets, * lifetimeBJetsL25, * lifetimeBJetsL3);

  if (correctedBJets.isValid() and softmuonBJetsL25.isValid() and softmuonBJetsL3.isValid())
    analyseSoftmuon(* correctedBJets, * softmuonBJetsL25, * softmuonBJetsL3);
  
  if (correctedBJets.isValid() and performanceBJetsL25.isValid() and performanceBJetsL3.isValid())
    analysePerformance(* correctedBJets, * performanceBJetsL25, * performanceBJetsL3);

}

void HLTBJet::clear

(

void

)

Definition at line 56 of file HLTBJet.cc.

References kMaxBJets, NohBJetL2, NohBJetL2Corrected, ohBJetIPL25Tag, ohBJetIPL3Tag, ohBJetL2CorrectedEnergy, ohBJetL2CorrectedEt, ohBJetL2CorrectedEta, ohBJetL2CorrectedPhi, ohBJetL2CorrectedPt, ohBJetL2Energy, ohBJetL2Et, ohBJetL2Eta, ohBJetL2Phi, ohBJetL2Pt, ohBJetMuL25Tag, ohBJetMuL3Tag, ohBJetPerfL25Tag, and ohBJetPerfL3Tag.

Referenced by analyze(), python.Vispa.Views.BoxDecayView.BoxDecayView::closeEvent(), python.Vispa.Views.LineDecayView.LineDecayView::setDataObjects(), and python.Vispa.Views.BoxDecayView.BoxDecayView::updateContent().

{
  NohBJetL2          = 0;
  NohBJetL2Corrected = 0;
  std::memset(ohBJetL2Energy,            '\0', kMaxBJets * sizeof(float));
  std::memset(ohBJetL2Et,                '\0', kMaxBJets * sizeof(float));
  std::memset(ohBJetL2Et,                '\0', kMaxBJets * sizeof(float));
  std::memset(ohBJetL2Pt,                '\0', kMaxBJets * sizeof(float));
  std::memset(ohBJetL2Eta,               '\0', kMaxBJets * sizeof(float));
  std::memset(ohBJetL2Phi,               '\0', kMaxBJets * sizeof(float));
  std::memset(ohBJetL2CorrectedEnergy,   '\0', kMaxBJets * sizeof(float));
  std::memset(ohBJetL2CorrectedEt,       '\0', kMaxBJets * sizeof(float));
  std::memset(ohBJetL2CorrectedPt,       '\0', kMaxBJets * sizeof(float));
  std::memset(ohBJetL2CorrectedEta,      '\0', kMaxBJets * sizeof(float));
  std::memset(ohBJetL2CorrectedPhi,      '\0', kMaxBJets * sizeof(float));

  std::memset(ohBJetIPL25Tag,            '\0', kMaxBJets * sizeof(float));
  std::memset(ohBJetIPL3Tag,             '\0', kMaxBJets * sizeof(float));

  std::memset(ohBJetMuL25Tag,            '\0', kMaxBJets * sizeof(int));
  std::memset(ohBJetMuL3Tag,             '\0', kMaxBJets * sizeof(float));
  std::memset(ohBJetPerfL25Tag,          '\0', kMaxBJets * sizeof(int));
  std::memset(ohBJetPerfL3Tag,           '\0', kMaxBJets * sizeof(int));
}

void HLTBJet::setup	(	const edm::ParameterSet &	config,
		TTree *	tree
	)

Definition at line 84 of file HLTBJet.cc.

References NohBJetL2, NohBJetL2Corrected, ohBJetIPL25Tag, ohBJetIPL3Tag, ohBJetL2CorrectedEnergy, ohBJetL2CorrectedEt, ohBJetL2CorrectedEta, ohBJetL2CorrectedPhi, ohBJetL2CorrectedPt, ohBJetL2Energy, ohBJetL2Et, ohBJetL2Eta, ohBJetL2Phi, ohBJetL2Pt, ohBJetMuL25Tag, ohBJetMuL3Tag, ohBJetPerfL25Tag, and ohBJetPerfL3Tag.

Referenced by HLTAnalyzer::HLTAnalyzer().

{
  // create the TTree branches
  if (tree) {
    tree->Branch("NohBJetL2",               & NohBJetL2,             "NohBJetL2/I");
    tree->Branch("ohBJetL2Energy",          ohBJetL2Energy,          "ohBJetL2Energy[NohBJetL2]/F");
    tree->Branch("ohBJetL2Et",              ohBJetL2Et,              "ohBJetL2Et[NohBJetL2]/F");
    tree->Branch("ohBJetL2Pt",              ohBJetL2Pt,              "ohBJetL2Pt[NohBJetL2]/F");
    tree->Branch("ohBJetL2Eta",             ohBJetL2Eta,             "ohBJetL2Eta[NohBJetL2]/F");
    tree->Branch("ohBJetL2Phi",             ohBJetL2Phi,             "ohBJetL2Phi[NohBJetL2]/F");
                
    tree->Branch("NohBJetL2Corrected",      & NohBJetL2Corrected,    "NohBJetL2Corrected/I");
    tree->Branch("ohBJetL2CorrectedEnergy", ohBJetL2CorrectedEnergy, "ohBJetL2CorrectedEnergy[NohBJetL2Corrected]/F");
    tree->Branch("ohBJetL2CorrectedEt",     ohBJetL2CorrectedEt,     "ohBJetL2CorrectedEt[NohBJetL2Corrected]/F");
    tree->Branch("ohBJetL2CorrectedPt",     ohBJetL2CorrectedPt,     "ohBJetL2CorrectedPt[NohBJetL2Corrected]/F");
    tree->Branch("ohBJetL2CorrectedEta",    ohBJetL2CorrectedEta,    "ohBJetL2CorrectedEta[NohBJetL2Corrected]/F");
    tree->Branch("ohBJetL2CorrectedPhi",    ohBJetL2CorrectedPhi,    "ohBJetL2CorrectedPhi[NohBJetL2Corrected]/F");
    
    tree->Branch("ohBJetIPL25Tag",          ohBJetIPL25Tag,          "ohBJetIPL25Tag[NohBJetL2]/F");
    tree->Branch("ohBJetIPL3Tag",           ohBJetIPL3Tag,           "ohBJetIPL3Tag[NohBJetL2]/F");

    tree->Branch("ohBJetMuL25Tag",          ohBJetMuL25Tag,          "ohBJetMuL25Tag[NohBJetL2]/I");
    tree->Branch("ohBJetMuL3Tag",           ohBJetMuL3Tag,           "ohBJetMuL3Tag[NohBJetL2]/F");
    tree->Branch("ohBJetPerfL25Tag",        ohBJetPerfL25Tag,        "ohBJetPerfL25Tag[NohBJetL2]/I");
    tree->Branch("ohBJetPerfL3Tag",         ohBJetPerfL3Tag,         "ohBJetPerfL3Tag[NohBJetL2]/I");
  }
}

Member Data Documentation

int HLTBJet::NohBJetL2

private

Definition at line 54 of file HLTBJet.h.

Referenced by analyseJets(), clear(), HLTBJet(), and setup().

int HLTBJet::NohBJetL2Corrected

private

Definition at line 62 of file HLTBJet.h.

Referenced by analyseCorrectedJets(), clear(), HLTBJet(), and setup().

float* HLTBJet::ohBJetIPL25Tag

private

Definition at line 70 of file HLTBJet.h.

Referenced by analyseLifetime(), clear(), HLTBJet(), and setup().

float* HLTBJet::ohBJetIPL3Tag

private

Definition at line 71 of file HLTBJet.h.

Referenced by analyseLifetime(), clear(), HLTBJet(), and setup().

float* HLTBJet::ohBJetL2CorrectedEnergy

private

Definition at line 63 of file HLTBJet.h.

Referenced by analyseCorrectedJets(), clear(), HLTBJet(), and setup().

float* HLTBJet::ohBJetL2CorrectedEt

private

Definition at line 64 of file HLTBJet.h.

Referenced by analyseCorrectedJets(), clear(), HLTBJet(), and setup().

float* HLTBJet::ohBJetL2CorrectedEta

private

Definition at line 66 of file HLTBJet.h.

Referenced by analyseCorrectedJets(), clear(), HLTBJet(), and setup().

float* HLTBJet::ohBJetL2CorrectedPhi

private

Definition at line 67 of file HLTBJet.h.

Referenced by analyseCorrectedJets(), clear(), HLTBJet(), and setup().

float* HLTBJet::ohBJetL2CorrectedPt

private

Definition at line 65 of file HLTBJet.h.

Referenced by analyseCorrectedJets(), clear(), HLTBJet(), and setup().

float* HLTBJet::ohBJetL2Energy

private

Definition at line 55 of file HLTBJet.h.

Referenced by analyseJets(), clear(), HLTBJet(), and setup().

float* HLTBJet::ohBJetL2Et

private

Definition at line 56 of file HLTBJet.h.

Referenced by analyseJets(), clear(), HLTBJet(), and setup().

float* HLTBJet::ohBJetL2Eta

private

Definition at line 58 of file HLTBJet.h.

Referenced by analyseJets(), clear(), HLTBJet(), and setup().

float* HLTBJet::ohBJetL2Phi

private

Definition at line 59 of file HLTBJet.h.

Referenced by analyseJets(), clear(), HLTBJet(), and setup().

float* HLTBJet::ohBJetL2Pt

private

Definition at line 57 of file HLTBJet.h.

Referenced by analyseJets(), clear(), HLTBJet(), and setup().

int* HLTBJet::ohBJetMuL25Tag

private

Definition at line 74 of file HLTBJet.h.

Referenced by analyseSoftmuon(), clear(), HLTBJet(), and setup().

float* HLTBJet::ohBJetMuL3Tag

private

Definition at line 75 of file HLTBJet.h.

Referenced by analyseSoftmuon(), clear(), HLTBJet(), and setup().

int* HLTBJet::ohBJetPerfL25Tag

private

Definition at line 78 of file HLTBJet.h.

Referenced by analysePerformance(), clear(), HLTBJet(), and setup().

int* HLTBJet::ohBJetPerfL3Tag

private

Definition at line 79 of file HLTBJet.h.

Referenced by analysePerformance(), clear(), HLTBJet(), and setup().