OccupancyMapMaker.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef PIXELCONDITIONSALGS_OCCUPANCYMAPMAKER_H
#define PIXELCONDITIONSALGS_OCCUPANCYMAPMAKER_H
 
#include "AthenaBaseComps/AthAlgorithm.h"
#include "GaudiKernel/ServiceHandle.h"
 
#include "InDetConditionsSummaryService/IInDetConditionsTool.h"
 
#include <string>
#include <sstream>
#include <vector>
#include <fstream>
#include <utility> // pair
 
class ITHistSvc;
class PixelID;
class TH2D;
class TH2C;
class TH1D;
 
class OccupancyMapMaker: public AthAlgorithm{
  public:
    OccupancyMapMaker (const std::string& name, ISvcLocator* pSvcLocator);
    ~OccupancyMapMaker();
 
    StatusCode initialize();
    StatusCode execute(const EventContext& ctx);
    StatusCode finalize();
 
  private:
    std::string getDCSIDFromPosition(int bec, int layer, int modPhi, int modEta);
 
    std::vector<std::string> &splitter(const std::string &str, char delim, std::vector<std::string> &elems);
 
    std::vector<std::string> splitter(const std::string &str, char delim); 
    StatusCode registerHistograms();
 
    const std::string histoSuffix(const int bec, const int layer);
 
  private:
    ServiceHandle <ITHistSvc> m_tHistSvc;
 
    ToolHandle<IInDetConditionsTool> m_pixelConditionsTool
    {this, "PixelConditionsSummaryTool", "PixelConditionsSummaryTool", "Tool to retrieve Pixel Conditions summary"};
 
    const PixelID *m_pixelID;
 
    // vector of modulename and vector(barrel/endcap, layer, phi, eta)
    std::vector< std::pair< std::string, std::vector<int> > > m_pixelMapping;
 
    std::string m_pixelRDOKey;
    std::vector<int> m_moduleHashList;
 
    double m_nEvents;
    //  double m_occupancyCut;
 
    TH1D* m_nEventsHist;
    TH1D* m_nEventsLBHist;
    std::vector<TH2D*> m_hitMaps;
    std::vector<TH2C*> m_noiseMaps;
    std::vector<TH1D*> m_LBdependence;
    std::vector<TH1D*> m_BCIDdependence;
    std::vector<TH1D*> m_TOTdistributions;
 
    TH1D *m_disabledModules;
    TH2D *m_overlayedPixelNoiseMap;
    TH2D *m_overlayedIBLDCNoiseMap; // Planar Double Chip
    TH2D *m_overlayedIBLSCNoiseMap; // 3D Single Chip
    //TH2D* m_overlayedDBMNoiseMap; // DBM
 
    // cuts ....
    double m_disk1ACut; // disk-1, A-side
    double m_disk2ACut; // disk-2, A-side
    double m_disk3ACut; // disk-3, A-side
    double m_disk1CCut; // disk-1, C-side
    double m_disk2CCut; // disk-2, C-side
    double m_disk3CCut; // disk-3, C-side
    double m_iblCut;    // IBL
    double m_bLayerCut; // B-Layer
    double m_layer1Cut; // Layer 1
    double m_layer2Cut; // Layer 2
    double m_dbmCut;    // DBM
 
    int m_hist_lbMax;   // max number of LB 
 
    double m_longPixelMultiplier;
    double m_gangedPixelMultiplier;
 
    //  int m_maxLVL1A;
 
    bool m_occupancyPerBC;
 
    unsigned int m_nBCReadout;
 
    unsigned int m_evt_lbMin; // lower limit for LB to be taken into account
    unsigned int m_evt_lbMax; // upper limit for LB to be taken into account
 
    bool m_calculateNoiseMaps;
 
    std::string m_mapFile;
};
 
#endif  // PIXELCONDITIONSALGS_OCCUPANCYMAPMAKERH_