SmearPerfectDetector.cxx

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file SmearPerfectDetector.cxx

//-----------------------
//Simple Detector: perfect resolution
//Acceptance from [-15,15] in eta
// Original author: Barak Schmookler
//-----------------------

#include "eicsmear/erhic/VirtualParticle.h"
#include "eicsmear/smear/Acceptance.h"
#include "eicsmear/smear/Device.h"
#include "eicsmear/smear/Detector.h"
#include "eicsmear/smear/Smearer.h"
#include "eicsmear/smear/ParticleMCS.h"
#include "eicsmear/smear/PerfectID.h"


//Convert pseudorapidity (eta) to polar angle (theta) in radians.
//Make use of TLorentzVector to do eta-to-theta conversion.
static double etaToTheta(const double eta) {
  TLorentzVector v;
  v.SetPtEtaPhiM(1., eta, 0., 0.);
  return v.Theta();
}

//Build Detector
Smear::Detector BuildPerfectDetector(){

  //Create Devices
  Smear::Device energy(Smear::kE,"0"); //"0":Perfect Resolution
  Smear::Device momentum(Smear::kP,"0");
  Smear::Device theta(Smear::kTheta,"0");
  Smear::Device phi(Smear::kPhi,"0");

  //Detector Acceptance
  Smear::Acceptance::Zone acceptall(etaToTheta(15.),etaToTheta(-15.));

  energy.Accept.AddZone(acceptall);
  momentum.Accept.AddZone(acceptall);
  theta.Accept.AddZone(acceptall);
  phi.Accept.AddZone(acceptall);

  //PID performance is unparameterised as of now
  Smear::PerfectID pid;
  pid.Accept.AddZone(acceptall);

  //Create the detector and add devices
  Smear::Detector det;
  det.AddDevice(energy);
  det.AddDevice(momentum);
  det.AddDevice(theta);
  det.AddDevice(phi);
  det.AddDevice(pid);
  det.SetEventKinematicsCalculator("NM JB DA");

  return det;

}