ATLAS Offline Software
Public Types | Public Member Functions | Private Member Functions | Private Attributes | List of all members
SiSmearedDigitizationTool Class Reference

#include <SiSmearedDigitizationTool.h>

Inheritance diagram for SiSmearedDigitizationTool:
Collaboration diagram for SiSmearedDigitizationTool:

Public Types

typedef std::multimap< IdentifierHash, InDet::PixelCluster * > Pixel_detElement_RIO_map
 
typedef std::multimap< IdentifierHash, InDet::SCT_Cluster * > SCT_detElement_RIO_map
 

Public Member Functions

 SiSmearedDigitizationTool (const std::string &type, const std::string &name, const IInterface *parent)
 Constructor with parameters. More...
 
StatusCode initialize ()
 
StatusCode prepareEvent (const EventContext &ctx, unsigned int)
 
StatusCode processBunchXing (int bunchXing, SubEventIterator bSubEvents, SubEventIterator eSubEvents)
 
StatusCode processAllSubEvents (const EventContext &ctx)
 
StatusCode mergeEvent (const EventContext &ctx)
 
StatusCode mergeClusters (Pixel_detElement_RIO_map *cluster_map)
 
StatusCode mergeClusters (SCT_detElement_RIO_map *cluster_map)
 
StatusCode digitize (const EventContext &ctx, TimedHitCollection< SiHit > &thpcsi)
 
StatusCode createAndStoreRIOs (const EventContext &ctx)
 
StatusCode retrieveTruth ()
 
StatusCode finalize ()
 
template<typename CLUSTER >
double calculateDistance (CLUSTER *clusterA, CLUSTER *clusterB)
 
template<typename CLUSTER >
double calculateSigma (CLUSTER *clusterA, CLUSTER *clusterB)
 
template<typename CLUSTER >
ClusterInfo calculateNewCluster (CLUSTER *clusterA, CLUSTER *clusterB)
 
template<typename CLUSTER >
StatusCode FillTruthMap (PRD_MultiTruthCollection *, CLUSTER *, const TimedHitPtr< SiHit > &)
 

Private Member Functions

 SiSmearedDigitizationTool ()
 
 SiSmearedDigitizationTool (const SiSmearedDigitizationTool &)
 
SiSmearedDigitizationTooloperator= (const SiSmearedDigitizationTool &)
 

Private Attributes

ServiceHandle< IAthRNGSvcm_rndmSvc {this, "RndmSvc", "AthRNGSvc", ""}
 Random number service. More...
 
SG::ReadCondHandleKey< InDetDD::SiDetectorElementCollectionm_pixelDetEleCollKey {this, "PixelDetEleCollKey", "PixelDetectorElementCollection", "Key of SiDetectorElementCollection for Pixel"}
 
SG::ReadCondHandleKey< InDetDD::SiDetectorElementCollectionm_SCTDetEleCollKey {this, "SCTDetEleCollKey", "SCT_DetectorElementCollection", "Key of SiDetectorElementCollection for SCT"}
 
const PixelIDm_pixel_ID
 Handle to the ID helper. More...
 
const SCT_IDm_sct_ID
 Handle to the ID helper. More...
 
std::string m_randomEngineName
 Name of the random number stream. More...
 
float m_pitch_X
 
float m_pitch_Y
 
bool m_merge
 
double m_nSigma
 
bool m_useDiscSurface
 
InDet::PixelClusterContainerm_pixelClusterContainer
 the PixelClusterContainer More...
 
InDet::SCT_ClusterContainerm_sctClusterContainer
 the SCT_ClusterContainer More...
 
ServiceHandle< PileUpMergeSvcm_mergeSvc
 PileUp Merge service. More...
 
int m_HardScatterSplittingMode
 Process all SiHit or just those from signal or background events. More...
 
bool m_HardScatterSplittingSkipper
 
PRD_MultiTruthCollectionm_pixelPrdTruth {}
 
std::string m_prdTruthNamePixel
 
PRD_MultiTruthCollectionm_SCTPrdTruth {}
 
std::string m_prdTruthNameSCT
 
SiHitCollectionm_simHitColl {}
 
std::string m_inputObjectName
 
std::vector< SiHitCollection * > m_siHitCollList
 name of the sub event hit collections. More...
 
Pixel_detElement_RIO_mapm_pixelClusterMap {}
 
SCT_detElement_RIO_mapm_sctClusterMap {}
 
bool m_SmearPixel
 
bool m_emulateAtlas
 
std::string m_pixel_SiClustersName
 
std::string m_Sct_SiClustersName
 
bool m_checkSmear
 
ServiceHandle< ITHistSvc > m_thistSvc {this, "THistSvc", "THistSvc"}
 
TFile * m_outputFile
 the root file More...
 
TTree * m_currentTree
 the tree to store information from pixel and SCT (before and after smearing) More...
 
double m_x_pixel
 
double m_y_pixel
 
double m_x_exit_pixel
 
double m_y_exit_pixel
 
double m_z_exit_pixel
 
double m_x_entry_pixel
 
double m_y_entry_pixel
 
double m_z_entry_pixel
 
double m_x_pixel_global
 
double m_y_pixel_global
 
double m_z_pixel_global
 
double m_x_SCT
 
double m_x_exit_SCT
 
double m_y_exit_SCT
 
double m_z_exit_SCT
 
double m_x_entry_SCT
 
double m_y_entry_SCT
 
double m_z_entry_SCT
 
double m_x_SCT_global
 
double m_y_SCT_global
 
double m_z_SCT_global
 
double m_x_pixel_smeared
 
double m_y_pixel_smeared
 
double m_x_SCT_smeared
 
double m_Err_x_pixel
 
double m_Err_y_pixel
 
double m_Err_x_SCT
 
double m_Err_y_SCT
 

structors and AlgTool implementation

virtual bool toProcess (int bunchXing) const override
 the method this base class helps implementing More...
 
virtual StatusCode processAllSubEvents (const EventContext &ctx)=0
 dummy implementation to allow compilation while all Digitization packages are migrated to use this new interface. More...
 
virtual bool filterPassed () const override
 dummy implementation of passing filter More...
 
virtual void resetFilter () override
 dummy implementation of filter reset More...
 
Gaudi::Property< int > m_firstXing
 
Gaudi::Property< int > m_lastXing
 
Gaudi::Property< int > m_vetoPileUpTruthLinks
 
bool m_filterPassed {true}
 

Detailed Description

Definition at line 61 of file SiSmearedDigitizationTool.h.

Member Typedef Documentation

◆ Pixel_detElement_RIO_map

Definition at line 82 of file SiSmearedDigitizationTool.h.

◆ SCT_detElement_RIO_map

Definition at line 83 of file SiSmearedDigitizationTool.h.

Constructor & Destructor Documentation

◆ SiSmearedDigitizationTool() [1/3]

SiSmearedDigitizationTool::SiSmearedDigitizationTool ( const std::string &  type,
const std::string &  name,
const IInterface *  parent 
)

Constructor with parameters.

Definition at line 66 of file SiSmearedDigitizationTool.cxx.

67  :
69  m_pixel_ID(nullptr),
70  m_sct_ID(nullptr),
71  m_randomEngineName("SiSmearedDigitization"),
72  m_pitch_X(0),
73  m_pitch_Y(0),
74  m_merge(false),
75  m_nSigma(0.),
76  m_useDiscSurface(false),
77  m_pixelClusterContainer(nullptr),
78  m_sctClusterContainer(nullptr),
79  m_mergeSvc("PileUpMergeSvc",name),
82  m_prdTruthNamePixel("PRD_MultiTruthPixel"),
83  m_prdTruthNameSCT("PRD_MultiTruthSCT"),
84  m_SmearPixel(true), //true: smear pixel --- false: smear SCT
85  m_emulateAtlas(true), // error rotation for endcap SCT
86  m_checkSmear(false),
87  m_outputFile(nullptr),
88  m_currentTree(nullptr),
89  m_x_pixel(0),
90  m_y_pixel(0),
91  m_x_exit_pixel(0),
92  m_y_exit_pixel(0),
93  m_z_exit_pixel(0),
94  m_x_entry_pixel(0),
95  m_y_entry_pixel(0),
96  m_z_entry_pixel(0),
100  m_x_SCT(0),
101  m_x_exit_SCT(0),
102  m_y_exit_SCT(0),
103  m_z_exit_SCT(0),
104  m_x_entry_SCT(0),
105  m_y_entry_SCT(0),
106  m_z_entry_SCT(0),
107  m_x_SCT_global(0),
108  m_y_SCT_global(0),
109  m_z_SCT_global(0),
112  m_x_SCT_smeared(0),
113  m_Err_x_pixel(0),
114  m_Err_y_pixel(0),
115  m_Err_x_SCT(0),
116  m_Err_y_SCT(0)
117 {
118  declareProperty("RndmEngine", m_randomEngineName, "Random engine name");
119  declareProperty("InputObjectName", m_inputObjectName="PixelHits", "Input Object name" );
120  declareProperty("pitch_X", m_pitch_X);
121  declareProperty("pitch_Y", m_pitch_Y);
122  declareProperty("MergeClusters", m_merge);
123  declareProperty("Nsigma", m_nSigma);
124  declareProperty("SmearPixel", m_SmearPixel, "Enable Pixel or SCT Smearing");
125  declareProperty("PixelClusterContainerName", m_pixel_SiClustersName="PixelClusters");
126  declareProperty("SCT_ClusterContainerName", m_Sct_SiClustersName="SCT_Clusters");
127  declareProperty("CheckSmear", m_checkSmear);
128 
129  declareProperty("HardScatterSplittingMode" , m_HardScatterSplittingMode, "Control pileup & signal splitting" );
130 
131 }

◆ SiSmearedDigitizationTool() [2/3]

SiSmearedDigitizationTool::SiSmearedDigitizationTool ( )
private

◆ SiSmearedDigitizationTool() [3/3]

SiSmearedDigitizationTool::SiSmearedDigitizationTool ( const SiSmearedDigitizationTool )
private

Member Function Documentation

◆ calculateDistance()

template<typename CLUSTER >
double SiSmearedDigitizationTool::calculateDistance ( CLUSTER clusterA,
CLUSTER clusterB 
)

Definition at line 503 of file SiSmearedDigitizationTool.cxx.

503  {
504 
505  // take needed information on the first clusters
506  Amg::Vector2D intersection_a = clusterA->localPosition();
507 
508  // take needed information on the second clusters
509  Amg::Vector2D intersection_b = clusterB->localPosition();
510 
511  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: intersection_a = " << intersection_a);
512  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: intersection_b = " << intersection_b);
513 
514  double distX = intersection_a.x() - intersection_b.x();
515  double distY = intersection_a.y() - intersection_b.y();
516 
517  return sqrt(distX*distX + distY*distY);
518 }

◆ calculateNewCluster()

template<typename CLUSTER >
ClusterInfo SiSmearedDigitizationTool::calculateNewCluster ( CLUSTER clusterA,
CLUSTER clusterB 
)

Definition at line 538 of file SiSmearedDigitizationTool.cxx.

