PndFsmSttPid Class Reference

#include <PndFsmSttPid.h>

Inheritance diagram for PndFsmSttPid:

Public Member Functions
	PndFsmSttPid ()
	PndFsmSttPid (ArgList &par)
virtual	~PndFsmSttPid ()
virtual PndFsmResponse *	respond (PndFsmTrack *t)
const std::string &	detName ()
void	setName (std::string &name)
Bool_t	doesPid () const
void	setStorePid (Bool_t doespid=kTRUE)

Protected Member Functions
virtual bool	setParameter (std::string &name, std::string &value)
void	parseParameterList (ArgList &par)

Protected Attributes
std::string	_detName
double	_efficiency
TRandom3 *	_rand
TDatabasePDG *	_fdbPDG
Bool_t	_doesPid

Private Member Functions
bool	detected (PndFsmTrack *t) const
double	compdEdx (double p, double M)
double	gauss (double x, double x0, double s)
void	initParameters ()
void	print (std::ostream &o)
bool	setParameter (std::string &name, double value)

Private Attributes
double	_thtMin
double	_thtMax
double	_ptmin
double	_dEdxRes

Detailed Description

Definition at line 41 of file PndFsmSttPid.h.

Constructor & Destructor Documentation

PndFsmSttPid::PndFsmSttPid

(

)

Default ctor.

Definition at line 57 of file PndFsmSttPid.cxx.

References _thtMax, _thtMin, and initParameters().

{
  initParameters();

  _thtMin=_thtMin*M_PI/180.0;
  _thtMax=_thtMax*M_PI/180.0;
  //print(std::cout);
}

PndFsmSttPid::PndFsmSttPid

(

ArgList &

par

)

Definition at line 66 of file PndFsmSttPid.cxx.

References _thtMax, _thtMin, initParameters(), and PndFsmAbsDet::parseParameterList().

{
  initParameters();
  //set default parameter values and parses a parameter list
  //i.e. std::list<std::string> of the form
  //"a=1" "b=2" "c=3"
  parseParameterList(par);

  _thtMin=_thtMin*M_PI/180.0;
  _thtMax=_thtMax*M_PI/180.0;
  //print(std::cout);
}

PndFsmSttPid::~PndFsmSttPid ( )

virtual

Destructor

Definition at line 83 of file PndFsmSttPid.cxx.

{
}

Member Function Documentation

double PndFsmSttPid::compdEdx ( double  p, double  M  )

private

Definition at line 149 of file PndFsmSttPid.cxx.

References a, C(), dEdX, log(), me, sqrt(), X, Z, and z.

Referenced by respond().

{
  double dEdX;

  p*=1000;
  M*=1000;

  const double Z=10;
  const double A=20;
  const double z=1;//charge of incident particle in unit of e

  double beta;
  beta=p/sqrt(M*M+p*p);//CalculateBeta(KE,M);

  double gamma;
  gamma=1./sqrt(1-beta*beta);

  const double I=10e-6*Z;//0.000188;//MeV
  const double me=0.511;//Mev/c2

  double Wmax;
  Wmax=(2*me*beta*beta*gamma*gamma) / (1 + 2*gamma*me/M + (me/M)*(me/M));

  //const double C1=0.1535;//MeV cm2/g

  double X,X0,X1;
  double kappa=0.307075;
  X0=0.201;
  X1=3;
  X=log10(beta*gamma);
  double delta;
  double C,a;
  C=-5.217;
  a=0.196;

  if(X<=X0)
    delta=0;
  else if(X<=X1)
    delta=2*log(10.)*X+C+a*(X1-X)*(X1-X)*(X1-X);
  else
    delta=2*log(10.)*X+C;

  dEdX= ( kappa * (Z/A)   * z*z /(beta*beta)) * (log(2*me*beta*beta*gamma*gamma*Wmax / (I*I)) - 2*beta*beta - delta); //-Dshell??

  return dEdX;
}

bool PndFsmSttPid::detected ( PndFsmTrack *  t ) const

