PndFsmCmpDet Class Reference

#include <PndFsmCmpDet.h>

Inheritance diagram for PndFsmCmpDet:

Public Member Functions
	PndFsmCmpDet ()
	PndFsmCmpDet (ArgList &par)
virtual	~PndFsmCmpDet ()
virtual PndFsmResponse *	respond (PndFsmTrack *t)
bool	AddDetector (std::string name, std::string params="")
bool	SubtractDetector (std::string name, std::string params="")
virtual void	print (std::ostream &o)
const std::string &	detName ()
void	setName (std::string &name)
Bool_t	doesPid () const
void	setStorePid (Bool_t doespid=kTRUE)

Protected Member Functions
void	parseParameterList (ArgList &par)

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

Private Types
typedef std::list< PndFsmAbsDet * >	FsmAbsDetList

Private Member Functions
bool	setParameter (std::string &name, std::string &value)
bool	setParameter (std::string &name, double value)
void	readParameters ()
void	initParameters ()
double	eval (TSpline3 *spline, double theta)

Private Attributes
FsmAbsDetList	fDetList
FsmAbsDetList	fSubtractDetList
PndFsmDetFactory	fDetFact
TString	_parFileName
bool	_parFile
double	_d0ResMulti
double	_z0ResMulti
double	_thtResMulti
double	_phiResMulti
double	_momResMulti
double	_pidLhMulti
std::map< int, TSpline3 * >	_d0
std::map< int, TSpline3 * >	_z0
std::map< int, TSpline3 * >	_tht
std::map< int, TSpline3 * >	_phi
std::map< int, TSpline3 * >	_mom
std::map< int, TParameter < double > * >	_mom0
TParameter< double > *	_tht0
TParameter< double > *	_tht1
std::map< int, TF1 * >	_d0Scale
std::map< int, TF1 * >	_z0Scale
std::map< int, TF1 * >	_thtScale
std::map< int, TF1 * >	_phiScale
std::map< int, TF1 * >	_momScale

Detailed Description

Definition at line 47 of file PndFsmCmpDet.h.

Member Typedef Documentation

typedef std::list<PndFsmAbsDet*> PndFsmCmpDet::FsmAbsDetList

private

Definition at line 69 of file PndFsmCmpDet.h.

Constructor & Destructor Documentation

PndFsmCmpDet::PndFsmCmpDet

(

)

Definition at line 93 of file PndFsmCmpDet.cxx.

References initParameters(), and readParameters().

                           {
  initParameters();
  readParameters();
}

PndFsmCmpDet::PndFsmCmpDet

(

ArgList &

par

)

Definition at line 62 of file PndFsmCmpDet.cxx.

References _parFileName, AddDetector(), initParameters(), PndFsmAbsDet::parseParameterList(), and readParameters().

                                       {
  initParameters();
  parseParameterList(par);
  std::cout<<" -I- (PndFsmCmpDet::PndFsmCmpDet) - Using parameter file '"<<_parFileName<<"'"<<endl;
  readParameters();

  // when PndFsmCmpDet has been created by the detector
  // factory, some arbitrary detector set is added

  std::cout<<" -I- (PndFsmCmpDet::PndFsmCmpDet) - Creating arbitrary detector set"<<endl;

  //EM Calorimeters w/ default parameters
  AddDetector("EmcBarrel");
  AddDetector("EmcFwCap");
  AddDetector("EmcBwCap");
  AddDetector("EmcFS");

  //Tracking, Vertexing, dE/dx
  AddDetector("Mvd2");
  AddDetector("Stt");
  AddDetector("MdcTS");
  AddDetector("MdcFS");
//AddDetector("Tpc");

  //PID detectors
  AddDetector("DrcBarrel");
  AddDetector("DrcDisc");
  AddDetector("Rich");
  AddDetector("Tof");
}

PndFsmCmpDet::~PndFsmCmpDet ( )

virtual

Definition at line 102 of file PndFsmCmpDet.cxx.

                            {
}

Member Function Documentation

bool PndFsmCmpDet::AddDetector	(	std::string	name,
		std::string	params = ""
	)

Definition at line 349 of file PndFsmCmpDet.cxx.

References PndFsmDetFactory::create(), fDetFact, and fDetList.

Referenced by PndFsmCmpDet().

                                                         {
  PndFsmAbsDet *det=fDetFact.create(name,params);
  if (det)
    fDetList.push_back(det);
  return det;
}

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

inlineinherited

