calib::AttenCurve Class Reference

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-16/Calibrator/func/AttenCurve.h"

Classes
struct	Pt

Public Member Functions
float	MeanPEPerCmAt (double w) const
	Mean response of this channel at this distance from detector centre. More...
float	MeanPEPerCmAt (double w, std::string version) const
bool	IsCalibrated () const
void	AddInterpPoint (float w, float factor)
void	WriteToCSVs (FILE *fConsts, FILE *fPoints, bool mc) const

Static Public Member Functions
static AttenCurve *	Uninitialized (int det, geo::OfflineChan chan)
	Return a new AttenCurve objects with fields uninitialized. More...

Public Attributes
std::vector< Pt >	interp_pts
float	atten_length
float	coeff_exp
float	background
float	edge_low
float	edge_high
float	coeff_low
float	coeff_high
float	chisq
float	cell_length
float	center_offset
	Nonzero in short muon-catcher cells. Positive is closer to readout. More...
int	det
geo::OfflineChan	chan
bool	initialized

Protected Member Functions
	AttenCurve ()

Detailed Description

Definition at line 18 of file AttenCurve.h.

Constructor & Destructor Documentation

calib::AttenCurve::AttenCurve ( )

inlineprotected

Definition at line 57 of file AttenCurve.h.

References calib::operator<<().

Referenced by Uninitialized().

{}

Member Function Documentation

void calib::AttenCurve::AddInterpPoint	(	float	w,
		float	factor
	)

Definition at line 122 of file AttenCurve.cxx.

References interp_pts, gen_hdf5record::pt, and w.

Referenced by calib::AttenFit::lowessFit(), and calib::AttenuationFit::lowessFit().

  {
    Pt pt = {w, factor};
    interp_pts.push_back(pt);
  }

bool calib::AttenCurve::IsCalibrated

(

)

const

Definition at line 97 of file AttenCurve.cxx.

References chisq.

Referenced by calib::Calibrator::GetAttenCurve().

  {
    // This isn't actually a chisq, but more an average fractional deviation of
    // the fit from the data. 0.2 is the eyeballed position of the tails of the
    // distribution in FD data.
    return chisq >= 0 && chisq < 0.2;

    // Checking the LOWESS points is another way to achieve something
    // similar. But this implementation tends to give to much weight to
    // problems that only occur in single locations.

    // for(unsigned int n = 0; n < interp_pts.size(); ++n){
    //   // Don't let the interpolated corrections pull more than 15% off the
    //   // original fit. 15% is arbitrary, based on a desire to calibrate a
    //   // reasonable number of channels.
    //   if(fabs(interp_pts[n].factor-1) > .15){
    //     // Don't count interpolations where just the final point is off as bad
    //     if(n > 0 && n < interp_pts.size()-1)
    //       return false;
    //   }
    // }
    // return true;
  }

float calib::AttenCurve::MeanPEPerCmAt

(

double

w

)

const

Mean response of this channel at this distance from detector centre.

Definition at line 37 of file AttenCurve.cxx.

Referenced by calib::FindOutliers::endRun(), calib::AttenFit::fit_channel_prof(), calib::AttenuationFit::fit_channel_prof(), and calib::Calibrator::GetAttenScale().

  {
    return MeanPEPerCmAt(w,"");
  }

float calib::AttenCurve::MeanPEPerCmAt	(	double	w,
		std::string	version
	)		const

Definition at line 43 of file AttenCurve.cxx.

