IDScanZFinder.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2019 CERN for the benefit of the ATLAS collaboration
*/
 
// IdScan: IDScanZFinder
// ( see header-file for details )
// -------------------------------
// ATLAS Collaboration
 
#include <cmath>
#include <vector>
 
#include "IDScanZFinder.h"
#include "ZFinderConstants.h"
 
#include "TrigInDetEvent/TrigVertex.h"
#include "TrigInDetEvent/TrigSiSpacePoint.h"
#include "IRegionSelector/IRoiDescriptor.h"
 
 
using CLHEP::millimeter;
 
 
 
IDScanZFinder::IDScanZFinder( const std::string& type, const std::string& name, const IInterface* parent) : 
  Run1::IDScanZFinderInternal<TrigSiSpacePoint>(type, name), AthAlgTool( type, name, parent ),
  m_numberingTool("TrigL2LayerNumberTool")
{
  declareInterface< ITrigRun1ZFinder >( this );
 
  declareProperty( "PhiBinSize",         m_phiBinSize       = 0.2   );
  declareProperty( "UseOnlyPixels",      m_pixOnly          = false );
  declareProperty( "MinZBinSize",        m_minZBinSize      = 0.2   );
  declareProperty( "nFirstLayers",       m_nFirstLayers     = 3     );
  declareProperty( "ZBinSizeEtaCoeff",   m_zBinSizeEtaCoeff = 0.1   );
  declareProperty( "DPhiDEta",           m_dphideta         = -0.02 );
  declareProperty( "NeighborMultiplier", m_neighborMultiplier = 1.);
  declareProperty( "NumberOfPeaks",      m_numberOfPeaks    = 1     );
  declareProperty( "ChargeAware",        m_chargeAware      = false );
  declareProperty( "ZHistoPerPhi",       m_zHistoPerPhi     = false );
  declareProperty( "VrtxDistCut",        m_vrtxDistCut      = 0.    );
  declareProperty( "nVrtxSeparation",    m_nvrtxSeparation  = 0     );
  declareProperty( "VrtxMixing",         m_vrtxMixing       = 0.    );
  declareProperty( "PreferCentralZ",     m_preferCentralZ   = true  );
  declareProperty( "TrustSPProvider",    m_trustSPprovider  = true  );
  declareProperty( "FullScanMode",       m_fullScanMode     = false );
  declareProperty( "TripletMode",        m_tripletMode      = 0     ); 
  declareProperty( "TripletDZ",          m_tripletDZ        = 25.   );
  declareProperty( "TripletDK",          m_tripletDK        = 0.005 );
  declareProperty( "TripletDP",          m_tripletDP        = 0.05  );
  declareProperty( "WeightThreshold",    m_weightThreshold  = 0     );
}
 
 
StatusCode IDScanZFinder::initialize()
{
  //   NB: This should go into the InitializeInternal !!!!!!!!!
  //       NO internal settings should be changed in here, this should just 
  //       be an athena wrapper !!!
 
  //   phiBinSize is expected in degrees (~0.2); make sure it is > ZFinder_MinPhiSliceSize
  //
 
  if ( m_phiBinSize < ZFinder_MinPhiSliceSize and ! m_forcePhiBinSize ){
    ATH_MSG_WARNING("Requested PhiBinSize of "  << m_phiBinSize
                    << " degrees is smaller than the minimum allowed (" << ZFinder_MinPhiSliceSize
                    << " degrees). Set to the minimum value.");
    //  m_phiBinSize = ZFinder_MinPhiSliceSize;
  }
 
  // NB: This should go into the InitializeInternal !!!!
  if ( m_dphideta > 0 ){
    ATH_MSG_WARNING("Requested DPhiDEta of "  << m_dphideta
                    << " is positive.  Setting to its negative!");
    // m_dphideta *= -1.;
  }
 
  ATH_CHECK( m_numberingTool.retrieve() );
 
  int offsetEndcapPixels = m_numberingTool->offsetEndcapPixels();
  int maxSiliconLayerNum = m_numberingTool->maxSiliconLayerNum(); 
  //  int offsetBarrelSCT    = m_numberingTool->offsetBarrelSCT(); 
 
 
  // std::cout << "ZFinder::initialise() offset pixels " << offsetEndcapPixels 
  //        << "\toffsetBarrelSCT " << offsetBarrelSCT 
  //        << "\tmaxlayers " << maxSiliconLayerNum << std::endl;
 
  initializeInternal(maxSiliconLayerNum,offsetEndcapPixels-1);
 
  ATH_MSG_INFO("IDScanZFinder constructed:     name()  "    << name()            );
  ATH_MSG_INFO("IDScanZFinderInternal version:         "    << getVersion()      );
  ATH_MSG_INFO("IDScanZFinder::PixOnly        set to   "    << m_pixOnly         );
  ATH_MSG_INFO("IDScanZFinder::FullScanMode            "    << m_fullScanMode    );
  ATH_MSG_INFO("IDScanZFinder::PhiBinSize     set to   "    << m_phiBinSize      );
  ATH_MSG_INFO("IDScanZFinder::# of peaks to consider: "    << m_numberOfPeaks   );
  ATH_MSG_INFO("IDScanZFinder::z bin size              "    << m_minZBinSize     );
  ATH_MSG_INFO("IDScanZFinder::eta coeff               "    << m_zBinSizeEtaCoeff);
  
  ATH_MSG_INFO("IDScanZFinder::m_nFirstLayers     = " << m_nFirstLayers    );
  ATH_MSG_INFO("IDScanZFinder::m_invPhiSliceSize  = " << m_invPhiSliceSize );
  ATH_MSG_INFO("IDScanZFinder::m_phiBinSize       = " << m_phiBinSize      );
  ATH_MSG_INFO("IDScanZFinder::m_dphideta         = " << m_dphideta        );
  ATH_MSG_INFO("IDScanZFinder::m_neighborMultiplier = " << m_neighborMultiplier);
  ATH_MSG_INFO("IDScanZFinder::m_minZBinSize      = " << m_minZBinSize     );
  ATH_MSG_INFO("IDScanZFinder::m_zBinSizeEtaCoeff = " << m_zBinSizeEtaCoeff);
  ATH_MSG_INFO("IDScanZFinder::m_chargeAware      = " << m_chargeAware     );
  ATH_MSG_INFO("IDScanZFinder::m_zHistoPerPhi     = " << m_zHistoPerPhi    );
 
  ATH_MSG_INFO("IDScanZFinder::m_nvrtxSeparation  = " << m_nvrtxSeparation );
  ATH_MSG_INFO("IDScanZFinder::m_vrtxDistCut      = " << m_vrtxDistCut     );
  ATH_MSG_INFO("IDScanZFinder::m_vrtxMixing       = " << m_vrtxMixing      );
  ATH_MSG_INFO("IDScanZFinder::m_preferCentralZ   = " << m_preferCentralZ  );
 
  ATH_MSG_INFO("IDScanZFinder::m_trustSPprovider  = " << m_trustSPprovider );
 
  ATH_MSG_INFO("IDScanZFinder::m_tripletMode      = " << m_tripletMode     );
 
  ATH_MSG_INFO("IDScanZFinder::m_weigthThreshold  = " << m_weightThreshold );
 
  return StatusCode::SUCCESS;
}
 
 
TrigVertexCollection* IDScanZFinder::findZ( const std::vector<const TrigSiSpacePoint *>& spVec, const IRoiDescriptor& roi)
{
 
  TrigVertexCollection* output = new TrigVertexCollection;
 
  std::vector<vertex>* vertices = findZInternal( spVec, roi);
 
  ATH_MSG_DEBUG("roi: "    << roi);
  ATH_MSG_DEBUG("m_NumPhiSlices: " << m_NumPhiSlices);
 
 
  if ( GetInternalStatus()==-1 ) { 
    ATH_MSG_WARNING("phi of spacepoint out of range! phi=" << GetReturnValue());
    ATH_MSG_WARNING("Exiting ZFinder...");
 
  }
 
  for ( unsigned int i=0 ; i<vertices->size() ; i++ ) { 
    output->push_back( new TrigVertex( (*vertices)[i]._z, (*vertices)[i]._weight, TrigVertex::NULLID ) );
  }
 
  delete vertices;
 
  return output;
}