538  {
539  // take needed information on the first clusters
540  Amg::Vector2D intersection_a = clusterA->localPosition();
541  const Amg::MatrixX& clusterErr_a = clusterA->localCovariance();
542 
543  // take needed information on the second clusters
544  Amg::Vector2D intersection_b = clusterB->localPosition();
545  const Amg::MatrixX& clusterErr_b = clusterB->localCovariance();
546 
547  double sigmaX = sqrt(Amg::error(clusterErr_a,Trk::locX) * Amg::error(clusterErr_a,Trk::locX) +
548  Amg::error(clusterErr_b,Trk::locX) * Amg::error(clusterErr_b,Trk::locX));
549 
550  double sigmaY = sqrt(Amg::error(clusterErr_a,Trk::locY) * Amg::error(clusterErr_a,Trk::locY) +
551  Amg::error(clusterErr_b,Trk::locY) * Amg::error(clusterErr_b,Trk::locY));
552 
553  double interX = 0.5*(intersection_a.x()+intersection_b.x());
554  double interY = 0.5*(intersection_a.y()+intersection_b.y());
555 
556  Amg::Vector2D intersection(interX, interY);
557 
558  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: intersection = " << intersection);
559 
560  const InDet::SiWidth& siWidth_a = clusterA->width();
561  const InDet::SiWidth& siWidth_b = clusterB->width();
562 
563  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: siWidth_a = " << siWidth_a);
564  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: siWidth_b = " << siWidth_b);
565 
566  Amg::Vector2D colRow = siWidth_a.colRow() + siWidth_b.colRow();
567  Amg::Vector2D phiRz = siWidth_a.widthPhiRZ() + siWidth_b.widthPhiRZ();
568 
569  InDet::SiWidth siWidth(colRow, phiRz);
570 
571  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: siWidth = " << siWidth);
572 
573  AmgSymMatrix(2) covariance;
574  covariance.setIdentity();
575  covariance(Trk::locX,Trk::locX) = sigmaX*sigmaX;
576  covariance(Trk::locY,Trk::locY) = sigmaY*sigmaY;
577  Amg::MatrixX clusterErr = covariance;
578 
579  return ClusterInfo( intersection, siWidth, clusterErr );
580 
581 }

◆ calculateSigma()

template<typename CLUSTER >
double SiSmearedDigitizationTool::calculateSigma ( CLUSTER clusterA,
CLUSTER clusterB 
)

Definition at line 521 of file SiSmearedDigitizationTool.cxx.

521  {
522  // take needed information on the first cluster
523  const Amg::MatrixX& clusterErr_a = clusterA->localCovariance();
524 
525  // take needed information on the second clusters
526  const Amg::MatrixX& clusterErr_b = clusterB->localCovariance();
527 
528  double sigmaX = sqrt(Amg::error(clusterErr_a,Trk::locX) * Amg::error(clusterErr_a,Trk::locX) +
529  Amg::error(clusterErr_b,Trk::locX) * Amg::error(clusterErr_b,Trk::locX));
530 
531  double sigmaY = sqrt(Amg::error(clusterErr_a,Trk::locY) * Amg::error(clusterErr_a,Trk::locY) +
532  Amg::error(clusterErr_b,Trk::locY) * Amg::error(clusterErr_b,Trk::locY));
533 
534  return sqrt(sigmaX*sigmaX + sigmaY*sigmaY);
535 }

◆ createAndStoreRIOs()

StatusCode SiSmearedDigitizationTool::createAndStoreRIOs ( const EventContext &  ctx)

Definition at line 1125 of file SiSmearedDigitizationTool.cxx.

1126 {
1127  // Get PixelDetectorElementCollection
1128  const InDetDD::SiDetectorElementCollection* elementsPixel = nullptr;
1129  if (m_SmearPixel) {
1131  elementsPixel = pixelDetEle.retrieve();
1132  if (elementsPixel==nullptr) {
1133  ATH_MSG_FATAL(m_pixelDetEleCollKey.fullKey() << " could not be retrieved");
1134  return StatusCode::FAILURE;
1135  }
1136  }
1137  // Get SCT_DetectorElementCollection
1138  const InDetDD::SiDetectorElementCollection* elementsSCT = nullptr;
1139  if (not m_SmearPixel) {
1141  elementsSCT = sctDetEle.retrieve();
1142  if (elementsSCT==nullptr) {
1143  ATH_MSG_FATAL(m_SCTDetEleCollKey.fullKey() << " could not be retrieved");
1144  return StatusCode::FAILURE;
1145  }
1146  }
1147 
1148  if ( m_SmearPixel ) { // Store Pixel RIOs
1149 
1150  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: in pixel createAndStoreRIOs() ---" );
1151 
1154 
1155  for ( ; i != e; i = m_pixelClusterMap->upper_bound(i->first) ) {
1156 
1157  std::pair <Pixel_detElement_RIO_map::iterator, Pixel_detElement_RIO_map::iterator> range;
1158  range = m_pixelClusterMap->equal_range(i->first);
1159 
1161  firstDetElem = range.first;
1162 
1163  IdentifierHash waferID;
1164  waferID = firstDetElem->first;
1165 
1166  const InDetDD::SiDetectorElement* detElement = elementsPixel->getDetectorElement(waferID);
1167 
1168  InDet::PixelClusterCollection *clusterCollection = new InDet::PixelClusterCollection(waferID);
1169  clusterCollection->setIdentifier(detElement->identify());
1170 
1171  for ( Pixel_detElement_RIO_map::iterator iter = range.first; iter != range.second; ++iter ) {
1172 
1173  InDet::PixelCluster* pixelCluster = (*iter).second;
1174  pixelCluster->setHashAndIndex(clusterCollection->identifyHash(),clusterCollection->size());
1175  clusterCollection->push_back(pixelCluster);
1176  }
1177 
1178  if ( m_pixelClusterContainer->addCollection( clusterCollection, waferID ).isFailure() ) {
1179  ATH_MSG_WARNING( "Could not add collection to Identifyable container !" );
1180  }
1181  } // end for
1182 
1183  m_pixelClusterMap->clear();
1184 
1185  }
1186  else { // Store SCT RIOs
1187 
1188  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: in SCT createAndStoreRIOs() ---" );
1189 
1192 
1193  for ( ; i != e; i = m_sctClusterMap->upper_bound(i->first) ) {
1194  std::pair <SCT_detElement_RIO_map::iterator, SCT_detElement_RIO_map::iterator> range;
1195  range = m_sctClusterMap->equal_range(i->first);
1196 
1197  SCT_detElement_RIO_map::iterator firstDetElem;
1198  firstDetElem = range.first;
1199 
1200  IdentifierHash waferID;
1201  waferID = firstDetElem->first;
1202  const InDetDD::SiDetectorElement* detElement = elementsSCT->getDetectorElement(waferID);
1203 
1204  InDet::SCT_ClusterCollection *clusterCollection = new InDet::SCT_ClusterCollection(waferID);
1205  clusterCollection->setIdentifier(detElement->identify());
1206 
1207 
1208  for ( SCT_detElement_RIO_map::iterator iter = range.first; iter != range.second; ++iter ) {
1209  InDet::SCT_Cluster* sctCluster = (*iter).second;
1210  sctCluster->setHashAndIndex(clusterCollection->identifyHash(),clusterCollection->size());
1211  clusterCollection->push_back(sctCluster);
1212  }
1213 
1214  if ( m_sctClusterContainer->addCollection( clusterCollection, clusterCollection->identifyHash() ).isFailure() ) {
1215  ATH_MSG_WARNING( "Could not add collection to Identifyable container !" );
1216  }
1217 
1218  } // end for
1219 
1220  m_sctClusterMap->clear();
1221  }
1222 
1223  return StatusCode::SUCCESS;
1224 }

◆ digitize()

StatusCode SiSmearedDigitizationTool::digitize ( const EventContext &  ctx,
TimedHitCollection< SiHit > &  thpcsi 
)

Definition at line 717 of file SiSmearedDigitizationTool.cxx.

