Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/FastJetGridMedianEstimator.cc

    r341014c r973b92a  
    1717 */
    1818
     19
    1920/** \class FastJetGridMedianEstimator
    2021 *
     
    3132#include "classes/DelphesFormula.h"
    3233
     34#include "ExRootAnalysis/ExRootResult.h"
     35#include "ExRootAnalysis/ExRootFilter.h"
    3336#include "ExRootAnalysis/ExRootClassifier.h"
    34 #include "ExRootAnalysis/ExRootFilter.h"
    35 #include "ExRootAnalysis/ExRootResult.h"
    3637
     38#include "TMath.h"
     39#include "TString.h"
     40#include "TFormula.h"
     41#include "TRandom3.h"
     42#include "TObjArray.h"
    3743#include "TDatabasePDG.h"
    38 #include "TFormula.h"
    3944#include "TLorentzVector.h"
    40 #include "TMath.h"
    41 #include "TObjArray.h"
    42 #include "TRandom3.h"
    43 #include "TString.h"
    4445
    4546#include <algorithm>
     47#include <stdexcept>
    4648#include <iostream>
    4749#include <sstream>
    48 #include <stdexcept>
     50#include <vector>
    4951#include <utility>
    50 #include <vector>
    5152
     53#include "fastjet/PseudoJet.hh"
     54#include "fastjet/JetDefinition.hh"
    5255#include "fastjet/ClusterSequence.hh"
     56#include "fastjet/Selector.hh"
     57#include "fastjet/RectangularGrid.hh"
    5358#include "fastjet/ClusterSequenceArea.hh"
    54 #include "fastjet/JetDefinition.hh"
    55 #include "fastjet/PseudoJet.hh"
    56 #include "fastjet/RectangularGrid.hh"
    57 #include "fastjet/Selector.hh"
    5859#include "fastjet/tools/JetMedianBackgroundEstimator.hh"
    5960
    6061#include "fastjet/tools/GridMedianBackgroundEstimator.hh"
    6162
     63#include "fastjet/plugins/SISCone/fastjet/SISConePlugin.hh"
     64#include "fastjet/plugins/CDFCones/fastjet/CDFMidPointPlugin.hh"
    6265#include "fastjet/plugins/CDFCones/fastjet/CDFJetCluPlugin.hh"
    63 #include "fastjet/plugins/CDFCones/fastjet/CDFMidPointPlugin.hh"
    64 #include "fastjet/plugins/SISCone/fastjet/SISConePlugin.hh"
    6566
    66 #include "fastjet/contribs/Nsubjettiness/ExtraRecombiners.hh"
     67#include "fastjet/contribs/Nsubjettiness/Nsubjettiness.hh"
    6768#include "fastjet/contribs/Nsubjettiness/Njettiness.hh"
    6869#include "fastjet/contribs/Nsubjettiness/NjettinessPlugin.hh"
    69 #include "fastjet/contribs/Nsubjettiness/Nsubjettiness.hh"
     70#include "fastjet/contribs/Nsubjettiness/ExtraRecombiners.hh"
    7071
    7172using namespace std;
    7273using namespace fastjet;
    7374using namespace fastjet::contrib;
     75
    7476
    7577//------------------------------------------------------------------------------
     
    7880  fItInputArray(0)
    7981{
     82
    8083}
    8184
     
    8487FastJetGridMedianEstimator::~FastJetGridMedianEstimator()
    8588{
     89
    8690}
    8791
     
    100104
    101105  fEstimators.clear();
    102   for(i = 0; i < size / 4; ++i)
     106  for(i = 0; i < size/4; ++i)
    103107  {
    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();
    108112    fEstimators.push_back(new GridMedianBackgroundEstimator(rapMin, rapMax, drap, dphi));
    109113  }
     
    121125void FastJetGridMedianEstimator::Finish()
    122126{
    123   vector<GridMedianBackgroundEstimator *>::iterator itEstimators;
     127  vector< GridMedianBackgroundEstimator * >::iterator itEstimators;
    124128
    125129  for(itEstimators = fEstimators.begin(); itEstimators != fEstimators.end(); ++itEstimators)
     
    140144  Double_t rho = 0;
    141145  PseudoJet jet;
    142   vector<PseudoJet> inputList, outputList;
     146  vector< PseudoJet > inputList, outputList;
    143147
    144   vector<GridMedianBackgroundEstimator *>::iterator itEstimators;
    145   ;
     148  vector< GridMedianBackgroundEstimator * >::iterator itEstimators;;
    146149
    147150  DelphesFactory *factory = GetFactory();
     
    152155  fItInputArray->Reset();
    153156  number = 0;
    154   while((candidate = static_cast<Candidate *>(fItInputArray->Next())))
     157  while((candidate = static_cast<Candidate*>(fItInputArray->Next())))
    155158  {
    156159    momentum = candidate->Momentum;
Note: See TracChangeset for help on using the changeset viewer.