#include <HitCreatorSilicon.h>

Inheritance diagram for iFatras::HitCreatorSilicon:

Collaboration diagram for iFatras::HitCreatorSilicon:

Public Member Functions
	HitCreatorSilicon (const std::string &, const std::string &, const IInterface *)
	Constructor. More...

virtual	~HitCreatorSilicon ()=default
	Destructor. More...

StatusCode	initialize ()
	AlgTool initailize method. More...

StatusCode	finalize ()
	AlgTool finalize method. More...

void	handle (const Incident &inc)
	handle for incident service More...

void	createSimHit (const ISF::ISFParticle &isp, const Trk::TrackParameters &, double) const
	Return nothing - store the HIT in hit collection. More...

void	createSimHit (const ISF::ISFParticle &isp, const Trk::TrackParameters &pars, double time, const InDetDD::SiDetectorElement &hitSiDetElement, bool isSiDetElement) const
	templated function Return nothing - store the HIT in hit collection More...

const ParametersROT *	createHit (const ISF::ISFParticle &, const Trk::TrackParameters &) const
	Return the cluster on Track – the PrepRawData is contained in this one. More...

const std::vector< ParametersROT > *	createHits (const ISF::ISFParticle &, const ParametersLayer &) const
	Return the cluster on Track – the PrepRawData is contained in this one. More...

Protected Member Functions
double	energyDeposit_fast (const ISF::ISFParticle &isp, bool &isPix, bool &isSCT) const
	Calculate Energyloss with simple Landau approximation. More...

double	energyDeposit_exact (const ISF::ISFParticle &isp, bool &isPix, bool &isSCT) const
	Calculate Energyloss with exact Landau*Gauss. More...

Protected Attributes
ServiceHandle< IIncidentSvc >	m_incidentSvc {this, "IncidentService", "IncidentSvc"}

SiHitCollection *	m_hitColl {}
	the SiHit collection More...

StringProperty	m_collectionName {this, "CollectionName", "PixelHits"}
	name of the collection on storegate More...

ServiceHandle< IAtRndmGenSvc >	m_randomSvc {this, "RandomNumberService", "AtRndmGenSvc"}
	Pointer to the random number generator service. More...

StringProperty	m_randomEngineName {this, "RandomStreamName", "FatrasRnd"}
	Name of the random number stream. More...

CLHEP::HepRandomEngine *	m_randomEngine {}
	Random Engine. More...

StringProperty	m_siIdHelperName {this, "IdHelperName", "PixelID"}
	where to find the Si helper More...

const PixelID *	m_pixIdHelper {}
	the Pixel ID helper More...

const SCT_ID *	m_sctIdHelper {}
	the SCT ID helper More...

ToolHandle< IInDetConditionsTool >	m_condSummaryTool {this, "ConditionsTool", "PixelConditionsSummaryTool"}
	ToolHandle to ClusterMaker. More...

BooleanProperty	m_useConditionsTool {this, "UseConditionsTool", true}

TF1 *	m_dEdX_function {}
	function to evaluate dEdx More...

DoubleProperty	m_siPathToCharge {this, "PathToChargeConversion", 500.}
	convert path in silicon to charge More...

BooleanProperty	m_fastEnergyDepositionModel {this, "FastEnergyDepositionModel", true}
	use fast energy deposition model (landau approximation ) More...

Detailed Description

RIOOnTrack creation, starting from intersection on an active surface

Author: Andreas.Salzburger -at- cern.ch

Definition at line 50 of file HitCreatorSilicon.h.

Constructor & Destructor Documentation

◆ HitCreatorSilicon()

iFatras::HitCreatorSilicon::HitCreatorSilicon	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p
	)

Constructor.

Definition at line 47 of file HitCreatorSilicon.cxx.

                                                                                                          :
   base_class(t,n,p)
 {
 }

◆ ~HitCreatorSilicon()

virtual iFatras::HitCreatorSilicon::~HitCreatorSilicon ( )

virtualdefault

Destructor.

Member Function Documentation

◆ createHit()

