PulserFromRawSummaryList.cxx

Go to the documentation of this file.

// $Id: PulserFromRawSummaryList.cxx,v 1.3 2003/09/03 09:31:16 rjn Exp $
//
// PulserFromRawSummaryList.cxx
//
// List of averages of pulser summaries
//
// Author: Ryan Nichol
//
//


#include "TNamed.h"
#include "MessageService/MsgService.h"
#include "Validity/VldContext.h"
#include "PulserCalibration/PulserFromRawSummaryList.h"
#include "PulserCalibration/PulserSummary.h"
#include "RawData/RawLIAdcSummaryBlock.h"
#include "RawData/RawLIAdcSummary.h"
#include "RawData/RawDigit.h"
#include "Plex/PlexHandle.h"
#include "Plex/PlexusABC.h"
#include "RawData/RawChannelId.h"
#include <map>
#include <set>
#include <algorithm>
#include <cassert>

ClassImp(PulserFromRawSummaryList)
  
  // ............

CVSID("$Id: PulserFromRawSummaryList.cxx,v 1.3 2003/09/03 09:31:16 rjn Exp $");

PulserFromRawSummaryList::PulserFromRawSummaryList() : 
  TNamed("PulserFromRawSummaryList","PulserFromRawSummaryList")
{
  MSG("Pulser",Msg::kVerbose) <<" PulserFromRawSummaryList ctor \n";
  fLowPin = new PulserSummary();
  fHighPin = new PulserSummary();
  fTriggers=0;
  fUsePlex = false;
}

PulserFromRawSummaryList::~PulserFromRawSummaryList()
{
  MSG("Pulser",Msg::kVerbose) <<" PulserFromRawSummaryList dtor \n";
  delete fLowPin;
  delete fHighPin;
}


PulserFromRawSummaryList::PulserFromRawSummaryList(const RawDigitDataBlock &rawddb, 
                                             const PulserFromRawLiBlock &rawsum) : 
  TNamed("PulserFromRawSummaryList","PulserFromRawSummaryList")
{
  MSG("Pulser",Msg::kVerbose) <<
    " PulserFromRawSummaryList constructing from RawDigitDataBlock \n";
  fLowPin = new PulserSummary();
  fHighPin = new PulserSummary();
  fTriggers=0;
  fNumPulses = rawsum.GetPulses();
  fLed=rawsum.GetLed();
  fPulserBox=rawsum.GetPulserBox();
  fHeight=rawsum.GetPulseHeight();
  fWidth=rawsum.GetPulseWidth();
  fPointIndex=rawsum.GetCalibPoint();
  fExpectedPoints=rawsum.GetCalibType(); // Will be 1 for dougpoint, 
                                         // more than 1 for gaincurve
  fAggregateNo = PlexLedId(rawsum.GetVldContext().GetDetector(),
                           fPulserBox, fLed);
  fIdent = fLed + 100 * fExpectedPoints;
  fUsePlex = false;
  fValidity = rawsum.GetVldContext();
  this->Add(rawddb,rawsum);
}

VldContext PulserFromRawSummaryList::GetVldContext() const
{
  return fValidity;
}



