TileCisDefaultCalibTool.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef TILECALIBALG_TILECISDEFAULTCALIBTOOL_H
#define TILECALIBALG_TILECISDEFAULTCALIBTOOL_H
 
// Gaudi includes
#include "GaudiKernel/ToolHandle.h"
#include "GaudiKernel/ServiceHandle.h"
 
// Athena includes
#include "AthenaBaseComps/AthAlgTool.h"
#include "Identifier/HWIdentifier.h"
 
// Tile includes
#include "TileCalibAlgs/ITileCalibTool.h"
#include "TileConditions/TileCablingService.h"
#include "TileEvent/TileDQstatus.h"
#include "StoreGate/ReadHandleKey.h"
#include "TileEvent/TileRawChannelContainer.h"
#include "TileEvent/TileDigitsContainer.h"
#include "TileCalibBlobObjs/TileCalibUtils.h"
#include "TileMonitoring/ITileStuckBitsProbsTool.h"
#include "TileConditions/TileInfo.h"
 
#include "TString.h"
#include <stdint.h>
#include <string> 
#include <map>
 
#define NBITS 10
#define NBSTATUS 4
 
// Forward declaration
class TileHWID;
class TileCablingSvc;
class TFile;
class TileRawChannelContainer;
class TMap;
 
class TileCisDefaultCalibTool: public AthAlgTool
                             , virtual public ITileCalibTool {
 
  public:
    TileCisDefaultCalibTool(const std::string& type, const std::string& name,
        const IInterface* pParent);
    virtual ~TileCisDefaultCalibTool();
 
    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:
 
    // Bit definitions for calibration quality flag
    enum QualityType {
      includedBit = 0,
      calibratedBit = 1,
      rangeBit = 2,
      probBit = 3,
      noiseBit = 4,
      injRMSBit = 5,
      digiErrorBit = 6,
      probChi2Bit = 7,
      edgeSamp = 8,
      nextToEdgeSamp = 9,
      stuckbitBit = 10
    };
    // typedefs for maps
    // dac map: uses dac value as key
    typedef std::map<uint32_t, int> TDACIntMap;
    typedef std::map<uint32_t, double> TDACDoubleMap;
    typedef std::map<uint32_t, double>::iterator TDACDoubleMapIter;
    // adc map: uses hwid as key
    // values of these maps are dac maps
    typedef std::map<HWIdentifier, TDACIntMap*> TAdcIntMap;
    typedef std::map<HWIdentifier, TDACDoubleMap*> TAdcDoubleMap;
    typedef std::map<HWIdentifier, TDACDoubleMap*>::iterator TAdcDoubleMapIter;
 
    // Maps for number of events, mean, and rms:
    // Fill these with "running" values, i.e. values are updated event-by-event,
    TAdcIntMap m_NEvtMap;
    TAdcIntMap m_NDigitalErrorsMap;
    TAdcDoubleMap m_MeanMap;
    TAdcDoubleMap m_MeanSqMap;
 
    // Tools / storegate info
    const TileHWID* m_tileHWID;
    const TileCablingService* m_cabling;
    ServiceHandle<TileCablingSvc> m_cablingSvc;
    ToolHandle<ITileStuckBitsProbsTool> m_stuckBitsProbs{this,
      "StuckBitsProbsTool","","Tile stuck bits probabilities tool"};
 
    SG::ReadHandleKey<TileDQstatus> m_dqStatusKey;
    SG::ReadHandleKey<TileDigitsContainer> m_digitsContainerKey{this,
      "TileDigitsContainer", "TileDigitsCnt", "Tile digits container"};
    SG::ReadHandleKey<TileRawChannelContainer> m_rawChannelContainerKey{this,
      "TileRawChannelContainer", "TileRawChannelFit", "Tile raw channel container"};
    // jobOptions
    std::string m_rawChannelContainerName;
    std::string m_ntupleID;
    std::string m_DigitsContainerName;
    // whether or not to use small capacitor, by default no
    bool m_useSmallCap;
 
    // remove pedestal events: should only be necessary for runs where
    // pedestal triggers are mixed with CIS triggers
    bool m_removePed;
 
    // expected maximal ped to set overflow limit (mostly needed for low gain)
    double m_maxPed;
 
    // phase range for calibration, by default use full range
    double m_phaseMin;
    double m_phaseMax;
 
    // Maximum charge (gain-dependent) for filling dac map: not equivalent to charge range used in linear fit
    // by default, fill low gain up to 800 pC, high gain up to 12.5 pC (25 pC for demonstrator)
    double m_chargeMaxHi;
    double m_chargeMinHi;
    double m_chargeMaxLo;
    double m_chargeMinLo;
 
    // Range for linear fits: by default do high gain fit from 3-10 pC, and low gain from 300 - 700 pC
    double m_linfitMaxHi;
    double m_linfitMinHi;
    double m_linfitMaxLo;
    double m_linfitMinLo;
 
    // Range for linear fits for demonstrator: by default do high gain fit from 6-20 pC, and low gain from 300 - 700 pC
    double m_linfitMaxHiDemo;
    double m_linfitMinHiDemo;
    double m_linfitMaxLoDemo;
    double m_linfitMinLoDemo;
 
    bool m_doSampleChecking;
 
    using Tile = TileCalibUtils;
    // Results
    float (*m_calib)[Tile::MAX_DRAWER][Tile::MAX_CHAN][Tile::MAX_GAIN];
    int (*m_qflag)[Tile::MAX_DRAWER][Tile::MAX_CHAN][Tile::MAX_GAIN];
    int (*m_nDAC)[Tile::MAX_DRAWER][Tile::MAX_CHAN][Tile::MAX_GAIN]; // This is now deprecated since you can get this form the TGraph
    int (*m_nDigitalErrors)[Tile::MAX_DRAWER][Tile::MAX_CHAN][Tile::MAX_GAIN];
    float (*m_chi2)[Tile::MAX_DRAWER][Tile::MAX_CHAN][Tile::MAX_GAIN];
 
    // These will not be saved to the ntuple
    // They are for the sample check that gets written to qflag
    // Mike Miller - 4 June 2009
    int (*m_edgeSample)[Tile::MAX_DRAWER][Tile::MAX_CHAN][Tile::MAX_GAIN];
    int (*m_nextToEdgeSample)[Tile::MAX_DRAWER][Tile::MAX_CHAN][Tile::MAX_GAIN];
 
    // these arrays are used to contain information about stuck bits in adc's
    // this array pertains to the "stuck bit" quality flag; it is not written to
    // the ntuple
    int (*m_sampleBit)[Tile::MAX_DRAWER][Tile::MAX_CHAN][Tile::MAX_GAIN][NBITS];
    // This array contains information about each bit in the adc
    // it IS written into the ntuple
    unsigned short (*m_bitStatus)[Tile::MAX_DRAWER][Tile::MAX_CHAN][Tile::MAX_GAIN][NBSTATUS];
    // used to count the number of injection samples analyzed for odd bit behavior
    int (*m_numSamp)[Tile::MAX_DRAWER][Tile::MAX_CHAN][Tile::MAX_GAIN];
 
    TMap* m_scanMap;
    TMap* m_scanMapRMS;
 
    double m_defaultCalib[4] = {0., 0., 0., 0.};
    double m_dac2Charge[4] = {0., 0., 0., 0.};
    double m_chargeMin[4] = {0., 0., 0., 0.};
    double m_chargeMax[4] = {0., 0., 0., 0.};
    double m_linfitMin[4] = {0., 0., 0., 0.};
    double m_linfitMax[4] = {0., 0., 0., 0.};
    double m_maxAmp[4] = {0., 0., 0., 0.};
    std::vector<int> m_fragIDsDemonstrators;
 
    // TileInfo
    std::string m_infoName;
    const TileInfo* m_tileInfo;
    
    // Functions
    inline int chanIsConnected(int ros, int chan) {
      if (m_cabling->channel2hole(ros, chan) < 0) return 0; //negative means not connected
      return 1;
    }
 
    inline void setBit(QualityType qb, int& bitflag) {
      bitflag |= (1 << qb);
    }
 
    inline void unsetBit(QualityType qb, int& bitflag) {
      bitflag &= ~(1 << qb);
    }
 
    inline TString arrayString(int ros, int drawer, int chan, int gain) {
      TString str = "";
      str += ros;
      str += "_";
      str += drawer;
      str += "_";
      str += chan;
      str += "_";
      str += gain;
      return str;
    }
 
};
 
#endif // #ifndef TILECALIBALG_TILECISDEFAULTCALIBTOOL_H