HcalPulseContainmentAlgo Class Reference

#include <HcalPulseContainmentAlgo.h>

Public Member Functions
std::pair< double, double >	calcpair (double)
	HcalPulseContainmentAlgo (const HcalPulseShape *shape, int num_samples, double fixedphase_ns, bool phaseAsInSim, const HcalTimeSlew *hcalTimeSlew_delay)
	HcalPulseContainmentAlgo (int num_samples, double fixedphase_ns, bool phaseAsInSim, const HcalTimeSlew *hcalTimeSlew_delay)

Private Member Functions
void	init (int num_samples)

Private Attributes
double	fixedphasens_
const HcalTimeSlew *	hcalTimeSlew_delay_
double	integrationwindowns_
HcalShapeIntegrator	integrator_
bool	phaseAsInSim_
double	time0shiftns_

Detailed Description

Definition at line 9 of file HcalPulseContainmentAlgo.h.

Constructor & Destructor Documentation

HcalPulseContainmentAlgo::HcalPulseContainmentAlgo	(	const HcalPulseShape *	shape,
		int	num_samples,
		double	fixedphase_ns,
		bool	phaseAsInSim,
		const HcalTimeSlew *	hcalTimeSlew_delay
	)

Definition at line 23 of file HcalPulseContainmentAlgo.cc.

References init().

    : fixedphasens_(fixedphase_ns),
      phaseAsInSim_(phaseAsInSim),
      integrator_(shape),
      hcalTimeSlew_delay_(hcalTimeSlew_delay) {
  init(num_samples);
}

HcalPulseContainmentAlgo::HcalPulseContainmentAlgo	(	int	num_samples,
		double	fixedphase_ns,
		bool	phaseAsInSim,
		const HcalTimeSlew *	hcalTimeSlew_delay
	)

Definition at line 12 of file HcalPulseContainmentAlgo.cc.

References init().

    : fixedphasens_(fixedphase_ns),
      phaseAsInSim_(phaseAsInSim),
      integrator_(&(HcalPulseShapes().hbShape())),
      hcalTimeSlew_delay_(hcalTimeSlew_delay) {
  init(num_samples);
}

Member Function Documentation

std::pair< double, double > HcalPulseContainmentAlgo::calcpair

(

double

truefc

)

Definition at line 70 of file HcalPulseContainmentAlgo.cc.

References HcalTimeSlew::delay(), fixedphasens_, hcalTimeSlew_delay_, funct::integral(), integrationwindowns_, integrator_, SiStripPI::max, HcalTimeSlew::Medium, phaseAsInSim_, alignCSCRings::s, time0shiftns_, and tmax.

                                                                        {
  double timeslew_ns = hcalTimeSlew_delay_->delay(std::max(0.0, (double)truefc), HcalTimeSlew::Medium);

  double tmin = 0;
  if (phaseAsInSim_) {  // timePhase as in hcalSimParameters, no time0shift
    tmin = fixedphasens_ - timeslew_ns;
  } else {  // Run 2: timePhase opposite to SIM, time0shift
    double shift_ns = fixedphasens_ - time0shiftns_ + timeslew_ns;
    //edm::LogPrint("HcalPulseContainmentAlgo") << "SHIFT " << fixedphasens_ << " " << time0shiftns_ << " " << timeslew_ns;
    tmin = -shift_ns;
  }
  double tmax = tmin + integrationwindowns_;

  //double integral  = shape_.integrate( tmin, tmax );
  double integral = integrator_(tmin, tmax);
  //edm::LogPrint("HcalPulseContainmentAlgo") << "INTEGRAL " << integral << " " << truefc << " " << tmin << " "  << tmax;
  double corfactor = 1.0 / integral;
  double recofc = (double)truefc * integral;

#if 0
  char s[80];
  sprintf (s, "%8.2f %8.4f %8.4f %8.5f %8.5f %8.5f ",
       truefc, tmin, tmax, integral, corfactor, recofc);
  edm::LogPrint("HcalPulseContainmentAlgo") << s;
#endif

  std::pair<double, double> thepair(recofc, corfactor);
  return thepair;
}

void HcalPulseContainmentAlgo::init ( int  num_samples )

private

Definition at line 35 of file HcalPulseContainmentAlgo.cc.

References integrationwindowns_, integrator_, alignCSCRings::s, and time0shiftns_.

Referenced by HcalPulseContainmentAlgo().

                                                   {
  const int binsize_ns = 25;

  // First set up controlling parameters for calculating the correction factor:
  // Integration window size...
  //
  integrationwindowns_ = (double)(binsize_ns * num_samples);

  // First find the point at which time bin "1" exceeds time bin "0",
  // and call that point "time 0".
  //
  for (int shift_ns = 0; shift_ns < binsize_ns; shift_ns++) {
    // Digitize by integrating to find all time sample
    // bin values for this shift.
    //
    double tmin = -(double)shift_ns;
    double bin0val = (double)integrator_(tmin, tmin + binsize_ns);
    double bin1val = (double)integrator_(tmin + binsize_ns, tmin + 2 * binsize_ns);

#if 0
    char s[80];
    sprintf (s, "%7.3f %8.5f %8.5f\n", tmin, bin0val, bin1val);
    edm::LogPrint("HcalPulseContainmentAlgo") << s;
#endif

    if (bin1val > bin0val) {
      time0shiftns_ = shift_ns;
      break;
    }
  }
#if 0
  edm::LogPrint("HcalPulseContainmentAlgo") << "time0shiftns_ = " << time0shiftns_;
#endif
}

Member Data Documentation

double HcalPulseContainmentAlgo::fixedphasens_

private

Definition at line 24 of file HcalPulseContainmentAlgo.h.

Referenced by calcpair().

const HcalTimeSlew* HcalPulseContainmentAlgo::hcalTimeSlew_delay_

private

Definition at line 29 of file HcalPulseContainmentAlgo.h.

Referenced by calcpair().

double HcalPulseContainmentAlgo::integrationwindowns_

private

Definition at line 25 of file HcalPulseContainmentAlgo.h.

Referenced by calcpair(), and init().

HcalShapeIntegrator HcalPulseContainmentAlgo::integrator_

private

Definition at line 28 of file HcalPulseContainmentAlgo.h.

Referenced by calcpair(), and init().

bool HcalPulseContainmentAlgo::phaseAsInSim_

private

Definition at line 27 of file HcalPulseContainmentAlgo.h.

Referenced by calcpair().

double HcalPulseContainmentAlgo::time0shiftns_

private

Definition at line 26 of file HcalPulseContainmentAlgo.h.

Referenced by calcpair(), and init().