#include <PixelFastDigitizationTool.h>

Inheritance diagram for PixelFastDigitizationTool:

Collaboration diagram for PixelFastDigitizationTool:

Public Member Functions
	PixelFastDigitizationTool (const std::string &type, const std::string &name, const IInterface *parent)
	Constructor with parameters. More...

	~PixelFastDigitizationTool ()
	Destructor. 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	digitize (const EventContext &ctx, TimedHitCollection< SiHit > &thpcsi)

StatusCode	createAndStoreRIOs (const EventContext &ctx)

Private Types
typedef std::multimap< IdentifierHash, InDet::PixelCluster * >	Pixel_detElement_RIO_map

Private Member Functions
	PixelFastDigitizationTool ()

	PixelFastDigitizationTool (const PixelFastDigitizationTool &)

PixelFastDigitizationTool &	operator= (const PixelFastDigitizationTool &)

Trk::DigitizationModule *	buildDetectorModule (const InDetDD::SiDetectorElement *) const

Static Private Member Functions
static bool	areNeighbours (const std::vector< Identifier > &group, const Identifier &rdoID, const InDetDD::SiDetectorElement *, const PixelID &pixelID)

static Amg::Vector3D	CalculateIntersection (const Amg::Vector3D &Point, const Amg::Vector3D &Direction, Amg::Vector2D PlaneBorder, double halfthickness)

static void	Diffuse (HepGeom::Point3D< double > &localEntry, HepGeom::Point3D< double > &localExit, double shiftX, double shiftY)

Private Attributes
TimedHitCollection< SiHit > *	m_thpcsi {}

ServiceHandle< IAthRNGSvc >	m_rndmSvc {this, "RndmSvc", "AthRNGSvc", ""}
	Random number service. More...

StringProperty	m_randomEngineName {this, "RndmEngine", "FastPixelDigitization"}
	Name of the random number stream. More...

const PixelID *	m_pixel_ID {}
	Handle to the ID helper. More...

PublicToolHandle< InDet::ClusterMakerTool >	m_clusterMaker {this, "ClusterMaker", "InDet::ClusterMakerTool/FatrasClusterMaker"}
	ToolHandle to ClusterMaker. More...

ToolHandle< ISiLorentzAngleTool >	m_lorentzAngleTool {this, "LorentzAngleTool", "PixelLorentzAngleTool", "Tool to retreive Lorentz angle"}

BooleanProperty	m_pixUseClusterMaker {this, "PixelUseClusterMaker", true}
	use the pixel cluster maker or not More...

InDet::PixelClusterContainer *	m_pixelClusterContainer {}
	the PixelClusterContainer More...

StringProperty	m_pixel_SiClustersName {this, "PixelClusterContainerName", "PixelClusters"}

ServiceHandle< PileUpMergeSvc >	m_mergeSvc {this, "MergeSvc", "PileUpMergeSvc"}
	PileUp Merge service. More...

IntegerProperty	m_HardScatterSplittingMode {this, "HardScatterSplittingMode", 0, "Control pileup & signal splitting"}
	Process all SiHit or just those from signal or background events. More...

bool	m_HardScatterSplittingSkipper {false}

Pixel_detElement_RIO_map *	m_pixelClusterMap {}

StringProperty	m_prdTruthNamePixel {this, "TruthNamePixel", "PRD_MultiTruthPixel"}

PRD_MultiTruthCollection *	m_pixPrdTruth {}
	the PRD truth map for Pixel measurements More...

PublicToolHandle< InDet::PixelGangedAmbiguitiesFinder >	m_gangedAmbiguitiesFinder {this, "gangedAmbiguitiesFinder", "InDet::PixelGangedAmbiguitiesFinder"}

StringProperty	m_inputObjectName {this, "InputObjectName", "PixelHits"}

std::vector< SiHitCollection * >	m_siHitCollList {}
	name of the sub event hit collections. More...

double	m_pixTanLorentzAngleScalor {1.}
	scale the lorentz angle effect More...

BooleanProperty	m_pixEmulateSurfaceCharge {this, "PixelEmulateSurfaceCharge", true}
	emulate the surface charge More...

DoubleProperty	m_pixSmearPathLength {this, "PixelSmearPathSigma", 0.01}
	the 2. More...

BooleanProperty	m_pixSmearLandau {this, "PixelSmearLandau", true}
	if true : landau else: gauss More...

DoubleProperty	m_pixMinimalPathCut {this, "PixelMinimalPathLength", 0.06}
	the 1. More...

DoubleProperty	m_pixPathLengthTotConv {this, "PixelPathLengthTotConversion", 125.}
	from path length to tot More...

BooleanProperty	m_pixModuleDistortion {this, "PixelEmulateModuleDistortion", true}
	simulationn of module bowing More...

DoubleArrayProperty	m_pixPhiError {this, "PixelErrorPhi", {} }
	phi error when not using the ClusterMaker More...

DoubleArrayProperty	m_pixEtaError {this, "PixelErrorEta", {} }
	eta error when not using the ClusterMaker More...

IntegerProperty	m_pixErrorStrategy {this, "PixelErrorStrategy", 2}
	error strategy for the ClusterMaker More...

DoubleProperty	m_pixDiffShiftBarrX {this, "PixDiffShiftBarrX", 0.005}

DoubleProperty	m_pixDiffShiftBarrY {this, "PixDiffShiftBarrY", 0.005}

DoubleProperty	m_pixDiffShiftEndCX {this, "PixDiffShiftEndCX", 0.008}

DoubleProperty	m_pixDiffShiftEndCY {this, "PixDiffShiftEndCY", 0.008}

DoubleProperty	m_ThrConverted {this, "ThrConverted", 50000}

bool	m_mergeCluster {true}
	enable the merging of neighbour Pixel clusters > More...

short	m_splitClusters {0}
	merging parameter used to define two clusters as neighbour > More...

bool	m_acceptDiagonalClusters {true}
	merging parameter used to define two clusters as neighbour > More...

StringProperty	m_pixelClusterAmbiguitiesMapName {this, "PixelClusterAmbiguitiesMapName", "PixelClusterAmbiguitiesMap"}

InDet::PixelGangedClusterAmbiguities *	m_ambiguitiesMap {}

ServiceHandle< InDetDD::IPixelReadoutManager >	m_pixelReadout {this, "PixelReadoutManager", "PixelReadoutManager", "Pixel readout manager" }

SG::ReadCondHandleKey< PixelChargeCalibCondData >	m_chargeDataKey {this, "PixelChargeCalibCondData", "PixelChargeCalibCondData", "Pixel charge calibration data"}

SG::ReadCondHandleKey< PixelCalib::PixelOfflineCalibData >	m_offlineCalibDataKey {this, "PixelOfflineCalibData", "PixelOfflineCalibData", "Pixel offline calibration data"}

SG::ReadCondHandleKey< PixelDistortionData >	m_distortionKey {this, "PixelDistortionData", "PixelDistortionData", "Output readout distortion data"}

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

PublicToolHandle< Trk::IModuleStepper >	m_digitizationStepper {this, "DigitizationStepper", "Trk::PlanarModuleStepper"}

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 59 of file PixelFastDigitizationTool.h.

Member Typedef Documentation

◆ Pixel_detElement_RIO_map

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

private

Definition at line 110 of file PixelFastDigitizationTool.h.

Constructor & Destructor Documentation

◆ PixelFastDigitizationTool() [1/3]

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

Constructor with parameters.

Definition at line 69 of file PixelFastDigitizationTool.cxx.

                                                                               :
  
   PileUpToolBase(type, name, parent)
 {
 }

◆ ~PixelFastDigitizationTool()

PixelFastDigitizationTool::~PixelFastDigitizationTool ( )

Destructor.

Definition at line 76 of file PixelFastDigitizationTool.cxx.

                                                       {
   if(m_pixelClusterMap) {
     delete m_pixelClusterMap;
   }
 }

◆ PixelFastDigitizationTool() [2/3]

PixelFastDigitizationTool::PixelFastDigitizationTool ( )

private

◆ PixelFastDigitizationTool() [3/3]

PixelFastDigitizationTool::PixelFastDigitizationTool ( const PixelFastDigitizationTool & )

private

Member Function Documentation

◆ areNeighbours()

