Fork me on GitHub

source: git/classes/DelphesFormula.cc@ 3db5282

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 3db5282 was 1fa50c2, checked in by Pavel Demin <pavel.demin@…>, 10 years ago

fix GPLv3 header

  • Property mode set to 100644
File size: 2.6 KB
Line 
1/*
2 * Delphes: a framework for fast simulation of a generic collider experiment
3 * Copyright (C) 2012-2014 Universite catholique de Louvain (UCL), Belgium
4 *
5 * This program is free software: you can redistribute it and/or modify
6 * it under the terms of the GNU General Public License as published by
7 * the Free Software Foundation, either version 3 of the License, or
8 * (at your option) any later version.
9 *
10 * This program is distributed in the hope that it will be useful,
11 * but WITHOUT ANY WARRANTY; without even the implied warranty of
12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13 * GNU General Public License for more details.
14 *
15 * You should have received a copy of the GNU General Public License
16 * along with this program. If not, see <http://www.gnu.org/licenses/>.
17 */
18
19
20#include "classes/DelphesFormula.h"
21
22#include "TString.h"
23
24#include <stdexcept>
25#include <string>
26
27using namespace std;
28
29//------------------------------------------------------------------------------
30
31DelphesFormula::DelphesFormula() :
32 TFormula()
33{
34}
35
36//------------------------------------------------------------------------------
37
38DelphesFormula::DelphesFormula(const char *name, const char *expression) :
39 TFormula()
40{
41}
42
43//------------------------------------------------------------------------------
44
45DelphesFormula::~DelphesFormula()
46{
47}
48
49//------------------------------------------------------------------------------
50
51Int_t DelphesFormula::Compile(const char *expression)
52{
53 string buffer;
54 const char *it;
55 for(it = expression; *it; ++it)
56 {
57 if(*it == ' ' || *it == '\t' || *it == '\r' || *it == '\n' || *it == '\\' ) continue;
58 buffer.push_back(*it);
59 }
60 if(TFormula::Compile(buffer.c_str()) != 0)
61 {
62 throw runtime_error("Invalid formula.");
63 }
64 return 0;
65}
66
67//------------------------------------------------------------------------------
68
69Double_t DelphesFormula::Eval(Double_t pt, Double_t eta, Double_t phi, Double_t energy)
70{
71 Double_t x[4] = {pt, eta, phi, energy};
72 return EvalPar(x);
73}
74
75//------------------------------------------------------------------------------
76
77Int_t DelphesFormula::DefinedVariable(TString &chaine, Int_t &action)
78{
79 action = kVariable;
80 if(chaine == "pt")
81 {
82 if(fNdim < 1) fNdim = 1;
83 return 0;
84 }
85 else if(chaine == "eta")
86 {
87 if(fNdim < 2) fNdim = 2;
88 return 1;
89 }
90 else if(chaine == "phi")
91 {
92 if(fNdim < 3) fNdim = 3;
93 return 2;
94 }
95 else if(chaine == "energy")
96 {
97 if(fNdim < 4) fNdim = 4;
98 return 3;
99 }
100 return -1;
101}
102
103//------------------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.