semilCorr Class Reference

#include <semilCorr.h>

Collaboration diagram for semilCorr:

Public Types
enum	Systematics {   TAGGINGWEIGHT = 1, FRAGMENTATION = 2, DECAY = 3, MSRESO = 4,   IDRESO = 5, DECAYREW = 6, MUONSPECTRUM = 7, ALL = 8 }

Public Member Functions
	semilCorr (const TString &fIn, const std::string &suffix="", bool DebugIn=false)
	~semilCorr ()
float	getSemilCorrToIncl (const TLorentzVector &jet, const TLorentzVector &mu)
float	getBjetCorrToIncl (const TLorentzVector &jet, const TLorentzVector &mu)
float	getSemilCorrToInclSyst (const TLorentzVector &jet, const TLorentzVector &mu, semilCorr::Systematics syst=semilCorr::ALL)

Private Member Functions
float	getResponse (float pt, float eta, const std::vector< TH1F * > &h)
float	getSemilCorrToIncl (const TLorentzVector &jet, const TLorentzVector &mu, const std::vector< TH1F * > &histos)
std::vector< int >	getHistoIndices (semilCorr::Systematics syst)

Private Attributes
std::vector< std::vector< TH1F * > >	m_histos
std::vector< float >	m_etas
TFile *	m_f
bool	m_Debug

Detailed Description

Definition at line 15 of file semilCorr.h.

Member Enumeration Documentation

Systematics

enum semilCorr::Systematics

Enumerator
TAGGINGWEIGHT
FRAGMENTATION
DECAY
MSRESO
IDRESO
DECAYREW
MUONSPECTRUM
ALL

Definition at line 18 of file semilCorr.h.

                  {
    TAGGINGWEIGHT = 1,
    FRAGMENTATION = 2,
    DECAY = 3,
    MSRESO = 4,
    IDRESO = 5,
    DECAYREW = 6,
    MUONSPECTRUM = 7,
    ALL = 8
  };

Constructor & Destructor Documentation

semilCorr()

semilCorr::semilCorr	(	const TString &	fIn,
		const std::string &	suffix = "",
		bool	DebugIn = false
	)

Definition at line 9 of file semilCorr.cxx.

                                                                     {
  m_Debug = DebugIn;
  m_f = TFile::Open(fIn);
  m_etas.push_back(0);
  m_etas.push_back(0.8);
  m_etas.push_back(1.2);
  m_etas.push_back(1.7);
  m_etas.push_back(2.1);
  m_etas.push_back(2.5);
 
  vector<string> etastr;
  etastr.push_back("_e0");
  etastr.push_back("_e1");
  etastr.push_back("_e2");
  etastr.push_back("_e3");
  etastr.push_back("_e4");
 
  vector<string> prefix;
  prefix.push_back("corr");
  prefix.push_back("tagSyst");
  prefix.push_back("fragSyst");
  prefix.push_back("decaySyst");
  prefix.push_back("msSyst");
  prefix.push_back("idSyst");
  prefix.push_back("decayRewSyst");
  prefix.push_back("muRewSyst");
  prefix.push_back("corrIncl");
 
  for(unsigned int j = 0; j<prefix.size(); j++){
    vector<TH1F*> corr;
    for(unsigned int i = 0; i<etastr.size(); i++){
      corr.push_back((TH1F*) m_f->Get((prefix[j]+etastr[i]).c_str())); 
    }
    m_histos.push_back(corr);
  }
}

~semilCorr()

semilCorr::~semilCorr

(

)

Definition at line 47 of file semilCorr.cxx.

                     {
  for(unsigned int i = 0; i<m_histos.size();i++){
    for(unsigned int j = 0; j<m_histos[i].size();j++){
      delete m_histos[i][j]; // TH2F's that store response info
    }
  }
  m_f->Close();
  delete m_f;
}

Member Function Documentation

getBjetCorrToIncl()

