Fork me on GitHub

source: git/classes/DelphesCscClusterFormula.cc@ 29b722a

Last change on this file since 29b722a was 29b722a, checked in by christinaw97 <christina.wang@…>, 3 years ago

Added DelphesCscClusterFormula

  • Property mode set to 100644
File size: 2.4 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#include "classes/DelphesCscClusterFormula.h"
20#include "classes/DelphesClasses.h"
21
22#include "TString.h"
23
24#include <stdexcept>
25#include <iostream>
26
27using namespace std;
28
29//------------------------------------------------------------------------------
30
31DelphesCscClusterFormula::DelphesCscClusterFormula() :
32 TFormula()
33{
34}
35
36//------------------------------------------------------------------------------
37
38DelphesCscClusterFormula::DelphesCscClusterFormula(const char *name, const char *expression) :
39 TFormula()
40{
41}
42
43//------------------------------------------------------------------------------
44
45DelphesCscClusterFormula::~DelphesCscClusterFormula()
46{
47}
48
49//------------------------------------------------------------------------------
50
51Int_t DelphesCscClusterFormula::Compile(const char *expression)
52{
53 TString 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.Append(*it);
59 }
60 buffer.ReplaceAll("decayR", "x");
61 buffer.ReplaceAll("decayZ", "y");
62 buffer.ReplaceAll("Ehad", "z");
63 buffer.ReplaceAll("Eem", "t");
64
65#if ROOT_VERSION_CODE < ROOT_VERSION(6, 3, 0)
66 TFormula::SetMaxima(100000, 1000, 1000000);
67#endif
68 if(TFormula::Compile(buffer) != 0)
69 {
70 throw runtime_error("Invalid formula.");
71 }
72 return 0;
73}
74
75//------------------------------------------------------------------------------
76
77Double_t DelphesCscClusterFormula::Eval(Double_t decayR, Double_t decayZ, Double_t Ehad, Double_t Eem)
78{
79 Double_t x[4] = {decayR, decayZ, Ehad, Eem};
80 return EvalPar(x);
81}
82
83//------------------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.