bool PixelFastDigitizationTool::areNeighbours	(	const std::vector< Identifier > &	group,
		const Identifier &	rdoID,
		const InDetDD::SiDetectorElement *	,
		const PixelID &	pixelID
	)

staticprivate

Definition at line 881 of file PixelFastDigitizationTool.cxx.

 {
   // note: in the PixelClusteringToolBase, m_splitClusters is a variable; here
   // splitClusters was explicitly set to zero, acceptDiagonalClusters = 1
   // so much of the original code is redundant and only one path through the code
   // is possible.
  
   std::vector<Identifier>::const_iterator groupBegin = group.begin();
   std::vector<Identifier>::const_iterator groupEnd = group.end();
  
   int row2 = pixelID.phi_index(rdoID);
   int col2 = pixelID.eta_index(rdoID);
  
   int rowmax = row2;
   bool match=false;
   while (groupBegin!=groupEnd)
     {
       Identifier id = *groupBegin;
       int row1 = pixelID.phi_index(id);
       int col1 = pixelID.eta_index(id);
       if(row1 > rowmax) rowmax = row1;
       int deltarow = abs(row2-row1);
       int deltacol = abs(col2-col1);
  
       // a side in common
       if(deltacol+deltarow < 2) match = true;
       //condition "if (acceptDiagonalClusters !=0") is redundant
       if(deltacol == 1 && deltarow == 1) match = true;
  
       ++groupBegin;
     }
  
  
   return match;
 }

◆ buildDetectorModule()

Trk::DigitizationModule * PixelFastDigitizationTool::buildDetectorModule ( const InDetDD::SiDetectorElement * hitSiDetElement ) const

private

Definition at line 921 of file PixelFastDigitizationTool.cxx.

                                                                                                                             {
  
   const InDetDD::PixelModuleDesign* design(dynamic_cast<const InDetDD::PixelModuleDesign*>(&hitSiDetElement->design()));
  
   if (!design) {
     ATH_MSG_DEBUG ( "Could not get design"<< design) ;
     return nullptr;
   }
  
   //Read from the SiDetectorElement information to build the digitization module
   const double halfThickness = hitSiDetElement->thickness() * 0.5;
   const double halfWidth     = design->width() * 0.5;
   const double halfLength    = design->length() * 0.5;
  
   int binsX = design->rows();
   int binsY = design->columns();
   double numberOfChip = design->numberOfCircuits();
  
   InDetDD::SiCellId cell(0,binsY/2);
   float LongPitch  =design->parameters(cell).width().xEta();
   //std::cout<<"numberOfChip "<<numberOfChip<<" LongPitch "<<LongPitch<<std::endl;
  
   ATH_MSG_VERBOSE("Retrieving infos: halfThickness = " << halfThickness << " --- halfWidth = " << halfWidth << " --- halfLength = " << halfLength );
   ATH_MSG_VERBOSE("Retrieving infos: binsX = " << binsX << " --- binsY = " << binsY << " --- numberOfChip = " << numberOfChip);
  
   int readoutDirection = design->readoutSide();
  
   const bool useLorentzAngle{true};
   const IdentifierHash detElHash = hitSiDetElement->identifyHash();
   float lorentzAngle   = useLorentzAngle ? hitSiDetElement->hitDepthDirection()*hitSiDetElement->hitPhiDirection()*std::atan(m_lorentzAngleTool->getTanLorentzAngle(detElHash)) : 0.;
  
   // added for degugging
  
   // std::cout << "barrel_ec  = " << m_detID->barrel_ec(hitSiDetElement->identify())
   //           << " --  layer_disk = " << m_detID->layer_disk(hitSiDetElement->identify())
   //           << " --  eta_module = " << m_detID->eta_module(hitSiDetElement->identify())
   //           << " --  lorentzAngle = " << lorentzAngle
   //        << " --  lorentzCorrection = " << shift
   //        << " --  readoutDirection = " << readoutDirection << std::endl;
   // std::cout << "element->hitDepthDirection() = " << hitSiDetElement->hitDepthDirection()
   //        << " --  element->hitPhiDirection() = " << hitSiDetElement->hitPhiDirection() << std::endl;
  
   // rectangle bounds
   auto rectangleBounds = std::make_shared<const Trk::RectangleBounds>(halfWidth,halfLength);
   ATH_MSG_VERBOSE("Initialized rectangle Bounds");
   // create the segmentation
   std::shared_ptr<const Trk::Segmentation> rectangleSegmentation(new Trk::RectangularSegmentation(std::move(rectangleBounds),(size_t)binsX,LongPitch,(size_t)binsY, numberOfChip));
   // build the module
   ATH_MSG_VERBOSE("Initialized rectangleSegmentation");
   Trk::DigitizationModule * digitizationModule = new Trk::DigitizationModule(std::move(rectangleSegmentation),
                                                                                halfThickness,
                                                                                readoutDirection,
                                                                                lorentzAngle);
   ATH_MSG_VERBOSE("Building Rectangle Segmentation with dimensions (halfX, halfY) = (" << halfWidth << ", " << halfLength << ")");
  
   // success return
   return digitizationModule;
 }

◆ CalculateIntersection()

Amg::Vector3D PixelFastDigitizationTool::CalculateIntersection	(	const Amg::Vector3D &	Point,
		const Amg::Vector3D &	Direction,
		Amg::Vector2D	PlaneBorder,
		double	halfthickness
	)

staticprivate

Definition at line 981 of file PixelFastDigitizationTool.cxx.

 {
   Amg::Vector3D Intersection(0.,0.,0.);
  
   //Using parameter definition of a line in 3d z=z_point + direction_z t
   std::vector<double> parameters;
   parameters.push_back((PlaneBorder.x() - Point.x())/Direction.x());
   parameters.push_back((-PlaneBorder.x() - Point.x())/Direction.x());
   parameters.push_back((PlaneBorder.y() - Point.y())/Direction.y());
   parameters.push_back((-PlaneBorder.y() - Point.y())/Direction.y());
  
   for(double parameter: parameters)
     {
       double z =  Point.z() + Direction.z() * parameter;
       if( std::abs(z) > halfthickness )
         continue;
  
  
       double x = Point.x() + Direction.x() * parameter;
       double y = Point.y() + Direction.y() * parameter;
  
       if(std::abs(x) > PlaneBorder.x() || std::abs(y) > PlaneBorder.y())
         continue;
  
  
       Intersection = Amg::Vector3D(x,y,z);
       break;
  
     }
  
   return Intersection;
 }

◆ createAndStoreRIOs()

StatusCode PixelFastDigitizationTool::createAndStoreRIOs ( const EventContext & ctx )

Definition at line 821 of file PixelFastDigitizationTool.cxx.

 {
   SG::ReadCondHandle<InDetDD::SiDetectorElementCollection> pixelDetEleHandle(m_pixelDetEleCollKey, ctx);
   const InDetDD::SiDetectorElementCollection* elements(*pixelDetEleHandle);
   if (not pixelDetEleHandle.isValid() or elements==nullptr) {
     ATH_MSG_FATAL(m_pixelDetEleCollKey.fullKey() << " is not available.");
     return StatusCode::FAILURE;
   }
  
   Pixel_detElement_RIO_map::iterator i = m_pixelClusterMap->begin();
   Pixel_detElement_RIO_map::iterator e = m_pixelClusterMap->end();
  
   InDet::PixelClusterCollection* clusterCollection = nullptr;
   IdentifierHash waferHash;
  
   for (; i != e; i = m_pixelClusterMap->upper_bound(i->first)){
  
     std::pair <Pixel_detElement_RIO_map::iterator, Pixel_detElement_RIO_map::iterator> range;
     range = m_pixelClusterMap->equal_range(i->first);
  
     Pixel_detElement_RIO_map::iterator firstDetElem;
     firstDetElem = range.first;
  
     waferHash = firstDetElem->first;
  
     const InDetDD::SiDetectorElement* detElement = elements->getDetectorElement(waferHash);
  
     clusterCollection = new InDet::PixelClusterCollection(waferHash);
     clusterCollection->setIdentifier(detElement->identify());
  
  
     for ( Pixel_detElement_RIO_map::iterator iter = range.first; iter != range.second; ++iter){
  
       InDet::PixelCluster* pixelCluster = (*iter).second;
       pixelCluster->setHashAndIndex(clusterCollection->identifyHash(),clusterCollection->size());
       clusterCollection->push_back(pixelCluster);
  
     }
  
  
     if (clusterCollection) {
       if (!clusterCollection->empty()) {
         ATH_MSG_DEBUG ( "Filling ambiguities map" );
         m_gangedAmbiguitiesFinder->execute(clusterCollection,*m_ambiguitiesMap);
         ATH_MSG_DEBUG ( "Ambiguities map: " << m_ambiguitiesMap->size() << " elements" );
         if ((m_pixelClusterContainer->addCollection(clusterCollection, waferHash)).isFailure()){
           ATH_MSG_WARNING( "Could not add collection to Identifyable container !" );
         }
       }
       else {delete clusterCollection;}
     }
  
  
   } // end for
  
   m_pixelClusterMap->clear();
  
   return StatusCode::SUCCESS;
 }

