LArBarrelCalculator.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
// LArBarrelCalculator.hh
// The Cell Identifier for the EM Barrel readout cells
 
// Adapted from code written by Gaston Parrour
// Adaptation by Sylvain Negroni
 
// 12-Jul-2002 WGS: Added LArG4Identifier.
 
// 18-03-2005  G.Unal: major revision to include new current maps
#ifndef LARBARRELCALCULATOR_H
#define LARBARRELCALCULATOR_H
 
#include "LArG4Code/LArG4Identifier.h"
#include "LArG4Code/LArCalculatorSvcImp.h"
#include "ILArBarrelGeometry.h"
#include "MapEta.h"
 
#include <stdexcept>
#include <vector>
#include <string>
#include <memory>
 
class G4Step;
class AccMap;
 
class LArG4BirksLaw;
 
class LArBarrelCalculator : public LArCalculatorSvcImp
{
public:
 
  LArBarrelCalculator(const std::string& name, ISvcLocator* pSvcLocator);
  virtual StatusCode initialize() override final;
  virtual StatusCode finalize() override final;
 
  LArBarrelCalculator (const LArBarrelCalculator&) = delete;
  LArBarrelCalculator& operator= (const LArBarrelCalculator&) = delete;
 
  void initializeForSDCreation() override final { m_geometry->initializeForSDCreation(); };
 
  virtual G4bool Process(const G4Step* a_step, std::vector<LArHitData>& hdata) const override final;
 
  // Check if the current hitTime is in-time
  virtual G4bool isInTime(G4double hitTime) const override final
  {
    return !(std::fabs(hitTime) > m_OOTcut);
  }
 
  //
  // get functions:
  //
  virtual G4float OOTcut() const override final { return m_OOTcut; }
 
private:
 
  ServiceHandle<ILArBarrelGeometry> m_geometry;
  const AccMap*   m_accmap;
  std::unique_ptr<MapEta>   m_etamap1;
  std::unique_ptr<MapEta>   m_etamap2;
  std::unique_ptr<MapEta>   m_etamap3;
 
  // RUN Options
  bool m_IflCur;
  bool m_IflMapTrans;
  bool m_IflXtalk;
 
  double m_dstep;
 
  const LArG4BirksLaw *m_birksLaw;
  bool  m_doHV;
 
  // global EMBarrel dimensions
  double m_etaMaxBarrel;
  double m_zMinBarrel;
  double m_zMaxBarrel;
  // global Accordion dimensions
  double m_rMinAccordion;
  double m_rMaxAccordion;
  // half thickness of absorber and electrode
  double m_ThickAbs;
  double m_ThickEle;
  // GU 11/06/2003  total number of cells in phi
  int m_NCellTot;
  int m_NCellMax;
  // to handle small difference (mostly phi wrapping and +-z symmetry)
  // between atlas and test beam
  bool m_testbeam;
 
 
  // Hv values
  //   0,1 = positive/negative barrel
  //   0->1023 = electrode number
  //   0->6 = eta region number (0.2 granularity) (1.2 to 1.475 in same eta bin)
  //   0,1 = below, above the electrode (according to phi in global Atlas frame)
 
  double m_hv[2][1024][7][2]{};
 
  G4bool FiducialCuts(G4double,G4double,G4double) const;
 
  void InitHV();
  double ScaleHV(double, double, double, double) const;
 
};
 
#endif // LARBARRELCALCULATOR_H