#include "hades.h"
#include "hdetpario.h"
#include "hpario.h"
#include "hspectrometer.h"
#include "hstart2calpar.h"
#include "hstart2detector.h"
#include <iomanip>
#include <iostream>
using namespace std;
ClassImp(HStart2Calpar)
ClassImp(HStart2CalparMod)
ClassImp(HStart2CalparCell)
HStart2CalparMod::HStart2CalparMod(Int_t n)
{
array = new TObjArray(n);
for (Int_t i = 0 ; i < n ; i++) {
array->AddAt(new HStart2CalparCell(), i);
}
}
HStart2CalparMod::~HStart2CalparMod()
{
array->Delete();
delete array;
}
HStart2Calpar::HStart2Calpar(const Char_t* name, const Char_t* title,
const Char_t* context)
: HParSet(name, title, context)
{
strcpy(detName, "Start");
HStart2Detector* det = (HStart2Detector*)(gHades->getSetup()->getDetector("Start"));
Int_t nMod = det->getMaxModInSetup();
Int_t nComp = det->getMaxComponents();
array = new TObjArray(nMod);
for (Int_t i = 0; i < nMod; i++) {
if (det->getModule(-1, i) > 0) {
array->AddAt(new HStart2CalparMod(nComp), i);
} else {
array->AddAt(new HStart2CalparMod(0), i);
}
}
}
HStart2Calpar::~HStart2Calpar()
{
array->Delete();
delete array;
}
Bool_t HStart2Calpar::init(HParIo* inp, Int_t* set)
{
HDetParIo* input = inp->getDetParIo("HStartParIo");
if (input) return (input->init(this, set));
return kFALSE;
}
Int_t HStart2Calpar::write(HParIo* output)
{
HDetParIo* out = output->getDetParIo("HStartParIo");
if (out) return out->write(this);
return -1;
}
void HStart2Calpar::clear()
{
for (Int_t i = 0; i < getSize(); i++) {
HStart2CalparMod& mod = (*this)[i];
for (Int_t j = 0; j < mod.getSize(); j++)
mod[j].clear();
}
status = kFALSE;
resetInputVersions();
}
void HStart2Calpar::printParam()
{
printf("Calibration parameters for the START2 detector\n");
printf("module strip TdcSlope TdcOffset AdcSlope AdcOffset\n");
Float_t data[4];
for (Int_t i = 0; i < getSize(); i++) {
HStart2CalparMod& mod = (*this)[i];
for (Int_t j = 0; j < mod.getSize(); j++) {
HStart2CalparCell& cell = mod[j];
cell.getData(data);
if (data[0] != 1.F || data[1] != 0.F || data[2] != 1.F || data[3] != 0.) {
printf("%4i %4i %10.5f %10.3f %10.5f %10.3f\n",
i, j, data[0], data[1], data[2], data[3]);
}
}
}
}
void HStart2Calpar::readline(const Char_t *buf, Int_t *set)
{
Int_t mod, strip;
Float_t data[4] = {0.F, 0.F, 0.F, 0.F};
sscanf(buf, "%i%i%f%f%f%f", &mod, &strip, &data[0], &data[1], &data[2], &data[3]);
if (!set[mod]) return;
(*this)[mod][strip].fill(data);
set[mod] = 999;
}
void HStart2Calpar::putAsciiHeader(TString& header)
{
header =
"# Calibration parameters for the START2 detector\n"
"# Format:\n"
"# module strip TdcSlope TdcOffset AdcSlope AdcOffset\n";
}
Bool_t HStart2Calpar::writeline(Char_t *buf, Int_t mod, Int_t strip)
{
HStart2CalparCell& cell = (*this)[mod][strip];
Float_t data[4];
cell.getData(data);
if (data[0] != 1.F || data[1] != 0.F || data[2] != 1.F || data[3] != 0.F) {
sprintf(buf, "%4i %4i %10.5f %10.3f %10.5f %10.3f\n",
mod, strip, data[0], data[1], data[2], data[3]);
return kTRUE;
}
return kFALSE;
}
Bool_t HStart2Calpar::isInTable(Int_t mod,Int_t strip, Bool_t silent)
{
Int_t maxMod = getSize()-1 ;
if(mod>maxMod){
if(!silent) Warning("isInTable()","Module number out of range! module = %i (max=%i)",mod,maxMod);
return kFALSE;
}
HStart2CalparMod* pMod = (&(*this)[mod]);
if(!pMod){
if(!silent) Warning("isInTable()","Module number not in HStart2Calpar ! module = %i",mod);
return kFALSE;
}
Int_t maxStrip = (*this)[mod].getSize()-1;
if(maxStrip==-1){
if(!silent) Warning("isInTable()","Module number not in HStart2Calpar ! module = %i",mod);
return kFALSE;
}
if(strip<0) return kTRUE;
if(strip>maxStrip){
if(!silent) Warning("isInTable()","Strip number out of range! module = %i , strip = %i (max=%i)",mod,strip,maxStrip);
return kFALSE;
}
return kTRUE;
}