◆ Diffuse()

void PixelFastDigitizationTool::Diffuse	(	HepGeom::Point3D< double > &	localEntry,
		HepGeom::Point3D< double > &	localExit,
		double	shiftX,
		double	shiftY
	)

staticprivate

Definition at line 1014 of file PixelFastDigitizationTool.cxx.

                                                                                                                                            {
  
   double localEntryX = localEntry.x();
   double localEntryY = localEntry.y();
   double localExitX = localExit.x();
   double localExitY = localExit.y();
  
   double signX = localExitX>localEntryX ? 1:-1 ;
   double signY = localExitY>localEntryY ? 1:-1 ;
  
   localEntryX += shiftX*signX*(-1);
   localExitX += shiftX*signX;
   localEntryY += shiftY*signY*(-1);
   localExitY += shiftY*signY;
  
   localEntry.setX(localEntryX);
   localEntry.setY(localEntryY);
   localExit.setX(localExitX);
   localExit.setY(localExitY);
  
 }

◆ digitize()

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

Definition at line 403 of file PixelFastDigitizationTool.cxx.

 {
   // Set the RNG to use for this event.
   ATHRNG::RNGWrapper* rngWrapper = m_rndmSvc->getEngine(this, m_randomEngineName);
   const std::string rngName = name()+m_randomEngineName;
   rngWrapper->setSeed( rngName, ctx );
   CLHEP::HepRandomEngine *rndmEngine = rngWrapper->getEngine(ctx);
  
   SG::ReadCondHandle<InDetDD::SiDetectorElementCollection> pixelDetEleHandle(m_pixelDetEleCollKey, ctx);
   const InDetDD::SiDetectorElementCollection* elements(*pixelDetEleHandle);
   if (not pixelDetEleHandle.isValid() or elements==nullptr) {
     ATH_MSG_FATAL(m_pixelDetEleCollKey.fullKey() << " is not available.");
     return StatusCode::FAILURE;
   }
  
   TimedHitCollection<SiHit>::const_iterator i, e;
  
   if(!m_pixelClusterMap) { m_pixelClusterMap = new Pixel_detElement_RIO_map; }
   else { m_pixelClusterMap->clear(); }
  
   SG::ReadCondHandle<PixelChargeCalibCondData> calibDataHandle(m_chargeDataKey, ctx);
   const PixelChargeCalibCondData *calibData = *calibDataHandle;
   SG::ReadCondHandle<PixelCalib::PixelOfflineCalibData> offlineCalibData(m_offlineCalibDataKey, ctx);
   std::vector<int> truthIdList;
   std::vector<Identifier> detEl;
  
   while (thpcsi.nextDetectorElement(i, e)) {
  
     Pixel_detElement_RIO_map PixelDetElClusterMap;
  
     truthIdList.clear();
     detEl.clear();
  
     while (i != e) {
  
       const TimedHitPtr<SiHit>& hit(*i++);
       // check the status of truth information for this SiHit
       // some Truth information is cut for pile up events
       const HepMcParticleLink currentLink = HepMcParticleLink::getRedirectedLink(hit->particleLink(), hit.eventId(), ctx); // This link should now correctly resolve to the TruthEvent McEventCollection in the main StoreGateSvc.
       const int barrelEC = hit->getBarrelEndcap();
       const int layerDisk = hit->getLayerDisk();
       const int phiModule = hit->getPhiModule();
       const int etaModule = hit->getEtaModule();
  
       const Identifier moduleID = m_pixel_ID->wafer_id(barrelEC, layerDisk, phiModule, etaModule);
       const IdentifierHash waferHash = m_pixel_ID->wafer_hash(moduleID);
       const InDetDD::SiDetectorElement* hitSiDetElement = elements->getDetectorElement(waferHash);
       if (!hitSiDetElement) {ATH_MSG_ERROR( " could not get detector element "); continue;}
  
       if (!(hitSiDetElement->isPixel())) {continue;}
  
  
  
       std::vector<HepMcParticleLink> hit_vector; //Store the hits in merged cluster
  
       const int truthID = (currentLink.barcode() !=0 && currentLink.id() == 0) ? 3 : currentLink.id(); // FIXME barcode-based Patch for reading in legacy barcode-based EDM - if the barcode is non-zero, but the id is zero then we must be looking at a particle linked to suppressed pile-up truth - such SiHits would be linked to the third GenParticle in their GenEvents (if they were present)
  
       const Identifier hitId = hitSiDetElement->identify(); // Isn't this is identical to moduleID?
       //const IdentifierHash hitIdHash = hitSiDetElement->identifyHash();
  
  
       bool isRep = false;
  
       for (int j : truthIdList) {
         for (auto & k : detEl) {
           if ((truthID > 0) && (truthID == j) && (hitId == k)) {isRep = true; break;}
         }
         if (isRep) break;
       }
  
       if (isRep) continue;
  
       truthIdList.push_back(truthID);
       detEl.push_back(hitId);
  
       HepGeom::Point3D<double> localStartPosition = hit->localStartPosition();
       HepGeom::Point3D<double> localEndPosition = hit->localEndPosition();
  
       localStartPosition = hitSiDetElement->hitLocalToLocal3D(localStartPosition);
       localEndPosition = hitSiDetElement->hitLocalToLocal3D(localEndPosition);
  
       int isEndcap = (barrelEC==0) ? 0:1;
  
       double shiftX = isEndcap ? m_pixDiffShiftEndCX : m_pixDiffShiftBarrX;
       double shiftY = isEndcap ? m_pixDiffShiftEndCY : m_pixDiffShiftBarrY;
  
       //       std::cout<<"Thr "<<m_ThrConverted<<std::endl;
  
       //New function to tune the cluster size
       Diffuse(localStartPosition, localEndPosition, shiftX, shiftY);
  
       const Amg::Vector3D localDirection(localEndPosition.x()-localStartPosition.x(), localEndPosition.y()-localStartPosition.y(), localEndPosition.z()-localStartPosition.z());
  
       Amg::Vector3D entryPoint(localStartPosition.x(),localStartPosition.y(),localStartPosition.z());
       Amg::Vector3D exitPoint(localEndPosition.x(),localEndPosition.y(),localEndPosition.z());
  
       const InDetDD::PixelModuleDesign* design(dynamic_cast<const InDetDD::PixelModuleDesign*>(&hitSiDetElement->design()));
       if (not design){
         ATH_MSG_FATAL("Failed to cast the hitSiDetElement::design pointer to a PixelModuleDesign*; aborting due to null pointer");
         return StatusCode::FAILURE;
       }
       const Amg::Vector2D localEntry(localStartPosition.x(),localStartPosition.y());
       const Amg::Vector2D localExit(localEndPosition.x(),localEndPosition.y());
  
       Identifier entryId = hitSiDetElement->identifierOfPosition(localEntry);
       Identifier exitId  = hitSiDetElement->identifierOfPosition(localExit);
  
       InDetDD::SiCellId entryCellId = hitSiDetElement->cellIdFromIdentifier(entryId);
       InDetDD::SiCellId exitCellId = hitSiDetElement->cellIdFromIdentifier(exitId);
  
       double halfthickness = hitSiDetElement->thickness()*0.5;
  
       bool EntryValid(entryCellId.isValid());
       bool ExitValid(exitCellId.isValid());
  
       double pixMinimalPathCut= 1. / m_pixPathLengthTotConv;
  
       Identifier diodeID = hitId;
       unsigned int FE = m_pixelReadout->getFE(diodeID, moduleID);
       InDetDD::PixelDiodeType type = m_pixelReadout->getDiodeType(diodeID);
  
       double th0 = calibData->getThresholds(type, waferHash, FE).value/ m_ThrConverted;
  
       //        if (old_th != th0) std::cout<<"converted threshold "<<th0<<std::endl, old_th= th0;
  
       //Avoid to store pixels with 0 ToT
       //if (m_pixMinimalPathCut > pixMinimalPathCut) pixMinimalPathCut=m_pixMinimalPathCut;
       //if (th0 > pixMinimalPathCut) ;
       pixMinimalPathCut=th0;
  
       if (!EntryValid || !ExitValid)
         {
           //If entry or exit aren't valid search for the valid entry/exit of the hit as Intersection with the module
  
           if ( !EntryValid && !ExitValid) continue;
  
           Amg::Vector3D Point = EntryValid ? entryPoint : exitPoint ;
  
           Amg::Vector3D Intersection =  CalculateIntersection(Point, localDirection, Amg::Vector2D(design->width() * 0.5,design->length() * 0.5),halfthickness);
  
           if( Intersection == Amg::Vector3D(0.,0.,0.))
             {
               ATH_MSG_WARNING("ATTENTION THE INTERSECTION COORDINATE IS OUT OF THE MODULE");
               continue;
             }
  
           const Amg::Vector2D Intersection_2d(Intersection.x(),Intersection.y());
           Identifier Intersection_2dId = hitSiDetElement->identifierOfPosition(Intersection_2d);
  
           InDetDD::SiCellId Intersection_2dCellId = hitSiDetElement->cellIdFromIdentifier(Intersection_2dId);
  
  
           if(!Intersection_2dCellId.isValid()) continue;
  
           if(EntryValid)
             exitPoint = Intersection;
           else
             entryPoint = Intersection;
  
  
         }
  
  
  
       Trk::DigitizationModule * digitizationModule = buildDetectorModule(hitSiDetElement);
       if(!digitizationModule){
         ATH_MSG_FATAL( " could not get build detector module ");
         return StatusCode::FAILURE;
       }
  
       // Getting the steps in the sensor
       std::vector<Trk::DigitizationStep> digitizationSteps = m_digitizationStepper->cellSteps(*digitizationModule,entryPoint,exitPoint);
  
  
  
       // the pixel positions and other needed stuff for the geometrical clustering
       std::unique_ptr<InDet::PixelCluster> pixelCluster = nullptr;
       Amg::Vector2D       clusterPosition(0.,0.);
  
       std::vector<Identifier>           rdoList;
       std::vector<int>                  totList;
  
       const bool   isGanged = false;
       int lvl1a = 0;
  
       double accumulatedPathLength=0.;
  
       //ATTENTION index max e min da rdo + manager
       int phiIndexMax = -999999;
       int phiIndexMin = 1000000;
       int etaIndexMax = -999999;
       int etaIndexMin = 1000000;
  
  
       for (auto& dStep : digitizationSteps){
  
         double pathlength = dStep.stepLength;
         // two options fro charge smearing: landau / gauss
         if ( m_pixSmearPathLength > 0. ) {
           // create the smdar parameter
           double sPar = m_pixSmearLandau ?
             m_pixSmearPathLength*CLHEP::RandLandau::shoot(rndmEngine) :
             m_pixSmearPathLength*CLHEP::RandGaussZiggurat::shoot(rndmEngine);
           pathlength *=  (1.+sPar);
         }
  
  
         if (pathlength < pixMinimalPathCut) continue;
  
         // position on the diode map
         Trk::DigitizationCell cell(dStep.stepCell.first,dStep.stepCell.second);
         Amg::Vector2D PositionOnModule = digitizationModule->segmentation().cellPosition(cell);
         InDetDD::SiCellId diode = hitSiDetElement->cellIdOfPosition(PositionOnModule);
  
  
         if (!diode.isValid())
           continue;
  
         Amg::Vector2D chargeCenterPosition = hitSiDetElement->rawLocalPositionOfCell(diode);
  
         const Identifier rdoId            =  hitSiDetElement->identifierOfPosition(chargeCenterPosition);
         clusterPosition += pathlength * chargeCenterPosition;
  
         int currentEtaIndex = diode.etaIndex();
         int currentPhiIndex = diode.phiIndex();
         if(currentEtaIndex > etaIndexMax) etaIndexMax = currentEtaIndex;
         if(currentEtaIndex < etaIndexMin) etaIndexMin = currentEtaIndex;
         if(currentPhiIndex > phiIndexMax) phiIndexMax = currentPhiIndex;
         if(currentPhiIndex < phiIndexMin) phiIndexMin = currentPhiIndex;
  
  
         // record - positions, rdoList and totList
         accumulatedPathLength += pathlength;
         //Fail
         rdoList.push_back(rdoId);
         totList.push_back(int(pathlength*m_pixPathLengthTotConv));
  
       }
  
       delete digitizationModule;
  
       // the col/row
       int siDeltaPhiCut = phiIndexMax-phiIndexMin+1;
       int siDeltaEtaCut = etaIndexMax-etaIndexMin+1;
  
       int totalToT=std::accumulate(totList.begin(), totList.end(), 0);;
  
       // bail out if 0 pixel or path length problem
       if (rdoList.empty() || accumulatedPathLength < pixMinimalPathCut || totalToT == 0) {
         if (totalToT == 0 && !rdoList.empty() ) ATH_MSG_WARNING("The total ToT of the cluster is 0, this should never happen");
         continue;
       }
  
       // weight the cluster position
       clusterPosition *= 1./accumulatedPathLength;
       Identifier clusterId = hitSiDetElement->identifierOfPosition(clusterPosition);
  
  
       // merging clusters
  
       bool merged = false;
       if(m_mergeCluster){ // merge to the current cluster "near" cluster in the cluster map, in the current detector element
  
         for(Pixel_detElement_RIO_map::iterator currentClusIter = PixelDetElClusterMap.begin(); currentClusIter != PixelDetElClusterMap.end();) {
           //make a temporary to use within the loop and possibly erase - increment the main interator at the same time.
           Pixel_detElement_RIO_map::iterator clusIter = currentClusIter++;
           InDet::PixelCluster* currentCluster = clusIter->second;
           const std::vector<Identifier> &currentRdoList = currentCluster->rdoList();
           bool areNb = false;
           for (auto rdoIter : rdoList) {
             areNb = PixelFastDigitizationTool::areNeighbours(currentRdoList, rdoIter, hitSiDetElement,*m_pixel_ID);
             if (areNb) { break; }
           }
           if (areNb) {
             const std::vector<int> &currentTotList = currentCluster->totList();
             rdoList.insert(rdoList.end(), currentRdoList.begin(), currentRdoList.end() );
             totList.insert(totList.end(), currentTotList.begin(), currentTotList.end() );
             Amg::Vector2D       currentClusterPosition(currentCluster->localPosition());
             float c1 = (float)currentRdoList.size();
             float c2 = (float)rdoList.size();
             clusterPosition = (clusterPosition*c2 + currentClusterPosition*c1)/((c1 + c2));
             clusterId = hitSiDetElement->identifierOfPosition(clusterPosition);
             merged = true;
             PixelDetElClusterMap.erase(clusIter);
  
             //Store HepMcParticleLink connected to the cluster removed from the collection
             std::pair<PRD_MultiTruthCollection::iterator,PRD_MultiTruthCollection::iterator> saved_hit = m_pixPrdTruth->equal_range(currentCluster->identify());
             for (PRD_MultiTruthCollection::iterator this_hit = saved_hit.first; this_hit != saved_hit.second; ++this_hit)
               {
                 hit_vector.push_back(this_hit->second);
               }
             //Delete all the occurency of the currentCluster from the multi map
             if (saved_hit.first != saved_hit.second) m_pixPrdTruth->erase(currentCluster->identify());
             delete currentCluster;
             //break; //commenting out this break statement allows for multiple existing clusters to be merged.
           }
         }
       }
  
       bool not_valid = false;
       for (auto & entry : rdoList) {
         if (!(entry.is_valid())) { not_valid = true; break;}
       }
  
       if (not_valid) continue;
  
       if(merged) {
         //Hacks for merged clusters
         for (auto rdoIter : rdoList) {
           const InDetDD::SiCellId& chargeCellId =  hitSiDetElement->cellIdFromIdentifier(rdoIter);
           // phi/ eta index
           int chargePhiIndex = chargeCellId.phiIndex();
           int chargeEtaIndex = chargeCellId.etaIndex();
           // set max min
           phiIndexMin = chargePhiIndex < phiIndexMin ?  chargePhiIndex : phiIndexMin;
           phiIndexMax = chargePhiIndex > phiIndexMax ?  chargePhiIndex : phiIndexMax;
           etaIndexMin = chargeEtaIndex < etaIndexMin ?  chargeEtaIndex : etaIndexMin;
           etaIndexMax = chargeEtaIndex > etaIndexMax ?  chargeEtaIndex : etaIndexMax;
         }
         siDeltaPhiCut = (phiIndexMax-phiIndexMin)+1;
         siDeltaEtaCut = (etaIndexMax-etaIndexMin)+1;
       }
  
       // ---------------------------------------------------------------------------------------------
       //  PART 2: Cluster && ROT creation
  
  
       //ATTENTION
       //       // correct shift implied by the scaling of the Lorentz angle
       //       double newshift = 0.5*thickness*tanLorAng;
       //       double corr = ( shift - newshift );
       // 2a) Cluster creation ------------------------------------
       if (m_pixUseClusterMaker) {
  
         // from InDetReadoutGeometry: width from eta
         const auto* pixModDesign =
             dynamic_cast<const InDetDD::PixelModuleDesign*>(
                 &hitSiDetElement->design());
         if (!pixModDesign) {
           return StatusCode::FAILURE;
         }
         double etaWidth =
             pixModDesign->widthFromColumnRange(etaIndexMin, etaIndexMax);
         // from InDetReadoutGeometry : width from phi
         double phiWidth =
             pixModDesign->widthFromRowRange(phiIndexMin, phiIndexMax);
  
         InDet::SiWidth siWidth(Amg::Vector2D(siDeltaPhiCut, siDeltaEtaCut),
                                Amg::Vector2D(phiWidth, etaWidth));
  
         // use the cluster maker from the offline software
         pixelCluster =
             std::make_unique<PixelCluster>(m_clusterMaker->pixelCluster(
                 clusterId, clusterPosition, std::vector<Identifier>(rdoList),
                 lvl1a, std::vector<int>(totList), siWidth,
                 hitSiDetElement, isGanged, m_pixErrorStrategy, *m_pixel_ID,
                 false, 0.0, 0.0, calibData, *offlineCalibData));
         if (isGanged) {
           pixelCluster->setGangedPixel(isGanged);
         }
       } else {
         ATH_MSG_WARNING(
             "[ cluster - pix ] No pixels errors provided, but configured to "
             "use them.");
         ATH_MSG_WARNING(
             "                  -> No pixels cluster will be created.");
         continue;
       }
  
       if (!(pixelCluster->identify().is_valid())) {
         pixelCluster.reset();
         continue;
       }
  
       if (!(m_pixel_ID->is_pixel(pixelCluster->identify()))) {
         pixelCluster.reset();
         continue;
       }
  
       const auto it =
           PixelDetElClusterMap.insert(Pixel_detElement_RIO_map::value_type(
               waferHash, pixelCluster.release()));
       const PixelCluster* insertedCluster = it->second;
  
       if (currentLink.isValid()) {
         if (!HepMC::ignoreTruthLink(currentLink, m_vetoPileUpTruthLinks)) {
           m_pixPrdTruth->insert(
               std::make_pair(insertedCluster->identify(), currentLink));
           ATH_MSG_DEBUG("Truth map filled with cluster"
                         << insertedCluster
                         << " and link = " << currentLink);
         }
       } else {
         ATH_MSG_DEBUG(
             "Particle link NOT valid!! Truth map NOT filled with cluster"
             << insertedCluster << " and link = " << currentLink);
       }
  
       // Add all hit that was connected to the cluster
       for (const HepMcParticleLink& p : hit_vector) {
  
         m_pixPrdTruth->insert(std::make_pair(insertedCluster->identify(), p ));
       }
  
       hit_vector.clear();
     } // end hit while
  
  
     (void) m_pixelClusterMap->insert(PixelDetElClusterMap.begin(), PixelDetElClusterMap.end());
  
  
   } // end nextDetectorElement while
  
  
   return StatusCode::SUCCESS;
 }

