PndMdtParamDigi Class Reference

#include <PndMdtParamDigi.h>

Inheritance diagram for PndMdtParamDigi:

Classes
struct	AvaBinType
struct	ClusInfo

Public Member Functions
	PndMdtParamDigi ()
	~PndMdtParamDigi ()
void	SetVerbose (Int_t v)
Bool_t	Init ()
PndMdtParamDigi &	SetParams (Int_t ptlType, TVector3 iniP, TVector3 iniPos, TVector3 finalPos, Double_t stripLen=100.)
void	UseNoise (Bool_t swith)
void	UseDetailedSim (Bool_t swith=kTRUE)
void	UsePlot (Bool_t swith=kTRUE)
void	UseGaussianAmp (Bool_t swith)
void	SetOptimization (Int_t val)
void	SetNoiseWidth (Double_t anode, Double_t strip)
std::vector< std::pair< Int_t, Double_t > >	GetFiredInfo ()
void	Compute (Bool_t useConvolution=kTRUE)
const std::vector< Double_t > &	GetWireSignal () const
const std::map< Int_t, std::vector< Double_t > > &	GetStripSignals () const
void	Draw ()

Private Member Functions
void	GetSignal (Bool_t useConvolution=kTRUE)
Bool_t	Digitize (Double_t &time, Double_t &amp)
Bool_t	Digitize (Int_t stripNo, Double_t &time, Double_t &amp)
void	GetRawSignalbySimAvalanche (Double_t fNoiseLevel=1.)
void	GetRawSignalbyWeightingAvalanche (Double_t fNoiseLevel=1.)
void	AddNoise (Double_t fNoiseLevel=1., Int_t isAnode=1)
void	ApplyTransferFunction (Double_t *fSignalData, Int_t nSize=fSamplingSize)
Double_t	GetElectronDriftTime (TVector2 iniPos)
Double_t	GetShiftTime (TVector2 fWireSurfacePos, TVector2 fProductionPos)
void	GetMPVofPrimaryIonization (Int_t particleType, const TVector3 &momentum, ValueErrorType &val) const
Int_t	GetAmplicationFactor (Int_t particleType, Double_t momentum) const
void	SamplingPosition (TVector2 fDirection, TVector2 &fIonProductionPos)
Double_t	GetMeanFreePath (Int_t ptlType, Double_t mom) const
	ClassDef (PndMdtParamDigi, 1)

Private Attributes
Int_t	fSamplingRate
Int_t	fSamplingInterval
std::vector< Double_t >	fSignalDataAnode
std::map< Int_t, std::vector < Double_t > >	fSignalDataStripM
TF1 *	fTrF
Bool_t	fUseNoise
Bool_t	fGaussianAmp
Double_t	fNoiseSigmaAnode
Double_t	fNoiseSigmaStrip
Int_t	fNumofTruncation
Int_t	fParticleType
TVector3	fParticleMomentum
TVector3	fInitPostion
TVector3	fExitPosition
Bool_t	fParamsChanged
Double_t	fStripLength
Double_t	fRestMass [5]
Bool_t	fDetailedSim
TGraphErrors *	gFreePath [5]
Bool_t	fUseAvaHist
TH1F *	hAvaSize
Double_t	cWIRERADIUS
Double_t	cCELLSIZE
Double_t	cSTRIPWIDTH
Double_t	cMAXRADIUS
Double_t	cMAPFACTOR
Double_t	cUNITLOGDR
Double_t	cUNITPHI
TH1F *	hAnodeI1d [NRAD]
TH1F *	hCathodeI2d [NRAD][NPHI]
Int_t	fNr [NRAD]
Int_t	fNrNphi [NRAD *NPHI]
Bool_t	fUsePlot
TH1F *	hIndex
TH2F *	h2dIndex
TH2F *	hTrack
TH1F *	hIonNum
TH2F *	hAvaPos
TH1F *	hAvaNum
TH1F *	hAva1D
TH2F *	hAva2D
TH1F *	hExp1
TH1F *	hExp2
TH1F *	hGen
TH1F *	hLan
Int_t	fVerbose
std::vector< AvaBinType >	fProbFunc1D
std::vector< AvaBinType >	fProbFunc2D
Double_t	SPEEDOFLIGHT
CLHEP::RandGeneral *	fRandAva

Static Private Attributes
static const Int_t	fSamplingSize = 20
static const Int_t	NRAD = 100
static const Int_t	NPHI = 100

Detailed Description

Definition at line 27 of file PndMdtParamDigi.h.

Constructor & Destructor Documentation

PndMdtParamDigi::PndMdtParamDigi

(

)

Default constructor

Definition at line 30 of file PndMdtParamDigi.cxx.

References cCELLSIZE, cMAPFACTOR, cMAXRADIUS, cSTRIPWIDTH, cUNITLOGDR, cUNITPHI, cWIRERADIUS, fRestMass, fSamplingInterval, fSamplingRate, CAMath::Log(), NPHI, NRAD, Pi, and SPEEDOFLIGHT.

  : TNamed("PndMdtParamDigi", "parameterized simulation of MDT")
  , fSamplingRate(1e8)
  , fTrF(0)
  , fUseNoise(kTRUE)
  , fGaussianAmp(kTRUE)
  , fNoiseSigmaAnode(1e-3)
  , fNoiseSigmaStrip(1e-5)
  , fNumofTruncation(1.)
  , fParamsChanged(kTRUE)
  , fDetailedSim(kFALSE)
  , fUseAvaHist(kFALSE)
  , fUsePlot(kTRUE)
  , fVerbose(0)
{
  cWIRERADIUS = 0.0025001;//unit: cm
  cCELLSIZE   = 0.5;//unit: cm
  cSTRIPWIDTH = 1.0;//unit: cm
  cMAXRADIUS = 0.00525;//not changeable
  cMAPFACTOR = 1000.;//not changeable

  cUNITLOGDR = TMath::Log(1. + cMAPFACTOR*(cMAXRADIUS - cWIRERADIUS))/NRAD;
  cUNITPHI = TMath::Pi()*2./NPHI;
  fSamplingInterval = 1e9/fSamplingRate;
  fRestMass[0] = 0.00510998928;
  fRestMass[1] = 0.1056583715;
  fRestMass[2] = 0.13957018;
  fRestMass[3] = 0.493667;
  fRestMass[4] = 0.938272046;// GeV

  SPEEDOFLIGHT = 299792458;//m/s
}

PndMdtParamDigi::~PndMdtParamDigi

(

)

Destructor

Definition at line 158 of file PndMdtParamDigi.cxx.

{ }

Member Function Documentation

void PndMdtParamDigi::AddNoise ( Double_t  fNoiseLevel = 1., Int_t  isAnode = 1  )

private

Definition at line 621 of file PndMdtParamDigi.cxx.

References Double_t, fNoiseSigmaAnode, fNoiseSigmaStrip, fSamplingSize, fSignalDataAnode, fSignalDataStripM, and i.

Referenced by GetRawSignalbySimAvalanche(), and GetRawSignalbyWeightingAvalanche().

{
  Double_t fSigmaAnode = fNoiseLevel*fNoiseSigmaAnode;//micro ampere
  Double_t fSigmaStrip = fNoiseLevel*fNoiseSigmaStrip;//micro ampere

  std::map<Int_t, std::vector<Double_t> >::iterator it = fSignalDataStripM.begin();
  std::map<Int_t, std::vector<Double_t> >::iterator end = fSignalDataStripM.end();
  while( it != end){
    std::vector<Double_t>& fSignalDataStrip = it->second;
    for(Int_t i=0; i< fSamplingSize; ++i){
      fSignalDataStrip[i] += gRandom->Gaus(0, fSigmaStrip);
    }
    ++ it;
  }
  for(Int_t i=0; i< fSamplingSize; ++i){
    fSignalDataAnode[i] += gRandom->Gaus(0, fSigmaAnode);
  }
}

