PndHypCalcStrip Class Reference

#include <PndHypCalcStrip.h>

Public Member Functions
	PndHypCalcStrip ()
	PndHypCalcStrip (Double_t pitch, Double_t orient, Int_t nrStrips, Int_t nrFeChannels, const TVector2 &firstStripAnchor, Double_t threshold, Double_t noise)
	PndHypCalcStrip (const PndHypStripDigiPar *digipar, SensorSide side=kTOP)
std::vector< PndHypStrip >	GetStrips (Double_t inx, Double_t iny, Double_t inz, Double_t outx, Double_t outy, Double_t outz, Double_t eLoss, int id)
Int_t	CalcFEfromStrip (Int_t stripNr) const
Int_t	CalcChannelfromStrip (Int_t stripNr) const
void	CalcFeChToStrip (Int_t fe, Int_t channel, Int_t &strip, enum SensorSide &side) const
void	CalcStripPointOnLine (const Double_t strip, TVector2 &point) const
const TVector2	GetStripDirection () const
void	SetVerboseLevel (Int_t level)
void	SetreStrip (Int_t reStp)
void	SetAnchor (TVector2 edge)
void	print () const

Private Member Functions
Double_t	CalcStripFromPoint (Double_t x, Double_t y)
Double_t	SmearCharge (Double_t charge)
Double_t	ChargeFromEloss (Double_t eloss) const

Private Attributes
Double_t	fPitch
Double_t	fOrient
	strip pitch (cm) More...
Int_t	fNrStrips
	strip orientation angle More...
Int_t	fNrFeChannels
	Nr. of strips on active area. More...
Double_t	fThreshold
	Nr of Channels per FE. More...
Double_t	fNoise
	charge threshold More...
TVector2	fAnchor
	ENC. More...
TVector2	fStripDir
	anchor point on first strip More...
TVector2	fOrthoDir
	vector perpendicular to strip direction More...
TRandom3 *	fRNG
	vector orthogonal to strip direction More...
Int_t	fVerboseLevel
	Random Number Generator. More...
Int_t	fsize
int	_id

Detailed Description

Definition at line 28 of file PndHypCalcStrip.h.

Constructor & Destructor Documentation

PndHypCalcStrip::PndHypCalcStrip

(

)

Definition at line 20 of file PndHypCalcStrip.cxx.

References fAnchor, fNoise, fNrStrips, fOrient, fPitch, fRNG, fThreshold, and fVerboseLevel.

                                {
  fPitch = 0.;
  fOrient = 0.;
  fAnchor = TVector2(0.,0.);
  fNrStrips = 0;
  fThreshold = 0.;
  fNoise = 0.;
  fVerboseLevel = 1;
  fRNG = new TRandom3();
}

PndHypCalcStrip::PndHypCalcStrip	(	Double_t	pitch,
		Double_t	orient,
		Int_t	nrStrips,
		Int_t	nrFeChannels,
		const TVector2 &	firstStripAnchor,
		Double_t	threshold,
		Double_t	noise
	)

Definition at line 31 of file PndHypCalcStrip.cxx.

References cos(), fOrient, fOrthoDir, fRNG, fStripDir, fVerboseLevel, and sin().

  : fPitch(pitch), fOrient(orient),
    fNrStrips(nrStrips), fNrFeChannels(nrFeChannels),
    fThreshold(threshold), fNoise(noise),
    fAnchor(firstStripAnchor)
{
  fStripDir.Set(cos(fOrient),sin(fOrient));
  fOrthoDir.Set(sin(fOrient),-cos(fOrient));
  fVerboseLevel = 1;
  fRNG = new TRandom3();
  //Print();
}

PndHypCalcStrip::PndHypCalcStrip	(	const PndHypStripDigiPar *	digipar,
		SensorSide	side = kTOP
	)

Definition at line 47 of file PndHypCalcStrip.cxx.

References cos(), fAnchor, fNoise, fNrFeChannels, fNrStrips, fOrient, fOrthoDir, fPitch, fRNG, fStripDir, fThreshold, fVerboseLevel, PndHypStripDigiPar::GetBotAnchor(), PndHypStripDigiPar::GetBotPitch(), PndHypStripDigiPar::GetNoise(), PndHypStripDigiPar::GetNrBotFE(), PndHypStripDigiPar::GetNrFECh(), PndHypStripDigiPar::GetNrTopFE(), PndHypStripDigiPar::GetOrient(), PndHypStripDigiPar::GetSkew(), PndHypStripDigiPar::GetThreshold(), PndHypStripDigiPar::GetTopAnchor(), PndHypStripDigiPar::GetTopPitch(), kBOTTOM, kTOP, and sin().