719 {
720  // Set the RNG to use for this event.
721  ATHRNG::RNGWrapper* rngWrapper = m_rndmSvc->getEngine(this, m_randomEngineName);
722  const std::string rngName = name()+m_randomEngineName;
723  rngWrapper->setSeed( rngName, ctx );
724  CLHEP::HepRandomEngine *rndmEngine = rngWrapper->getEngine(ctx);
725 
726  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: in SiSmearedDigizationTool::digitize() ---" );
727 
728  // Get PixelDetectorElementCollection
729  const InDetDD::SiDetectorElementCollection* elementsPixel = nullptr;
730  if (m_SmearPixel) {
732  elementsPixel = pixelDetEle.retrieve();
733  if (elementsPixel==nullptr) {
734  ATH_MSG_FATAL(m_pixelDetEleCollKey.fullKey() << " could not be retrieved");
735  return StatusCode::FAILURE;
736  }
737  }
738  // Get SCT_DetectorElementCollection
739  const InDetDD::SiDetectorElementCollection* elementsSCT = nullptr;
740  if (not m_SmearPixel) {
742  elementsSCT = sctDetEle.retrieve();
743  if (elementsSCT==nullptr) {
744  ATH_MSG_FATAL(m_SCTDetEleCollKey.fullKey() << " could not be retrieved");
745  return StatusCode::FAILURE;
746  }
747  }
748 
750 
751  if(m_SmearPixel) { // Smear Pixel
753  } else {// Smear SCT
755  }
756 
757  while (thpcsi.nextDetectorElement(i, e)) {
758 
759  while (i != e) {
760  m_useDiscSurface = false;
761 
762  const TimedHitPtr<SiHit>& hit(*i++);
763  int barrelEC = hit->getBarrelEndcap();
764  int layerDisk = hit->getLayerDisk();
765  int phiModule = hit->getPhiModule();
766  int etaModule = hit->getEtaModule();
767  int side = 0;
768 
769  const InDetDD::SiDetectorElement* hitSiDetElement = nullptr;
770 
771  if(m_SmearPixel) { // Smear Pixel
772  Identifier wafer_id = m_pixel_ID->wafer_id(barrelEC,layerDisk,phiModule,etaModule);
773  IdentifierHash wafer_hash = m_pixel_ID->wafer_hash(wafer_id);
774  const InDetDD::SiDetectorElement* hitSiDetElement_temp = elementsPixel->getDetectorElement(wafer_hash);
775  ATH_MSG_DEBUG("Pixel SiDetectorElement --> barrel_ec " << barrelEC << ", layer_disk " << layerDisk << ", phi_module " << phiModule << ", eta_module " << etaModule );
776  hitSiDetElement = hitSiDetElement_temp;
777  } else { // Smear SCT
778  side = hit->getSide();
779  Identifier idwafer = m_sct_ID->wafer_id(barrelEC,layerDisk,phiModule,etaModule,side);
780  IdentifierHash idhash = m_sct_ID->wafer_hash(m_sct_ID->wafer_id(idwafer));
781  const InDetDD::SiDetectorElement* hitSiDetElement_temp = elementsSCT->getDetectorElement(idhash);
782  ATH_MSG_DEBUG("SCT SiDetectorElement --> barrel_ec " << barrelEC << ", layer_disk " << layerDisk << ", phi_module " << phiModule << ", eta_module " << etaModule << ", side " << side);
783  hitSiDetElement = hitSiDetElement_temp;
784  }
785 
786  // untangling the logic: if not using custom geometry, hitSiDetElement
787  // gets dereferenced (and should be checked)
788  //
789  // if using custom geometry, hitPlanarDetElement gets dereferenced
790  // (and should be checked)
791  if (not hitSiDetElement) {
792  ATH_MSG_FATAL("hitSiDetElement is null in SiSmearedDigitizationTool:"<<__LINE__);
793  throw std::runtime_error(std::string("hitSiDetElement is null in SiSmearedDigitizationTool::digitize() "));
794  }
795 
796  if (m_SmearPixel && !(hitSiDetElement->isPixel())) continue;
797  if (!m_SmearPixel && !(hitSiDetElement->isSCT())) continue;
798 
799  IdentifierHash waferID;
800 
801  if(m_SmearPixel) { // Smear Pixel
802  waferID = m_pixel_ID->wafer_hash(hitSiDetElement->identify());
803  } else { // Smear SCT
804  waferID = m_sct_ID->wafer_hash(hitSiDetElement->identify());
805  }
806 
807  HepGeom::Point3D<double> pix_localStartPosition = hit->localStartPosition();
808  HepGeom::Point3D<double> pix_localEndPosition = hit->localEndPosition();
809 
810  pix_localStartPosition = hitSiDetElement->hitLocalToLocal3D(pix_localStartPosition);
811  pix_localEndPosition = hitSiDetElement->hitLocalToLocal3D(pix_localEndPosition);
812 
813  double localEntryX = pix_localStartPosition.x();
814  double localEntryY = pix_localStartPosition.y();
815  double localEntryZ = pix_localStartPosition.z();
816  double localExitX = pix_localEndPosition.x();
817  double localExitY = pix_localEndPosition.y();
818  double localExitZ = pix_localEndPosition.z();
819 
820  double thickness = 0.0;
821  thickness = hitSiDetElement->thickness();
822 
823  // Transform to reconstruction local coordinates (different x,y,z ordering and sign conventions compared to simulation coords)
824  if (!m_SmearPixel) { // Smear SCT
825  HepGeom::Point3D<double> sct_localStartPosition = hit->localStartPosition();
826  HepGeom::Point3D<double> sct_localEndPosition = hit->localEndPosition();
827 
828  sct_localStartPosition = hitSiDetElement->hitLocalToLocal3D(sct_localStartPosition);
829  sct_localEndPosition = hitSiDetElement->hitLocalToLocal3D(sct_localEndPosition);
830 
831  localEntryX = sct_localStartPosition.x();
832  localEntryY = sct_localStartPosition.y();
833  localEntryZ = sct_localStartPosition.z();
834  localExitX = sct_localEndPosition.x();
835  localExitY = sct_localEndPosition.y();
836  localExitZ = sct_localEndPosition.z();
837  }
838 
839  double distX = std::abs(std::abs(localExitX)-std::abs(localEntryX));
840  double distY = std::abs(std::abs(localExitY)-std::abs(localEntryY));
841 
842  if(m_SmearPixel) { // Smear Pixel
843  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: pixel start position --- " << localEntryX << ", " << localEntryY << ", " << localEntryZ );
844  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: pixel exit position --- " << localExitX << ", " << localExitY << ", " << localExitZ );
845  m_x_entry_pixel = localEntryX;
846  m_y_entry_pixel = localEntryY;
847  m_z_entry_pixel = localEntryZ;
848  m_x_exit_pixel = localExitX;
849  m_y_exit_pixel = localExitY;
850  m_z_exit_pixel = localExitZ;
851  } else { // Smear SCT
852  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: SCT start position --- " << localEntryX << ", " << localEntryY << ", " << localEntryZ );
853  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: SCT exit position --- " << localExitX << ", " << localExitY << ", " << localExitZ );
854  m_x_entry_SCT = localEntryX;
855  m_y_entry_SCT = localEntryY;
856  m_z_entry_SCT = localEntryZ;
857  m_x_exit_SCT = localExitX;
858  m_y_exit_SCT = localExitY;
859  m_z_exit_SCT = localExitZ;
860  }
861 
862  Amg::Vector2D localEntry(localEntryX,localEntryY);
863  Amg::Vector2D localExit(localExitX,localExitY);
864 
865  // the pixel positions and other needed stuff for the geometrical clustering
866  std::vector<Identifier> rdoList;
867 
868  Amg::Vector3D localDirection(localExitX-localEntryX, localExitY-localEntryY, localExitZ-localEntryZ);
869 
870  InDetDD::SiCellId entryCellId;
871  InDetDD::SiCellId exitCellId;
872 
873  // get the identifier of the entry and the exit
874  Identifier entryId = hitSiDetElement->identifierOfPosition(localEntry);
875  Identifier exitId = hitSiDetElement->identifierOfPosition(localExit);
876 
877  // now get the cellIds and check whether they're valid
878  entryCellId = hitSiDetElement->cellIdFromIdentifier(entryId);
879  exitCellId = hitSiDetElement->cellIdFromIdentifier(exitId);
880 
881  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: entryId " << entryId << " --- exitId " << exitId );
882  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: entryCellId " << entryCellId << " --- exitCellId " << exitCellId );
883 
884  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: surface " << hitSiDetElement->surface());
885 
886  // entry / exit validity
887  bool entryValid = entryCellId.isValid();
888  bool exitValid = exitCellId.isValid();
889 
890  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: entryValid? " << entryValid << " --- exitValid? " << exitValid );
891 
892  if (!entryValid && !exitValid) continue;
893 
894  // the intersecetion id and cellId of it
895  double interX = 0.5*(localEntryX+localExitX);
896  double interY = 0.5*(localEntryY+localExitY);
897 
898  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: " << (m_SmearPixel ? "pixel" : "SCT") << " inter X --- " << interX );
899  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: " << (m_SmearPixel ? "pixel" : "SCT") << " inter Y --- " << interY );
900  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: " << (m_SmearPixel ? "pixel" : "SCT") << " dist X --- " << distX );
901  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: " << (m_SmearPixel ? "pixel" : "SCT") << " dist Y --- " << distY );
902 
903  //the particle crosses at least n pixels (in x direction you have timesX, in y direction you have timesY)
904  double timesX = (m_pitch_X) ? floor(distX/m_pitch_X) : 0;
905  double timesY = (m_pitch_Y) ? floor(distY/m_pitch_Y) : 0;
906 
907  double newdistX = distX - (timesX*m_pitch_X);
908  double newdistY = distY - (timesY*m_pitch_Y);
909 
910  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: times X --- " << timesX );
911  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: times Y --- " << timesY );
912  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: new dist X --- " << newdistX );
913  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: new dist Y --- " << newdistY );
914  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: thickness --- " << thickness );
915 
916  //Probability
917 
918  double ProbY = 2*newdistY/(m_pitch_Y+newdistY);
919  double ProbX = 2*newdistX/(m_pitch_X+newdistX);
920 
921  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: ProbX --- " << ProbX );
922  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: ProbY --- " << ProbY );
923 
924  // create the errors
925  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: pitch X --- " << m_pitch_X );
926  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: pitch Y --- " << m_pitch_Y );
927 
928  double sigmaX = m_pitch_X/sqrt(12.);
929  double sigmaY = m_pitch_Y/sqrt(12.);
930 
931  int elementX = timesX+1;
932  int elementY = timesY+1;
933 
934  if(m_SmearPixel) {
935  if (CLHEP::RandFlat::shoot(rndmEngine, 0.0, 1.0) < ProbY) { // number of crossed pixel is (timesY+1)+1
936  sigmaY = (double)(timesY+2)*m_pitch_Y/sqrt(12.);
937  elementY++;
938  } else // number of crossed pixel is (timesY+1)
939  sigmaY = (double)(timesY+1)*m_pitch_Y/sqrt(12.);
940 
941  if (CLHEP::RandFlat::shoot(rndmEngine, 0.0, 1.0) < ProbX) { // number of crossed pixel is (timesY+1)+1
942  sigmaX = (double)(timesX+2)*m_pitch_X/sqrt(12.);
943  elementX++;
944  } else // number of crossed pixel is (timesY+1)
945  sigmaX = (double)(timesX+1)*m_pitch_X/sqrt(12.);
946  }
947 
948  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool " << (m_SmearPixel ? "pixel" : "SCT") << " sigma X --- " << sigmaX);
949  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool " << (m_SmearPixel ? "pixel" : "SCT") << " sigma Y --- " << sigmaY);
950 
951 
952  double temp_X = interX;
953  double temp_Y = interY;
954 
955  Amg::Vector2D intersection(interX, interY);
956 
957  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: Intersection after smearing: " << intersection);
958 
959  Identifier intersectionId;
960  intersectionId = hitSiDetElement->identifierOfPosition(intersection);
961 
962  rdoList.push_back(intersectionId);
963  InDetDD::SiCellId currentCellId = hitSiDetElement->cellIdFromIdentifier(intersectionId);
964 
965  if (!currentCellId.isValid()) continue;
966 
967  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: Intersection Id = " << intersectionId << " --- currentCellId = " << currentCellId );
968 
969  if(m_SmearPixel) { // Smear Pixel --> Create Pixel Cluster
970 
971  double lengthX = (m_pitch_X) ? elementX*m_pitch_X : 1.;
972  double lengthY = (m_pitch_Y) ? elementY*m_pitch_Y : 1.;
973 
974  if (m_pitch_X == 0. || m_pitch_Y == 0.)
975  ATH_MSG_WARNING( "--- SiSmearedDigitizationTool: pitchX and/or pitchY are 0. Cluster length is forced to be 1. mm");
976 
977  InDet::SiWidth siWidth(Amg::Vector2D(elementX,elementY), Amg::Vector2D(lengthX, lengthY));
978 
980 
981  AmgSymMatrix(2) covariance;
982  covariance.setIdentity();
983  covariance(Trk::locX,Trk::locX) = sigmaX*sigmaX;
984  covariance(Trk::locY,Trk::locY) = sigmaY*sigmaY;
985 
986  // create the cluster
987  pixelCluster = new InDet::PixelCluster(intersectionId,
988  intersection,
989  std::vector<Identifier>(rdoList),
990  siWidth,
991  hitSiDetElement,
992  Amg::MatrixX(covariance));
993  m_pixelClusterMap->insert(std::pair<IdentifierHash, InDet::PixelCluster* >(waferID, pixelCluster));
994 
995  if (FillTruthMap(m_pixelPrdTruth, pixelCluster, hit).isFailure()) {
996  ATH_MSG_FATAL ( "FillTruthMap() for pixel failed!" );
997  return StatusCode::FAILURE;
998  }
999 
1000  if (m_checkSmear) {
1001 
1002  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: pixelCluster --> " << *pixelCluster);
1003  //Take info to store in the tree
1004  m_x_pixel = temp_X;
1005  m_y_pixel = temp_Y;
1006 
1007  m_x_pixel_smeared = (pixelCluster->localPosition()).x();
1008  m_y_pixel_smeared = (pixelCluster->localPosition()).y();
1009 
1010  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: BEFORE SMEARING LocalPosition --> X = " << m_x_pixel << " Y = " << m_y_pixel);
1011  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: LocalPosition --> X = " << m_x_pixel_smeared << " Y = " << m_y_pixel_smeared);
1012 
1013  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: GlobalPosition --> X = " << (pixelCluster->globalPosition()).x() << " Y = " << (pixelCluster->globalPosition()).y());
1014  m_x_pixel_global = (pixelCluster->globalPosition()).x();
1015  m_y_pixel_global = (pixelCluster->globalPosition()).y();
1016  m_z_pixel_global = (pixelCluster->globalPosition()).z();
1017 
1018  m_Err_x_pixel = Amg::error(pixelCluster->localCovariance(), Trk::locX);
1019  m_Err_y_pixel = Amg::error(pixelCluster->localCovariance(), Trk::locY);
1020 
1021  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: Error --> X = " << m_Err_x_pixel << " Y = " << m_Err_y_pixel);
1022 
1023  m_currentTree -> Fill();
1024  } // End of Pix smear check
1025 
1026  } else { // Smear SCT --> Create SCT Cluster
1027 
1028  // prepare the clusters
1029  InDet::SCT_Cluster * sctCluster = nullptr;
1030 
1031  // Pixel Design needed -------------------------------------------------------------
1032  const InDetDD::SCT_ModuleSideDesign* design_sct;
1033 
1034  design_sct = dynamic_cast<const InDetDD::SCT_ModuleSideDesign*>(&hitSiDetElement->design());
1035 
1036  if (!design_sct) {
1037  ATH_MSG_INFO ( "Could not get design"<< design_sct) ;
1038  continue;
1039  }
1040 
1041  // Find length of strip at centre
1042  double clusterWidth = rdoList.size()*hitSiDetElement->phiPitch(intersection);
1043  const std::pair<InDetDD::SiLocalPosition, InDetDD::SiLocalPosition> ends(design_sct->endsOfStrip(intersection));
1044  double stripLength = std::abs(ends.first.xEta()-ends.second.xEta());
1045 
1046  InDet::SiWidth siWidth(Amg::Vector2D(int(rdoList.size()),1),
1047  Amg::Vector2D(clusterWidth,stripLength) );
1048 
1049  const Amg::Vector2D& colRow = siWidth.colRow();
1050 
1051  AmgSymMatrix(2) mat;
1052  mat.setIdentity();
1053  mat(Trk::locX,Trk::locX) = sigmaX*sigmaX;
1054  mat(Trk::locY,Trk::locY) = hitSiDetElement->length()*hitSiDetElement->length()/12.;
1055 
1056  // single strip - resolution close to pitch/sqrt(12)
1057  // two-strip hits: better resolution, approx. 40% lower
1058 
1059  InDetDD::DetectorShape elShape = hitSiDetElement->design().shape();
1060  if(m_emulateAtlas && elShape == InDetDD::Trapezoid)
1061  { // rotation for endcap SCT
1062 
1063  if(colRow.x() == 1) {
1064  mat(Trk::locX,Trk::locX) = pow(siWidth.phiR(),2)/12;
1065  }
1066  else if(colRow.x() == 2) {
1067  mat(Trk::locX,Trk::locX) = pow(0.27*siWidth.phiR(),2)/12;
1068  }
1069  else {
1070  mat(Trk::locX,Trk::locX) = pow(siWidth.phiR(),2)/12;
1071  }
1072 
1073  mat(Trk::locY,Trk::locY) = pow(siWidth.z()/colRow.y(),2)/12;
1074  double sn = hitSiDetElement->sinStereoLocal(intersection);
1075  double sn2 = sn*sn;
1076  double cs2 = 1.-sn2;
1077  double w = hitSiDetElement->phiPitch(intersection)/hitSiDetElement->phiPitch();
1078  double v0 = mat(Trk::locX, Trk::locX)*w*w;
1079  double v1 = mat(Trk::locY, Trk::locY);
1080  mat(Trk::locX, Trk::locX) = (cs2*v0+sn2*v1);
1081  mat(Trk::locY, Trk::locX) = (sn*sqrt(cs2)*(v0-v1));
1082  mat(Trk::locY, Trk::locY) = (sn2*v0+cs2*v1);
1083  } // End of rotation endcap SCT
1084 
1085 
1086  sctCluster = new InDet::SCT_Cluster(intersectionId,
1087  intersection,
1088  std::vector<Identifier>(rdoList),
1089  siWidth,
1090  hitSiDetElement,
1091  Amg::MatrixX(mat));
1092 
1093  m_sctClusterMap->insert(std::pair<IdentifierHash, InDet::SCT_Cluster* >(waferID, sctCluster));
1094 
1095  if (FillTruthMap(m_SCTPrdTruth, sctCluster, hit).isFailure()) {
1096  ATH_MSG_FATAL ( "FillTruthMap() for SCT failed!" );
1097  return StatusCode::FAILURE;
1098  }
1099 
1100  if (m_checkSmear) {
1101 
1102  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: SCT_Cluster --> " << *sctCluster);
1103 
1104  //Take info to store in the tree
1105  m_x_SCT = m_useDiscSurface ? temp_Y : temp_X;
1106  m_x_SCT_smeared = (sctCluster->localPosition()).x();
1107 
1108  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: BEFORE SMEARING LocalPosition --> X = " << m_x_SCT );
1109  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: LocalPosition --> X = " << m_x_SCT_smeared );
1110 
1111  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: GlobalPosition --> X = " << (sctCluster->globalPosition()).x() << " Y = " << (sctCluster->globalPosition()).y());
1112  m_x_SCT_global = (sctCluster->globalPosition()).x();
1113  m_y_SCT_global = (sctCluster->globalPosition()).y();
1114  m_z_SCT_global = (sctCluster->globalPosition()).z();
1115 
1116  m_currentTree -> Fill();
1117  } // End smear checking
1118  } // End SCT
1119  } // while
1120  } //while
1121  return StatusCode::SUCCESS;
1122 }