◆ filterPassed()

virtual bool PileUpToolBase::filterPassed ( ) const

inlineoverridevirtualinherited

dummy implementation of passing filter

Definition at line 49 of file PileUpToolBase.h.

49 { return m_filterPassed; }

◆ initialize()

StatusCode PixelFastDigitizationTool::initialize ( )

virtual

Reimplemented from PileUpToolBase.

Definition at line 83 of file PixelFastDigitizationTool.cxx.

 {
  
   ATH_MSG_DEBUG ( "PixelDigitizationTool::initialize()" );
  
   ATH_CHECK(m_pixelReadout.retrieve());
   ATH_CHECK(m_chargeDataKey.initialize());
   ATH_CHECK(m_offlineCalibDataKey.initialize());
   ATH_CHECK(m_pixelDetEleCollKey.initialize());
  
   //locate the AtRndmGenSvc and initialize our local ptr
   if (!m_rndmSvc.retrieve().isSuccess())
     {
       ATH_MSG_ERROR ( "Could not find given RndmSvc" );
       return StatusCode::FAILURE;
     }
  
   if (detStore()->retrieve(m_pixel_ID, "PixelID").isFailure()) {
     ATH_MSG_ERROR ( "Could not get Pixel ID helper" );
     return StatusCode::FAILURE;
   }
  
   if (m_inputObjectName.empty())
     {
       ATH_MSG_FATAL ( "Property InputObjectName not set !" );
       return StatusCode::FAILURE;
     }
   else
     {
       ATH_MSG_DEBUG ( "Input objects: '" << m_inputObjectName << "'" );
     }
  
   // retrieve the offline cluster maker : for pixel and/or sct
   if ( m_pixUseClusterMaker) {
     if (m_clusterMaker.retrieve().isFailure()){
       ATH_MSG_WARNING( "Could not retrieve " << m_clusterMaker );
       ATH_MSG_WARNING( "-> Switching to simplified cluster creation!" );
       m_pixUseClusterMaker = false;
       m_clusterMaker.disable();
     }
   } else {
     m_clusterMaker.disable();
   }
  
   if (m_pixModuleDistortion) {
     ATH_CHECK(m_distortionKey.initialize());
   }
  
   //locate the PileUpMergeSvc and initialize our local ptr
   if (!m_mergeSvc.retrieve().isSuccess()) {
     ATH_MSG_ERROR ( "Could not find PileUpMergeSvc" );
     return StatusCode::FAILURE;
   }
  
  
   // get the InDet::PixelGangedAmbiguitiesFinder
   if ( m_gangedAmbiguitiesFinder.retrieve().isFailure() ) {
     ATH_MSG_FATAL( m_gangedAmbiguitiesFinder.propertyName() << ": Failed to retrieve tool " << m_gangedAmbiguitiesFinder.type() );
     return StatusCode::FAILURE;
   } else {
     ATH_MSG_DEBUG ( m_gangedAmbiguitiesFinder.propertyName() << ": Retrieved tool " << m_gangedAmbiguitiesFinder.type() );
   }
  
   ATH_CHECK(m_lorentzAngleTool.retrieve());
  
   return StatusCode::SUCCESS;
 }