const ParametersROT* iFatras::HitCreatorSilicon::createHit	(	const ISF::ISFParticle &	,
		const Trk::TrackParameters &
	)		const

inline

Return the cluster on Track – the PrepRawData is contained in this one.

Definition at line 76 of file HitCreatorSilicon.h.

76 { return 0; }

◆ createHits()

const std::vector< ParametersROT >* iFatras::HitCreatorSilicon::createHits	(	const ISF::ISFParticle &	,
		const ParametersLayer &
	)		const

inline

Return the cluster on Track – the PrepRawData is contained in this one.

Definition at line 79 of file HitCreatorSilicon.h.

79 { return 0; }

◆ createSimHit() [1/2]

void iFatras::HitCreatorSilicon::createSimHit	(	const ISF::ISFParticle &	isp,
		const Trk::TrackParameters &	pars,
		double	time
	)		const

Return nothing - store the HIT in hit collection.

Definition at line 326 of file HitCreatorSilicon.cxx.

                                                                                                                          {
  
   // get surface and DetElement base
   const Trk::Surface &hitSurface = pars.associatedSurface();
   const Trk::TrkDetElementBase* detElementBase = hitSurface.associatedDetectorElement();
  
   // the detector element cast -> needed
   const InDetDD::SiDetectorElement* SiDetElement = dynamic_cast<const InDetDD::SiDetectorElement*>((detElementBase));
  
   // return triggered if no siDetElement or no current particle link to the stack
   if ( !SiDetElement ){
     ATH_MSG_WARNING("[ sihit ] No Silicon Detector element available. Ignore this one.");
     return;
   }
  
   return createSimHit( isp, pars, time, *SiDetElement, 1);
  
 }

◆ createSimHit() [2/2]

void iFatras::HitCreatorSilicon::createSimHit	(	const ISF::ISFParticle &	isp,
		const Trk::TrackParameters &	pars,
		double	time,
		const InDetDD::SiDetectorElement &	hitSiDetElement,
		bool	isSiDetElement
	)		const

templated function Return nothing - store the HIT in hit collection

<

Todo:: : fix edge effects