◆ FillTruthMap()

template<typename CLUSTER >
StatusCode SiSmearedDigitizationTool::FillTruthMap ( PRD_MultiTruthCollection map,
CLUSTER cluster,
const TimedHitPtr< SiHit > &  hit 
)

Definition at line 479 of file SiSmearedDigitizationTool.cxx.

479  {
480 
481  ATH_MSG_DEBUG("Truth map filling with cluster " << *cluster << " and link = " << hit->particleLink());
482  if (hit->particleLink().isValid()){
483  if (!HepMC::ignoreTruthLink(hit->particleLink(), m_vetoPileUpTruthLinks)) {
484  map->insert(std::make_pair(cluster->identify(), hit->particleLink()));
485  ATH_MSG_DEBUG("Truth map filled with cluster " << *cluster << " and link = " << hit->particleLink());
486  }
487  }else{
488  ATH_MSG_DEBUG("Particle link NOT valid!! Truth map NOT filled with cluster" << cluster << " and link = " << hit->particleLink());
489  }
490 
491  return StatusCode::SUCCESS;
492 }

◆ filterPassed()

virtual bool PileUpToolBase::filterPassed ( ) const
inlineoverridevirtualinherited

dummy implementation of passing filter

Definition at line 49 of file PileUpToolBase.h.

49 { return m_filterPassed; }

◆ finalize()

StatusCode SiSmearedDigitizationTool::finalize ( )

Definition at line 243 of file SiSmearedDigitizationTool.cxx.

244 {
245 
246  if (m_checkSmear){
247  m_outputFile->cd();
248  m_currentTree->Write();
249  m_outputFile->Close();
250  ATH_MSG_DEBUG ( "SiSmearedDigitizationTool : Writing Tree" );
251 
252  }
253 
254  ATH_MSG_DEBUG ( "SiSmearedDigitizationTool : finalize()" );
255 
256 
257  return StatusCode::SUCCESS;
258 
259 }

◆ initialize()

StatusCode SiSmearedDigitizationTool::initialize ( )
virtual

Reimplemented from PileUpToolBase.

Definition at line 135 of file SiSmearedDigitizationTool.cxx.

136 {
137 
138  ATH_MSG_DEBUG ( "SiSmearedDigitizationTool::initialize()" );
139 
140  //locate the AtRndmGenSvc and initialize our local ptr
141  if (!m_rndmSvc.retrieve().isSuccess())
142  {
143  ATH_MSG_ERROR ( "Could not find given RndmSvc" );
144  return StatusCode::FAILURE;
145  }
146 
147  if (detStore()->retrieve(m_pixel_ID, "PixelID").isFailure()) {
148  ATH_MSG_ERROR ( "Could not get Pixel ID helper" );
149  return StatusCode::FAILURE;
150  }
151 
152  if (not m_SmearPixel){ // Smear SCT
153  if (detStore()->retrieve(m_sct_ID, "SCT_ID").isFailure()) {
154  ATH_MSG_ERROR ( "Could not get SCT ID helper" );
155  return StatusCode::FAILURE;
156  }
157 
158  m_inputObjectName="SCT_Hits"; // Set the input object name
159  }
160 
161  // Initialize ReadCondHandleKeys
164 
165  if (m_inputObjectName.empty())
166  {
167  ATH_MSG_FATAL ( "Property InputObjectName not set !" );
168  return StatusCode::FAILURE;
169  }
170  else
171  {
172  ATH_MSG_DEBUG ( "Input objects: '" << m_inputObjectName << "'" );
173  }
174 
175  //locate the PileUpMergeSvc and initialize our local ptr
176  if (!m_mergeSvc.retrieve().isSuccess()) {
177  ATH_MSG_ERROR ( "Could not find PileUpMergeSvc" );
178  return StatusCode::FAILURE;
179  }
180 
181  if (m_checkSmear){
182 
183  // get THistSvc
184  ATH_CHECK(m_thistSvc.retrieve());
185 
186  if (m_SmearPixel){
187  m_outputFile = new TFile("CheckSmearing_Pixel.root","RECREATE");
188  m_currentTree = new TTree("PixelTree","Check smearing Pixel");
189  m_currentTree->Branch("pixel_X" , &m_x_pixel , "m_x_pixel/D");
190  m_currentTree->Branch("pixel_Y" , &m_y_pixel , "m_y_pixel/D");
191  m_currentTree->Branch("pixel_X_exit" , &m_x_exit_pixel , "m_x_exit_pixel/D");
192  m_currentTree->Branch("pixel_Y_exit" , &m_y_exit_pixel , "m_y_exit_pixel/D");
193  m_currentTree->Branch("pixel_Z_exit" , &m_z_exit_pixel , "m_z_exit_pixel/D");
194  m_currentTree->Branch("pixel_X_entry" , &m_x_entry_pixel , "m_x_entry_pixel/D");
195  m_currentTree->Branch("pixel_Y_entry" , &m_y_entry_pixel , "m_y_entry_pixel/D");
196  m_currentTree->Branch("pixel_Z_entry" , &m_z_entry_pixel , "m_z_entry_pixel/D");
197  m_currentTree->Branch("pixel_X_global" , &m_x_pixel_global , "m_x_pixel_global/D");
198  m_currentTree->Branch("pixel_Y_global" , &m_y_pixel_global , "m_y_pixel_global/D");
199  m_currentTree->Branch("pixel_Z_global" , &m_z_pixel_global , "m_z_pixel_global/D");
200  m_currentTree->Branch("pixel_X_smear" , &m_x_pixel_smeared , "m_x_pixel_smeared/D");
201  m_currentTree->Branch("pixel_Y_smear" , &m_y_pixel_smeared , "m_y_pixel_smeared/D");
202  m_currentTree->Branch("pixel_Err_X" , &m_Err_x_pixel , "m_Err_x_pixel/D");
203  m_currentTree->Branch("pixel_Err_Y" , &m_Err_y_pixel , "m_Err_y_pixel/D");
204 
205  if( m_thistSvc->regTree("PixelTree",m_currentTree).isFailure()) {
206  ATH_MSG_ERROR("Cannot register Ttree");
207  return StatusCode::FAILURE;
208  }else{
209  ATH_MSG_DEBUG ( "Ttree registered" );
210  }
211  }
212  else{
213  m_outputFile = new TFile("CheckSmearing_SCT.root","RECREATE");
214  m_currentTree = new TTree("SCT_Tree","Check smearing SCT");
215  m_currentTree->Branch("SCT_X" , &m_x_SCT , "m_x_SCT/D");
216  m_currentTree->Branch("SCT_exit_X" , &m_x_exit_SCT , "m_x_exit_SCT/D");
217  m_currentTree->Branch("SCT_exit_Y" , &m_y_exit_SCT , "m_y_exit_SCT/D");
218  m_currentTree->Branch("SCT_exit_Z" , &m_z_exit_SCT , "m_z_exit_SCT/D");
219  m_currentTree->Branch("SCT_entry_X" , &m_x_entry_SCT , "m_x_entry_SCT/D");
220  m_currentTree->Branch("SCT_entry_Y" , &m_y_entry_SCT , "m_y_entry_SCT/D");
221  m_currentTree->Branch("SCT_entry_Z" , &m_z_entry_SCT , "m_z_entry_SCT/D");
222  m_currentTree->Branch("SCT_X_global" , &m_x_SCT_global , "m_x_SCT_global/D");
223  m_currentTree->Branch("SCT_Y_global" , &m_y_SCT_global , "m_y_SCT_global/D");
224  m_currentTree->Branch("SCT_Z_global" , &m_z_SCT_global , "m_z_SCT_global/D");
225  m_currentTree->Branch("SCT_X_smear" , &m_x_SCT_smeared , "m_x_SCT_smeared/D");
226  m_currentTree->Branch("SCT_Err_X" , &m_Err_x_SCT , "m_Err_x_SCT/D");
227 
228  if( m_thistSvc->regTree("SCT_Tree",m_currentTree).isFailure()) {
229  ATH_MSG_ERROR("Cannot register Ttree");
230  return StatusCode::FAILURE;
231  }else{
232  ATH_MSG_DEBUG ( "Ttree registered" );
233  }
234  }
235 
236 
237  }
238 
239  return StatusCode::SUCCESS;
240 }

◆ mergeClusters() [1/2]

StatusCode SiSmearedDigitizationTool::mergeClusters ( Pixel_detElement_RIO_map cluster_map)

Definition at line 583 of file SiSmearedDigitizationTool.cxx.

