LArBarrelCalibrationCalculator.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
*/
 
/********************************************************************
 
NAME:     CalibrationCalculator.cxx
PACKAGE:  offline/LArCalorimeter/LArG4/LArG4Barrel
 
AUTHORS:  G. Unal, L. Carminati (on a template from Bill Selingman)
CREATED:  September, 2004
 
PURPOSE:  This class calculates the values needed for calibration hits
          in the barrel LAr calorimeter. This calculator is called
          in calibration runs (see LArBarrelSDConsultant) for calibration
          hits in the accordion (no presampler).
 
UPDATES:
 
********************************************************************/
 
//#define DEBUG_HITS
 
#include "LArBarrelCalibrationCalculator.h"
 
#include "LArG4Code/LArG4Identifier.h"
 
#include "G4Step.hh"
#include "globals.hh"
#include <algorithm>
 
namespace LArG4 {
 
  namespace Barrel {
 
    CalibrationCalculator::CalibrationCalculator(const std::string& name, ISvcLocator *pSvcLocator)
      : LArCalibCalculatorSvcImp(name, pSvcLocator)
      , m_geometryCalculator("LArBarrelGeometry", name)
    {
      declareProperty("GeometryCalculator", m_geometryCalculator);
    }
 
    StatusCode CalibrationCalculator::initialize() {
      // Initialize the geometry calculator.
      ATH_CHECK(m_geometryCalculator.retrieve());
      return StatusCode::SUCCESS;
    }
 
    CalibrationCalculator::~CalibrationCalculator()
    {
    }
 
 
    G4bool CalibrationCalculator::Process(const G4Step* step, LArG4Identifier & identifier,
                                          std::vector<G4double> & energies,
                                          const eCalculatorProcessing process) const
    {
      // Use the calculators to determine the energies and the
      // identifier associated with this G4Step.  Note that the
      // default is to process both the energy and the ID.
 
      if ( process == kEnergyAndID  ||  process == kOnlyEnergy )
        {
#ifdef DEBUG_HITS
          std::cout << "LArG4::Barrel::CalibrationCalculator::Process"
                    << " calling SimulationEnergies" << std::endl;
#endif
          m_energyCalculator.Energies( step, energies );
 
          // First, get the energy.
          //      m_energy = step->GetTotalEnergyDeposit();
        }
      else
        for (unsigned int i=0; i != 4; i++) energies.push_back( 0. );
 
 
      if ( process == kEnergyAndID  ||  process == kOnlyID )
        {
          // Calculate the identifier.
          identifier = m_geometryCalculator->CalculateIdentifier( step );
        }
      else
        identifier = LArG4Identifier();
 
 
#ifdef DEBUG_HITS
      G4double energy = accumulate(energies.begin(),energies.end(),0.);
      std::cout << "LArG4::Barrel::CalibrationCalculator::Process"
                << " ID=" << std::string(identifier)
                << " energy=" << energy
                << " energies=(" << energies[0]
                << "," << energies[1]
                << "," << energies[2]
                << "," << energies[3] << ")"
                << std::endl;
#endif
 
      // Check for bad result.
      if ( identifier == LArG4Identifier() )
        return false;
 
      return true;
    }
 
  } // namespace Barrel
 
} // namespace LAr