PndSdsCalcStripDif Class Reference

Class for calculating strip indices from wafer hits. More...

#include <PndSdsCalcStripDif.h>

Public Member Functions
	PndSdsCalcStripDif ()
	PndSdsCalcStripDif (Double_t pitch, Double_t orient, Int_t nrStrips, Int_t nrFeChannels, const TVector2 &firstStripAnchor, Double_t threshold, Double_t noise, Double_t sigma)
	PndSdsCalcStripDif (const PndSdsStripDigiPar *digipar, SensorSide side=kTOP)
std::vector< PndSdsStrip >	GetStrips (Double_t inx, Double_t iny, Double_t inz, Double_t outx, Double_t outy, Double_t outz, Double_t eLoss)
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)
std::ostream &	operator<< (std::ostream &out)
Double_t	ChargeDiffusion (Double_t u, Double_t a, Double_t path, Double_t dir, Double_t Q)
void	Print () const

Private Member Functions
Double_t	CalcStripFromPoint (Double_t x, Double_t y)
Double_t	SmearCharge (Double_t charge)
const 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 to x axis More...
Int_t	fNrFeChannels
	Nr. of strips on active area. More...
TVector2	fAnchor
	Nr of Channels per FE. More...
Double_t	fThreshold
	anchor point on first strip More...
Double_t	fNoise
	charge threshold More...
TVector2	fStripDir
	ENC. More...
TVector2	fOrthoDir
	vector perpendicular to strip direction More...
TRandom3 *	fRNG
	vector orthogonal to strip direction More...
Int_t	fVerboseLevel
	Random Number Generator. More...
Double_t	fSigma

Detailed Description

Class for calculating strip indices from wafer hits.

Author

HG Zaunick hg.za.nosp@m.unic.nosp@m.k@phy.nosp@m.sik..nosp@m.tu-dr.nosp@m.esde.nosp@m.n.de

Definition at line 32 of file PndSdsCalcStripDif.h.

Constructor & Destructor Documentation

PndSdsCalcStripDif::PndSdsCalcStripDif

(

)

Default constructor No reason to use it

Definition at line 17 of file PndSdsCalcStripDif.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 = 0;
  fRNG = new TRandom3();
}

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

Constructor Create Object with all necessary parameters

Parameters

pitch	Pitch between Strips (=width of strips)
orient	Orientation Angle of Strips in x-y-Plane from x-Axis
nrStrips	Number of Strips in the active area
nrFeChannels	Number of Channels per Frontend
firstStripAnchor	Coordinate in x-y-Plane which belongs to the first Strip (Index 0). Usually this is the readout point of the strip
threshold	Charge-Threshold applied to the signal (electrons)
noise	Noise superimposed to the signal (electrons)

Definition at line 28 of file PndSdsCalcStripDif.cxx.

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

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

PndSdsCalcStripDif::PndSdsCalcStripDif	(	const PndSdsStripDigiPar *	digipar,
		SensorSide	side = kTOP
	)

Definition at line 45 of file PndSdsCalcStripDif.cxx.

References cos(), fAnchor, fNoise, fNrFeChannels, fNrStrips, fOrient, fOrthoDir, fPitch, fRNG, fStripDir, fThreshold, fVerboseLevel, PndSdsStripDigiPar::GetBotAnchor(), PndSdsStripDigiPar::GetBotPitch(), PndSdsStripDigiPar::GetNoise(), PndSdsStripDigiPar::GetNrBotFE(), PndSdsStripDigiPar::GetNrFECh(), PndSdsStripDigiPar::GetNrTopFE(), PndSdsStripDigiPar::GetOrient(), PndSdsStripDigiPar::GetSkew(), PndSdsStripDigiPar::GetThreshold(), PndSdsStripDigiPar::GetTopAnchor(), PndSdsStripDigiPar::GetTopPitch(), kBOTTOM, kTOP, Print(), 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 = 0;
  fRNG = new TRandom3();
  if (fVerboseLevel > 0) Print();
}

