Changes in modules/FastJetGridMedianEstimator.cc [341014c:973b92a] in git
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/FastJetGridMedianEstimator.cc
r341014c r973b92a 17 17 */ 18 18 19 19 20 /** \class FastJetGridMedianEstimator 20 21 * … … 31 32 #include "classes/DelphesFormula.h" 32 33 34 #include "ExRootAnalysis/ExRootResult.h" 35 #include "ExRootAnalysis/ExRootFilter.h" 33 36 #include "ExRootAnalysis/ExRootClassifier.h" 34 #include "ExRootAnalysis/ExRootFilter.h"35 #include "ExRootAnalysis/ExRootResult.h"36 37 38 #include "TMath.h" 39 #include "TString.h" 40 #include "TFormula.h" 41 #include "TRandom3.h" 42 #include "TObjArray.h" 37 43 #include "TDatabasePDG.h" 38 #include "TFormula.h"39 44 #include "TLorentzVector.h" 40 #include "TMath.h"41 #include "TObjArray.h"42 #include "TRandom3.h"43 #include "TString.h"44 45 45 46 #include <algorithm> 47 #include <stdexcept> 46 48 #include <iostream> 47 49 #include <sstream> 48 #include < stdexcept>50 #include <vector> 49 51 #include <utility> 50 #include <vector>51 52 53 #include "fastjet/PseudoJet.hh" 54 #include "fastjet/JetDefinition.hh" 52 55 #include "fastjet/ClusterSequence.hh" 56 #include "fastjet/Selector.hh" 57 #include "fastjet/RectangularGrid.hh" 53 58 #include "fastjet/ClusterSequenceArea.hh" 54 #include "fastjet/JetDefinition.hh"55 #include "fastjet/PseudoJet.hh"56 #include "fastjet/RectangularGrid.hh"57 #include "fastjet/Selector.hh"58 59 #include "fastjet/tools/JetMedianBackgroundEstimator.hh" 59 60 60 61 #include "fastjet/tools/GridMedianBackgroundEstimator.hh" 61 62 63 #include "fastjet/plugins/SISCone/fastjet/SISConePlugin.hh" 64 #include "fastjet/plugins/CDFCones/fastjet/CDFMidPointPlugin.hh" 62 65 #include "fastjet/plugins/CDFCones/fastjet/CDFJetCluPlugin.hh" 63 #include "fastjet/plugins/CDFCones/fastjet/CDFMidPointPlugin.hh"64 #include "fastjet/plugins/SISCone/fastjet/SISConePlugin.hh"65 66 66 #include "fastjet/contribs/Nsubjettiness/ ExtraRecombiners.hh"67 #include "fastjet/contribs/Nsubjettiness/Nsubjettiness.hh" 67 68 #include "fastjet/contribs/Nsubjettiness/Njettiness.hh" 68 69 #include "fastjet/contribs/Nsubjettiness/NjettinessPlugin.hh" 69 #include "fastjet/contribs/Nsubjettiness/ Nsubjettiness.hh"70 #include "fastjet/contribs/Nsubjettiness/ExtraRecombiners.hh" 70 71 71 72 using namespace std; 72 73 using namespace fastjet; 73 74 using namespace fastjet::contrib; 75 74 76 75 77 //------------------------------------------------------------------------------ … … 78 80 fItInputArray(0) 79 81 { 82 80 83 } 81 84 … … 84 87 FastJetGridMedianEstimator::~FastJetGridMedianEstimator() 85 88 { 89 86 90 } 87 91 … … 100 104 101 105 fEstimators.clear(); 102 for(i = 0; i < size /4; ++i)106 for(i = 0; i < size/4; ++i) 103 107 { 104 rapMin = param[i *4].GetDouble();105 rapMax = param[i *4 + 1].GetDouble();106 drap = param[i *4 + 2].GetDouble();107 dphi = param[i *4 + 3].GetDouble();108 rapMin = param[i*4].GetDouble(); 109 rapMax = param[i*4 + 1].GetDouble(); 110 drap = param[i*4 + 2].GetDouble(); 111 dphi = param[i*4 + 3].GetDouble(); 108 112 fEstimators.push_back(new GridMedianBackgroundEstimator(rapMin, rapMax, drap, dphi)); 109 113 } … … 121 125 void FastJetGridMedianEstimator::Finish() 122 126 { 123 vector< GridMedianBackgroundEstimator *>::iterator itEstimators;127 vector< GridMedianBackgroundEstimator * >::iterator itEstimators; 124 128 125 129 for(itEstimators = fEstimators.begin(); itEstimators != fEstimators.end(); ++itEstimators) … … 140 144 Double_t rho = 0; 141 145 PseudoJet jet; 142 vector< PseudoJet> inputList, outputList;146 vector< PseudoJet > inputList, outputList; 143 147 144 vector<GridMedianBackgroundEstimator *>::iterator itEstimators; 145 ; 148 vector< GridMedianBackgroundEstimator * >::iterator itEstimators;; 146 149 147 150 DelphesFactory *factory = GetFactory(); … … 152 155 fItInputArray->Reset(); 153 156 number = 0; 154 while((candidate = static_cast<Candidate 157 while((candidate = static_cast<Candidate*>(fItInputArray->Next()))) 155 158 { 156 159 momentum = candidate->Momentum;
Note:
See TracChangeset
for help on using the changeset viewer.