 |
ATLAS Offline Software
|
#include <ClusterMakerTool.h>
|
| | ClusterMakerTool (const std::string &type, const std::string &name, const IInterface *parent) |
| |
| | ~ClusterMakerTool ()=default |
| |
| StatusCode | initialize () |
| |
| PixelCluster | pixelCluster (const Identifier &clusterID, const Amg::Vector2D &localPos, std::vector< Identifier > &&rdoList, const int lvl1a, std::vector< int > &&totList, const SiWidth &width, const InDetDD::SiDetectorElement *element, bool ganged, int errorStrategy, const PixelID &pixelID, bool split, double splitProb1, double splitProb2, const PixelChargeCalibCondData *calibData, const PixelCalib::PixelOfflineCalibData *offlineCalibData, const EventContext &ctx) const |
| |
| xAOD::PixelCluster * | xAODpixelCluster (xAOD::PixelCluster &cluster, const Amg::Vector2D &localPos, const std::vector< Identifier > &rdoList, const int lvl1a, const std::vector< int > &totList, const SiWidth &width, const InDetDD::SiDetectorElement *element, bool ganged, int errorStrategy, const PixelID &pixelID, bool split, double splitProb1, double splitProb2, const PixelChargeCalibCondData *calibData, const PixelCalib::PixelOfflineCalibData *offlineCalibData, const EventContext &ctx) const |
| |
| SCT_Cluster | sctCluster (const Identifier &clusterID, const Amg::Vector2D &localPos, std::vector< Identifier > &&rdoList, const SiWidth &width, const InDetDD::SiDetectorElement *element, int errorStrategy) const |
| |
| ServiceHandle< StoreGateSvc > & | evtStore () |
| | The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
|
| |
| const ServiceHandle< StoreGateSvc > & | evtStore () const |
| | The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
|
| |
| const ServiceHandle< StoreGateSvc > & | detStore () const |
| | The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc. More...
|
| |
| virtual StatusCode | sysInitialize () override |
| | Perform system initialization for an algorithm. More...
|
| |
| virtual StatusCode | sysStart () override |
| | Handle START transition. More...
|
| |
| virtual std::vector< Gaudi::DataHandle * > | inputHandles () const override |
| | Return this algorithm's input handles. More...
|
| |
| virtual std::vector< Gaudi::DataHandle * > | outputHandles () const override |
| | Return this algorithm's output handles. More...
|
| |
| Gaudi::Details::PropertyBase & | declareProperty (Gaudi::Property< T, V, H > &t) |
| |
| Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &) |
| | Declare a new Gaudi property. More...
|
| |
| Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &) |
| | Declare a new Gaudi property. More...
|
| |
| Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &) |
| |
| Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &) |
| | Declare a new Gaudi property. More...
|
| |
| Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, T &property, const std::string &doc="none") |
| | Declare a new Gaudi property. More...
|
| |
| void | updateVHKA (Gaudi::Details::PropertyBase &) |
| |
| MsgStream & | msg () const |
| |
| MsgStream & | msg (const MSG::Level lvl) const |
| |
| bool | msgLvl (const MSG::Level lvl) const |
| |
|
| template<typename ClusterType , typename IdentifierVec , typename ToTList > |
| ClusterType | makePixelCluster (const Identifier &clusterID, const Amg::Vector2D &localPos, IdentifierVec &&rdoList, const int lvl1a, ToTList &&totList, const SiWidth &width, const InDetDD::SiDetectorElement *element, bool ganged, int errorStrategy, const PixelID &pixelID, bool split, double splitProb1, double splitProb2, const PixelChargeCalibCondData *calibData, const PixelCalib::PixelOfflineCalibData *offlineCalibData, const EventContext &ctx, xAOD::PixelCluster *cluster=nullptr) const |
| |
| double | getPixelCTBPhiError (int layer, int phi, int PhiClusterSize) const |
| |
| Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &) |
| | specialization for handling Gaudi::Property<SG::VarHandleKey> More...
|
| |
| Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyArrayType &) |
| | specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
|
| |
| Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleType &) |
| | specialization for handling Gaudi::Property<SG::VarHandleBase> More...
|
| |
| Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T, V, H > &t, const SG::NotHandleType &) |
| | specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
|
| |
Definition at line 55 of file ClusterMakerTool.h.
◆ StoreGateSvc_t
◆ ClusterMakerTool()
| InDet::ClusterMakerTool::ClusterMakerTool |
( |
const std::string & |
type, |
|
|
const std::string & |
name, |
|
|
const IInterface * |
parent |
|
) |
| |
◆ ~ClusterMakerTool()
| InDet::ClusterMakerTool::~ClusterMakerTool |
( |
| ) |
|
|
default |
◆ declareGaudiProperty() [1/4]
specialization for handling Gaudi::Property<SG::VarHandleKeyArray>
Definition at line 170 of file AthCommonDataStore.h.
175 hndl.documentation());
◆ declareGaudiProperty() [2/4]
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition at line 156 of file AthCommonDataStore.h.
161 hndl.documentation());
◆ declareGaudiProperty() [3/4]
specialization for handling Gaudi::Property<SG::VarHandleBase>
Definition at line 184 of file AthCommonDataStore.h.
189 hndl.documentation());
◆ declareGaudiProperty() [4/4]
◆ declareProperty() [1/6]
Declare a new Gaudi property.
- Parameters
-
| name | Name of the property. |
| hndl | Object holding the property value. |
| doc | Documentation string for the property. |
This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.
Definition at line 245 of file AthCommonDataStore.h.
250 this->declare(hndl.
vhKey());
251 hndl.
vhKey().setOwner(
this);
253 return PBASE::declareProperty(
name,hndl,
doc);
◆ declareProperty() [2/6]
Declare a new Gaudi property.
- Parameters
-
| name | Name of the property. |
| hndl | Object holding the property value. |
| doc | Documentation string for the property. |
This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.
Definition at line 221 of file AthCommonDataStore.h.
229 return PBASE::declareProperty(
name,hndl,
doc);
◆ declareProperty() [3/6]
◆ declareProperty() [4/6]
Declare a new Gaudi property.
- Parameters
-
| name | Name of the property. |
| property | Object holding the property value. |
| doc | Documentation string for the property. |
This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.
Definition at line 333 of file AthCommonDataStore.h.
338 return PBASE::declareProperty(
name, property,
doc);
◆ declareProperty() [5/6]
Declare a new Gaudi property.
- Parameters
-
| name | Name of the property. |
| property | Object holding the property value. |
| doc | Documentation string for the property. |
This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.
Definition at line 352 of file AthCommonDataStore.h.
◆ declareProperty() [6/6]
◆ detStore()
◆ evtStore() [1/2]
◆ evtStore() [2/2]
◆ extraDeps_update_handler()
Add StoreName to extra input/output deps as needed.
use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given
◆ getPixelCTBPhiError()
| double InDet::ClusterMakerTool::getPixelCTBPhiError |
( |
int |
layer, |
|
|
int |
phi, |
|
|
int |
PhiClusterSize |
|
) |
| const |
|
private |
Definition at line 587 of file ClusterMakerTool.cxx.
597 if(phiClusterSize > 3)
return 14.6*
micrometer;
599 if(
layer == 0 &&
phi == 0)
return sigmaL0Phi0[phiClusterSize-1];
600 if(
layer == 1 &&
phi == 0)
return sigmaL1Phi0[phiClusterSize-1];
601 if(
layer == 2 &&
phi == 0)
return sigmaL2Phi0[phiClusterSize-1];
602 if(
layer == 0 &&
phi == 1)
return sigmaL0Phi1[phiClusterSize-1];
603 if(
layer == 1 &&
phi == 1)
return sigmaL1Phi1[phiClusterSize-1];
604 if(
layer == 2 &&
phi == 1)
return sigmaL2Phi1[phiClusterSize-1];
◆ initialize()
| StatusCode InDet::ClusterMakerTool::initialize |
( |
| ) |
|
◆ inputHandles()
Return this algorithm's input handles.
We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.
◆ interfaceID()
| static const InterfaceID& InDet::ClusterMakerTool::interfaceID |
( |
| ) |
|
|
inlinestatic |
◆ makePixelCluster()
template<typename ClusterType , typename IdentifierVec , typename ToTList >
| ClusterType InDet::ClusterMakerTool::makePixelCluster |
( |
const Identifier & |
clusterID, |
|
|
const Amg::Vector2D & |
localPos, |
|
|
IdentifierVec && |
rdoList, |
|
|
const int |
lvl1a, |
|
|
ToTList && |
totList, |
|
|
const SiWidth & |
width, |
|
|
const InDetDD::SiDetectorElement * |
element, |
|
|
bool |
ganged, |
|
|
int |
errorStrategy, |
|
|
const PixelID & |
pixelID, |
|
|
bool |
split, |
|
|
double |
splitProb1, |
|
|
double |
splitProb2, |
|
|
const PixelChargeCalibCondData * |
calibData, |
|
|
const PixelCalib::PixelOfflineCalibData * |
offlineCalibData, |
|
|
const EventContext & |
ctx, |
|
|
xAOD::PixelCluster * |
cluster = nullptr |
|
) |
| const |
|
private |
Definition at line 191 of file ClusterMakerTool.cxx.
216 throw std::runtime_error(
"Dynamic cast failed for design in ClusterMakerTool.cxx");
218 int rowMin = design->rows();
220 int colMin = design->columns();
222 float qRowMin = 0;
float qRowMax = 0;
223 float qColMin = 0;
float qColMax = 0;
224 std::vector<float> chargeList;
225 int nRDO=rdoList.size();
227 chargeList.reserve(nRDO);
229 for (
int i=0;
i<nRDO;
i++) {
235 std::array<InDetDD::PixelDiodeTree::CellIndexType,2> diode_idx
240 auto diode_type = design->getDiodeType(si_param);
242 && design->numberOfConnectedCells( design->readoutIdOfCell(
InDetDD::SiCellId(diode_idx[0],diode_idx[1])))>1) {
248 charge = ToT/8.0*(8000.0-1200.0)+1200.0;
250 chargeList.push_back(
charge);
254 for (
int i=0;
i<nRDO;
i++) {
259 if (calibData) {
charge=chargeList[
i]; }
275 if (col == colMin) qColMin +=
charge;
281 if (col == colMax) qColMax +=
charge;
293 if(qRowMin+qRowMax > 0) omegax = qRowMax/
float(qRowMin+qRowMax);
294 if(qColMin+qColMax > 0) omegay = qColMax/
float(qColMin+qColMax);
303 double Ax[3] = {
T(0,0),
T(1,0),
T(2,0)};
304 double Ay[3] = {
T(0,1),
T(1,1),
T(2,1)};
305 double R [3] = {
T(0,3),
T(1,3),
T(2,3)};
308 Amg::Vector3D globalPos(M[0]*Ax[0]+M[1]*Ay[0]+R[0],M[0]*Ax[1]+M[1]*Ay[1]+R[1],M[0]*Ax[2]+M[1]*Ay[2]+R[2]);
314 errorMatrix.setIdentity();
319 double eta = std::abs(globalPos.eta());
320 double zPitch =
width.z()/colRow.y();
325 throw std::runtime_error(
"Wrong helper type in ClusterMakerTool.cxx.");
328 int layer =
pid->layer_disk(clusterID);
330 switch (errorStrategy){
344 if (offlineCalibData) {
349 errorMatrix.fillSymmetric(0,0,
pow(phiError,2));
356 errorMatrix.fillSymmetric(0,0,square(phiError));
357 errorMatrix.fillSymmetric(1,1,square(
etaError));
362 errorMatrix.fillSymmetric(1,1,square(zPitch)*
ONE_TWELFTH);
381 static_assert(std::is_same_v<ClusterType, PixelCluster> ||
382 std::is_same_v<std::remove_pointer_t<ClusterType>,
384 "Not an InDet::PixelCluster or xAOD::PixelCluster");
386 return newInDetpixelCluster(newClusterID, locpos,
388 std::forward<IdentifierVec>(rdoList),
390 std::forward<ToTList>(totList),
391 std::move(chargeList),
394 std::move(errorMatrix),
401 return AddNewxAODpixelCluster(*cluster)(newClusterID,
404 std::forward<IdentifierVec>(rdoList),
406 std::forward<ToTList>(totList),
◆ msg() [1/2]
◆ msg() [2/2]
◆ msgLvl()
◆ outputHandles()
Return this algorithm's output handles.
We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.
◆ pixelCluster()
| PixelCluster InDet::ClusterMakerTool::pixelCluster |
( |
const Identifier & |
clusterID, |
|
|
const Amg::Vector2D & |
localPos, |
|
|
std::vector< Identifier > && |
rdoList, |
|
|
const int |
lvl1a, |
|
|
std::vector< int > && |
totList, |
|
|
const SiWidth & |
width, |
|
|
const InDetDD::SiDetectorElement * |
element, |
|
|
bool |
ganged, |
|
|
int |
errorStrategy, |
|
|
const PixelID & |
pixelID, |
|
|
bool |
split, |
|
|
double |
splitProb1, |
|
|
double |
splitProb2, |
|
|
const PixelChargeCalibCondData * |
calibData, |
|
|
const PixelCalib::PixelOfflineCalibData * |
offlineCalibData, |
|
|
const EventContext & |
ctx |
|
) |
| const |
◆ renounce()
◆ renounceArray()
◆ sctCluster()
Definition at line 508 of file ClusterMakerTool.cxx.
525 errorMatrix.setIdentity();
531 switch (errorStrategy){
543 else if(colRow.x() == 2){
557 else if(colRow.x() == 2){
574 double v0 = (errorMatrix)(0,0)*
w*
w;
575 double v1 = (errorMatrix)(1,1);
576 errorMatrix.fillSymmetric(0,0,cs2*
v0+sn2*v1);
577 errorMatrix.fillSymmetric(0,1,sn*sqrt(cs2)*(
v0-v1));
578 errorMatrix.fillSymmetric(1,1,sn2*
v0+cs2*v1);
581 return SCT_Cluster(clusterID, locpos, std::move(rdoList),
width, element, std::move(errorMatrix));
◆ sysInitialize()
◆ sysStart()
Handle START transition.
We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.
◆ updateVHKA()
◆ xAODpixelCluster()
| xAOD::PixelCluster * InDet::ClusterMakerTool::xAODpixelCluster |
( |
xAOD::PixelCluster & |
cluster, |
|
|
const Amg::Vector2D & |
localPos, |
|
|
const std::vector< Identifier > & |
rdoList, |
|
|
const int |
lvl1a, |
|
|
const std::vector< int > & |
totList, |
|
|
const SiWidth & |
width, |
|
|
const InDetDD::SiDetectorElement * |
element, |
|
|
bool |
ganged, |
|
|
int |
errorStrategy, |
|
|
const PixelID & |
pixelID, |
|
|
bool |
split, |
|
|
double |
splitProb1, |
|
|
double |
splitProb2, |
|
|
const PixelChargeCalibCondData * |
calibData, |
|
|
const PixelCalib::PixelOfflineCalibData * |
offlineCalibData, |
|
|
const EventContext & |
ctx |
|
) |
| const |
◆ m_detStore
◆ m_evtStore
◆ m_forceErrorStrategy1B
| bool InDet::ClusterMakerTool::m_forceErrorStrategy1B {false} |
|
private |
◆ m_pixelLorentzAngleTool
◆ m_sctLorentzAngleTool
◆ m_varHandleArraysDeclared
◆ m_vhka
The documentation for this class was generated from the following files:
int phi_index(const Identifier &id) const
Eigen::Matrix< double, Eigen::Dynamic, Eigen::Dynamic > MatrixX
Dynamic Matrix - dynamic allocation.
Scalar phi() const
phi method
float getPixelBarrelEtaError(int ibin) const
Scalar eta() const
pseudorapidity method
Eigen::Matrix< double, 2, 1 > Vector2D
virtual DetectorShape shape() const
Shape of element.
PixelClusterErrorData * getPixelClusterErrorData()
Trk::Surface & surface()
Element Surface.
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
#define ATH_MSG_VERBOSE(x)
virtual IdentifierHash identifyHash() const override final
identifier hash (inline)
Identifier wafer_id(int barrel_ec, int layer_disk, int phi_module, int eta_module) const
For a single crystal.
float getPixelEndcapPhiError(int ibin) const
const AtlasDetectorID * getIdHelper() const
Returns the id helper (inline)
virtual void setOwner(IDataHandleHolder *o)=0
float getPixelBarrelPhiError(int ibin) const
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
double phiPitch() const
Pitch (inline methods)
IdentifierHash wafer_hash(Identifier wafer_id) const
wafer hash from id
def phimod(flags, cells_name, *args, **kw)
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
Eigen::Affine3d Transform3D
int eta_index(const Identifier &id) const
float getPixelEndcapRError(int ibin) const
double R(const INavigable4Momentum *p1, const double v_eta, const double v_phi)
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
virtual void renounce()=0
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
virtual HelperType helper() const
Type of helper, defaulted to 'Unimplemented'.
float etaError(const U &p)
int getBarrelBin(double eta, int etaClusterSize, int phiClusterSize) const
double charge(const T &p)
Eigen::Matrix< double, 3, 1 > Vector3D
int getEndcapBin(int etaClusterSize, int phiClusterSize) const
Helper class to access parameters of a diode.
#define ATH_MSG_WARNING(x)
Identifier pixel_id(int barrel_ec, int layer_disk, int phi_module, int eta_module, int phi_index, int eta_index) const
For an individual pixel.
double sinStereoLocal(const Amg::Vector2D &localPos) const
Angle of strip in local frame with respect to the etaAxis.
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
constexpr double ONE_TWELFTH
float getCharge(InDetDD::PixelDiodeType type, unsigned int moduleHash, unsigned int FE, float ToT) const
static constexpr std::array< PixelDiodeTree::CellIndexType, 2 > makeCellIndex(T local_x_idx, T local_y_idx)
Create a 2D cell index from the indices in local-x (phi, row) and local-y (eta, column) direction.
This is a "hash" representation of an Identifier. This encodes a 32 bit index which can be used to lo...
virtual const SiDetectorDesign & design() const override final
access to the local description (inline):
constexpr int pow(int base, int exp) noexcept
const Amg::Transform3D & transform() const
Returns HepGeom::Transform3D by reference.
This class provides an interface to generate or decode an identifier for the upper levels of the dete...