Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Changing linkage of the "timeshift_ns_hbheho" function to external #14957

Merged
merged 4 commits into from Jun 29, 2016
Merged

Changing linkage of the "timeshift_ns_hbheho" function to external #14957

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
491caf9
Changing linkage of timeshift_ns_hbheho function to external
Jun 23, 2016
0f09376
turns out, we need to change the linkage of one more function
Jun 23, 2016
a11b65c
renamig a file
Jun 26, 2016
a3397a0
updating the guard name
Jun 26, 2016
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
14 changes: 14 additions & 0 deletions RecoLocalCalo/HcalRecAlgos/interface/HcalCorrectionFunctions.h
View file
Open in desktop
@@ -0,0 +1,14 @@
#ifndef RecoLocalCalo_HcalRecAlgos_HcalCorrectionFunctions_h_
#define RecoLocalCalo_HcalRecAlgos_HcalCorrectionFunctions_h_

///Timeshift correction for HPDs based on the position of the peak ADC measurement.
/// Allows for an accurate determination of the relative phase of the pulse shape from
/// the HPD. Calculated based on a weighted sum of the -1,0,+1 samples relative to the peak
/// as follows: wpksamp = (0*sample[0] + 1*sample[1] + 2*sample[2]) / (sample[0] + sample[1] + sample[2])
/// where sample[1] is the maximum ADC sample value.
float timeshift_ns_hbheho(float wpksamp);

/// Special energy correction for some HB- cells
float hbminus_special_ecorr(int ieta, int iphi, double energy, int runnum);

#endif // RecoLocalCalo_HcalRecAlgos_HcalCorrectionFunctions_h_
27 changes: 9 additions & 18 deletions RecoLocalCalo/HcalRecAlgos/src/HcalSimpleRecAlgo.cc
View file
Open in desktop
Expand Up @@ -3,6 +3,7 @@
#include "CalibCalorimetry/HcalAlgos/interface/HcalTimeSlew.h"
#include "RecoLocalCalo/HcalRecAlgos/src/HcalTDCReco.h"
#include "RecoLocalCalo/HcalRecAlgos/interface/rawEnergy.h"
#include "RecoLocalCalo/HcalRecAlgos/interface/HcalCorrectionFunctions.h"
#include "DataFormats/METReco/interface/HcalCaloFlagLabels.h"

#include <algorithm>
Expand Down Expand Up @@ -107,19 +108,9 @@ void HcalSimpleRecAlgo::setBXInfo(const BunchXParameter* info,
lenBunchCrossingInfo_ = lenInfo;
}

///Timeshift correction for HPDs based on the position of the peak ADC measurement.
/// Allows for an accurate determination of the relative phase of the pulse shape from
/// the HPD. Calculated based on a weighted sum of the -1,0,+1 samples relative to the peak
/// as follows: wpksamp = (0*sample[0] + 1*sample[1] + 2*sample[2]) / (sample[0] + sample[1] + sample[2])
/// where sample[1] is the maximum ADC sample value.
static float timeshift_ns_hbheho(float wpksamp);

///Same as above, but for the HF PMTs.
///Timeshift correction for the HF PMTs.
static float timeshift_ns_hf(float wpksamp);

/// Ugly hack to apply energy corrections to some HB- cells
static float eCorr(int ieta, int iphi, double ampl, int runnum);

/// Leak correction
static float leakCorr(double energy);

Expand Down Expand Up @@ -396,9 +387,9 @@ namespace HcalSimpleRecAlgoImpl {
if (cell.subdet() == HcalBarrel) {
const int ieta = cell.ieta();
const int iphi = cell.iphi();
uncorr_ampl *= eCorr(ieta, iphi, uncorr_ampl, runnum);
ampl *= eCorr(ieta, iphi, ampl, runnum);
m3_ampl *= eCorr(ieta, iphi, m3_ampl, runnum);
uncorr_ampl *= hbminus_special_ecorr(ieta, iphi, uncorr_ampl, runnum);
ampl *= hbminus_special_ecorr(ieta, iphi, ampl, runnum);
m3_ampl *= hbminus_special_ecorr(ieta, iphi, m3_ampl, runnum);
}
}

Expand Down Expand Up @@ -504,9 +495,9 @@ namespace HcalSimpleRecAlgoImpl {
if (cell.subdet() == HcalBarrel) {
const int ieta = cell.ieta();
const int iphi = cell.iphi();
uncorr_ampl *= eCorr(ieta, iphi, uncorr_ampl, runnum);
ampl *= eCorr(ieta, iphi, ampl, runnum);
m3_ampl *= eCorr(ieta, iphi, m3_ampl, runnum);
uncorr_ampl *= hbminus_special_ecorr(ieta, iphi, uncorr_ampl, runnum);
ampl *= hbminus_special_ecorr(ieta, iphi, ampl, runnum);
m3_ampl *= hbminus_special_ecorr(ieta, iphi, m3_ampl, runnum);
}
}

Expand Down Expand Up @@ -665,7 +656,7 @@ HFRecHit HcalSimpleRecAlgo::reconstructHFUpgrade(const HcalUpgradeDataFrame& dig


/// Ugly hack to apply energy corrections to some HB- cells
float eCorr(int ieta, int iphi, double energy, int runnum) {
float hbminus_special_ecorr(int ieta, int iphi, double energy, int runnum) {
// return energy correction factor for HBM channels
// iphi=6 ieta=(-1,-15) and iphi=32 ieta=(-1,-7)
// I.Vodopianov 28 Feb. 2011
Expand Down