private

Definition at line 197 of file PndFsmSttPid.cxx.

References PndFsmAbsDet::_efficiency, _ptmin, PndFsmAbsDet::_rand, _thtMin, PndFsmTrack::charge(), PndFsmTrack::p4(), and theta.

Referenced by respond().

{
    double theta = t->p4().Theta();
    double p_t=t->p4().Vect().Perp(TVector3(0.,0.,1.));
    double charge=t->charge();

    return ( charge!=0.0 && theta>=_thtMin && theta<=_thtMax && p_t>_ptmin  && _rand->Rndm()<=_efficiency);
}

const std::string& PndFsmAbsDet::detName ( )

inlineinherited

Definition at line 74 of file PndFsmAbsDet.h.

References PndFsmAbsDet::_detName.

Referenced by PndFastSim::AddDetector(), PndFsmIdealPid::print(), PndFsmSimpleVtx::print(), PndFsmMvdPid::print(), PndFsmSimpleTracker::print(), PndFsmDetTemplate::print(), PndFsmMdtPid::print(), PndFsmTof::print(), PndFsmAbsDet::print(), print(), PndFsmMdcFS::print(), PndFsmRich::print(), PndFsmMdcTS::print(), PndFsmMvd2::print(), PndFsmDrcDisc::print(), PndFsmDrcBarrel::print(), PndFsmStt::print(), PndFsmMvd::print(), PndFsmEffTracker::print(), PndFsmResponse::print(), PndFastSim::Register(), PndFastSim::smearTrack(), and PndFastSim::sumResponse().

{return _detName;}

Bool_t PndFsmAbsDet::doesPid ( ) const

inlineinherited

Definition at line 76 of file PndFsmAbsDet.h.

References PndFsmAbsDet::_doesPid.

Referenced by PndFastSim::Register(), and PndFastSim::smearTrack().

{return _doesPid;};

double PndFsmSttPid::gauss ( double  x, double  x0, double  s  )

private

Definition at line 142 of file PndFsmSttPid.cxx.

References exp(), and sqrt().

Referenced by respond().

{
  return (1.0/(sqrt(2.0*M_PI)*s))*
          exp(-(x-x0)*(x-x0)/(2.0*s*s));
}

void PndFsmSttPid::initParameters ( )

privatevirtual

Reimplemented from PndFsmAbsDet.

Definition at line 218 of file PndFsmSttPid.cxx.

References _dEdxRes, PndFsmAbsDet::_detName, PndFsmAbsDet::_efficiency, _ptmin, _thtMax, and _thtMin.

Referenced by PndFsmSttPid().

{
  _detName         = "SttPid";
  _thtMin          = 7.765;
  _thtMax          = 159.44;
  _ptmin            = 0.15;
  _dEdxRes         = 0.2; // 20% dEdx resolution
  _efficiency      = 1.0;

}

void PndFsmAbsDet::parseParameterList ( ArgList &  par )

protectedinherited

Definition at line 119 of file PndFsmAbsDet.cxx.

References CStrTok::GetFirst(), CStrTok::GetNext(), name, and PndFsmAbsDet::setParameter().

Referenced by PndFsmAbsDet::PndFsmAbsDet(), PndFsmCmpDet::PndFsmCmpDet(), PndFsmCombiDet::PndFsmCombiDet(), PndFsmDetTemplate::PndFsmDetTemplate(), PndFsmDrcBarrel::PndFsmDrcBarrel(), PndFsmDrcDisc::PndFsmDrcDisc(), PndFsmEffTracker::PndFsmEffTracker(), PndFsmEmcBarrel::PndFsmEmcBarrel(), PndFsmEmcBwCap::PndFsmEmcBwCap(), PndFsmEmcFS::PndFsmEmcFS(), PndFsmEmcFwCap::PndFsmEmcFwCap(), PndFsmEmcPid::PndFsmEmcPid(), PndFsmIdealPid::PndFsmIdealPid(), PndFsmMdcFS::PndFsmMdcFS(), PndFsmMdcTS::PndFsmMdcTS(), PndFsmMdtPid::PndFsmMdtPid(), PndFsmMvd::PndFsmMvd(), PndFsmMvd2::PndFsmMvd2(), PndFsmMvdPid::PndFsmMvdPid(), PndFsmRich::PndFsmRich(), PndFsmSimpleTracker::PndFsmSimpleTracker(), PndFsmSimpleVtx::PndFsmSimpleVtx(), PndFsmStt::PndFsmStt(), PndFsmSttPid(), and PndFsmTof::PndFsmTof().

