#include <PndFsmEmcPid.h>

Inheritance diagram for PndFsmEmcPid:

Public Member Functions
	PndFsmEmcPid ()

	PndFsmEmcPid (ArgList &par)

virtual	~PndFsmEmcPid ()

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
TH2F *	_emcPidPdf [6]

double	_pmin

double	_ptmin

double	_thtMin

double	_thtMax

double	_phiMin

double	_phiMax

std::string	_parFileName

Detailed Description

Definition at line 42 of file PndFsmEmcPid.h.

Constructor & Destructor Documentation

PndFsmEmcPid::PndFsmEmcPid ( )

Default ctor.

Definition at line 58 of file PndFsmEmcPid.cxx.

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

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

PndFsmEmcPid::PndFsmEmcPid ( ArgList & par )

Definition at line 71 of file PndFsmEmcPid.cxx.

References _phiMax, _phiMin, _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;
   _phiMin=_phiMin*M_PI/180.0;
   _phiMax=_phiMax*M_PI/180.0;
   readParameters();
 
   //print(std::cout);
 }

PndFsmEmcPid::~PndFsmEmcPid ( )

virtual

Destructor

Definition at line 95 of file PndFsmEmcPid.cxx.

96 {

97 }

Member Function Documentation

bool PndFsmEmcPid::detected ( PndFsmTrack * t ) const

private

Definition at line 175 of file PndFsmEmcPid.cxx.

References PndFsmAbsDet::_efficiency, _phiMin, _pmin, _ptmin, PndFsmAbsDet::_rand, _thtMin, PndFsmTrack::charge(), fabs(), p, PndFsmTrack::p4(), phi, pt(), and theta.