Member Function Documentation

Int_t PndSdsCalcStripDif::CalcChannelfromStrip

(

Int_t

stripNr

)

const

Calculate Channel index (on Frontend) from strip-index

Parameters

stripNr

strip index

Returns

channel on frontend

Definition at line 438 of file PndSdsCalcStripDif.cxx.

References fNrFeChannels.

Referenced by PndSdsStripHitProducerDif::Exec().

{ return (stripNr%fNrFeChannels); }

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

Calculate global strip index from Frontend number and -channel

Parameters

fe	frontend number
channel	frontend channel
strip	global strip index
side

Definition at line 440 of file PndSdsCalcStripDif.cxx.

References fNrFeChannels, fNrStrips, kBOTTOM, and kTOP.

{
  //Caution! The top side s always the reference side!
  Int_t nr = fe * fNrFeChannels + channel;
  if (nr < fNrStrips) {
    side = kTOP;
  } else {
    nr -= fNrStrips;
    side = kBOTTOM;
  }
  strip = nr;
}

Int_t PndSdsCalcStripDif::CalcFEfromStrip

(

Int_t

stripNr

)

const

Calculate Frontend number from strip-index

Parameters

stripNr

strip index

Returns

Frontend number

Definition at line 436 of file PndSdsCalcStripDif.cxx.

References fNrFeChannels.

Referenced by PndSdsStripHitProducerDif::Exec().

{ return (stripNr/fNrFeChannels); }

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

private

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

Parameters

x	x-coordinate
y	y-coordinate

Returns

strip-parameter nu

Definition at line 425 of file PndSdsCalcStripDif.cxx.

References fAnchor, fOrthoDir, and fPitch.

Referenced by GetStrips().

{
  // fOrthoDir is already set to magnitude == 1., makes it cheaper here.
  return ((x-fAnchor.X())*fOrthoDir.X() + (y-fAnchor.Y())*fOrthoDir.Y())/fPitch;
}

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

Definition at line 431 of file PndSdsCalcStripDif.cxx.

References fAnchor, fOrthoDir, and fPitch.

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

Double_t PndSdsCalcStripDif::ChargeDiffusion	(	Double_t	u,
		Double_t	a,
		Double_t	path,
		Double_t	dir,
		Double_t	Q
	)

Definition at line 457 of file PndSdsCalcStripDif.cxx.

References Double_t, fSigma, n, and sqrt().

Referenced by GetStrips().

{  Int_t N = 10;
  //Int_t n = 0;
  Double_t Npath = path/N;
  Double_t NQ = Q/N;
  Double_t Delta;         //the area of Gaussian from (-infinite,u)
  Double_t dQ=0;   //the Q in the given range
  //define sigma of charge diffusion as a function of sensor thickness
  // suppose the function is linear the sigma woulde be decreased with the thickness
  // if the sigma = 8um when thickness is 300um, then the sigma = 4um for thickness = 150um
  // the sigma shall be changed in the cycle calculation.
  //  Double_t sigma = fSigma
  // cout<<endl;
  //  cout<<"the parameters for charge diffusion are: u, a, path, dir, Q "<< u<<","<<a<<","<<path<<","<<dir<<","<<Q<<endl;
  
  if(path==0.)
  {
    // Delta = 0.5 *(1+TMath::Erf((u-a)/fSigma*std::sqrt(2.)));
    Delta = 0.5 *(1+TMath::Erf((u-a)/fSigma*std::sqrt(2.)));
    dQ = NQ*Delta;
    // cout<<"when path ==0 ,the dQ = "<<dQ<<endl;
  }
  else {
    if(dir>0.)
    {
      for(Double_t da =a,n = 0; n<N; n++)
      { //da = a;
        da +=Npath;
        //  Delta = 0.5*(1+TMath::Erf((u-da)/(fSigma*std::sqrt(2.))));
        //here the sigma square is linear as thickness
        Delta = 0.5 *(1+TMath::Erf((u-a)/(std::sqrt(((N-n)/N))*(fSigma/std::sqrt(2.))*std::sqrt(2.))));//fSigma is from para list and for 300um sensor,so
        dQ +=NQ * Delta;
        //   cout<<"when dir = 1 the delta = "<<Delta<<endl;
        //   cout<<"when dir = 1 the dQ = "<<dQ<<endl;
      }
    } // cout<<"when dir = 1 the dQ = "<<dQ<<endl;
    
    if(dir<0)
    {
      for(Double_t da=a-path,n=0;n<N; n++)
      {
        da+=Npath;
        //  Delta = 0.5*(1+TMath::Erf((u-da)/(fSigma*std::sqrt(2.))));
        Delta = 0.5*(1+TMath::Erf((u-da)/(std::sqrt(((N-n)/N))*(fSigma/std::sqrt(2.))*std::sqrt(2.))));
        dQ +=NQ * Delta;
        //    cout<<"when dir = -1 the delta = "<<Delta<<endl;
        //    cout<<"when dir = -1 the dQ = "<<dQ<<endl;
      }
      //  cout<<"when dir = -1 the dQ = "<<dQ<<endl;
    }//dir<<0
    
  }//path>0
  return dQ;
}