Definition at line 345 of file HitCreatorSilicon.cxx.

                                                                                                                                                                                               {
  
   // get surface and DetElement base
   const Trk::Surface &hitSurface = pars.associatedSurface();
  
   // get the identifier and hash identifier
   Identifier hitId         = hitSurface.associatedDetectorElementIdentifier();
   IdentifierHash hitIdHash = hitSiDetElement.identifyHash();
   // check conditions of the intersection
   if ( m_useConditionsTool ) {
     const EventContext& ctx = Gaudi::Hive::currentContext();
     bool isActive = m_condSummaryTool->isActive(hitIdHash, hitId, ctx);                   // active = "element returns data"
     bool isGood   = isActive ? m_condSummaryTool->isGood(hitIdHash, hitId, ctx) : false;  // good   = "data are reliable"
     if (!isActive)
       ATH_MSG_VERBOSE("[ sihit ]  ID " << hitId << ", hash " << hitIdHash << " is not active. ");
     else if (!isGood)
       ATH_MSG_VERBOSE("[ sihit ]  ID " << hitId << ", hash " << hitIdHash << " is active but not good. ");
     else
       ATH_MSG_VERBOSE("[ sihit ]  ID " << hitId << ", hash " << hitIdHash << " is active and good.");
     if (!isActive || !isGood) return;
   }
  
   // intersection
   const Amg::Vector2D& intersection = pars.localPosition();
   double interX = intersection.x();
   double interY = intersection.y();
   // thickness of the module
   const double thickness = hitSiDetElement.thickness();
   // get the momentum direction into the local frame
   Amg::Vector3D particleDir = pars.momentum().unit();
   const Amg::Transform3D& sTransform = hitSurface.transform();
   Amg::Vector3D localDirection = sTransform.inverse().linear() * particleDir;
   localDirection *= thickness/cos(localDirection.theta());
   // moving direction
   int movingDirection = localDirection.z() > 0. ? 1 : -1;
   // get he local distance of the intersection in x,y
   double distX = localDirection.x();
   double distY = localDirection.y();
   // local entries in x,y
   double localEntryX = interX-0.5*distX;
   double localEntryY = interY-0.5*distY;
   double localExitX  = interX+0.5*distX;
   double localExitY  = interY+0.5*distY;
   const Amg::Transform3D &hitTransform = hitSiDetElement.transformHit().inverse();
   // transform into the hit frame
   Amg::Vector3D localEntry(hitTransform*(sTransform*Amg::Vector3D(localEntryX,localEntryY,-0.5*movingDirection*thickness)));
   Amg::Vector3D localExit(hitTransform*(sTransform*Amg::Vector3D(localExitX,localExitY,0.5*movingDirection*thickness)));
  
   bool isPix = hitSiDetElement.isPixel();
   bool isSCT = hitSiDetElement.isSCT();
  
   // Landau approximation
   const double dEdX = m_fastEnergyDepositionModel ?
     energyDeposit_fast(isp,isPix,isSCT)
     : energyDeposit_exact(isp,isPix,isSCT);
  
   const double energyDeposit = dEdX * (localExit - localEntry).mag();
  
   // create the silicon hit
   const HepGeom::Point3D<double> localEntryHep( localEntry.x(), localEntry.y(), localEntry.z() );
   const HepGeom::Point3D<double> localExitHep( localExit.x(), localExit.y(), localExit.z() );
  
   if ( isSiDetElement )
     ATH_MSG_VERBOSE("[ sihit ] Adding an SiHit SiDetElement to the SiHitCollection.");
   else
     ATH_MSG_VERBOSE("[ sihit ] SiHit SiDetElement not found to add to the SiHitCollection.");
  
   HepMcParticleLink partLink(HepMC::uniqueID(isp), 0,
                              HepMcParticleLink::IS_POSITION,
                              HepMcParticleLink::IS_ID);
   m_hitColl->Emplace(localEntryHep,
                      localExitHep,
                      energyDeposit,
                      time,
                      partLink,
                      m_pixIdHelper ? 0 : 1,
                      m_pixIdHelper ? m_pixIdHelper->barrel_ec(hitId)  : m_sctIdHelper->barrel_ec(hitId),
                      m_pixIdHelper ? m_pixIdHelper->layer_disk(hitId) : m_sctIdHelper->layer_disk(hitId),
                      m_pixIdHelper ? m_pixIdHelper->eta_module(hitId) : m_sctIdHelper->eta_module(hitId),
                      m_pixIdHelper ? m_pixIdHelper->phi_module(hitId) : m_sctIdHelper->phi_module(hitId),
                      m_pixIdHelper ? 0 : m_sctIdHelper->side(hitId));
  
 }

◆ energyDeposit_exact()

double iFatras::HitCreatorSilicon::energyDeposit_exact	(	const ISF::ISFParticle &	isp,
		bool &	isPix,
		bool &	isSCT
	)		const

protected

Calculate Energyloss with exact Landau*Gauss.