{
  // ArgList is a std::list<std::string> of the form ("a=1","b=2","c=3","d=4")

  // Default values for the parameters

  CStrTok tokenizer;
  char csrc[200];

  if (par.size() != 0) {

    //cout <<"Parameters for detector <"<<_detName<<">"<<endl;
    for(ArgList::const_iterator argIt=par.begin(); argIt!=par.end();argIt++) {

      const char *src=argIt->data();
      strcpy(csrc,src);

      char* token = tokenizer.GetFirst(csrc,"=");
      string name(token);
      //cout <<"variable: -" << name << "-   ";
      token = tokenizer.GetNext("=");

      double value = atof(token);
      string strvalue(token);

      //cout <<"value: -"<<value<<"-   "<<endl;

      if (value!=0 || strvalue=="0" || strvalue=="0.0" || strvalue=="0.")
      {
        if (!setParameter(name,value))
          cout  <<" -W- (PndFsmAbsDet::parseParameterList) Unknown Parameter: <"<<name<< endl;
      }
      else
      {
        if (!setParameter(name,strvalue))
          cout  <<" -W- (PndFsmAbsDet::parseParameterList) Unknown Parameter: <"<<name<< endl;
      }

    }
  }

  /*
  if (par.size() != 0) {

    //cout <<"Parameters for detector <"<<_detName<<">"<<endl;
    for(ArgList::const_iterator argIt=par.begin(); argIt!=par.end();argIt++) {
      typedef boost::tokenizer<boost::char_separator<char> > tokenizer;
      boost::char_separator<char> sep("=");
      tokenizer tokens(*argIt, sep);
      tokenizer::iterator tok_iter = tokens.begin();
      string name(*tok_iter);
      ++tok_iter;
      double value=atof((*tok_iter).c_str());

      //cout  <<"<"<<name<<"> = "<<value<<endl;

      string name("efficiency");
      double value=0.99;
      if (!setParameter(name,value))
          cout  <<" -W- Unknown Parameter: <"<<name<<"> in PndFsmAbsDet::parseParameterList"  << endl;
    }
  }
  */
  //print(std::cout);

}

void PndFsmSttPid::print ( std::ostream &  o )

privatevirtual

Reimplemented from PndFsmAbsDet.

Definition at line 207 of file PndFsmSttPid.cxx.

References _dEdxRes, PndFsmAbsDet::_efficiency, _ptmin, _thtMax, _thtMin, and PndFsmAbsDet::detName().

{
  o <<"Parameters for detector <"<<detName()<<">"<<endl;
  o  <<"  _thtMin = "<<_thtMin<<endl;
  o  <<"  _thtMax = "<<_thtMax<<endl;
  o  <<"  _ptmin = "<<_ptmin<<endl;
  o  <<"  _dEdxRes = "<<_dEdxRes << " (rel)"<< endl;
  o  <<"  _efficiency = "<<_efficiency<<endl;
}

PndFsmResponse * PndFsmSttPid::respond ( PndFsmTrack *  t )

virtual

Accessors to contained information

Implements PndFsmAbsDet.

Definition at line 92 of file PndFsmSttPid.cxx.

