PEGainSummary.cxx

Go to the documentation of this file.

//#include <cstdio>
#include "PEGain/PEGainSummary.h"
#include "MessageService/MsgService.h"     // MSG text output

#include <cassert>
#include <cmath>

ClassImp(PEGainSummary)

CVSID("$Id: PEGainSummary.cxx,v 1.3 2006/05/27 09:49:04 rhatcher Exp $");

//......................................................................

PEGainSummary::PEGainSummary()
   :fGotNear(0),fGotFar(0),fEncodedStrip(-1)//,fBestSide(-1)
{

   MSG("PEGain",Msg::kVerbose) << "PEGainSummary::PEGainSummary()\n";

}
//......................................................................

PEGainSummary::~PEGainSummary()
{

   MSG("PEGain",Msg::kVerbose) << "PEGainSummary::~PEGainSummary()\n";

}
//......................................................................

void PEGainSummary::AddNearEnd(const PulserSummary &ps, const PlexStripEndId &seid)
{
   MSG("PEGain",Msg::kVerbose) << "PEGainSummary::AddNearEnd()\n";
   
   if(fGotNear==1) {
      MSG("PEGain",Msg::kWarning) << "I've already got a near end for : "
                                  << seid << endl;
   }
   else {
      if(fEncodedStrip==-1) {
         fSeid=seid;
         fEncodedStrip=seid.BuildPlnStripEndKey();
      }
      else assert(seid.BuildPlnStripEndKey()==fEncodedStrip);
      fNearEnd=ps;
      fGotNear=1;
   }
}

void PEGainSummary::AddFarEnd(const PulserSummary &ps, const PlexStripEndId &seid)
{
   MSG("PEGain",Msg::kVerbose) << "PEGainSummary::AddFarEnd()\n";
   
   if(fGotFar==1) {
      MSG("PEGain",Msg::kWarning) << "I've already got a far end for : "
                                  << seid << endl;
   }
   else {
      if(fEncodedStrip==-1) {
         fSeid=seid;
         fEncodedStrip=seid.BuildPlnStripEndKey();
      }
      else assert(seid.BuildPlnStripEndKey()==fEncodedStrip);
      fFarEnd=ps;
      fGotFar=1;
   }
}

Float_t PEGainSummary::GetNearGain() const
{
   MSG("PEGain",Msg::kVerbose) << "PEGainSummary::GetNearGain()\n";
   if(fGotNear) {
      Float_t tempRms=fNearEnd.GetRMS();
      Float_t tempMean=fNearEnd.GetMean();

      return (tempRms*tempRms/tempMean);
   }
   else return -1;
}

Float_t PEGainSummary::GetFarGain() const
{
   MSG("PEGain",Msg::kVerbose) << "PEGainSummary::GetFarGain()\n";
   if(fGotFar) {
      Float_t tempRms=fFarEnd.GetRMS();
      Float_t tempMean=fFarEnd.GetMean();

      return (tempRms*tempRms/tempMean);
   }
   return -1;
}

Float_t PEGainSummary::GetNearError() const
{
   MSG("PEGain",Msg::kVerbose) << "PEGainSummary::GetNearError()\n";
   if(fGotNear) {
      Float_t tempRms=fNearEnd.GetRMS();
      Float_t tempMean=fNearEnd.GetMean();
      Float_t tempNent=(float)fNearEnd.GetEntries();
      if(tempNent>1) {
         Float_t tempRmsErr=tempRms/sqrt(2.0*tempNent);
         Float_t tempMeanErr=tempRms/sqrt(tempNent-1);
         
         return sqrt(pow(2.0*tempRms/tempMean,2)*tempRmsErr*tempRmsErr
                     +pow(tempRms/tempMean,4)*tempMeanErr*tempMeanErr);
      }
      
   }
   return -1;
}

Float_t PEGainSummary::GetFarError() const
{
   MSG("PEGain",Msg::kVerbose) << "PEGainSummary::GetFarError()\n";
   if(fGotFar) {
      Float_t tempRms=fFarEnd.GetRMS();
      Float_t tempMean=fFarEnd.GetMean();
      Float_t tempNent=(float)fFarEnd.GetEntries();
      if(tempNent>1) {
         Float_t tempRmsErr=tempRms/sqrt(2.0*tempNent);
         Float_t tempMeanErr=tempRms/sqrt(tempNent-1);
         
         return sqrt(pow(2.0*tempRms/tempMean,2)*tempRmsErr*tempRmsErr
                     +pow(tempRms/tempMean,4)*tempMeanErr*tempMeanErr);
      }
      
   }
   return -1;
}

Int_t PEGainSummary::GetBest() const {
   MSG("PEGain",Msg::kVerbose) << "PEGainSummary::SetBest()\n";
   if(fGotNear) {
      if(fNearEnd.GetEntries()>200) {
         if(fGotFar) {
            if(fFarEnd.GetEntries()>200) {
               // Got two chances;
               Int_t nearNent=fNearEnd.GetEntries();
               Int_t farNent=fFarEnd.GetEntries();
               Float_t nentRatio=float(nearNent)/float(farNent);
               if(nentRatio>1.1) {
                  //fBestSide=1;
                  return 1;
               }
               if(nentRatio<0.9) {
                  //fBestSide=2;
                  return 2;
               }
               Float_t nearMean=fNearEnd.GetMean();
               Float_t farMean=fFarEnd.GetMean();
               if(nearMean<8000 && nearMean>800) {
                   return 1;
               }
               if(farMean<8000 && farMean>800) {
               }
               if(fabs(nearMean-5000.0)<fabs(farMean-5000.0)) {
                  //fBestSide=1;
                  return 1;
               }
               else {
                  //fBestSide=2;
                  return 2;
               }
            }
         }
         // Only Near
         //fBestSide=1;
         return 1;
      }
      
   }
   else if(fGotFar) {
      if(fFarEnd.GetEntries()>200) {
         //fBestSide=2;
         return 2;
      }
      // Only Far
      
   }
   //fBestSide=0;
   return 0;
}



Float_t PEGainSummary::GetBestGain() const {
   MSG("PEGain",Msg::kVerbose) << "PEGainSummary::GetBestGain()\n";
   Int_t fBestSide = GetBest();
   //   if(fBestSide==-1) SetBest();
   if(fBestSide==1) {
      return GetNearGain();
   }
   else if(fBestSide==2) {
      return GetFarGain();
   }
   return -1;
}


Float_t PEGainSummary::GetBestError() const
{
   MSG("PEGain",Msg::kVerbose) << "PEGainSummary::GetBestError()\n";
   Int_t fBestSide=GetBest();
   //   if(fBestSide==-1) SetBest();
   if(fBestSide==1) {
      return GetNearError();
   }
   else if(fBestSide==2) {
      return GetFarError();
   }
   return -1;
   
}
   
Float_t PEGainSummary::GetProperNearGain() const
{
   return 0.8*GetNearGain();
}
Float_t PEGainSummary::GetProperFarGain() const
{
   return 0.8*GetFarGain();
}
Float_t PEGainSummary::GetProperBestGain() const
{
   return 0.8*GetBestGain();
}
Float_t PEGainSummary::GetProperNearError() const 
{
   return 0.8*GetNearError();
}
Float_t PEGainSummary::GetProperFarError() const
{
   return 0.8*GetFarError();
}
Float_t PEGainSummary::GetProperBestError() const
{
   return 0.8*GetBestError();
}

---