References ana::assert(), atten_length, background, cell_length, coeff_exp, coeff_high, coeff_low, edge_high, edge_low, stan::math::exp(), calib::FOURTH(), initialized, interp_pts, w, and y1.

  {
    assert(initialized);

    //Now check just once in Calibrator instead
    //if(!IsCalibrated()) return -1;

    // TODO: take account of center_offset here for muon catcher planes. Need
    // to fix the todo below at the same time.

    //The attenuation formulae has a minor correction relative to tag v6 and earlier. 
    //Backwards compatibility is ensured in calibrator by passing the tag version to 
    //this function. Users can also call this function directly without specifying a
    //tag in which case they get the latest, corrected, version of the attenuation
    //formulae.

    // Distance to travel the long way round the fibre
    double wlong;
    if(version=="v1"||version=="v2"||version=="v3"||version=="v3.1"||version=="v4"||version=="v5"||version=="v5.1"||version=="v6"){
      wlong = 1.5*cell_length+w;
    }
    else{
      wlong = cell_length+w;
    }
    
    const double expPart = coeff_exp*(exp(w/atten_length)+
                                      exp(-wlong/atten_length))+background;

    double edgePart = 1;
    if(w < -edge_low) edgePart -= coeff_low *FOURTH(w+edge_low);
    if(w > edge_high) edgePart -= coeff_high*FOURTH(w-edge_high);

    double interpPart = 1;
    if(interp_pts.size() >= 2){
      typedef std::vector<Pt>::const_iterator it_t;

      // Find the point after w
      it_t itNext = std::lower_bound(interp_pts.begin(), interp_pts.end(), w);
      if(itNext == interp_pts.end()) --itNext;
      if(itNext == interp_pts.begin()) ++itNext;
      // And the point before
      it_t itPrev = itNext; --itPrev;

      // Interpolate
      const double w0 = itPrev->w; const double y0 = itPrev->factor;
      const double w1 = itNext->w; const double y1 = itNext->factor;

      interpPart = ((w1-w)*y0+(w-w0)*y1)/(w1-w0);
    }

    return expPart*edgePart*interpPart;
  }

AttenCurve * calib::AttenCurve::Uninitialized ( int  det, geo::OfflineChan  chan  )

static

Return a new AttenCurve objects with fields uninitialized.

Definition at line 25 of file AttenCurve.cxx.

References AttenCurve(), center_offset, chan, det, and initialized.

Referenced by calib::AttenFit::endJob(), calib::AttenCache::LoadFromDatabase(), and calib::AttenuationFit::writeResults().

  {
    AttenCurve* res = new AttenCurve;
    res->det = det;
    res->chan = chan;
    res->initialized = false;
    res->center_offset = 0;
    return res;
  }

void calib::AttenCurve::WriteToCSVs	(	FILE *	fConsts,
		FILE *	fPoints,
		bool	mc
	)		const

Definition at line 129 of file AttenCurve.cxx.

References ana::assert(), atten_length, background, chan, chisq, coeff_exp, coeff_high, coeff_low, edge_high, edge_low, initialized, interp_pts, getGoodRuns4SAM::n, geo::OfflineChan::ToDBValidityChan(), and w.

Referenced by calib::AttenFit::endJob(), and calib::AttenuationFit::writeResults().

  {
    const int vldTime = 1; // 1970


    assert(initialized);
    fprintf(fConsts, "%g,%g,%g,%g,%g,%g,%g,%g,%d,%d\n",
            coeff_exp, atten_length, background,
            edge_low, edge_high, coeff_low, coeff_high, chisq,
            chan.ToDBValidityChan(), vldTime);

    const int N = interp_pts.size();
    for(int n = 0; n < N; ++n){
      fprintf(fPoints, "%g,%g,%d,%d\n",
              interp_pts[n].w, interp_pts[n].factor,
              1000*chan.ToDBValidityChan()+n, vldTime);
    }
  }

Member Data Documentation

float calib::AttenCurve::atten_length

Definition at line 45 of file AttenCurve.h.

Referenced by calib::AttenFit::expFit(), calib::AttenuationFit::expFit(), calib::AttenFit::fit_channel_prof(), calib::AttenuationFit::fit_channel_prof(), calib::AttenCache::LoadFromDatabase(), MeanPEPerCmAt(), calib::operator<<(), and WriteToCSVs().

float calib::AttenCurve::background

Definition at line 45 of file AttenCurve.h.

Referenced by calib::AttenFit::expFit(), calib::AttenuationFit::expFit(), calib::AttenFit::fit_channel_prof(), calib::AttenuationFit::fit_channel_prof(), calib::AttenCache::LoadFromDatabase(), MeanPEPerCmAt(), calib::operator<<(), and WriteToCSVs().

float calib::AttenCurve::cell_length

Definition at line 49 of file AttenCurve.h.

Referenced by calib::AttenFit::expFit(), calib::AttenuationFit::expFit(), calib::AttenFit::fit_channel_prof(), calib::AttenuationFit::fit_channel_prof(), calib::AttenCache::LoadFromDatabase(), and MeanPEPerCmAt().