void PndMdtParamDigi::ApplyTransferFunction ( Double_t *  fSignalData, Int_t  nSize = fSamplingSize  )

private

Definition at line 590 of file PndMdtParamDigi.cxx.

References Double_t, fSamplingSize, fTrF, i, and t.

Referenced by GetSignal().

{
  Double_t dt=0.01;//microsecond
  Double_t t(0.);
  Double_t fIntegral(0.);
  Double_t dtau(0.);
  Double_t fF(0.);
  Double_t fG(0.);
  Double_t tprime(0.);
  Int_t iBin(1);
  dtau = (fTrF->GetXmax() - fTrF->GetXmin())/1.e3;
  Double_t fSigCopy[fSamplingSize];
  if(fTrF){
    memcpy(fSigCopy, fSignalData, fSamplingSize*sizeof(Double_t));
    //std::vector<Double_t> fSigCopy(fSignalData.begin(),fSignalData.end());
    for(Int_t i=0; i< fSamplingSize; ++i){
      t = i*dt;//microsecond
      fIntegral = 0.;
      tprime = fTrF->GetXmin();
      while( tprime  < fTrF->GetXmax()){
        iBin = (Int_t)(tprime/dt);
        fF = (iBin <0 || iBin > fSamplingSize) ? 0 : fSigCopy[iBin];
        fG = (t-tprime < fTrF->GetXmin()) || (t-tprime > fTrF->GetXmax()) ? 0 :fTrF->Eval((t-tprime));//convent to microsecond
        fIntegral += fF*fG*dtau;
        tprime += dtau;
      }
      fSignalData[i] = fIntegral;
    }
  }
}

PndMdtParamDigi::ClassDef ( PndMdtParamDigi  , 1    )

private

void PndMdtParamDigi::Compute

(

Bool_t

useConvolution = kTRUE

)

Definition at line 242 of file PndMdtParamDigi.cxx.

References GetSignal().

Referenced by PndMdtPointsToWaveform::exec_e(), and PndMdtPointsToWaveform::exec_t().

{ 
  GetSignal(useConvolution); 
}

Bool_t PndMdtParamDigi::Digitize ( Double_t &  time, Double_t &  amp  )

private

Definition at line 640 of file PndMdtParamDigi.cxx.

References CAMath::Abs(), Bool_t, Double_t, fNoiseSigmaAnode, fSamplingInterval, fSamplingSize, fSignalDataAnode, and is.

Referenced by GetFiredInfo().

{
  time = -1;
  amp = -9999;
  Double_t fPeak = 0.;
  Bool_t NotFound = kTRUE;
  for(size_t is=0; is < fSamplingSize; ++ is){
    if(TMath::Abs(fSignalDataAnode[is]) > TMath::Abs(fPeak)) fPeak = fSignalDataAnode[is];
    if(TMath::Abs(fSignalDataAnode[is]) > 5*fNoiseSigmaAnode && NotFound) {
      time = is*fSamplingInterval;//nano seconds
      NotFound = kFALSE;
    }
  }
  amp = fPeak;
  return time > 0;
}

Bool_t PndMdtParamDigi::Digitize ( Int_t  stripNo, Double_t &  time, Double_t &  amp  )

private

Definition at line 656 of file PndMdtParamDigi.cxx.

References CAMath::Abs(), Bool_t, Double_t, fNoiseSigmaStrip, fSamplingInterval, fSamplingSize, fSignalDataStripM, and is.

{
  time = -1;
  amp = -9999;
  Double_t fPeak = 0.;
  Bool_t NotFound = kTRUE;
  std::vector<Double_t>& fSignalDataStrip = fSignalDataStripM[id];
  for(size_t is=0; is < fSamplingSize; ++ is){
    if(TMath::Abs(fSignalDataStrip[is]) > TMath::Abs(fPeak)) fPeak = fSignalDataStrip[is];
    if(TMath::Abs(fSignalDataStrip[is]) > 5*fNoiseSigmaStrip && NotFound) {
      time = is*fSamplingInterval;//nano seconds
      NotFound = kFALSE;
    }
  }
  amp = fPeak;
  return time > 0;
}

void PndMdtParamDigi::Draw

(

)

Definition at line 182 of file PndMdtParamDigi.cxx.

References c, c1, c2, fSamplingInterval, fSamplingSize, fSignalDataAnode, fSignalDataStripM, fUsePlot, h2dIndex, hAva1D, hAva2D, hAvaPos, hExp1, hGen, hIndex, hIonNum, hLan, hTrack, i, name, and TString.

{
  if(fUsePlot)
  {
    TCanvas* c = new TCanvas("cPndMdtParamDigi", "PndMdtParamDigi", 800, 600);
    c->Divide(2,3);
    c->cd(1);
    hTrack->Draw("colz");
    c->cd(2);
    //hAva2D->Draw("colz");
    hIonNum->Draw();
    c->cd(3);
    hAvaPos->Draw("colz");
    c->cd(5);
    hIndex->DrawNormalized();
    c->cd(5);
    hAva1D->SetLineColor(kRed);
    hAva1D->SetLineWidth(2);
    hAva1D->Draw("same");
    c->cd(6);
    hAva2D->Draw("colz");
    c->cd(4);
    h2dIndex->Draw("colz");
    TCanvas* c2 = new TCanvas("dist", "PndMdtParamDigi", 800, 600);
    c2->Divide(2,2);
    c2->cd(1);
    hExp1->Draw();
    c2->cd(2);
    c2->cd(3);
    hGen->Draw();
    c2->cd(4);
    hLan->Draw();
  }
  TCanvas* c1 = new TCanvas("cSignal", "Signals", 900, 900);
  c1->Divide(3, (fSignalDataStripM.size() + 1)/3);
  c1->cd(1);
  TH1F* hAI = new TH1F("hAI","induced current @anode", fSamplingSize, 0, fSamplingSize*fSamplingInterval/1.e3);
  hAI->GetYaxis()->SetTitle("I(#muA)");
  for(Int_t i=0; i < fSamplingSize; ++i){
    hAI->SetBinContent(i+1, fSignalDataAnode[i]);
  }
  hAI->Draw();
  std::map<Int_t, std::vector<Double_t> >::iterator it = fSignalDataStripM.begin();
  std::map<Int_t, std::vector<Double_t> >::iterator end = fSignalDataStripM.end();
  Int_t ipad=1;
  for(; it != end; ++it, ++ipad){
    c1->cd(1+ipad);
    TString name="hCI";
    name += it->first;
    TString title="induced current @strip ";
    title += it->first;
    TH1F* hCI = new TH1F(name, title, fSamplingSize, 0, fSamplingSize*fSamplingInterval/1.e3);
    hCI->GetYaxis()->SetTitle("I(#muA)");
    std::vector<Double_t>& fSignalDataStrip = it->second;
    for(Int_t i=0; i < fSamplingSize; ++i){
      hCI->SetBinContent(i+1, fSignalDataStrip[i]);
    }
    hCI->Draw();
  }
}

Int_t PndMdtParamDigi::GetAmplicationFactor ( Int_t  particleType, Double_t  momentum  ) const

private

Definition at line 684 of file PndMdtParamDigi.cxx.

References Double_t.