{
  if(side == kTOP)
    {
      fPitch = digipar->GetTopPitch();
      fOrient = digipar->GetOrient();
      fAnchor = digipar->GetTopAnchor();
      fNrStrips = digipar->GetNrTopFE()*digipar->GetNrFECh();
    }
  else if(side == kBOTTOM)
    {
      fPitch = digipar->GetBotPitch();
      fOrient = digipar->GetOrient() + digipar->GetSkew();
      fAnchor = digipar->GetBotAnchor();
      fNrStrips = digipar->GetNrBotFE()*digipar->GetNrFECh();
    }
  fNrFeChannels = digipar->GetNrFECh();
  fThreshold = digipar->GetThreshold();
  fNoise = digipar->GetNoise();
  
  fStripDir.Set(cos(fOrient),sin(fOrient));
  fOrthoDir.Set(sin(fOrient),-cos(fOrient));
  fVerboseLevel = 1;
  fRNG = new TRandom3();
  //Print();
  
}

Member Function Documentation

Int_t PndHypCalcStrip::CalcChannelfromStrip

(

Int_t

stripNr

)

const

84 * Calculate Channel index (on Frontend) from strip-index 85 *

Parameters

stripNr

strip index 86 *

Returns

channel on frontend 87

Definition at line 240 of file PndHypCalcStrip.cxx.

References fNrFeChannels.

Referenced by PndHypStripHitProducer::Exec().

{ return (stripNr%fNrFeChannels); }

void PndHypCalcStrip::CalcFeChToStrip	(	Int_t	fe,
		Int_t	channel,
		Int_t &	strip,
		enum SensorSide &	side
	)		const

91 * Calculate global strip index from Frontend number and -channel 92 *

Parameters

fe	frontend number 93 *
channel	frontend channel 94 *
strip	global strip index 95 *
side	96

Definition at line 242 of file PndHypCalcStrip.cxx.

References fNrFeChannels, fNrStrips, kBOTTOM, and kTOP.

Referenced by PndHypStripClusterTask::Exec().

  {
    Int_t nr = fe * fNrFeChannels + channel;
   //  if (nr > fNrStrips) {
//       nr -= fNrStrips;
//       side = kBOTTOM;
//     } else side = kTOP;

    if (nr < fNrStrips) {
      side = kTOP;
    } else {
      nr -= fNrStrips;
      side = kBOTTOM;
    }
    strip = nr;
    //std::cout<< "strip "<<strip<<"side "<<side<<" chan "<<channel<<std::endl;
    
  }

Int_t PndHypCalcStrip::CalcFEfromStrip

(

Int_t

stripNr

)

const

Calculate Frontend number from strip-index

Parameters

stripNr

strip index

Returns

Frontend number

Definition at line 238 of file PndHypCalcStrip.cxx.

References fNrFeChannels.

Referenced by PndHypStripHitProducer::Exec().

{ return (stripNr/fNrFeChannels); }

Double_t PndHypCalcStrip::CalcStripFromPoint ( Double_t  x, Double_t  y  )

private

125 * Calculate continuous strip-number parameter from given Point in local system 126 *

Parameters

x	x-coordinate 127 *
y	y-coordinate 128 *

Returns

strip-parameter nu 129

Definition at line 229 of file PndHypCalcStrip.cxx.

References fAnchor, fOrthoDir, and fPitch.

Referenced by GetStrips().

  {
    return ( (x-fAnchor.X())*fOrthoDir.Y()
             -(y-fAnchor.Y())*fOrthoDir.X() )/fPitch;
    //std::cout<<" nrstrip "
    //   <<( (x-fAnchor.X())*fStripDir.Y()-(y-fAnchor.Y())*fStripDir.X() )/fPitch<<std::endl;
    
  }

void PndHypCalcStrip::CalcStripPointOnLine	(	const Double_t	strip,
		TVector2 &	point
	)		const

Definition at line 262 of file PndHypCalcStrip.cxx.

References fAnchor, fOrthoDir, and fPitch.

Referenced by PndHypStripClusterTask::Exec().

{
  point =  fAnchor + ( fPitch * (strip+0.5) * fOrthoDir ) ;
}

