PulserSummaryList.cxx

Go to the documentation of this file.

// $Id: PulserSummaryList.cxx,v 1.7 2006/04/19 20:00:55 rhatcher Exp $
//
// PulserSummaryList.cxx
//
// List of averages of pulser summaries
//
// Author: Phil Adamson 4/2001
//
// $Log: PulserSummaryList.cxx,v $
// Revision 1.7  2006/04/19 20:00:55  rhatcher
// convert DetectorType:: to Detector:: everywhere.
//
// Revision 1.6  2004/11/30 07:56:25  rhatcher
// react to modified RawChannelId ctor
//
// Revision 1.5  2004/09/08 21:50:51  pa
// Checkin of new PulserCalibration version. MIght even work
//
// Revision 1.4  2004/09/03 17:31:12  phil.adamson
//
//
// Fix breakage introduced by change to RawLIAdcSummaryBlock.h
//
// Revision 1.3  2002/10/11 10:38:12  phil.adamson
// Add MakeLin, some minor fixups
//
// Revision 1.2  2002/08/07 16:16:52  phil.adamson
// Minor cosmetic details, plus provide ZeroCorr() methods in PulserSummary.
// Phil.
//
//


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

ClassImp(PulserSummaryList)
  
  // ............

CVSID("$Id: PulserSummaryList.cxx,v 1.7 2006/04/19 20:00:55 rhatcher Exp $");

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

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


PulserSummaryList::PulserSummaryList(const RawLIAdcSummaryBlock & rawsum) : 
  TNamed("PulserSummaryList","PulserSummaryList")
{
  MSG("Pulser",Msg::kVerbose) <<
    " PulserSummaryList constructing from RawLIAdcSummary \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(rawsum);
}

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



void PulserSummaryList::Add(const RawLIAdcSummaryBlock &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
  fHighPin->SetPinDiodeId(myPins.first);
  fLowPin->SetPinDiodeId(myPins.second);
  //myGetPinDiodeIds(plex,fAggregateNo);

  MSG("Pulser",Msg::kVerbose)<<"We have "<<rawsum.GetNumberOfSummaries() <<" summaries"<<endl;
  TIter rawit = rawsum.GetDatumIter(); 
  while (RawLIAdcSummary *rawch = 
         dynamic_cast<RawLIAdcSummary *>(rawit.Next())) {
      RawChannelId rcid = rawch->GetChannel();
      ReadoutType::Readout_t rt = plex.GetReadoutType(rcid);
      if (rawsum.GetNumberOfSummaries()<10) 
        MSG("Pulser",Msg::kVerbose) << "Channel "<<rcid<<endl; 
      if (rt & ReadoutType::kScintStrip) {
        MSG("Pulser",Msg::kVerbose) << "Got LI in strip at "<<rawch->GetChannel()<<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;
             // near end
             // MSG("Pulser",Msg::kVerbose) << "Near "<<rcid<<endl;
             fNearMap[seid.BuildPlnStripEndKey()].Add(*rawch,seid);
          }
          else if (plex.GetLedId(seid.GetOppositeSEId())==fAggregateNo) {
             MSG("Pulser",Msg::kVerbose) << "Got LI in strip "
                                         << seid << endl;
             // far end
             // MSG("Pulser",Msg::kVerbose) << "Far "<<rcid<<endl; 
             fFarMap[seid.BuildPlnStripEndKey()].Add(*rawch,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 "<<rawch->GetChannel()<<endl;
         PlexPinDiodeId pinid = plex.GetPinDiodeId(rcid);
         //     MSG("Pulser",Msg::kDebug)<<"PIN "<<rcid<<" is "<<pinid<<endl;
         if (pinid==myPins.first) {
            //High gain PIN
            //MSG("Pulser",Msg::kDebug)<<"High Gain PIN "<<rcid<<" is "<<pinid<<endl;
            fHighPin->Add(*rawch,pinid);
         }
         else if (pinid==myPins.second) {
            //Low Gain PIN
            //MSG("Pulser",Msg::kDebug)<<"Low Gain PIN "<<rcid<<" is "<<pinid<<endl;
            fLowPin->Add(*rawch,pinid);
         }
         //MSG("Pulser",Msg::kDebug)<<"Other PIN "<<pinid<<" with "<<rawch->GetEntries()<<endl;
         // Else not this LED's PIN
      }
      else if (rt & ReadoutType::kFlashTrigPMT) {
         MSG("Pulser",Msg::kVerbose) << "Got LI in trigger at "<<rawch->GetChannel()<<endl;
         numTriggers = max(numTriggers,int(rawch->GetEntries()));
      }
  }
  
  
  fTriggers +=numTriggers;
  MSG("Pulser",Msg::kDebug) << "Trig PMT: "<<numTriggers<<" of "<<fTriggers<<" triggers\n";
  
}

void PulserSummaryList::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);
}


//  std::pair<PlexPinDiodeId,PlexPinDiodeId> PulserSummaryList::myGetPinDiodeIds(PlexHandle plex, PlexLedId led)
//  {
//     Int_t highpinchads[6]={4146,4274,4402,4786,4914,5042};
//     Int_t lowpinchads[6]={4114,4242,4370,4626,4882,5010};
//     PlexPinDiodeId high = 
//       plex.GetPinDiodeId(RawChannelId(Detector::kCalDet,ElecType::kVA,
//                                   0,highpinchads[led.GetLedInBox()-1]));
//     PlexPinDiodeId low = 
//       plex.GetPinDiodeId(RawChannelId(Detector::kCalDet,ElecType::kVA,
//                                   0,lowpinchads[led.GetLedInBox()-1]));
//      return std::make_pair(high,low);
//  }

//  PlexLedId PulserSummaryList::myGetLedId(PlexHandle plex, PlexStripEndId seid)
//  {
//    int led = getLed(plex.GetRawChannelId(seid).GetCrate(),seid.GetStrip(),seid.GetPlane());
//    return PlexLedId(Detector::kCalDet, 0 , led);

//  }

//  int PulserSummaryList::getLed(int crate,int strip, int plane)
//  {
//    switch (crate)
//      {
//      case 0:
//        if (strip<4) {
//      if (plane==40) return 2;
//      if ((plane%4)==0||(plane%4)==1) return 1;
//      return 2;
//        } else if (strip<14) {
//      if (plane==50) return 3;
//      if (plane==40) return 1;
//      if ((plane%4)==0||(plane%4)==1) return 2;
//      return 1;
//        } else {
//      if (plane==50) return 1;
//      return 3;
//        }
//      case 1:
//        if (strip<10) {
//      if (plane==56||plane==58) return 4;
//      return 6;
//        } else if (strip<20) {
//      if (plane==58) return 6;
//      if ((plane%4)==0||(plane%4)==1) return 5;
//      return 4;
//        } else {
//      if ((plane%4)==0||(plane%4)==1) return 4;
//      return 5;
//        }
//      default:
//        return 0;
//      }

//  }

---