TRTNoise.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef TRT_DIGITIZATION_TRTNOISE_H
#define TRT_DIGITIZATION_TRTNOISE_H
 
#include "TRTDigit.h"
 
#include "AthenaBaseComps/AthMessaging.h"
#include "TRT_ConditionsServices/ITRT_StrawNeighbourSvc.h"
#include "TRT_ConditionsServices/ITRT_StrawStatusSummaryTool.h"
 
#include "GaudiKernel/ServiceHandle.h"
#include "GaudiKernel/ToolHandle.h"
 
#include "CLHEP/Random/RandomEngine.h"
 
#include <set>
#include <vector>
 
class Identifier;
class TRTDigCondBase;
class TRTElectronicsNoise;
class TRTElectronicsProcessing;
class TRT_ID;
 
namespace InDetDD {
  class TRT_DetectorManager;
}
 
class TRTDigSettings;

class TRTNoise : public AthMessaging {
public:
  TRTNoise( const TRTDigSettings*,
            const InDetDD::TRT_DetectorManager*,
            CLHEP::HepRandomEngine* noiseRndmEngine,
            CLHEP::HepRandomEngine* elecNoiseRndmEngine,
            CLHEP::HepRandomEngine* elecProcRndmEngine,
            CLHEP::HepRandomEngine* elecNoiseResetRndmEngine,
            TRTDigCondBase* digcond,
            TRTElectronicsProcessing * ep,
            TRTElectronicsNoise * electronicsnoise,
            const TRT_ID* trt_id,
            int UseGasMix,
            ToolHandle<ITRT_StrawStatusSummaryTool> sumTool
            );

  void appendPureNoiseToProperDigits( std::vector<TRTDigit>& digitVect,
                                      const std::set<int>& sim_hitids,
                                      CLHEP::HepRandomEngine* noiseRndmEngine) ;
 
 
  void appendCrossTalkNoiseToProperDigits(std::vector<TRTDigit>& digitVect,
                                          const std::set<Identifier>& simhitsIdentifiers,
                                          const ServiceHandle<ITRT_StrawNeighbourSvc>& m_TRTStrawNeighbourSvc,
                                          CLHEP::HepRandomEngine* noiseRndmEngine);
 
  static void sortDigits(std::vector<TRTDigit>& digitVect);
 
  ~TRTNoise();
 
private:
 
  const TRTDigSettings* m_settings{};
  const InDetDD::TRT_DetectorManager* m_detmgr{};
 
  TRTDigCondBase * m_pDigConditions{};
  TRTElectronicsProcessing * m_pElectronicsProcessing{};
  TRTElectronicsNoise * m_pElectronicsNoise{};
  CLHEP::HepRandomEngine* m_noise_randengine{};
 
  //---  Initialization:

  void InitThresholdsAndNoiseAmplitudes_and_ProduceNoiseDigitPool(CLHEP::HepRandomEngine* noiseRndmEngine,
                                                                  CLHEP::HepRandomEngine* elecNoiseRndmEngine,
                                                                  CLHEP::HepRandomEngine* elecProcRndmEngine);

  void ProduceNoiseDigitPool( const std::vector<float>& lowthresholds,
                              const std::vector<float>& noiseamps,
                              const std::vector<int>& strawType,
                              CLHEP::HepRandomEngine* noiseRndmEngine,
                              CLHEP::HepRandomEngine* elecNoiseRndmEngine,
                              CLHEP::HepRandomEngine* elecProcRndmEngine
                              );
 
  const TRT_ID* m_id_helper;
  unsigned int m_digitPoolLength;
  unsigned int m_digitPoolLength_nextaccessindex;
  std::vector<unsigned int> m_digitPool;
 
  double simulateNoiseFrequency( const unsigned int& totalNumberOfTests,
                                 const double& electronicsNoiseAmplitude,
                                 double lowthreshold = -1.0 );

  static float useLookupTable(const float& x,
                       const std::vector<float>& y_given_x,
                       const float & min_x,
                       const float & max_x);

  static void makeInvertedLookupTable( const std::vector<float>& y_given_x,
                                const float & min_x,
                                const float & max_x,
                                std::vector<float>& x_given_y,
                                float & min_y,
                                float & max_y );
  static void evolve_LT2AmpVsNL_to_include_LTfluct( std::vector<float>& nl_given_lt2na,
                                             const float & min_lt2na,
                                             const float & max_lt2na,
                                             const float relativeLTFluct,
                                             float & new_min_lt2na,
                                             float & new_max_lt2na,
                                             const unsigned int& number_new_bins );
 
  Identifier getStrawIdentifier (int hitID);
 
  int m_UseGasMix{};
  ToolHandle<ITRT_StrawStatusSummaryTool> m_sumTool; // need for Argon
 
};
 
#endif