icarus::calo::NormalizeTPC Class Reference

Inheritance diagram for icarus::calo::NormalizeTPC:

Classes
class	ScaleInfo

Public Member Functions
	NormalizeTPC (fhicl::ParameterSet const &pset)
void	configure (const fhicl::ParameterSet &pset) override
double	Normalize (double dQdx, const art::Event &e, const recob::Hit &h, const geo::Point_t &location, const geo::Vector_t &direction, double t0) override
Public Member Functions inherited from INormalizeCharge
virtual	~INormalizeCharge () noexcept=default
	Virtual Destructor. More...

Private Member Functions
ScaleInfo	GetScaleInfo (uint64_t timestamp)
std::string	URL (uint64_t timestamp)

Private Attributes
int	fTimeout
std::string	fURL
bool	fVerbose
std::map< uint64_t, ScaleInfo >	fScaleInfos

Detailed Description

Definition at line 27 of file NormalizeTPC_tool.cc.

Constructor & Destructor Documentation

icarus::calo::NormalizeTPC::NormalizeTPC

(

fhicl::ParameterSet const &

pset

)

Definition at line 61 of file NormalizeTPC_tool.cc.

                                                                  {
  this->configure(pset);
}

Member Function Documentation

void icarus::calo::NormalizeTPC::configure ( const fhicl::ParameterSet &  pset )

overridevirtual

Implements INormalizeCharge.

Definition at line 65 of file NormalizeTPC_tool.cc.

                                                                    {
  fURL = pset.get<std::string>("URL");
  fTimeout = pset.get<unsigned>("Timeout");
  fVerbose = pset.get<bool>("Verbose", false);
}

icarus::calo::NormalizeTPC::ScaleInfo icarus::calo::NormalizeTPC::GetScaleInfo ( uint64_t  timestamp )

private

Definition at line 75 of file NormalizeTPC_tool.cc.

                                                                                       {
  // check the cache
  if (fScaleInfos.count(timestamp)) {
    return fScaleInfos.at(timestamp);
  }

  // Otherwise, look it up
  int error = 0;
  std::string url = URL(timestamp);

  if (fVerbose) std::cout << "NormalizeTPC Tool -- New Scale info, requesting data from url:\n" << url << std::endl;

  Dataset d = getDataWithTimeout(url.c_str(), "", fTimeout, &error);
  if (error) {
    throw cet::exception("NormalizeTPC") << "Calibration Database access failed. URL: (" << url << ") Error Code: " << error;
  }

  if (fVerbose) std::cout << "NormalizeTPC Tool -- Received HTTP response:\n" << getHTTPmessage(d) << std::endl;

  if (getHTTPstatus(d) != 200) {
    throw cet::exception("NormalizeTPC") 
      << "Calibration Database access failed. URL: (" << url
      << "). HTTP error status: " << getHTTPstatus(d) << ". HTTP error message: " << getHTTPmessage(d);
  }

  // Collect the timestamp info
  ScaleInfo thisscale;

  // Number of rows
  int n_tuple = getNtuples(d);
  if (n_tuple < 0) {
    throw cet::exception("NormalizeTPC") << "Calibration Database access failed. URL: (" << url << ") Bad Tuple Number: " << n_tuple;
  }

  // Iterate over the rows
  // The first 4 are metadata
  for (unsigned row = 4; row < (unsigned)n_tuple; row++) {
    Tuple tup = getTuple(d, row);

    int err = 0;
    // Get the itpc number
    int ch = getLongValue(tup, 0, &err);
    if (error) {
      throw cet::exception("NormalizeTPC") << "Calibration Database access failed. URL: (" << url << ") Failed on tuple access, row: " << row << ", col 0. Error Code: " << error;
    }

    // and the scale
    double scale = getDoubleValue(tup, 2, &err);
    if (error) {
      throw cet::exception("NormalizeTPC") << "Calibration Database access failed. URL: (" << url << ") Failed on tuple access, row: " << row << ", col 1. Error Code: " << error;
    }

    thisscale.scale[ch] = scale;
  }

  // Set the cache
  fScaleInfos[timestamp] = thisscale;

  return thisscale;
}

double icarus::calo::NormalizeTPC::Normalize ( double  dQdx, const art::Event &  e, const recob::Hit &  h, const geo::Point_t &  location, const geo::Vector_t &  direction, double  t0  )

overridevirtual

Implements INormalizeCharge.

Definition at line 136 of file NormalizeTPC_tool.cc.

                                                                                              {
  // Get the info
  ScaleInfo i = GetScaleInfo(e.time().timeHigh());

  // Lookup the TPC, cryo
  unsigned tpc = hit.WireID().TPC;
  unsigned cryo = hit.WireID().Cryostat;
  // Get the TPC index
  unsigned itpc = 2*cryo + tpc/2;

  double scale = 1;

  // TODO: what to do if no scale is found? throw an exception??
  if (i.scale.count(itpc)) scale = i.scale.at(itpc);

  if (fVerbose) std::cout << "NormalizeTPC Tool -- Data at itpc: " << itpc << " scale: " << scale << std::endl;

  return dQdx * scale;
}

std::string icarus::calo::NormalizeTPC::URL ( uint64_t  timestamp )

private

Definition at line 71 of file NormalizeTPC_tool.cc.

                                                        {
  return fURL + std::to_string(timestamp);
}

Member Data Documentation

std::map<uint64_t, ScaleInfo> icarus::calo::NormalizeTPC::fScaleInfos

private

Definition at line 52 of file NormalizeTPC_tool.cc.

int icarus::calo::NormalizeTPC::fTimeout

private

Definition at line 37 of file NormalizeTPC_tool.cc.

std::string icarus::calo::NormalizeTPC::fURL

private

Definition at line 38 of file NormalizeTPC_tool.cc.

bool icarus::calo::NormalizeTPC::fVerbose

private

Definition at line 39 of file NormalizeTPC_tool.cc.

---

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

• NormalizeTPC_tool.cc