ATLAS Offline Software
LArBarrelCalibrationCalculator.cxx
Go to the documentation of this file.
1 /*
2  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
3 */
4 
5 /********************************************************************
6 
7 NAME: CalibrationCalculator.cxx
8 PACKAGE: offline/LArCalorimeter/LArG4/LArG4Barrel
9 
10 AUTHORS: G. Unal, L. Carminati (on a template from Bill Selingman)
11 CREATED: September, 2004
12 
13 PURPOSE: This class calculates the values needed for calibration hits
14  in the barrel LAr calorimeter. This calculator is called
15  in calibration runs (see LArBarrelSDConsultant) for calibration
16  hits in the accordion (no presampler).
17 
18 UPDATES:
19 
20 ********************************************************************/
21 
22 //#define DEBUG_HITS
23 
25 
27 
28 #include "G4Step.hh"
29 #include "globals.hh"
30 #include <algorithm>
31 
32 namespace LArG4 {
33 
34  namespace Barrel {
35 
36  CalibrationCalculator::CalibrationCalculator(const std::string& name, ISvcLocator *pSvcLocator)
37  : LArCalibCalculatorSvcImp(name, pSvcLocator)
38  {
39  }
40 
42  // Initialize the geometry calculator.
43  ATH_CHECK(m_geometryCalculator.retrieve());
44  return StatusCode::SUCCESS;
45  }
46 
47  G4bool CalibrationCalculator::Process(const G4Step* step, LArG4Identifier & identifier,
48  std::vector<G4double> & energies,
49  const eCalculatorProcessing process) const
50  {
51  // Use the calculators to determine the energies and the
52  // identifier associated with this G4Step. Note that the
53  // default is to process both the energy and the ID.
54 
55  if ( process == kEnergyAndID || process == kOnlyEnergy ) {
56 #ifdef DEBUG_HITS
57  std::cout << "LArG4::Barrel::CalibrationCalculator::Process"
58  << " calling SimulationEnergies" << std::endl;
59 #endif
60  m_energyCalculator.Energies( step, energies );
61  }
62  else {
63  for (unsigned int i=0; i != 4; i++) energies.push_back( 0. );
64  }
65 
66 
67  if ( process == kEnergyAndID || process == kOnlyID ) {
68  // Calculate the identifier.
69  identifier = m_geometryCalculator->CalculateIdentifier( step );
70  }
71  else {
73  }
74 
75 
76 #ifdef DEBUG_HITS
77  G4double energy = accumulate(energies.begin(),energies.end(),0.);
78  std::cout << "LArG4::Barrel::CalibrationCalculator::Process"
79  << " ID=" << std::string(identifier)
80  << " energy=" << energy
81  << " energies=(" << energies[0]
82  << "," << energies[1]
83  << "," << energies[2]
84  << "," << energies[3] << ")"
85  << std::endl;
86 #endif
87 
88  // Check for bad result.
89  return ( identifier != LArG4Identifier() );
90  }
91 
92  } // namespace Barrel
93 
94 } // namespace LAr
LArG4Identifier
Definition: LArG4Identifier.h:121
LArCalibCalculatorSvcImp
Definition: LArCalibCalculatorSvcImp.h:12
initialize
void initialize()
Definition: run_EoverP.cxx:894
LArG4::kOnlyEnergy
@ kOnlyEnergy
Definition: LArG4EnumDefs.h:10
xAOD::identifier
identifier
Definition: UncalibratedMeasurement_v1.cxx:15
LArG4
Definition: LArWheelCalculatorEnums.h:8
LArBarrelCalibrationCalculator.h
SUSY_SimplifiedModel_PostInclude.process
string process
Definition: SUSY_SimplifiedModel_PostInclude.py:43
DetType::Barrel
@ Barrel
Definition: DetType.h:14
ParticleGun_FastCalo_ChargeFlip_Config.energy
energy
Definition: ParticleGun_FastCalo_ChargeFlip_Config.py:78
lumiFormat.i
int i
Definition: lumiFormat.py:85
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
LArG4::kOnlyID
@ kOnlyID
Definition: LArG4EnumDefs.h:10
LArG4Identifier.h
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:240
runIDAlign.accumulate
accumulate
Update flags based on parser line args.
Definition: runIDAlign.py:60
LArG4::Barrel::CalibrationCalculator::CalibrationCalculator
CalibrationCalculator(const std::string &name, ISvcLocator *pSvcLocator)
Definition: LArBarrelCalibrationCalculator.cxx:66
LArCellBinning.step
step
Definition: LArCellBinning.py:158
LArG4::kEnergyAndID
@ kEnergyAndID
Definition: LArG4EnumDefs.h:10
LArG4::eCalculatorProcessing
eCalculatorProcessing
Definition: LArG4EnumDefs.h:10