◆ mergeEvent()

StatusCode PixelFastDigitizationTool::mergeEvent ( const EventContext & ctx )

Definition at line 320 of file PixelFastDigitizationTool.cxx.

 {
  
   m_pixelClusterContainer = new InDet::PixelClusterContainer(m_pixel_ID->wafer_hash_max());
  
   if(!m_pixelClusterContainer) {
     ATH_MSG_FATAL( "[ --- ] Could not create PixelClusterContainer");
     return StatusCode::FAILURE;
   }
  
   InDet::SiClusterContainer* symSiContainer=nullptr;
  
   // --------------------------------------
   // Pixel_Cluster container registration
   m_pixelClusterContainer->cleanup();
   if ((evtStore()->record(m_pixelClusterContainer, "PixelClusters")).isFailure())   {
     ATH_MSG_FATAL("[ hitproc ] Error while registering PixelCluster container");
     return StatusCode::FAILURE;
   }
  
   // symlink the SCT Container
   if ((evtStore()->symLink(m_pixelClusterContainer,symSiContainer)).isFailure()) {
     ATH_MSG_FATAL( "[ --- ] PixelClusterContainer could not be symlinked to SiClusterContainter in StoreGate !" );
     return StatusCode::FAILURE;
   } else {
     ATH_MSG_DEBUG( "[ hitproc ] PixelClusterContainer symlinked to SiClusterContainer in StoreGate" );
   }
  
   // truth info
   m_pixPrdTruth = new PRD_MultiTruthCollection;
  
   if ((evtStore()->contains<PRD_MultiTruthCollection>(m_prdTruthNamePixel))){
     if((evtStore()->retrieve(m_pixPrdTruth, m_prdTruthNamePixel)).isFailure()){
       ATH_MSG_FATAL("Could not retrieve collection " << m_prdTruthNamePixel);
       return StatusCode::FAILURE;
     }
   }else{
     if((evtStore()->record(m_pixPrdTruth, m_prdTruthNamePixel)).isFailure()){
       ATH_MSG_FATAL("Could not record collection " << m_prdTruthNamePixel);
       return StatusCode::FAILURE;
     }
   }
  
   m_ambiguitiesMap =new PixelGangedClusterAmbiguities();
  
   if (m_thpcsi != nullptr) {
     if(digitize(ctx, *m_thpcsi).isFailure()) {
       ATH_MSG_FATAL ( "Pixel digitize method failed!" );
       return StatusCode::FAILURE;
     }
   }
  
   delete m_thpcsi;
   for(SiHitCollection* ptr : m_siHitCollList) delete ptr;
   m_siHitCollList.clear();
  
  
   if (createAndStoreRIOs(ctx).isFailure()) {
     ATH_MSG_FATAL ( "createAndStoreRIOs() failed!" );
     return StatusCode::FAILURE;
   }
   else {
     ATH_MSG_DEBUG ( "createAndStoreRIOs() succeeded" );
   }
  
   if ((evtStore()->setConst(m_pixelClusterContainer)).isFailure()) {
     ATH_MSG_ERROR("[ ---- ] Could not set Pixel ROT container ");
   }
   ATH_MSG_DEBUG ("Ambiguities map has " << m_ambiguitiesMap->size() << " elements" );
   StatusCode sc = evtStore()->record(m_ambiguitiesMap,m_pixelClusterAmbiguitiesMapName,false);
   if (sc.isFailure()){
     ATH_MSG_FATAL ( "PixelClusterAmbiguitiesMap could not be recorded in StoreGate !" );
     return StatusCode::FAILURE;
   }else{
     ATH_MSG_DEBUG ( "PixelClusterAmbiguitiesMap recorded in StoreGate" );
   }
  
  
  
   return StatusCode::SUCCESS;
 }

