calib::Calibrator Class Reference

#include <Calibrator.h>

List of all members.

Public Member Functions
	Calibrator ()
TTimeStamp	GetTimeStamp () const
rawdata::det_id_	GetDetector () const
void	MakeCellHit (const RawDigit *rawdigit, CellHit *cellhit)
float	GetPE (const RawDigit *rawdigit)
float	GetPECorrFactor (const RawDigit *rawdigit)
float	GetTNS (const RawDigit *rawdigit)
void	MakeRecoHit (const CellHit *cellhit, RecoHit *recohit, float w)
float	GetT (const CellHit *cellhit, float dist_to_readout)
float	GetPEAtt (const CellHit *cellhit, float dist_to_readout)
float	GetMIPScale (const CellHit *cellhit)
float	GetGeVScale (const CellHit *cellhit)
Static Public Member Functions
Calibrator &	Instance (TTimeStamp timestamp, rawdata::det_id_ detector, const geo::Geometry &geometry)
Public Attributes
float	fPEACperMIP
float	fPEACperMIPErr
std::vector< float >	fMIPtoGeV
std::vector< float >	fMIPtoGeVErr
std::map< int, float >	fADCtoPE
std::map< int, float >	fPEtoSigCor
std::map< int, attlen >	fSCtoPEAC
std::vector< float >	fFTDistNDX
std::vector< float >	fFTDistNDY
std::vector< float >	fFTDistFDX
std::vector< float >	fFTDistFDY
Private Member Functions
void	SetContext (TTimeStamp timestamp, rawdata::det_id_ detector, const geo::Geometry &geometry)
void	GetXYZD (const CellHit *cellhit, float w, double *xyzd)
bool	LoadADCtoPE ()
bool	LoadPEtoSigCor ()
bool	LoadSCtoPEAC ()
bool	LoadPEACtoMIP ()
bool	LoadMIPtoGeV ()
bool	LoadFineTimeDistros ()
Private Attributes
TTimeStamp	fCurrTimeStamp
rawdata::det_id_	fCurrDetector
const geo::Geometry *	fCurrGeo
bool	fLoadedFineTimeDistros
Static Private Attributes
Calibrator *	fInstance = 0

---

Constructor & Destructor Documentation

Calibrator::Calibrator (   )

Definition at line 26 of file Calibrator.cxx. References fADCtoPE, fFTDistFDX, fFTDistFDY, fFTDistNDX, fFTDistNDY, fMIPtoGeV, fMIPtoGeVErr, and rawdata::kNear. : fCurrTimeStamp(0), fCurrDetector(rawdata::kNear), fCurrGeo(0), fPEACperMIP(0), fPEACperMIPErr(0), fLoadedFineTimeDistros(0) { fADCtoPE.clear(); fMIPtoGeV.clear(); fMIPtoGeVErr.clear(); fFTDistNDX.clear(); fFTDistNDY.clear(); fFTDistFDX.clear(); fFTDistFDY.clear(); }

---

Member Function Documentation

rawdata::det_id_ calib::Calibrator::GetDetector (   )  const [inline]

Definition at line 50 of file Calibrator.h. { return fCurrDetector; }

float Calibrator::GetGeVScale (  const CellHit *  cellhit  )

