TestStand.cxx

Go to the documentation of this file.

#include "AtlfastAlgs/MuonSmearer.h"
#include <iostream>

int main (/*int argc, char* args[]*/) {

  int randSeed = 12345;
  int lumi = 1;
  // Muon smearing key = 1 assumes standalone spectrometer muons
  // Need to give a GenParticle for ID (muSmearKey = 2) or combined 
  // (muSmearKey = 3) muons
  int muSmearKey = 1;
  std::string muonResFile = "atlfastDatafiles/MuonResolutionTable.xml";
  std::vector<double> smearParamArray;
  smearParamArray.push_back(0.065);
  smearParamArray.push_back(0.050);
  smearParamArray.push_back(0.040);
  smearParamArray.push_back(0.10);
  smearParamArray.push_back(0.245);
  smearParamArray.push_back(0.007);
  smearParamArray.push_back(0.306);
  smearParamArray.push_back(2.4);
  smearParamArray.push_back(0.228);
  smearParamArray.push_back(0.007);
  int smearParamSchema = 1;

  Atlfast::MuonSmearer smearer(randSeed, lumi, muSmearKey, muonResFile);
  smearer.setSmearParameters(smearParamArray);
  smearer.setSmearParamSchema(smearParamSchema);

  HepLorentzVector hlv(10000.0,10000.0,10000.0,20000.0);
  HepLorentzVector hlv_smear = smearer.smear(hlv);

  std::cout << "Unsmeared: " << hlv << ", smeared: " << hlv_smear << std::endl;
  
}

    

---