FiltSGateModule Class Reference

#include <FiltSGateModule.h>

Inheritance diagram for FiltSGateModule:

List of all members.

Public Member Functions
	FiltSGateModule ()
JobCResult	Ana (const MomNavigator *mom)
const Registry &	DefaultConfig () const
void	Config (const Registry &r)
void	Report ()
void	EndJob ()
void	ProcessVtmTimeInfoBlks (std::vector< const RawVtmTimeInfoBlock * > &vtmBlks)
bool	ProcessDigitDataBlks (std::vector< const RawDigitDataBlock * > &digitBlks, Int_t tfnum)
Double_t	near_tdc_sec (Int_t tdc)
void	DrawHistgrams ()
Private Attributes
Int_t	fTicksPerSec
Int_t	fQieErrorMask
Int_t	fPedSubAdcMin
Int_t	fCheckTimeFrame
Double_t	fTminWindow
Double_t	fTmaxWindow
Int_t	fTriggerMaskAllow
Int_t	fDrawHists
bool	fWriteHists
bool	fKeepCanvas
TCanvas *	fCanvas
TH1I	fLastRFClockHist
TH1I	fSpillTimeRFHist
TH1D	fAdcVsTimeRFHist1
Double_t	fHist1TScale
TH1D	fSnarlVsTimeRFHist
TH1D	fAdcVsTimeRFHist2
Double_t	fHist2TScale
TNtuple	fNtuple
Int_t	fLatestRun
Int_t	fLatestSubRun
Int_t	fLatestSpillTimeTF
Int_t	fLatestSpillTimeRF
Int_t	fLastRFClock
Int_t	fLastRFClockPrevTF
Int_t	fLatchedSpillTimeTF

---

Detailed Description

Definition at line 24 of file FiltSGateModule.h.

---

Constructor & Destructor Documentation

FiltSGateModule::FiltSGateModule

(

)

Definition at line 80 of file FiltSGateModule.cxx.

References fAdcVsTimeRFHist1, fAdcVsTimeRFHist2, fLastRFClockHist, fSnarlVsTimeRFHist, and fSpillTimeRFHist.

                                 :
  // Configuration data
  fTicksPerSec(kDefaultTicksPerSec),
  fQieErrorMask(0),
  fPedSubAdcMin(0),
  fCheckTimeFrame(0),
  fTminWindow(-20.0),  // SGate window (usec) based on digitTDC - VtmSpillRF
  fTmaxWindow( 20.0),
  fTriggerMaskAllow(0),
  fDrawHists(1),
  fWriteHists(true),
  fKeepCanvas(true),
  // Classes for drawing
  fCanvas(new TCanvas(std::tmpnam(0),"SGate Filter",800,800)),
  // Histograms
  //           Name of object  Histogram title             bins Low  High
  //           ==============  ==========================  ==== ==== ====
  fLastRFClockHist("fLastRFClock", 
                    Form("Last RF Clock - %d",kDefaultTicksPerSec),
                    150,-25,125),
  fSpillTimeRFHist("fSpillTimeRF", "Spill Time RF",            50,  0,53104000),
  fAdcVsTimeRFHist1("fAdcVsTimeRFHist1", "Adc vs. Time rel RF", 500,-20.0,20.0),
  fSnarlVsTimeRFHist("fSnarlVsTimeRFHist1", "Snarl vs. Time rel RF", 500,-20.0,20.0),
  fAdcVsTimeRFHist2("fAdcVsTimeRFHist2", "Adc vs. Time rel RF", 500,-100.0,100.0),
  // Ntuple
  fNtuple("fNtuple","Event Summary",gsChTags),
  // Module data
  fLatestRun(-1),
  fLatestSubRun(-1),
  fLatestSpillTimeTF(-1),
  fLatestSpillTimeRF(-1),
  fLastRFClock(-1),
  fLastRFClockPrevTF(-1),
  fLatchedSpillTimeTF(-1)

{
//======================================================================
// Set the initial state of your module. The code above initializes the
// variables listed to the values in parentheses. In this case that's
// all I need to do so the method itself is empty...
//======================================================================

  fLastRFClockHist.SetXTitle("Ticks");
  fLastRFClockHist.SetYTitle("VtmTimeInfo blocks");

  fSpillTimeRFHist.SetXTitle("Ticks");
  fSpillTimeRFHist.SetYTitle("VtmTimeInfo blocks");

  fAdcVsTimeRFHist1.SetXTitle("TDC-SpillTimeRF (sec)");
  fAdcVsTimeRFHist1.SetYTitle("#Sigma Adc");

  fSnarlVsTimeRFHist.SetXTitle("TDC-SpillTimeRF (sec)");
  fSnarlVsTimeRFHist.SetYTitle("#Snarls");

  fAdcVsTimeRFHist2.SetXTitle("TDC-SpillTimeRF (sec)");
  fAdcVsTimeRFHist2.SetYTitle("#Sigma Adc");

}