Double_t PndHypCalcStrip::ChargeFromEloss ( Double_t  eloss ) const

inlineprivate

Definition at line 116 of file PndHypCalcStrip.h.

Referenced by GetStrips().

{return eloss/(3.61e-9);}

const TVector2 PndHypCalcStrip::GetStripDirection ( ) const

inline

Get Strip Direction (strip orientation angle)

Definition at line 72 of file PndHypCalcStrip.h.

References fStripDir.

Referenced by PndHypStripClusterTask::Exec().

{return fStripDir;}

std::vector< PndHypStrip > PndHypCalcStrip::GetStrips	(	Double_t	inx,
		Double_t	iny,
		Double_t	inz,
		Double_t	outx,
		Double_t	outy,
		Double_t	outz,
		Double_t	eLoss,
		int	id
	)

Definition at line 77 of file PndHypCalcStrip.cxx.

References CalcStripFromPoint(), ChargeFromEloss(), Double_t, fabs(), fNrStrips, fThreshold, fVerboseLevel, n, out, and SmearCharge().

Referenced by PndHypStripHitProducer::Exec().

{//1
  if (fVerboseLevel > 2) std::cout<<"-I-  PndHypCalcStrip::GetStrips "<<std::endl;
  
  // 2d-Projection of trajectory
  //mvd local c.s(x,y);hyp local c.s(x-z);3.03.08 new hyp cs(xy)
  TVector2 in(inx,iny);//TVector2 in(inx,inz);
  TVector2 out(outx,outy);//TVector2 out(outx,outz);
  TVector2 path=out-in;
  
  if (fVerboseLevel > 2){
    std::cout<<" InPoint: ("<<in.X()<<","<<in.Y()<<std::endl;
    std::cout<<" OutPoint: ("<<out.X()<<","<<out.Y()<<std::endl;
  }
  
  std::vector<PndHypStrip> strips;
  Double_t SmearedQ;
  
  if (path.Mod()<1E-18) {
    std::cout<<"-W-  PndHypCalcStrip::GetStrips : No Trajectory inside Sensor!"<<path.Mod()<<std::endl;
    return strips;
  }    
  
  if (fVerboseLevel > 1) std::cout<<" pathlength: "<<path.Mod()<<std::endl;
  
  Double_t nuIn  = CalcStripFromPoint(inx,iny);
  Double_t nuOut = CalcStripFromPoint(outx,outy);
  
  if (fVerboseLevel > 2) std::cout<<" nuIn = "<<nuIn<<" ; nuOut = "<<nuOut<<std::endl;
  
  //Double_t dir = 0.;
  //if (nuOut>nuIn) dir=1.;
  //else dir = -1.;
  
  Double_t Q = ChargeFromEloss(eLoss);//*1E9/3.61; // 3.6 Electrons/eV in Silicon
  if (fVerboseLevel > 1) std::cout<<" integral charge = "<<Q<<std::endl;
  
  // did we hit the active area ?
  //if ( (nuIn<0.5 && nuOut<0.5) || 
  //   (((nuIn+0.5) >Double_t(fNrStrips-1)) && ((nuOut+0.5) > Double_t(fNrStrips-1))))
   
  if ( (nuIn<0. && nuOut<0.) || ((nuIn >Double_t(fNrStrips)) && (nuOut > Double_t(fNrStrips))))
    {
      if (fVerboseLevel > 1) std::cout<<" Hit outside active area."<<std::endl;
      return strips;
    }
  
  // is the In-Point inside active area ?// only charge fraction inside active area taken
  if (nuIn<0.){
    Q *= (nuOut)/(nuOut-nuIn);
    nuIn = 0.;
  } else if (nuIn > (Double_t(fNrStrips))){
    Q *= ((Double_t)fNrStrips-nuOut)/(-nuOut+nuIn);
    nuIn = Double_t(fNrStrips);
  }
  
  // is the Out-Point inside active area ?
  if (nuOut<0.){
    Q *= (nuIn)/(-nuOut+nuIn);
    nuOut = 0.;
    
  } else if (nuOut > (Double_t(fNrStrips))){
    Q *= ((Double_t)fNrStrips-nuIn)/(nuOut-nuIn);
    nuOut = Double_t(fNrStrips);
  }
  
  // only one strip hit ?
  if (Int_t(nuIn) == Int_t(nuOut))
    {
      // this strip collected the entire charge
      SmearedQ = SmearCharge(Q);
      if (SmearedQ >= fThreshold)strips.push_back(PndHypStrip(Int_t(nuOut),SmearedQ));
      
      if (fVerboseLevel > 1) std::cout<<" -> 1 strip hit."<<std::endl;
      //return strips;
    } 
  else {
    
    Int_t nrHits = 0;
    Double_t dQ=Q/std::fabs(nuOut-nuIn);
    Double_t dir = (nuOut>nuIn) ? 1. : -1.;
    // calculate portion of track in first strip
    Int_t nextIn = Int_t(nuIn + 0.5+0.5*dir);
    Double_t Q1 = dQ*std::fabs(nextIn-nuIn);
    
    if (fVerboseLevel > 2)
      {
        std::cout<<" part of first strip : "<<nextIn-nuIn<<std::endl ;
        std::cout<<" charge : "<<Q1<<std::endl ;
        std::cout<<" next strip : "<<nextIn<<std::endl ;
      }              
    SmearedQ = SmearCharge(Q1);
    if (SmearedQ >= fThreshold)strips.push_back(PndHypStrip(Int_t(nuIn),SmearedQ));
    nrHits++;
    Q -= Q1;
    
    
    
    // calculate portion of track in last strip
    Int_t prevOut = Int_t(nuOut + 0.5-0.5*dir);
    Double_t Q2 = dQ*std::fabs(nuOut-prevOut);
    
    if (fVerboseLevel > 2){
      std::cout<<" part of last strip : "<<(nuOut-prevOut)<<std::endl ;
      std::cout<<" charge : "<<Q2<<std::endl ;
      std::cout<<" end of previous strip : "<<prevOut<<std::endl ;
    }  
    SmearedQ = SmearCharge(Q2);
    
    if (SmearedQ >= fThreshold)strips.push_back(PndHypStrip(Int_t(nuOut),SmearedQ));
    nrHits++;
    Q -= Q2;
    
    
    
    // Distribute the charge among the intermediate strips
    nextIn = Int_t(nextIn - 0.5 + 0.5*dir);
    prevOut = Int_t(prevOut - 0.5 + 0.5*dir);
    
    
    if (fVerboseLevel > 2) {
      std::cout<<" dir="<<Int_t((dir))<<std::endl;
      std::cout<<" begin="<<nextIn<<" end="<<prevOut<<std::endl;
    }
    
    for (Int_t n = nextIn ; n != prevOut; n += Int_t(dir) )
      {
        
        if (fVerboseLevel > 2) std::cout<<" n = "<<n<<std::endl;
        SmearedQ = SmearCharge(dQ);
        if (SmearedQ >= fThreshold)strips.push_back(PndHypStrip(n,SmearedQ));
        nrHits++;
        Q -= dQ;        //std::cout<<" loop over MORE strips : charge "<<Q<<" dQ "<<dQ<<std::endl;
      }
    if (fVerboseLevel > 2)if (fabs(Q)>1.) std::cout<<" charge Q = "<<Q<<" not detected!"<<std::endl;
    if (fVerboseLevel > 1) std::cout<<" -> "<<nrHits<<" strips hit."<<std::endl;
  }
  
  
  
  return strips;
}

