HarpoDccMap.cxx

Go to the documentation of this file.

// 
// File HarpoDccMap.cxx
//
#include "HarpoDccMap.h"
#include <iostream>

ClassImp(HarpoDccMap) 

HarpoDccMap::HarpoDccMap() :  map(NALL,NADC) , fEmptyCell(-1.0)  
{
  map = fEmptyCell;
  fPedestal = NULL;
}

// void  HarpoDccMap::SetEmptyCell(Short_t val)
// {
//   fEmptyCell = (Double_t) val;
// }

void HarpoDccMap::SetEmptyCell(Double_t val)
{
  fEmptyCell = val;
}


// void HarpoDccMap::SetPedestal(UShort_t ped)
// {

//   for(Int_t i = 0; i<NALL; i++)
//     fPedestal[i] = (Double_t) ped;
// }

// void HarpoDccMap::SetPedestal(Int_t i, Double_t ped){
//   if(i<0 || i>=NALL) Warning("SetPedestal","index out of bounds (%d)",i);
//   else fPedestal[i] = ped;
// }

void HarpoDccMap::SetPedestal(Double_t ped)
{
    fPedestal->SetMean(ped);
}

Double_t HarpoDccMap::GetPedestal(Int_t i){
  if(i<0 || i>=NALL){
    Warning("GetPedestal","index out of bounds (%d)",i);
    return -1;
  }else return fPedestal->GetPed((UInt_t) i);
}


void HarpoDccMap::SetSigmaNoise(Double_t val)
{
  fPedestal->SetSigma(val);
}

// void HarpoDccMap::SetSigmaNoise(Int_t i, Double_t val){
//   if(i<0 || i>=NALL) Warning("SetSigmaNoise","index out of bounds (%d)",i);
//   else fSigmaNoise[i] = val;
// }

Double_t HarpoDccMap::GetSigmaNoise(Int_t i){
  if(i<0 || i>=NALL){
    Warning("GetSigmaNoise","index out of bounds (%d)",i);
    return -1;
  }else return  fPedestal->GetSigma((UInt_t) i);
}

Double_t HarpoDccMap::GetRawData(Int_t i, Int_t j){ 
  if(i>=0 && i<NALL && j>=0 && j<NADC) {
    return map(i,j);
  } else {
    Warning("GetData","index out of matix boundd (%d,%d)",i,j);
    return -2000.0;
  }
}

Double_t HarpoDccMap::GetData(Int_t i, Int_t j){ 
  if(i>=0 && i<NALL && j>=0 && j<NADC) {
    if(map(i,j) == -1) return -1000.0; 
    else return map(i,j) - fPedestal->GetPed((UInt_t) i);
  } else {
    Warning("GetData","index out of matix boundd (%d,%d)",i,j);
    return -2000.0;
  }
}

void HarpoDccMap::SetRawData(Int_t i, Int_t j, Double_t val){ 
  if(i>=0 && i<NALL && j>=0 && j<NADC) {
    map(i,j) = val;
  } else {
    Warning("SetRawData","index out of matix boundd (%d,%d)",i,j);
  }
}

void HarpoDccMap::SetData(Int_t i, Int_t j, Double_t val){ 
  if(i>=0 && i<NALL && j>=0 && j<NADC) {
    map(i,j) = val + fPedestal->GetPed((UInt_t) i);
  } else {
    Warning("SetData","index out of matix boundd (%d,%d)",i,j);
  }
}

void HarpoDccMap::FillIf(HarpoDccMap *m, Int_t row, Int_t col, UShort_t word)
{
  //  std::cout << "row " << row << " col " << col << " val " << word << std::endl;
  if ( ( row >= 0 ) && (row < m->map.GetNrows() ) 
       && (col >= 0 ) && (col < m->map.GetNcols())) {
    //    m->map(row,col) = (Double_t) word;
    m->SetRawData(row,col,(Double_t) word);
  } 
  else
 {
   std::cout << "---DccMap--- Bad Row Col " << row << " " << col << std::endl;
 }    
}

void          HarpoDccMap::print()
{
  std::cout << "---DccMap--- : fEmptyCell " << fEmptyCell  << std::endl;
  //  map.Print("f=  %6.2f  ");
  map.Print(); // "%11.4g"
}