---

Member Function Documentation

JobCResult FiltSGateModule::Ana

(

const MomNavigator *

mom

)

[virtual]

Implement this for read only access to the MomNavigator

Reimplemented from JobCModule.

Definition at line 326 of file FiltSGateModule.cxx.

References DrawHistgrams(), fDrawHists, fLatestRun, fLatestSpillTimeTF, fLatestSubRun, MomNavigator::FragmentIter(), fTriggerMaskAllow, RecMinos::GetHeader(), RawRecord::GetRawBlockIter(), RawDaqHeader::GetRun(), RawDaqSnarlHeader::GetSnarl(), RawDaqHeader::GetSubRun(), RawDaqHeader::GetTimeFrameNum(), RawDaqSnarlHeader::GetTrigSrc(), JobCResult::kAOK, Msg::kDebug, JobCResult::kFailed, JobCResult::kPassed, Msg::kSynopsis, MSG, ProcessDigitDataBlks(), and ProcessVtmTimeInfoBlks().

{
//======================================================================
// Given the object to which the data is attached perform some
// analysis of the data
//======================================================================

  vector<const RawVtmTimeInfoBlock*> vtmBlks;
  vector<const RawDigitDataBlock*> digitBlks;
  Int_t snarlTfNum = -1;

  int    nrawrec=0;
  int    nrawblk=0;
  bool   acceptBecauseTrigMask = false;
  UInt_t trigSrc = 0;
  Int_t  snarlNum = -1;
  bool   sawMonitorRecord = false;

  TObject* recObj = 0;
  TIter recItr = mom->FragmentIter();
  while ( ( recObj = recItr.Next() ) ) {
    RawRecord* rawrec = dynamic_cast<RawRecord*>(recObj);
    if (!rawrec) continue;  // skip things in mom that aren't RawRecords
    nrawrec++;

    // look for the DaqSnarl records with TrigSrc having SGate bit set
    const RawDaqSnarlHeader* snarlHdr = 
      dynamic_cast<const RawDaqSnarlHeader*>(rawrec->GetHeader());
    if ( snarlHdr ) {
      fLatestRun    = snarlHdr->GetRun();
      fLatestSubRun = snarlHdr->GetSubRun();
      snarlTfNum    = snarlHdr->GetTimeFrameNum();
      snarlNum      = snarlHdr->GetSnarl();
      if ( snarlHdr->GetTrigSrc() & fTriggerMaskAllow ) {
        // rawrec->Print();
        acceptBecauseTrigMask = true;
        trigSrc = snarlHdr->GetTrigSrc();
      }
    } // snarl header
    else
      sawMonitorRecord = true;

    TObject* blkObj = 0;
    TIter blkItr = rawrec->GetRawBlockIter();
    while ( ( blkObj = blkItr.Next() ) ) {
      const RawDataBlock* rawBlk = dynamic_cast<const RawDataBlock*>(blkObj);
      if (!rawBlk) continue;  // odd!
      nrawblk++;

      const RawVtmTimeInfoBlock* vtmblk = 
        dynamic_cast<const RawVtmTimeInfoBlock*>(rawBlk);
      if (vtmblk) vtmBlks.push_back(vtmblk);
      
      const RawDigitDataBlock* digitblk = 
        dynamic_cast<const RawDigitDataBlock*>(rawBlk);
      if (digitblk) digitBlks.push_back(digitblk);

    } // loop over blocks
  }  // loop over records

  MSG("Filt",Msg::kDebug)
   << "Saw " << (vtmBlks.size()) << " RawVtmTimeInfoBlocks "
   << "and " << (digitBlks.size()) << " RawDigitDataBlocks "
   << " from " << nrawblk << " blocks in " << nrawrec << " records"
   << endl;

  // first process all available VtmTimeInfo blocks
  ProcessVtmTimeInfoBlks(vtmBlks);

  // then process digits
  bool digitsInWindow = ProcessDigitDataBlks(digitBlks,snarlTfNum);

  // If graphics are enabled (fDrawHist!=0), draw the histograms if:
  //    * there were digits seen in the time window
  //    * fDrawHist is < 0 (e.g. ever snarl record --> very slow!)
  //    * fDrawHist is > 0 and that many TimeFrames have passed 
  //        since the histograms were last drawn
  if (fDrawHists!=0) {
    static Int_t lastTF_drawn = -9999; // when we see N-th new time frame ...
    if ( ( digitsInWindow ) ||
         ( fDrawHists < 0 ) || 
         ( fLatestSpillTimeTF >= lastTF_drawn+fDrawHists  ) )
      {
        this->DrawHistgrams();
        lastTF_drawn = fLatestSpillTimeTF;
      }
  }

  if ( acceptBecauseTrigMask || digitsInWindow ) {
    MSG("Filt",Msg::kSynopsis )
      << std::right
      << "FiltSGate Snarl=" << setw(6) << snarlNum 
      << " in TF=" << setw(5) << snarlTfNum
      << " TrigMask="
      << (acceptBecauseTrigMask?"yes":"no ")
      << " (0x" << std::hex <<  std::setfill('0') << std::right << setw(8)
      << trigSrc
      << std::setfill(' ') << std::left << std::dec
      << ") digitsInWindow="
      << (digitsInWindow?"yes":"no ")
      << endl;
    return JobCResult::kPassed; // Passed cut
  }
  if ( sawMonitorRecord ) {
    MSG("Filt",Msg::kDebug )
      << "FiltSGateModule AOK because DaqMonitor record in recset."
      << endl;
    return JobCResult::kAOK;   // Passed
  }

  return JobCResult::kFailed; // no reason to keep any of these records
}

