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