Gryfindor/detector/channel.cxx

Go to the documentation of this file.

/***************************************************************************
                          channel.cpp  -  description
                             -------------------
    begin                : Fri Feb 8 2002
    copyright            : (C) 2002 by Robert Morse
    email                : r.morse@sussex.ac.uk
 ***************************************************************************/


#include "channel.h"
void Channel::Clean(){
  no_entries=0;
  mean=0;
  rms=0;
  sigma1=0;
  sigma2=0;
  channel_no =-1;
};

void Channel::SetMean(Double_t value){
  mean = value;
};

void Channel::SetEntries(Int_t values){
  no_entries = values;
};

void Channel::SetRMS(Double_t value){
  rms = value;
};

Double_t Channel::GetMean(){
  return mean;
};

Int_t Channel::GetEntries(){
  return no_entries;
};

Double_t Channel::GetRMS(){
        if(sigma1 ==0)
                return rms;
         return   rms = pow(( sigma1/Double_t(no_entries) - pow( sigma2/Double_t(no_entries),2)),0.5);
};

void Channel::SetEntry(Int_t entries, Double_t new_mean, Double_t new_rms){
  no_entries = entries;
  mean = new_mean;
  rms = new_rms;
};

void Channel::AddEntry(Int_t new_entries, Double_t new_mean, Double_t new_rms){
  if(new_entries > 0){

         sigma1 += Double_t(new_entries)*(pow(new_rms,2) + pow(new_mean,2));
         sigma2 += new_mean * Double_t(new_entries);

    mean = mean * ( Double_t(no_entries) / (Double_t(no_entries + new_entries)) )
    + new_mean*( Double_t(new_entries) / (Double_t(new_entries+no_entries)) );
        no_entries += new_entries;
        }
};

---