const Double_t PndSdsCalcStripDif::ChargeFromEloss ( Double_t  eloss ) const

inlineprivate

Definition at line 139 of file PndSdsCalcStripDif.h.

Referenced by GetStrips().

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

const TVector2 PndSdsCalcStripDif::GetStripDirection ( ) const

inline

Get Strip Direction (strip orientation angle)

Definition at line 102 of file PndSdsCalcStripDif.h.

References fStripDir.

{return fStripDir;}

std::vector< PndSdsStrip > PndSdsCalcStripDif::GetStrips	(	Double_t	inx,
		Double_t	iny,
		Double_t	inz,
		Double_t	outx,
		Double_t	outy,
		Double_t	outz,
		Double_t	eLoss
	)

Get List of hit channels from given Wafer hit.

Parameters

inx	x-coordinate of entry point
iny	y-coordinate of entry point
inz	z-coordinate of entry point
outx	x-coordinate of exit point
outy	y-coordinate of exit point
outz	z-coordinate of exit point
eLoss	Energy Loss inside the detector element (GeV)

Returns

List of fired strips

the path portion on the first and last strip

Definition at line 74 of file PndSdsCalcStripDif.cxx.

References CalcStripFromPoint(), ChargeDiffusion(), ChargeFromEloss(), Double_t, fabs(), fNrStrips, fPitch, fThreshold, fVerboseLevel, PndSdsStrip::GetIndex(), i, n, out, genfit::tools::QR(), and SmearCharge().

Referenced by PndSdsStripHitProducerDif::Exec().