float semilCorr::getBjetCorrToIncl	(	const TLorentzVector &	jet,
		const TLorentzVector &	mu
	)

Definition at line 63 of file semilCorr.cxx.

{
  return getSemilCorrToIncl(jet,mu,m_histos[8]);
}

getHistoIndices()

vector< int > semilCorr::getHistoIndices ( semilCorr::Systematics  syst )

private

Definition at line 79 of file semilCorr.cxx.

{
  vector<int> indices;
  if(syst == semilCorr::ALL){
    indices.push_back(((int)semilCorr::TAGGINGWEIGHT));
    indices.push_back(((int)semilCorr::FRAGMENTATION));
    indices.push_back(((int)semilCorr::DECAY));
    indices.push_back(((int)semilCorr::MSRESO));
    indices.push_back(((int)semilCorr::IDRESO));
    indices.push_back(((int)semilCorr::DECAYREW));
    indices.push_back(((int)semilCorr::MUONSPECTRUM));
  }else
    indices.push_back((int) syst);
  return indices;
}

getResponse()

float semilCorr::getResponse ( float  pt, float  eta, const std::vector< TH1F * > &  h  )

private

Definition at line 108 of file semilCorr.cxx.

{
  float usePt = pt;
  int histbin = -1;
  for(unsigned int i = 0; i<m_etas.size()-1; i++){
    if(fabs(eta)>=m_etas[i] && fabs(eta)<m_etas[i+1])
      histbin = i;
  }
  
  if(histbin == -1) histbin = h.size()-1;
  float min = h[histbin]->GetBinCenter(1);
  float max = h[histbin]->GetBinCenter(h[histbin]->GetNbinsX());
  if(pt>max)
    usePt = max-0.001;
  if(pt<min) usePt = min+0.001;
  float result = h[histbin]->Interpolate(usePt);
  return result;
}

getSemilCorrToIncl() [1/2]

float semilCorr::getSemilCorrToIncl	(	const TLorentzVector &	jet,
		const TLorentzVector &	mu
	)

Definition at line 57 of file semilCorr.cxx.

{
  return getSemilCorrToIncl(jet,mu,m_histos[0]);
}

getSemilCorrToIncl() [2/2]

float semilCorr::getSemilCorrToIncl ( const TLorentzVector &  jet, const TLorentzVector &  mu, const std::vector< TH1F * > &  histos  )

private

Definition at line 69 of file semilCorr.cxx.

{
  TLorentzVector jetmu = jet+mu;
  //correction to get things to 1 (or to pttruth), not to reference
  float corr = getResponse(jetmu.Pt(), jetmu.Eta(), histos);  
  return corr;
}

getSemilCorrToInclSyst()

float semilCorr::getSemilCorrToInclSyst	(	const TLorentzVector &	jet,
		const TLorentzVector &	mu,
		semilCorr::Systematics	syst = semilCorr::ALL
	)

Definition at line 95 of file semilCorr.cxx.

{
  //vector<int> indices = getHistoIndices(up,syst);
  vector<int> indices1 = getHistoIndices(syst);
  float systr = 0;
  for(unsigned int i = 0; i<indices1.size(); i++){
    systr += pow(getSemilCorrToIncl(jet,mu,m_histos[indices1[i]]),2);
  }
  return sqrt(systr);
}

Member Data Documentation

m_Debug

bool semilCorr::m_Debug

private

Definition at line 35 of file semilCorr.h.

m_etas

std::vector<float> semilCorr::m_etas

private

Definition at line 31 of file semilCorr.h.

m_f

TFile* semilCorr::m_f

private

Definition at line 33 of file semilCorr.h.

m_histos

std::vector<std::vector<TH1F*> > semilCorr::m_histos

private

Definition at line 30 of file semilCorr.h.

---

The documentation for this class was generated from the following files:

• semilCorr.h
• semilCorr.cxx