Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/AngularSmearing.cc

    r341014c r95aa610  
    1717 */
    1818
     19
    1920/** \class AngularSmearing
    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
    45 #include <algorithm>
     46#include <algorithm>
     47#include <stdexcept>
    4648#include <iostream>
    4749#include <sstream>
    48 #include <stdexcept>
    4950
    5051using namespace std;
     
    7677  fFormulaPhi->Compile(GetString("PhiResolutionFormula", "0.0"));
    7778
     79
    7880  // import input array
    7981
     
    101103
    102104  fItInputArray->Reset();
    103   while((candidate = static_cast<Candidate *>(fItInputArray->Next())))
     105  while((candidate = static_cast<Candidate*>(fItInputArray->Next())))
    104106  {
    105107    const TLorentzVector &candidatePosition = candidate->Position;
     
    114116    eta = gRandom->Gaus(eta, fFormulaEta->Eval(pt, eta, phi, e));
    115117    phi = gRandom->Gaus(phi, fFormulaPhi->Eval(pt, eta, phi, e));
    116 
     118   
    117119    if(pt <= 0.0) continue;
    118120
    119121    mother = candidate;
    120     candidate = static_cast<Candidate *>(candidate->Clone());
     122    candidate = static_cast<Candidate*>(candidate->Clone());
    121123    eta = candidateMomentum.Eta();
    122124    phi = candidateMomentum.Phi();
    123     candidate->Momentum.SetPtEtaPhiE(pt, eta, phi, pt * TMath::CosH(eta));
     125    candidate->Momentum.SetPtEtaPhiE(pt, eta, phi, pt*TMath::CosH(eta));
    124126    candidate->AddCandidate(mother);
    125 
     127       
    126128    fOutputArray->Add(candidate);
    127129  }
Note: See TracChangeset for help on using the changeset viewer.