{
  if (fVerboseLevel > 2) std::cout<<"-I- PndSdsCalcStrip::GetStrips "<<std::endl;
  
  // 2d-Projection of trajectory
  TVector2 in(inx,iny);
  TVector2 out(outx,outy);
  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<PndSdsStrip> strips;
  //Double_t smearedQ;
  
  if (path.Mod()<1E-18) {
    std::cout<<"-W- PndSdsCalcStrip::GetStrips : No Trajectory inside Sensor! (out-in).Mod() = "<<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 Q = ChargeFromEloss(eLoss);//*1E9/3.61; // 3.6 eV/Electron 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<<"-W- PndSdsCalcStrip::GetStrips: 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);
  }
  
  //path direction
  Double_t dir = (nuIn<nuOut)? 1 : -1;
  
  //Strip number for FE channel
  Int_t inStrip = Int_t(nuIn);
  Int_t outStrip = Int_t(nuOut);
  
  Double_t Path1 = (nuIn-inStrip) * fPitch;
  Double_t Path2 = (nuOut-outStrip) * fPitch;
  Double_t pathO = std::fabs(nuOut-nuIn) * fPitch;
  
  // only one strip hit ?
  if (inStrip == outStrip){
    // this strip collected the entire charge  
    // smearedQ = SmearCharge(Q);
    //if (smearedQ >= fThreshold)
    //  strips.push_back(PndSdsStrip(Int_t(nuOut),smearedQ));
    
    //--------------------------Tsito Digi with diffusion------------------------------------------------
    //Double_t path_real = path.Mod();
    Int_t index = inStrip;
    
    Double_t charge1, charge2, dQ, dQR, dQL, QL, QR, QM;
    Int_t id;
    Double_t q_head,q_tail;
    //charge diffusion according to the gaussian smearing
    //between the 3 strips of index-1, index, and index+1
    
    QL = ChargeDiffusion(0,Path1,pathO,dir,Q);//the diffusion charge from -infinite to index
    QR = Q - ChargeDiffusion(fPitch,Path1,pathO,dir,Q);//the diffusion charge from index to infinite
    QM = Q - QR -QL;// the charge left on the strip index
    
    if (fVerboseLevel > 1)
      std::cout<<"the QL,QR,QM are: "<<QL<<","<<QR<<","<<QM<<std::endl;
    
    //noise smearing for the QL,QR and QM
    dQR = SmearCharge(QR);
    dQL = SmearCharge(QL);
    
    //is it only 1 strip involved
    if(dQL<fThreshold && dQR<fThreshold){
      dQ = SmearCharge(QM);
      if(dQ>fThreshold)
        strips.push_back(PndSdsStrip(index,dQ));
      
    }else{
      Int_t i=0;
      Int_t j=1;
      
      if(QL>=fThreshold){
        //charge distribution at the head of the track
        charge1 = ChargeDiffusion((i-1)*fPitch,Path1,pathO,dir,Q);//charge distribution from -infinite to strip of index-1
        if (fVerboseLevel > 1)
          std::cout<<"if QL>fThrehold, the charge on the strip(i-1)*fPitch is :"<<charge1<<",path1= "<<Path1<<std::endl;
        charge2 = ChargeDiffusion(i*fPitch,Path1,pathO,dir,Q);
        if (fVerboseLevel > 1)    
          std::cout<<"if QL>fThrehold, the charge on the strip(i)*fPitch is :"<<charge2<<std::endl;
        q_head = charge2 - charge1;
        id = index + i-1;
        dQ = SmearCharge(q_head);
        if(dQ>fThreshold)
          strips.push_back(PndSdsStrip(id,dQ));
      }
      
      if(QR>=fThreshold){
        //charge distribution at the head of the track
        charge1 = ChargeDiffusion(j*fPitch,Path1,pathO,dir,Q);//charge distribution from -infinite to strip of index-1
        if (fVerboseLevel > 1)      
          std::cout<<"if dQR>fThrehold, the charge on the strip(j)*fPitch is :"<<charge1<<std::endl;
        charge2 = ChargeDiffusion((j+1)*fPitch,Path1,pathO,dir,Q);
        if (fVerboseLevel > 1)      
          std::cout<<"if dQR>fThrehold, the charge on the strip(j+1)*fPitch is :"<<charge2<<std::endl;
        q_tail = charge2 - charge1;
        id = index + j;
        dQ = SmearCharge(q_tail);
        if(dQ>fThreshold)
          strips.push_back(PndSdsStrip(id,dQ));
      }
      
      //the charge left on the center strip
      Double_t restQ = ChargeDiffusion(fPitch,Path1,pathO,dir,Q) - ChargeDiffusion(0,Path1,pathO,dir,Q);
      restQ = SmearCharge(restQ);
      
      if(restQ>fThreshold)//added on Sep.17,2009
        strips.push_back(PndSdsStrip(index,restQ));
      if (fVerboseLevel > 1)
        std::cout<<"the rest Q on the middle strip is :"<<restQ<<std::endl;
    }
    
    //-----------------------End Tsito Digi with diffusion-------------------------------------------------
    
    if (fVerboseLevel > 1) std::cout<<" -> "<<strips.size()<<" strips hit."<<std::endl;
  } else {
    Double_t dQ=Q/std::fabs(nuOut-nuIn);
    //Double_t dPath=path.Mod()/std::fabs(nuOut-nuIn);
    //Double_t dir = (nuOut>nuIn) ? 1. : -1.;
    Int_t nrHits = 0;
    // 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 (Q1 >= fThreshold) strips.push_back(PndSdsStrip(Int_t(nuIn),Q1));
    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 (Q2 >= fThreshold) strips.push_back(PndSdsStrip(Int_t(nuOut),Q2));
    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 (dQ >= fThreshold) strips.push_back(PndSdsStrip(n,dQ));
      nrHits++;
      Q -= dQ;
    }
    
    
    
    //--------------------------Tsito Digi with diffusion-------------------------------------------------
    
    if(strips.size()<1)
      return strips;
    
    //for the strips
    Double_t head_Path;  //, tail_Path;
    std::vector<PndSdsStrip> strip_head;
    std::vector<PndSdsStrip> strip_tail;
    std::vector<PndSdsStrip> strip_middle;
    PndSdsStrip str_head;
    PndSdsStrip str_tail;
    std::vector<PndSdsStrip> str_middle;
    Int_t size = strips.size();
    Double_t total_Q = ChargeFromEloss(eLoss);
    
    if(inStrip<outStrip){
        str_head = strips[0];
        //if(!str_head)cout<<"str_head is available"<<endl;
        str_tail = strips[size-1];
        head_Path = Path1;
        //tail_Path = Path2+(size-1)*fPitch;
      //     cout<<"when inStrip!=outStrip, the head Path is "<<head_Path<<endl;
        for(Int_t i = 1; i<size-1;i++)
          str_middle.push_back(strips[i]);
      
      dir=-1;
      
    }
    else{
        str_tail = strips[0];
        //if(!str_tail)cout<<"str_tail is available, inStrip<outStrip"<<endl;
        str_head = strips[size-1];
        head_Path = Path2;
        //tail_Path = Path1+(size-1)*fPitch;
        //      cout<<"when inStrip>outStrip, the head Path is "<<head_Path<<endl;
      //                cout<<"the dir and size is "<<dir<<","<<size<<endl;
        for(Int_t i = size-2; i>0;i--)
          str_middle.push_back(strips[i]);
      
        dir=1;
      
    }
    
    strips.erase(strips.begin(),strips.begin()+size);
    
    //strips initialization
    Int_t index_head = str_head.GetIndex();
    //Double_t path_head = str_head.GetPath();
    
    Int_t index_tail = str_tail.GetIndex();
    //Double_t path_tail = str_tail.GetPath();
    
    //for the left strip of the head strip
    Double_t charge1, charge2;
    Int_t id;
    //Double_t ipath;
    //Double_t q_head,q_tail;
    Double_t dq_head = ChargeDiffusion(0,head_Path,path.Mod(),dir,total_Q);//the diffusion charge from -infinite to index
    //  cout<<"the dir is "<<dir<<endl;
    Double_t dq_tail = total_Q - ChargeDiffusion(size*fPitch,head_Path,path.Mod(),dir,total_Q);//the diffusion charge from index+size to infinite
    //    cout<<"the dq between -infinite to 0.01 is"<<ChargeDiffusion(size*fPitch,head_Path,path,dir,total_Q)<<endl;
    Double_t q_head = ChargeDiffusion(fPitch,head_Path,path.Mod(),dir,total_Q) - dq_head;// the charge left on the strip index
    //    cout<<"the dq between -infinite to 0.005 is"<<ChargeDiffusion(fPitch,head_Path,path,dir,total_Q)<<endl;
    Double_t q_tail = ChargeDiffusion(size*fPitch,head_Path,path.Mod(),dir,total_Q)
    -ChargeDiffusion((size-1)*fPitch,head_Path,path.Mod(),dir,total_Q);
    //  cout<<"the dir is "<<dir<<endl;
    if (fVerboseLevel > 2)
      std::cout<<"the dq_head,dq_tail,q_head,q_tail are :"<<dq_head<<","<<dq_tail<<","<<q_head<<","<<q_tail<<std::endl;
    
    Double_t dQ_head = SmearCharge(dq_head);
    Double_t dQ_tail = SmearCharge(dq_tail);
    Double_t Q_head = SmearCharge(q_head);
    Double_t Q_tail = SmearCharge(q_tail);
    
    if(dQ_head>fThreshold) {
      id = index_head-1;
        strip_head.push_back(PndSdsStrip(id,dQ_head));
    }
    if(Q_head>fThreshold) {
      id = index_head;
        strip_head.push_back(PndSdsStrip(id,Q_head));
    }
    strips.insert(strips.begin(),strip_head.begin(),strip_head.end());
    
    if(dQ_tail>fThreshold) {
      id = index_tail+1;
        strip_tail.push_back(PndSdsStrip(id,dQ_tail));
    }
    if(Q_tail>fThreshold) {
        id = index_tail;
        strip_head.push_back(PndSdsStrip(id,Q_tail));
    }
    
    strips.insert(strips.end(),strip_tail.begin(),strip_tail.end());
    
    //for the strips between the head and the tail
    //condition for middle strips
    nextIn = Int_t(nuIn+0.5+0.5*dir);
    prevOut = Int_t(nuOut+0.5-0.5*dir);
    
    nextIn=Int_t(nextIn-0.5+0.5*dir);
    prevOut=Int_t(prevOut-0.5+0.5*dir);
    if(nextIn != prevOut){
      for(Int_t i= 1;i< size-1;i++){
          charge1 = ChargeDiffusion(i*fPitch,head_Path,path.Mod(),dir,total_Q);
          charge2 = ChargeDiffusion((i+1)*fPitch,head_Path,path.Mod(),dir,total_Q);
          dQ = charge2 - charge1;
          //std::vector<PndSdsStrip>::iterator iter_ii=strip_middle.begin();
        if(dir==-1)
          id = i;
        else
          id = size-1-i;
          //id = iter_ii->GetIndex();
          //ipath = iter_ii->GetPath();
        Double_t tmpCharge = SmearCharge(dQ);
          if(tmpCharge>fThreshold) strip_middle.push_back(PndSdsStrip(id,tmpCharge));
          //iter_ii++;
        }//loop for all strips between head and tail
        strips.insert(strips.begin()+strip_head.size(),strip_middle.begin(),strip_middle.begin()+strip_middle.size());
      
      /*for(Int_t i=0; i<strip_head.size(); i++)
       strips.push_back(strip_head[i]);
       for(Int_t i=0; i<strip_middle.size(); i++)
       strips.push_back(strip_middle.[i]);
       for(Int_t i=0; i<strip_tail.size(); i++)
       strips.push_back(strip_tail[i]);*/
    }//if strips fired exceed 3
    
    //-----------------------End Tsito Digi with diffusion-------------------------------------------------
    
    if (fVerboseLevel > 2) if (fabs(Q)>1.) std::cout<<" charge Q = "<<Q<<" not detected!"<<std::endl;
    if (fVerboseLevel > 1) std::cout<<" -> "<<strips.size()<<" strips hit."<<std::endl;
  }
  
  return strips;
}

