PulserSummary.cxx

Go to the documentation of this file.

// $Id: PulserSummary.cxx,v 1.8 2007/03/01 17:06:40 rhatcher Exp $
//
// PulserSummary.cxx
//
// List of averages of pulser summaries
//
// Author: Phil Adamson 4/2001
//
// $Log: PulserSummary.cxx,v $
// Revision 1.8  2007/03/01 17:06:40  rhatcher
// fix a bunch of CVSID lines.  The quoted string must be "$Id: $" (especially
// common is to leave off the trailing $).
//
// Revision 1.7  2004/09/08 21:50:51  pa
// Checkin of new PulserCalibration version. MIght even work
//
// Revision 1.6  2003/02/06 20:02:59  rjn
// Modified PulserSummary.h and added classes to make PulserSummary objects from Raw Data.
//
// Revision 1.5  2002/10/11 10:38:12  phil.adamson
// Add MakeLin, some minor fixups
//
// Revision 1.4  2002/10/07 10:07:43  phil.adamson
// Try to fix the occasional nan...
//
// Revision 1.3  2002/09/25 08:54:08  phil.adamson
// Add Ryan's bits to make PulserSummary objects from Raw Data
//
// Revision 1.2  2002/08/07 16:16:52  phil.adamson
// Minor cosmetic details, plus provide ZeroCorr() methods in PulserSummary.
// Phil.
//
//

#include "MessageService/MsgService.h"
#include "PulserCalibration/PulserSummary.h"
#include "RawData/RawLIAdcSummary.h"
#include <cmath>
#include <cassert>

CVSID("$Id: PulserSummary.cxx,v 1.8 2007/03/01 17:06:40 rhatcher Exp $");

ClassImp(PulserSummary)
  
  // ............

PulserSummary::PulserSummary()
{
  MSG("Pulser",Msg::kVerbose) << "PulserSummary ctor\n";
  fMean=0;
  fMeanSquare=0;
  fRMS=0;
  fNumEntries=0;
  fChannel=0; // empty object has channel 0
}

PulserSummary::~PulserSummary()
{

}

void PulserSummary::Add(RawLIAdcSummary &from, const PlexPinDiodeId &pinid)
{
  fPinId = pinid;
  fReadoutType = ReadoutType::kPinDiode;
  MSG("Pulser",Msg::kVerbose) << "PIN "<<pinid<<" has "<<fMean<<" "<<fNumEntries<<" "<<fNumTriggers<<endl;
  Add(from);
}

void PulserSummary::AddDigit(RawDigit &from, 
                             const PlexPinDiodeId &pinid)
{
   fPinId = pinid;
   fReadoutType = ReadoutType::kPinDiode;
   MSG("Pulser",Msg::kVerbose) << "PIN "<<pinid<<" has "<<fMean<<" "<<fNumEntries<<" "<<fNumTriggers<<endl;
   AddDigit(from);
}

void PulserSummary::Add(RawLIAdcSummary &from, const PlexStripEndId &seid)
{
  fSEId = seid;
  fReadoutType = ReadoutType::kScintStrip;
  Add(from);
}

void PulserSummary::AddDigit(RawDigit &from, const PlexStripEndId &seid)
{
  fSEId = seid;
  fReadoutType = ReadoutType::kScintStrip;
  AddDigit(from);
}

void PulserSummary::AddDigit(PulserNearDigit &from, const PlexStripEndId &seid)
{
  fSEId = seid;
  fReadoutType = ReadoutType::kScintStrip;
  AddDigit(from);
}