{
  TGraphErrors* gAmp(0);
  if(!gAmp){
    gAmp= new TGraphErrors(9);
    gAmp->SetPoint(0, 0.05, 5440.5);
    gAmp->SetPointError(0, 0, 2402.1);
    gAmp->SetPoint(1, 0.10, 5683.1);
    gAmp->SetPointError(1, 0, 2661.3);
    gAmp->SetPoint(2, 0.20, 5774.8);
    gAmp->SetPointError(2, 0, 2803.5);
    gAmp->SetPoint(3, 0.30, 5629.1);
    gAmp->SetPointError(3, 0, 2590.4);
    gAmp->SetPoint(4, 0.50, 5708.1);
    gAmp->SetPointError(4, 0, 2734.2);
    gAmp->SetPoint(5, 1.00, 5667.6);
    gAmp->SetPointError(5, 0, 2755.5);
    gAmp->SetPoint(6, 2.00, 5958.2);
    gAmp->SetPointError(7, 0, 2712.1);
    gAmp->SetPoint(7, 5.00, 5643.8);
    gAmp->SetPointError(7, 0, 2389.5);
    gAmp->SetPoint(8, 10.00, 5749.3);
    gAmp->SetPointError(8, 0, 2552.6);
  }
  const Double_t fScaleFactor = 1.0;//simulation calibration to experiments 
  Double_t fMeanAmp = gAmp->Eval(momentum);
  Double_t fAmpErr =  gAmp->GetErrorY(momentum);
  return  (Int_t) fScaleFactor*gRandom->Gaus(fMeanAmp, fAmpErr);
}

Double_t PndMdtParamDigi::GetElectronDriftTime ( TVector2  iniPos )

inlineprivate

Definition at line 102 of file PndMdtParamDigi.h.

References cWIRERADIUS, Double_t, log(), p1, p2, and x.

Referenced by GetRawSignalbySimAvalanche(), and GetRawSignalbyWeightingAvalanche().

    {
      const Double_t  p0 =  5.35513e+03;
      const Double_t  p1 =  3.92749e+02;
      const Double_t  p2 = -5.41433e+03;
      const Double_t  p3 = -3.31810e+03;
      Double_t x = iniPos.Mod() - cWIRERADIUS;
      Double_t x1 = log(1.+x);
      return (x*(p0 + p1*pow(x, 0.25)) + x1*(p2 + p3*x1));
    }

std::vector< std::pair< Int_t, Double_t > > PndMdtParamDigi::GetFiredInfo

(

)

Definition at line 545 of file PndMdtParamDigi.cxx.

References Digitize(), Double_t, fSignalDataStripM, and GetSignal().

{
  std::vector<std::pair<Int_t, Double_t> > fFiredInfo;
  Double_t mTime;
  Double_t mAmplitude;
  GetSignal();

  if(Digitize(mTime, mAmplitude))
    fFiredInfo.push_back(std::pair<Int_t, Double_t>(-1, mTime));

  std::map<Int_t, std::vector<Double_t> >::iterator it = fSignalDataStripM.begin();
  std::map<Int_t, std::vector<Double_t> >::iterator end = fSignalDataStripM.end();
  for(; it != end; ++it){
    if(Digitize(it->first, mTime, mAmplitude)){
      fFiredInfo.push_back(std::pair<Int_t, Double_t>(it->first, mTime));
    }
  }
  return fFiredInfo;
}

Double_t PndMdtParamDigi::GetMeanFreePath ( Int_t  ptlType, Double_t  mom  ) const

inlineprivate

Definition at line 143 of file PndMdtParamDigi.h.

References Double_t, fRestMass, and gFreePath.

    {
      Double_t mMass = fRestMass[ptlType];
      //Double_t mEnergy = sqrt(mom*mom + mMass*mMass);
      Double_t mBeta = mom/mMass;
      Double_t mPath = gFreePath[ptlType]->Eval(mBeta);
      return gRandom->Exp(mPath);
    }

void PndMdtParamDigi::GetMPVofPrimaryIonization ( Int_t  particleType, const TVector3 &  momentum, ValueErrorType &  val  ) const

private

Definition at line 720 of file PndMdtParamDigi.cxx.

References Double_t, exp(), fRestMass, log(), p1, and p2.

Referenced by GetRawSignalbySimAvalanche(), and GetRawSignalbyWeightingAvalanche().

{
  Double_t p0 = 9.83401e-01;
  Double_t p1 = 7.42822e-02;
  Double_t p2 = 5.62250e-01;
  Double_t p3 = 2.22743e-01;
  //index: electron 0, muon 1, pion 2, kaon, 3, proton 4;
  Double_t fMass = fRestMass[particleType];
  Double_t betagamma = fMomentum.Mag()/fMass;
  val.first = p0*exp(p1*log(1+betagamma)) + p2/(pow(betagamma, p3));
  val.second = betagamma < 0.3 ? 0.234 : 0.334;
}

void PndMdtParamDigi::GetRawSignalbySimAvalanche ( Double_t  fNoiseLevel = 1. )

private

Definition at line 247 of file PndMdtParamDigi.cxx.

References AddNoise(), cMAPFACTOR, cSTRIPWIDTH, cUNITLOGDR, cUNITPHI, cWIRERADIUS, Double_t, fExitPosition, fInitPostion, fNumofTruncation, fParamsChanged, fParticleMomentum, fParticleType, fRestMass, fSamplingInterval, fSamplingSize, fSignalDataAnode, fSignalDataStripM, fUseNoise, fUsePlot, GetElectronDriftTime(), GetMPVofPrimaryIonization(), gFreePath, h2dIndex, hAnodeI1d, hAvaNum, hAvaPos, hAvaSize, hCathodeI2d, hIndex, hIonNum, hTrack, jj, CAMath::Log(), NPHI, NRAD, SamplingPosition(), and SPEEDOFLIGHT.