std::ostream& PndSdsCalcStripDif::operator<<

(

std::ostream &

out

)

void PndSdsCalcStripDif::Print

(

)

const

Definition at line 512 of file PndSdsCalcStripDif.cxx.

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

Referenced by PndSdsCalcStripDif().

{
  std::cout<<"-I- PndSdsCalcStrip Info :"<<std::endl;
  std::cout<<"     pitch                  = "<<fPitch*10000.<<" um"<<std::endl;
  std::cout<<"     orientation angle      = "<<fOrient/TMath::Pi()*180.<<" deg"<<std::endl;
  std::cout<<"     nr of 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-direction vector = ("<<fStripDir.X()<<","<<fStripDir.Y()<<")"<<std::endl;
}

void PndSdsCalcStripDif::SetVerboseLevel ( Int_t  level )

inline

Definition at line 104 of file PndSdsCalcStripDif.h.

References fVerboseLevel.

Referenced by PndSdsStripHitProducerDif::Init().

{ fVerboseLevel = level;}

Double_t PndSdsCalcStripDif::SmearCharge ( Double_t  charge )

private

Apply gaussian distributed noise to given charge

Parameters

charge

Equivalent Noise Charge

Returns

Smeared Charge

Definition at line 418 of file PndSdsCalcStripDif.cxx.

