Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/ImpactParameterSmearing.cc

    r341014c r632da30  
    2525 */
    2626
     27
    2728#include "modules/ImpactParameterSmearing.h"
    2829
     
    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
    4546#include <algorithm>
     47#include <stdexcept>
    4648#include <iostream>
    4749#include <sstream>
    48 #include <stdexcept>
    4950
    5051using namespace std;
     
    99100
    100101  fItInputArray->Reset();
    101   while((candidate = static_cast<Candidate *>(fItInputArray->Next())))
     102  while((candidate = static_cast<Candidate*>(fItInputArray->Next())))
    102103  {
    103104
    104105    // take momentum before smearing (otherwise apply double smearing on d0)
    105     particle = static_cast<Candidate *>(candidate->GetCandidates()->At(0));
     106    particle = static_cast<Candidate*>(candidate->GetCandidates()->At(0));
    106107
    107108    const TLorentzVector &candidateMomentum = particle->Momentum;
     
    111112    phi = candidateMomentum.Phi();
    112113    e = candidateMomentum.E();
    113 
     114   
    114115    px = candidateMomentum.Px();
    115116    py = candidateMomentum.Py();
    116117
    117118    // calculate coordinates of closest approach to track circle in transverse plane xd, yd, zd
    118     xd = candidate->Xd;
    119     yd = candidate->Yd;
    120     zd = candidate->Zd;
     119    xd =  candidate->Xd;
     120    yd =  candidate->Yd;
     121    zd =  candidate->Zd;
    121122
    122123    // calculate smeared values
     
    130131
    131132    // calculate impact parameter (after-smearing)
    132     d0 = (xd * py - yd * px) / pt;
     133    d0 = (xd*py - yd*px)/pt;
    133134
    134135    dd0 = gRandom->Gaus(0.0, fFormula->Eval(pt, eta, phi, e));
     
    137138    mother = candidate;
    138139
    139     candidate = static_cast<Candidate *>(candidate->Clone());
     140    candidate = static_cast<Candidate*>(candidate->Clone());
    140141    candidate->Xd = xd;
    141142    candidate->Yd = yd;
Note: See TracChangeset for help on using the changeset viewer.