TileLaserLinearityCalibTool.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef TILECALIBALG_TILELASERLINEARITYCALIBTOOL_H
#define TILECALIBALG_TILELASERLINEARITYCALIBTOOL_H
 
#include "AthenaBaseComps/AthAlgTool.h"
 
#include "TileCalibAlgs/ITileCalibTool.h"
#include "TileLaserDefaultCalibTool.h"
 
#include <string> 
#include <map>
 
#include "TileConditions/TileCablingService.h"
#include "TileConditions/TileCondToolEmscale.h"
#include "TileEvent/TileDQstatus.h"
#include "TileEvent/TileRawChannelContainer.h"
#include "TileEvent/TileLaserObject.h"
#include "StoreGate/ReadHandleKey.h"
 
#define NFILTERS 10
#define NDIODES 10
#define NDIODES_LASER1 4
#define NMONITORS 4
#define NPMTS 2
 
#define NGAINS 2
#define NPARTITIONS 4 
#define NDRAWERS 64
#define NCOUPLES 22 
#define NCHANNELS 48
#define NSLICES 100
#define NFIBERS 2
 
class TileRawChannelContainer;
class TileLaserObject;
class TileHWID;
class TFile;
 
class TileLaserLinearityCalibTool : public AthAlgTool, virtual public ITileCalibTool
{
 
 public:
  TileLaserLinearityCalibTool(const std::string& type, const std::string& name,const IInterface* pParent);
  virtual ~TileLaserLinearityCalibTool();
 
  virtual StatusCode initialize() override;
  virtual StatusCode initNtuple(int runNumber, int runType, TFile * rootfile) override;
  virtual StatusCode execute() override;
  virtual StatusCode finalizeCalculations() override;
  virtual StatusCode writeNtuple(int runNumber, int runType, TFile * rootfile) override;
  virtual StatusCode finalize() override;
 
 private:
 
  // Name of ntuple of the Tool
  std::string m_toolNtuple;
 
  std::string m_rawChannelContainerName;
  std::string m_laserContainerName;
  SG::ReadHandleKey<TileDQstatus> m_dqStatusKey;
 
  SG::ReadHandleKey<TileRawChannelContainer> m_rawChannelContainerKey{this,
      "TileRawChannelContainer", "TileRawChannelOpt2", "Input Tile raw channel container"};
  SG::ReadHandleKey<TileLaserObject> m_laserContainerKey{this,
      "TileLaserObject", "TileLaserObject", "Input Tile laser object"};
 
  const TileHWID* m_tileHWID;    
  const TileCablingService* m_cabling;
  ToolHandle<TileCondToolEmscale>  m_tileToolEmscale{this,
    "TileCondToolEmscale", "TileCondToolEmscale", "Tile em scale tool"};
 
  // Parameter which will end up in the ROOTuple
  //
  // Details concerning this parameters could be found on the LASER webpage :
  // 
  // http://atlas-tile-laser.web.cern.ch/
  //
  // Here we are working on linearity runs, it explains the [8] in the arrays
  // we store one value per filter position
 
  int m_toolRunNo;                   // The run number
  int m_ADC_problem;                 // Expert tag for LastROD problem noticfication
  int m_las_filter;                  // Filter wheel position
  double m_las_requ_amp;             // Requested amplitude
  double m_hrate;                    // Humidity rate (in %)
  double m_flow;                     // Gas flow in diodes box (in L/h)
  double m_head_temp;                // Temperature of the LASER head
  double m_las_time;                 // Event time
  double (*m_LG_PMT)[NPMTS];                // Mean value for box PMTs
  double (*m_LG_PMT_S)[NPMTS];              // Corresponding RMS
  double (*m_LG_diode)[NDIODES];              // Mean value for box Photodiodes
  double (*m_LG_diode_S)[NDIODES];            // Corresponding RMS
  double (*m_HG_PMT)[NPMTS];                // Mean value for box PMTs
  double (*m_HG_PMT_S)[NPMTS];              // Corresponding RMS
  double (*m_HG_diode)[NDIODES];              // Mean value for box Photodiodes
  double (*m_HG_diode_S)[NDIODES];            // Corresponding RMS
  double (*m_mean)[NPARTITIONS][NDRAWERS][NCHANNELS][NGAINS];       // Mean signal computed for this run
  double (*m_mean_S)[NPARTITIONS][NDRAWERS][NCHANNELS][NGAINS];     // Corresponding RMS
  double (*m_LG_ratio)[NPARTITIONS][NDRAWERS][NCHANNELS][NGAINS];      // Calib coefficients computed for this run (for all diodes)
  double (*m_LG_ratio_S)[NPARTITIONS][NDRAWERS][NCHANNELS][NGAINS];    // Corresponding RMS
  double (*m_LG_ratio2)[NPARTITIONS][NDRAWERS][NCHANNELS][NGAINS];      // Calib coefficients computed for this run (for all diodes)
  double (*m_LG_ratio2_S)[NPARTITIONS][NDRAWERS][NCHANNELS][NGAINS];    // Corresponding RMS
  double (*m_HG_ratio)[NPARTITIONS][NDRAWERS][NCHANNELS][NGAINS];      // Calib coefficients computed for this run (for all diodes)
  double (*m_HG_ratio_S)[NPARTITIONS][NDRAWERS][NCHANNELS][NGAINS];    // Corresponding RMS
  double (*m_HG_ratio2)[NPARTITIONS][NDRAWERS][NCHANNELS][NGAINS];      // Calib coefficients computed for this run (for all diodes)
  double (*m_HG_ratio2_S)[NPARTITIONS][NDRAWERS][NCHANNELS][NGAINS];    // Corresponding RMS
  int    (*m_entries)[NPARTITIONS][NDRAWERS][NCHANNELS][NGAINS];    // Number of LASER events collected for one channel (and a particular gain)
 
 
  // Local results - not sent to ROOTuple
 
  int m_PMT1_ADC_prev[2];
  int m_PMT2_ADC_prev[2];
 
  int m_first_filter;      // Value of the first filter position
 
  int m_last_evt_filter;   // Filter position of the previous event
  int m_n_same_filt_evts;  // How many events in a row with the same filter setting
 
  bool m_complete_turn;
 
  RunningStat* (*m_HG_diode_signal)[NDIODES];
  RunningStat* (*m_HG_PMT_signal)[NPMTS];
  RunningStat* (*m_LG_diode_signal)[NDIODES];
  RunningStat* (*m_LG_PMT_signal)[NPMTS];
  RunningStat* (*m_signal)[NPARTITIONS][NDRAWERS][NCHANNELS][NGAINS];
  RunningStat* (*m_LG_ratio_stat)[NPARTITIONS][NDRAWERS][NCHANNELS][NGAINS];
  RunningStat* (*m_LG_ratio2_stat)[NPARTITIONS][NDRAWERS][NCHANNELS][NGAINS];
  RunningStat* (*m_HG_ratio_stat)[NPARTITIONS][NDRAWERS][NCHANNELS][NGAINS];
  RunningStat* (*m_HG_ratio2_stat)[NPARTITIONS][NDRAWERS][NCHANNELS][NGAINS];
 
  // Functions
  inline int chanIsConnected(int ros, int chan) {
    if(m_cabling->channel2hole(ros,chan)<0) return 0; //negative means not connected
    return 1;
  }
};
 
 
 
#endif // #ifndef TILECALIBALG_TILELASERLINEARITYCALIBTOOL_H