References Double_t, fNoise, fRNG, and fVerboseLevel.

Referenced by GetStrips().

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

Member Data Documentation

TVector2 PndSdsCalcStripDif::fAnchor

private

Nr of Channels per FE.

Definition at line 116 of file PndSdsCalcStripDif.h.

Referenced by CalcStripFromPoint(), CalcStripPointOnLine(), PndSdsCalcStripDif(), and Print().

Double_t PndSdsCalcStripDif::fNoise

private

charge threshold

Definition at line 118 of file PndSdsCalcStripDif.h.

Referenced by PndSdsCalcStripDif(), and SmearCharge().

Int_t PndSdsCalcStripDif::fNrFeChannels

private

Nr. of strips on active area.

Definition at line 115 of file PndSdsCalcStripDif.h.

Referenced by CalcChannelfromStrip(), CalcFeChToStrip(), CalcFEfromStrip(), PndSdsCalcStripDif(), and Print().

Int_t PndSdsCalcStripDif::fNrStrips

private

strip orientation angle to x axis

Definition at line 114 of file PndSdsCalcStripDif.h.

Referenced by CalcFeChToStrip(), GetStrips(), PndSdsCalcStripDif(), and Print().

Double_t PndSdsCalcStripDif::fOrient

private

strip pitch (cm)

