PndFsmMvd2 Class Reference

#include <PndFsmMvd2.h>

Inheritance diagram for PndFsmMvd2:

Public Member Functions
	PndFsmMvd2 ()
	PndFsmMvd2 (ArgList &par)
virtual	~PndFsmMvd2 ()
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 Types
enum	{ nPidType = 5 }
enum	PidType {   electron, muon, pion, kaon,   proton }

Private Member Functions
bool	detected (PndFsmTrack *t) const
double	dp (PndFsmTrack *t) const
double	dphi (PndFsmTrack *t) const
double	dtheta (PndFsmTrack *t) const
void	initParameters ()
void	print (std::ostream &o)
bool	setParameter (std::string &name, double value)
double	MeanEnergyLoss (PidType particle)
double	LandauGaus (double s_mpv, double width1, double width2)
double	Likelihood (PidType particle)
double	mpv (PidType particle)
double	width1 (PidType particle)
double	width2 (PidType particle)

Private Attributes
double	_thtMin
double	_thtMax
double	_radiationLength
double	_pmin
double	_vtxRes
double	_pRes
double	_dEdxResMulti
double	_phiRes
double	_thetaRes
double	_momentum
double	_energyloss

Detailed Description

Definition at line 41 of file PndFsmMvd2.h.

Member Enumeration Documentation

anonymous enum

private

Enumerator
nPidType

Definition at line 91 of file PndFsmMvd2.h.

{nPidType = 5};

enum PndFsmMvd2::PidType

private

Enumerator
electron
muon
pion
kaon
proton

Definition at line 92 of file PndFsmMvd2.h.

               {
    electron,
    muon,
    pion,
    kaon,
    proton
  };

Constructor & Destructor Documentation

PndFsmMvd2::PndFsmMvd2

(

)

Default ctor.

Definition at line 58 of file PndFsmMvd2.cxx.

References _thtMax, _thtMin, and initParameters().

{
  initParameters();

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

PndFsmMvd2::PndFsmMvd2

(

ArgList &

par

)

Definition at line 67 of file PndFsmMvd2.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);
}

PndFsmMvd2::~PndFsmMvd2 ( )

virtual

Destructor

Definition at line 84 of file PndFsmMvd2.cxx.

{
}

Member Function Documentation

bool PndFsmMvd2::detected ( PndFsmTrack *  t ) const

private

Definition at line 146 of file PndFsmMvd2.cxx.

References PndFsmAbsDet::_efficiency, _pmin, PndFsmAbsDet::_rand, _thtMin, PndFsmTrack::charge(), PndFsmTrack::hitMapResponse(), PndFsmTrack::hitMapValid(), FsmDetEnum::Mvd, p, PndFsmTrack::p4(), and theta.

Referenced by respond().