float calib::AttenCurve::center_offset

Nonzero in short muon-catcher cells. Positive is closer to readout.

Definition at line 51 of file AttenCurve.h.

Referenced by calib::AttenCache::LoadFromDatabase(), and Uninitialized().

geo::OfflineChan calib::AttenCurve::chan

Definition at line 53 of file AttenCurve.h.

Referenced by calib::AttenCache::Add(), calib::AttenFit::endJob(), calib::Calibrator::GetAttenScale(), calib::AttenCache::LoadFromDatabase(), calib::operator<<(), Uninitialized(), calib::AttenuationFit::writeResults(), and WriteToCSVs().

float calib::AttenCurve::chisq

Definition at line 48 of file AttenCurve.h.

Referenced by calib::FindOutliers::endRun(), calib::AttenFit::expFit(), calib::AttenuationFit::expFit(), IsCalibrated(), calib::AttenCache::LoadFromDatabase(), calib::operator<<(), calib::AttenFit::rolloffFit(), calib::AttenuationFit::rolloffFit(), calib::AttenuationFit::writeResults(), and WriteToCSVs().

float calib::AttenCurve::coeff_exp

Definition at line 45 of file AttenCurve.h.

Referenced by calib::AttenFit::expFit(), calib::AttenuationFit::expFit(), calib::AttenFit::fit_channel_prof(), calib::AttenuationFit::fit_channel_prof(), calib::AttenCache::LoadFromDatabase(), MeanPEPerCmAt(), calib::operator<<(), and WriteToCSVs().

float calib::AttenCurve::coeff_high

Definition at line 47 of file AttenCurve.h.

Referenced by calib::AttenFit::expFit(), calib::AttenuationFit::expFit(), calib::AttenCache::LoadFromDatabase(), MeanPEPerCmAt(), calib::operator<<(), calib::AttenFit::rolloffFit(), calib::AttenuationFit::rolloffFit(), and WriteToCSVs().

float calib::AttenCurve::coeff_low

Definition at line 47 of file AttenCurve.h.

Referenced by calib::AttenFit::expFit(), calib::AttenuationFit::expFit(), calib::AttenCache::LoadFromDatabase(), MeanPEPerCmAt(), calib::operator<<(), calib::AttenFit::rolloffFit(), calib::AttenuationFit::rolloffFit(), and WriteToCSVs().

int calib::AttenCurve::det

Definition at line 52 of file AttenCurve.h.

Referenced by calib::AttenCache::Add(), calib::operator<<(), and Uninitialized().

float calib::AttenCurve::edge_high

Definition at line 46 of file AttenCurve.h.

Referenced by calib::AttenFit::expFit(), calib::AttenuationFit::expFit(), calib::AttenFit::fit_channel_prof(), calib::AttenuationFit::fit_channel_prof(), calib::AttenCache::LoadFromDatabase(), MeanPEPerCmAt(), calib::operator<<(), calib::AttenFit::rolloffFit(), calib::AttenuationFit::rolloffFit(), and WriteToCSVs().

float calib::AttenCurve::edge_low

Definition at line 46 of file AttenCurve.h.

Referenced by calib::AttenFit::expFit(), calib::AttenuationFit::expFit(), calib::AttenFit::fit_channel_prof(), calib::AttenuationFit::fit_channel_prof(), calib::AttenCache::LoadFromDatabase(), MeanPEPerCmAt(), calib::operator<<(), calib::AttenFit::rolloffFit(), calib::AttenuationFit::rolloffFit(), and WriteToCSVs().

bool calib::AttenCurve::initialized

Definition at line 54 of file AttenCurve.h.

Referenced by calib::AttenFit::fit_channel_prof(), calib::AttenuationFit::fit_channel_prof(), calib::AttenCache::LoadFromDatabase(), MeanPEPerCmAt(), calib::operator<<(), Uninitialized(), and WriteToCSVs().

std::vector<Pt> calib::AttenCurve::interp_pts

Definition at line 43 of file AttenCurve.h.

Referenced by AddInterpPoint(), calib::AttenCache::LoadFromDatabase(), MeanPEPerCmAt(), calib::operator<<(), and WriteToCSVs().

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-16/Calibrator/func/AttenCurve.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-16/Calibrator/func/AttenCurve.cxx