#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, std::vector< G4double > &energies, const eCalculatorProcessing process=kEnergyAndID) 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.

43 {

44 }

Member Function Documentation

◆ Process()

G4bool LArG4::DM::CalibrationCalculator::Process	(	const G4Step *	step,
		LArG4Identifier &	identifier,
		std::vector< G4double > &	energies,
		const eCalculatorProcessing	process = `kEnergyAndID`
	)		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 56 of file DMCalibrationCalculator.h.

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

Public Member Functions

Private Attributes