{

  if(!fParamsChanged) return;

  fParamsChanged = kFALSE;
  memset(&fSignalDataAnode[0], 0, fSamplingSize*sizeof(Double_t));
  //memset(&fSignalDataStrip[0], 0, fSamplingSize*sizeof(Double_t));
  fSignalDataStripM.clear();

  TVector3 fUnitMotion = (fExitPosition-fInitPostion).Unit();
  Double_t fMaxPath = (fExitPosition-fInitPostion).Mag();
  Double_t mBeta = fParticleMomentum.Mag()/fRestMass[fParticleType];
  if(fParticleType == 4 && mBeta < 0.2) return;
  Double_t mPath = gFreePath[fParticleType]->Eval(mBeta);
  if(mPath < 0.) return;
  const Double_t fScaleFactor = 40.;

  Int_t jj;
  Int_t ii;

  ValueErrorType NumofPrimaryIonizaionMpv; 
  GetMPVofPrimaryIonization(fParticleType, fParticleMomentum, NumofPrimaryIonizaionMpv);


  TVector3 fCurrentPosofCluster = fInitPostion;
  Int_t    fNumofCluster(0);
  Int_t    fNumofSingleIonization  =0;
  Int_t    fNumofTotalIonizaion = 0;
  Int_t    fNumofIonsofThisCluster = 0;
  Int_t    fNumofTotalIons = 0;
  TVector2 fIonProductionPos;
  Double_t fStepLen = 0.;
  Double_t fTrackLen = 0.;
  Double_t fTime = 0.;

  TStopwatch sw;
  sw.Start();
  do
  {
    fStepLen = gRandom->Exp(mPath);
    fTime += fStepLen/(mBeta*SPEEDOFLIGHT*1.e7);//nano second
    fTrackLen += fStepLen; 
    fCurrentPosofCluster += fStepLen*fUnitMotion;

    if(fCurrentPosofCluster.XYvector().Mod() <= cWIRERADIUS)
      continue;

    fNumofSingleIonization = gRandom->Landau(NumofPrimaryIonizaionMpv.first, NumofPrimaryIonizaionMpv.second);
    fNumofTotalIonizaion  += fNumofSingleIonization;

    if(fUsePlot){
      hIonNum->Fill(fNumofSingleIonization);
      hTrack->Fill(fCurrentPosofCluster.X(), fCurrentPosofCluster.Y());
    }

    Int_t iSign = fCurrentPosofCluster.Z() >= 0 ? 1 : -1;
    Int_t zIndex = Int_t(fCurrentPosofCluster.Z()/cSTRIPWIDTH + iSign*0.5);
    zIndex = zIndex >= 0 ? zIndex*2 : -2*zIndex -1;


    //cout<<"fSignalDataStrip.size()#"<<fSignalDataStrip.size()<<endl;
    std::vector<Double_t>& fSignalDataStrip = fSignalDataStripM[zIndex];
    if(fSignalDataStrip.size() != fSamplingSize)
      fSignalDataStrip.resize(fSamplingSize);

    if(fNumofSingleIonization > 0)
    {
      fNumofIonsofThisCluster = hAvaSize->GetRandom();
      fNumofTotalIons  += fNumofIonsofThisCluster;
      if(fUsePlot)
        hAvaNum->Fill(fNumofIonsofThisCluster);

      if( fNumofIonsofThisCluster > 10 ) 
      {
        //Double_t lx = fCurrentPosofCluster.XYvector().Mod() - cWIRERADIUS;
        //Double_t lx1 = log(1.+lx);
        //clusInfo.fStartTime = fTime +  (lx*(p0 + p1*pow(lx, 0.25)) + lx1*(p2 + p3*lx1));

        Double_t fStartTime = fTime + GetElectronDriftTime(fCurrentPosofCluster.XYvector());//micro second
        Int_t tShift = Int_t(fStartTime/fSamplingInterval);
        if(tShift >= fSamplingSize) tShift = fSamplingSize-1;
        //Int_t refPhiIndex = Int_t(fCurrentPosofCluster.Phi()/cUNITPHI); //[R.K. 01/2017] unused variable?
        //if(refPhiIndex>=100) refPhiIndex = 99;

        TVector2 unitV = fCurrentPosofCluster.XYvector().Unit();
        std::map<Int_t, Int_t> fAnodeIndex;//fired region
        std::map<Int_t, Int_t> fStripIndex;//fired region

        for(ii = 0; ii < fNumofIonsofThisCluster; ++ii)
        {//avalanche profile
          //get position of this ion
          SamplingPosition(unitV, fIonProductionPos);
          //cout<<"sampling done, r = "<<fIonProductionPos.Mod()<<endl;
          Int_t rIndex = Int_t(TMath::Log(1. + cMAPFACTOR*(fIonProductionPos.Mod() - cWIRERADIUS))/cUNITLOGDR);
          Int_t phiIndex = Int_t(fIonProductionPos.Phi()/cUNITPHI);
          if(fUsePlot){
            hIndex->Fill(rIndex);
            //h2dIndex->Fill(rIndex, phiIndex-refPhiIndex < 0 ? 100 + phiIndex-refPhiIndex : phiIndex-refPhiIndex);
            h2dIndex->Fill(rIndex, phiIndex);
            hAvaPos->Fill(fIonProductionPos.X(), fIonProductionPos.Y());
          }
          //cout<<"plot done"<<endl;
          //if(rIndex>=NRAD) rIndex = NRAD-1;
          //if(phiIndex >= NPHI) phiIndex = NPHI -1;
          if(rIndex>=NRAD) continue;
          if(phiIndex >= NPHI) continue;
          //cout<<"rIndex = "<<rIndex<<", phiIndex = "<<phiIndex<<endl;
          ++ fAnodeIndex[rIndex];
          ++ fStripIndex[rIndex + phiIndex*NRAD];
        }
        std::map<Int_t, Int_t>::const_iterator mit = fAnodeIndex.begin();
        std::map<Int_t, Int_t>::const_iterator mend = fAnodeIndex.end();
        for(; mit != mend; ++mit){
          if(mit->second > fNumofTruncation){
            for(ii = tShift, jj=1; ii < fSamplingSize; ++ii, ++jj){
              fSignalDataAnode[ii] += mit->second * hAnodeI1d[mit->first]->GetBinContent(jj);
            }
          }
        }
        mit = fStripIndex.begin();
        mend = fStripIndex.end();
        for(; mit != mend; ++mit){
          Double_t fN = mit->second;
          if(fN > fNumofTruncation)
          {
            Int_t rIdx = mit->first%NRAD;
            Int_t phiIdx = mit->first/NRAD;
            //cout<<"Index ("<<rIdx<<", "<<phiIdx<<") = "<<fN<<endl;
            for(ii = tShift, jj=1; ii < fSamplingSize; ++ii, ++jj){
              fSignalDataStrip[ii] += fScaleFactor* mit->second * hCathodeI2d[rIdx][phiIdx]->GetBinContent(jj);
            }
          }
        }
      }
    }
    ++ fNumofCluster;
  }while(fTrackLen < fMaxPath);

  if(fUseNoise){
    AddNoise(fNoiseLevel);
  }
  sw.Stop();
  cout<<"PndMdtParamDigi::GetRawSignal, cost "<<sw.CpuTime()<<" s"
    <<" clusters#"<<fNumofCluster
    <<" primary ionization#"<<fNumofTotalIonizaion
    <<" total ions#"<<fNumofTotalIons<<endl;
}

void PndMdtParamDigi::GetRawSignalbyWeightingAvalanche ( Double_t  fNoiseLevel = 1. )

private

Definition at line 395 of file PndMdtParamDigi.cxx.

References AddNoise(), cMAPFACTOR, CAMath::Cos(), cSTRIPWIDTH, cUNITLOGDR, cUNITPHI, cWIRERADIUS, Double_t, exp(), fExitPosition, fInitPostion, fNumofTruncation, fParamsChanged, fParticleMomentum, fParticleType, fProbFunc1D, fProbFunc2D, fRestMass, fSamplingInterval, fSamplingSize, fSignalDataAnode, fSignalDataStripM, fUseNoise, fUsePlot, GetElectronDriftTime(), GetMPVofPrimaryIonization(), gFreePath, h2dIndex, hAnodeI1d, hAvaNum, hAvaPos, hAvaSize, hCathodeI2d, hIndex, hIonNum, hTrack, PndMdtParamDigi::AvaBinType::Index, jj, NPHI, NRAD, PndMdtParamDigi::AvaBinType::Probabilty, CAMath::Sin(), SPEEDOFLIGHT, x, and y.

Referenced by GetSignal().

