EvtMultiChannelParser Class Reference

#include <EvtMultiChannelParser.hh>

List of all members.

Public Member Functions
	EvtMultiChannelParser ()
	~EvtMultiChannelParser ()
void	parse (const char *file, const char *model)
void	parse (const std::vector< std::string > &v)
double	pdfMax () const
int	nScan () const
double	dm () const
std::vector< std::string >	amp (int i) const
std::vector< std::string >	ampConj (int i) const
EvtComplex	ampCoef (int i) const
EvtComplex	ampConjCoef (int i) const
int	coefFormat (int i) const
int	coefConjFormat (int i) const
int	getNAmp () const
int	getNAmpConj () const
Static Public Member Functions
static EvtDecayMode	getDecayMode (const char *file)
static void	parseComplexCoef (int &i, const std::vector< std::string > &v, EvtComplex &c, int &format)
static double	parseRealCoef (int &i, const std::vector< std::string > &v)
static bool	isKeyword (const std::string &s)
Private Attributes
double	_pdfMax
int	_nScan
double	_dm
std::vector< std::vector< std::string > >	_amp
std::vector< std::vector< std::string > >	_ampConj
std::vector< EvtComplex >	_ampCoef
std::vector< int >	_coefFormat
std::vector< EvtComplex >	_ampConjCoef
std::vector< int >	_coefConjFormat

---

Detailed Description

Definition at line 31 of file EvtMultiChannelParser.hh.

---

Constructor & Destructor Documentation

EvtMultiChannelParser::EvtMultiChannelParser

(

)

[inline]

Definition at line 35 of file EvtMultiChannelParser.hh.

: _pdfMax(-1.), _nScan(0), _dm(0.) {}

EvtMultiChannelParser::~EvtMultiChannelParser

(

)

[inline]

Definition at line 36 of file EvtMultiChannelParser.hh.

{}

---

Member Function Documentation

std::vector<std::string> EvtMultiChannelParser::amp

(

int

i

)

const [inline]

Definition at line 52 of file EvtMultiChannelParser.hh.

References _amp.

Referenced by EvtAmpFactory< EvtDalitzPoint >::build().

{ return _amp[i]; }

EvtComplex EvtMultiChannelParser::ampCoef

(

int

i

)

const [inline]

Definition at line 54 of file EvtMultiChannelParser.hh.

References _ampCoef.

Referenced by EvtAmpFactory< EvtDalitzPoint >::build().

{ return _ampCoef[i]; }

std::vector<std::string> EvtMultiChannelParser::ampConj

(

int

i

)

const [inline]

Definition at line 53 of file EvtMultiChannelParser.hh.

References _ampConj.

Referenced by EvtAmpFactory< EvtDalitzPoint >::build().

{ return _ampConj[i]; }

EvtComplex EvtMultiChannelParser::ampConjCoef

(

int

i

)

const [inline]

Definition at line 55 of file EvtMultiChannelParser.hh.

References _ampConjCoef.

Referenced by EvtAmpFactory< EvtDalitzPoint >::build().

{ return _ampConjCoef[i]; }

int EvtMultiChannelParser::coefConjFormat

(

int

i

)

const [inline]

Definition at line 58 of file EvtMultiChannelParser.hh.

References _coefConjFormat.

{ return _coefConjFormat[i]; }

int EvtMultiChannelParser::coefFormat

(

int

i

)

const [inline]

Definition at line 57 of file EvtMultiChannelParser.hh.

References _coefFormat.

{ return _coefFormat[i]; }

double EvtMultiChannelParser::dm

(

)

const [inline]

Definition at line 50 of file EvtMultiChannelParser.hh.

References _dm.

{ return _dm; }

EvtDecayMode EvtMultiChannelParser::getDecayMode

(

const char *

file

)

[static]

Definition at line 32 of file EvtMultiChannelParser.cc.

References EvtPDL::getId(), EvtParser::getNToken(), EvtPDL::getStdHep(), EvtParser::getToken(), genRecEmupikp::i, EvtDecayMode::mode(), TrigConf::N, EvtParser::Read(), and deljobs::string.

{
  // Open file, read tokens

  EvtParser parser;
  parser.Read(file);

  // Seek Decay

  int i = 0;
  int N = parser.getNToken();
  while(i<N) {
    
    std::string tok = parser.getToken(i++);
    if(tok == std::string("Decay")) break;
  }
  
  // Get mother

  string mother = string(parser.getToken(i++).c_str());
  std::string bf = parser.getToken(i++);

  vector<string> dauV;
  // Get daughters
    
  while(1) {

    std::string d = parser.getToken(i++);

    if(EvtPDL::getStdHep(EvtPDL::getId(d.c_str())) == 0) break;
    
    dauV.push_back(string(d.c_str()));
  }
  
  EvtDecayMode mode(mother,dauV);
  printf("Decay File defines mode %s\n",mode.mode());

  return mode;
}