void PndHypCalcStrip::print

(

)

const

Definition at line 268 of file PndHypCalcStrip.cxx.

References fAnchor, fNrFeChannels, fNrStrips, fOrient, fPitch, fStripDir, and Pi.

{
  std::cout<<"-I- PndHypCalcStrip Info :"<<std::endl;
  std::cout<<"     Pitch        = "<<fPitch*10000.<<" um"<<std::endl;
  std::cout<<"     Orientation  = "<<fOrient/TMath::Pi()*180.<<" deg"<<std::endl;
  std::cout<<"     Nr Strips    = "<<fNrStrips<<std::endl;
  std::cout<<"     nr of channels per FE  = "<<fNrFeChannels<<std::endl;
  std::cout<<"     nr of frontends        = "<<fNrStrips/fNrFeChannels<<std::endl;
  std::cout<<"     Anchor Point = "<<fAnchor.X()<<","<<fAnchor.Y()<<")"<<std::endl;
  std::cout<<"     Strip Vector = ("<<fStripDir.X()<<","<<fStripDir.Y()<<")"<<std::endl;
}

void PndHypCalcStrip::SetAnchor ( TVector2  edge )

inline

Definition at line 77 of file PndHypCalcStrip.h.

References fAnchor.

Referenced by PndHypStripHitProducer::Exec().