void FiltSGateModule::Config

(

const Registry &

r

)

[virtual]

Return the actual configuration. If your module directly pulls its configuration from the fConfig Registry, you don't need to override this. Override if you have local config variables.

Reimplemented from JobCModule.

Definition at line 204 of file FiltSGateModule.cxx.

References chooseTimeScale(), fAdcVsTimeRFHist1, fAdcVsTimeRFHist2, fCheckTimeFrame, fDrawHists, fHist1TScale, fHist2TScale, fKeepCanvas, fLastRFClockHist, Form(), fPedSubAdcMin, fQieErrorMask, fSnarlVsTimeRFHist, fTicksPerSec, fTmaxWindow, fTminWindow, fTriggerMaskAllow, fWriteHists, Registry::Get(), Msg::kSynopsis, and MSG.

{
//======================================================================
// Configure the module given the registry r
//======================================================================
//  char   tmpb;
  double tmpd;
  int    tmpi;

  if (r.Get("TicksPerSec",     tmpi)) { fTicksPerSec      = tmpi; }
  if (r.Get("QieErrorMask",    tmpi)) { fQieErrorMask     = tmpi; }
  if (r.Get("PedSubAdcMin",    tmpi)) { fPedSubAdcMin     = tmpi; }
  if (r.Get("CheckTimeFrame",  tmpi)) { fCheckTimeFrame   = tmpi; }

  if (r.Get("TminWindow",      tmpd)) { fTminWindow       = tmpd; }
  if (r.Get("TmaxWindow",      tmpd)) { fTmaxWindow       = tmpd; }

  if (r.Get("TriggerMaskAllow",tmpi)) { fTriggerMaskAllow = tmpi; }

  if (r.Get("DrawHists",       tmpi)) { fDrawHists        = tmpi; }
  if (r.Get("WriteHists",      tmpi)) { fWriteHists       = tmpi; }
  if (r.Get("KeepCanvas",      tmpi)) { fKeepCanvas       = tmpi; }

  // histogram parameters
  bool   reBin;
  int    nbin;
  double xmin, xmax;
  const char* timeUnitsTag;

  reBin = false;
  nbin = fAdcVsTimeRFHist1.GetNbinsX();
  xmin = fAdcVsTimeRFHist1.GetXaxis()->GetXmin();
  xmax = fAdcVsTimeRFHist1.GetXaxis()->GetXmax();
  if (r.Get("Nbin1",tmpi)) { reBin = true; nbin = tmpi; } 
  if (r.Get("Tmin1",tmpd)) { reBin = true; xmin = tmpd; } 
  if (r.Get("Tmax1",tmpd)) { reBin = true; xmax = tmpd; } 
  if (reBin) {
    MSG("Filt",Msg::kSynopsis)
      << "Rebin Hist1 nbins=" << nbin 
      << " [ " << xmin << " : " << xmax << " ]" 
      << endl;
    chooseTimeScale(xmax-xmin,fHist1TScale,timeUnitsTag);
    fAdcVsTimeRFHist1.SetBins(nbin,xmin*fHist1TScale,xmax*fHist1TScale);
    fAdcVsTimeRFHist1.SetXTitle(Form("TDC-SpillTimeRF (%s)",timeUnitsTag));
    fSnarlVsTimeRFHist.SetBins(nbin,xmin*fHist1TScale,xmax*fHist1TScale);
    fSnarlVsTimeRFHist.SetXTitle(Form("TDC-SpillTimeRF (%s)",timeUnitsTag));
  }

  reBin = false;
  nbin = fAdcVsTimeRFHist2.GetNbinsX();
  xmin = fAdcVsTimeRFHist2.GetXaxis()->GetXmin();
  xmax = fAdcVsTimeRFHist2.GetXaxis()->GetXmax();
  if (r.Get("Nbin2",tmpi)) { reBin = true; nbin = tmpi; } 
  if (r.Get("Tmin2",tmpd)) { reBin = true; xmin = tmpd; } 
  if (r.Get("Tmax2",tmpd)) { reBin = true; xmax = tmpd; } 
  if (reBin) {
    MSG("Filt",Msg::kSynopsis)
      << "Rebin Hist2 nbins=" << nbin 
      << " [ " << xmin << " : " << xmax << " ]" 
      << endl;
    chooseTimeScale(xmax-xmin,fHist2TScale,timeUnitsTag);
    fAdcVsTimeRFHist2.SetBins(nbin,xmin*fHist2TScale,xmax*fHist2TScale);
    fAdcVsTimeRFHist2.SetXTitle(Form("TDC-SpillTimeRF (%s)",timeUnitsTag));
  }

  reBin = false;
  nbin = fLastRFClockHist.GetNbinsX();
  xmin = fLastRFClockHist.GetXaxis()->GetXmin();
  xmax = fLastRFClockHist.GetXaxis()->GetXmax();
  if (r.Get("NbinLastRF",tmpi)) { reBin = true; nbin = tmpi; } 
  if (r.Get("TminLastRF",tmpd)) { reBin = true; xmin = tmpd; } 
  if (r.Get("TmaxLastRF",tmpd)) { reBin = true; xmax = tmpd; } 
  if (reBin) {
    MSG("Filt",Msg::kSynopsis)
      << "Rebin RFClock nbins=" << nbin 
      << " [ " << xmin << " : " << xmax << " ]" 
      << endl;
    fLastRFClockHist.SetBins(nbin,xmin,xmax);
  }
}