{
  static Double_t fScaleFactor = 40.;

  if(!fParamsChanged) return;

  fParamsChanged = kFALSE;
  memset(&fSignalDataAnode[0], 0, fSamplingSize*sizeof(Double_t));
  //memset(&fSignalDataStrip[0], 0, fSamplingSize*sizeof(Double_t));
  fSignalDataStripM.clear();

  TVector3 fUnitMotion = (fExitPosition-fInitPostion).Unit();
  Double_t fMaxPath = (fExitPosition-fInitPostion).Mag();
  Double_t mBeta = fParticleMomentum.Mag()/fRestMass[fParticleType];
  Double_t mPath = gFreePath[fParticleType]->Eval(mBeta);
  if(mPath < 0.) return;

  Int_t jj;
  Int_t ii;

  ValueErrorType NumofPrimaryIonizaionMpv; 
  GetMPVofPrimaryIonization(fParticleType, fParticleMomentum, NumofPrimaryIonizaionMpv);


  TVector3 fCurrentPosofCluster = fInitPostion;
  Int_t    fNumofCluster(0);
  Int_t    fNumofSingleIonization  =0;
  Int_t    fNumofTotalIonizaion = 0;
  Int_t    fNumofIonsofThisCluster = 0;
  Int_t    fNumofTotalIons = 0;
  Double_t fStepLen = 0.;
  Double_t fTrackLen = 0.;
  Double_t fTime = 0.;

  //TStopwatch sw;
  //sw.Start();
  do
  {
    fStepLen = gRandom->Exp(mPath);
    fTime += fStepLen/(mBeta*SPEEDOFLIGHT*1.e7);//nano second
    fTrackLen += fStepLen; 
    fCurrentPosofCluster += fStepLen*fUnitMotion;

    if(fCurrentPosofCluster.XYvector().Mod() <= cWIRERADIUS)
      continue;

    fNumofSingleIonization = gRandom->Landau(NumofPrimaryIonizaionMpv.first, NumofPrimaryIonizaionMpv.second);
    fNumofTotalIonizaion  += fNumofSingleIonization;

    if(fUsePlot){
      hIonNum->Fill(fNumofSingleIonization);
      hTrack->Fill(fCurrentPosofCluster.X(), fCurrentPosofCluster.Y());
    }

    Int_t iSign = fCurrentPosofCluster.Z() >= 0 ? 1 : -1;
    Int_t zIndex = Int_t(fCurrentPosofCluster.Z()/cSTRIPWIDTH + iSign*0.5);
    zIndex = zIndex >= 0 ? zIndex*2 : -2*zIndex -1;


    //cout<<"fSignalDataStrip.size()#"<<fSignalDataStrip.size()<<endl;
    std::vector<Double_t>& fSignalDataStrip = fSignalDataStripM[zIndex];
    if(fSignalDataStrip.size() != fSamplingSize)
      fSignalDataStrip.resize(fSamplingSize);

    if(fNumofSingleIonization > 0)
    {
      fNumofIonsofThisCluster = hAvaSize->GetRandom();
      fNumofTotalIons  += fNumofIonsofThisCluster;
      if(fUsePlot)
        hAvaNum->Fill(fNumofIonsofThisCluster);

      if( fNumofIonsofThisCluster > 1 ) 
      {
        //Double_t lx = fCurrentPosofCluster.XYvector().Mod() - cWIRERADIUS;
        //Double_t lx1 = log(1.+lx);
        //clusInfo.fStartTime = fTime +  (lx*(p0 + p1*pow(lx, 0.25)) + lx1*(p2 + p3*lx1));

        Double_t fStartTime = fTime + GetElectronDriftTime(fCurrentPosofCluster.XYvector());//micro second
        Int_t tShift = Int_t(fStartTime/fSamplingInterval);
        if(tShift >= fSamplingSize) tShift = fSamplingSize-1;
        Int_t refPhiIndex = Int_t(fCurrentPosofCluster.Phi()/cUNITPHI);
        //wire signal
        //cout<<"1D, number of fired bins #"<<fProbFunc1D.size()<<endl;
        std::vector<AvaBinType>::const_iterator vit = fProbFunc1D.begin();
        std::vector<AvaBinType>::const_iterator end = fProbFunc1D.end();
        for(; vit != end; ++ vit)
        {
          const AvaBinType& bin = *vit;
          Double_t fN = fNumofIonsofThisCluster* bin.Probabilty;
          if(fN > fNumofTruncation){
            for(ii = tShift, jj=1; ii < fSamplingSize; ++ ii, ++jj){
              fSignalDataAnode[ii] += fN*hAnodeI1d[bin.Index]->GetBinContent(jj);
            }
          }
          if(fUsePlot){
            hIndex->SetBinContent(bin.Index+1, hIndex->GetBinContent(bin.Index+1) + fN);
          }
        }
        //cout<<"wire signal done"<<endl;
        //strip signal
        //cout<<"2D, number of fired bins #"<<fProbFunc2D.size()<<endl;
        vit = fProbFunc2D.begin();
        end = fProbFunc2D.end();
        std::map<Int_t, Double_t> fInfo;
        for(; vit != end; ++ vit)
        {
          const AvaBinType& bin = *vit;
          Int_t fN = fNumofIonsofThisCluster* bin.Probabilty;
          if(fN > fNumofTruncation){
            //retrieve the rIndex, phiIndex;
            Int_t rIdx = bin.Index%NRAD;
            Int_t phiIdx = bin.Index/NRAD;
            phiIdx += refPhiIndex;
            phiIdx %= NPHI;
            fInfo[rIdx + phiIdx*NRAD] = fN;
            for(ii = tShift, jj=1; ii < fSamplingSize; ++ ii, ++jj){
              fSignalDataStrip[ii] += fScaleFactor* fN *hCathodeI2d[rIdx][phiIdx]->GetBinContent(jj);
            }
            if(fUsePlot){
              Int_t irr = rIdx;
              Int_t ipp = phiIdx;
              h2dIndex->SetBinContent(irr+1, ipp+1, h2dIndex->GetBinContent(irr+1, ipp+1) + fN);
              //Int_t rIndex = Int_t(TMath::Log(1. + cMAPFACTOR*(fIonProductionPos.Mod() - cWIRERADIUS))/cUNITLOGDR);
              //Int_t phiIndex = Int_t(fIonProductionPos.Phi()/cUNITPHI);
              Double_t PHI = (ipp+1)*cUNITPHI;
              Double_t RADIUS = (exp(irr+1) - 1.)*cUNITLOGDR/cMAPFACTOR + cWIRERADIUS;
              Double_t x = RADIUS*TMath::Cos(PHI);
              Double_t y = RADIUS*TMath::Sin(PHI);
              //Int_t ix = (x + 0.01)/0.0002;
              //Int_t iy = (y + 0.01)/0.0002;
              for(Int_t iv = 0; iv < fN; ++ iv)
                hAvaPos->Fill(x,y);
              //hAvaPos->SetBinContent(ix, iy, hAvaPos->GetBinContent(ix, iy) + fNrNphi[ir][ipp]);
            }
          }
        }
      }
    }
    ++ fNumofCluster;
  }while(fTrackLen < fMaxPath);

  if(fUseNoise){
    AddNoise(fNoiseLevel);
  }
  //sw.Stop();
  //cout<<"PndMdtParamDigi::GetRawSignal, cost "<<sw.CpuTime()<<" s"
  //  <<" clusters#"<<fNumofCluster
  //  <<" primary ionization#"<<fNumofTotalIonizaion
  //  <<" total ions#"<<fNumofTotalIons<<endl;
}

Double_t PndMdtParamDigi::GetShiftTime ( TVector2  fWireSurfacePos, TVector2  fProductionPos  )

inlineprivate

Definition at line 113 of file PndMdtParamDigi.h.

References Double_t, p1, and p2.

    {
      const Double_t p0 = 1.81940e-8;//unit: micro sencond
      const Double_t p1 = 1.80248;
      const Double_t p2 = 3.68652e+2;
      const Double_t p3 = 1.67353e+3;
      //const Double_t mean = 0; //[R.K. 01/2017] unused variable?
      //const Double_t sigma = 2.73477e-3;//micro second //[R.K. 01/2017] unused variable?
      Double_t dr = (fProductionPos-fWireSurfacePos).Mod();
      return (((p3*dr+p2)*dr+p1)*dr+p0 /*+ gRandom->Gaus(mean, sigma)*/)*1.e3;//nano sencod
    }

