Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/TimeSmearing.cc

    r043e5b3 r341014c  
    1717 */
    1818
    19 
    2019/** \class TimeSmearing
    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
    4645#include <algorithm>
    47 #include <stdexcept>
    4846#include <iostream>
    4947#include <sstream>
     48#include <stdexcept>
    5049
    5150using namespace std;
     
    5453
    5554TimeSmearing::TimeSmearing() :
    56 fItInputArray(0)
     55  fItInputArray(0)
    5756{
    5857}
     
    9796
    9897  fItInputArray->Reset();
    99   while((candidate = static_cast<Candidate*>(fItInputArray->Next())))
     98  while((candidate = static_cast<Candidate *>(fItInputArray->Next())))
    10099  {
    101100    const TLorentzVector &candidateInitialPosition = candidate->InitialPosition;
    102101    const TLorentzVector &candidateFinalPosition = candidate->Position;
    103102
    104     ti = candidateInitialPosition.T()*1.0E-3/c_light;
    105     tf = candidateFinalPosition.T()*1.0E-3/c_light;
     103    ti = candidateInitialPosition.T() * 1.0E-3 / c_light;
     104    tf = candidateFinalPosition.T() * 1.0E-3 / c_light;
    106105
    107106    // apply smearing formula
     
    109108    ti = ti + tf_smeared - tf;
    110109    tf = tf_smeared;
    111    
     110
    112111    mother = candidate;
    113     candidate = static_cast<Candidate*>(candidate->Clone());
    114     candidate->InitialPosition.SetT(ti*1.0E3*c_light);
    115     candidate->Position.SetT(tf*1.0E3*c_light);
     112    candidate = static_cast<Candidate *>(candidate->Clone());
     113    candidate->InitialPosition.SetT(ti * 1.0E3 * c_light);
     114    candidate->Position.SetT(tf * 1.0E3 * c_light);
    116115
    117     candidate->ErrorT = fTimeResolution*1.0E3*c_light;
     116    candidate->ErrorT = fTimeResolution * 1.0E3 * c_light;
    118117
    119118    candidate->AddCandidate(mother);
Note: See TracChangeset for help on using the changeset viewer.