void PulserSummary::Add(RawLIAdcSummary &from)
{
   MSG("Pulser",Msg::kVerbose) << "fChannel " << RawChannelId(fChannel)
                               << " and Channel is: " << from.GetChannel()
                               << endl; 
   if (fChannel==0) {
      fChannel=from.GetChannel().GetEncoded();
   } 
   MSG("Pulser",Msg::kVerbose) << "fChannel " << RawChannelId(fChannel)
                               << " and Channel is: " << from.GetChannel()
                              << endl; 
  assert(from.GetChannel().GetEncoded()==fChannel);
  Int_t tempEnd=(int)fSEId.GetEnd();
  if(fSEId.GetPlane()==10 && fSEId.GetStrip()==4 && tempEnd==1) {
     MSG("Pulser",Msg::kVerbose) << "My Channel " 
                                 << "\tMean: " << from.GetMean()
                                 << "\tRMS: " <<  from.GetRms()
                                 << "\tEntries: " << from.GetEntries() << endl;
  }

 
    MSG("Pulser",Msg::kVerbose) << "PulserSummary::Add() Got channel "<<fChannel<<endl;
    fMean = 
      fNumEntries*(1.0/(fNumEntries + from.GetEntries())) * fMean 
      + from.GetEntries()*(1.0/(fNumEntries + from.GetEntries())) 
      * from.GetMean();
    fMeanSquare = 
      fNumEntries*(1.0/(fNumEntries + from.GetEntries())) *fMeanSquare
      + from.GetEntries()*(1.0/(fNumEntries + from.GetEntries())) 
      * (float(from.GetRms())*from.GetRms()+float(from.GetMean())*from.GetMean());
    fNumEntries +=from.GetEntries();
}

void PulserSummary::AddDigit(RawDigit &from)
{
   MSG("Pulser",Msg::kVerbose) << "fChannel " << RawChannelId(fChannel)
                               << " and Channel is: " << from.GetChannel()
                               << endl; 
   if (fChannel==0) {
      fChannel=from.GetChannel().GetEncoded();
   } 
   MSG("Pulser",Msg::kVerbose) << "fChannel " << RawChannelId(fChannel)
                               << " and Channel is: " << from.GetChannel()
                               << endl; 
   assert(from.GetChannel().GetEncoded()==fChannel);
   //Int_t tempEnd=(int)fSEId.GetEnd();
   MSG("Pulser",Msg::kVerbose) 
      << "PulserSummary::Add() Got channel "<<fChannel<<endl;
   fMean = (float(fNumEntries)/(float(fNumEntries)+1.0))*fMean +
      (1.0/(float(fNumEntries)+1.0))*from.GetADC();
   fMeanSquare = (float(fNumEntries)/(float(fNumEntries)+1.0))*fMeanSquare +
      (1.0/(float(fNumEntries)+1.0))*(from.GetADC()*from.GetADC());
   fNumEntries ++;
}

void PulserSummary::AddDigit(PulserNearDigit &from)
{
   MSG("Pulser",Msg::kVerbose) << "fChannel " << RawChannelId(fChannel)
                               << " and Channel is: " << from.GetChannel()
                               << endl; 
   if (fChannel==0) {
      fChannel=from.GetChannel().GetEncoded();
   } 
   MSG("Pulser",Msg::kVerbose) << "fChannel " << RawChannelId(fChannel)
                               << " and Channel is: " << from.GetChannel()
                               << endl; 
   assert(from.GetChannel().GetEncoded()==fChannel);
   //Int_t tempEnd=(int)fSEId.GetEnd();
   MSG("Pulser",Msg::kVerbose) 
      << "PulserSummary::Add() Got channel "<<fChannel<<endl;
   fMean = (float(fNumEntries)/(float(fNumEntries)+1.0))*fMean +
      (1.0/(float(fNumEntries)+1.0))*from.GetADC();
   fMeanSquare = (float(fNumEntries)/(float(fNumEntries)+1.0))*fMeanSquare +
      (1.0/(float(fNumEntries)+1.0))*(from.GetADC()*from.GetADC());
   fNumEntries ++;
}

void PulserSummary::Finish(int numTriggers)
{
   
  if ((fMeanSquare - fMean*fMean)>1) {
    fRMS = std::sqrt(fMeanSquare - fMean*fMean);
  }
  else {
    fRMS = 1;
  }
   fNumTriggers = numTriggers;
   fMeanZC = fMean * float(fNumEntries)/float(numTriggers);
   float ms = fMeanSquare * float(fNumEntries)/float(numTriggers) -
     fMeanZC*fMeanZC;
   if (ms>1) {
     fRMSZC = std::sqrt(ms);
   }
   else {
     fRMSZC = 1;
   }
   MSG("Pulser",Msg::kVerbose) 
      << "PulserSummary::Finish() finished summary in channel: "
      << fChannel << " Mean:  " << fMean << " RMS: " 
      << fRMS << " NumEntries: " << fNumEntries << endl;
}





void PulserSummary::SetPinDiodeId(PlexPinDiodeId pin)
{
  fPinId = pin;

}

---