Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/Merger.cc

    r341014c rcab38f6  
    1717 */
    1818
     19
    1920/** \class Merger
    2021 *
     
    3233#include "classes/DelphesFormula.h"
    3334
     35#include "ExRootAnalysis/ExRootResult.h"
     36#include "ExRootAnalysis/ExRootFilter.h"
    3437#include "ExRootAnalysis/ExRootClassifier.h"
    35 #include "ExRootAnalysis/ExRootFilter.h"
    36 #include "ExRootAnalysis/ExRootResult.h"
    3738
     39#include "TMath.h"
     40#include "TString.h"
     41#include "TFormula.h"
     42#include "TRandom3.h"
     43#include "TObjArray.h"
    3844#include "TDatabasePDG.h"
    39 #include "TFormula.h"
    4045#include "TLorentzVector.h"
    41 #include "TMath.h"
    42 #include "TObjArray.h"
    43 #include "TRandom3.h"
    44 #include "TString.h"
    4546
    46 #include <algorithm>
     47#include <algorithm>
     48#include <stdexcept>
    4749#include <iostream>
    4850#include <sstream>
    49 #include <stdexcept>
    5051
    5152using namespace std;
     
    8889
    8990  fMomentumOutputArray = ExportArray(GetString("MomentumOutputArray", "momentum"));
    90 
     91 
    9192  fEnergyOutputArray = ExportArray(GetString("EnergyOutputArray", "energy"));
    9293}
     
    9697void Merger::Finish()
    9798{
    98   vector<TIterator *>::iterator itInputList;
     99  vector< TIterator * >::iterator itInputList;
    99100  TIterator *iterator;
    100101
     
    112113  Candidate *candidate;
    113114  TLorentzVector momentum;
    114   Double_t sumPT, sumE;
    115   vector<TIterator *>::iterator itInputList;
     115  Double_t sumPT, sumE; 
     116  vector< TIterator * >::iterator itInputList;
    116117  TIterator *iterator;
    117118
    118119  DelphesFactory *factory = GetFactory();
    119 
     120 
    120121  momentum.SetPxPyPzE(0.0, 0.0, 0.0, 0.0);
    121122  sumPT = 0;
     
    129130    // loop over all candidates
    130131    iterator->Reset();
    131     while((candidate = static_cast<Candidate *>(iterator->Next())))
     132    while((candidate = static_cast<Candidate*>(iterator->Next())))
    132133    {
    133134      const TLorentzVector &candidateMomentum = candidate->Momentum;
     
    142143
    143144  candidate = factory->NewCandidate();
    144 
     145 
    145146  candidate->Position.SetXYZT(0.0, 0.0, 0.0, 0.0);
    146147  candidate->Momentum = momentum;
    147 
     148 
    148149  fMomentumOutputArray->Add(candidate);
    149150
    150151  candidate = factory->NewCandidate();
    151 
     152 
    152153  candidate->Position.SetXYZT(0.0, 0.0, 0.0, 0.0);
    153154  candidate->Momentum.SetPtEtaPhiE(sumPT, 0.0, 0.0, sumE);
    154 
     155 
    155156  fEnergyOutputArray->Add(candidate);
    156157}
Note: See TracChangeset for help on using the changeset viewer.