Definition at line 74 of file PndFsmAbsDet.h.

References PndFsmAbsDet::_detName.

Referenced by PndFastSim::AddDetector(), PndFsmIdealPid::print(), PndFsmSimpleTracker::print(), PndFsmSimpleVtx::print(), PndFsmMvdPid::print(), PndFsmDetTemplate::print(), PndFsmMdtPid::print(), PndFsmTof::print(), PndFsmAbsDet::print(), PndFsmSttPid::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 PndFsmCmpDet::eval ( TSpline3 *  spline, double  theta  )

private

Definition at line 535 of file PndFsmCmpDet.cxx.

References xmax, and xmin.

Referenced by respond().

                                               {
  double xmin=s->GetXmin();
  double xmax=s->GetXmax();
  if (x<xmin)
    return s->Eval(xmin)+(x-xmin)*s->Derivative(xmin);
  else if (x>xmax)
    return s->Eval(xmax)+(x-xmax)*s->Derivative(xmax);
  else
    return s->Eval(x);
}

void PndFsmCmpDet::initParameters ( )

privatevirtual

Reimplemented from PndFsmAbsDet.

Definition at line 514 of file PndFsmCmpDet.cxx.

References _d0ResMulti, PndFsmAbsDet::_detName, _mom0, _momResMulti, _parFile, _parFileName, _phiResMulti, _pidLhMulti, _tht0, _tht1, _thtResMulti, and _z0ResMulti.

