PndFsmMdtPid Class Reference

#include <PndFsmMdtPid.h>

Inheritance diagram for PndFsmMdtPid:

Public Member Functions
	PndFsmMdtPid ()
	PndFsmMdtPid (ArgList &par)
virtual	~PndFsmMdtPid ()
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
void	initParameters ()
void	print (std::ostream &o)
bool	setParameter (std::string &name, double value)
bool	readParameters ()

Private Attributes
TH3F *	_mdtPidPdf [6]
double	_thtMin
double	_thtMax
double	_pmin
double	_misId
bool	_useFlat
std::string	_parFileName

Detailed Description

Definition at line 42 of file PndFsmMdtPid.h.

Constructor & Destructor Documentation

PndFsmMdtPid::PndFsmMdtPid

(

)

Default ctor.

Definition at line 60 of file PndFsmMdtPid.cxx.

References _thtMax, _thtMin, initParameters(), and readParameters().

{
  initParameters();

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

PndFsmMdtPid::PndFsmMdtPid

(

ArgList &

par

)

Definition at line 70 of file PndFsmMdtPid.cxx.

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

{
  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;
  readParameters();
  //print(std::cout);
}

PndFsmMdtPid::~PndFsmMdtPid ( )

virtual

Destructor

Definition at line 88 of file PndFsmMdtPid.cxx.

{
}

Member Function Documentation

bool PndFsmMdtPid::detected ( PndFsmTrack *  t ) const

private

Definition at line 181 of file PndFsmMdtPid.cxx.

References PndFsmAbsDet::_efficiency, _misId, _pmin, PndFsmAbsDet::_rand, _thtMin, _useFlat, PndFsmTrack::charge(), p, PndFsmTrack::p4(), PndFsmTrack::pdt(), and theta.

Referenced by respond().

{
    double theta = t->p4().Theta();
    double p=t->p4().Vect().Mag();
    double charge=t->charge();
    int type=abs(t->pdt());

        if (_useFlat) return ( charge!=0.0 && theta>=_thtMin && theta<=_thtMax && p>_pmin && ((type==13 && _rand->Rndm()<=_efficiency) || (type==211 && _rand->Rndm()<=_misId)) );

        return ( charge!=0.0 && theta>=_thtMin && theta<=_thtMax && p>_pmin );
}

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

inlineinherited

Definition at line 74 of file PndFsmAbsDet.h.

References PndFsmAbsDet::_detName.

Referenced by PndFastSim::AddDetector(), PndFsmIdealPid::print(), PndFsmMvdPid::print(), PndFsmSimpleTracker::print(), PndFsmSimpleVtx::print(), PndFsmDetTemplate::print(), 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;};

void PndFsmMdtPid::initParameters ( )

privatevirtual

Reimplemented from PndFsmAbsDet.

Definition at line 205 of file PndFsmMdtPid.cxx.

References PndFsmAbsDet::_detName, PndFsmAbsDet::_efficiency, _misId, _parFileName, _pmin, _thtMax, _thtMin, and _useFlat.

Referenced by PndFsmMdtPid().

{
  _detName         = "MdtPid";
  _thtMin          = 5.;
  _thtMax          = 160.;
  _pmin            = 1.0;
  _misId                   = 0.01;
  _efficiency      = 1.0;
  _parFileName     = "$VMCWORKDIR/fastsim/MdtPidPdf3D.root";
  _useFlat                 = false;
}

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(), 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 PndFsmMdtPid::print ( std::ostream &  o )

privatevirtual

Reimplemented from PndFsmAbsDet.

Definition at line 194 of file PndFsmMdtPid.cxx.

References PndFsmAbsDet::_efficiency, _misId, _pmin, _thtMax, _thtMin, and PndFsmAbsDet::detName().

{
  o <<"Parameters for detector <"<<detName()<<">"<<endl;
  o  <<"  _thtMin = "<<_thtMin<<endl;
  o  <<"  _thtMax = "<<_thtMax<<endl;
  o  <<"  _pmin = "<<_pmin<<endl;
  o  <<"  _misId = "<<_misId<<endl;
  o  <<"  _efficiency = "<<_efficiency<<endl;
}

bool PndFsmMdtPid::readParameters ( )

private

Definition at line 246 of file PndFsmMdtPid.cxx.

References _mdtPidPdf, _parFileName, _useFlat, f, and i.

Referenced by PndFsmMdtPid().