584 {
585  // The idea is first to check how many cluster we have to merge and then merge them.
586  // The loop is done until there aren't any clusters to merge
587 
588  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: in mergeClusters() using PixelClusters --- ");
589 
590  Pixel_detElement_RIO_map::iterator i = cluster_map->begin();
591  Pixel_detElement_RIO_map::iterator e = cluster_map->end();
592 
593  for (; i != e; i = cluster_map->upper_bound(i->first)){
594  IdentifierHash current_id = i->first;
595  // Check if clusters with current_id have been already considered
596 
597  bool NewMerge = true;
598 
599  while (NewMerge) {
600  NewMerge = false;
601  std::pair <Pixel_detElement_RIO_map::iterator, Pixel_detElement_RIO_map::iterator> range = cluster_map->equal_range(current_id);
602 
603  for ( Pixel_detElement_RIO_map::iterator iter = range.first; iter != range.second; ++iter){
604  for (Pixel_detElement_RIO_map::iterator inner_iter = std::next(iter); inner_iter != range.second; ++inner_iter){
605 
606  double dist = calculateDistance((*iter).second,(*inner_iter).second);
607  double sigma = calculateSigma((*iter).second,(*inner_iter).second);
608 
609  if( dist <= m_nSigma * sigma) {
610 
611  std::vector<Identifier> rdoList;
612 
614  InDet::SiWidth siWidth;
615  Amg::MatrixX clusterErr;
616  std::tie( intersection, siWidth, clusterErr ) = calculateNewCluster( iter->second, inner_iter->second );
617 
618  const InDetDD::SiDetectorElement* hitSiDetElement = (((*inner_iter).second)->detectorElement());
619  Identifier intersectionId = hitSiDetElement->identifierOfPosition(intersection);
620 
621  rdoList.push_back(intersectionId);
622 
623  InDetDD::SiCellId currentCellId = hitSiDetElement->cellIdFromIdentifier(intersectionId);
624 
625  if ( !currentCellId.isValid() ) {
626  continue;
627  }
628 
630  intersection,
631  std::move(rdoList),
632  siWidth,
633  hitSiDetElement,
634  std::move(clusterErr));
635  ((*inner_iter).second) = pixelCluster;
636 
637  cluster_map->erase(iter);
638  NewMerge = true;
639  goto REPEAT_LOOP;
640  }
641  }
642  }
643  REPEAT_LOOP: ;
644  }
645  }
646 
647  return StatusCode::SUCCESS;
648 }

◆ mergeClusters() [2/2]

StatusCode SiSmearedDigitizationTool::mergeClusters ( SCT_detElement_RIO_map cluster_map)

◆ mergeEvent()

StatusCode SiSmearedDigitizationTool::mergeEvent ( const EventContext &  ctx)

Definition at line 494 of file SiSmearedDigitizationTool.cxx.

494  {
495 
496  if (m_SmearPixel)
498  else
500 }

◆ operator=()

SiSmearedDigitizationTool& SiSmearedDigitizationTool::operator= ( const SiSmearedDigitizationTool )
private

◆ prepareEvent()

StatusCode SiSmearedDigitizationTool::prepareEvent ( const EventContext &  ctx,
unsigned int   
)

Definition at line 261 of file SiSmearedDigitizationTool.cxx.

262 {
263 
264  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: in pixel prepareEvent() ---" );
265 
266  m_siHitCollList.clear();
268 
269  return StatusCode::SUCCESS;
270 }

◆ processAllSubEvents() [1/2]

StatusCode SiSmearedDigitizationTool::processAllSubEvents ( const EventContext &  ctx)
virtual

Reimplemented from PileUpToolBase.

Definition at line 314 of file SiSmearedDigitizationTool.cxx.

314  {
315 
316  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: in pixel processAllSubEvents() ---" );
317 
318  InDet::SiClusterContainer* symSiContainer=nullptr;
319 
320  if(m_SmearPixel){ // Smear Pixel
322 
324  ATH_MSG_FATAL( "[ --- ] Could not create PixelClusterContainer");
325  return StatusCode::FAILURE;
326  }
327 
328  // --------------------------------------
329  // PixelCluster container registration
330 
331  m_pixelClusterContainer->cleanup();
332  if ((evtStore()->record(m_pixelClusterContainer, m_pixel_SiClustersName)).isFailure()) {
333  if ((evtStore()->retrieve(m_pixelClusterContainer, m_pixel_SiClustersName)).isFailure()) {
334  ATH_MSG_FATAL("[ hitproc ] Error while registering PixelCluster container");
335  return StatusCode::FAILURE;
336  }
337  }
338 
339  // symlink the Pixel Container
340  // Pixel
341 
342  if ((evtStore()->symLink(m_pixelClusterContainer,symSiContainer)).isFailure()) {
343  ATH_MSG_FATAL( "[ --- ] PixelClusterContainer could not be symlinked to SiClusterContainter in StoreGate !" );
344  return StatusCode::FAILURE;
345  } else {
346  ATH_MSG_INFO( "[ hitproc ] PixelClusterContainer symlinked to SiClusterContainer in StoreGate" );
347  }
348 
349  }else{ // Smear SCT
351 
352  if(!m_sctClusterContainer) {
353  ATH_MSG_FATAL( "[ --- ] Could not create SCT_ClusterContainer");
354  return StatusCode::FAILURE;
355  }
356 
357  // --------------------------------------
358  // SCT_Cluster container registration
359  m_sctClusterContainer->cleanup();
360  if ((evtStore()->record(m_sctClusterContainer, m_Sct_SiClustersName)).isFailure()) {
361  ATH_MSG_FATAL("[ hitproc ] Error while registering SCT_Cluster container");
362  return StatusCode::FAILURE;
363  }
364  // symlink the SCT Container
365  // SCT
366 
367  if ((evtStore()->symLink(m_sctClusterContainer,symSiContainer)).isFailure()) {
368  ATH_MSG_FATAL( "[ --- ] SCT_ClusterContainer could not be symlinked to SiClusterContainter in StoreGate !" );
369  return StatusCode::FAILURE;
370  } else {
371  ATH_MSG_DEBUG( "[ hitproc ] SCT_ClusterContainer symlinked to SiClusterContainer in StoreGate" );
372  }
373 
374  }
375 
376  if (retrieveTruth().isFailure()) {
377  ATH_MSG_FATAL ( "retrieveTruth() failed!" );
378  return StatusCode::FAILURE;
379  }
380 
381  // get the container(s)
382  using TimedHitCollList = PileUpMergeSvc::TimedList<SiHitCollection>::type;
383 
385 
386  //this is a list<pair<time_t, DataLink<SCTUncompressedHitCollection> > >
387  TimedHitCollList hitCollList;
388  unsigned int numberOfSimHits(0);
389  if ( !(m_mergeSvc->retrieveSubEvtsData(m_inputObjectName, hitCollList, numberOfSimHits).isSuccess()) && hitCollList.empty() ) {
390  ATH_MSG_ERROR ( "Could not fill TimedHitCollList" );
391  return StatusCode::FAILURE;
392  } else {
393  ATH_MSG_DEBUG ( hitCollList.size() << " SiHitCollections with key " << m_inputObjectName << " found" );
394  }
395 
396  // Define Hit Collection
397  TimedHitCollection<SiHit> thpcsi(numberOfSimHits);
398 
399  //now merge all collections into one
400  TimedHitCollList::iterator iColl(hitCollList.begin());
401  TimedHitCollList::iterator endColl(hitCollList.end() );
402 
404  // loop on the hit collections
405  while ( iColl != endColl ) {
407  if (m_HardScatterSplittingMode == 1 && m_HardScatterSplittingSkipper ) { ++iColl; continue; }
409  const SiHitCollection* p_collection(iColl->second);
410  thpcsi.insert(iColl->first, p_collection);
411  ATH_MSG_DEBUG ( "SiHitCollection found with " << p_collection->size() << " hits" );
412  ++iColl;
413  }
414 
415  // Process the Hits straw by straw: get the iterator pairs for given straw
416  if(this->digitize(ctx, thpcsi).isFailure()) {
417  ATH_MSG_FATAL ( "digitize method failed!" );
418  return StatusCode::FAILURE;
419  }
420 
421  if (m_merge)
422  if(this->mergeEvent(ctx).isFailure()) {
423  ATH_MSG_FATAL ( "merge method failed!" );
424  return StatusCode::FAILURE;
425  }
426 
427  if (createAndStoreRIOs(ctx).isFailure()) {
428  ATH_MSG_FATAL ( "createAndStoreRIOs() failed!" );
429  return StatusCode::FAILURE;
430  }
431  else {
432  ATH_MSG_DEBUG ( "createAndStoreRIOs() succeeded" );
433  }
434 
435  return StatusCode::SUCCESS;
436 
437 }

◆ processAllSubEvents() [2/2]

virtual StatusCode IPileUpTool::processAllSubEvents
inherited

dummy implementation to allow compilation while all Digitization packages are migrated to use this new interface.

◆ processBunchXing()

StatusCode SiSmearedDigitizationTool::processBunchXing ( int  bunchXing,
SubEventIterator  bSubEvents,
SubEventIterator  eSubEvents 
)
virtual

Reimplemented from PileUpToolBase.

Definition at line 273 of file SiSmearedDigitizationTool.cxx.

276 {
277  ATH_MSG_DEBUG( "--- SiSmearedDigitizationTool: in pixel processBunchXing() ---" );
278  //decide if this event will be processed depending on HardScatterSplittingMode & bunchXing
279  if (m_HardScatterSplittingMode == 2 && !m_HardScatterSplittingSkipper ) { m_HardScatterSplittingSkipper = true; return StatusCode::SUCCESS; }
280  if (m_HardScatterSplittingMode == 1 && m_HardScatterSplittingSkipper ) { return StatusCode::SUCCESS; }
282 
283  using TimedHitCollList = PileUpMergeSvc::TimedList<SiHitCollection>::type;
284  TimedHitCollList hitCollList;
285 
286  if (!(m_mergeSvc->retrieveSubSetEvtData(m_inputObjectName, hitCollList, bunchXing,
287  bSubEvents, eSubEvents).isSuccess()) &&
288  hitCollList.empty()) {
289  ATH_MSG_ERROR("Could not fill TimedHitCollList");
290  return StatusCode::FAILURE;
291  } else {
292  ATH_MSG_VERBOSE(hitCollList.size() << " SiHitCollections with key " <<
293  m_inputObjectName << " found");
294  }
295 
296  TimedHitCollList::iterator iColl(hitCollList.begin());
297  TimedHitCollList::iterator endColl(hitCollList.end());
298 
299  for( ; iColl != endColl; ++iColl) {
300  SiHitCollection *siHitColl = new SiHitCollection(*iColl->second);
301  PileUpTimeEventIndex timeIndex(iColl->first);
302  ATH_MSG_DEBUG("SiHitCollection found with " << siHitColl->size() <<
303  " hits");
304  ATH_MSG_VERBOSE("time index info. time: " << timeIndex.time()
305  << " index: " << timeIndex.index()
306  << " type: " << timeIndex.type());
307  m_siHitCollList.push_back(siHitColl);
308  }
309 
310  return StatusCode::SUCCESS;
311 }

◆ resetFilter()

virtual void PileUpToolBase::resetFilter ( )
inlineoverridevirtualinherited

dummy implementation of filter reset

Reimplemented in MergeTruthJetsTool.

Definition at line 51 of file PileUpToolBase.h.

51 { m_filterPassed=true; }

◆ retrieveTruth()

StatusCode SiSmearedDigitizationTool::retrieveTruth ( )

Definition at line 439 of file SiSmearedDigitizationTool.cxx.

439  {
440 
441 
442  if(m_SmearPixel){ // Smear Pixel
443 
445 
446  if ((evtStore()->contains<PRD_MultiTruthCollection>(m_prdTruthNamePixel))){
447  if((evtStore()->retrieve(m_pixelPrdTruth, m_prdTruthNamePixel)).isFailure()){
448  ATH_MSG_FATAL("Could not retrieve collection " << m_prdTruthNamePixel);
449  return StatusCode::FAILURE;
450  }
451  }else{
452  if((evtStore()->record(m_pixelPrdTruth, m_prdTruthNamePixel)).isFailure()){
453  ATH_MSG_FATAL("Could not record collection " << m_prdTruthNamePixel);
454  return StatusCode::FAILURE;
455  }
456  }
457  }else{ // Smear SCT
459 
460  if ((evtStore()->contains<PRD_MultiTruthCollection>(m_prdTruthNameSCT))){
461  if((evtStore()->retrieve(m_SCTPrdTruth, m_prdTruthNameSCT)).isFailure()){
462  ATH_MSG_FATAL("Could not retrieve collection " << m_prdTruthNameSCT);
463  return StatusCode::FAILURE;
464  }
465  }else{
466  if((evtStore()->record(m_SCTPrdTruth, m_prdTruthNameSCT)).isFailure()){
467  ATH_MSG_FATAL("Could not record collection " << m_prdTruthNameSCT);
468  return StatusCode::FAILURE;
469  }
470  }
471  }
472 
473 
474  return StatusCode::SUCCESS;
475 
476 }

