HcalTimeSlewSim Class Reference

#include <HcalTimeSlewSim.h>

Public Member Functions
void	delay (CaloSamples &samples, CLHEP::HepRandomEngine *, const HcalTimeSlew *hcalTimeSlew_delay) const
	HcalTimeSlewSim (const CaloVSimParameterMap *parameterMap, double minFCToDelay)

Private Member Functions
double	charge (const CaloSamples &samples) const

Private Attributes
double	minFCToDelay_
const CaloVSimParameterMap *	theParameterMap

Detailed Description

Definition at line 19 of file HcalTimeSlewSim.h.

Constructor & Destructor Documentation

HcalTimeSlewSim()

HcalTimeSlewSim::HcalTimeSlewSim	(	const CaloVSimParameterMap *	parameterMap,
		double	minFCToDelay
	)

Definition at line 13 of file HcalTimeSlewSim.cc.

    : theParameterMap(parameterMap), minFCToDelay_(minFCToDelay) {}

Member Function Documentation

charge()

double HcalTimeSlewSim::charge ( const CaloSamples &  samples ) const

private

Definition at line 17 of file HcalTimeSlewSim.cc.

                                                               {
  double totalCharge = 0.;
  for (int i = 0; i < 4; ++i) {
    int bin = i + samples.presamples();
    if (bin < samples.size()) {
      totalCharge += samples[bin];
    }
  }
  return totalCharge;
}

References newFWLiteAna::bin, mps_fire::i, and EgammaValidation_cff::samples.

delay()

void HcalTimeSlewSim::delay	(	CaloSamples &	samples,
		CLHEP::HepRandomEngine *	engine,
		const HcalTimeSlew *	hcalTimeSlew_delay
	)		const

Definition at line 28 of file HcalTimeSlewSim.cc.

                                                                          {
  // HO goes slow, HF shouldn't be used at all
  //ZDC not used for the moment
 
  DetId detId(cs.id());
  if (detId.det() == DetId::Calo && detId.subdetId() == HcalZDCDetId::SubdetectorId)
    return;
  HcalDetId hcalDetId(detId);
 
  if (hcalDetId.subdet() == HcalBarrel || hcalDetId.subdet() == HcalEndcap || hcalDetId.subdet() == HcalOuter) {
    HcalTimeSlew::BiasSetting biasSetting =
        (hcalDetId.subdet() == HcalOuter) ? HcalTimeSlew::Slow : HcalTimeSlew::Medium;
 
    // double totalCharge = charge(cs); // old TS...
 
    int maxbin = cs.size();
    CaloSamples data(detId, maxbin);  // for a temporary copy
    data = cs;
 
    // smearing
    int soi = cs.presamples();
    double eps = 1.e-6;
    double scale_factor = 0.6;
    double scale = cs[soi] / scale_factor;
    double smearns = 0.;
    double cut = minFCToDelay_;  //5. fC (above mean) for signal to be delayed
 
    const HcalSimParameters& params = static_cast<const HcalSimParameters&>(theParameterMap->simParameters(detId));
    if (params.doTimeSmear()) {
      double rms = params.timeSmearRMS(scale);
      smearns = CLHEP::RandGaussQ::shoot(engine) * rms;
      LogDebug("HcalTimeSlewSim") << "TimeSmear charge " << scale << " rms " << rms << " smearns " << smearns;
    }
 
    // cycle over TS',  it - current TS index
    for (int it = 0; it < maxbin - 1;) {
      double datai = cs[it];
      int nts = 0;
      double tshift = smearns;
      double totalCharge = datai / scale_factor;
 
      // until we get more precise/reliable QIE8 simulation...
      if (totalCharge <= 0.)
        totalCharge = eps;  // protecion against negaive v.
      tshift += hcalTimeSlew_delay->delay(totalCharge, biasSetting);
      if (tshift <= 0.)
        tshift = eps;
 
      if (cut > -999.) {  //preserve compatibility
        if ((datai > cut) && (it == 0 || (datai > cs[it - 1]))) {
          // number of TS affected by current move depends on the signal shape:
          // rising or peaking
          nts = 2;  // rising
          if (datai > cs[it + 1])
            nts = 3;  // peaking
 
          // 1 or 2 TS to move from here,
          // 2d or 3d TS gets leftovers to preserve the sum
          for (int j = it; j < it + nts && j < maxbin; ++j) {
            // snippet borrowed from  CaloSamples::offsetTime(offset)
            // CalibFormats/CaloObjects/src/CaloSamples.cc
            double t = j * 25. - tshift;
            int firstbin = floor(t / 25.);
            double f = t / 25. - firstbin;
            int nextbin = firstbin + 1;
            double v1 = (firstbin < 0 || firstbin >= maxbin) ? 0. : cs[firstbin];
            double v2 = (nextbin < 0 || nextbin >= maxbin) ? 0. : cs[nextbin];
 
            // Keeping the overal sum intact
            if (nts == 2) {
              if (j == it)
                data[j] = v2 * f;
              else
                data[j] = v1 * (1. - f) + v2;
            } else {  // nts = 3
              if (j == it)
                data[j] = v2 * f;
              if (j == it + 1)
                data[j] = v1 * (1. - f) + v2 * f;
              if (j == it + nts - 1)
                data[j] = v1 * (1. - f) + v2;
            }
 
          }  // end of local move of TS', now update...
          cs = data;
 
        }  // end of rising edge or peak finding
      } else {
        double t = it * 25. - tshift;
        int firstbin = floor(t / 25.);
        double f = t / 25. - firstbin;
        int nextbin = firstbin + 1;
        double v2 = (nextbin < 0 || nextbin >= maxbin) ? 0. : data[nextbin];
        data[it] = v2 * f;
        data[it + 1] += (v2 - data[it]);
        cs = data;
      }
 
      if (nts < 3)
        it++;
      else
        it += 2;
    }
 
    // final update of the shifted TS array
    cs = data;
  }
}

References DetId::Calo, fwrapper::cs, TkAlMuonSelectors_cfi::cut, data, HcalTimeSlew::delay(), f, HcalBarrel, HcalEndcap, HcalOuter, dqmiolumiharvest::j, LogDebug, HcalTimeSlew::Medium, minFCToDelay_, CalibrationSummaryClient_cfi::params, SiStripPI::rms, Scenarios_cff::scale, CaloVSimParameterMap::simParameters(), HcalTimeSlew::Slow, HcalDetId::subdet(), HcalZDCDetId::SubdetectorId, OrderedSet::t, and theParameterMap.

Referenced by HcalAmplifier::amplify().

Member Data Documentation

minFCToDelay_

double HcalTimeSlewSim::minFCToDelay_

private

Definition at line 29 of file HcalTimeSlewSim.h.

Referenced by delay().

theParameterMap

const CaloVSimParameterMap* HcalTimeSlewSim::theParameterMap

private

Definition at line 28 of file HcalTimeSlewSim.h.

Referenced by delay().