SimDaqPlaneTrigger.cxx

Go to the documentation of this file.

#include "SimDaqTrigMaker.h"
#include "SimDaqPlaneTrigger.h"
#include "MessageService/MsgService.h"
#include <vector>

CVSID("$Id: SimDaqPlaneTrigger.cxx,v 1.5 2007/03/01 17:22:11 rhatcher Exp $");

SimDaqTrigMakerProxy<SimDaqPlaneTrigger> gSimDaqPlaneTriggerProxy("PlaneTrigger");

ClassImp(SimDaqPlaneTrigger)


Bool_t SimDaqPlaneTrigger::ApplyTriggerAt(PlexHandle    plex,
                                          SimDigitList& sortedDigitList, // Digit list
                                          UInt_t         startAt,         // Which digit to start looking at
                                          UInt_t         stopAt           // Which digit is the end of the burst
                                          )
{
  MSG("DetSim",Msg::kDebug) << "Plane Trigger Running." << endl;

  // Quick check: do we have enough hits at all?
  if(Int_t(stopAt-startAt) < fN) return false; // Not enough digits total.

  // Add the hits to a plane map.
  std::vector<Int_t> planeMap(600);
  for(UInt_t i = 0;i<planeMap.size();i++) planeMap[i]=0;


  Int_t totplanes = 0;
  Int_t firstPlane = 9999;
  Int_t lastPlane = -9999;

  for(UInt_t idig = startAt; idig<stopAt; idig++) {
    // Find the plane.
    if(sortedDigitList[idig].GetErrors()) continue;
    PlexPixelSpotId psid = sortedDigitList[idig].GetPixelSpotId();
    PlexStripEndId  seid = plex.GetStripEndId(psid);
    if(!seid.IsValid()) {
      psid.SetSpot(1);
      seid = plex.GetStripEndId(psid);
    }
    Int_t plane = seid.GetPlane();

    if(plane > (Int_t)planeMap.size()) continue;

    // Add plane to map.
    if(planeMap[plane] ==0){
      totplanes++;
      if(plane<firstPlane) firstPlane = plane;
      if(plane>lastPlane) lastPlane = plane;
    }

    planeMap[plane]++;
  }
  
  // Second quick check: enough planes were hit in total.
  if(totplanes < fN) return false;

  for(int i = firstPlane; i <= lastPlane; i++) {
    Int_t nPlane = 0;

    for (Int_t j=0; (j<fM) && (i+j<=lastPlane); j++) {
      if ( planeMap[i+j]>0 ) nPlane++;
    }
    if (nPlane>=fN) {
      return true;
    }
    if (i+fM>lastPlane) break; // We're done... stop looking.
  }
  
  return false;
}

---