int EvtMultiChannelParser::getNAmp

(

)

const [inline]

Definition at line 60 of file EvtMultiChannelParser.hh.

References _amp.

Referenced by EvtAmpFactory< EvtDalitzPoint >::build().

{ return _amp.size(); } 

int EvtMultiChannelParser::getNAmpConj

(

)

const [inline]

Definition at line 61 of file EvtMultiChannelParser.hh.

References _ampConj.

Referenced by EvtAmpFactory< EvtDalitzPoint >::build().

{ return _ampConj.size(); } 

bool EvtMultiChannelParser::isKeyword

(

const std::string &

s

)

[static]

Definition at line 267 of file EvtMultiChannelParser.cc.

References deljobs::string.

Referenced by parse().

{
  if(s == std::string("AMPLITUDE")) return true;
  if(s == std::string("CONJUGATE")) return true;
  if(s == std::string("COEFFICIENT")) return true;
  return false;
}

int EvtMultiChannelParser::nScan

(

)

const [inline]

Definition at line 49 of file EvtMultiChannelParser.hh.

References _nScan.

Referenced by EvtIntervalDecayAmp< EvtDalitzPoint >::init().

{ return _nScan; }

void EvtMultiChannelParser::parse

(

const std::vector< std::string > &

v

)

Definition at line 121 of file EvtMultiChannelParser.cc.

References _amp, _ampCoef, _ampConj, _ampConjCoef, _coefConjFormat, _coefFormat, _dm, _nScan, _pdfMax, EvtCyclic3::c_str(), genRecEmupikp::i, isKeyword(), parseComplexCoef(), and deljobs::string.

{
  // place holder for strtod
  char** tc = 0;


  // Get PDF maximum or number of points to 
  // use in the scan.

  if(v[0] == std::string("MAXPDF")) {
    
    _pdfMax = strtod(v[1].c_str(),tc);
    if(_pdfMax <= 0) { printf("Bad pdfMax=%f\n",_pdfMax); assert(0); }
  }
  else 
    if(v[0] == std::string("SCANPDF")) {

      _nScan = atoi(v[1].c_str());
    }
    else {

      printf("Error parsing decay file\n");
      assert(0);
    }


  // Now parse the rest of file for amplitude specifications.

  bool conjugate = false;
  int i = 2;
  assert(isKeyword(v[2]));
  
  while(i < v.size()) {
  
    int i0 = i;
  
    // Switch to conjugate amplitudes after keyword
    if(v[i] == std::string("CONJUGATE")) {
      
      assert(conjugate == false);
      conjugate = true;
      assert(!isKeyword(v[++i]));
      _dm =  strtod(v[i].c_str(),tc);
      i++;
    }

    std::vector<std::string> params;
    EvtComplex c;
    int format;

    if(!conjugate && v[i] == std::string("AMPLITUDE")) {

      while(!isKeyword(v[++i])) params.push_back(v[i]);
      _amp.push_back(params);
      
      parseComplexCoef(i,v,c,format);
      _ampCoef.push_back(c);
      _coefFormat.push_back(format);
      continue;
    }
    else 
      if(conjugate && v[i] == std::string("AMPLITUDE")) {

        while(!isKeyword(v[++i])) params.push_back(v[i]);
        _ampConj.push_back(params);
        parseComplexCoef(i,v,c,format);
        _ampConjCoef.push_back(c);
        _coefConjFormat.push_back(format);
        continue;
      }
      else {
        
        printf("Expect keyword, found parameter %s\n",v[i].c_str());
        assert(0);
      }
    

    assert(i > i0);
  }

  printf("PARSING SUCCESSFUL\n");
  printf("%d amplitude terms\n",_amp.size());
  printf("%d conj amplitude terms\n",_ampConj.size());
}

void EvtMultiChannelParser::parse	(	const char *	file,
		const char *	model
	)

Definition at line 74 of file EvtMultiChannelParser.cc.

References EvtParser::getNToken(), EvtParser::getToken(), genRecEmupikp::i, TrigConf::N, EvtParser::Read(), deljobs::string, and v.

Referenced by EvtIntervalDecayAmp< EvtDalitzPoint >::init().