const Registry & FiltSGateModule::DefaultConfig

(

)

const [virtual]

Get the default configuration registry. This should normally be overridden. One useful idiom is to implement it like:

const Registry& MyModule::DefaultConfig() const { static Registry cfg; // never is destroyed if (cfg.Size()) return cfg; // already filled it // set defaults: cfg.Set("TheAnswer",42); cfg.Set("Units","unknown"); return cfg; }

Reimplemented from JobCModule.

Definition at line 141 of file FiltSGateModule.cxx.

References MuELoss::e, JobCModule::GetName(), RawQieDigit::kCapIdError, RawQieDigit::kCapIdSameMenuError, RawQieDigit::kCapIdSharedMinderError, kDefaultTicksPerSec, RawQieDigit::kMisCountError, RawQieDigit::kParityError, RawQieDigit::kPriorityTruncate, Registry::LockValues(), Registry::Set(), TRIGGER_BIT_SPILL, TRIGGER_BIT_SPILL_IP, and Registry::UnLockValues().

{
//======================================================================
// Create a registry which holds the default configuration and return it
//======================================================================
  static Registry r;

  // Set name of config
  std::string name = this->GetName();
  name += ".config.default";
  r.SetName(name.c_str());

  Int_t qieErrorMask = 0;
  qieErrorMask |= RawQieDigit::kParityError;             // 0x0001
  qieErrorMask |= RawQieDigit::kCapIdError;              // 0x0002
  qieErrorMask |= RawQieDigit::kMisCountError;           // 0x0004
  qieErrorMask |= RawQieDigit::kPriorityTruncate;        // 0x0008
  qieErrorMask |= RawQieDigit::kCapIdSameMenuError;      // 0x0010
  qieErrorMask |= RawQieDigit::kCapIdSharedMinderError;  // 0x0020

  Int_t trigMaskAllow = 0;
  //trigMaskAllow |= TRIGGER_BIT_NULL;       //     0x1
  //trigMaskAllow |= TRIGGER_BIT_RAWLI;      //     0x2
  //trigMaskAllow |= TRIGGER_BIT_PLANE;      //     0x4
  //trigMaskAllow |= TRIGGER_BIT_TOF;        //     0x8
  //trigMaskAllow |= TRIGGER_BIT_ACTIVITY;   //    0x10
  //trigMaskAllow |= TRIGGER_BIT_E4;         //    0x20
  //trigMaskAllow |= TRIGGER_BIT_MINBIAS;    //  0x8000
  trigMaskAllow |= TRIGGER_BIT_SPILL;      // 0x10000
  trigMaskAllow |= TRIGGER_BIT_SPILL_IP;   // 0x20000

  // Set values of config
  r.UnLockValues();
  r.Set("TicksPerSec",      kDefaultTicksPerSec);
  r.Set("QieErrorMask",     qieErrorMask);
  r.Set("PedSubAdcMin",     0);
  r.Set("CheckTimeFrame",   0);
  // Tmin & Tmax values are in seconds!!
  r.Set("TminWindow",      -100.0e-6);
  r.Set("TmaxWindow",       100.0e-6);
  r.Set("TriggerMaskAllow", trigMaskAllow);

  r.Set("DrawHists",        1);
  r.Set("WriteHists",       true);

  r.Set("Nbin1",            200);
  r.Set("Tmin1",           -20.0e-6);
  r.Set("Tmax1",            20.0e-6);
  r.Set("Nbin2",            200);
  r.Set("Tmin2",           -100.0e-6);
  r.Set("Tmax2",            100.0e-6);

  r.Set("NbinLastRF",       600);
  r.Set("TminLastRF",      -25);
  r.Set("TmaxLastRF",       575);

  r.LockValues();

  return r;
}

void FiltSGateModule::DrawHistgrams

(

)

Definition at line 661 of file FiltSGateModule.cxx.

References fAdcVsTimeRFHist1, fCanvas, fLastRFClockHist, fSnarlVsTimeRFHist, and fSpillTimeRFHist.

Referenced by Ana(), and Report().

{
//======================================================================
// Draw the histograms
//======================================================================
  
  static bool first = true;
  if (first) {
    first = false;
    fCanvas->Divide(1,4); // Like PAW's zone command
    fCanvas->cd(4);
    gPad->Divide(2,1);
    gStyle->SetOptStat(111111);
  }

  fCanvas->cd(1); 
  gPad->SetLogy(0); 
  fAdcVsTimeRFHist1.DrawCopy();

  fCanvas->cd(2); 
  gPad->SetLogy(0); 
  fSnarlVsTimeRFHist.DrawCopy();

  fCanvas->cd(3); 
  gPad->SetLogy(0); 
  fAdcVsTimeRFHist1.DrawCopy();

  fCanvas->cd(4); 
  gPad->cd(1);
  //gPad->SetLogy(1); 
  fLastRFClockHist.DrawCopy();

  fCanvas->cd(4); 
  gPad->cd(2);
  //gPad->SetLogy(1); 
  fSpillTimeRFHist.DrawCopy();

  fCanvas->cd(1);
  fCanvas->Update();
}

void FiltSGateModule::EndJob

(

)

[virtual]

Implement for notification of end of job

Reimplemented from JobCModule.

Definition at line 302 of file FiltSGateModule.cxx.

References fAdcVsTimeRFHist1, fAdcVsTimeRFHist2, fCanvas, fKeepCanvas, fLastRFClockHist, fNtuple, fSnarlVsTimeRFHist, fSpillTimeRFHist, and fWriteHists.

{
//======================================================================
// At the end of the job print some stuff out and save the histogram
// to a file
//======================================================================
  
  // Create a file to write the histogram to. By default it becomes
  // the current directory. Write uses current directory by default
  if (fWriteHists) {
    TFile userFile("sgate_filter.root","RECREATE");
    fLastRFClockHist.Write(); 
    fSpillTimeRFHist.Write();
    fAdcVsTimeRFHist1.Write();
    fSnarlVsTimeRFHist.Write();
    fAdcVsTimeRFHist2.Write();
    fNtuple.Write();
  }

  if ( ! fKeepCanvas ) delete fCanvas;
}

Double_t FiltSGateModule::near_tdc_sec

(

Int_t

tdc

)

Definition at line 651 of file FiltSGateModule.cxx.

References fTicksPerSec.

Referenced by ProcessDigitDataBlks().

{
  return tdc/((double)fTicksPerSec);

  // Peter's version
  // return tdc*1.5625*16.0*58.0/77.0;
}

bool FiltSGateModule::ProcessDigitDataBlks	(	std::vector< const RawDigitDataBlock * > &	digitBlks,
		Int_t	tfnum
	)

Definition at line 523 of file FiltSGateModule.cxx.

References RawQieDigit::AsString(), RawChannelId::ClearModeBits(), digit(), fAdcVsTimeRFHist1, fAdcVsTimeRFHist2, fCheckTimeFrame, fHist1TScale, fHist2TScale, fLastRFClockPrevTF, fLatchedSpillTimeTF, fLatestRun, fLatestSpillTimeRF, fLatestSubRun, fPedSubAdcMin, fQieErrorMask, fSnarlVsTimeRFHist, fTmaxWindow, fTminWindow, Msg::kDebug, Detector::kNear, kPedAdc, ElecType::kQIE, RawQieDigit::kSGate, Msg::kSynopsis, Msg::kWarning, MSG, and near_tdc_sec().

Referenced by Ana().

{
//======================================================================
// Process the DigitDataBlocks
//   Determine whether any digit TDC's were in the user specified
//   time window.  Histogram TDC distribution.
//
// Inputs:  cdlh - Pointer to a candidate digit list handle
//
// Outputs: true is any TDC's were in the time window
//======================================================================

  size_t n = digitBlks.size();
  if ( ! n ) return false;   // nothing to process

  bool sawWindowDigits = false;

  // real bad hack !!!!!
#ifdef RAWCHANNELID_NEW_CTOR
  RawChannelId rcid_bad(Detector::kNear,ElecType::kQIE,5,18,0,12);
#else
  RawChannelId rcid_bad(Detector::kNear,ElecType::kQIE,5,18,0,12,false,false);
#endif

  for (size_t i = 0; i < n; ++i) {
    const RawDigitDataBlock* rddb = digitBlks[i];
    
    TIter rditr = rddb->GetDatumIter();
    const RawDigit* digit = 0;
    Int_t nd = 0;
    //Long64_t sumtdc = 0;
    Int_t mintdc = 0x7FFFFFF;
    Int_t maxtdc = -(mintdc-1);
    Char_t dtlast = -1;
    Msg::LogLevel_t dumpLevel = Msg::kDebug;

    Double_t sumAdcTime = 0;
    Double_t sumAdc     = 0;

    // What we're doing here for calculating times isn't really
    // right for events that span the TF boundary ... we should
    // be making use of the crate T0 and not the TF #.
    Double_t tfOffset =
      (snarlTfNum-fLatchedSpillTimeTF) * near_tdc_sec(fLastRFClockPrevTF);
    if ( fCheckTimeFrame == 0 ) tfOffset = 0;

    while ( ( digit = dynamic_cast<RawDigit*>(rditr.Next()) ) ) {

      // skip digits that have error codes on that we don't want
      if ( digit->GetErrorCode() & fQieErrorMask ) continue;

      // real bad hack !!!!!
      RawChannelId rcid = digit->GetChannel();
      rcid.ClearModeBits();
      if ( rcid == rcid_bad ) continue;

      Int_t tdc = digit->GetTDC();
      Double_t subTime_sec = 
        near_tdc_sec(tdc-fLatestSpillTimeRF) + tfOffset;

      // don't forget to subtract off the base ADC value of 50
      Int_t pedSubAdc = digit->GetADC()-kPedAdc;
      // skip digits that are less than the minimum ADC
      if ( pedSubAdc < fPedSubAdcMin ) continue;

      fAdcVsTimeRFHist1.Fill(subTime_sec*fHist1TScale,pedSubAdc);
      fAdcVsTimeRFHist2.Fill(subTime_sec*fHist2TScale,pedSubAdc);

      sumAdcTime += pedSubAdc * subTime_sec*fHist1TScale;
      sumAdc     += pedSubAdc;
      
      if ( subTime_sec >= fTminWindow && 
           subTime_sec <= fTmaxWindow    ) {
        sawWindowDigits = true;
        dumpLevel = Msg::kSynopsis;
      }

      nd++;
      //sumtdc += tdc;
      if ( tdc < mintdc ) mintdc = tdc;
      if ( tdc > maxtdc ) maxtdc = tdc;
      Char_t dt = digit->GetDataType();
      if ( dtlast != dt && dtlast != -1 ) {
        MSG("Filt",Msg::kWarning)
          << " RawDigitDataBlock has mixed DataTypes!" 
          << (Int_t)dtlast << " -> " << (Int_t)dt 
          << endl;
      }
      dtlast = dt;
    }  // loop over digits

    if ( sumAdc > 0.0 )  fSnarlVsTimeRFHist.Fill(sumAdcTime/sumAdc,1.0);

    //Int_t avgtdc = sumtdc/nd;
    const char* tagrf = " ";
    if ( mintdc <= fLatestSpillTimeRF && 
         fLatestSpillTimeRF <= maxtdc     ) tagrf = "<== ";

    MSG("Filt",dumpLevel)
      << std::right
      << " Run " << setw(5) << fLatestRun
      << "/" << setw(2) <<  fLatestSubRun
      << " TF " << setw(5) << snarlTfNum
      << " TDC"
      << " min " << setw(9) << mintdc 
      << " max " << setw(9) << maxtdc 
      //<< " avg " << setw(9) << avgtdc 
      //<< " usec " << setw(10) << setprecision(7)
      //<< near_tdc_sec(avgtdc)/Munits::microsecond
      << " nd " << setw(5) << nd << " "
      << RawQieDigit::AsString((RawQieDigit::EQieDataTypes)dtlast)
      << (dtlast==RawQieDigit::kSGate?"  ##":"   ")
      << tagrf
      //<< endl
      //<< "   " << mintdc 
      //<< ((mintdc<=fLatestSpillTimeRF)?" < ":" !< ") 
      //<< fLatestSpillTimeRF 
      //<< ((fLatestSpillTimeRF<=maxtdc)?" < ":" !< ") 
      //<< maxtdc 
      << endl;

  }  // loop over blocks

  return sawWindowDigits;
}

void FiltSGateModule::ProcessVtmTimeInfoBlks

(

std::vector< const RawVtmTimeInfoBlock * > &

vtmBlks

)

Definition at line 441 of file FiltSGateModule.cxx.

References fLastRFClock, fLastRFClockHist, fLastRFClockPrevTF, fLatchedSpillTimeTF, fLatestRun, fLatestSpillTimeRF, fLatestSpillTimeTF, fLatestSubRun, fSpillTimeRFHist, Msg::kDebug, kDefaultTicksPerSec, Msg::kSynopsis, and MSG.

Referenced by Ana().

{
//======================================================================
// Calculate the total charge and number of hit strips in an event
//
// Inputs:  cdlh - Pointer to a candidate digit list handle
//
// Outputs: nHit - total number of digits in the list
//======================================================================

  size_t n = vtmBlks.size();
  if ( ! n ) return;   // nothing to process

  const size_t maxCrates = 16;
  int crate[maxCrates];
  int tfnum[maxCrates];
  int spillTimeRF[maxCrates];
  int lastRFClock[maxCrates];
  Msg::LogLevel_t dumpLevel = Msg::kDebug;

  Int_t  sum = 0;

  for (size_t i = 0; i < n; ++i) {
    const RawVtmTimeInfoBlock* vti = vtmBlks[i];

    crate[i]       = vti->GetCrate();
    tfnum[i]       = vti->GetTimeFrameNum();
    spillTimeRF[i] = vti->GetSpillTimeRF();
    lastRFClock[i] = vti->GetLastTimeStamp();
    // Run & SubRun should be the same in all the blocks
    if ( 0 == i ) {
      fLatestRun     = vti->GetRun();
      fLatestSubRun  = vti->GetSubRun();
    }

    sum += spillTimeRF[i];

    fLastRFClockHist.Fill(lastRFClock[i]-kDefaultTicksPerSec);

    MSG("Filt",Msg::kDebug)
      << " TF " << setw(5) << tfnum[i] 
      << " VtmTimeInfo crate " << setw(2) << crate[i] 
      << " RF " << setw(9) << spillTimeRF[i] 
      << " GPS " << setw(9) << vti->GetSpillTimeGPS()
      << endl;
  }

  Int_t fSpillTimeRFAvg = sum/n;
  fLatestSpillTimeTF      = (Int_t)TMath::Median(n,tfnum);
  Int_t medianSpillTimeRF = (Int_t)TMath::Median(n,spillTimeRF);
  // We recognize a new latching of a spill time by it having
  // a different value -- and thus we know in which TimeFrame
  // the latch occurred.  Or so we assume -- and yes it is broken
  // if the new spill time happens *exactly* at the same time.
  if ( medianSpillTimeRF != fLatestSpillTimeRF ) {
    fLatchedSpillTimeTF   = fLatestSpillTimeTF;
    // fill only for new latched value (and not initial junk value)
    if (fLatestSpillTimeRF>=0) fSpillTimeRFHist.Fill(fLatestSpillTimeRF); 
    dumpLevel = Msg::kSynopsis;
  }
  fLatestSpillTimeRF      = medianSpillTimeRF;
  fLastRFClockPrevTF      = fLastRFClock;
  fLastRFClock            = (Int_t)TMath::Median(n,lastRFClock);

  MSG("Filt",dumpLevel) 
    << std::right
    << " Run " << setw(5) << fLatestRun
    << "/" << setw(2) <<  fLatestSubRun
    << " TF " << setw(5) << fLatestSpillTimeTF
    << " (" << setw(5) << fLatchedSpillTimeTF << ")"
    << " SpillTimeRF"
    << " median " << setw(9) << fLatestSpillTimeRF
    << ( (fLatestSpillTimeRF == fSpillTimeRFAvg) ? " consistent":" median!=avg")
    //<< " avg " << setw(9) << fSpillTimeRFAvg 
    //<< " usec " << setw(10) << setprecision(7)
    //<< near_tdc_sec(fSpillTimeRFAvg)/Munits::microsecond
    //<< " *******"
    << endl;
}

void FiltSGateModule::Report

(

)

[virtual]

Implement to spew end of running report

Reimplemented from JobCModule.

Definition at line 287 of file FiltSGateModule.cxx.

References DrawHistgrams(), fDrawHists, Msg::kInfo, and MSG.

{
//======================================================================
// ** No longer needed. Use Module::Report() to see list of parameters
// ** which can be set using Set("name=value").
//
// Print help for this module
//======================================================================
  MSG("Demo",Msg::kInfo) 
    << "Report For FiltSGateModule Module:\n" << endl;
  if (fDrawHists!=0) this->DrawHistgrams();
}

---

Member Data Documentation

TH1D FiltSGateModule::fAdcVsTimeRFHist1 [private]

Definition at line 59 of file FiltSGateModule.h.

Referenced by Config(), DrawHistgrams(), EndJob(), FiltSGateModule(), and ProcessDigitDataBlks().

TH1D FiltSGateModule::fAdcVsTimeRFHist2 [private]

Definition at line 62 of file FiltSGateModule.h.

Referenced by Config(), EndJob(), FiltSGateModule(), and ProcessDigitDataBlks().

TCanvas* FiltSGateModule::fCanvas [private]

Definition at line 56 of file FiltSGateModule.h.

Referenced by DrawHistgrams(), and EndJob().

Int_t FiltSGateModule::fCheckTimeFrame [private]

Definition at line 47 of file FiltSGateModule.h.

Referenced by Config(), and ProcessDigitDataBlks().

Int_t FiltSGateModule::fDrawHists [private]

Definition at line 52 of file FiltSGateModule.h.

Referenced by Ana(), Config(), and Report().

Double_t FiltSGateModule::fHist1TScale [private]

Definition at line 60 of file FiltSGateModule.h.

Referenced by Config(), and ProcessDigitDataBlks().

Double_t FiltSGateModule::fHist2TScale [private]

Definition at line 63 of file FiltSGateModule.h.

Referenced by Config(), and ProcessDigitDataBlks().

bool FiltSGateModule::fKeepCanvas [private]

Definition at line 54 of file FiltSGateModule.h.

Referenced by Config(), and EndJob().

Int_t FiltSGateModule::fLastRFClock [private]

Definition at line 70 of file FiltSGateModule.h.

Referenced by ProcessVtmTimeInfoBlks().

TH1I FiltSGateModule::fLastRFClockHist [private]

Definition at line 57 of file FiltSGateModule.h.

Referenced by Config(), DrawHistgrams(), EndJob(), FiltSGateModule(), and ProcessVtmTimeInfoBlks().

Int_t FiltSGateModule::fLastRFClockPrevTF [private]

Definition at line 71 of file FiltSGateModule.h.

Referenced by ProcessDigitDataBlks(), and ProcessVtmTimeInfoBlks().

Int_t FiltSGateModule::fLatchedSpillTimeTF [private]

Definition at line 72 of file FiltSGateModule.h.

Referenced by ProcessDigitDataBlks(), and ProcessVtmTimeInfoBlks().

Int_t FiltSGateModule::fLatestRun [private]

Definition at line 66 of file FiltSGateModule.h.

Referenced by Ana(), ProcessDigitDataBlks(), and ProcessVtmTimeInfoBlks().

Int_t FiltSGateModule::fLatestSpillTimeRF [private]

Definition at line 69 of file FiltSGateModule.h.

Referenced by ProcessDigitDataBlks(), and ProcessVtmTimeInfoBlks().

Int_t FiltSGateModule::fLatestSpillTimeTF [private]

Definition at line 68 of file FiltSGateModule.h.

Referenced by Ana(), and ProcessVtmTimeInfoBlks().

Int_t FiltSGateModule::fLatestSubRun [private]

Definition at line 67 of file FiltSGateModule.h.

Referenced by Ana(), ProcessDigitDataBlks(), and ProcessVtmTimeInfoBlks().

TNtuple FiltSGateModule::fNtuple [private]

Definition at line 64 of file FiltSGateModule.h.

Referenced by EndJob().

Int_t FiltSGateModule::fPedSubAdcMin [private]

Definition at line 46 of file FiltSGateModule.h.

Referenced by Config(), and ProcessDigitDataBlks().

Int_t FiltSGateModule::fQieErrorMask [private]

Definition at line 45 of file FiltSGateModule.h.

Referenced by Config(), and ProcessDigitDataBlks().

TH1D FiltSGateModule::fSnarlVsTimeRFHist [private]

Definition at line 61 of file FiltSGateModule.h.

Referenced by Config(), DrawHistgrams(), EndJob(), FiltSGateModule(), and ProcessDigitDataBlks().

TH1I FiltSGateModule::fSpillTimeRFHist [private]

Definition at line 58 of file FiltSGateModule.h.

Referenced by DrawHistgrams(), EndJob(), FiltSGateModule(), and ProcessVtmTimeInfoBlks().

Int_t FiltSGateModule::fTicksPerSec [private]

Definition at line 44 of file FiltSGateModule.h.

Referenced by Config(), and near_tdc_sec().

Double_t FiltSGateModule::fTmaxWindow [private]

Definition at line 49 of file FiltSGateModule.h.

Referenced by Config(), and ProcessDigitDataBlks().

Double_t FiltSGateModule::fTminWindow [private]

Definition at line 48 of file FiltSGateModule.h.

Referenced by Config(), and ProcessDigitDataBlks().

Int_t FiltSGateModule::fTriggerMaskAllow [private]

Definition at line 50 of file FiltSGateModule.h.

Referenced by Ana(), and Config().

bool FiltSGateModule::fWriteHists [private]

Definition at line 53 of file FiltSGateModule.h.

Referenced by Config(), and EndJob().

---

The documentation for this class was generated from the following files:

• FiltSGateModule.h
• FiltSGateModule.cxx

---