{

  TFile *f=new TFile(_parFileName.c_str());

  for (int i=0;i<6;i++)
  {
    _mdtPidPdf[i]=0;
  }

  if (f->IsZombie())
  {
    cout <<" -W-  (PndFsmMdtPid::readParameters) - file "<<_parFileName.c_str()
         <<" doesn't exist. Using simple likelihood."<<endl;
        _useFlat = true;
  }
  else
  {
    _mdtPidPdf[0]=(TH3F*)f->Get("hpdf_e");
    _mdtPidPdf[1]=(TH3F*)f->Get("hpdf_mu");
    _mdtPidPdf[2]=(TH3F*)f->Get("hpdf_pi");
    _mdtPidPdf[3]=(TH3F*)f->Get("hpdf_k");
    _mdtPidPdf[4]=(TH3F*)f->Get("hpdf_p");
    _mdtPidPdf[5]=(TH3F*)f->Get("hpdf_pb");

        cout << _mdtPidPdf[0]<<endl;

    for (int i=0;i<6;i++) _mdtPidPdf[i]->SetDirectory(0);

    f->Close();
  }
  delete f;

  return true;
}

PndFsmResponse * PndFsmMdtPid::respond ( PndFsmTrack *  t )

virtual

Accessors to contained information

Implements PndFsmAbsDet.

Definition at line 97 of file PndFsmMdtPid.cxx.

References _mdtPidPdf, _misId, _useFlat, PndFsmTrack::charge(), detected(), idx, P, p, PndFsmTrack::p4(), PndFsmTrack::pdt(), PndFsmResponse::setDetected(), PndFsmResponse::setDetector(), PndFsmResponse::setLHElectron(), PndFsmResponse::setLHKaon(), PndFsmResponse::setLHMuon(), PndFsmResponse::setLHPion(), PndFsmResponse::setLHProton(), PndFsmResponse::setMuoIron(), and tht().

{
  PndFsmResponse *result=new PndFsmResponse();

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

  if (wasDetected)
  {
          if (_useFlat)  // just use very simple likelihoods for muons and pions
          {
                        result->setLHElectron(0.);
                        result->setLHMuon( 1-_misId );
                        result->setLHPion( _misId );
                        result->setLHKaon( 0. );
                        result->setLHProton( 0. );
          }
          else   // use sophisticated methode with histogram lookup likelihoods
          {
                        int    type   = t->pdt();
                        double p      = t->p4().Vect().Mag();
                        double tht    = t->p4().Theta()*57.296;
                        double charge = t->charge();

                        // convert pdg code in type: 0=e+-, 1=mu+-, 2=pi+-, 3=K+-, 4=p, 5=p-bar
                        int idx;
                        if (abs(type)==11) idx=0;
                        else if (abs(type)==13) idx=1;
                        else if (abs(type)==211) idx=2;
                        else if (abs(type)==321) idx=3;
                        else if (type==2212) idx=4;
                        else idx=5;

                        // histogram max p = 8.0
                        if (p>7.99) p=7.99;

                        // get the bin corresponding to p
                        int currbinT = _mdtPidPdf[0]->GetYaxis()->FindBin(tht);
                        int currbinP = _mdtPidPdf[0]->GetZaxis()->FindBin(p);

                        // get the slice containing the pdf emcecal(p) for true particle type
                        TH1D *hpdf = _mdtPidPdf[idx]->ProjectionX("_tmppdf",currbinT, currbinT, currbinP, currbinP);

                        // get a random emcecal value from the true distribution
                        // if distribution is empty, choose 0
                        double xsig = 0.;
                        if (hpdf->Integral()>0) xsig = hpdf->GetRandom();

                        // store the value in the response object
                        result->setMuoIron(xsig);

                        // find the bin corresponding to the value
                        int xsigbin = hpdf->FindBin(xsig);

                        // get the probability values for the different particle types for this emcecal value; sum needed for normalization
                        double P[6], Psum=0.;

                        // sum up e ... K
                        for (int k=0;k<6;++k)
                        {
                                P[k] = _mdtPidPdf[k]->GetBinContent(xsigbin,currbinT,currbinP);
                                if (k<4) Psum += P[k];
                        }
                        // add P_p or P_pbar depending on charge of particle
                        if (charge>0) Psum+=P[4];
                        else Psum+=P[5];

                        if (Psum<=0.) Psum=1.;

                        result->setLHElectron(P[0]/Psum);
                        result->setLHMuon(P[1]/Psum);
                        result->setLHPion(P[2]/Psum);
                        result->setLHKaon(P[3]/Psum);
                        if (charge>0) result->setLHProton(P[4]/Psum);
                        else result->setLHProton(P[5]/Psum);

          }
  }

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

privatevirtual

Reimplemented from PndFsmAbsDet.

Definition at line 218 of file PndFsmMdtPid.cxx.

References PndFsmAbsDet::_efficiency, _misId, _pmin, _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 == "pmin")
        _pmin=value;
  else
  if (name == "misId")
        _misId=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

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(), 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(), PndFsmEmcFwCap::print(), PndFsmEmcBwCap::print(), PndFsmEmcFS::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(), PndFsmMvd2::detected(), PndFsmEmcFwCap::detected(), PndFsmMdcTS::detected(), PndFsmMdcFS::detected(), PndFsmTof::detected(), PndFsmSttPid::detected(), PndFsmMvd::detected(), PndFsmEmcBwCap::detected(), PndFsmEmcFS::detected(), PndFsmStt::detected(), PndFsmSimpleVtx::detected(), PndFsmMvdPid::detected(), PndFsmRich::detected(), PndFsmEmcBarrel::detected(), PndFsmDrcDisc::detected(), detected(), PndFsmDrcBarrel::detected(), PndFsmEmcPid::detected(), PndFsmEffTracker::detected(), PndFsmSimpleTracker::initParameters(), PndFsmSimpleVtx::initParameters(), PndFsmDetTemplate::initParameters(), PndFsmMvdPid::initParameters(), PndFsmTof::initParameters(), initParameters(), PndFsmEmcPid::initParameters(), PndFsmMvd2::initParameters(), PndFsmEmcFwCap::initParameters(), PndFsmRich::initParameters(), PndFsmMdcFS::initParameters(), PndFsmSttPid::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(), print(), PndFsmEmcPid::print(), PndFsmSttPid::print(), PndFsmMdcFS::print(), PndFsmMdcTS::print(), PndFsmEmcBwCap::print(), PndFsmMvd2::print(), PndFsmEmcFS::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(), setParameter(), PndFsmTof::setParameter(), PndFsmEmcPid::setParameter(), PndFsmEmcFS::setParameter(), PndFsmEmcBwCap::setParameter(), PndFsmMdcTS::setParameter(), PndFsmMvd2::setParameter(), PndFsmSttPid::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(), PndFsmSttPid::respond(), PndFsmRich::respond(), PndFsmDrcDisc::respond(), PndFsmDrcBarrel::respond(), and PndFsmEffTracker::respond().