◆ toProcess()

virtual bool PileUpToolBase::toProcess ( int  bunchXing) const
inlineoverridevirtualinherited

the method this base class helps implementing

Reimplemented in MergeHijingParsTool, and MergeTrackRecordCollTool.

Definition at line 32 of file PileUpToolBase.h.

32  {
33  //closed interval [m_firstXing,m_lastXing]
34  return !((m_firstXing > bunchXing) || (bunchXing > m_lastXing));
35  }

Member Data Documentation

◆ m_checkSmear

bool SiSmearedDigitizationTool::m_checkSmear
private

Definition at line 155 of file SiSmearedDigitizationTool.h.

◆ m_currentTree

TTree* SiSmearedDigitizationTool::m_currentTree
private

the tree to store information from pixel and SCT (before and after smearing)

Definition at line 159 of file SiSmearedDigitizationTool.h.

◆ m_emulateAtlas

bool SiSmearedDigitizationTool::m_emulateAtlas
private

Definition at line 150 of file SiSmearedDigitizationTool.h.

◆ m_Err_x_pixel

double SiSmearedDigitizationTool::m_Err_x_pixel
private

Definition at line 188 of file SiSmearedDigitizationTool.h.

◆ m_Err_x_SCT

double SiSmearedDigitizationTool::m_Err_x_SCT
private

Definition at line 190 of file SiSmearedDigitizationTool.h.

◆ m_Err_y_pixel

double SiSmearedDigitizationTool::m_Err_y_pixel
private

Definition at line 189 of file SiSmearedDigitizationTool.h.

◆ m_Err_y_SCT

double SiSmearedDigitizationTool::m_Err_y_SCT
private

Definition at line 191 of file SiSmearedDigitizationTool.h.

◆ m_filterPassed

bool PileUpToolBase::m_filterPassed {true}
protectedinherited

Definition at line 60 of file PileUpToolBase.h.

◆ m_firstXing

Gaudi::Property<int> PileUpToolBase::m_firstXing
protectedinherited
Initial value:
{this, "FirstXing", -999,
"First bunch-crossing in which det is live"}

Definition at line 54 of file PileUpToolBase.h.

◆ m_HardScatterSplittingMode

int SiSmearedDigitizationTool::m_HardScatterSplittingMode
private

Process all SiHit or just those from signal or background events.

Definition at line 130 of file SiSmearedDigitizationTool.h.

◆ m_HardScatterSplittingSkipper

bool SiSmearedDigitizationTool::m_HardScatterSplittingSkipper
private

Definition at line 131 of file SiSmearedDigitizationTool.h.

◆ m_inputObjectName

std::string SiSmearedDigitizationTool::m_inputObjectName
private

Definition at line 140 of file SiSmearedDigitizationTool.h.

◆ m_lastXing

Gaudi::Property<int> PileUpToolBase::m_lastXing
protectedinherited
Initial value:
{this, "LastXing", 999,
"Last bunch-crossing in which det is live"}

Definition at line 56 of file PileUpToolBase.h.

◆ m_merge

bool SiSmearedDigitizationTool::m_merge
private

Definition at line 120 of file SiSmearedDigitizationTool.h.

◆ m_mergeSvc

ServiceHandle<PileUpMergeSvc> SiSmearedDigitizationTool::m_mergeSvc
private

PileUp Merge service.

Definition at line 129 of file SiSmearedDigitizationTool.h.

◆ m_nSigma

double SiSmearedDigitizationTool::m_nSigma
private

Definition at line 121 of file SiSmearedDigitizationTool.h.

◆ m_outputFile

TFile* SiSmearedDigitizationTool::m_outputFile
private

the root file

Definition at line 158 of file SiSmearedDigitizationTool.h.

◆ m_pitch_X

float SiSmearedDigitizationTool::m_pitch_X
private

Definition at line 117 of file SiSmearedDigitizationTool.h.

◆ m_pitch_Y

float SiSmearedDigitizationTool::m_pitch_Y
private

Definition at line 118 of file SiSmearedDigitizationTool.h.

◆ m_pixel_ID

const PixelID* SiSmearedDigitizationTool::m_pixel_ID
private

Handle to the ID helper.

Definition at line 112 of file SiSmearedDigitizationTool.h.

◆ m_pixel_SiClustersName

std::string SiSmearedDigitizationTool::m_pixel_SiClustersName
private

Definition at line 152 of file SiSmearedDigitizationTool.h.

◆ m_pixelClusterContainer

InDet::PixelClusterContainer* SiSmearedDigitizationTool::m_pixelClusterContainer
private

the PixelClusterContainer

Definition at line 125 of file SiSmearedDigitizationTool.h.

◆ m_pixelClusterMap

Pixel_detElement_RIO_map* SiSmearedDigitizationTool::m_pixelClusterMap {}
private

Definition at line 144 of file SiSmearedDigitizationTool.h.

◆ m_pixelDetEleCollKey

SG::ReadCondHandleKey<InDetDD::SiDetectorElementCollection> SiSmearedDigitizationTool::m_pixelDetEleCollKey {this, "PixelDetEleCollKey", "PixelDetectorElementCollection", "Key of SiDetectorElementCollection for Pixel"}
private

Definition at line 109 of file SiSmearedDigitizationTool.h.

◆ m_pixelPrdTruth

PRD_MultiTruthCollection* SiSmearedDigitizationTool::m_pixelPrdTruth {}
private

Definition at line 133 of file SiSmearedDigitizationTool.h.

◆ m_prdTruthNamePixel

std::string SiSmearedDigitizationTool::m_prdTruthNamePixel
private

Definition at line 134 of file SiSmearedDigitizationTool.h.

◆ m_prdTruthNameSCT

std::string SiSmearedDigitizationTool::m_prdTruthNameSCT
private

Definition at line 137 of file SiSmearedDigitizationTool.h.

◆ m_randomEngineName

std::string SiSmearedDigitizationTool::m_randomEngineName
private

Name of the random number stream.

Definition at line 115 of file SiSmearedDigitizationTool.h.

◆ m_rndmSvc

ServiceHandle<IAthRNGSvc> SiSmearedDigitizationTool::m_rndmSvc {this, "RndmSvc", "AthRNGSvc", ""}
private

Random number service.

Definition at line 107 of file SiSmearedDigitizationTool.h.

◆ m_sct_ID

const SCT_ID* SiSmearedDigitizationTool::m_sct_ID
private

Handle to the ID helper.

Definition at line 113 of file SiSmearedDigitizationTool.h.

◆ m_Sct_SiClustersName

std::string SiSmearedDigitizationTool::m_Sct_SiClustersName
private

Definition at line 153 of file SiSmearedDigitizationTool.h.

◆ m_sctClusterContainer

InDet::SCT_ClusterContainer* SiSmearedDigitizationTool::m_sctClusterContainer
private

the SCT_ClusterContainer

Definition at line 127 of file SiSmearedDigitizationTool.h.

◆ m_sctClusterMap

SCT_detElement_RIO_map* SiSmearedDigitizationTool::m_sctClusterMap {}
private

Definition at line 146 of file SiSmearedDigitizationTool.h.

◆ m_SCTDetEleCollKey

SG::ReadCondHandleKey<InDetDD::SiDetectorElementCollection> SiSmearedDigitizationTool::m_SCTDetEleCollKey {this, "SCTDetEleCollKey", "SCT_DetectorElementCollection", "Key of SiDetectorElementCollection for SCT"}
private

Definition at line 110 of file SiSmearedDigitizationTool.h.

◆ m_SCTPrdTruth

PRD_MultiTruthCollection* SiSmearedDigitizationTool::m_SCTPrdTruth {}
private

Definition at line 136 of file SiSmearedDigitizationTool.h.

◆ m_siHitCollList

std::vector<SiHitCollection*> SiSmearedDigitizationTool::m_siHitCollList
private

name of the sub event hit collections.

Definition at line 142 of file SiSmearedDigitizationTool.h.

◆ m_simHitColl

SiHitCollection* SiSmearedDigitizationTool::m_simHitColl {}
private

Definition at line 139 of file SiSmearedDigitizationTool.h.

◆ m_SmearPixel

bool SiSmearedDigitizationTool::m_SmearPixel
private

Definition at line 148 of file SiSmearedDigitizationTool.h.

◆ m_thistSvc

ServiceHandle<ITHistSvc> SiSmearedDigitizationTool::m_thistSvc {this, "THistSvc", "THistSvc"}
private

Definition at line 157 of file SiSmearedDigitizationTool.h.

◆ m_useDiscSurface

bool SiSmearedDigitizationTool::m_useDiscSurface
private

Definition at line 123 of file SiSmearedDigitizationTool.h.

◆ m_vetoPileUpTruthLinks

Gaudi::Property<int> PileUpToolBase::m_vetoPileUpTruthLinks
protectedinherited
Initial value:
{this, "VetoPileUpTruthLinks", true,
"Ignore links to suppressed pile-up truth"}

Definition at line 58 of file PileUpToolBase.h.

◆ m_x_entry_pixel

double SiSmearedDigitizationTool::m_x_entry_pixel
private

Definition at line 166 of file SiSmearedDigitizationTool.h.

◆ m_x_entry_SCT

double SiSmearedDigitizationTool::m_x_entry_SCT
private

Definition at line 177 of file SiSmearedDigitizationTool.h.

◆ m_x_exit_pixel

double SiSmearedDigitizationTool::m_x_exit_pixel
private

Definition at line 163 of file SiSmearedDigitizationTool.h.

◆ m_x_exit_SCT

double SiSmearedDigitizationTool::m_x_exit_SCT
private

Definition at line 174 of file SiSmearedDigitizationTool.h.

◆ m_x_pixel

double SiSmearedDigitizationTool::m_x_pixel
private

Definition at line 161 of file SiSmearedDigitizationTool.h.

◆ m_x_pixel_global

double SiSmearedDigitizationTool::m_x_pixel_global
private

Definition at line 169 of file SiSmearedDigitizationTool.h.

◆ m_x_pixel_smeared

double SiSmearedDigitizationTool::m_x_pixel_smeared
private

Definition at line 184 of file SiSmearedDigitizationTool.h.

◆ m_x_SCT

double SiSmearedDigitizationTool::m_x_SCT
private

Definition at line 173 of file SiSmearedDigitizationTool.h.

◆ m_x_SCT_global

double SiSmearedDigitizationTool::m_x_SCT_global
private

Definition at line 180 of file SiSmearedDigitizationTool.h.

◆ m_x_SCT_smeared

double SiSmearedDigitizationTool::m_x_SCT_smeared
private

Definition at line 186 of file SiSmearedDigitizationTool.h.

◆ m_y_entry_pixel

double SiSmearedDigitizationTool::m_y_entry_pixel
private

Definition at line 167 of file SiSmearedDigitizationTool.h.

◆ m_y_entry_SCT

double SiSmearedDigitizationTool::m_y_entry_SCT
private

Definition at line 178 of file SiSmearedDigitizationTool.h.

◆ m_y_exit_pixel

double SiSmearedDigitizationTool::m_y_exit_pixel
private

Definition at line 164 of file SiSmearedDigitizationTool.h.

◆ m_y_exit_SCT

double SiSmearedDigitizationTool::m_y_exit_SCT
private

Definition at line 175 of file SiSmearedDigitizationTool.h.

◆ m_y_pixel

double SiSmearedDigitizationTool::m_y_pixel
private

Definition at line 162 of file SiSmearedDigitizationTool.h.

◆ m_y_pixel_global

double SiSmearedDigitizationTool::m_y_pixel_global
private

Definition at line 170 of file SiSmearedDigitizationTool.h.

◆ m_y_pixel_smeared

double SiSmearedDigitizationTool::m_y_pixel_smeared
private

Definition at line 185 of file SiSmearedDigitizationTool.h.

