LArG4::DM::CalibrationCalculator Class Reference

#include <DMCalibrationCalculator.h>

Inheritance diagram for LArG4::DM::CalibrationCalculator:

Collaboration diagram for LArG4::DM::CalibrationCalculator:

Public Member Functions
	CalibrationCalculator (const std::string &name, ISvcLocator *pSvcLocator)
virtual	~CalibrationCalculator ()
virtual G4bool	Process (const G4Step *step, LArG4Identifier &identifier, LArG4Identifier &identifier_sr, std::vector< double > &energies, const LArG4::eCalculatorProcessing process) const override final

Private Attributes
CaloG4::SimulationEnergies	m_energyCalculator

Detailed Description

Definition at line 34 of file DMCalibrationCalculator.h.

Constructor & Destructor Documentation

CalibrationCalculator()

LArG4::DM::CalibrationCalculator::CalibrationCalculator	(	const std::string &	name,
		ISvcLocator *	pSvcLocator )

Definition at line 36 of file DMCalibrationCalculator.cxx.

      : LArCalibCalculatorSvcImp(name,pSvcLocator)
    {
    }

~CalibrationCalculator()

LArG4::DM::CalibrationCalculator::~CalibrationCalculator ( )

virtual

Definition at line 42 of file DMCalibrationCalculator.cxx.

    {
    }

Member Function Documentation

Process()

G4bool LArG4::DM::CalibrationCalculator::Process ( const G4Step * step, LArG4Identifier & identifier, LArG4Identifier & identifier_sr, std::vector< double > & energies, const LArG4::eCalculatorProcessing process ) const

finaloverridevirtual

Definition at line 46 of file DMCalibrationCalculator.cxx.

      {
      // 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 )
        {
          m_energyCalculator.Energies( step, energies );
        }
      else
        for (unsigned int i=0; i != 4; i++) energies.push_back( 0. );
 
 
      identifier.clear();
      if ( process == kEnergyAndID  ||  process == kOnlyID )
        {
          // Calculate the identifier.
 
          // static ?
          static const double oneOverDeta = 10.;       //   1/Deta = 1./0.1 = 10.
          static const double oneOverDphi = 32./M_PI;  //   1/Dphi
 
          // Calculate the mid-point of the step, and the simple geometry variables.
 
          G4StepPoint* pre_step_point = step->GetPreStepPoint();
          G4StepPoint* post_step_point = step->GetPostStepPoint();
 
          G4ThreeVector startPoint = pre_step_point->GetPosition();
          G4ThreeVector endPoint   = post_step_point->GetPosition();
          G4ThreeVector p = (startPoint + endPoint) * 0.5;
 
          //G4double rho = p.perp();
          G4double eta = fabs( p.pseudoRapidity() );
          G4double phi = p.phi();
          if ( phi < 0. ) phi += 2.*M_PI; // Normalize for phiBin calculation
 
#ifdef DEBUG_HITS
          ATH_MSG_DEBUG("Process: rho,eta,phi " << p.perp() << " " << eta << " " << phi);
#endif
 
          // subdet = +/-4    "+" or " -" according to sign of Z in World coorinate
          G4int subdet = ( p.z() > 0.) ? 4 : -4;
          G4int phiBin = (int) ( phi * oneOverDphi );
          if (phiBin>63) phiBin=63;
          G4int type = 1;
 
          // Initialize identifier variables with (invalid) default
          // values (INT_MIN is defined in <climits>).
          G4int sampling = INT_MIN;
          G4int region = INT_MIN;
          G4int etaBin = INT_MIN;
 
          if (eta < 1.0) {
            sampling = 2;
            region   = 1;
            etaBin = (int) ( eta * oneOverDeta );
          }
          else if (eta < 1.5) {
            sampling = 2;
            region   = 2;
            etaBin = (int) ( (eta-1.) * oneOverDeta );
          }
          else if (eta < 1.6) {
            sampling = 1;
            region   = 4;
            etaBin   = 0;
          }
          else if (eta < 1.8) {
            sampling = 1;
            region   = 5;
            etaBin = (int) ( (eta-1.5) * oneOverDeta );
          }
          else if (eta < 3.2) {
            sampling = 1;
            region   = 6;
            etaBin = (int) ( (eta-1.3) * oneOverDeta );
          }
          else {
            ATH_MSG_WARNING("hit above 3.2 in eta !!!");
            return false;
          }
 
 
          // Append the cell ID to the (empty) identifier.
          identifier    << 10         // Calorimeter
                        << subdet     // LAr +/-4 where "+" or " -" according to
            // the sign of Z in World coorinate
                        << type
                        << sampling
                        << region
                        << etaBin
                        << phiBin;
 
        }
 
#ifdef DEBUG_HITS
      G4double energy = accumulate(energies.begin(),energies.end(),0.);
      ATH_MSG_DEBUG("Process:"
                << " ID=" << std::string(identifier)
                << " energy=" << energy
                << " energies=(" << energies[0]
                << "," << energies[1]
                << "," << energies[2]
                    << "," << energies[3] << ")");
#endif
 
      // Check for bad result.
      if ( identifier == LArG4Identifier() )
        return false;
 
      return true;
    }

Member Data Documentation

m_energyCalculator

CaloG4::SimulationEnergies LArG4::DM::CalibrationCalculator::m_energyCalculator

private

Definition at line 58 of file DMCalibrationCalculator.h.

---

The documentation for this class was generated from the following files:

• DMCalibrationCalculator.h
• DMCalibrationCalculator.cxx