Referenced by PndFsmCmpDet().

                                  {
  _detName = "CmpDet";
  _parFileName = "$VMCWORKDIR/fastsim/cmpdetparams.root";
  _parFile=false;
  _d0ResMulti=1.0;
  _z0ResMulti=1.0;
  _thtResMulti=1.0;
  _phiResMulti=1.0;
  _momResMulti=1.0;
  _pidLhMulti=1.0;
  _mom0[11]=0;
  _mom0[13]=0;
  _mom0[211]=0;
  _mom0[321]=0;
  _mom0[2212]=0;
  _tht0=0;
  _tht1=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(), 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::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 PndFsmAbsDet::print ( std::ostream &  o )

virtualinherited

Reimplemented in PndFsmEffTracker, PndFsmMvd, PndFsmStt, PndFsmDrcBarrel, PndFsmDrcDisc, PndFsmEmcBarrel, PndFsmEmcBwCap, PndFsmEmcFS, PndFsmEmcFwCap, PndFsmMdcFS, PndFsmMdcTS, PndFsmMvd2, PndFsmRich, PndFsmSttPid, PndFsmEmcPid, PndFsmMdtPid, PndFsmTof, PndFsmDetTemplate, PndFsmMvdPid, PndFsmSimpleTracker, PndFsmSimpleVtx, and PndFsmIdealPid.

Definition at line 112 of file PndFsmAbsDet.cxx.

References PndFsmAbsDet::detName().

Referenced by PndFastSim::Init().

{
  o<<"Detector:"<<detName()<<endl;
}

void PndFsmCmpDet::readParameters ( )

private

Definition at line 414 of file PndFsmCmpDet.cxx.

References _d0, _d0Scale, _mom, _mom0, _momScale, _parFile, _parFileName, _phi, _phiScale, _tht, _tht0, _tht1, _thtScale, _z0, _z0Scale, f, for(), i, and if().

Referenced by PndFsmCmpDet().

                                  {

  TFile f(_parFileName);

  if (f.IsZombie()) {
    cout <<" -W-  (PndFsmCmpDet::readParameters) - file "<<_parFileName
         <<" doesn't exist. Using constant vertex reso"<<endl;
  } else {
    // electrons
    _d0[11]=(TSpline3*)f.Get("d0E-");
    _z0[11]=(TSpline3*)f.Get("z0E-");
    _tht[11]=(TSpline3*)f.Get("thtE-");
    _phi[11]=(TSpline3*)f.Get("phiE-");
    _mom[11]=(TSpline3*)f.Get("momE-");
    if (!_mom0[11]) _mom0[11]=(TParameter<double>*)f.Get("mom0E-");
    _d0Scale[11]=(TF1*)f.Get("d0ScaleE-");
    _z0Scale[11]=(TF1*)f.Get("z0ScaleE-");
    _thtScale[11]=(TF1*)f.Get("thtScaleE-");
    _phiScale[11]=(TF1*)f.Get("phiScaleE-");
    _momScale[11]=(TF1*)f.Get("momScaleE-");
    // muons
    _d0[13]=(TSpline3*)f.Get("d0Mu+");
    _z0[13]=(TSpline3*)f.Get("z0Mu+");
    _tht[13]=(TSpline3*)f.Get("thtMu+");
    _phi[13]=(TSpline3*)f.Get("phiMu+");
    if (!_mom0[13]) _mom[13]=(TSpline3*)f.Get("momMu+");
    _mom0[13]=(TParameter<double>*)f.Get("mom0Mu+");
    _d0Scale[13]=(TF1*)f.Get("d0ScaleMu+");
    _z0Scale[13]=(TF1*)f.Get("z0ScaleMu+");
    _thtScale[13]=(TF1*)f.Get("thtScaleMu+");
    _phiScale[13]=(TF1*)f.Get("phiScaleMu+");
    _momScale[13]=(TF1*)f.Get("momScaleMu+");
    // pions
    _d0[211]=(TSpline3*)f.Get("d0Pi-");
    _z0[211]=(TSpline3*)f.Get("z0Pi-");
    _tht[211]=(TSpline3*)f.Get("thtPi-");
    _phi[211]=(TSpline3*)f.Get("phiPi-");
    _mom[211]=(TSpline3*)f.Get("momPi-");
    if (!_mom0[211]) _mom0[211]=(TParameter<double>*)f.Get("mom0Pi-");
    _d0Scale[211]=(TF1*)f.Get("d0ScalePi-");
    _z0Scale[211]=(TF1*)f.Get("z0ScalePi-");
    _thtScale[211]=(TF1*)f.Get("thtScalePi-");
    _phiScale[211]=(TF1*)f.Get("phiScalePi-");
    _momScale[211]=(TF1*)f.Get("momScalePi-");
    // kaons
    _d0[321]=(TSpline3*)f.Get("d0K-");
    _z0[321]=(TSpline3*)f.Get("z0K-");
    _tht[321]=(TSpline3*)f.Get("thtK-");
    _phi[321]=(TSpline3*)f.Get("phiK-");
    _mom[321]=(TSpline3*)f.Get("momK-");
    if (!_mom0[321]) _mom0[321]=(TParameter<double>*)f.Get("mom0K-");
    _d0Scale[321]=(TF1*)f.Get("d0ScaleK-");
    _z0Scale[321]=(TF1*)f.Get("z0ScaleK-");
    _thtScale[321]=(TF1*)f.Get("thtScaleK-");
    _phiScale[321]=(TF1*)f.Get("phiScaleK-");
    _momScale[321]=(TF1*)f.Get("momScaleK-");
    // protons
    _d0[2212]=(TSpline3*)f.Get("d0P+");
    _z0[2212]=(TSpline3*)f.Get("z0P+");
    _tht[2212]=(TSpline3*)f.Get("thtP+");
    _phi[2212]=(TSpline3*)f.Get("phiP+");
    _mom[2212]=(TSpline3*)f.Get("momP+");
    if (!_mom0[2212]) _mom0[2212]=(TParameter<double>*)f.Get("mom0P+");
    _d0Scale[2212]=(TF1*)f.Get("d0ScaleP+");
    _z0Scale[2212]=(TF1*)f.Get("z0ScaleP+");
    _thtScale[2212]=(TF1*)f.Get("thtScaleP+");
    _phiScale[2212]=(TF1*)f.Get("phiScaleP+");
    _momScale[2212]=(TF1*)f.Get("momScaleP+");

    if (!_tht0) _tht0=(TParameter<double>*)f.Get("tht0");
    if (!_tht1) _tht1=(TParameter<double>*)f.Get("tht1");

    f.Close();
    // now check that all parameters have been loaded correctly
    _parFile=true;
    std::map<int, TSpline3*>::iterator i;
    std::map<int, TParameter<double>*>::iterator j;
    std::map<int, TF1*>::iterator k;

    for (i=_d0.begin(); i!=_d0.end(); i++) _parFile = _parFile && i->second;
    for (i=_z0.begin(); i!=_z0.end(); i++) _parFile = _parFile && i->second;
    for (i=_tht.begin(); i!=_tht.end(); i++) _parFile = _parFile && i->second;
    for (i=_phi.begin(); i!=_phi.end(); i++) _parFile = _parFile && i->second;
    for (i=_mom.begin(); i!=_mom.end(); i++) _parFile = _parFile && i->second;

    for (j=_mom0.begin(); j!=_mom0.end(); j++) _parFile = _parFile && j->second;
    _parFile = _parFile && _tht0 && _tht1;

    for (k=_d0Scale.begin(); k!=_d0Scale.end(); k++) _parFile = _parFile && k->second;
    for (k=_z0Scale.begin(); k!=_z0Scale.end(); k++) _parFile = _parFile && k->second;
    for (k=_thtScale.begin(); k!=_thtScale.end(); k++) _parFile = _parFile && k->second;
    for (k=_phiScale.begin(); k!=_phiScale.end(); k++) _parFile = _parFile && k->second;
    for (k=_momScale.begin(); k!=_momScale.end(); k++) _parFile = _parFile && k->second;

    if (!_parFile)
      cout <<" -W-  (PndFsmCmpDet::readParameters) - file "<<_parFileName
           <<" is defective. Using constant vertex reso"<<endl;
  }
}

PndFsmResponse * PndFsmCmpDet::respond ( PndFsmTrack *  t )

virtual

Implements PndFsmAbsDet.

Definition at line 110 of file PndFsmCmpDet.cxx.

References _d0, _d0ResMulti, _d0Scale, _mom, _mom0, _momResMulti, _momScale, _parFile, _phi, _phiResMulti, _phiScale, _pidLhMulti, _tht, _tht0, _tht1, _thtResMulti, _thtScale, _z0, _z0ResMulti, _z0Scale, PndFsmTrack::charge(), dE, PndFsmResponse::dE(), PndFsmResponse::detected(), PndFsmResponse::dm(), PndFsmResponse::dp(), PndFsmResponse::dphi(), PndFsmResponse::DrcBarrelThtc(), PndFsmResponse::DrcBarrelThtcErr(), PndFsmResponse::DrcDiscThtc(), PndFsmResponse::DrcDiscThtcErr(), PndFsmResponse::dt(), dtheta, PndFsmResponse::dtheta(), PndFsmResponse::dV(), eval(), fabs(), fDetList, fSubtractDetList, PndFsmResponse::LHElectron(), PndFsmResponse::LHKaon(), PndFsmResponse::LHMuon(), PndFsmResponse::LHPion(), PndFsmResponse::LHProton(), PndFsmResponse::m2(), m2(), PndFsmResponse::m2Err(), PndFsmResponse::MvddEdx(), PndFsmResponse::MvddEdxErr(), p, PndFsmTrack::p4(), PndFsmTrack::pdt(), pid(), PndFsmResponse::RichThtc(), PndFsmResponse::RichThtcErr(), PndFsmResponse::setdE(), PndFsmResponse::setDetected(), PndFsmResponse::setDetector(), PndFsmResponse::setdm(), PndFsmResponse::setdp(), PndFsmResponse::setdphi(), PndFsmResponse::setDrcBarrelThtc(), PndFsmResponse::setDrcDiscThtc(), PndFsmResponse::setdt(), PndFsmResponse::setdtheta(), PndFsmResponse::setdV(), PndFsmResponse::setLHElectron(), PndFsmResponse::setLHKaon(), PndFsmResponse::setLHMuon(), PndFsmResponse::setLHPion(), PndFsmResponse::setLHProton(), PndFsmResponse::setm2(), PndFsmResponse::setMvddEdx(), PndFsmResponse::setRichThtc(), PndFsmResponse::setSttdEdx(), PndFsmResponse::setTpcdEdx(), sqrt(), PndFsmResponse::SttdEdx(), PndFsmResponse::SttdEdxErr(), theta, PndFsmResponse::TpcdEdx(), PndFsmResponse::TpcdEdxErr(), and val.

                                    {
  PndFsmResponse *result=new PndFsmResponse();

  result->setDetector(this);

  bool detected=false;

  double dE=0.0;
  double dp=0.0;
  double dtheta=0.0;
  double dphi=0.0;
  double dt=0.0;
  double dm=0.0;

  double m2=0;
  double MvddEdx=0;
  double TpcdEdx=0;
  double SttdEdx=0;
  double DrcDiscThtc=0;
  double DrcBarrelThtc=0;
  double RichThtc=0;

  double m2Err=0;
  double MvddEdxErr=0;
  double TpcdEdxErr=0;
  double SttdEdxErr=0;
  double DrcDiscThtcErr=0;
  double DrcBarrelThtcErr=0;
  double RichThtcErr=0;

  double dVx=0;
  double dVy=0;
  double dVz=0;

  double LH_e=1.0;
  double LH_mu=1.0;
  double LH_pi=1.0;
  double LH_K=1.0;
  double LH_p=1.0;

  double val=0.0;

  for (FsmAbsDetList::iterator iter=fDetList.begin();iter!=fDetList.end(); iter++) {
    PndFsmResponse* resp=(*iter)->respond(t);

    detected |= resp->detected();

    if (resp->detected()) {
      if (fabs(val = resp->dE()) > 1e-8)    dE+=1/(val*val);
      if (fabs(val = resp->dt()) > 1e-8)     dt += val*val;
      if (fabs(val = resp->dm()) > 1e-8)     dm +=val;
      if (fabs (val = resp->m2()) > 1e-11)    m2+=val;
      if (fabs (val = resp->MvddEdx()) > 1e-11)    MvddEdx+=val;
      if (fabs (val = resp->TpcdEdx()) > 1e-11)    TpcdEdx+=val;
      if (fabs (val = resp->SttdEdx()) > 1e-11)    SttdEdx+=val;
      if (fabs (val = resp->DrcDiscThtc()) > 1e-11)   DrcDiscThtc+=val;
      if (fabs (val = resp->DrcBarrelThtc()) > 1e-11)    DrcBarrelThtc+=val;
      if (fabs (val = resp->RichThtc()) > 1e-11)    RichThtc+=val;

      if (fabs (val = resp->m2Err()) > 1e-11)    m2Err+=val;
      if (fabs (val = resp->MvddEdxErr()) > 1e-11)    MvddEdxErr+=val;
      if (fabs (val = resp->TpcdEdxErr()) > 1e-11)    TpcdEdxErr+=val;
      if (fabs (val = resp->SttdEdxErr()) > 1e-11)    SttdEdxErr+=val;
      if (fabs (val = resp->DrcDiscThtcErr()) > 1e-11)   DrcDiscThtcErr+=val;
      if (fabs (val = resp->DrcBarrelThtcErr()) > 1e-11)    DrcBarrelThtcErr+=val;
      if (fabs (val = resp->RichThtcErr()) > 1e-11)    RichThtcErr+=val;

      double rawLHe  = resp->LHElectron();
      double rawLHmu = resp->LHMuon();
      double rawLHpi = resp->LHPion();
      double rawLHK  = resp->LHKaon();
      double rawLHp  = resp->LHProton();

      double sumRaw = rawLHe+rawLHmu+rawLHpi+rawLHK+rawLHp;

      if (sumRaw>0) {
        rawLHe  /= sumRaw;
        rawLHmu /= sumRaw;
        rawLHpi /= sumRaw;
        rawLHK  /= sumRaw;
        rawLHp  /= sumRaw;
        LH_e  *= rawLHe;
        LH_mu *= rawLHmu;
        LH_pi *= rawLHpi;
        LH_K  *= rawLHK;
        LH_p  *= rawLHp;
      } else {
        LH_e  *= 0.2;
        LH_mu *= 0.2;
        LH_pi *= 0.2;
        LH_K  *= 0.2;
        LH_p  *= 0.2;
      }
    }
  }

  // invoke parameterised vertex/momentum resolution
  // (this will overwrite dp, dtheta, dphi and dV)
    if ( _parFile ) {
        // haven't done neutral particles yet
        if ( fabs(t->charge())>1e-8 ) {
            int pid=abs(t->pdt());
            if ( pid==11 || pid==13 || pid==211 || pid==321 || pid==2212) {
                double p=t->p4().Vect().Mag();
                double theta=t->p4().Vect().Theta()*180/M_PI;
                // cut off slow and out of theta range particles
                // (this avoids floating point exceptions)
                detected &= ( p > _mom0[pid]->GetVal() );
                detected &= ( theta >= _tht0->GetVal() );
                detected &= ( theta <= _tht1->GetVal() );

                if (detected) {
                    dtheta = eval(_tht[pid], theta) * _thtScale[pid]->Eval(p) * _thtResMulti * 3.1416/180;
                    dphi = eval(_phi[pid], theta) * _phiScale[pid]->Eval(p) * _phiResMulti * 3.1416/180;
                    dp = eval(_mom[pid], theta) * _momScale[pid]->Eval(p) * _momResMulti * p;

                    dVx = eval(_d0[pid], theta) * _d0Scale[pid]->Eval(p) * _d0ResMulti;
                    dVy = eval(_d0[pid], theta) * _d0Scale[pid]->Eval(p) * _d0ResMulti;
                    dVz = eval(_z0[pid], theta) * _z0Scale[pid]->Eval(p) * _z0ResMulti;

                    dtheta=1/(dtheta*dtheta);
                    dphi=1/(dphi*dphi);
                    dp=1/(dp*dp);
                }
            }
        }
    }

  for (FsmAbsDetList::iterator iter=fSubtractDetList.begin();iter!=fSubtractDetList.end(); iter++) {
    PndFsmResponse* resp=(*iter)->respond(t);

    if (resp->detected()) {
      if (fabs(val = resp->dE()) > 1e-8)    dE-=1/(val*val);
      if (fabs(val = resp->dp()) > 1e-8)     dp-=1/(val*val);
      if (fabs(val = resp->dtheta())> 1e-8)    dtheta-=1/(val*val);
      if (fabs(val = resp->dphi()) > 1e-8)  dphi-=1/(val*val);
      if (fabs(val = resp->dt()) > 1e-8)     dt -= val*val;
      if (fabs(val = resp->dm()) > 1e-8)     dm -=val;
      if (fabs (val = resp->m2()) > 1e-11)    m2-=val;
      if (fabs (val = resp->MvddEdx()) > 1e-11)    MvddEdx-=val;
      if (fabs (val = resp->TpcdEdx()) > 1e-11)    TpcdEdx-=val;
      if (fabs (val = resp->SttdEdx()) > 1e-11)    SttdEdx-=val;
      if (fabs (val = resp->DrcDiscThtc()) > 1e-11)   DrcDiscThtc-=val;
      if (fabs (val = resp->DrcBarrelThtc()) > 1e-11)    DrcBarrelThtc-=val;
      if (fabs (val = resp->RichThtc()) > 1e-11)    RichThtc-=val;

      if (fabs (val = resp->m2Err()) > 1e-11)    m2Err-=val;
      if (fabs (val = resp->MvddEdxErr()) > 1e-11)    MvddEdxErr-=val;
      if (fabs (val = resp->TpcdEdxErr()) > 1e-11)    TpcdEdxErr-=val;
      if (fabs (val = resp->SttdEdxErr()) > 1e-11)    SttdEdxErr-=val;
      if (fabs (val = resp->DrcDiscThtcErr()) > 1e-11)   DrcDiscThtcErr-=val;
      if (fabs (val = resp->DrcBarrelThtcErr()) > 1e-11)    DrcBarrelThtcErr-=val;
      if (fabs (val = resp->RichThtcErr()) > 1e-11)    RichThtcErr-=val;

      if (fabs (val = resp->dV().X()) > 1e-11) dVx-=1/(val*val);
      if (fabs (val = resp->dV().Y()) > 1e-11) dVy-=1/(val*val);
      if (fabs (val = resp->dV().Z()) > 1e-11) dVz-=1/(val*val);

      double rawLHe  = resp->LHElectron();
      double rawLHmu = resp->LHMuon();
      double rawLHpi = resp->LHPion();
      double rawLHK  = resp->LHKaon();
      double rawLHp  = resp->LHProton();

      double sumRaw = rawLHe+rawLHmu+rawLHpi+rawLHK+rawLHp;

      if (sumRaw>0) {
        rawLHe  /= sumRaw;
        rawLHmu /= sumRaw;
        rawLHpi /= sumRaw;
        rawLHK  /= sumRaw;
        rawLHp  /= sumRaw;
        LH_e  /= rawLHe;
        LH_mu /= rawLHmu;
        LH_pi /= rawLHpi;
        LH_K  /= rawLHK;
        LH_p  /= rawLHp;
      }
    }
  }

  double sumLH = LH_e + LH_mu + LH_pi + LH_K + LH_p;

  if (sumLH>0) {
    LH_e  /= sumLH;
    LH_mu /= sumLH;
    LH_pi /= sumLH;
    LH_K  /= sumLH;
    LH_p  /= sumLH;
  } else {
    LH_e  = 0.2;
    LH_mu = 0.2;
    LH_pi = 0.2;
    LH_K  = 0.2;
    LH_p  = 0.2;
  }

  // this shifts all likelihoods linearly to the state
  // of no pid information at all i. e. all lhs are 0.2
  LH_e  = 0.2*(1-_pidLhMulti) + LH_e*_pidLhMulti;
  LH_mu = 0.2*(1-_pidLhMulti) + LH_mu*_pidLhMulti;
  LH_pi = 0.2*(1-_pidLhMulti) + LH_pi*_pidLhMulti;
  LH_K  = 0.2*(1-_pidLhMulti) + LH_K*_pidLhMulti;
  LH_p  = 0.2*(1-_pidLhMulti) + LH_p*_pidLhMulti;

  result->setdE( dE>0. ? 1/sqrt(dE) : 0.0 );
  result->setdp( dp>0. ? 1/sqrt(dp) : 0.0 );
  result->setdtheta( dtheta>0. ? 1/sqrt(dtheta) : 0.0 );
  result->setdphi( dphi>0. ? 1/sqrt(dphi) : 0.0 );
  result->setdt( sqrt(dt) );
  result->setdm(dm);

  result->setm2(m2, m2Err);
  result->setMvddEdx(MvddEdx,MvddEdxErr);
  result->setTpcdEdx(TpcdEdx,TpcdEdxErr);
  result->setSttdEdx(SttdEdx,SttdEdxErr);

  result->setDrcDiscThtc(DrcDiscThtc,DrcDiscThtcErr);
  result->setDrcBarrelThtc(DrcBarrelThtc,DrcBarrelThtcErr);
  result->setRichThtc(RichThtc,RichThtcErr);

  if (dVx > 0.) dVx=1./sqrt(dVx); else dVx = 0.0;
  if (dVy > 0.) dVy=1./sqrt(dVy); else dVy = 0.0;
  if (dVz > 0.) dVz=1./sqrt(dVz); else dVz = 0.0;

  result->setdV( dVx , dVy , dVz );

  result->setLHElectron(LH_e);
  result->setLHMuon(LH_mu);
  result->setLHPion(LH_pi);
  result->setLHKaon(LH_K);
  result->setLHProton(LH_p);

  result->setDetected(detected);

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

privatevirtual

Reimplemented from PndFsmAbsDet.

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

privatevirtual

Reimplemented from PndFsmAbsDet.

Definition at line 378 of file PndFsmCmpDet.cxx.

References _d0ResMulti, _momResMulti, _phiResMulti, _pidLhMulti, _tht0, _tht1, _thtResMulti, and _z0ResMulti.

                                                             {
  // *****************
  // include here all float parameters which should be settable
  // *****************

  bool knownName=true;

  if (name == "d0ResMulti")
    _d0ResMulti=value;
  else
  if (name == "z0ResMulti")
    _z0ResMulti=value;
  else
  if (name == "thtResMulti")
    _thtResMulti=value;
  else
  if (name == "phiResMulti")
    _phiResMulti=value;
  else
  if (name == "momResMulti")
    _momResMulti=value;
  else
  if (name == "thtMin")
    _tht0 = new TParameter<double>("tht0", value);
  else
  if (name == "thtMax")
    _tht1 = new TParameter<double>("tht1", value);
  else
  if (name == "pidLhMulti")
    _pidLhMulti=value;
  else
    knownName=false;

  return knownName;
}

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;};

bool PndFsmCmpDet::SubtractDetector	(	std::string	name,
		std::string	params = ""
	)

Definition at line 356 of file PndFsmCmpDet.cxx.

References PndFsmDetFactory::create(), fDetFact, and fSubtractDetList.

                                                              {
  PndFsmAbsDet *det=fDetFact.create(name,params);
  if (det)
    fSubtractDetList.push_back(det);
  return det;
}

Member Data Documentation

std::map<int, TSpline3*> PndFsmCmpDet::_d0

private

Definition at line 86 of file PndFsmCmpDet.h.

Referenced by readParameters(), and respond().

double PndFsmCmpDet::_d0ResMulti

private

Definition at line 76 of file PndFsmCmpDet.h.

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

std::map<int, TF1*> PndFsmCmpDet::_d0Scale

private

Definition at line 94 of file PndFsmCmpDet.h.

Referenced by readParameters(), and respond().

std::string PndFsmAbsDet::_detName

protectedinherited

Definition at line 92 of file PndFsmAbsDet.h.

Referenced by PndFsmAbsDet::detName(), 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(), 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(), PndFsmEmcBwCap::print(), PndFsmEmcFS::print(), PndFsmEmcFwCap::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(), PndFsmSttPid::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(), PndFsmMvd2::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(), PndFsmEmcFS::print(), PndFsmSttPid::print(), PndFsmMdcTS::print(), PndFsmEmcBwCap::print(), PndFsmMvd2::print(), PndFsmMdcFS::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(), PndFsmMvd2::setParameter(), PndFsmMdcTS::setParameter(), PndFsmSttPid::setParameter(), PndFsmEmcFwCap::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(), PndFsmSttPid::respond(), PndFsmRich::respond(), PndFsmDrcDisc::respond(), PndFsmDrcBarrel::respond(), and PndFsmEffTracker::respond().

std::map<int, TSpline3*> PndFsmCmpDet::_mom

private

Definition at line 90 of file PndFsmCmpDet.h.

Referenced by readParameters(), and respond().

std::map<int, TParameter<double>*> PndFsmCmpDet::_mom0

private

Definition at line 91 of file PndFsmCmpDet.h.

Referenced by initParameters(), readParameters(), and respond().

double PndFsmCmpDet::_momResMulti

private

Definition at line 80 of file PndFsmCmpDet.h.

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

std::map<int, TF1*> PndFsmCmpDet::_momScale

private

Definition at line 98 of file PndFsmCmpDet.h.

Referenced by readParameters(), and respond().

bool PndFsmCmpDet::_parFile

private

Definition at line 75 of file PndFsmCmpDet.h.

Referenced by initParameters(), readParameters(), and respond().

TString PndFsmCmpDet::_parFileName

private

Definition at line 74 of file PndFsmCmpDet.h.

Referenced by initParameters(), PndFsmCmpDet(), and readParameters().

std::map<int, TSpline3*> PndFsmCmpDet::_phi

private

Definition at line 89 of file PndFsmCmpDet.h.

Referenced by readParameters(), and respond().

double PndFsmCmpDet::_phiResMulti

private

Definition at line 79 of file PndFsmCmpDet.h.

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

std::map<int, TF1*> PndFsmCmpDet::_phiScale

private

Definition at line 97 of file PndFsmCmpDet.h.

Referenced by readParameters(), and respond().

double PndFsmCmpDet::_pidLhMulti

private

Definition at line 84 of file PndFsmCmpDet.h.

Referenced by initParameters(), respond(), 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(), PndFsmSttPid::detected(), PndFsmMvdPid::detected(), PndFsmTof::detected(), PndFsmMdtPid::detected(), PndFsmEffTracker::detected(), PndFsmDrcBarrel::detected(), PndFsmEmcPid::detected(), PndFsmAbsDet::PndFsmAbsDet(), PndFsmRich::respond(), PndFsmTof::respond(), PndFsmMvd::respond(), PndFsmSttPid::respond(), PndFsmMvd2::respond(), PndFsmStt::respond(), PndFsmMvdPid::respond(), PndFsmDrcDisc::respond(), PndFsmDrcBarrel::respond(), and PndFsmEffTracker::respond().

std::map<int, TSpline3*> PndFsmCmpDet::_tht

private

Definition at line 88 of file PndFsmCmpDet.h.

Referenced by readParameters(), and respond().

TParameter<double>* PndFsmCmpDet::_tht0

private

Definition at line 92 of file PndFsmCmpDet.h.

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

TParameter<double>* PndFsmCmpDet::_tht1

private

Definition at line 93 of file PndFsmCmpDet.h.

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

double PndFsmCmpDet::_thtResMulti

private

Definition at line 78 of file PndFsmCmpDet.h.

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

std::map<int, TF1*> PndFsmCmpDet::_thtScale

private

Definition at line 96 of file PndFsmCmpDet.h.

Referenced by readParameters(), and respond().

std::map<int, TSpline3*> PndFsmCmpDet::_z0

private

Definition at line 87 of file PndFsmCmpDet.h.

Referenced by readParameters(), and respond().

double PndFsmCmpDet::_z0ResMulti

private

Definition at line 77 of file PndFsmCmpDet.h.

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

std::map<int, TF1*> PndFsmCmpDet::_z0Scale

private

Definition at line 95 of file PndFsmCmpDet.h.

Referenced by readParameters(), and respond().

PndFsmDetFactory PndFsmCmpDet::fDetFact

private

Definition at line 72 of file PndFsmCmpDet.h.

Referenced by AddDetector(), and SubtractDetector().

FsmAbsDetList PndFsmCmpDet::fDetList

private

Definition at line 70 of file PndFsmCmpDet.h.

Referenced by AddDetector(), and respond().

FsmAbsDetList PndFsmCmpDet::fSubtractDetList

private

Definition at line 71 of file PndFsmCmpDet.h.

Referenced by respond(), and SubtractDetector().

---

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

• PndFsmCmpDet.h
• PndFsmCmpDet.cxx