Definition at line 245 of file HitCreatorSilicon.cxx.

                                                                                                                {
   int pdg_id= isp.pdgCode();
   Amg::Vector3D mom = isp.momentum();
   double Momentum = sqrt((mom.x()*mom.x()) +( mom.y()*mom.y())+(mom.z()*mom.z()));
   m_dEdX_function->SetNpx(4);
   double xmin = 0;
   double xmax = 0.;
   if(abs(pdg_id)==11){ // electrons
     m_dEdX_function->SetParameters(0.009807, 0.3 ,453.7,0.03733);
     xmax = 1.;
   }
   else{ //  Kaonen, Protonen, Pionen, Myonen
  
     double rest_mass= isp.mass();
     double para0 = -1.12409e-02;
     double para1 = 1.42786e-11;
     double beta = Momentum/ sqrt((Momentum*Momentum)+(rest_mass*rest_mass));
     double MPV = (para0/pow(beta,2))*(log(para1*(pow(beta,2)/(1-pow(beta,2))))-pow(beta,2));
     if(Momentum < 1000){ // momentum < 1GeV
       if(abs(pdg_id)==211){  // Pionen
         if(isSCT){
           m_dEdX_function->SetParameters(0.0250, MPV ,453.7,0.0109);
         }
         if(isPix){
           m_dEdX_function->SetParameters(0.0205, MPV ,453.7,0.0219);
         }
       }
       else if(abs(pdg_id)==321){ //Kaonen
         if(isSCT){
           m_dEdX_function->SetParameters(0.0458, MPV ,453.7,0.0064);
         }
         if(isPix){
           m_dEdX_function->SetParameters(0.0465, MPV,453.7,0.000);
         }
       }
       else if(abs(pdg_id)==13){ // Muonen
         if(isSCT){
           m_dEdX_function->SetParameters(0.0140, MPV,453.7,0.0245);
         }
         if(isPix){
           m_dEdX_function->SetParameters(0.0064, MPV ,453.7,0.0480);
         }
       }
       else if(abs(pdg_id)==2212){ // Protonen
         if(isSCT){
           m_dEdX_function->SetParameters(0.0458, MPV ,453.7,0.0064);
         }
         if(isPix){
           m_dEdX_function->SetParameters(0.0465, MPV ,453.7,0.000);
         }
       }
  
     }
     else if(Momentum>= 1000){
       if(isSCT){
         m_dEdX_function->SetParameters(0.0135, MPV ,453.7,0.0245);
       }
       if(isPix){
         m_dEdX_function->SetParameters(0.0102, MPV ,453.7,0.0297);
       }
     }
     xmax = MPV*10.;
   } //  Parametrisierung fuer Kaonen, Protonen, Pionen, Myonen ENDE
  
   double ymin = 0;
   double ymax = m_dEdX_function -> GetMaximum(xmin, xmax);
   double energyLoss = 0;
  
   for (int count=0; count <1000; count++) {
     double ry = CLHEP::RandFlat::shoot(m_randomEngine)*(ymax-ymin) + ymin;
     double rx = CLHEP::RandFlat::shoot(m_randomEngine)*(xmax-xmin) + xmin;
     double y = m_dEdX_function->Eval(rx);
     if ( ry <= y ){ energyLoss = rx; break;}
   }
  
  
  
   return energyLoss;
 }

◆ energyDeposit_fast()

double iFatras::HitCreatorSilicon::energyDeposit_fast	(	const ISF::ISFParticle &	isp,
		bool &	isPix,
		bool &	isSCT
	)		const

protected

Calculate Energyloss with simple Landau approximation.

Definition at line 177 of file HitCreatorSilicon.cxx.

                                                                                                               {
  
   int pdg_id= isp.pdgCode();
   Amg::Vector3D mom = isp.momentum();
   double Momentum = sqrt((mom.x()*mom.x()) +( mom.y()*mom.y())+(mom.z()*mom.z()));
   double randLand = CLHEP::RandLandau::shoot(m_randomEngine);
   double energyLoss=0;
   if(abs(pdg_id)==11){ // electrons
     energyLoss = (randLand*0.032)+0.32;
   }
   else{ // kaons, protons, pions, myons
     double rest_mass= isp.mass();
     double para0 = -1.12409e-02;
     double para1 = 1.42786e-11;
     // ST : TEMPORARY remove beta dependence ( fixed to 1 GeV )
     Momentum = 1000.;
     // ST
     double beta = Momentum/ sqrt((Momentum*Momentum)+(rest_mass*rest_mass));
     double MPV = (para0/pow(beta,2))*(log(para1*(pow(beta,2)/(1-pow(beta,2))))-pow(beta,2));
  
     if( Momentum< 1000){ //momentum < 1GeV
       if(abs(pdg_id)==211){  // Pionen
         if(isSCT){
           energyLoss = (randLand*0.023)+(MPV+0.008);
         }
         if(isPix){
           energyLoss = (randLand*0.0225)+(MPV+0.008);
         }
       }
       else if(abs(pdg_id)==321){ //Kaonen
         if(isSCT){
           energyLoss = (randLand*0.048)+(MPV+0.01);
         }
         if(isPix){
           energyLoss = (randLand*0.0465)+(MPV+0.01);
         }
       }
       else if(abs(pdg_id)==13){ // Muonen
         if(isSCT){
           energyLoss = (randLand*0.025)+(MPV+0.011);
         }
         if(isPix){
           energyLoss = (randLand*0.048)+(MPV+0.011);
         }
       }
       else if(abs(pdg_id)==2212){ // Protonen
         if(isSCT){
           energyLoss = (randLand*0.0458)+(MPV+0.008);
         }
         if(isPix){
           energyLoss = (randLand*0.0465)+(MPV+0.008);
         }
       }
  
     }//end momentum < 1GeV
     else if(Momentum >= 1000){
       if(isSCT){
         energyLoss = (randLand*0.025)+(MPV+0.013);
       }
       if(isPix){
         energyLoss = (randLand*0.025)+(MPV+0.014);
       }
     }//end Momentum >= 1GeV
   } // End kaons, protons, pions, myons
  
   return energyLoss;
 }