Definition at line 113 of file PndSdsCalcStripDif.h.

Referenced by PndSdsCalcStripDif(), and Print().

TVector2 PndSdsCalcStripDif::fOrthoDir

private

vector perpendicular to strip direction

Definition at line 120 of file PndSdsCalcStripDif.h.

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

Double_t PndSdsCalcStripDif::fPitch

private

Definition at line 112 of file PndSdsCalcStripDif.h.

Referenced by CalcStripFromPoint(), CalcStripPointOnLine(), GetStrips(), PndSdsCalcStripDif(), and Print().

TRandom3* PndSdsCalcStripDif::fRNG

private

vector orthogonal to strip direction

Definition at line 121 of file PndSdsCalcStripDif.h.

Referenced by PndSdsCalcStripDif(), and SmearCharge().

Double_t PndSdsCalcStripDif::fSigma

private

Definition at line 123 of file PndSdsCalcStripDif.h.

Referenced by ChargeDiffusion(), and PndSdsCalcStripDif().

TVector2 PndSdsCalcStripDif::fStripDir

private

ENC.

Definition at line 119 of file PndSdsCalcStripDif.h.

Referenced by GetStripDirection(), PndSdsCalcStripDif(), and Print().

Double_t PndSdsCalcStripDif::fThreshold

private

anchor point on first strip

Definition at line 117 of file PndSdsCalcStripDif.h.

Referenced by GetStrips(), and PndSdsCalcStripDif().

Int_t PndSdsCalcStripDif::fVerboseLevel

private

Random Number Generator.

Definition at line 122 of file PndSdsCalcStripDif.h.

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

---

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

• PndSdsCalcStripDif.h
• PndSdsCalcStripDif.cxx