SimDaqEnergyTrigger.cxx

Go to the documentation of this file.

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

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

SimDaqTrigMakerProxy<SimDaqEnergyTrigger> gSimDaqEnergyTriggerProxy("EnergyTrigger");

ClassImp(SimDaqEnergyTrigger)


Bool_t SimDaqEnergyTrigger::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
                                          )
{
  //
  // This trigger looks if the energy sum of any M consecutive planes 
  // is larger than Threshold.
  // It has the further requirement that N planes out of M must be hit,
  // and Nhits digits must have occoured.
  //
  MSG("DetSim",Msg::kDebug) << "Energy Trigger Running: "                          
                            << fThreshold << " ADC and " 
                            << fNhits << " hits and " 
                            << fN << " planes hit in a window of " 
                            << fM  << " planes.\n";

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

  // Add the hits to a plane map.
  std::vector<Int_t> planeEnergy(600,0);
  std::vector<Int_t> planeHits(600,0);

  Int_t totADC = 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();
    Int_t adc = sortedDigitList[idig].GetADC();
    totADC+=adc;

    if( (UInt_t)plane >= planeEnergy.size()) continue;

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

    planeEnergy[plane]+= adc;
    planeHits[plane]+= 1;
  }
  
  // Second quick check: enough ADCs in total to fire trigger
  if(totADC < fThreshold) return false;



  for(Int_t i = firstPlane; i <= lastPlane; i++) {
    Int_t local_energy = 0;
    Int_t local_nhit = 0;
    Int_t local_planeshit = 0;

    for (Int_t j=0; (j<fM) && (i+j<=lastPlane); j++) {
      local_energy += planeEnergy[i+j];
      local_nhit   += planeHits[i+j];
      if(planeHits[i+j]>0) local_planeshit++;
    }
    if ((local_energy >= fThreshold) && (local_nhit >=fNhits) && (local_planeshit >= fN)) {
      return true;
    }

    if (i+fM>lastPlane) break; // We're done... stop looking.
  }
  
  return false;
}

---