Skimmers/Trigger/TriggerDSTSkimmer.cc

0001
0002 #include "TriggerDSTSkimmer.h"
0003
0004 #include <fun4all/Fun4AllReturnCodes.h>
0005
0006 #include <ffarawobjects/Gl1Packet.h>
0007
0008 #include <phool/getClass.h>
0009
0010 #include <iostream>
0011 #include <map>
0012 #include <utility>
0013
0014 //____________________________________________________________________________..
0015 TriggerDSTSkimmer::TriggerDSTSkimmer(const std::string &name)
0016   : SubsysReco(name)
0017 {
0018 }
0019
0020 //____________________________________________________________________________..
0021 int TriggerDSTSkimmer::process_event(PHCompositeNode *topNode)
0022 {
0023   if (Verbosity() > 0)
0024   {
0025     if (ievent % 1000 == 0)
0026     {
0027       std::cout << "Processing event " << ievent << std::endl;
0028     }
0029     ievent++;
0030   }
0031
0032   // trigger cut
0033   int n_trigger_index = m_trigger_index.size();
0034   for (int it = 0; it < n_trigger_index; ++it)
0035   {
0036     if (m_trigger_index[it] < 0 || m_trigger_index[it] > 63)
0037     {
0038       std::cerr << "ERROR in TriggerDSTSkimmer: trigger index "
0039                 << m_trigger_index[it] << " is out of range [0..63]. "
0040                 << "Aborting run." << std::endl;
0041       return Fun4AllReturnCodes::ABORTRUN;
0042     }
0043   }
0044
0045   if (n_trigger_index != 0)
0046   {
0047     bool trigger_fired = false;
0048     Gl1Packet *_gl1PacketInfo = findNode::getClass<Gl1Packet>(topNode, "GL1Packet");
0049     int gl1_trigger_vector_scaled[64] = {0};
0050     if (_gl1PacketInfo)
0051     {
0052       uint64_t scaledtriggervec = _gl1PacketInfo->lValue(0, "ScaledVector");
0053       for (auto &val : gl1_trigger_vector_scaled)
0054       {
0055         val = ((scaledtriggervec & 0x1U) == 0x1U);
0056         scaledtriggervec = (scaledtriggervec >> 1U) & 0xffffffffU;
0057       }
0058     }
0059     else
0060     {
0061       std::cout << "TriggerDSTSkimmer::process_event - Error - Can't find Trigger Node Gl1Packet therefore no selection can be made" << std::endl;
0062       return Fun4AllReturnCodes::ABORTEVENT;
0063     }
0064     for (int it = 0; it < n_trigger_index; ++it)
0065     {
0066       if (gl1_trigger_vector_scaled[m_trigger_index[it]] == 1)
0067       {
0068         trigger_fired = true;
0069         break;
0070       }
0071     }
0072     if (!trigger_fired)
0073     {
0074       return Fun4AllReturnCodes::ABORTEVENT;
0075     }
0076   }
0077   return Fun4AllReturnCodes::EVENT_OK;
0078 }