◆ operator=()

PixelFastDigitizationTool& PixelFastDigitizationTool::operator= ( const PixelFastDigitizationTool & )

private

◆ prepareEvent()

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

Definition at line 153 of file PixelFastDigitizationTool.cxx.

 {
  
   m_siHitCollList.clear();
   m_thpcsi = new TimedHitCollection<SiHit>();
   m_HardScatterSplittingSkipper = false;
   return StatusCode::SUCCESS;
 }

◆ processAllSubEvents() [1/2]

StatusCode PixelFastDigitizationTool::processAllSubEvents ( const EventContext & ctx )

virtual

Reimplemented from PileUpToolBase.

Definition at line 204 of file PixelFastDigitizationTool.cxx.

                                                                                  {
  
   m_pixelClusterContainer = new InDet::PixelClusterContainer(m_pixel_ID->wafer_hash_max());
  
   if(!m_pixelClusterContainer) {
     ATH_MSG_FATAL( "[ --- ] Could not create PixelClusterContainer");
     return StatusCode::FAILURE;
   }
  
   InDet::SiClusterContainer* symSiContainer=nullptr;
  
   // --------------------------------------
   // Pixel Cluster container registration
   m_pixelClusterContainer->cleanup();
   if ((evtStore()->record(m_pixelClusterContainer, m_pixel_SiClustersName)).isFailure())   {
     ATH_MSG_FATAL("[ hitproc ] Error while registering PixelCluster container");
     return StatusCode::FAILURE;
   }
  
   // symlink the Pixel Cluster Container
   if ((evtStore()->symLink(m_pixelClusterContainer,symSiContainer)).isFailure()) {
     ATH_MSG_FATAL( "[ --- ] PixelClusterContainer could not be symlinked to SiClusterContainter in StoreGate !" );
     return StatusCode::FAILURE;
   } else {
     ATH_MSG_DEBUG( "[ hitproc ] PixelClusterContainer symlinked to SiClusterContainer in StoreGate" );
   }
  
   // truth info
  
   m_pixPrdTruth = new PRD_MultiTruthCollection;
  
   if ((evtStore()->contains<PRD_MultiTruthCollection>(m_prdTruthNamePixel))){
     if((evtStore()->retrieve(m_pixPrdTruth, m_prdTruthNamePixel)).isFailure()){
       ATH_MSG_FATAL("Could not retrieve collection " << m_prdTruthNamePixel);
       return StatusCode::FAILURE;
     }
   }else{
     if((evtStore()->record(m_pixPrdTruth, m_prdTruthNamePixel)).isFailure()){
       ATH_MSG_FATAL("Could not record collection " << m_prdTruthNamePixel);
       return StatusCode::FAILURE;
     }
   }
  
  
   m_ambiguitiesMap =new PixelGangedClusterAmbiguities();
  
  
   //  get the container(s)
   using TimedHitCollList = PileUpMergeSvc::TimedList<SiHitCollection>::type;
  
   //this is a list<pair<time_t, DataLink<SCTUncompressedHitCollection> > >
   TimedHitCollList hitCollList;
   unsigned int numberOfSimHits(0);
   if ( !(m_mergeSvc->retrieveSubEvtsData(m_inputObjectName.value(), hitCollList, numberOfSimHits).isSuccess()) && hitCollList.empty() ) {
     ATH_MSG_ERROR ( "Could not fill TimedHitCollList" );
     return StatusCode::FAILURE;
   } else {
     ATH_MSG_DEBUG ( hitCollList.size() << " SiHitCollections with key " << m_inputObjectName << " found" );
   }
  
   // Define Hit Collection
   TimedHitCollection<SiHit> thpcsi(numberOfSimHits);
  
   //now merge all collections into one
   TimedHitCollList::iterator   iColl(hitCollList.begin());
   TimedHitCollList::iterator endColl(hitCollList.end()  );
  
   m_HardScatterSplittingSkipper = false;
   // loop on the hit collections
   while ( iColl != endColl ) {
     //decide if this event will be processed depending on HardScatterSplittingMode & bunchXing
     if (m_HardScatterSplittingMode == 2 && !m_HardScatterSplittingSkipper ) { m_HardScatterSplittingSkipper = true; ++iColl; continue; }
     if (m_HardScatterSplittingMode == 1 && m_HardScatterSplittingSkipper )  { ++iColl; continue; }
     if (m_HardScatterSplittingMode == 1 && !m_HardScatterSplittingSkipper ) { m_HardScatterSplittingSkipper = true; }
     const SiHitCollection* p_collection(iColl->second);
     thpcsi.insert(iColl->first, p_collection);
     ATH_MSG_DEBUG ( "SiHitCollection found with " << p_collection->size() << " hits" );
     ++iColl;
   }
  
   // Process the Hits straw by straw: get the iterator pairs for given straw
   if(this->digitize(ctx, thpcsi).isFailure()) {
     ATH_MSG_FATAL ( "digitize method failed!" );
     return StatusCode::FAILURE;
   }
  
   if (createAndStoreRIOs(ctx).isFailure()) {
     ATH_MSG_FATAL ( "createAndStoreRIOs() failed!" );
     return StatusCode::FAILURE;
   }
   else {
     ATH_MSG_DEBUG ( "createAndStoreRIOs() succeeded" );
   }
  
  
   if ((evtStore()->setConst(m_pixelClusterContainer)).isFailure()) {
     ATH_MSG_ERROR("[ ---- ] Could not set Pixel ROT container ");
   }
   ATH_MSG_DEBUG ("Ambiguities map has " << m_ambiguitiesMap->size() << " elements" );
   StatusCode sc = evtStore()->record(m_ambiguitiesMap,m_pixelClusterAmbiguitiesMapName,false);
   if (sc.isFailure()){
     ATH_MSG_FATAL ( "PixelClusterAmbiguitiesMap could not be recorded in StoreGate !" );
     return StatusCode::FAILURE;
   }else{
     ATH_MSG_DEBUG ( "PixelClusterAmbiguitiesMap recorded in StoreGate" );
   }
  
  
  
  
   return StatusCode::SUCCESS;
  
 }