References _dEdxRes, PndFsmAbsDet::_fdbPDG, PndFsmAbsDet::_rand, PndFsmTrack::charge(), compdEdx(), detected(), fabs(), gauss(), p, PndFsmTrack::p4(), PndFsmTrack::pdt(), PndFsmResponse::setDetected(), PndFsmResponse::setDetector(), PndFsmResponse::setLHElectron(), PndFsmResponse::setLHKaon(), PndFsmResponse::setLHMuon(), PndFsmResponse::setLHPion(), PndFsmResponse::setLHProton(), and PndFsmResponse::setSttdEdx().

{
  PndFsmResponse *result=new PndFsmResponse();

  result->setDetector(this);
  bool wasDetected=detected(t);
  result->setDetected(wasDetected);

  if (wasDetected && fabs(t->charge())>1e-8)
  {
     // now the dEdx information
    TParticlePDG* part = _fdbPDG->GetParticle(t->pdt());
    double mass = (part) ? part->Mass() : t->p4().M();
    double p=t->p4().Vect().Mag();

    // overall resolution for the dEdx measurement
    double dEdx = compdEdx(p,mass);
    double sig = _dEdxRes*dEdx;

    double m_e  = _fdbPDG->GetParticle(11)->Mass();
    double m_mu = _fdbPDG->GetParticle(13)->Mass();
    double m_pi = _fdbPDG->GetParticle(211)->Mass();
    double m_K  = _fdbPDG->GetParticle(321)->Mass();
    double m_p  = _fdbPDG->GetParticle(2212)->Mass();

    // compute the expected dEdx values for all particle types
    // we need these to determine the pdf's (gaussian around nominal dEdx with res sig)

    double dEdx_e  = compdEdx(p,m_e);
    double dEdx_mu = compdEdx(p,m_mu);
    double dEdx_pi = compdEdx(p,m_pi);
    double dEdx_K  = compdEdx(p,m_K);
    double dEdx_p  = compdEdx(p,m_p);

    double measdEdx = _rand->Gaus(dEdx,sig);

    result ->setSttdEdx(measdEdx,sig);

    if (dEdx_e)  result->setLHElectron( gauss(measdEdx,dEdx_e,sig) );
    if (dEdx_mu) result->setLHMuon( gauss(measdEdx,dEdx_mu,sig) );
    if (dEdx_pi) result->setLHPion( gauss(measdEdx,dEdx_pi,sig) );
    if (dEdx_K)  result->setLHKaon( gauss(measdEdx,dEdx_K,sig) );
    if (dEdx_p)  result->setLHProton(gauss(measdEdx,dEdx_p,sig) );
  }


  return result;
}

void PndFsmAbsDet::setName ( std::string &  name )

inlineinherited

Definition at line 75 of file PndFsmAbsDet.h.

References PndFsmAbsDet::_detName, and name.

Referenced by PndFsmDetFactory::create().

{_detName = name;};

bool PndFsmSttPid::setParameter ( std::string &  name, double  value  )

privatevirtual

Reimplemented from PndFsmAbsDet.

Definition at line 230 of file PndFsmSttPid.cxx.

References _dEdxRes, PndFsmAbsDet::_efficiency, _ptmin, _thtMax, and _thtMin.

{
  // *****************
  // include here all parameters which should be settable via tcl
  // *****************

  bool knownName=true;

  if (name == "thtMin")
    _thtMin=value;
  else
  if (name == "thtMax")
    _thtMax=value;
  else
  if (name == "ptmin")
        _ptmin=value;
  else
  if (name == "dEdxRes")
    _dEdxRes=value;
  else
  if (name == "efficiency")
    _efficiency=value;
  else
    knownName=false;

  return knownName;
}

bool PndFsmAbsDet::setParameter ( std::string &  name, std::string &  value  )

protectedvirtualinherited

Reimplemented in PndFsmDrcBarrel, PndFsmDrcDisc, PndFsmCombiDet, and PndFsmCmpDet.

Definition at line 99 of file PndFsmAbsDet.cxx.

{
  return false;
}

void PndFsmAbsDet::setStorePid ( Bool_t  doespid = kTRUE )

inlineinherited

