Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/FastJetGridMedianEstimator.cc

    r973b92a r341014c  
    1717 */
    1818
    19 
    2019/** \class FastJetGridMedianEstimator
    2120 *
     
    3231#include "classes/DelphesFormula.h"
    3332
     33#include "ExRootAnalysis/ExRootClassifier.h"
     34#include "ExRootAnalysis/ExRootFilter.h"
    3435#include "ExRootAnalysis/ExRootResult.h"
    35 #include "ExRootAnalysis/ExRootFilter.h"
    36 #include "ExRootAnalysis/ExRootClassifier.h"
    3736
     37#include "TDatabasePDG.h"
     38#include "TFormula.h"
     39#include "TLorentzVector.h"
    3840#include "TMath.h"
     41#include "TObjArray.h"
     42#include "TRandom3.h"
    3943#include "TString.h"
    40 #include "TFormula.h"
    41 #include "TRandom3.h"
    42 #include "TObjArray.h"
    43 #include "TDatabasePDG.h"
    44 #include "TLorentzVector.h"
    4544
    4645#include <algorithm>
    47 #include <stdexcept>
    4846#include <iostream>
    4947#include <sstream>
     48#include <stdexcept>
     49#include <utility>
    5050#include <vector>
    51 #include <utility>
    5251
     52#include "fastjet/ClusterSequence.hh"
     53#include "fastjet/ClusterSequenceArea.hh"
     54#include "fastjet/JetDefinition.hh"
    5355#include "fastjet/PseudoJet.hh"
    54 #include "fastjet/JetDefinition.hh"
    55 #include "fastjet/ClusterSequence.hh"
     56#include "fastjet/RectangularGrid.hh"
    5657#include "fastjet/Selector.hh"
    57 #include "fastjet/RectangularGrid.hh"
    58 #include "fastjet/ClusterSequenceArea.hh"
    5958#include "fastjet/tools/JetMedianBackgroundEstimator.hh"
    6059
    6160#include "fastjet/tools/GridMedianBackgroundEstimator.hh"
    6261
     62#include "fastjet/plugins/CDFCones/fastjet/CDFJetCluPlugin.hh"
     63#include "fastjet/plugins/CDFCones/fastjet/CDFMidPointPlugin.hh"
    6364#include "fastjet/plugins/SISCone/fastjet/SISConePlugin.hh"
    64 #include "fastjet/plugins/CDFCones/fastjet/CDFMidPointPlugin.hh"
    65 #include "fastjet/plugins/CDFCones/fastjet/CDFJetCluPlugin.hh"
    6665
    67 #include "fastjet/contribs/Nsubjettiness/Nsubjettiness.hh"
     66#include "fastjet/contribs/Nsubjettiness/ExtraRecombiners.hh"
    6867#include "fastjet/contribs/Nsubjettiness/Njettiness.hh"
    6968#include "fastjet/contribs/Nsubjettiness/NjettinessPlugin.hh"
    70 #include "fastjet/contribs/Nsubjettiness/ExtraRecombiners.hh"
     69#include "fastjet/contribs/Nsubjettiness/Nsubjettiness.hh"
    7170
    7271using namespace std;
    7372using namespace fastjet;
    7473using namespace fastjet::contrib;
    75 
    7674
    7775//------------------------------------------------------------------------------
     
    8078  fItInputArray(0)
    8179{
    82 
    8380}
    8481
     
    8784FastJetGridMedianEstimator::~FastJetGridMedianEstimator()
    8885{
    89 
    9086}
    9187
     
    104100
    105101  fEstimators.clear();
    106   for(i = 0; i < size/4; ++i)
     102  for(i = 0; i < size / 4; ++i)
    107103  {
    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();
    112108    fEstimators.push_back(new GridMedianBackgroundEstimator(rapMin, rapMax, drap, dphi));
    113109  }
     
    125121void FastJetGridMedianEstimator::Finish()
    126122{
    127   vector< GridMedianBackgroundEstimator * >::iterator itEstimators;
     123  vector<GridMedianBackgroundEstimator *>::iterator itEstimators;
    128124
    129125  for(itEstimators = fEstimators.begin(); itEstimators != fEstimators.end(); ++itEstimators)
     
    144140  Double_t rho = 0;
    145141  PseudoJet jet;
    146   vector< PseudoJet > inputList, outputList;
     142  vector<PseudoJet> inputList, outputList;
    147143
    148   vector< GridMedianBackgroundEstimator * >::iterator itEstimators;;
     144  vector<GridMedianBackgroundEstimator *>::iterator itEstimators;
     145  ;
    149146
    150147  DelphesFactory *factory = GetFactory();
     
    155152  fItInputArray->Reset();
    156153  number = 0;
    157   while((candidate = static_cast<Candidate*>(fItInputArray->Next())))
     154  while((candidate = static_cast<Candidate *>(fItInputArray->Next())))
    158155  {
    159156    momentum = candidate->Momentum;
Note: See TracChangeset for help on using the changeset viewer.