{
  if (t->hitMapValid()) {
    return t->hitMapResponse(FsmDetEnum::Mvd);
  } else {
    double theta = t->p4().Theta();
    double p=t->p4().Vect().Mag();
    double charge=t->charge();
    return ( charge!=0.0 && theta>=_thtMin && theta<=_thtMax && p>_pmin  && _rand->Uniform()<=_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(), PndFsmSttPid::print(), PndFsmMdcFS::print(), PndFsmRich::print(), PndFsmMdcTS::print(), 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 PndFsmMvd2::dp ( PndFsmTrack *  t ) const

private

Definition at line 159 of file PndFsmMvd2.cxx.

References _pRes, p, and PndFsmTrack::p4().

Referenced by respond().

{
  double p=t->p4().Vect().Mag();
  double Dp=_pRes*p;   //10% smearing
  return ( Dp );  //to be refined
}

double PndFsmMvd2::dphi ( PndFsmTrack *  t ) const

private

Definition at line 167 of file PndFsmMvd2.cxx.

References _phiRes.

Referenced by respond().

{
  double Dphi=_phiRes*M_PI/180.0;   
  return Dphi; //to be refined
}

double PndFsmMvd2::dtheta ( PndFsmTrack *  t ) const

private

Definition at line 174 of file PndFsmMvd2.cxx.

References _thetaRes.

Referenced by respond().

{
  double dt=_thetaRes*M_PI/180.0;   
  return dt; //to be refined 
}

void PndFsmMvd2::initParameters ( )

privatevirtual

Reimplemented from PndFsmAbsDet.

Definition at line 197 of file PndFsmMvd2.cxx.

References _dEdxResMulti, PndFsmAbsDet::_detName, PndFsmAbsDet::_efficiency, _phiRes, _pmin, _pRes, _radiationLength, _thetaRes, _thtMax, _thtMin, _vtxRes, FsmDetEnum::Mvd2, and FsmDetName::name().

Referenced by PndFsmMvd2().

{
  _detName         = FsmDetName::name(FsmDetEnum::Mvd2);
  _thtMin          = 5.0;
  _thtMax          = 160.0;
  _radiationLength = 0.0;
  _pmin            = 0.0;
  _vtxRes          = 100e-4;  //vtx res = 100 um 
  _pRes            =0.1;//10% smearing
  _dEdxResMulti    =1.00; //normal energy loss resolution
  _phiRes          =0.01; //0.01 degree smearing
  _thetaRes        =0.01; //0.01 degree smearing
  _efficiency      =1.0; 
}

double PndFsmMvd2::LandauGaus ( double  s_mpv, double  width1, double  width2  )

private

Definition at line 270 of file PndFsmMvd2.cxx.

References CAMath::Abs(), c1, c2, f1, f2, h, i, and x.

Referenced by Likelihood().

                                                                            {
  // this is the adapted TF1::Integral function from
  // ROOT 5.14. GOTOs have been removed and interval 
  // division has been modified to maximize performance

  if (widthone<=0) 
    return TMath::Gaus(s_mpv, 0, widthtwo, true);
  else if (widthtwo<=0)
    return TMath::Landau(s_mpv, 0, widthone, true);
  else {
    static double x[12] = { 0.96028985649753623,  0.79666647741362674,
                            0.52553240991632899,  0.18343464249564980,
                            0.98940093499164993,  0.94457502307323258,
                            0.86563120238783174,  0.75540440835500303,
                            0.61787624440264375,  0.45801677765722739,
                            0.28160355077925891,  0.09501250983763744};

    static double w[12] = { 0.10122853629037626,  0.22238103445337447,
                            0.31370664587788729,  0.36268378337836198,
                            0.02715245941175409,  0.06225352393864789,
                            0.09515851168249278,  0.12462897125553387,
                            0.14959598881657673,  0.16915651939500254,
                            0.18260341504492359,  0.18945061045506850};

    double h, bb, aa, c1, c2, u, s8, s16, f1, f2;
    double xx;
    bool redo=true;
    int i;

    h = 0;
    aa = -5.0*widthtwo;
    bb = -2.5*widthtwo;

    do {
      c1 = 0.5*(bb+aa);
      c2 = 0.5*(bb-aa);

      s8 = 0;
      for (i=0;i<4;i++) {
        u  = c2*x[i];
        xx = c1+u;
        f1 = TMath::Landau(s_mpv+xx, 0, widthone, true)*TMath::Gaus(xx, 0, widthtwo, true);
        xx = c1-u;
        f2 = TMath::Landau(s_mpv+xx, 0, widthone, true)*TMath::Gaus(xx, 0, widthtwo, true);
        s8+= w[i]*(f1 + f2);
      }
      s16 = 0;
      for (i=4;i<12;i++) {
        u   = c2*x[i];
        xx  = c1+u;
        f1  = TMath::Landau(s_mpv+xx, 0, widthone, true)*TMath::Gaus(xx, 0, widthtwo, true);
        xx  = c1-u;
        f2  = TMath::Landau(s_mpv+xx, 0, widthone, true)*TMath::Gaus(xx, 0, widthtwo, true);
        s16+= w[i]*(f1 + f2);
      }
      s16 = c2*s16;
      if (TMath::Abs(s16-c2*s8) <= 1e-12*(s16+1) ) {
        aa =bb;
        bb+=2*c2*1.5;
        if (bb>=5*widthtwo) {
          bb=5*widthtwo;
          redo=false;
        }
        h += s16;
      } else 
        bb = c1;
    } while (redo);
    return h;
  }
}

double PndFsmMvd2::Likelihood ( PidType  particle )

private

Definition at line 341 of file PndFsmMvd2.cxx.

References _dEdxResMulti, _energyloss, LandauGaus(), MeanEnergyLoss(), mpv(), width1(), and width2().

Referenced by respond().

                                          {
  return LandauGaus( _energyloss - MeanEnergyLoss(part) - mpv(part), width1(part), width2(part)*_dEdxResMulti );
}

double PndFsmMvd2::MeanEnergyLoss ( PidType  particle )

private

Definition at line 256 of file PndFsmMvd2.cxx.

References _momentum, c, electron, log(), Mass, and nPidType.

Referenced by Likelihood(), and respond().

                                              {
  //Calculate the lower boundary of the energy loss distribution.

  //[GeV]
  static float eb=0.14e-6;
  //[m/s]
  static float c=2.99792458e8;
  //[GeV/c**2]
  static float Mass[nPidType]={ 0.511e-3, 0.1058, 0.1396, 0.4937, 0.9383 };

  double sqrBeta=1/(1+pow(Mass[part]/_momentum,2));
  return 4.9312e-05 * (log(2*Mass[electron]*c*c/eb*sqrBeta/(1-sqrBeta))-sqrBeta)/sqrBeta;
};

double PndFsmMvd2::mpv ( PidType  particle )

private

Definition at line 345 of file PndFsmMvd2.cxx.

References _momentum, c1, electron, kaon, muon, pion, proton, x, and x0.

Referenced by Likelihood(), and respond().

                                   {
  double x=_momentum;
  double x0=0;
  double x1=0;
  double c0=0;
  double c1=0;
  double d1=0;
  double a3=0;
  double a4=0;
  double a5=0;

  switch(part) {
  case proton:
    x0 = 0.45;
    x1 = 1.3;
    c0 = 0.383451e-03;
    c1 = -0.126986e-03;
    d1 = -5.21351e-06;
    a3 = -3.05821;
    a4 = 24.6356;
    a5 = -68.632;
    break;
  case kaon:
    x0 = 0.25;
    x1 = 1.05;
    c0 = 0.326259e-03;
    c1 = -7.68052e-05;
    d1 = 1.51033e-05;
    a3 = -19.6615;
    a4 = 264.382;
    a5 = -1238.09;
    break;
  case pion:
    x0 = 0.1;
    x1 = 1.0;
    c0 = 0.274692e-03;
    c1 = 3.2571e-05;
    d1 = 9.16527e-06;
    a5 = -6624.05;
    break;
  case muon:
    x0 = 0.15;
    x1 = 1.15;
    a3 = 4.33244;
    a4 = -107.686;
    a5 = 699.522;
    c0 = 0.248749e-03;
    c1 = 6.57118e-05;
    d1 = -4.09447e-06;
    break;
  case electron:
    x1 = 1.20;
    c0 = 2.93999e-03;
    c1 = 1.76792e-05;
    break;
  }
  if (x>=x1)
    return c0+c1*(x1-x0)+d1*(x-x1);
  if (x>=x0)
    return c0+c1*(x-x0);
  else
    return c0+c1*(x-x0)+pow(x0-x,3)*(a3+(x0-x)*(a4+(x0-x)*a5));
}

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(), PndFsmMvdPid::PndFsmMvdPid(), PndFsmRich::PndFsmRich(), PndFsmSimpleTracker::PndFsmSimpleTracker(), PndFsmSimpleVtx::PndFsmSimpleVtx(), PndFsmStt::PndFsmStt(), PndFsmSttPid::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 PndFsmMvd2::print ( std::ostream &  o )

privatevirtual

Reimplemented from PndFsmAbsDet.

Definition at line 181 of file PndFsmMvd2.cxx.

References _dEdxResMulti, PndFsmAbsDet::_efficiency, _phiRes, _pmin, _pRes, _radiationLength, _thetaRes, _thtMax, _thtMin, _vtxRes, and PndFsmAbsDet::detName().

{
  o <<"Parameters for detector <"<<detName()<<">"<<endl;
  o  <<"  _thtMin = "<<_thtMin<<endl; 
  o  <<"  _thtMax = "<<_thtMax<<endl; 
  o  <<"  _radiationLength = "<<_radiationLength<<endl; 
  o  <<"  _pmin = "<<_pmin<<endl; 
  o  <<"  _pRes = "<<100*_pRes<< "\%" << endl; 
  o  <<"  _dEdxResMulti = "<<100*_dEdxResMulti<< "\%" << endl; 
  o  <<"  _vtxRes = "<<_vtxRes*1e5<<"um "<<endl;
  o  <<"  _phiRes = "<<_phiRes << " degree" << endl;
  o  <<"  _thetaRes = "<<_thetaRes <<" degree" << endl; 
  o  <<"  _efficiency = "<<_efficiency<<endl; 
}

PndFsmResponse * PndFsmMvd2::respond ( PndFsmTrack *  t )

virtual

Accessors to contained information

Implements PndFsmAbsDet.

Definition at line 93 of file PndFsmMvd2.cxx.

References _dEdxResMulti, _energyloss, _momentum, PndFsmAbsDet::_rand, _vtxRes, PndFsmTrack::charge(), detected(), dp(), dphi(), dtheta(), electron, fabs(), kaon, Likelihood(), MeanEnergyLoss(), mpv(), muon, PndFsmTrack::p4(), PndFsmTrack::pdt(), pion, proton, PndFsmResponse::setDetected(), PndFsmResponse::setDetector(), PndFsmResponse::setdp(), PndFsmResponse::setdphi(), PndFsmResponse::setdtheta(), PndFsmResponse::setdV(), PndFsmResponse::setLHElectron(), PndFsmResponse::setLHKaon(), PndFsmResponse::setLHMuon(), PndFsmResponse::setLHPion(), PndFsmResponse::setLHProton(), PndFsmResponse::setMvddEdx(), width1(), and width2().

{
  PndFsmResponse *result=new PndFsmResponse();
  
  result->setDetector(this);
  bool wasDetected=detected(t);
  result->setDetected(wasDetected);
  
  if (wasDetected && fabs(t->charge())>1e-8)
  {
    //select particle
    PidType part;
    switch(abs(t->pdt())) {
      case 2212:part=proton; break;
      case 321: part=kaon; break;
      case 211: part=pion; break;
      case 13:  part=muon; break;
      case 11:  part=electron; break;
    }

    //build random energy loss
    _momentum = t->p4().Vect().Mag();
    _energyloss = MeanEnergyLoss(part) + mpv(part) + _rand->Landau(0, width1(part)) + _rand->Gaus(0, width2(part))*_dEdxResMulti;

    //store random energy loss
    result->setMvddEdx(_energyloss/1e-3);

    //store energy loss 
    if (_momentum>=0 && _momentum<=2.5) {
      result->setLHElectron(Likelihood(electron));
      result->setLHMuon(Likelihood(muon));
      result->setLHPion(Likelihood(pion));
      result->setLHKaon(Likelihood(kaon));
      result->setLHProton(Likelihood(proton));
    } else {
      result->setLHElectron(1);
      result->setLHMuon(1);
      result->setLHPion(1);
      result->setLHKaon(1);
      result->setLHProton(1);
    }

    result->setdV(_vtxRes, _vtxRes, _vtxRes);
    result->setdp(dp(t));
    result->setdphi(dphi(t));
    result->setdtheta(dtheta(t));
  }

  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 PndFsmMvd2::setParameter ( std::string &  name, double  value  )

privatevirtual

Reimplemented from PndFsmAbsDet.

Definition at line 213 of file PndFsmMvd2.cxx.

References _dEdxResMulti, PndFsmAbsDet::_efficiency, _phiRes, _pmin, _pRes, _radiationLength, _thetaRes, _thtMax, _thtMin, and _vtxRes.

{
  // *****************
  // 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 == "radiationLength")
    _radiationLength=value;
  else
  if (name == "pmin")
    _pmin=value;
  else
  if (name == "vtxRes")
    _vtxRes=value;
  else
  if (name == "pRes")
    _pRes=value;
  else
  if (name == "dEdxResMulti")
    _dEdxResMulti=value;
  else
  if (name == "thetaRes")
    _thetaRes=value;
  else
  if (name == "phiRes")
    _phiRes=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;};

double PndFsmMvd2::width1 ( PidType  particle )

private

Definition at line 409 of file PndFsmMvd2.cxx.

References _momentum, electron, kaon, muon, pion, proton, and x.

Referenced by Likelihood(), and respond().

                                      {
  double x=_momentum;
  switch(part) {
  case proton:
    if (x>=1.10)
      return +3.81174e-04+x*(-2.25108e-04+x*+5.45154e-05);
    else
      return -5.28145e-05+x*(+8.29883e-04+x*-5.35972e-04);
    break;
  case kaon:
    if (x>=1.05)
      return +2.61134e-04+x*(-1.30818e-04+x*+3.44165e-05);
    else
      return +3.41858e-04+x*(-3.21115e-04+x*+1.37459e-04);
    break;
  case pion:
    if (x>=1.00)
      return +1.88718e-04+x*(-6.38948e-05+x*+1.78590e-05);
    else
      return +1.82872e-04+x*(-1.28373e-04+x*+8.01459e-05);
    break;
  case muon:
    if (x>=1.20)
      return +1.06142e-04+x*(+3.68777e-05+x*-1.00190e-05);
    else
      return +1.89374e-04+x*(-1.46441e-04+x*+9.10813e-05);
    break;
  case electron:
    if (x>1.2) x=1.2;
    // electrons are constant for momentum > 1.2GeV
    return +1.27955e-04+x*(-3.15732e-06+x*+9.64736e-06);
    break;
  }
  return 0;
}

double PndFsmMvd2::width2 ( PidType  particle )

private

Definition at line 445 of file PndFsmMvd2.cxx.

References _momentum, electron, kaon, muon, pion, proton, and x.

Referenced by Likelihood(), and respond().

                                      {
  double x=_momentum;
  switch(part) {
  case proton:
    if (x>=1.10)
      return +6.41067e-04+x*(-3.82507e-04+x*+9.03732e-05);
    else
      return +6.40328e-04-3.21725e-04*x+3.17708e-05*pow(x,-3);
    break;
  case kaon:
    if (x>=1.05)
      return +2.22504e-04+x*(-6.40051e-06+x*+2.14434e-06);
    else
      return +3.86684e-04-1.61873e-04*x+7.76586e-06*pow(x,-3);
    break;
  case pion:
    if (x>=1.00)
      return +1.32999e-04+x*(+1.19714e-04+x*-3.53302e-05);
    else
      return +2.21603e-04-3.21357e-06*x+4.64793e-06*pow(x,-2);
    break;
  case muon:
   if (x>=1.20)
      return +7.84582e-05+x*(+1.88988e-04+x*-5.49637e-05);
    else
      return +1.67388e-04+5.67991e-05*x+3.42702e-06*pow(x,-2);
    break;
  case electron:
    if (x>1.2) x=1.2;
    // electrons are constant for momentum > 1.2GeV
    return +4.08849e-04-3.56548e-05*x+1.84825e-08*pow(x,-3);
    break;
  }
  return 0;
}

Member Data Documentation

double PndFsmMvd2::_dEdxResMulti

private

Definition at line 86 of file PndFsmMvd2.h.

Referenced by initParameters(), Likelihood(), 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(), initParameters(), PndFsmEmcBwCap::initParameters(), PndFsmEmcFS::initParameters(), PndFsmSttPid::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(), PndFsmSttPid::detected(), PndFsmEmcFS::detected(), detected(), PndFsmMvd::detected(), PndFsmMdcFS::detected(), PndFsmStt::detected(), PndFsmEmcFwCap::detected(), PndFsmMdcTS::detected(), PndFsmSimpleVtx::detected(), PndFsmTof::detected(), PndFsmEmcBwCap::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(), PndFsmEmcBwCap::initParameters(), PndFsmEmcFS::initParameters(), initParameters(), PndFsmEmcFwCap::initParameters(), PndFsmRich::initParameters(), PndFsmMdcFS::initParameters(), PndFsmSttPid::initParameters(), PndFsmMdcTS::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(), print(), PndFsmSttPid::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(), PndFsmRich::setParameter(), PndFsmMdcFS::setParameter(), PndFsmEmcFS::setParameter(), PndFsmEmcBwCap::setParameter(), setParameter(), PndFsmMdcTS::setParameter(), PndFsmSttPid::setParameter(), PndFsmEmcFwCap::setParameter(), PndFsmEmcBarrel::setParameter(), PndFsmDrcDisc::setParameter(), PndFsmDrcBarrel::setParameter(), PndFsmStt::setParameter(), PndFsmMvd::setParameter(), and PndFsmEffTracker::setParameter().

double PndFsmMvd2::_energyloss

private

Definition at line 108 of file PndFsmMvd2.h.

Referenced by Likelihood(), and respond().

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(), PndFsmSttPid::respond(), PndFsmRich::respond(), PndFsmDrcDisc::respond(), PndFsmDrcBarrel::respond(), and PndFsmEffTracker::respond().

double PndFsmMvd2::_momentum

private

Definition at line 107 of file PndFsmMvd2.h.

Referenced by MeanEnergyLoss(), mpv(), respond(), width1(), and width2().

double PndFsmMvd2::_phiRes

private

Definition at line 87 of file PndFsmMvd2.h.

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

double PndFsmMvd2::_pmin

private

Definition at line 83 of file PndFsmMvd2.h.

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

double PndFsmMvd2::_pRes

private

Definition at line 85 of file PndFsmMvd2.h.

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

double PndFsmMvd2::_radiationLength

private

Definition at line 82 of file PndFsmMvd2.h.

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

TRandom3* PndFsmAbsDet::_rand

protectedinherited

Definition at line 94 of file PndFsmAbsDet.h.

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

double PndFsmMvd2::_thetaRes

private

Definition at line 88 of file PndFsmMvd2.h.

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

double PndFsmMvd2::_thtMax

private

Definition at line 81 of file PndFsmMvd2.h.

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

double PndFsmMvd2::_thtMin

private

Definition at line 80 of file PndFsmMvd2.h.

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

double PndFsmMvd2::_vtxRes

private

Definition at line 84 of file PndFsmMvd2.h.

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

---

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

• PndFsmMvd2.h
• PndFsmMvd2.cxx