Definition at line 77 of file PndFsmAbsDet.h.

References PndFsmAbsDet::_doesPid.

Referenced by PndFsmDetFactory::create().

{_doesPid=doespid;};

Member Data Documentation

double PndFsmSttPid::_dEdxRes

private

Definition at line 79 of file PndFsmSttPid.h.

Referenced by initParameters(), print(), respond(), and setParameter().

std::string PndFsmAbsDet::_detName

protectedinherited

Definition at line 92 of file PndFsmAbsDet.h.

Referenced by PndFsmAbsDet::detName(), PndFsmCmpDet::initParameters(), PndFsmIdealPid::initParameters(), PndFsmDetTemplate::initParameters(), PndFsmMvdPid::initParameters(), PndFsmSimpleTracker::initParameters(), PndFsmSimpleVtx::initParameters(), PndFsmTof::initParameters(), PndFsmMdtPid::initParameters(), PndFsmEmcPid::initParameters(), PndFsmCombiDet::initParameters(), PndFsmMvd2::initParameters(), PndFsmEmcBwCap::initParameters(), PndFsmEmcFS::initParameters(), initParameters(), PndFsmEmcFwCap::initParameters(), PndFsmMdcFS::initParameters(), PndFsmMdcTS::initParameters(), PndFsmRich::initParameters(), PndFsmDrcBarrel::initParameters(), PndFsmDrcDisc::initParameters(), PndFsmEmcBarrel::initParameters(), PndFsmMvd::initParameters(), PndFsmStt::initParameters(), PndFsmEffTracker::initParameters(), PndFsmAbsDet::initParameters(), PndFsmEmcPid::print(), PndFsmEmcFS::print(), PndFsmEmcFwCap::print(), PndFsmEmcBwCap::print(), PndFsmEmcBarrel::print(), and PndFsmAbsDet::setName().

Bool_t PndFsmAbsDet::_doesPid

protectedinherited

Definition at line 96 of file PndFsmAbsDet.h.

Referenced by PndFsmAbsDet::doesPid(), PndFsmAbsDet::PndFsmAbsDet(), and PndFsmAbsDet::setStorePid().

double PndFsmAbsDet::_efficiency

protectedinherited

Definition at line 93 of file PndFsmAbsDet.h.

Referenced by PndFsmDetTemplate::detected(), PndFsmSimpleTracker::detected(), PndFsmMdcFS::detected(), PndFsmMdcTS::detected(), PndFsmTof::detected(), PndFsmSimpleVtx::detected(), PndFsmEmcFS::detected(), PndFsmMvd::detected(), PndFsmEmcBwCap::detected(), PndFsmStt::detected(), PndFsmMvd2::detected(), PndFsmEmcFwCap::detected(), detected(), PndFsmMvdPid::detected(), PndFsmRich::detected(), PndFsmEmcBarrel::detected(), PndFsmDrcDisc::detected(), PndFsmMdtPid::detected(), PndFsmDrcBarrel::detected(), PndFsmEmcPid::detected(), PndFsmEffTracker::detected(), PndFsmSimpleTracker::initParameters(), PndFsmSimpleVtx::initParameters(), PndFsmDetTemplate::initParameters(), PndFsmMvdPid::initParameters(), PndFsmTof::initParameters(), PndFsmMdtPid::initParameters(), PndFsmEmcPid::initParameters(), PndFsmMvd2::initParameters(), PndFsmEmcFwCap::initParameters(), PndFsmRich::initParameters(), PndFsmMdcFS::initParameters(), initParameters(), PndFsmMdcTS::initParameters(), PndFsmEmcBwCap::initParameters(), PndFsmEmcFS::initParameters(), PndFsmEmcBarrel::initParameters(), PndFsmDrcDisc::initParameters(), PndFsmDrcBarrel::initParameters(), PndFsmStt::initParameters(), PndFsmMvd::initParameters(), PndFsmEffTracker::initParameters(), PndFsmSimpleTracker::print(), PndFsmSimpleVtx::print(), PndFsmDetTemplate::print(), PndFsmMvdPid::print(), PndFsmTof::print(), PndFsmMdtPid::print(), PndFsmEmcPid::print(), PndFsmMdcFS::print(), PndFsmEmcFS::print(), PndFsmMdcTS::print(), PndFsmEmcBwCap::print(), PndFsmMvd2::print(), print(), PndFsmEmcFwCap::print(), PndFsmRich::print(), PndFsmDrcDisc::print(), PndFsmEmcBarrel::print(), PndFsmDrcBarrel::print(), PndFsmStt::print(), PndFsmMvd::print(), PndFsmEffTracker::print(), PndFsmSimpleTracker::setParameter(), PndFsmMvdPid::setParameter(), PndFsmDetTemplate::setParameter(), PndFsmSimpleVtx::setParameter(), PndFsmMdtPid::setParameter(), PndFsmTof::setParameter(), PndFsmEmcPid::setParameter(), PndFsmEmcFS::setParameter(), PndFsmEmcBwCap::setParameter(), PndFsmMdcTS::setParameter(), PndFsmMvd2::setParameter(), setParameter(), PndFsmEmcFwCap::setParameter(), PndFsmMdcFS::setParameter(), PndFsmRich::setParameter(), PndFsmEmcBarrel::setParameter(), PndFsmDrcDisc::setParameter(), PndFsmDrcBarrel::setParameter(), PndFsmStt::setParameter(), PndFsmMvd::setParameter(), and PndFsmEffTracker::setParameter().