void PndMdtParamDigi::GetSignal ( Bool_t  useConvolution = kTRUE )

private

Definition at line 570 of file PndMdtParamDigi.cxx.

References ApplyTransferFunction(), fSignalDataAnode, fSignalDataStripM, and GetRawSignalbyWeightingAvalanche().

Referenced by Compute(), and GetFiredInfo().

{

  GetRawSignalbyWeightingAvalanche(1.);
  //GetRawSignalbySimAvalanche(1.);
  //convolution with readout circuit
  if(useConvolution){
    ApplyTransferFunction(&fSignalDataAnode[0]);
    //for_each(fSignalDataAnode.begin(), fSignalDataAnode.end(), print());
    //cout<<endl;

    std::map<Int_t, std::vector<Double_t> >::iterator it = fSignalDataStripM.begin();
    std::map<Int_t, std::vector<Double_t> >::iterator end = fSignalDataStripM.end();
    while( it != end){
      std::vector<Double_t>& fSignalDataStrip = it->second;
      ApplyTransferFunction(&fSignalDataStrip[0]);
      ++ it;
    }
  }
}

const std::map<Int_t, std::vector<Double_t> >& PndMdtParamDigi::GetStripSignals ( ) const

inline

Definition at line 66 of file PndMdtParamDigi.h.

References fSignalDataStripM.

Referenced by PndMdtPointsToWaveform::exec_e(), and PndMdtPointsToWaveform::exec_t().

{ return fSignalDataStripM; }

const std::vector<Double_t>& PndMdtParamDigi::GetWireSignal ( ) const

inline

Definition at line 65 of file PndMdtParamDigi.h.

References fSignalDataAnode.

Referenced by PndMdtPointsToWaveform::exec_e(), and PndMdtPointsToWaveform::exec_t().

{ return fSignalDataAnode; }

Bool_t PndMdtParamDigi::Init

(

)

Definition at line 62 of file PndMdtParamDigi.cxx.

References file, fProbFunc1D, fProbFunc2D, fSamplingSize, fSignalDataAnode, fTrF, fUseAvaHist, fUsePlot, gFreePath, h2dIndex, hAnodeI1d, hAva1D, hAva2D, hAvaNum, hAvaPos, hAvaSize, hCathodeI2d, hExp1, hGen, hIndex, hIonNum, hLan, hTrack, ir, NPHI, NRAD, and TString.

Referenced by PndMdtPointsToWaveform::Init().

{
  fTrF = new TF1("fTrF", "(x*6./0.015)**6*exp(-x*6./0.015)", 0., 2.);
  if(! fTrF) return kFALSE;

  TString dir(gSystem->Getenv("VMCWORKDIR"));
  TString mdtParamsFile(dir + "/macro/params/MdtDigiParams.root");
  TFile* file = TFile::Open(mdtParamsFile);
  if( ! file->IsOpen()) return kFALSE;

  gFreePath[0] = (TGraphErrors*)file->Get("gElec");
  gFreePath[1] = (TGraphErrors*)file->Get("gMuon");
  gFreePath[2] = (TGraphErrors*)file->Get("gPion");
  gFreePath[3] = (TGraphErrors*)file->Get("gKaon");
  gFreePath[4] = (TGraphErrors*)file->Get("gProton");
  hAvaSize = (TH1F*)file->Get("hava");
  //hAnodeI = (TH1F*)file->Get("hAnodeI");
  //hCathodeIx= (TH1F*)file->Get("hCathodeIx");
  //hCathodeIy= (TH1F*)file->Get("hCathodeIy");
  cout<<"PndMdtParamDigi::PndMdtParamDigi(), test"<<endl;
  cout<<"mean frem path @1 GeV of electron "<<gFreePath[0]->Eval(1.)<<endl;
  cout<<"mean frem path @1 GeV of muon "<<gFreePath[1]->Eval(1.)<<endl;
  cout<<"mean frem path @1 GeV of pion "<<gFreePath[2]->Eval(1.)<<endl;
  cout<<"mean frem path @1 GeV of kaon "<<gFreePath[3]->Eval(1.)<<endl;
  cout<<"mean frem path @1 GeV of proton "<<gFreePath[4]->Eval(1.)<<endl;
  //cout<<"avalanche size #"<<hAvaSize->GetRandom()<<endl;
  //cout<<"induced current on anode @t = 0.05us, "<<hAnodeI->GetBinContent(5)<<endl;
  //cout<<"induced current on cathode @t = 0.05us, theta 0, "<<hCathodeIx->GetBinContent(5)<<endl;
  //cout<<"induced current on cathode @t = 0.05us, theta PI/2, "<<hCathodeIy->GetBinContent(5)<<endl;
  //file->Close();
  fUseAvaHist = kTRUE;
  //double* lProbData = new double[hAvaSize->GetNbinsX()];
  //for(Int_t ib=0; ib < hAvaSize->GetNbinsX(); ++ ib)
  //  lProbData[ib] = hAvaSize->GetBinContent(ib+1);
  //fRandAva = new CLHEP::RandGeneral(lProbData, hAvaSize->GetNbinsX());

  TString mdtParamsFileI(dir + "/macro/params/MdtDigiParamsI.root");
  TFile* file2 = TFile::Open(mdtParamsFileI);
  if(! file2->IsOpen()) return kFALSE;

  for(Int_t ir=0; ir < NRAD; ++ ir)
  {
    TString hisName("Ir");
    hisName += ir;
    hAnodeI1d[ir] = (TH1F*)file2->Get(hisName);
  }
  for(Int_t ir=0; ir < NRAD; ++ ir)
  {
    for(Int_t ip=0; ip <NPHI; ++ ip){
      TString hisName("Ir");
      hisName += ir;
      hisName += "t";
      hisName += ip;
      hCathodeI2d[ir][ip] = (TH1F*)file2->Get(hisName);
    }
  }
  hAva1D = (TH1F*)file2->Get("hAva1D");
  hAva2D = (TH2F*)file2->Get("hAva2D");
  cout<<"1D avalanche PDF "<<hAva1D->GetTitle()<<endl;
  cout<<"2D avalanche PDF "<<hAva2D->GetTitle()<<endl;

  for(Int_t ir=0; ir < hAva1D->GetNbinsX(); ++ ir)
    if(hAva1D->GetBinContent(ir+1) > 0){
      fProbFunc1D.push_back(AvaBinType(ir, hAva1D->GetBinContent(ir+1)));
    }

  for(Int_t ir=0; ir < hAva2D->GetNbinsX(); ++ ir)
    for(Int_t ip=0; ip < hAva2D->GetNbinsY(); ++ ip)
      if(hAva2D->GetBinContent(ir+1, ip+1) > 0){
        fProbFunc2D.push_back(AvaBinType(ir+ip*NRAD, hAva2D->GetBinContent(ir+1, ip+1)));
      }


  if(fUsePlot)
  {
    hIonNum = new TH1F("hIonizationNum", "number of primary ionization", 20, 0, 20);
    hTrack  = new TH2F("hTrack", "trajectory of primary track", 100, -0.5, 0.5, 100, -0.5, 0.5);
    hAvaPos = new TH2F("hAvaPos", "avalanche profile", 100, -0.01, 0.01, 100, -0.01, 0.01);
    hAvaNum = new TH1F("hAvaNum", "avalanche", 20000, 0, 20000);
    hIndex  = new TH1F("hIndex", "raidal avlanche structure", 100,0, 100);
    hIndex->GetXaxis()->SetTitle("Start bin of Drift");
    h2dIndex  = new TH2F("h2dIndex", "raidal avlanche structure", 100,0, 100, 100, 0, 100);
    h2dIndex->GetXaxis()->SetTitle("Start bin of Drift");

    hExp1 = new TH1F("expf", "exponential dis", 100,0,1);
    //hExp2;
    hGen = new TH1F("gen", "general dist", 100, 0, 10000);
    hLan = new TH1F("lan", "landau dis", 20, 0, 20);
  }
  fSignalDataAnode.resize(fSamplingSize);

  return kTRUE;
}