◆ 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 PixelFastDigitizationTool::processBunchXing	(	int	bunchXing,
		SubEventIterator	bSubEvents,
		SubEventIterator	eSubEvents
	)

virtual

Reimplemented from PileUpToolBase.

Definition at line 163 of file PixelFastDigitizationTool.cxx.

 {
   //decide if this event will be processed depending on HardScatterSplittingMode & bunchXing
   if (m_HardScatterSplittingMode == 2 && !m_HardScatterSplittingSkipper ) { m_HardScatterSplittingSkipper = true; return StatusCode::SUCCESS; }
   if (m_HardScatterSplittingMode == 1 && m_HardScatterSplittingSkipper )  { return StatusCode::SUCCESS; }
   if (m_HardScatterSplittingMode == 1 && !m_HardScatterSplittingSkipper ) { m_HardScatterSplittingSkipper = true; }
  
   using TimedHitCollList = PileUpMergeSvc::TimedList<SiHitCollection>::type;
   TimedHitCollList hitCollList;
  
   if (!(m_mergeSvc->retrieveSubSetEvtData(m_inputObjectName.value(), hitCollList, bunchXing,
                                           bSubEvents, eSubEvents).isSuccess()) &&
       hitCollList.empty()) {
     ATH_MSG_ERROR("Could not fill TimedHitCollList");
     return StatusCode::FAILURE;
   } else {
     ATH_MSG_VERBOSE(hitCollList.size() << " SiHitCollections with key " <<
                     m_inputObjectName << " found");
   }
  
   TimedHitCollList::iterator iColl(hitCollList.begin());
   TimedHitCollList::iterator endColl(hitCollList.end());
  
   for( ; iColl != endColl; ++iColl) {
     SiHitCollection *siHitColl = new SiHitCollection(*iColl->second);
     PileUpTimeEventIndex timeIndex(iColl->first);
     ATH_MSG_DEBUG("SiHitCollection found with " << siHitColl->size() <<
                   " hits");
     ATH_MSG_VERBOSE("time index info. time: " << timeIndex.time()
                     << " index: " << timeIndex.index()
                     << " type: " << timeIndex.type());
     m_thpcsi->insert(timeIndex, siHitColl);
     m_siHitCollList.push_back(siHitColl);
   }
  
   return StatusCode::SUCCESS;
 }

◆ 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; }

◆ 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.

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

Member Data Documentation

◆ m_acceptDiagonalClusters

bool PixelFastDigitizationTool::m_acceptDiagonalClusters {true}

private

merging parameter used to define two clusters as neighbour >

Definition at line 140 of file PixelFastDigitizationTool.h.

◆ m_ambiguitiesMap

InDet::PixelGangedClusterAmbiguities* PixelFastDigitizationTool::m_ambiguitiesMap {}

private

Definition at line 142 of file PixelFastDigitizationTool.h.

◆ m_chargeDataKey

SG::ReadCondHandleKey<PixelChargeCalibCondData> PixelFastDigitizationTool::m_chargeDataKey {this, "PixelChargeCalibCondData", "PixelChargeCalibCondData", "Pixel charge calibration data"}

private

Definition at line 146 of file PixelFastDigitizationTool.h.

◆ m_clusterMaker

PublicToolHandle<InDet::ClusterMakerTool> PixelFastDigitizationTool::m_clusterMaker {this, "ClusterMaker", "InDet::ClusterMakerTool/FatrasClusterMaker"}

private

ToolHandle to ClusterMaker.

Definition at line 99 of file PixelFastDigitizationTool.h.

◆ m_digitizationStepper

PublicToolHandle<Trk::IModuleStepper> PixelFastDigitizationTool::m_digitizationStepper {this, "DigitizationStepper", "Trk::PlanarModuleStepper"}

private

Definition at line 164 of file PixelFastDigitizationTool.h.

◆ m_distortionKey

SG::ReadCondHandleKey<PixelDistortionData> PixelFastDigitizationTool::m_distortionKey {this, "PixelDistortionData", "PixelDistortionData", "Output readout distortion data"}

private

Definition at line 151 of file PixelFastDigitizationTool.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_gangedAmbiguitiesFinder

PublicToolHandle< InDet::PixelGangedAmbiguitiesFinder > PixelFastDigitizationTool::m_gangedAmbiguitiesFinder {this, "gangedAmbiguitiesFinder", "InDet::PixelGangedAmbiguitiesFinder"}

private

Definition at line 116 of file PixelFastDigitizationTool.h.

◆ m_HardScatterSplittingMode

IntegerProperty PixelFastDigitizationTool::m_HardScatterSplittingMode {this, "HardScatterSplittingMode", 0, "Control pileup & signal splitting"}

private

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

Definition at line 107 of file PixelFastDigitizationTool.h.

◆ m_HardScatterSplittingSkipper

bool PixelFastDigitizationTool::m_HardScatterSplittingSkipper {false}

private

Definition at line 108 of file PixelFastDigitizationTool.h.

◆ m_inputObjectName

StringProperty PixelFastDigitizationTool::m_inputObjectName {this, "InputObjectName", "PixelHits"}

private

Definition at line 118 of file PixelFastDigitizationTool.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_lorentzAngleTool

ToolHandle<ISiLorentzAngleTool> PixelFastDigitizationTool::m_lorentzAngleTool {this, "LorentzAngleTool", "PixelLorentzAngleTool", "Tool to retreive Lorentz angle"}

private

Definition at line 100 of file PixelFastDigitizationTool.h.

◆ m_mergeCluster

bool PixelFastDigitizationTool::m_mergeCluster {true}

private

enable the merging of neighbour Pixel clusters >

Definition at line 138 of file PixelFastDigitizationTool.h.

◆ m_mergeSvc

ServiceHandle<PileUpMergeSvc> PixelFastDigitizationTool::m_mergeSvc {this, "MergeSvc", "PileUpMergeSvc"}

private

PileUp Merge service.

Definition at line 106 of file PixelFastDigitizationTool.h.

◆ m_offlineCalibDataKey

SG::ReadCondHandleKey<PixelCalib::PixelOfflineCalibData> PixelFastDigitizationTool::m_offlineCalibDataKey {this, "PixelOfflineCalibData", "PixelOfflineCalibData", "Pixel offline calibration data"}

private

Definition at line 149 of file PixelFastDigitizationTool.h.

◆ m_pixDiffShiftBarrX

DoubleProperty PixelFastDigitizationTool::m_pixDiffShiftBarrX {this, "PixDiffShiftBarrX", 0.005}

private

Definition at line 132 of file PixelFastDigitizationTool.h.

◆ m_pixDiffShiftBarrY

DoubleProperty PixelFastDigitizationTool::m_pixDiffShiftBarrY {this, "PixDiffShiftBarrY", 0.005}

private

Definition at line 133 of file PixelFastDigitizationTool.h.

◆ m_pixDiffShiftEndCX

DoubleProperty PixelFastDigitizationTool::m_pixDiffShiftEndCX {this, "PixDiffShiftEndCX", 0.008}

private

Definition at line 134 of file PixelFastDigitizationTool.h.

◆ m_pixDiffShiftEndCY

DoubleProperty PixelFastDigitizationTool::m_pixDiffShiftEndCY {this, "PixDiffShiftEndCY", 0.008}

private

Definition at line 135 of file PixelFastDigitizationTool.h.

◆ m_pixel_ID

const PixelID* PixelFastDigitizationTool::m_pixel_ID {}

private

Handle to the ID helper.

Definition at line 97 of file PixelFastDigitizationTool.h.

◆ m_pixel_SiClustersName

StringProperty PixelFastDigitizationTool::m_pixel_SiClustersName {this, "PixelClusterContainerName", "PixelClusters"}

private

Definition at line 104 of file PixelFastDigitizationTool.h.

◆ m_pixelClusterAmbiguitiesMapName

StringProperty PixelFastDigitizationTool::m_pixelClusterAmbiguitiesMapName {this, "PixelClusterAmbiguitiesMapName", "PixelClusterAmbiguitiesMap"}

