SimDaqMultTrigger.cxx

Go to the documentation of this file.

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

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

SimDaqTrigMakerProxy<SimDaqMultTrigger> gSimDaqMultTriggerProxy("MultTrigger");

ClassImp(SimDaqMultTrigger)


SimDaqMultTrigger::SimDaqMultTrigger(void):
  SimDaqTrigger()
{
}


Bool_t SimDaqMultTrigger::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 is the standard  PlaneTrigger just working on significant hits
  // which are above some ADC threshold. it furthermore requires
  // a total number of hits to be present in the M plane window
  //
  MSG("DetSim",Msg::kDebug) << "Multiplicity Trigger Running." 
                            << "    Calculating " << fNhits << " hits in " 
                            << fN << " out 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.

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

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

    if(adc>=fThreshold) {
      // 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;
    Int_t nDigit = 0;

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

---