◆ m_y_SCT_global

double SiSmearedDigitizationTool::m_y_SCT_global
private

Definition at line 181 of file SiSmearedDigitizationTool.h.

◆ m_z_entry_pixel

double SiSmearedDigitizationTool::m_z_entry_pixel
private

Definition at line 168 of file SiSmearedDigitizationTool.h.

◆ m_z_entry_SCT

double SiSmearedDigitizationTool::m_z_entry_SCT
private

Definition at line 179 of file SiSmearedDigitizationTool.h.

◆ m_z_exit_pixel

double SiSmearedDigitizationTool::m_z_exit_pixel
private

Definition at line 165 of file SiSmearedDigitizationTool.h.

◆ m_z_exit_SCT

double SiSmearedDigitizationTool::m_z_exit_SCT
private

Definition at line 176 of file SiSmearedDigitizationTool.h.

◆ m_z_pixel_global

double SiSmearedDigitizationTool::m_z_pixel_global
private

Definition at line 171 of file SiSmearedDigitizationTool.h.

◆ m_z_SCT_global

double SiSmearedDigitizationTool::m_z_SCT_global
private

Definition at line 182 of file SiSmearedDigitizationTool.h.


The documentation for this class was generated from the following files:
python.PyKernel.retrieve
def retrieve(aClass, aKey=None)
Definition: PyKernel.py:110
xAOD::iterator
JetConstituentVector::iterator iterator
Definition: JetConstituentVector.cxx:68
InDetDD::SolidStateDetectorElementBase::identifierOfPosition
Identifier identifierOfPosition(const Amg::Vector2D &localPos) const
Full identifier of the cell for a given position: assumes a raw local position (no Lorentz shift)
Definition: SolidStateDetectorElementBase.cxx:217
SiSmearedDigitizationTool::m_y_entry_pixel
double m_y_entry_pixel
Definition: SiSmearedDigitizationTool.h:167
AllowedVariables::e
e
Definition: AsgElectronSelectorTool.cxx:37
ATHRNG::RNGWrapper::setSeed
void setSeed(const std::string &algName, const EventContext &ctx)
Set the random seed using a string (e.g.
Definition: RNGWrapper.h:169
SiSmearedDigitizationTool::m_z_exit_SCT
double m_z_exit_SCT
Definition: SiSmearedDigitizationTool.h:176
SiSmearedDigitizationTool::m_pitch_Y
float m_pitch_Y
Definition: SiSmearedDigitizationTool.h:118
SiSmearedDigitizationTool::m_x_pixel_global
double m_x_pixel_global
Definition: SiSmearedDigitizationTool.h:169
SiSmearedDigitizationTool::m_HardScatterSplittingSkipper
bool m_HardScatterSplittingSkipper
Definition: SiSmearedDigitizationTool.h:131
ATH_MSG_FATAL
#define ATH_MSG_FATAL(x)
Definition: AthMsgStreamMacros.h:34
pdg_comparison.sigma
sigma
Definition: pdg_comparison.py:324
Amg::MatrixX
Eigen::Matrix< double, Eigen::Dynamic, Eigen::Dynamic > MatrixX
Dynamic Matrix - dynamic allocation.
Definition: EventPrimitives.h:27
PixelCluster
Definition: Trigger/TrigAccel/TrigCudaFitter/src/PixelCluster.h:8
SiSmearedDigitizationTool::m_pitch_X
float m_pitch_X
Definition: SiSmearedDigitizationTool.h:117
InDetDD::SiDetectorElementCollection
Definition: SiDetectorElementCollection.h:30
SG::ReadCondHandle
Definition: ReadCondHandle.h:44
PileUpToolBase::m_filterPassed
bool m_filterPassed
Definition: PileUpToolBase.h:60
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
Trk::locX
@ locX
Definition: ParamDefs.h:37
Trk::locY
@ locY
local cartesian
Definition: ParamDefs.h:38
SiSmearedDigitizationTool::retrieveTruth
StatusCode retrieveTruth()
Definition: SiSmearedDigitizationTool.cxx:439
SiSmearedDigitizationTool::m_x_exit_SCT
double m_x_exit_SCT
Definition: SiSmearedDigitizationTool.h:174
SiSmearedDigitizationTool::m_checkSmear
bool m_checkSmear
Definition: SiSmearedDigitizationTool.h:155
SiSmearedDigitizationTool::m_SmearPixel
bool m_SmearPixel
Definition: SiSmearedDigitizationTool.h:148
SiSmearedDigitizationTool::m_z_exit_pixel
double m_z_exit_pixel
Definition: SiSmearedDigitizationTool.h:165
SiSmearedDigitizationTool::FillTruthMap
StatusCode FillTruthMap(PRD_MultiTruthCollection *, CLUSTER *, const TimedHitPtr< SiHit > &)
Definition: SiSmearedDigitizationTool.cxx:479
Amg::Vector2D
Eigen::Matrix< double, 2, 1 > Vector2D
Definition: GeoPrimitives.h:48
SiSmearedDigitizationTool::m_pixelClusterContainer
InDet::PixelClusterContainer * m_pixelClusterContainer
the PixelClusterContainer
Definition: SiSmearedDigitizationTool.h:125
SiSmearedDigitizationTool::m_merge
bool m_merge
Definition: SiSmearedDigitizationTool.h:120
InDetDD::SCT_ModuleSideDesign
Definition: SCT_ModuleSideDesign.h:40
InDet::SiWidth::widthPhiRZ
const Amg::Vector2D & widthPhiRZ() const
Definition: SiWidth.h:121
InDet
Primary Vertex Finder.
Definition: VP1ErrorUtils.h:36
PRD_MultiTruthCollection
A PRD is mapped onto all contributing particles.
Definition: PRD_MultiTruthCollection.h:24
SiSmearedDigitizationTool::m_pixelClusterMap
Pixel_detElement_RIO_map * m_pixelClusterMap
Definition: SiSmearedDigitizationTool.h:144
mat
GeoMaterial * mat
Definition: LArDetectorConstructionTBEC.cxx:55
SiSmearedDigitizationTool::m_thistSvc
ServiceHandle< ITHistSvc > m_thistSvc
Definition: SiSmearedDigitizationTool.h:157
CSV_InDetExporter.new
new
Definition: CSV_InDetExporter.py:145
AtlasHitsVector
Definition: AtlasHitsVector.h:33
InDetDD::SiCellId::isValid
bool isValid() const
Test if its in a valid state.
Definition: SiCellId.h:136
SiSmearedDigitizationTool::calculateSigma
double calculateSigma(CLUSTER *clusterA, CLUSTER *clusterB)
Definition: SiSmearedDigitizationTool.cxx:521
SiSmearedDigitizationTool::m_mergeSvc
ServiceHandle< PileUpMergeSvc > m_mergeSvc
PileUp Merge service.
Definition: SiSmearedDigitizationTool.h:129
SiSmearedDigitizationTool::m_z_SCT_global
double m_z_SCT_global
Definition: SiSmearedDigitizationTool.h:182
InDet::SCT_ClusterContainer
Trk::PrepRawDataContainer< SCT_ClusterCollection > SCT_ClusterContainer
Definition: SCT_ClusterContainer.h:27
SiSmearedDigitizationTool::m_pixelDetEleCollKey
SG::ReadCondHandleKey< InDetDD::SiDetectorElementCollection > m_pixelDetEleCollKey
Definition: SiSmearedDigitizationTool.h:109
SiSmearedDigitizationTool::calculateDistance
double calculateDistance(CLUSTER *clusterA, CLUSTER *clusterB)
Definition: SiSmearedDigitizationTool.cxx:503
SiSmearedDigitizationTool::calculateNewCluster
ClusterInfo calculateNewCluster(CLUSTER *clusterA, CLUSTER *clusterB)
Definition: SiSmearedDigitizationTool.cxx:538
SiSmearedDigitizationTool::m_y_SCT_global
double m_y_SCT_global
Definition: SiSmearedDigitizationTool.h:181
TimedHitPtr< SiHit >
InDet::PixelClusterContainer
Trk::PrepRawDataContainer< PixelClusterCollection > PixelClusterContainer
Definition: InnerDetector/InDetRecEvent/InDetPrepRawData/InDetPrepRawData/PixelClusterContainer.h:28
SiSmearedDigitizationTool::m_Err_y_SCT
double m_Err_y_SCT
Definition: SiSmearedDigitizationTool.h:191
ATH_MSG_VERBOSE
#define ATH_MSG_VERBOSE(x)
Definition: AthMsgStreamMacros.h:28
PileUpToolBase::m_lastXing
Gaudi::Property< int > m_lastXing
Definition: PileUpToolBase.h:56
beamspotman.sigmaX
sigmaX
Definition: beamspotman.py:1625
HepMC::ignoreTruthLink
bool ignoreTruthLink(const T &p, bool vetoPileUp)
Helper function for SDO creation in PileUpTools.
Definition: MagicNumbers.h:345
beamspotman.sigmaY
sigmaY
Definition: beamspotman.py:1625
x
#define x
TimedHitCollection::nextDetectorElement
bool nextDetectorElement(const_iterator &b, const_iterator &e)
sets an iterator range with the hits of current detector element returns a bool when done
intersection
std::vector< std::string > intersection(std::vector< std::string > &v1, std::vector< std::string > &v2)
Definition: compareFlatTrees.cxx:25
SiSmearedDigitizationTool::mergeClusters
StatusCode mergeClusters(Pixel_detElement_RIO_map *cluster_map)
Definition: SiSmearedDigitizationTool.cxx:583
AmgSymMatrix
#define AmgSymMatrix(dim)
Definition: EventPrimitives.h:50
PixelID::wafer_id
Identifier wafer_id(int barrel_ec, int layer_disk, int phi_module, int eta_module) const
For a single crystal.
Definition: PixelID.h:364
SiSmearedDigitizationTool::m_sctClusterContainer
InDet::SCT_ClusterContainer * m_sctClusterContainer
the SCT_ClusterContainer
Definition: SiSmearedDigitizationTool.h:127
PileUpToolBase::m_vetoPileUpTruthLinks
Gaudi::Property< int > m_vetoPileUpTruthLinks
Definition: PileUpToolBase.h:58
PileUpMergeSvc::TimedList::type
std::list< value_t > type
type of the collection of timed data object
Definition: PileUpMergeSvc.h:75
TRT::Hit::side
@ side
Definition: HitInfo.h:83
SiSmearedDigitizationTool::m_Err_y_pixel
double m_Err_y_pixel
Definition: SiSmearedDigitizationTool.h:189
SiSmearedDigitizationTool::m_x_pixel
double m_x_pixel
Definition: SiSmearedDigitizationTool.h:161
InDetDD::SiDetectorElement::cellIdFromIdentifier
virtual SiCellId cellIdFromIdentifier(const Identifier &identifier) const override final
SiCellId from Identifier.
Definition: SiDetectorElement.cxx:120
PixelID::wafer_hash
IdentifierHash wafer_hash(Identifier wafer_id) const
wafer hash from id
Definition: PixelID.h:387
SiSmearedDigitizationTool::m_randomEngineName
std::string m_randomEngineName
Name of the random number stream.
Definition: SiSmearedDigitizationTool.h:115
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
SiSmearedDigitizationTool::m_x_SCT_smeared
double m_x_SCT_smeared
Definition: SiSmearedDigitizationTool.h:186
fillPileUpNoiseLumi.next
next
Definition: fillPileUpNoiseLumi.py:52
SiSmearedDigitizationTool::m_useDiscSurface
bool m_useDiscSurface
Definition: SiSmearedDigitizationTool.h:123
parseMapping.v0
def v0
Definition: parseMapping.py:149
lumiFormat.i
int i
Definition: lumiFormat.py:85
SiSmearedDigitizationTool::m_sctClusterMap
SCT_detElement_RIO_map * m_sctClusterMap
Definition: SiSmearedDigitizationTool.h:146
SiSmearedDigitizationTool::m_z_entry_SCT
double m_z_entry_SCT
Definition: SiSmearedDigitizationTool.h:179
z
#define z
SiSmearedDigitizationTool::digitize
StatusCode digitize(const EventContext &ctx, TimedHitCollection< SiHit > &thpcsi)
Definition: SiSmearedDigitizationTool.cxx:717
Trk::PrepRawData::setHashAndIndex
void setHashAndIndex(unsigned short collHash, unsigned short objIndex)
TEMP for testing: might make some classes friends later ...
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
plotBeamSpotVxVal.range
range
Definition: plotBeamSpotVxVal.py:195
test_pyathena.parent
parent
Definition: test_pyathena.py:15
PileUpToolBase::PileUpToolBase
PileUpToolBase(const std::string &type, const std::string &name, const IInterface *parent)
Definition: PileUpToolBase.cxx:7
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
SiSmearedDigitizationTool::m_y_pixel_smeared
double m_y_pixel_smeared
Definition: SiSmearedDigitizationTool.h:185
SCT_ID::wafer_hash
IdentifierHash wafer_hash(const Identifier &wafer_id) const
wafer hash from id - optimized
Definition: SCT_ID.h:492
SiSmearedDigitizationTool::SCT_detElement_RIO_map
std::multimap< IdentifierHash, InDet::SCT_Cluster * > SCT_detElement_RIO_map
Definition: SiSmearedDigitizationTool.h:83
xAOD::double
double
Definition: CompositeParticle_v1.cxx:159
InDet::SCT_Cluster
Definition: InnerDetector/InDetRecEvent/InDetPrepRawData/InDetPrepRawData/SCT_Cluster.h:34
SiSmearedDigitizationTool::m_z_entry_pixel
double m_z_entry_pixel
Definition: SiSmearedDigitizationTool.h:168
SiSmearedDigitizationTool::m_x_entry_SCT
double m_x_entry_SCT
Definition: SiSmearedDigitizationTool.h:177
SiSmearedDigitizationTool::m_y_pixel_global
double m_y_pixel_global
Definition: SiSmearedDigitizationTool.h:170
SiSmearedDigitizationTool::m_y_exit_pixel
double m_y_exit_pixel
Definition: SiSmearedDigitizationTool.h:164
SiSmearedDigitizationTool::m_SCTDetEleCollKey
SG::ReadCondHandleKey< InDetDD::SiDetectorElementCollection > m_SCTDetEleCollKey
Definition: SiSmearedDigitizationTool.h:110
SiSmearedDigitizationTool::createAndStoreRIOs
StatusCode createAndStoreRIOs(const EventContext &ctx)
Definition: SiSmearedDigitizationTool.cxx:1125
SiSmearedDigitizationTool::m_prdTruthNamePixel
std::string m_prdTruthNamePixel
Definition: SiSmearedDigitizationTool.h:134
SiSmearedDigitizationTool::m_pixelPrdTruth
PRD_MultiTruthCollection * m_pixelPrdTruth
Definition: SiSmearedDigitizationTool.h:133
python.PyKernel.detStore
detStore
Definition: PyKernel.py:41
Trk
Ensure that the ATLAS eigen extensions are properly loaded.
Definition: FakeTrackBuilder.h:9
SiSmearedDigitizationTool::m_simHitColl
SiHitCollection * m_simHitColl
Definition: SiSmearedDigitizationTool.h:139
Amg
Definition of ATLAS Math & Geometry primitives (Amg)
Definition: AmgStringHelpers.h:19
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:228
SCT_ID::wafer_hash_max
size_type wafer_hash_max(void) const
Definition: SCT_ID.cxx:645
ATHRNG::RNGWrapper
A wrapper class for event-slot-local random engines.
Definition: RNGWrapper.h:56
Amg::error
double error(const Amg::MatrixX &mat, int index)
return diagonal error of the matrix caller should ensure the matrix is symmetric and the index is in ...
Definition: EventPrimitivesHelpers.h:40
Trk::PrepRawData::localPosition
const Amg::Vector2D & localPosition() const
return the local position reference
SiSmearedDigitizationTool::m_Err_x_pixel
double m_Err_x_pixel
Definition: SiSmearedDigitizationTool.h:188
InDetDD::SiDetectorElement
Definition: SiDetectorElement.h:109
SG::CondHandleKey::initialize
StatusCode initialize(bool used=true)
SiSmearedDigitizationTool::m_HardScatterSplittingMode
int m_HardScatterSplittingMode
Process all SiHit or just those from signal or background events.
Definition: SiSmearedDigitizationTool.h:130
PixelID::wafer_hash_max
size_type wafer_hash_max(void) const
Definition: PixelID.cxx:831
SiSmearedDigitizationTool::m_siHitCollList
std::vector< SiHitCollection * > m_siHitCollList
name of the sub event hit collections.
Definition: SiSmearedDigitizationTool.h:142
InDet::SiClusterContainer
Trk::PrepRawDataContainer< SiClusterCollection > SiClusterContainer
Definition: SiClusterContainer.h:26
SiSmearedDigitizationTool::m_emulateAtlas
bool m_emulateAtlas
Definition: SiSmearedDigitizationTool.h:150
Amg::Vector3D
Eigen::Matrix< double, 3, 1 > Vector3D
Definition: GeoPrimitives.h:47
InDetDD::DetectorShape
DetectorShape
Definition: DetectorDesign.h:41
SiSmearedDigitizationTool::m_z_pixel_global
double m_z_pixel_global
Definition: SiSmearedDigitizationTool.h:171
ATHRNG::RNGWrapper::getEngine
CLHEP::HepRandomEngine * getEngine(const EventContext &ctx) const
Retrieve the random engine corresponding to the provided EventContext.
Definition: RNGWrapper.h:134
SiSmearedDigitizationTool::m_currentTree
TTree * m_currentTree
the tree to store information from pixel and SCT (before and after smearing)
Definition: SiSmearedDigitizationTool.h:159
SiSmearedDigitizationTool::m_inputObjectName
std::string m_inputObjectName
Definition: SiSmearedDigitizationTool.h:140
InDetDD::SiCellId
Definition: SiCellId.h:29
InDet::PixelCluster
Definition: InnerDetector/InDetRecEvent/InDetPrepRawData/InDetPrepRawData/PixelCluster.h:49
Trk::pixelCluster
@ pixelCluster
Definition: MeasurementType.h:22
SiSmearedDigitizationTool::m_pixel_SiClustersName
std::string m_pixel_SiClustersName
Definition: SiSmearedDigitizationTool.h:152
SiHitCollection
AtlasHitsVector< SiHit > SiHitCollection
Definition: SiHitCollection.h:14
InDet::SiCluster::globalPosition
const Amg::Vector3D & globalPosition() const
return global position reference
y
#define y
TRT::Hit::phiModule
@ phiModule
Definition: HitInfo.h:80
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
InDetDD
Message Stream Member.
Definition: FakeTrackBuilder.h:8
SiSmearedDigitizationTool::m_rndmSvc
ServiceHandle< IAthRNGSvc > m_rndmSvc
Random number service.
Definition: SiSmearedDigitizationTool.h:107
SiSmearedDigitizationTool::m_y_exit_SCT
double m_y_exit_SCT
Definition: SiSmearedDigitizationTool.h:175
python.CaloScaleNoiseConfig.type
type
Definition: CaloScaleNoiseConfig.py:78
PileUpToolBase::m_firstXing
Gaudi::Property< int > m_firstXing
Definition: PileUpToolBase.h:54
SiSmearedDigitizationTool::m_Err_x_SCT
double m_Err_x_SCT
Definition: SiSmearedDigitizationTool.h:190
SiSmearedDigitizationTool::m_x_SCT
double m_x_SCT
Definition: SiSmearedDigitizationTool.h:173
if
if(febId1==febId2)
Definition: LArRodBlockPhysicsV0.cxx:567
InDet::SiWidth
Definition: SiWidth.h:25
InDet::SiWidth::colRow
const Amg::Vector2D & colRow() const
Definition: SiWidth.h:115
AtlasHitsVector::size
size_type size() const
Definition: AtlasHitsVector.h:143
SiSmearedDigitizationTool::m_Sct_SiClustersName
std::string m_Sct_SiClustersName
Definition: SiSmearedDigitizationTool.h:153
SiSmearedDigitizationTool::m_x_SCT_global
double m_x_SCT_global
Definition: SiSmearedDigitizationTool.h:180
SiSmearedDigitizationTool::m_x_entry_pixel
double m_x_entry_pixel
Definition: SiSmearedDigitizationTool.h:166
SiSmearedDigitizationTool::m_pixel_ID
const PixelID * m_pixel_ID
Handle to the ID helper.
Definition: SiSmearedDigitizationTool.h:112
SiSmearedDigitizationTool::m_y_pixel
double m_y_pixel
Definition: SiSmearedDigitizationTool.h:162
SiSmearedDigitizationTool::m_x_pixel_smeared
double m_x_pixel_smeared
Definition: SiSmearedDigitizationTool.h:184
SiSmearedDigitizationTool::m_sct_ID
const SCT_ID * m_sct_ID
Handle to the ID helper.
Definition: SiSmearedDigitizationTool.h:113
SCT_ID::wafer_id
Identifier wafer_id(int barrel_ec, int layer_disk, int phi_module, int eta_module, int side) const
For a single side of module.
Definition: SCT_ID.h:464
python.IoTestsLib.w
def w
Definition: IoTestsLib.py:200
InDet::SCT_ClusterCollection
Trk::PrepRawDataCollection< SCT_Cluster > SCT_ClusterCollection
Definition: SCT_ClusterCollection.h:26
pow
constexpr int pow(int base, int exp) noexcept
Definition: ap_fixedTest.cxx:15
PileUpTimeEventIndex
a struct encapsulating the identifier of a pile-up event
Definition: PileUpTimeEventIndex.h:12
SiSmearedDigitizationTool::mergeEvent
StatusCode mergeEvent(const EventContext &ctx)
Definition: SiSmearedDigitizationTool.cxx:494
SiSmearedDigitizationTool::m_nSigma
double m_nSigma
Definition: SiSmearedDigitizationTool.h:121
InDetDD::SCT_ModuleSideDesign::endsOfStrip
virtual std::pair< SiLocalPosition, SiLocalPosition > endsOfStrip(const SiLocalPosition &position) const override=0
give the ends of strips
length
double length(const pvec &v)
Definition: FPGATrackSimLLPDoubletHoughTransformTool.cxx:26
SiSmearedDigitizationTool::Pixel_detElement_RIO_map
std::multimap< IdentifierHash, InDet::PixelCluster * > Pixel_detElement_RIO_map
Definition: SiSmearedDigitizationTool.h:82
SiSmearedDigitizationTool::m_y_entry_SCT
double m_y_entry_SCT
Definition: SiSmearedDigitizationTool.h:178
InDetDD::SiDetectorElementCollection::getDetectorElement
const SiDetectorElement * getDetectorElement(const IdentifierHash &hash) const
Definition: SiDetectorElementCollection.cxx:15
TimedHitCollection
Definition: TimedHitCollection.h:15
InDetDD::SolidStateDetectorElementBase::identify
virtual Identifier identify() const override final
identifier of this detector element (inline)
SiSmearedDigitizationTool::m_outputFile
TFile * m_outputFile
the root file
Definition: SiSmearedDigitizationTool.h:158
InDet::PixelClusterCollection
Trk::PrepRawDataCollection< PixelCluster > PixelClusterCollection
Definition: PixelClusterCollection.h:26
SiSmearedDigitizationTool::m_SCTPrdTruth
PRD_MultiTruthCollection * m_SCTPrdTruth
Definition: SiSmearedDigitizationTool.h:136
InDetDD::Trapezoid
@ Trapezoid
Definition: DetectorDesign.h:42
ClusterInfo
std::tuple< Amg::Vector2D, InDet::SiWidth, Amg::MatrixX > ClusterInfo
Definition: SiSmearedDigitizationTool.h:44
SiSmearedDigitizationTool::m_prdTruthNameSCT
std::string m_prdTruthNameSCT
Definition: SiSmearedDigitizationTool.h:137
Identifier
Definition: IdentifierFieldParser.cxx:14
SiSmearedDigitizationTool::m_x_exit_pixel
double m_x_exit_pixel
Definition: SiSmearedDigitizationTool.h:163