Fork me on GitHub

Changeset 341014c in git for modules/MomentumSmearing.cc


Ignore:
Timestamp:
Feb 12, 2019, 9:29:17 PM (5 years ago)
Author:
Pavel Demin <pavel-demin@…>
Branches:
ImprovedOutputFile, Timing, llp, master
Children:
6455202
Parents:
45e58be
Message:

apply .clang-format to all .h, .cc and .cpp files

File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/MomentumSmearing.cc

    r45e58be r341014c  
    1717 */
    1818
    19 
    2019/** \class MomentumSmearing
    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;
     
    9998
    10099  fItInputArray->Reset();
    101   while((candidate = static_cast<Candidate*>(fItInputArray->Next())))
     100  while((candidate = static_cast<Candidate *>(fItInputArray->Next())))
    102101  {
    103102    const TLorentzVector &candidatePosition = candidate->Position;
     
    108107    e = candidateMomentum.E();
    109108    res = fFormula->Eval(pt, eta, phi, e);
    110  
     109
    111110    // apply smearing formula
    112111    //pt = gRandom->Gaus(pt, fFormula->Eval(pt, eta, phi, e) * pt);
    113    
    114     res = ( res > 1.0 ) ? 1.0 : res;
    115112
    116     pt = LogNormal(pt, res * pt );
    117    
     113    res = (res > 1.0) ? 1.0 : res;
     114
     115    pt = LogNormal(pt, res * pt);
     116
    118117    //if(pt <= 0.0) continue;
    119118
    120119    mother = candidate;
    121     candidate = static_cast<Candidate*>(candidate->Clone());
     120    candidate = static_cast<Candidate *>(candidate->Clone());
    122121    eta = candidateMomentum.Eta();
    123122    phi = candidateMomentum.Phi();
    124     candidate->Momentum.SetPtEtaPhiE(pt, eta, phi, pt*TMath::CosH(eta));
     123    candidate->Momentum.SetPtEtaPhiE(pt, eta, phi, pt * TMath::CosH(eta));
    125124    //candidate->TrackResolution = fFormula->Eval(pt, eta, phi, e);
    126125    candidate->TrackResolution = res;
    127126    candidate->AddCandidate(mother);
    128        
     127
    129128    fOutputArray->Add(candidate);
    130129  }
     
    138137  if(mean > 0.0)
    139138  {
    140     b = TMath::Sqrt(TMath::Log((1.0 + (sigma*sigma)/(mean*mean))));
    141     a = TMath::Log(mean) - 0.5*b*b;
     139    b = TMath::Sqrt(TMath::Log((1.0 + (sigma * sigma) / (mean * mean))));
     140    a = TMath::Log(mean) - 0.5 * b * b;
    142141
    143     return TMath::Exp(a + b*gRandom->Gaus(0.0, 1.0));
     142    return TMath::Exp(a + b * gRandom->Gaus(0.0, 1.0));
    144143  }
    145144  else
     
    149148}
    150149
    151 
    152150//------------------------------------------------------------------------------
Note: See TracChangeset for help on using the changeset viewer.