CacherEM.cxx

Go to the documentation of this file.

#include "CacherEM.h"
#include "MessageService/MsgService.h"

CVSID("$Id: CacherEM.cxx,v 1.2 2004/09/29 00:07:29 cbs Exp $");

//********************************
CacherEM::CacherEM(){
  fn = 0;
  fpars = NULL;
  fpscache.clear();
  frtcache.clear();
  fsumE = 0;
}

//********************************
CacherEM::~CacherEM(){
  delete [] fpars;
}

//********************************
map<Int_t,spef> CacherEM::GetPSMap(){
  return fpscache;
}

//********************************
vector<rte> CacherEM::GetRTMap(){
  return frtcache;
}

Double_t CacherEM::GetSumE(){
  return fsumE;
}

//********************************
Double_t *CacherEM::GetPars(){
  return fpars;
}

//********************************
void CacherEM::StorePS(Int_t n,Double_t* pars,map<Int_t,spef> cache)
{  
  for(Int_t i=0;i<n;i++){
    fpars[i] = pars[i];
  }
  fpscache = cache;
}

//********************************
void CacherEM::StoreRT(Int_t n,Double_t* pars,vector<rte> cache,Double_t sumE)
{
  fn = n;
  fpars = new Double_t[n];
  for(Int_t i=0;i<n;i++){
    fpars[i] = pars[i];
  }
  frtcache = cache;
  fsumE = sumE;
}

//********************************
Bool_t CacherEM::ValidPS(Double_t* pars){

  Bool_t valid = true;
  if(fn==0) return false;
  for(Int_t i=0;i<fn;i++){
    if(pars[i]!=fpars[i]) {
      valid = false;
      break;
    }
  }
  return valid;
}

//********************************
Bool_t CacherEM::ValidRT(Double_t* pars){

  if(fn==0) return false;
  if(pars[0]!=fpars[0]) return false;
  return true;

}

//********************************
void CacherEM::PrintPSMap()
{

  map<Int_t,spef>::iterator beg = fpscache.begin();
  map<Int_t,spef>::iterator end = fpscache.end();

  while(beg!=end){
    MSG("FitShowerEM", Msg::kVerbose) << beg->first << " " 
                                      << beg->second.en << " " 
                                      << beg->second.fl << " " 
                                      << endl;
    beg++;
  }
  
}

//********************************
void CacherEM::PrintRTMap()
{

  vector<rte>::iterator beg = frtcache.begin();
  vector<rte>::iterator end = frtcache.end();

  while(beg!=end){
    MSG("FitShowerEM", Msg::kVerbose) << beg->t << " " 
                                      << beg->r << " " 
                                      << beg->e << " " 
                                      << beg->f << " " 
                                      << endl;
    beg++;
  }
  
}

//********************************
void CacherEM::ClearCache(){
  
  fn = 0;
  Double_t *temp = fpars;
  fpars = NULL;
  delete [] temp;
  fpscache.clear();
  frtcache.clear();

}

//********************************
void CacherEM::ClearPSCache(){
  
  fpscache.clear();

}

---