TDatabasePDG* PndFsmAbsDet::_fdbPDG

protectedinherited

Definition at line 95 of file PndFsmAbsDet.h.

Referenced by PndFsmRich::detected(), PndFsmDrcDisc::detected(), PndFsmDrcBarrel::detected(), PndFsmAbsDet::PndFsmAbsDet(), PndFsmStt::respond(), PndFsmTof::respond(), PndFsmMvd::respond(), respond(), PndFsmRich::respond(), PndFsmDrcDisc::respond(), PndFsmDrcBarrel::respond(), and PndFsmEffTracker::respond().

double PndFsmSttPid::_ptmin

private

Definition at line 78 of file PndFsmSttPid.h.

Referenced by detected(), initParameters(), print(), and setParameter().

TRandom3* PndFsmAbsDet::_rand

protectedinherited

Definition at line 94 of file PndFsmAbsDet.h.

Referenced by PndFsmDetTemplate::detected(), PndFsmSimpleTracker::detected(), PndFsmMvd::detected(), PndFsmEmcFwCap::detected(), PndFsmSimpleVtx::detected(), PndFsmMvd2::detected(), PndFsmMdcFS::detected(), PndFsmMdcTS::detected(), PndFsmEmcFS::detected(), PndFsmEmcBwCap::detected(), PndFsmStt::detected(), detected(), PndFsmMvdPid::detected(), PndFsmTof::detected(), PndFsmMdtPid::detected(), PndFsmEffTracker::detected(), PndFsmDrcBarrel::detected(), PndFsmEmcPid::detected(), PndFsmAbsDet::PndFsmAbsDet(), PndFsmRich::respond(), PndFsmMvdPid::respond(), PndFsmTof::respond(), respond(), PndFsmMvd2::respond(), PndFsmStt::respond(), PndFsmMvd::respond(), PndFsmDrcDisc::respond(), PndFsmDrcBarrel::respond(), and PndFsmEffTracker::respond().

double PndFsmSttPid::_thtMax

private

Definition at line 77 of file PndFsmSttPid.h.

Referenced by initParameters(), PndFsmSttPid(), print(), and setParameter().

double PndFsmSttPid::_thtMin

private

Definition at line 76 of file PndFsmSttPid.h.

Referenced by detected(), initParameters(), PndFsmSttPid(), print(), and setParameter().

---

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

• PndFsmSttPid.h
• PndFsmSttPid.cxx