Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/AngularSmearing.cc

    r95aa610 r341014c  
    1717 */
    1818
    19 
    2019/** \class AngularSmearing
    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
    46 #include <algorithm>
    47 #include <stdexcept>
     45#include <algorithm>
    4846#include <iostream>
    4947#include <sstream>
     48#include <stdexcept>
    5049
    5150using namespace std;
     
    7776  fFormulaPhi->Compile(GetString("PhiResolutionFormula", "0.0"));
    7877
    79 
    8078  // import input array
    8179
     
    103101
    104102  fItInputArray->Reset();
    105   while((candidate = static_cast<Candidate*>(fItInputArray->Next())))
     103  while((candidate = static_cast<Candidate *>(fItInputArray->Next())))
    106104  {
    107105    const TLorentzVector &candidatePosition = candidate->Position;
     
    116114    eta = gRandom->Gaus(eta, fFormulaEta->Eval(pt, eta, phi, e));
    117115    phi = gRandom->Gaus(phi, fFormulaPhi->Eval(pt, eta, phi, e));
    118    
     116
    119117    if(pt <= 0.0) continue;
    120118
    121119    mother = candidate;
    122     candidate = static_cast<Candidate*>(candidate->Clone());
     120    candidate = static_cast<Candidate *>(candidate->Clone());
    123121    eta = candidateMomentum.Eta();
    124122    phi = candidateMomentum.Phi();
    125     candidate->Momentum.SetPtEtaPhiE(pt, eta, phi, pt*TMath::CosH(eta));
     123    candidate->Momentum.SetPtEtaPhiE(pt, eta, phi, pt * TMath::CosH(eta));
    126124    candidate->AddCandidate(mother);
    127        
     125
    128126    fOutputArray->Add(candidate);
    129127  }
Note: See TracChangeset for help on using the changeset viewer.