DemoNPlaneCut.cxx

Go to the documentation of this file.

// $Id: DemoNPlaneCut.cxx,v 1.1 2004/01/28 18:42:01 messier Exp $
//
// Select events with at least n continuous planes hit
//
// messier@huhepl.harvard.edu
#include "Demo/DemoNPlaneCut.h"
#include <cstdio>
// ROOT includes
#include "TFile.h"
#include "CandData/CandRecord.h"           // DigitList handling
#include "CandDigit/CandDigitListHandle.h" // "                "
#include "CandDigit/CandDigitHandle.h"     // "                "
#include "DataUtil/GetCandidate.h"         // Unpack data from Mom
#include "JobControl/JobCModuleRegistry.h" // JOBMODULE registration macro

JOBMODULE(DemoNPlaneCut,"NPlaneCut","Cut short events");

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

DemoNPlaneCut::DemoNPlaneCut() :
  fNplaneCut(15),
  fMaxGap(2)
{ }

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

int DemoNPlaneCut::CountPlanes(const CandDigitListHandle* cdlh)
{
  // The loop is done with a class that iterates over the list
  int iPlane = 0;
  memset(fPlane, 0, PLANEARRAYSIZE*sizeof(fPlane[0]));
  CandDigitHandleItr cdhItr(cdlh->GetDaughterIterator());
  for (; cdhItr.IsValid(); cdhItr.Next()) {
    const PlexSEIdAltL&   altList = (*cdhItr)->GetPlexSEIdAltL();
    const PlexStripEndId& strEnd  = altList.GetBestSEId();
    iPlane = strEnd.GetPlane();
    if (iPlane>=0 && iPlane<PLANEARRAYSIZE) fPlane[iPlane] = 1;
  }
  int nPlane    = 0;
  int gapSz     = 0;
  int nPlaneMax = 0;
  for (int i=0; i<PLANEARRAYSIZE; ++i) {
    if (fPlane[i]==1) {
      // We have a hit on a plane. Increase the plane count and zero
      // out the gap counter
      ++nPlane;  
      gapSz = 0;
      if (nPlane > nPlaneMax) nPlaneMax = nPlane;
    }
    else {

      // No hit on the plane. Increase the gap size. If we reach the
      // maximum gap size, start the plane count again

      ++gapSz;
      if (gapSz > fMaxGap) nPlane = 0;
    }
  }
  return nPlaneMax;
}

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

JobCResult DemoNPlaneCut::Ana(const MomNavigator *mom) 
{
//======================================================================
// Given the object to which the data is attached perform some
// analysis of the data
//======================================================================
  // Access to the event is given through handles. In this case I want
  // the list of candidate digits
  const CandDigitListHandle* canddigit = 
    DataUtil::GetCandidate<CandDigitListHandle>(mom,
                                                "CandDigitListHandle", 
                                                "canddigitlist");
  
  // If the GetDigitList is successful, loop over the elements in the
  // list to add up the charge
  int nplane = 0;
  if (canddigit) nplane = this->CountPlanes(canddigit);
  if (nplane >= fNplaneCut) return JobCResult::kPassed;
  return JobCResult::kFailed;
}

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

const Registry& DemoNPlaneCut::DefaultConfig() const
{
//======================================================================
// 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());

  // Set values of config
  r.UnLockValues();
  r.Set("NplaneCut", 15);
  r.Set("MaxGap",    2);
  r.LockValues();

  return r;
}

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

void DemoNPlaneCut::Config(const Registry& r)
{
//======================================================================
// Configure the module given the registry r
//======================================================================
  int    tmpi;
  if (r.Get("NplaneCut", tmpi)) { fNplaneCut = tmpi; }
  if (r.Get("MaxGap",    tmpi)) { fMaxGap    = tmpi; }
}

---