void PulserFromRawSummaryList::Add(const RawDigitDataBlock &rawddb, 
                                const PulserFromRawLiBlock &rawsum)
{
  // If we get here, we should be adding the same kinds of pulse
  // but just in case...
  MSG("Pulser",Msg::kVerbose) <<"LED "<<fLed<<"  "<<rawsum.GetLed()<<endl;
  MSG("Pulser",Msg::kVerbose) <<"H   "<<fHeight<<"  "<<rawsum.GetPulseHeight()<<endl;
  MSG("Pulser",Msg::kVerbose) <<"W   "<<fWidth<<"  "<<rawsum.GetPulseWidth()<<endl;

  assert(fLed==rawsum.GetLed());
  assert(fPulserBox==rawsum.GetPulserBox());
  assert(fHeight==rawsum.GetPulseHeight());
  assert(fWidth==rawsum.GetPulseWidth());
  assert(fPointIndex==rawsum.GetCalibPoint());
  assert(fExpectedPoints==rawsum.GetCalibType());

  //Get validity context
  VldContext vldc(rawsum.GetVldContext());
  PlexHandle plex(vldc);

  int numTriggers=0;
  std::pair<PlexPinDiodeId,PlexPinDiodeId> myPins = 
     plex.GetPinDiodeIds(fAggregateNo); // hi, low
  MSG("Pulser",Msg::kVerbose)<< "We have "
                         << rawddb.GetNumberOfDigits() 
                         << " digits"<<endl;
  TIter rawit = rawddb.GetDatumIter();

  uIntNearMap tempNearMap;
  uIntNearMap::iterator tempNearIter;
  while (RawDigit *rd = (RawDigit *) rawit()) {
     
     RawChannelId rcid = rd->GetChannel(); 
     ReadoutType::Readout_t rt = plex.GetReadoutType(rcid);
     switch (rcid.GetElecType()) {
         case (ElecType::kQIE):
             if(tempNearMap.find(rcid.GetEncoded())==
                tempNearMap.end()) {
                 tempNearMap[rcid.GetEncoded()].New(*rd);
             }
             else {
                 tempNearIter=tempNearMap.find(rcid.GetEncoded());
                 tempNearIter->second.Add(*rd);
             }
             break;
         case (ElecType::kVA): 
             if (rt & ReadoutType::kScintStrip) {
                 MSG("Pulser",Msg::kVerbose) 
                     << "Got LI in strip at "<< rcid<<endl;
                 
                 // Find out which strip-end we are. At CalDet, this is
                 // a bit pointless, but...
                 PlexSEIdAltL altlist = plex.GetSEIdAltL(rcid);
                 for (PlexSEIdAltL::iterator it = altlist.begin(); 
                      it!=altlist.end(); it++) {
                     PlexStripEndId seid = (*it).GetSEId();
                     if (plex.GetLedId(seid) == fAggregateNo) {
                         MSG("Pulser",Msg::kVerbose) << "Got LI in strip "
                                                   << seid << endl;
            
                         fNearMap[seid.BuildPlnStripEndKey()].AddDigit(*rd,seid);
                     }
                     else if (plex.GetLedId(seid.GetOppositeSEId())==
                              fAggregateNo) {
                         MSG("Pulser",Msg::kVerbose) << "Got LI in strip "
                                                   << seid << endl;
                         
                         fFarMap[seid.BuildPlnStripEndKey()].AddDigit(*rd,seid);
                     }
                     else {
                         // Crosstalk or such
                         //     MSG("Pulser",Msg::kVerbose) << "Xtalk "<<rcid<<endl; 
                     }
                 }
             }
             else if (rt & ReadoutType::kPinDiode) {
                 MSG("Pulser",Msg::kVerbose) << "Got LI in PIN at "
                                           << rd->GetChannel()<< endl;
                 
                 PlexPinDiodeId pinid = plex.GetPinDiodeId(rcid);
                 if (pinid==myPins.first) {
                     //High gain PIN
                     fHighPin->AddDigit(*rd,pinid);
                 }
                 else if (pinid==myPins.second) {
                     //Low Gain PIN
                     fLowPin->AddDigit(*rd,pinid);
                 }
                 // Else not this LED's PIN
             }
             else if (rt & ReadoutType::kFlashTrigPMT) {
                 MSG("Pulser",Msg::kVerbose) 
                     << "Got LI in trigger at "<<rd->GetChannel()<<endl;
                 numTriggers++;
             }
             
             break;
         default:
             break;
     }
     
  }
  for(tempNearIter = tempNearMap.begin(); 
      tempNearIter != tempNearMap.end(); ++tempNearIter) {

      PulserNearDigit nd = tempNearIter->second;
      RawChannelId rcid = nd.GetChannel();
 
      ReadoutType::Readout_t rt2 = plex.GetReadoutType(rcid); 
      if (rt2 & ReadoutType::kScintStrip) {
          MSG("Pulser",Msg::kVerbose) 
              << "Got LI in strip at "<< rcid<<endl;
          
          // Find out which strip-end we are. At CalDet, this is
          // a bit pointless, but...
          PlexSEIdAltL altlist = plex.GetSEIdAltL(rcid);
          for (PlexSEIdAltL::iterator it = altlist.begin(); 
               it!=altlist.end(); it++) {
              PlexStripEndId seid = (*it).GetSEId();
              if (plex.GetLedId(seid) == fAggregateNo) {
                  MSG("Pulser",Msg::kVerbose) << "Got LI in strip "
                                            << seid << endl;
                  
                  fNearMap[seid.BuildPlnStripEndKey()].AddDigit(nd,seid);
              }
              else if (plex.GetLedId(seid.GetOppositeSEId())==
                       fAggregateNo) {
                  MSG("Pulser",Msg::kVerbose) << "Got LI in strip "
                                            << seid << endl;
                  
                  fFarMap[seid.BuildPlnStripEndKey()].AddDigit(nd,seid);
              }
              else {
                  // Crosstalk or such
                  //    MSG("Pulser",Msg::kVerbose) << "Xtalk "<<rcid<<endl; 
              }
          }
      }
     
  }
  
  fTriggers +=numTriggers;
  MSG("Pulser",Msg::kDebug) << "Trig PMT: "<<numTriggers<<" of "<<fTriggers<<" triggers\n";
  
}

void PulserFromRawSummaryList::Finish()
{
  // Do zero correction here.
  for (std::map<Int_t,PulserSummary>::iterator it = fNearMap.begin();
       it!=fNearMap.end(); it++) (*it).second.Finish(fTriggers);
  for (std::map<Int_t,PulserSummary>::iterator it = fFarMap.begin();
       it!=fFarMap.end(); it++) (*it).second.Finish(fTriggers);
  fLowPin->Finish(fTriggers);
  fHighPin->Finish(fTriggers);
}

---