void PndMdtParamDigi::SamplingPosition ( TVector2  fDirection, TVector2 &  fIonProductionPos  )

private

Definition at line 732 of file PndMdtParamDigi.cxx.

References cos(), cWIRERADIUS, Double_t, and sin().

Referenced by GetRawSignalbySimAvalanche().

{
  const Double_t cThetaSigma = 2.67859e-01;
  //sampling theta
  Double_t mTheta = gRandom->Gaus(0., cThetaSigma);
  //sampling radial
  const Double_t mTauR1 = 1.51900e-3;
  const Double_t mTauR2 = 5.53950e-4;
  const Double_t mTauR3 = 2.94876e-4;
  Double_t mRadius(0);
  Double_t prob = gRandom->Rndm();
  if(prob < 2.38210768816632112e-02){
    mRadius = gRandom->Exp(mTauR1);
  }else if(prob < 4.34694439053706305e-01){
    mRadius = gRandom->Exp(mTauR2);
  } else{
    mRadius = gRandom->Exp(mTauR3);
  }
  mRadius += cWIRERADIUS;
  Double_t mX = mRadius*cos(mTheta);
  Double_t mY = mRadius*sin(mTheta);
  //unitV.Print();
  TVector2 unitV = fDirection.Unit();

  fIonProductionPos.Set( unitV.X()*mX - unitV.Y()*mY, unitV.Y()*mX + unitV.X()*mY);
  //fIonProductionPos.Print();
}

void PndMdtParamDigi::SetNoiseWidth ( Double_t  anode, Double_t  strip  )

inline

Definition at line 56 of file PndMdtParamDigi.h.

References fNoiseSigmaAnode, fNoiseSigmaStrip, and strip.

Referenced by PndMdtPointsToWaveform::Init().

{ fNoiseSigmaAnode = anode; fNoiseSigmaStrip = strip; }

void PndMdtParamDigi::SetOptimization ( Int_t  val )

inline

Definition at line 55 of file PndMdtParamDigi.h.

References fNumofTruncation, and val.

Referenced by PndMdtPointsToWaveform::Init().

{ fNumofTruncation = val; }

PndMdtParamDigi & PndMdtParamDigi::SetParams	(	Int_t	ptlType,
		TVector3	iniP,
		TVector3	iniPos,
		TVector3	finalPos,
		Double_t	stripLen = 100.
	)

Definition at line 160 of file PndMdtParamDigi.cxx.

References fExitPosition, fInitPostion, fParamsChanged, fParticleMomentum, fParticleType, fSignalDataStripM, fStripLength, fVerbose, and TString.

Referenced by PndMdtPointsToWaveform::exec_e(), and PndMdtPointsToWaveform::exec_t().

{
  fParticleType = ptlType;
  fParticleMomentum =  iniP;
  fInitPostion = iniPos;
  fExitPosition = finalPos;
  fStripLength = stripLen;
  fParamsChanged = kTRUE;
  fSignalDataStripM.clear();

  if(fVerbose > 1){
    cout<<"==============================================="<<endl;
    TString partileName[5] = {"Elec", "Muon", "Pion", "Kaon", "Proton"};
    cout<<"PndMdtParamDigi::SetParams, "<<partileName[ptlType]<<endl
      <<" initial momentum#"<<iniP.X()<<", "<<iniP.Y()<<", "<<iniP.Z()<<", "<<iniP.Mag()<<"GeV"<<endl
      <<" initial position#"<<iniPos.X()<<", "<<iniPos.Y()<<", "<<iniPos.Z()<<", "<<iniPos.Mag()<<"cm"<<endl
      <<" final position#"<<finalPos.X()<<", "<<finalPos.Y()<<", "<<finalPos.Z()<<", "<<finalPos.Mag()<<"cm"<<endl;
    cout<<"==============================================="<<endl;
  }
  return *this;
}

void PndMdtParamDigi::SetVerbose ( Int_t  v )

inline

Definition at line 37 of file PndMdtParamDigi.h.

References fVerbose, and v.

{ fVerbose = v; }

void PndMdtParamDigi::UseDetailedSim ( Bool_t  swith = kTRUE )

inline

Definition at line 52 of file PndMdtParamDigi.h.

References fDetailedSim.

Referenced by PndMdtPointsToWaveform::Init().

{ fDetailedSim = swith; }

void PndMdtParamDigi::UseGaussianAmp ( Bool_t  swith )

inline

Definition at line 54 of file PndMdtParamDigi.h.

References fGaussianAmp.

Referenced by PndMdtPointsToWaveform::Init().

{ fGaussianAmp = swith; }

void PndMdtParamDigi::UseNoise ( Bool_t  swith )

inline

Definition at line 51 of file PndMdtParamDigi.h.

References fUseNoise.

Referenced by PndMdtPointsToWaveform::Init().

{ fUseNoise = swith; }

void PndMdtParamDigi::UsePlot ( Bool_t  swith = kTRUE )

inline

Definition at line 53 of file PndMdtParamDigi.h.

References fUsePlot.

Referenced by PndMdtPointsToWaveform::Init().

{ fUsePlot = swith; }

Member Data Documentation

Double_t PndMdtParamDigi::cCELLSIZE

private

Definition at line 164 of file PndMdtParamDigi.h.

Referenced by PndMdtParamDigi().

Double_t PndMdtParamDigi::cMAPFACTOR

private

Definition at line 167 of file PndMdtParamDigi.h.

Referenced by GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and PndMdtParamDigi().

Double_t PndMdtParamDigi::cMAXRADIUS

private

Definition at line 166 of file PndMdtParamDigi.h.

Referenced by PndMdtParamDigi().

Double_t PndMdtParamDigi::cSTRIPWIDTH

private

Definition at line 165 of file PndMdtParamDigi.h.

Referenced by GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and PndMdtParamDigi().

Double_t PndMdtParamDigi::cUNITLOGDR

private

Definition at line 168 of file PndMdtParamDigi.h.

Referenced by GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and PndMdtParamDigi().

Double_t PndMdtParamDigi::cUNITPHI

private

Definition at line 169 of file PndMdtParamDigi.h.

Referenced by GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and PndMdtParamDigi().

Double_t PndMdtParamDigi::cWIRERADIUS

private

Definition at line 163 of file PndMdtParamDigi.h.

Referenced by GetElectronDriftTime(), GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), PndMdtParamDigi(), and SamplingPosition().

Bool_t PndMdtParamDigi::fDetailedSim

private

Definition at line 138 of file PndMdtParamDigi.h.

Referenced by UseDetailedSim().

TVector3 PndMdtParamDigi::fExitPosition

private

Definition at line 129 of file PndMdtParamDigi.h.

Referenced by GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and SetParams().

Bool_t PndMdtParamDigi::fGaussianAmp

private

Definition at line 94 of file PndMdtParamDigi.h.

Referenced by UseGaussianAmp().

TVector3 PndMdtParamDigi::fInitPostion

private

Definition at line 128 of file PndMdtParamDigi.h.

Referenced by GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and SetParams().

