KKCorrectionFactors Class Reference

#include <KKCorrectionFactors.h>

Public Member Functions
float	getScale (float genEnergy, float genEta, float simEnergy) const
	KKCorrectionFactors (const edm::ParameterSet &pset)
	~KKCorrectionFactors ()

Private Attributes
TH3F *	h3_
bool	interpolate3D_

Detailed Description

Definition at line 27 of file KKCorrectionFactors.h.

Constructor & Destructor Documentation

KKCorrectionFactors::KKCorrectionFactors

(

const edm::ParameterSet &

pset

)

Definition at line 3 of file KKCorrectionFactors.cc.

References Exception, edm::ParameterSet::exists(), MillePedeFileConverter_cfg::fileName, edm::ParameterSet::getUntrackedParameter(), h3_, interpolate3D_, NewMuonAssociatorByHits_cfi::obj, and AlCaHLTBitMon_QueryRunRegistry::string.

{

  interpolate3D_ = pset.exists("interpolate3D") && pset.getUntrackedParameter<bool>("interpolate3D");

  // Get the filename and histogram name
  std::string fileName = pset.getUntrackedParameter<std::string>("fileName");
  std::string histogramName = pset.getUntrackedParameter<std::string>("histogramName");

  // Read histo
  edm::FileInPath myDataFile( fileName );
  TFile * myFile = TFile::Open( myDataFile.fullPath().c_str() );

  gROOT->cd(); // create histogram in memory
  auto obj = myFile->Get( histogramName.c_str() );

  // Check if histogram exists in file
  if(!obj) {
    throw cms::Exception( "FileReadError" )
      << "Histogram \"" << histogramName
      << "\" not found in file \"" << fileName
      << "\".\n";
  }
  h3_ = new TH3F( *((TH3F*)obj) );

  delete myFile;

}

KKCorrectionFactors::~KKCorrectionFactors ( )

inline

Definition at line 34 of file KKCorrectionFactors.h.

References getScale(), and h3_.

{ delete h3_; };

Member Function Documentation

float KKCorrectionFactors::getScale	(	float	genEnergy,
		float	genEta,
		float	simEnergy
	)		const

Definition at line 33 of file KKCorrectionFactors.cc.

References h3_, interpolate3D_, alignCSCRings::r, and Scenarios_cff::scale.

Referenced by ~KKCorrectionFactors().

{

  float r = simE / genE;
  float scale = 1.;

  if( interpolate3D_
      // TH3::Interpolate can only interpolate inside the bondaries of the histogram
      && genE > h3_->GetXaxis()->GetXmin()
      && genE < h3_->GetXaxis()->GetXmax()
      && genEta > h3_->GetYaxis()->GetXmin()
      && genEta < h3_->GetYaxis()->GetXmax()
      && r < h3_->GetZaxis()->GetXmax()
      && r > h3_->GetZaxis()->GetXmax() ) {

    scale = h3_->Interpolate( genE, genEta, r );

  } else { // interpolation in r is mandatory

    int binE   = h3_->GetXaxis()->FindFixBin( genE );
    int binEta = h3_->GetYaxis()->FindFixBin( genEta );

    // find the two bins which are closest to the actual value
    auto binWidthR = h3_->GetZaxis()->GetBinWidth(0);
    int binRup = h3_->GetZaxis()->FindFixBin( r + binWidthR/2. );
    int binRdn = h3_->GetZaxis()->FindFixBin( r - binWidthR/2. );

    auto scaleUp = h3_->GetBinContent( binE, binEta, binRup );
    auto scaleDn = h3_->GetBinContent( binE, binEta, binRdn );

    // make a linear extrapolation between neighbour bins if they are not zero
    auto Rdn = h3_->GetZaxis()->GetBinCenter( binRdn );
    scale = scaleDn + (scaleUp-scaleDn) * ( r - Rdn ) / binWidthR;

  }
  return scale;

}

Member Data Documentation

TH3F* KKCorrectionFactors::h3_

private

Definition at line 40 of file KKCorrectionFactors.h.

Referenced by getScale(), KKCorrectionFactors(), and ~KKCorrectionFactors().

bool KKCorrectionFactors::interpolate3D_

private

Definition at line 41 of file KKCorrectionFactors.h.

Referenced by getScale(), and KKCorrectionFactors().