◆ finalize()

StatusCode iFatras::HitCreatorSilicon::finalize ( )

AlgTool finalize method.

Definition at line 87 of file HitCreatorSilicon.cxx.

 {
   delete m_dEdX_function;
   return StatusCode::SUCCESS;
 }

◆ handle()

void iFatras::HitCreatorSilicon::handle ( const Incident & inc )

handle for incident service

Definition at line 94 of file HitCreatorSilicon.cxx.

                                                            {
   // check the incident type
   if ( inc.type() == IncidentType::BeginEvent ){
     // check if the hit collection already contains:
     // (a)     if yes ... try to retrieve it
     if ( evtStore()->contains<SiHitCollection>(m_collectionName) ){
       if ( (evtStore()->retrieve(m_hitColl , m_collectionName)).isFailure() )
         ATH_MSG_ERROR( "[ --- ] Unable to retrieve SiHitCollection " << m_collectionName);
       // (b)     if no ... try to create it
     } else {
       m_hitColl = new SiHitCollection( m_collectionName);
       if ( (evtStore()->record(m_hitColl, m_collectionName, true)).isFailure() ) {
         ATH_MSG_ERROR( "[ --- ] Unable to record SiHitCollection " << m_collectionName);
         delete m_hitColl; m_hitColl=0;
       }
     }
   }
   return;
 }

◆ initialize()

StatusCode iFatras::HitCreatorSilicon::initialize ( )

AlgTool initailize method.

Definition at line 54 of file HitCreatorSilicon.cxx.

 {
   // Random number service
   ATH_CHECK ( m_randomSvc.retrieve() );
  
   //Get own engine with own seeds:
   m_randomEngine = m_randomSvc->GetEngine(m_randomEngineName);
   if (!m_randomEngine) {
     ATH_MSG_ERROR( "[ --- ] Could not get random engine '" << m_randomEngineName << "'" );
     return StatusCode::FAILURE;
   }
  
   ATH_CHECK ( m_condSummaryTool.retrieve( DisableTool{ !m_useConditionsTool } ) );
  
   // Get the Pixel Identifier-helper:
   if ( m_siIdHelperName == "PixelID" ) {
     ATH_CHECK ( detStore()->retrieve(m_pixIdHelper, m_siIdHelperName) );
   }
   else if ( m_siIdHelperName == "SCT_ID" ) {
     ATH_CHECK ( detStore()->retrieve(m_sctIdHelper, m_siIdHelperName) );
   }
  
   // Athena/Gaudi framework
   ATH_CHECK (m_incidentSvc.retrieve());
  
   // register to the incident service: BeginEvent for TrackCollection
   m_incidentSvc->addListener( this, IncidentType::BeginEvent);
   m_dEdX_function = new TF1 ("fitfunc2", langaufun_fast, 0.,10.,4);
   ATH_MSG_VERBOSE( "[ sihit ]  initialize() successful." );
   return StatusCode::SUCCESS;
 }

Member Data Documentation

◆ m_collectionName

StringProperty iFatras::HitCreatorSilicon::m_collectionName {this, "CollectionName", "PixelHits"}

protected

name of the collection on storegate

Definition at line 88 of file HitCreatorSilicon.h.