Definition at line 526 of file Calibrator.cxx. Referenced by MakeRecoHit(). { // This needs to be calibrated and put in the DB, obviously! // Super-rough number based on 1 MIP = ~2 MeV/cm. return 0.013; }

float Calibrator::GetMIPScale (  const CellHit *  cellhit  )

Definition at line 520 of file Calibrator.cxx. Referenced by MakeRecoHit(). { // This needs to be calibrated and put in the DB, obviously! return 0.033; }

float Calibrator::GetPE (  const RawDigit *  rawdigit  )

Definition at line 355 of file Calibrator.cxx. References rawdata::RawDigit::ADC(). Referenced by MakeCellHit(). { // Look up gain from DB. For now, must hard code the // number. (This should not be put in any XML file. // That'll just confuse the issue.) float gain = 1.43; return rawdigit->ADC(0)/gain; }

float Calibrator::GetPEAtt (  const CellHit *  cellhit, float  dist_to_readout )

Definition at line 504 of file Calibrator.cxx. References recobase::CellHit::PECorr(). Referenced by MakeRecoHit(). { // These constants belong in the DB. Don't put them in an XML file // in the meantime, as this will just confuse the issue... Double_t F0 = 5.744; // P. Shanahan tuning from old CalHit.xml Double_t L0 = 315.3; // . Double_t F1 = 1.546; // . Double_t L1 = 3329.; // . Double_t attenuation = F0*TMath::Exp(-dist_to_readout/L0) + F1*TMath::Exp(-dist_to_readout/L1); return cellhit->PECorr() / attenuation; }

float Calibrator::GetPECorrFactor (  const RawDigit *  rawdigit  )

Definition at line 364 of file Calibrator.cxx. Referenced by MakeCellHit(). { // Look up channel-to-channel correction factors in DB. // For now, must hard code the factor. return 1.0; }

float Calibrator::GetT (  const CellHit *  cellhit, float  dist_to_readout )

Definition at line 494 of file Calibrator.cxx. References recobase::CellHit::TNS(). Referenced by MakeRecoHit(). { // This needs to be calibrated and put in a DB table Double_t speedOfFiberTransport = 15.3; // cm/ns, "first principles" calib. // Differs from c/n due to zigzag // paths down fiber. But, this // value may be way off (10%? 30%?). return (cellhit->TNS()-dist_to_readout/speedOfFiberTransport); }

TTimeStamp calib::Calibrator::GetTimeStamp (   )  const [inline]

Definition at line 49 of file Calibrator.h. { return fCurrTimeStamp; }

float Calibrator::GetTNS (  const RawDigit *  rawdigit  )

Definition at line 371 of file Calibrator.cxx. References rawdata::RawDigit::ADC(), fCurrDetector, rawdata::RawDigit::NADC(), and rawdata::RawDigit::TDC(). Referenced by MakeCellHit(). { // You'd think this would be a simple function... Not anymore! // The new readout scheme allows for super time resolution, and // this is the code that gets it. Double_t t = 0; // First, see if we're using the new-style readout. This // condition should be unnecessary in the future. if (rawdigit->NADC()<=1) { // Old style readout. Return old style answer: t = rawdigit->TDC(0)*62.5; } else { // New style readout. Use ADC info to determine fine timing. // See doc-???? for details. // ADC readings before and after peak, scaled to peak, // yield "x" and "y": Double_t peakadc = rawdigit->ADC(0); Double_t x = (rawdigit->ADC(1))/peakadc; Double_t y = (rawdigit->ADC(2))/peakadc; // These constants belong in the the DB. For now, hard code them. // Tuned up by RBP, 2009-10-07 Double_t Xa = 0.0; Double_t Ya = 1.0; Double_t Xb = 0.0; Double_t Yb = (fCurrDetector==rawdata::kNear?0.90:0.885); Double_t Xc = 1.0; Double_t Yc = 0.800; // These could be evaluated offline and stuck in the DB. // For now, I'll evaluate them explicitly here so the math // isn't lost: Double_t L1 = sqrt((Yb-Ya)*(Yb-Ya)+(Xb-Xa)*(Xb-Xa)); Double_t L2 = sqrt((Yc-Yb)*(Yc-Yb)+(Xc-Xb)*(Xc-Xb)); // distance to first line segment: Double_t U1 = ((x-Xa)*(Xb-Xa)+(y-Ya)*(Yb-Ya))/L1/L1; if (U1<0) U1=0; if (U1>0.99999) U1=0.99999; Double_t PX1 = Xa + U1*(Xb-Xa); Double_t PY1 = Ya + U1*(Yb-Ya); Double_t D1 = sqrt((PX1-x)*(PX1-x)+(PY1-y)*(PY1-y)); // distance to second line segment: Double_t U2 = ((x-Xb)*(Xc-Xb)+(y-Yb)*(Yc-Yb))/L2/L2; if (U2<0) U2=0; if (U2>0.99999) U2=0.99999; Double_t PX2 = Xb + U2*(Xc-Xb); Double_t PY2 = Yb + U2*(Yc-Yb); Double_t D2 = sqrt((PX2-x)*(PX2-x)+(PY2-y)*(PY2-y)); // get the better of the two Double_t U; if (D1<D2) U = U1*L1/(L1+L2); else U = (L1+U2*L2)/(L1+L2); // map best U value to a fine time: vector<float>& ftx = fFTDistNDX; vector<float>& fty = fFTDistNDY; if (fCurrDetector == rawdata::kFar) { ftx = fFTDistFDX; fty = fFTDistFDY; } unsigned int bin=1; if (U < ftx[0]) bin = 0; else for (; bin<ftx.size()-1; ++bin) { if (U >= ftx[bin-1] && U < ftx[bin]) break; } if (bin == ftx.size()) { // should never get in here std::cout << "In GetTNS(): bin=" << bin << " u=" << U << std::endl; bin=ftx.size()-1; } Double_t timefrac = 0.5*( fty[bin] + fty[bin-1]); // Gap between digitizations Double_t gap = (fCurrDetector==rawdata::kNear?125.0:500.0); t = rawdigit->TDC(0)*62.5 - timefrac*gap; } return t; }

void Calibrator::GetXYZD (  const CellHit *  cellhit, float  w, double *  xyzd )  [private]

Definition at line 533 of file Calibrator.cxx. References geo::PlaneGeo::Cell(), recobase::CellHit::Cell(), geo::CellGeo::DistToReadOut(), fCurrGeo, geo::CellGeo::GetCenter(), geo::PlaneGeo::Ncells(), geo::Geometry::NPlanes(), geo::Geometry::Plane(), recobase::CellHit::Plane(), recobase::CellHit::View(), and geo::CellGeo::WorldToLocal(). Referenced by MakeRecoHit(). { // Return (via 'xyzd' argument) position of deposition in world // coordinates and the distance to the readout, assuming deposition // has cell-local-z position = w. UInt_t plane = cellhit->Plane(); UInt_t cell = cellhit->Cell(); if (plane>=fCurrGeo->NPlanes()) { throw geo::Exception(__FILE__,__LINE__,geo::Exception::BAD_PLANE_NUMBER); } const geo::PlaneGeo& p = fCurrGeo->Plane(plane); if (cell>=p.Ncells()) { throw geo::Exception(__FILE__,__LINE__,geo::Exception::BAD_CELL_NUMBER); } const geo::CellGeo& c = p.Cell(cell); // Should be replaced with alignment tables or whatnot. Okay for now. // w is not localZ, which makes this a lot messier than just // c.GetCenter(xyzd,w) and c.DistToReadOut(w) c.GetCenter(xyzd); if (cellhit->View()==geo::kX) xyzd[1] = w; else xyzd[0] = w; Double_t localxyz[3]; c.WorldToLocal(xyzd,localxyz); xyzd[3] = c.DistToReadOut(localxyz[2]); }

Calibrator & Calibrator::Instance (  TTimeStamp  timestamp, rawdata::det_id_  detector, const geo::Geometry &  geometry )  [static]

Definition at line 63 of file Calibrator.cxx. References fInstance, and SetContext(). Referenced by util::EDMUtils::GetCalibrator(), and main(). { if (!fInstance) fInstance = new Calibrator; fInstance->SetContext(timestamp, detector, geometry); return *fInstance; }

bool Calibrator::LoadADCtoPE (   )  [private]

Definition at line 76 of file Calibrator.cxx. References db::Table::Connect(), db::Table::Disconnect(), fADCtoPE, fCurrTimeStamp, and db::Table::Init(). { fADCtoPE.clear(); std::string dbFile = "ADCtoPE.xml"; db::Table* table = new db::Table(dbFile); if (!table) return false; if (table->Connect()) { table->Init(fCurrTimeStamp.GetSec()); if (table->empty()) return false; table->Disconnect(); } // now fill STL map... std::string str1 = "Channel"; std::string str2 = "PEperADC"; int chan; float val; for (unsigned int i=0; i<table->size(); ++i) { (*table)[i].Get(str1,chan); (*table)[i].Get(str2,val); fADCtoPE[chan] = val; } delete table; return true; }

bool Calibrator::LoadFineTimeDistros (   )  [private]

Definition at line 270 of file Calibrator.cxx. References db::Table::Connect(), db::Table::Disconnect(), fCurrTimeStamp, fFTDistFDX, fFTDistFDY, fFTDistNDX, fFTDistNDY, fLoadedFineTimeDistros, and db::Table::Init(). Referenced by SetContext(). { fFTDistNDX.clear(); fFTDistNDY.clear(); fFTDistFDX.clear(); fFTDistFDY.clear(); // get values for ND first std::string dbFile = "FTDistroND.xml"; db::Table* table = new db::Table(dbFile); if (!table) return false; if (table->Connect()) { table->Init(fCurrTimeStamp.GetSec()); if (table->empty()) return false; table->Disconnect(); } float x,y; std::string strx = "x"; std::string stry = "y"; for (unsigned int i=0; i<table->size(); ++i) { (*table)[i].Get(strx,x); (*table)[i].Get(stry,y); fFTDistNDX.push_back(x); fFTDistNDY.push_back(y); } delete table; // now get values for FD dbFile = "FTDistroFD.xml"; table = new db::Table(dbFile); if (!table) return false; if (table->Connect()) { table->Init(fCurrTimeStamp.GetSec()); if (table->empty()) return false; table->Disconnect(); } for (unsigned int i=0; i<table->size(); ++i) { (*table)[i].Get(strx,x); (*table)[i].Get(stry,y); fFTDistFDX.push_back(x); fFTDistFDY.push_back(y); } delete table; fLoadedFineTimeDistros = 1; return true; }

bool Calibrator::LoadMIPtoGeV (   )  [private]

Definition at line 231 of file Calibrator.cxx. References db::Table::Connect(), db::Table::Disconnect(), fCurrTimeStamp, fMIPtoGeV, and db::Table::Init(). { fMIPtoGeV.clear(); std::string dbFile = "MIPtoGeV.xml"; db::Table* table = new db::Table(dbFile); if (!table) return false; if (table->Connect()) { table->Init(fCurrTimeStamp.GetSec()); if (table->empty()) return false; table->Disconnect(); } float val; char buf[20]; std::string str1; for (int i=0; i<5; ++i) { sprintf(buf,"par%d",i); str1 = buf; (*table)[0].Get(str1,val); fMIPtoGeV.push_back(val); } delete table; if (fMIPtoGeV.empty()) return false; return true; }

bool Calibrator::LoadPEACtoMIP (   )  [private]

Definition at line 201 of file Calibrator.cxx. References db::Table::Connect(), db::Table::Disconnect(), fCurrTimeStamp, fPEACperMIP, and db::Table::Init(). { fPEACperMIP = 1.; std::string dbFile = "PEACtoMIP.xml"; db::Table* table = new db::Table(dbFile); if (!table) return false; if (table->Connect()) { table->Init(fCurrTimeStamp.GetSec()); if (table->empty()) return false; table->Disconnect(); } std::string str1 = "PEACperMIP"; (*table)[0].Get(str1,fPEACperMIP); delete table; return true; }

bool Calibrator::LoadPEtoSigCor (   )  [private]

Definition at line 114 of file Calibrator.cxx. References db::Table::Connect(), db::Table::Disconnect(), fCurrTimeStamp, fPEtoSigCor, and db::Table::Init(). { fPEtoSigCor.clear(); std::string dbFile = "PEtoSigCor.xml"; db::Table* table = new db::Table(dbFile); if (!table) return false; if (table->Connect()) { table->Init(fCurrTimeStamp.GetSec()); if (table->empty()) return false; table->Disconnect(); } // now fill STL map... std::string str1 = "Channel"; std::string str2 = "SCperPE"; int chan; float val; for (unsigned int i=0; i<table->size(); ++i) { (*table)[i].Get(str1,chan); (*table)[i].Get(str2,val); fPEtoSigCor[chan] = val; } delete table; return true; }

bool Calibrator::LoadSCtoPEAC (   )  [private]

Definition at line 152 of file Calibrator.cxx. References calib::attlen::a1, calib::attlen::a1sig, calib::attlen::a2, calib::attlen::a2sig, db::Table::Connect(), db::Table::Disconnect(), fCurrTimeStamp, fSCtoPEAC, db::Table::Init(), calib::attlen::l1, calib::attlen::l1sig, calib::attlen::l2, and calib::attlen::l2sig. { fSCtoPEAC.clear(); std::string dbFile = "SigCortoPEAC.xml"; db::Table* table = new db::Table(dbFile); if (!table) return false; if (table->Connect()) { table->Init(fCurrTimeStamp.GetSec()); if (table->empty()) return false; table->Disconnect(); } // now fill STL map... std::string str1 = "Channel"; std::string str2 = "A1"; std::string str3 = "A2"; std::string str4 = "L1"; std::string str5 = "L2"; int chan; for (unsigned int i=0; i<table->size(); ++i) { attlen a; (*table)[i].Get(str1,chan); (*table)[i].Get(str2,a.a1); (*table)[i].Get(str2+"Err",a.a1sig); (*table)[i].Get(str3,a.a2); (*table)[i].Get(str3+"Err",a.a2sig); (*table)[i].Get(str4,a.l1); (*table)[i].Get(str4+"Err",a.l1sig); (*table)[i].Get(str5,a.l2); (*table)[i].Get(str5+"Err",a.l2sig); fSCtoPEAC[chan] = a; } delete table; return true; }

void Calibrator::MakeCellHit (  const RawDigit *  rawdigit, CellHit *  cellhit )

Definition at line 337 of file Calibrator.cxx. References fCurrDetector, fCurrGeo, fCurrTimeStamp, cmap::CMap::GetCell(), GetPE(), GetPECorrFactor(), cmap::CMap::GetPlane(), GetTNS(), cmap::CMap::Instance(), recobase::CellHit::PE(), recobase::CellHit::Plane(), geo::Geometry::Plane(), recobase::CellHit::SetCell(), recobase::CellHit::SetPE(), recobase::CellHit::SetPECorr(), recobase::CellHit::SetPlane(), recobase::CellHit::SetStatus(), recobase::CellHit::SetTNS(), and recobase::CellHit::SetView(). Referenced by calhit::CalHit::Reco(). { // copy inherited data members (*((RawDigit*)cellhit)) = (*rawdigit); cmap::CMap& cmap = cmap::CMap::Instance(fCurrTimeStamp,fCurrDetector); // new info cellhit->SetStatus(0); // for now, it always works! cellhit->SetCell (cmap.GetCell(rawdigit)); cellhit->SetPlane (cmap.GetPlane(rawdigit)); cellhit->SetView (fCurrGeo->Plane(cellhit->Plane()).View()); cellhit->SetTNS (GetTNS(rawdigit)); cellhit->SetPE (GetPE(rawdigit)); cellhit->SetPECorr(cellhit->PE()*GetPECorrFactor(rawdigit)); }

void Calibrator::MakeRecoHit (  const CellHit *  cellhit, RecoHit *  recohit, float  w )

Definition at line 470 of file Calibrator.cxx. References GetGeVScale(), GetMIPScale(), GetPEAtt(), GetT(), GetXYZD(), recobase::RecoHit::MIP(), recobase::RecoHit::PEAtt(), recobase::RecoHit::SetGeV(), recobase::RecoHit::SetMIP(), recobase::RecoHit::SetPEAtt(), recobase::RecoHit::SetT(), recobase::RecoHit::SetW(), recobase::RecoHit::SetX(), recobase::RecoHit::SetY(), and recobase::RecoHit::SetZ(). Referenced by calhit::CalHit::RoughReco(). { // copy inherited data members (*((CellHit*)recohit)) = (*cellhit); // new info: recohit->SetW(w); // Need position info to continue... Double_t xyzd[4]; GetXYZD(cellhit,w,xyzd); recohit->SetX (xyzd[0]); recohit->SetY (xyzd[1]); recohit->SetZ (xyzd[2]); recohit->SetT (GetT (cellhit,xyzd[3])); recohit->SetPEAtt(GetPEAtt(cellhit,xyzd[3])); recohit->SetMIP (recohit->PEAtt() * GetMIPScale(cellhit)); recohit->SetGeV (recohit->MIP() * GetGeVScale(cellhit)); }

void Calibrator::SetContext (  TTimeStamp  timestamp, rawdata::det_id_  detector, const geo::Geometry &  geometry )  [private]

Definition at line 45 of file Calibrator.cxx. References fCurrDetector, fCurrGeo, fCurrTimeStamp, and LoadFineTimeDistros(). Referenced by Instance(). { fCurrTimeStamp = timestamp; fCurrDetector = detector; fCurrGeo = &geometry; LoadFineTimeDistros(); /* LoadADCtoPE(); LoadPEtoSigCor(); LoadSCtoPEAC(); LoadPEACtoMIP(); LoadMIPtoGeV(); */ }

---

Member Data Documentation

std::map<int,float> calib::Calibrator::fADCtoPE

Definition at line 95 of file Calibrator.h. Referenced by Calibrator(), and LoadADCtoPE().

rawdata::det_id_ calib::Calibrator::fCurrDetector [private]

Definition at line 60 of file Calibrator.h. Referenced by GetTNS(), MakeCellHit(), and SetContext().

const geo::Geometry* calib::Calibrator::fCurrGeo [private]

Definition at line 61 of file Calibrator.h. Referenced by GetXYZD(), MakeCellHit(), and SetContext().

TTimeStamp calib::Calibrator::fCurrTimeStamp [private]

Definition at line 59 of file Calibrator.h. Referenced by LoadADCtoPE(), LoadFineTimeDistros(), LoadMIPtoGeV(), LoadPEACtoMIP(), LoadPEtoSigCor(), LoadSCtoPEAC(), MakeCellHit(), and SetContext().

std::vector<float> calib::Calibrator::fFTDistFDX

Definition at line 100 of file Calibrator.h. Referenced by Calibrator(), and LoadFineTimeDistros().

std::vector<float> calib::Calibrator::fFTDistFDY

Definition at line 101 of file Calibrator.h. Referenced by Calibrator(), and LoadFineTimeDistros().

std::vector<float> calib::Calibrator::fFTDistNDX

Definition at line 98 of file Calibrator.h. Referenced by Calibrator(), LoadFineTimeDistros(), and main().

std::vector<float> calib::Calibrator::fFTDistNDY

Definition at line 99 of file Calibrator.h. Referenced by Calibrator(), LoadFineTimeDistros(), and main().

Calibrator * Calibrator::fInstance = 0 [static, private]

Definition at line 22 of file Calibrator.cxx. Referenced by Instance().

bool calib::Calibrator::fLoadedFineTimeDistros [private]

Definition at line 104 of file Calibrator.h. Referenced by LoadFineTimeDistros().

std::vector<float> calib::Calibrator::fMIPtoGeV

Definition at line 93 of file Calibrator.h. Referenced by Calibrator(), and LoadMIPtoGeV().

std::vector<float> calib::Calibrator::fMIPtoGeVErr

Definition at line 94 of file Calibrator.h. Referenced by Calibrator().

float calib::Calibrator::fPEACperMIP

Definition at line 91 of file Calibrator.h. Referenced by LoadPEACtoMIP().

float calib::Calibrator::fPEACperMIPErr

Definition at line 92 of file Calibrator.h.

std::map<int,float> calib::Calibrator::fPEtoSigCor

Definition at line 96 of file Calibrator.h. Referenced by LoadPEtoSigCor().

std::map<int,attlen> calib::Calibrator::fSCtoPEAC

Definition at line 97 of file Calibrator.h. Referenced by LoadSCtoPEAC().

---

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

• Calibrator.h
• Calibrator.cxx

---