{
  // Open file, read tokens

  EvtParser parser;
  parser.Read(file);


  // Get parameters (tokens between the model name and ;)

  int i = 0;
  int N = parser.getNToken();

  // Seek the model name

  while(i<N) {
   
    std::string tok = parser.getToken(i++);
    if(tok == std::string(model)) break;
  }
  if(i == N) {

    printf("No model %s found in decay file %s",model,file);
    exit(0);
  }  

  
  // Add all tokens up to a semicolon to vector 

  std::vector<std::string> v;
  while(i<N) {

    std::string tok = parser.getToken(i++);
    if(tok == std::string(";")) break;
    else v.push_back(tok);
  }
  
  if(i == N) {

    printf("No terminating ; found in decay file %s",file);
    assert(0);
  }

  parse(v);
}

void EvtMultiChannelParser::parseComplexCoef	(	int &	i,
		const std::vector< std::string > &	v,
		EvtComplex &	c,
		int &	format
	)			[static]

Definition at line 208 of file EvtMultiChannelParser.cc.

References EvtCyclic3::c_str(), CARTESIAN, cos(), EvtConst::pi, POLAR_DEG, POLAR_RAD, sin(), and deljobs::string.

Referenced by parse().

{
  // place holder for strtod
  char** tc = 0;

  assert(v[i++] == std::string("COEFFICIENT"));

  if(v[i] == std::string("POLAR_DEG")) {

      double mag = strtod(v[i+1].c_str(),tc);
      double phaseRad = strtod(v[i+2].c_str(),tc)*EvtConst::pi/180.0;
      i += 3;
      c = EvtComplex(mag*cos(phaseRad),mag*sin(phaseRad));  
      format = POLAR_DEG;
  }
  else if(v[i] == std::string("POLAR_RAD")) {
    
    double mag = strtod(v[i+1].c_str(),tc);
    double phaseRad = strtod(v[i+2].c_str(),tc);
    i += 3;
    c = EvtComplex(mag*cos(phaseRad),mag*sin(phaseRad));  
    format = POLAR_RAD;    
  }
  else if(v[i] == std::string("CARTESIAN")) {

    double re = strtod(v[i+1].c_str(),tc);
    double im = strtod(v[i+2].c_str(),tc);
    i += 3;
    c = EvtComplex(re,im);  
    format = CARTESIAN;
  }
  else {
    
    printf("Invalid format %s for complex coefficient\n",v[i].c_str());
    exit(0);
  }
}

double EvtMultiChannelParser::parseRealCoef	(	int &	i,
		const std::vector< std::string > &	v
	)			[static]

Definition at line 248 of file EvtMultiChannelParser.cc.

References EvtCyclic3::c_str(), and deljobs::string.

{
  // place holder for strtod
  char** tc = 0;
  double value = 0;

  if(v[i] == std::string("COEFFICIENT")) {

    value = strtod(v[i+1].c_str(),tc);
  }
  else assert(0);

  i += 2;

  assert(value > 0.);
  return value;
}

double EvtMultiChannelParser::pdfMax

(

)

const [inline]

Definition at line 48 of file EvtMultiChannelParser.hh.

References _pdfMax.

Referenced by EvtIntervalDecayAmp< EvtDalitzPoint >::init().

{ return _pdfMax; }

---

Member Data Documentation

std::vector<std::vector<std::string> > EvtMultiChannelParser::_amp [private]

Definition at line 69 of file EvtMultiChannelParser.hh.

Referenced by amp(), getNAmp(), and parse().

std::vector<EvtComplex> EvtMultiChannelParser::_ampCoef [private]

Definition at line 71 of file EvtMultiChannelParser.hh.

Referenced by ampCoef(), and parse().

std::vector<std::vector<std::string> > EvtMultiChannelParser::_ampConj [private]

Definition at line 70 of file EvtMultiChannelParser.hh.

Referenced by ampConj(), getNAmpConj(), and parse().

std::vector<EvtComplex> EvtMultiChannelParser::_ampConjCoef [private]

Definition at line 73 of file EvtMultiChannelParser.hh.

Referenced by ampConjCoef(), and parse().

std::vector<int> EvtMultiChannelParser::_coefConjFormat [private]

Definition at line 74 of file EvtMultiChannelParser.hh.

Referenced by coefConjFormat(), and parse().

std::vector<int> EvtMultiChannelParser::_coefFormat [private]

Definition at line 72 of file EvtMultiChannelParser.hh.

Referenced by coefFormat(), and parse().

double EvtMultiChannelParser::_dm [private]

Definition at line 67 of file EvtMultiChannelParser.hh.

Referenced by dm(), and parse().

int EvtMultiChannelParser::_nScan [private]

Definition at line 66 of file EvtMultiChannelParser.hh.

Referenced by nScan(), and parse().

double EvtMultiChannelParser::_pdfMax [private]

Definition at line 65 of file EvtMultiChannelParser.hh.

Referenced by parse(), and pdfMax().

---