Referenced by respond().

 {
     double theta = t->p4().Theta();
     double phi = t->p4().Phi();
     double p=t->p4().Vect().Mag();
     double pt=t->p4().Pt();
     double charge=t->charge();
     //int type=abs(t->pdt()); //[R.K. 01/2017] unused variable
 
         return ( fabs(charge)>1e-6 && theta>=_thtMin && theta<=_thtMax && phi>=_phiMin && phi<=_phiMax && p>_pmin && pt> _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(), 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().

74 {return _detName;}

PndFsmAbsDet::_detName

std::string _detName

Definition: PndFsmAbsDet.h:92

Bool_t PndFsmAbsDet::doesPid ( ) const

inlineinherited

Definition at line 76 of file PndFsmAbsDet.h.

References PndFsmAbsDet::_doesPid.

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

76 {return _doesPid;};

PndFsmAbsDet::_doesPid

Bool_t _doesPid

Definition: PndFsmAbsDet.h:96

void PndFsmEmcPid::initParameters ( )

privatevirtual

Reimplemented from PndFsmAbsDet.

Definition at line 203 of file PndFsmEmcPid.cxx.

References PndFsmAbsDet::_detName, PndFsmAbsDet::_efficiency, _parFileName, _phiMax, _phiMin, _pmin, _ptmin, _thtMax, and _thtMin.

Referenced by PndFsmEmcPid().

 {
   _detName = "EmcPid";
   _pmin   = 0.00;
   _ptmin  = 0.00;
   _thtMin = 22.0;
   _thtMax = 140.0;
   _phiMin=-360.;
   _phiMax=360.;
   _efficiency=1.0;
   _parFileName     = "$VMCWORKDIR/fastsim/EmcPidPdf.root";
 }

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

privatevirtual

Reimplemented from PndFsmAbsDet.

Definition at line 189 of file PndFsmEmcPid.cxx.

References PndFsmAbsDet::_detName, PndFsmAbsDet::_efficiency, _parFileName, _phiMax, _phiMin, _pmin, _ptmin, _thtMax, and _thtMin.

 {
   o <<"Detector <"<<_detName<<">"<<endl;
   o  <<"  _pmin = "<<_pmin<<endl;
   o  <<"  _ptmin = "<<_ptmin<<endl;
   o  <<"  _thtMin = "<<_thtMin<<endl;
   o  <<"  _thtMax = "<<_thtMax<<endl;
   o  <<"  _phiMin = "<<_phiMin<<endl;
   o  <<"  _phiMax = "<<_phiMax<<endl;
   o  <<"  _efficiency = "<<_efficiency<<endl;
   o  <<"  _parFileName     = "<<_parFileName<<endl;
 }

bool PndFsmEmcPid::readParameters ( )

private

Definition at line 251 of file PndFsmEmcPid.cxx.

References _emcPidPdf, _parFileName, exit(), f, and i.

Referenced by PndFsmEmcPid().

 {
 
   TFile *f=new TFile(_parFileName.c_str());
 
   for (int i=0;i<6;i++)
   {
     _emcPidPdf[i]=0;
   }
 
   if (f->IsZombie())
   {
     cout <<" -W-  (PndFsmEmcPid::readParameters) - file "<<_parFileName.c_str()
          <<" doesn't exist."<<endl;
         exit(0);
   }
   else
   {
     _emcPidPdf[0]=(TH2F*)f->Get("hpdf_e");
     _emcPidPdf[1]=(TH2F*)f->Get("hpdf_mu");
     _emcPidPdf[2]=(TH2F*)f->Get("hpdf_pi");
     _emcPidPdf[3]=(TH2F*)f->Get("hpdf_k");
     _emcPidPdf[4]=(TH2F*)f->Get("hpdf_p");
     _emcPidPdf[5]=(TH2F*)f->Get("hpdf_pb");
 
     for (int i=0;i<6;i++) _emcPidPdf[i]->SetDirectory(0);
 
     f->Close();
   }
   delete f;
 
   return true;
 }

PndFsmResponse * PndFsmEmcPid::respond ( PndFsmTrack * t )

virtual

Accessors to contained information

Implements PndFsmAbsDet.

Definition at line 104 of file PndFsmEmcPid.cxx.

References _emcPidPdf, PndFsmTrack::charge(), detected(), idx, P, p, PndFsmTrack::p4(), PndFsmTrack::pdt(), PndFsmResponse::setDetected(), PndFsmResponse::setDetector(), PndFsmResponse::setEmcEcal(), PndFsmResponse::setLHElectron(), PndFsmResponse::setLHKaon(), PndFsmResponse::setLHMuon(), PndFsmResponse::setLHPion(), and PndFsmResponse::setLHProton().

 {
   PndFsmResponse *result=new PndFsmResponse();
 
   result->setDetector(this);
   bool wasDetected=detected(t);
   result->setDetected(wasDetected);
 
   if (wasDetected)
   {
     int    type   = t->pdt();
     double p      = t->p4().Vect().Mag();
     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 = 10.0
         double hpmax = _emcPidPdf[0]->GetXaxis()->GetXmax();
         if (p>hpmax) p=hpmax;
 
         // get the bin corresponding to p
         int currbin = _emcPidPdf[0]->GetXaxis()->FindBin(p);
 
         // get the slice containing the pdf emcecal(p) for true particle type
         TH1D *hpdf = _emcPidPdf[idx]->ProjectionY("_tmppdf",currbin, currbin);
 
         // 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->setEmcEcal(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.;
         // find P up e ... pbar and sum
         for (int k=0;k<6;++k)
         {
                 P[k] = _emcPidPdf[k]->GetBinContent(currbin,xsigbin);
                 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().

75 {_detName = name;};

name

TString name

Definition: createRootGeoFile.C:82

PndFsmAbsDet::_detName

std::string _detName

Definition: PndFsmAbsDet.h:92

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

privatevirtual

Reimplemented from PndFsmAbsDet.

Definition at line 217 of file PndFsmEmcPid.cxx.

References PndFsmAbsDet::_efficiency, _phiMax, _phiMin, _pmin, _ptmin, _thtMax, and _thtMin.

 {
   // *****************
   // include here all parameters which should be settable via tcl
   // *****************
 
   bool knownName=true;
 
   if (name == "pmin")
     _pmin=value;
   else
   if (name == "ptmin")
     _ptmin=value;
   else
   if (name == "thtMin")
     _thtMin=value;
   else
   if (name == "thtMax")
     _thtMax=value;
   else
   if (name == "phiMin")
     _phiMin=value;
   else
   if (name == "phiMax")
     _phiMax=value;
   else
   if (name == "efficiency")
     _efficiency=value;
   else
     knownName=false;
 
   return knownName;
 }