eEmSort.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
//  eEmSort.cxx
//  TopoCore
//  algorithm to create sorted lists for eEMs, et order applied
//
#include "L1TopoAlgorithms/eEmSort.h"
#include "L1TopoEvent/TOBArray.h"
#include "L1TopoEvent/eEmTOBArray.h"
#include "L1TopoEvent/GenericTOB.h"
#include <algorithm>
 
REGISTER_ALG_TCS(eEmSort)
 
bool SortByEtLargesteEm(TCS::GenericTOB* tob1, TCS::GenericTOB* tob2)
{
   return tob1->Et() > tob2->Et();
}
 
// constructor
TCS::eEmSort::eEmSort(const std::string & name) : SortingAlg(name) {
   defineParameter( "InputWidth", 120 ); // for FW
   defineParameter( "InputWidth1stStage", 30 ); // for FW
   defineParameter( "OutputWidth", 6 );
   defineParameter( "REtaMin", 0 );
   defineParameter( "RHadMin", 0 );
   defineParameter( "WsTotMin", 0 );
   defineParameter( "MinEta", 0 );
   defineParameter( "MaxEta", 196 );
}
 
 
// destructor
TCS::eEmSort::~eEmSort() {}
 
 
TCS::StatusCode
TCS::eEmSort::initialize() {
   m_numberOfeEms = parameter("OutputWidth").value();
   m_minREta = parameter("REtaMin").value();
   m_minRHad = parameter("RHadMin").value();
   m_minWsTot = parameter("WsTotMin").value();
   return TCS::StatusCode::SUCCESS;
}
 
 
TCS::StatusCode
TCS::eEmSort::sort(const InputTOBArray & input, TOBArray & output) {
 
   const eEmTOBArray & eems = dynamic_cast<const eEmTOBArray&>(input);
 
   // fill output array with GenericTOB buildt from eEms
   for(eEmTOBArray::const_iterator eem = eems.begin(); eem!= eems.end(); ++eem ) {
 
      // isolation cut
      if ( !isocut(m_minREta, (*eem)-> Reta()) ) {continue;}
      if ( !isocut(m_minRHad, (*eem)-> Rhad()) ) {continue;}
      if ( !isocut(m_minWsTot, (*eem)-> Wstot()) ) {continue;}
      
      const GenericTOB gtob(**eem);
      output.push_back( gtob );
   }
 
   // sort
   output.sort(SortByEtLargesteEm);
   
 
   // keep only max number of eEms
   int par = m_numberOfeEms;
   unsigned int maxNumberOfeEms = std::clamp(par, 0, std::abs(par));
   if(maxNumberOfeEms>0) {
      while( output.size()> maxNumberOfeEms ) {
         if (output.size() == (maxNumberOfeEms+1)) {
            bool isAmbiguous = output[maxNumberOfeEms-1].EtDouble() == output[maxNumberOfeEms].EtDouble();
            if (isAmbiguous) { output.setAmbiguityFlag(true); }
         }
         output.pop_back();
      }
   }
   return TCS::StatusCode::SUCCESS;
}