Double_t PndMdtParamDigi::fNoiseSigmaAnode

private

Definition at line 95 of file PndMdtParamDigi.h.

Referenced by AddNoise(), Digitize(), and SetNoiseWidth().

Double_t PndMdtParamDigi::fNoiseSigmaStrip

private

Definition at line 96 of file PndMdtParamDigi.h.

Referenced by AddNoise(), Digitize(), and SetNoiseWidth().

Int_t PndMdtParamDigi::fNr[NRAD]

private

Definition at line 172 of file PndMdtParamDigi.h.

Int_t PndMdtParamDigi::fNrNphi[NRAD *NPHI]

private

Definition at line 173 of file PndMdtParamDigi.h.

Int_t PndMdtParamDigi::fNumofTruncation

private

Definition at line 97 of file PndMdtParamDigi.h.

Referenced by GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and SetOptimization().

Bool_t PndMdtParamDigi::fParamsChanged

private

Definition at line 130 of file PndMdtParamDigi.h.

Referenced by GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and SetParams().

TVector3 PndMdtParamDigi::fParticleMomentum

private

Definition at line 127 of file PndMdtParamDigi.h.

Referenced by GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and SetParams().

Int_t PndMdtParamDigi::fParticleType

private

Definition at line 126 of file PndMdtParamDigi.h.

Referenced by GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and SetParams().

std::vector<AvaBinType> PndMdtParamDigi::fProbFunc1D

private

Definition at line 210 of file PndMdtParamDigi.h.

Referenced by GetRawSignalbyWeightingAvalanche(), and Init().

std::vector<AvaBinType> PndMdtParamDigi::fProbFunc2D

private

Definition at line 211 of file PndMdtParamDigi.h.

Referenced by GetRawSignalbyWeightingAvalanche(), and Init().

CLHEP::RandGeneral* PndMdtParamDigi::fRandAva

private

Definition at line 215 of file PndMdtParamDigi.h.

Double_t PndMdtParamDigi::fRestMass[5]

private

Definition at line 132 of file PndMdtParamDigi.h.

Referenced by GetMeanFreePath(), GetMPVofPrimaryIonization(), GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and PndMdtParamDigi().

Int_t PndMdtParamDigi::fSamplingInterval

private

Definition at line 89 of file PndMdtParamDigi.h.

Referenced by Digitize(), Draw(), GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and PndMdtParamDigi().

Int_t PndMdtParamDigi::fSamplingRate

private

Definition at line 88 of file PndMdtParamDigi.h.

Referenced by PndMdtParamDigi().

const Int_t PndMdtParamDigi::fSamplingSize = 20

staticprivate

Definition at line 87 of file PndMdtParamDigi.h.

Referenced by AddNoise(), ApplyTransferFunction(), Digitize(), Draw(), GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and Init().

std::vector<Double_t> PndMdtParamDigi::fSignalDataAnode

private

Definition at line 90 of file PndMdtParamDigi.h.

Referenced by AddNoise(), Digitize(), Draw(), GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), GetSignal(), GetWireSignal(), and Init().

std::map<Int_t, std::vector<Double_t> > PndMdtParamDigi::fSignalDataStripM

private

Definition at line 91 of file PndMdtParamDigi.h.

Referenced by AddNoise(), Digitize(), Draw(), GetFiredInfo(), GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), GetSignal(), GetStripSignals(), and SetParams().

Double_t PndMdtParamDigi::fStripLength

private

Definition at line 131 of file PndMdtParamDigi.h.

Referenced by SetParams().

TF1* PndMdtParamDigi::fTrF

private

Definition at line 92 of file PndMdtParamDigi.h.

Referenced by ApplyTransferFunction(), and Init().

Bool_t PndMdtParamDigi::fUseAvaHist

private

Definition at line 153 of file PndMdtParamDigi.h.

Referenced by Init().

Bool_t PndMdtParamDigi::fUseNoise

private

Definition at line 93 of file PndMdtParamDigi.h.

Referenced by GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and UseNoise().

Bool_t PndMdtParamDigi::fUsePlot

private

Definition at line 187 of file PndMdtParamDigi.h.

Referenced by Draw(), GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), Init(), and UsePlot().

Int_t PndMdtParamDigi::fVerbose

private

Definition at line 201 of file PndMdtParamDigi.h.

Referenced by SetParams(), and SetVerbose().

TGraphErrors* PndMdtParamDigi::gFreePath[5]

private

Definition at line 141 of file PndMdtParamDigi.h.

Referenced by GetMeanFreePath(), GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and Init().

TH2F* PndMdtParamDigi::h2dIndex

private

Definition at line 189 of file PndMdtParamDigi.h.

Referenced by Draw(), GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and Init().

TH1F* PndMdtParamDigi::hAnodeI1d[NRAD]

private

Definition at line 170 of file PndMdtParamDigi.h.

Referenced by GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and Init().

TH1F* PndMdtParamDigi::hAva1D

private

Definition at line 194 of file PndMdtParamDigi.h.

Referenced by Draw(), and Init().

TH2F* PndMdtParamDigi::hAva2D

private

Definition at line 195 of file PndMdtParamDigi.h.

Referenced by Draw(), and Init().

TH1F* PndMdtParamDigi::hAvaNum

private

Definition at line 193 of file PndMdtParamDigi.h.

Referenced by GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and Init().

TH2F* PndMdtParamDigi::hAvaPos

private

Definition at line 192 of file PndMdtParamDigi.h.

Referenced by Draw(), GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and Init().

TH1F* PndMdtParamDigi::hAvaSize

private

Definition at line 154 of file PndMdtParamDigi.h.

Referenced by GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and Init().

TH1F* PndMdtParamDigi::hCathodeI2d[NRAD][NPHI]

private

Definition at line 171 of file PndMdtParamDigi.h.

Referenced by GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and Init().

TH1F* PndMdtParamDigi::hExp1

private

Definition at line 197 of file PndMdtParamDigi.h.

Referenced by Draw(), and Init().

TH1F* PndMdtParamDigi::hExp2

private

Definition at line 198 of file PndMdtParamDigi.h.

TH1F* PndMdtParamDigi::hGen

private

Definition at line 199 of file PndMdtParamDigi.h.

Referenced by Draw(), and Init().

TH1F* PndMdtParamDigi::hIndex

private

Definition at line 188 of file PndMdtParamDigi.h.

Referenced by Draw(), GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and Init().

TH1F* PndMdtParamDigi::hIonNum

private

Definition at line 191 of file PndMdtParamDigi.h.

Referenced by Draw(), GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and Init().

TH1F* PndMdtParamDigi::hLan

private

Definition at line 200 of file PndMdtParamDigi.h.

Referenced by Draw(), and Init().

TH2F* PndMdtParamDigi::hTrack

private

Definition at line 190 of file PndMdtParamDigi.h.

Referenced by Draw(), GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and Init().

const Int_t PndMdtParamDigi::NPHI = 100

staticprivate

Definition at line 161 of file PndMdtParamDigi.h.

Referenced by GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), Init(), and PndMdtParamDigi().

const Int_t PndMdtParamDigi::NRAD = 100

staticprivate

Definition at line 160 of file PndMdtParamDigi.h.

Referenced by GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), Init(), and PndMdtParamDigi().

Double_t PndMdtParamDigi::SPEEDOFLIGHT

private

Definition at line 214 of file PndMdtParamDigi.h.

Referenced by GetRawSignalbySimAvalanche(), GetRawSignalbyWeightingAvalanche(), and PndMdtParamDigi().

---

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

• PndMdtParamDigi.h
• PndMdtParamDigi.cxx