◆ m_condSummaryTool

ToolHandle<IInDetConditionsTool> iFatras::HitCreatorSilicon::m_condSummaryTool {this, "ConditionsTool", "PixelConditionsSummaryTool"}

protected

ToolHandle to ClusterMaker.

Handle to Pixel/SCT conditions tool

Definition at line 100 of file HitCreatorSilicon.h.

◆ m_dEdX_function

TF1* iFatras::HitCreatorSilicon::m_dEdX_function {}

protected

function to evaluate dEdx

Definition at line 102 of file HitCreatorSilicon.h.

◆ m_fastEnergyDepositionModel

BooleanProperty iFatras::HitCreatorSilicon::m_fastEnergyDepositionModel {this, "FastEnergyDepositionModel", true}

protected

use fast energy deposition model (landau approximation )

Definition at line 105 of file HitCreatorSilicon.h.

◆ m_hitColl

SiHitCollection* iFatras::HitCreatorSilicon::m_hitColl {}

protected

the SiHit collection

Definition at line 87 of file HitCreatorSilicon.h.

◆ m_incidentSvc

ServiceHandle<IIncidentSvc> iFatras::HitCreatorSilicon::m_incidentSvc {this, "IncidentService", "IncidentSvc"}

protected

Definition at line 84 of file HitCreatorSilicon.h.

◆ m_pixIdHelper

const PixelID* iFatras::HitCreatorSilicon::m_pixIdHelper {}

protected

the Pixel ID helper

Definition at line 96 of file HitCreatorSilicon.h.

◆ m_randomEngine

CLHEP::HepRandomEngine* iFatras::HitCreatorSilicon::m_randomEngine {}

protected

Random Engine.

Definition at line 93 of file HitCreatorSilicon.h.

◆ m_randomEngineName

StringProperty iFatras::HitCreatorSilicon::m_randomEngineName {this, "RandomStreamName", "FatrasRnd"}

protected

Name of the random number stream.

Definition at line 92 of file HitCreatorSilicon.h.

◆ m_randomSvc

ServiceHandle<IAtRndmGenSvc> iFatras::HitCreatorSilicon::m_randomSvc {this, "RandomNumberService", "AtRndmGenSvc"}

protected

Pointer to the random number generator service.

Random Svc

Definition at line 91 of file HitCreatorSilicon.h.

◆ m_sctIdHelper

const SCT_ID* iFatras::HitCreatorSilicon::m_sctIdHelper {}

protected

the SCT ID helper

Definition at line 97 of file HitCreatorSilicon.h.

◆ m_siIdHelperName

StringProperty iFatras::HitCreatorSilicon::m_siIdHelperName {this, "IdHelperName", "PixelID"}

protected

where to find the Si helper

Definition at line 95 of file HitCreatorSilicon.h.

◆ m_siPathToCharge

DoubleProperty iFatras::HitCreatorSilicon::m_siPathToCharge {this, "PathToChargeConversion", 500.}

protected

convert path in silicon to charge

Definition at line 104 of file HitCreatorSilicon.h.

◆ m_useConditionsTool

BooleanProperty iFatras::HitCreatorSilicon::m_useConditionsTool {this, "UseConditionsTool", true}

protected

Definition at line 101 of file HitCreatorSilicon.h.

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

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ HitCreatorSilicon()

◆ ~HitCreatorSilicon()

Member Function Documentation

◆ createHit()

◆ createHits()

◆ createSimHit() [1/2]

◆ createSimHit() [2/2]

◆ energyDeposit_exact()

◆ energyDeposit_fast()

◆ finalize()

◆ handle()

◆ initialize()

Member Data Documentation

◆ m_collectionName

◆ m_condSummaryTool

◆ m_dEdX_function

◆ m_fastEnergyDepositionModel

◆ m_hitColl

◆ m_incidentSvc

◆ m_pixIdHelper

◆ m_randomEngine

◆ m_randomEngineName

◆ m_randomSvc

◆ m_sctIdHelper

◆ m_siIdHelperName

◆ m_siPathToCharge

◆ m_useConditionsTool