iFatras::HitCreatorSilicon Class Reference

RIOOnTrack creation, starting from intersection on an active surface. More...

#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.
virtual	~HitCreatorSilicon ()=default
	Destructor.
StatusCode	initialize ()
	AlgTool initailize method.
StatusCode	finalize ()
	AlgTool finalize method.
void	handle (const Incident &inc)
	handle for incident service
void	createSimHit (const ISF::ISFParticle &isp, const Trk::TrackParameters &, double) const
	Return nothing - store the HIT in hit collection.
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
const ParametersROT *	createHit (const ISF::ISFParticle &, const Trk::TrackParameters &) const
	Return the cluster on Track – the PrepRawData is contained in this one.
const std::vector< ParametersROT > *	createHits (const ISF::ISFParticle &, const ParametersLayer &) const
	Return the cluster on Track – the PrepRawData is contained in this one.

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

Protected Attributes
ServiceHandle< IIncidentSvc >	m_incidentSvc {this, "IncidentService", "IncidentSvc"}
SiHitCollection *	m_hitColl {}
	the SiHit collection
StringProperty	m_collectionName {this, "CollectionName", "PixelHits"}
	name of the collection on storegate
ServiceHandle< IAtRndmGenSvc >	m_randomSvc {this, "RandomNumberService", "AtRndmGenSvc"}
	Pointer to the random number generator service.
StringProperty	m_randomEngineName {this, "RandomStreamName", "FatrasRnd"}
	Name of the random number stream.
CLHEP::HepRandomEngine *	m_randomEngine {}
	Random Engine.
StringProperty	m_siIdHelperName {this, "IdHelperName", "PixelID"}
	where to find the Si helper
const PixelID *	m_pixIdHelper {}
	the Pixel ID helper
const SCT_ID *	m_sctIdHelper {}
	the SCT ID helper
ToolHandle< IInDetConditionsTool >	m_condSummaryTool {this, "ConditionsTool", "PixelConditionsSummaryTool"}
	ToolHandle to ClusterMaker.
BooleanProperty	m_useConditionsTool {this, "UseConditionsTool", true}
TF1 *	m_dEdX_function {}
	function to evaluate dEdx
DoubleProperty	m_siPathToCharge {this, "PathToChargeConversion", 500.}
	convert path in silicon to charge
BooleanProperty	m_fastEnergyDepositionModel {this, "FastEnergyDepositionModel", true}
	use fast energy deposition model (landau approximation )

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.

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

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

{this, "CollectionName", "PixelHits"};           

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.

{this, "ConditionsTool", "PixelConditionsSummaryTool"};          

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.

{this, "FastEnergyDepositionModel", true}; 

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.

{this, "IncidentService", "IncidentSvc"};

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.

{this, "RandomStreamName", "FatrasRnd"};         

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.

{this, "RandomNumberService", "AtRndmGenSvc"};                

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.

{this, "IdHelperName", "PixelID"};          

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.

{this, "PathToChargeConversion", 500.};           

m_useConditionsTool

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

protected

Definition at line 101 of file HitCreatorSilicon.h.

{this, "UseConditionsTool", true};

---

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

• HitCreatorSilicon.h
• HitCreatorSilicon.cxx