private

Definition at line 141 of file PixelFastDigitizationTool.h.

◆ m_pixelClusterContainer

InDet::PixelClusterContainer* PixelFastDigitizationTool::m_pixelClusterContainer {}

private

the PixelClusterContainer

Definition at line 103 of file PixelFastDigitizationTool.h.

◆ m_pixelClusterMap

Pixel_detElement_RIO_map* PixelFastDigitizationTool::m_pixelClusterMap {}

private

Definition at line 111 of file PixelFastDigitizationTool.h.

◆ m_pixelDetEleCollKey

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

private

Definition at line 154 of file PixelFastDigitizationTool.h.

◆ m_pixelReadout

ServiceHandle<InDetDD::IPixelReadoutManager> PixelFastDigitizationTool::m_pixelReadout {this, "PixelReadoutManager", "PixelReadoutManager", "Pixel readout manager" }

private

Definition at line 143 of file PixelFastDigitizationTool.h.

◆ m_pixEmulateSurfaceCharge

BooleanProperty PixelFastDigitizationTool::m_pixEmulateSurfaceCharge {this, "PixelEmulateSurfaceCharge", true}

private

emulate the surface charge

Definition at line 123 of file PixelFastDigitizationTool.h.

◆ m_pixErrorStrategy

IntegerProperty PixelFastDigitizationTool::m_pixErrorStrategy {this, "PixelErrorStrategy", 2}

private

error strategy for the ClusterMaker

Definition at line 131 of file PixelFastDigitizationTool.h.

◆ m_pixEtaError

DoubleArrayProperty PixelFastDigitizationTool::m_pixEtaError {this, "PixelErrorEta", {} }

private

eta error when not using the ClusterMaker

Definition at line 130 of file PixelFastDigitizationTool.h.

◆ m_pixMinimalPathCut

DoubleProperty PixelFastDigitizationTool::m_pixMinimalPathCut {this, "PixelMinimalPathLength", 0.06}

private

the 1.

model parameter: minimal 3D path in pixel

Definition at line 126 of file PixelFastDigitizationTool.h.

◆ m_pixModuleDistortion

BooleanProperty PixelFastDigitizationTool::m_pixModuleDistortion {this, "PixelEmulateModuleDistortion", true}

private

simulationn of module bowing

Definition at line 128 of file PixelFastDigitizationTool.h.

◆ m_pixPathLengthTotConv

DoubleProperty PixelFastDigitizationTool::m_pixPathLengthTotConv {this, "PixelPathLengthTotConversion", 125.}

private

from path length to tot

Definition at line 127 of file PixelFastDigitizationTool.h.

◆ m_pixPhiError

DoubleArrayProperty PixelFastDigitizationTool::m_pixPhiError {this, "PixelErrorPhi", {} }

private

phi error when not using the ClusterMaker

Definition at line 129 of file PixelFastDigitizationTool.h.

◆ m_pixPrdTruth

PRD_MultiTruthCollection* PixelFastDigitizationTool::m_pixPrdTruth {}

private

the PRD truth map for Pixel measurements

Definition at line 114 of file PixelFastDigitizationTool.h.

◆ m_pixSmearLandau

BooleanProperty PixelFastDigitizationTool::m_pixSmearLandau {this, "PixelSmearLandau", true}

private

if true : landau else: gauss

Definition at line 125 of file PixelFastDigitizationTool.h.

◆ m_pixSmearPathLength

DoubleProperty PixelFastDigitizationTool::m_pixSmearPathLength {this, "PixelSmearPathSigma", 0.01}

private

the 2.

model parameter: smear the path

Definition at line 124 of file PixelFastDigitizationTool.h.

◆ m_pixTanLorentzAngleScalor

double PixelFastDigitizationTool::m_pixTanLorentzAngleScalor {1.}

private

scale the lorentz angle effect

Definition at line 122 of file PixelFastDigitizationTool.h.

◆ m_pixUseClusterMaker

BooleanProperty PixelFastDigitizationTool::m_pixUseClusterMaker {this, "PixelUseClusterMaker", true}

private

use the pixel cluster maker or not

Definition at line 101 of file PixelFastDigitizationTool.h.

◆ m_prdTruthNamePixel

StringProperty PixelFastDigitizationTool::m_prdTruthNamePixel {this, "TruthNamePixel", "PRD_MultiTruthPixel"}

private

Definition at line 113 of file PixelFastDigitizationTool.h.

◆ m_randomEngineName

StringProperty PixelFastDigitizationTool::m_randomEngineName {this, "RndmEngine", "FastPixelDigitization"}

private

Name of the random number stream.

Definition at line 95 of file PixelFastDigitizationTool.h.

◆ m_rndmSvc

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

private

Random number service.

Definition at line 93 of file PixelFastDigitizationTool.h.

◆ m_siHitCollList

std::vector<SiHitCollection*> PixelFastDigitizationTool::m_siHitCollList {}

private

name of the sub event hit collections.

Definition at line 120 of file PixelFastDigitizationTool.h.

◆ m_splitClusters

short PixelFastDigitizationTool::m_splitClusters {0}

private

merging parameter used to define two clusters as neighbour >

Definition at line 139 of file PixelFastDigitizationTool.h.

◆ m_thpcsi

TimedHitCollection<SiHit>* PixelFastDigitizationTool::m_thpcsi {}

private

Definition at line 91 of file PixelFastDigitizationTool.h.

◆ m_ThrConverted

DoubleProperty PixelFastDigitizationTool::m_ThrConverted {this, "ThrConverted", 50000}

private

Definition at line 136 of file PixelFastDigitizationTool.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.

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

Public Member Functions

Private Types

Private Member Functions

Static Private Member Functions

Private Attributes

structors and AlgTool implementation

Detailed Description

Member Typedef Documentation

◆ Pixel_detElement_RIO_map

Constructor & Destructor Documentation

◆ PixelFastDigitizationTool() [1/3]

◆ ~PixelFastDigitizationTool()

◆ PixelFastDigitizationTool() [2/3]

◆ PixelFastDigitizationTool() [3/3]

Member Function Documentation

◆ areNeighbours()

◆ buildDetectorModule()

◆ CalculateIntersection()

◆ createAndStoreRIOs()

◆ Diffuse()

◆ digitize()

◆ filterPassed()

◆ initialize()

◆ mergeEvent()

◆ operator=()

◆ prepareEvent()

◆ processAllSubEvents() [1/2]

◆ processAllSubEvents() [2/2]

◆ processBunchXing()

◆ resetFilter()

◆ toProcess()

Member Data Documentation

◆ m_acceptDiagonalClusters

◆ m_ambiguitiesMap

◆ m_chargeDataKey

◆ m_clusterMaker

◆ m_digitizationStepper

◆ m_distortionKey

◆ m_filterPassed

◆ m_firstXing

◆ m_gangedAmbiguitiesFinder

◆ m_HardScatterSplittingMode

◆ m_HardScatterSplittingSkipper

◆ m_inputObjectName

◆ m_lastXing

◆ m_lorentzAngleTool

◆ m_mergeCluster

◆ m_mergeSvc

◆ m_offlineCalibDataKey

◆ m_pixDiffShiftBarrX

◆ m_pixDiffShiftBarrY

◆ m_pixDiffShiftEndCX

◆ m_pixDiffShiftEndCY

◆ m_pixel_ID

◆ m_pixel_SiClustersName

◆ m_pixelClusterAmbiguitiesMapName

◆ m_pixelClusterContainer

◆ m_pixelClusterMap

◆ m_pixelDetEleCollKey

◆ m_pixelReadout

◆ m_pixEmulateSurfaceCharge

◆ m_pixErrorStrategy

◆ m_pixEtaError

◆ m_pixMinimalPathCut

◆ m_pixModuleDistortion

◆ m_pixPathLengthTotConv

◆ m_pixPhiError

◆ m_pixPrdTruth

◆ m_pixSmearLandau

◆ m_pixSmearPathLength

◆ m_pixTanLorentzAngleScalor

◆ m_pixUseClusterMaker

◆ m_prdTruthNamePixel

◆ m_randomEngineName

◆ m_rndmSvc

◆ m_siHitCollList

◆ m_splitClusters

◆ m_thpcsi

◆ m_ThrConverted

◆ m_vetoPileUpTruthLinks