ProjectiveCylinder.cpp

Go to the documentation of this file.

//==========================================================================
//  AIDA Detector description implementation 
//--------------------------------------------------------------------------
// Copyright (C) Organisation europeenne pour la Recherche nucleaire (CERN)
// All rights reserved.
//
// For the licensing terms see $DD4hepINSTALL/LICENSE.
// For the list of contributors see $DD4hepINSTALL/doc/CREDITS.
//
//==========================================================================
/*
 * ProjectiveCylinder.cpp
 *
 *  Created on: Oct 31, 2013
 *      Author: Christian Grefe, CERN
 */
 
#include <DDSegmentation/ProjectiveCylinder.h>
#include <DDSegmentation/SegmentationUtil.h>
 
#define _USE_MATH_DEFINES
#include <cmath>
 
namespace dd4hep {
namespace DDSegmentation {
 
using Util::thetaFromXYZ;
using Util::phiFromXYZ;
 
ProjectiveCylinder::ProjectiveCylinder(const std::string& cellEncoding) :
    CylindricalSegmentation(cellEncoding) {
    // define type and description
    _type = "ProjectiveCylinder";
    _description = "Projective segmentation in the global coordinates";
 
    // register all necessary parameters
    registerParameter("theta_bins", "Number of bins theta", _thetaBins, 1);
    registerParameter("phi_bins", "Number of bins phi", _phiBins, 1);
    registerParameter("offset_theta", "Angular offset in theta", _offsetTheta, 0., SegmentationParameter::AngleUnit, true);
    registerParameter("offset_phi", "Angular offset in phi", _offsetPhi, 0., SegmentationParameter::AngleUnit, true);
    registerIdentifier("identifier_theta", "Cell ID identifier for theta", _thetaID, "theta");
    registerIdentifier("identifier_phi", "Cell ID identifier for phi", _phiID, "phi");
}
 
 
ProjectiveCylinder::ProjectiveCylinder(const BitFieldCoder* decode) :   CylindricalSegmentation(decode) {
    // define type and description
    _type = "ProjectiveCylinder";
    _description = "Projective segmentation in the global coordinates";
 
    // register all necessary parameters
    registerParameter("theta_bins", "Number of bins theta", _thetaBins, 1);
    registerParameter("phi_bins", "Number of bins phi", _phiBins, 1);
    registerParameter("offset_theta", "Angular offset in theta", _offsetTheta, 0., SegmentationParameter::AngleUnit, true);
    registerParameter("offset_phi", "Angular offset in phi", _offsetPhi, 0., SegmentationParameter::AngleUnit, true);
    registerIdentifier("identifier_theta", "Cell ID identifier for theta", _thetaID, "theta");
    registerIdentifier("identifier_phi", "Cell ID identifier for phi", _phiID, "phi");
}
 
ProjectiveCylinder::~ProjectiveCylinder() {
 
}
 
Vector3D ProjectiveCylinder::position(const CellID& cID) const {
    return Util::positionFromRThetaPhi(1.0, theta(cID), phi(cID));
}
 
CellID ProjectiveCylinder::cellID(const Vector3D& /* localPosition */, const Vector3D& globalPosition, const VolumeID& vID) const {
        CellID cID = vID ;
    double lTheta = thetaFromXYZ(globalPosition);
    double lPhi = phiFromXYZ(globalPosition);
    _decoder->set(cID,_thetaID, positionToBin(lTheta, M_PI / (double) _thetaBins, _offsetTheta));
    _decoder->set(cID,_phiID  , positionToBin(lPhi, 2 * M_PI / (double) _phiBins, _offsetPhi));
    return cID;
}
 
double ProjectiveCylinder::theta(const CellID& cID) const {
        CellID thetaIndex = _decoder->get(cID,_thetaID);
    return M_PI * ((double) thetaIndex + 0.5) / (double) _thetaBins;
}
double ProjectiveCylinder::phi(const CellID& cID) const {
        CellID phiIndex = _decoder->get(cID,_phiID);
    return 2. * M_PI * ((double) phiIndex + 0.5) / (double) _phiBins;
}
 
 
} /* namespace DDSegmentation */
} /* namespace dd4hep */