TH3F* PndFsmMdtPid::_mdtPidPdf[6]

private

Definition at line 75 of file PndFsmMdtPid.h.

Referenced by readParameters(), and respond().

double PndFsmMdtPid::_misId

private

Definition at line 80 of file PndFsmMdtPid.h.

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

std::string PndFsmMdtPid::_parFileName

private

Definition at line 84 of file PndFsmMdtPid.h.

Referenced by initParameters(), and readParameters().

double PndFsmMdtPid::_pmin

private

Definition at line 79 of file PndFsmMdtPid.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(), PndFsmStt::detected(), PndFsmSttPid::detected(), PndFsmEmcFwCap::detected(), PndFsmMvd::detected(), PndFsmSimpleVtx::detected(), PndFsmMvdPid::detected(), PndFsmMdcTS::detected(), PndFsmEmcFS::detected(), PndFsmMdcFS::detected(), PndFsmEmcBwCap::detected(), PndFsmMvd2::detected(), PndFsmTof::detected(), detected(), PndFsmEffTracker::detected(), PndFsmDrcBarrel::detected(), PndFsmEmcPid::detected(), PndFsmAbsDet::PndFsmAbsDet(), PndFsmRich::respond(), PndFsmMvd::respond(), PndFsmMvdPid::respond(), PndFsmSttPid::respond(), PndFsmMvd2::respond(), PndFsmStt::respond(), PndFsmTof::respond(), PndFsmDrcDisc::respond(), PndFsmDrcBarrel::respond(), and PndFsmEffTracker::respond().

double PndFsmMdtPid::_thtMax

private

Definition at line 78 of file PndFsmMdtPid.h.

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

double PndFsmMdtPid::_thtMin

private

Definition at line 77 of file PndFsmMdtPid.h.

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

bool PndFsmMdtPid::_useFlat

private

Definition at line 82 of file PndFsmMdtPid.h.

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

---

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

• PndFsmMdtPid.h
• PndFsmMdtPid.cxx