{fAnchor = edge;}

void PndHypCalcStrip::SetreStrip ( Int_t  reStp )

inline

Definition at line 76 of file PndHypCalcStrip.h.

References fNrStrips.

Referenced by PndHypStripHitProducer::Exec().

{fNrStrips = reStp;}

void PndHypCalcStrip::SetVerboseLevel ( Int_t  level )

inline

Definition at line 74 of file PndHypCalcStrip.h.

References fVerboseLevel.

Referenced by PndHypStripHitProducer::Exec(), and PndHypStripHitProducer::Init().

{ fVerboseLevel = level;}

Double_t PndHypCalcStrip::SmearCharge ( Double_t  charge )

private

133 * Apply gaussian distributed noise to given charge 134 *

Parameters

charge

Equivalent Noise Charge 135 *

Returns

Smeared Charge 136

Definition at line 222 of file PndHypCalcStrip.cxx.

References Double_t, fNoise, fRNG, and fVerboseLevel.

Referenced by GetStrips().

{
  Double_t smeared = fRNG->Gaus(charge,fNoise);
  if (fVerboseLevel > 2) std::cout<<" charge = "<<charge<<", smeared = "<<smeared<<std::endl;
    return smeared;//fRNG->Gaus(charge,fNoise);//smeared;
}

Member Data Documentation

int PndHypCalcStrip::_id

private

Definition at line 100 of file PndHypCalcStrip.h.

TVector2 PndHypCalcStrip::fAnchor

private

ENC.

Definition at line 92 of file PndHypCalcStrip.h.

Referenced by CalcStripFromPoint(), CalcStripPointOnLine(), PndHypCalcStrip(), print(), and SetAnchor().

Double_t PndHypCalcStrip::fNoise

private

charge threshold

Definition at line 91 of file PndHypCalcStrip.h.

Referenced by PndHypCalcStrip(), and SmearCharge().

Int_t PndHypCalcStrip::fNrFeChannels

private

Nr. of strips on active area.

Definition at line 89 of file PndHypCalcStrip.h.

Referenced by CalcChannelfromStrip(), CalcFeChToStrip(), CalcFEfromStrip(), PndHypCalcStrip(), and print().

Int_t PndHypCalcStrip::fNrStrips

private

strip orientation angle

Definition at line 88 of file PndHypCalcStrip.h.

Referenced by CalcFeChToStrip(), GetStrips(), PndHypCalcStrip(), print(), and SetreStrip().

Double_t PndHypCalcStrip::fOrient

private

strip pitch (cm)

Definition at line 87 of file PndHypCalcStrip.h.

Referenced by PndHypCalcStrip(), and print().

TVector2 PndHypCalcStrip::fOrthoDir

private

vector perpendicular to strip direction

Definition at line 94 of file PndHypCalcStrip.h.

Referenced by CalcStripFromPoint(), CalcStripPointOnLine(), and PndHypCalcStrip().

Double_t PndHypCalcStrip::fPitch

private

Definition at line 86 of file PndHypCalcStrip.h.

Referenced by CalcStripFromPoint(), CalcStripPointOnLine(), PndHypCalcStrip(), and print().

TRandom3* PndHypCalcStrip::fRNG

private

vector orthogonal to strip direction

Definition at line 95 of file PndHypCalcStrip.h.

Referenced by PndHypCalcStrip(), and SmearCharge().

Int_t PndHypCalcStrip::fsize

private

Definition at line 97 of file PndHypCalcStrip.h.

TVector2 PndHypCalcStrip::fStripDir

private

anchor point on first strip

Definition at line 93 of file PndHypCalcStrip.h.

Referenced by GetStripDirection(), PndHypCalcStrip(), and print().

Double_t PndHypCalcStrip::fThreshold

private

Nr of Channels per FE.

Definition at line 90 of file PndHypCalcStrip.h.

Referenced by GetStrips(), and PndHypCalcStrip().

Int_t PndHypCalcStrip::fVerboseLevel

private

Random Number Generator.

Definition at line 96 of file PndHypCalcStrip.h.

Referenced by GetStrips(), PndHypCalcStrip(), SetVerboseLevel(), and SmearCharge().

---

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

• PndHypCalcStrip.h
• PndHypCalcStrip.cxx