CTPPSDirectProtonSimulation Class Reference

Inheritance diagram for CTPPSDirectProtonSimulation:

Public Member Functions
	CTPPSDirectProtonSimulation (const edm::ParameterSet &)
	~CTPPSDirectProtonSimulation () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Member Functions
void	processProton (const HepMC::GenVertex *in_vtx, const HepMC::GenParticle *in_trk, const CTPPSGeometry &geometry, const LHCInfo &lhcInfo, const CTPPSBeamParameters &beamParameters, const LHCInterpolatedOpticalFunctionsSetCollection &opticalFunctions, CLHEP::HepRandomEngine *rndEngine, std::vector< CTPPSLocalTrackLite > &out_tracks, edm::DetSetVector< TotemRPRecHit > &out_strip_hits, edm::DetSetVector< CTPPSPixelRecHit > &out_pixel_hits, edm::DetSetVector< CTPPSDiamondRecHit > &out_diamond_hits, std::map< int, edm::DetSetVector< TotemRPRecHit >> &out_strip_hits_per_particle, std::map< int, edm::DetSetVector< CTPPSPixelRecHit >> &out_pixel_hits_per_particle, std::map< int, edm::DetSetVector< CTPPSDiamondRecHit >> &out_diamond_hits_per_particle) const
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
bool	checkApertures_
	simulation parameters More...
bool	checkIsHit_
std::unique_ptr< TF2 >	empiricalAperture45_
std::unique_ptr< TF2 >	empiricalAperture56_
edm::EDGetTokenT< edm::HepMCProduct >	hepMCToken_
double	insensitiveMarginStrips_
	size of insensitive margin at sensor's edge facing the beam, in mm More...
std::string	lhcInfoLabel_
	input More...
std::string	opticsLabel_
double	pitchPixelsHor_
double	pitchPixelsVer_
double	pitchStrips_
	strip pitch in mm More...
bool	produceHitsRelativeToBeam_
bool	produceRecHits_
bool	produceScoringPlaneHits_
	flags what output to be produced More...
bool	roundToPitch_
double	stripZeroPosition_
	internal variable: v position of strip 0, in mm More...
std::unique_ptr< TF1 >	timeResolutionDiamonds45_
std::unique_ptr< TF1 >	timeResolutionDiamonds56_
	x-dependent time resolution of diamonds (per rec hit) in ns More...
bool	useEmpiricalApertures_
unsigned int	verbosity_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Definition at line 57 of file CTPPSDirectProtonSimulation.cc.

Constructor & Destructor Documentation

CTPPSDirectProtonSimulation()

CTPPSDirectProtonSimulation::CTPPSDirectProtonSimulation ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 128 of file CTPPSDirectProtonSimulation.cc.

    : lhcInfoLabel_(iConfig.getParameter<std::string>("lhcInfoLabel")),
      opticsLabel_(iConfig.getParameter<std::string>("opticsLabel")),
      hepMCToken_(consumes<edm::HepMCProduct>(iConfig.getParameter<edm::InputTag>("hepMCTag"))),
 
      produceScoringPlaneHits_(iConfig.getParameter<bool>("produceScoringPlaneHits")),
      produceRecHits_(iConfig.getParameter<bool>("produceRecHits")),
 
      useEmpiricalApertures_(iConfig.getParameter<bool>("useEmpiricalApertures")),
      empiricalAperture45_(
          new TF2("empiricalAperture45", iConfig.getParameter<std::string>("empiricalAperture45").c_str())),
      empiricalAperture56_(
          new TF2("empiricalAperture56", iConfig.getParameter<std::string>("empiricalAperture56").c_str())),
 
      produceHitsRelativeToBeam_(iConfig.getParameter<bool>("produceHitsRelativeToBeam")),
      roundToPitch_(iConfig.getParameter<bool>("roundToPitch")),
      checkIsHit_(iConfig.getParameter<bool>("checkIsHit")),
 
      pitchStrips_(iConfig.getParameter<double>("pitchStrips")),
      insensitiveMarginStrips_(iConfig.getParameter<double>("insensitiveMarginStrips")),
 
      pitchPixelsHor_(iConfig.getParameter<double>("pitchPixelsHor")),
      pitchPixelsVer_(iConfig.getParameter<double>("pitchPixelsVer")),
 
      timeResolutionDiamonds45_(
          new TF1("timeResolutionDiamonds45", iConfig.getParameter<std::string>("timeResolutionDiamonds45").c_str())),
      timeResolutionDiamonds56_(
          new TF1("timeResolutionDiamonds56", iConfig.getParameter<std::string>("timeResolutionDiamonds56").c_str())),
 
      verbosity_(iConfig.getUntrackedParameter<unsigned int>("verbosity", 0)) {
  if (produceScoringPlaneHits_)
    produces<std::vector<CTPPSLocalTrackLite>>();
 
  if (produceRecHits_) {
    produces<edm::DetSetVector<TotemRPRecHit>>();
    produces<edm::DetSetVector<CTPPSDiamondRecHit>>();
    produces<edm::DetSetVector<CTPPSPixelRecHit>>();
 
    produces<std::map<int, edm::DetSetVector<TotemRPRecHit>>>();
    produces<std::map<int, edm::DetSetVector<CTPPSDiamondRecHit>>>();
    produces<std::map<int, edm::DetSetVector<CTPPSPixelRecHit>>>();
  }
 
  // v position of strip 0
  stripZeroPosition_ = RPTopology::last_strip_to_border_dist_ + (RPTopology::no_of_strips_ - 1) * RPTopology::pitch_ -
                       RPTopology::y_width_ / 2.;

References RPTopology::last_strip_to_border_dist_, RPTopology::no_of_strips_, RPTopology::pitch_, produceRecHits_, produceScoringPlaneHits_, stripZeroPosition_, and RPTopology::y_width_.

~CTPPSDirectProtonSimulation()

CTPPSDirectProtonSimulation::~CTPPSDirectProtonSimulation ( )

inlineoverride

Definition at line 63 of file CTPPSDirectProtonSimulation.cc.

:
  void produce(edm::Event &, const edm::EventSetup &) override;

Member Function Documentation

fillDescriptions()

void CTPPSDirectProtonSimulation::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 178 of file CTPPSDirectProtonSimulation.cc.

                                                                                             {
  edm::ParameterSetDescription desc;
  desc.addUntracked<unsigned int>("verbosity", 0);
 
  desc.add<std::string>("lhcInfoLabel", "")->setComment("label of the LHCInfo record");
  desc.add<std::string>("opticsLabel", "")->setComment("label of the optics records");
  desc.add<edm::InputTag>("hepMCTag", edm::InputTag("generator", "unsmeared"));
 
  desc.add<bool>("produceScoringPlaneHits", true);
  desc.add<bool>("produceRecHits", true);
 
  desc.add<bool>("useEmpiricalApertures", false);
  desc.add<std::string>("empiricalAperture45", "0")->setComment("2D function of xi and xangle for cutoff on sec 45");
  desc.add<std::string>("empiricalAperture56", "0")->setComment("2D function of xi and xangle for cutoff on sec 56");
 
  desc.add<bool>("produceHitsRelativeToBeam", false);
  desc.add<bool>("roundToPitch", true);
  desc.add<bool>("checkIsHit", true);
  desc.add<double>("pitchStrips", 66.e-3);              // in mm
  desc.add<double>("insensitiveMarginStrips", 34.e-3);  // in mm
 
  desc.add<double>("pitchPixelsHor", 100.e-3);
  desc.add<double>("pitchPixelsVer", 150.e-3);
 
  desc.add<std::string>("timeResolutionDiamonds45", "0.200")
      ->setComment("time resolution of single diamond sensor in sector 45, in ns");
  desc.add<std::string>("timeResolutionDiamonds56", "0.200")
      ->setComment("time resolution of single diamond sensor in sector 56, in ns");
 
  descriptions.add("ctppsDirectProtonSimulation", desc);

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, HLT_FULL_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

processProton()

void CTPPSDirectProtonSimulation::processProton ( const HepMC::GenVertex *  in_vtx, const HepMC::GenParticle *  in_trk, const CTPPSGeometry &  geometry, const LHCInfo &  lhcInfo, const CTPPSBeamParameters &  beamParameters, const LHCInterpolatedOpticalFunctionsSetCollection &  opticalFunctions, CLHEP::HepRandomEngine *  rndEngine, std::vector< CTPPSLocalTrackLite > &  out_tracks, edm::DetSetVector< TotemRPRecHit > &  out_strip_hits, edm::DetSetVector< CTPPSPixelRecHit > &  out_pixel_hits, edm::DetSetVector< CTPPSDiamondRecHit > &  out_diamond_hits, std::map< int, edm::DetSetVector< TotemRPRecHit >> &  out_strip_hits_per_particle, std::map< int, edm::DetSetVector< CTPPSPixelRecHit >> &  out_pixel_hits_per_particle, std::map< int, edm::DetSetVector< CTPPSDiamondRecHit >> &  out_diamond_hits_per_particle  ) const

private

xi is positive for diffractive protons, thus proton momentum p = (1-xi) * p_nom horizontal component of proton momentum: p_x = th_x * (1-xi) * p_nom

Definition at line 294 of file CTPPSDirectProtonSimulation.cc.

                                                                                           {
 
  std::stringstream ssLog;
 
  // vectors in CMS convention
  const HepMC::FourVector &vtx_cms = in_vtx->position();  // in mm
  const HepMC::FourVector &mom_cms = in_trk->momentum();
 
  // transformation to LHC/TOTEM convention
  HepMC::FourVector vtx_lhc(-vtx_cms.x(), vtx_cms.y(), -vtx_cms.z(), vtx_cms.t());
  HepMC::FourVector mom_lhc(-mom_cms.x(), mom_cms.y(), -mom_cms.z(), mom_cms.t());
 
  // determine the LHC arm and related parameters
  unsigned int arm = 3;
  double z_sign;
  double beamMomentum = 0.;
  double xangle = 0.;
  const std::unique_ptr<TF2> *empiricalAperture;
  if (mom_lhc.z() < 0)  // sector 45
  {
    arm = 0;
    z_sign = -1;
    beamMomentum = beamParameters.getBeamMom45();
    xangle = beamParameters.getHalfXangleX45();
    empiricalAperture = &empiricalAperture45_;
  } else {  // sector 56
    arm = 1;
    z_sign = +1;
    beamMomentum = beamParameters.getBeamMom56();
    xangle = beamParameters.getHalfXangleX56();
    empiricalAperture = &empiricalAperture56_;
  }
 
  // calculate effective RP arrival time
  // effective time mimics the timing calibration -> effective times are distributed about 0
  // units:
  //    vertex: all components in mm
  //    c_light: in mm/ns
  //    time_eff: in ns
  const double time_eff = (vtx_lhc.t() - z_sign * vtx_lhc.z()) / CLHEP::c_light;
 
  // calculate kinematics for optics parametrisation
  const double p = mom_lhc.rho();
  const double xi = 1. - p / beamMomentum;
  const double th_x_phys = mom_lhc.x() / p;
  const double th_y_phys = mom_lhc.y() / p;
  const double vtx_lhc_eff_x = vtx_lhc.x() - vtx_lhc.z() * (mom_lhc.x() / mom_lhc.z() + xangle);
  const double vtx_lhc_eff_y = vtx_lhc.y() - vtx_lhc.z() * (mom_lhc.y() / mom_lhc.z());
 
  if (verbosity_) {
    ssLog << "simu: xi = " << xi << ", th_x_phys = " << th_x_phys << ", th_y_phys = " << th_y_phys
          << ", vtx_lhc_eff_x = " << vtx_lhc_eff_x << ", vtx_lhc_eff_y = " << vtx_lhc_eff_y << std::endl;
  }
 
  // check empirical aperture
  if (useEmpiricalApertures_) {
    const auto &xangle = lhcInfo.crossingAngle();
    (*empiricalAperture)->SetParameter("xi", xi);
    (*empiricalAperture)->SetParameter("xangle", xangle);
    const double th_x_th = (*empiricalAperture)->EvalPar(nullptr);
 
    if (th_x_th > th_x_phys) {
      if (verbosity_) {
        ssLog << "stop because of empirical appertures";
        edm::LogInfo("CTPPSDirectProtonSimulation") << ssLog.str();
      }
 
      return;
    }
  }
 
  // transport the proton into each pot/scoring plane
  for (const auto &ofp : opticalFunctions) {
    CTPPSDetId rpId(ofp.first);
    const unsigned int rpDecId = rpId.arm() * 100 + rpId.station() * 10 + rpId.rp();
 
    // first check the arm
    if (rpId.arm() != arm)
      continue;
 
    if (verbosity_)
      ssLog << "  RP " << rpDecId << std::endl;
 
    // transport proton
    LHCInterpolatedOpticalFunctionsSet::Kinematics k_in = {
        vtx_lhc_eff_x * 1E-1, th_x_phys, vtx_lhc_eff_y * 1E-1, th_y_phys, xi};  // conversions: mm -> cm
    LHCInterpolatedOpticalFunctionsSet::Kinematics k_out;
    ofp.second.transport(k_in, k_out, true);
 
    double b_x = k_out.x * 1E1, b_y = k_out.y * 1E1;  // conversions: cm -> mm
    double a_x = k_out.th_x, a_y = k_out.th_y;
 
    // if needed, subtract beam position and angle
    if (produceHitsRelativeToBeam_) {
      // determine beam position
      LHCInterpolatedOpticalFunctionsSet::Kinematics k_be_in = {0., 0., 0., 0., 0.};
      LHCInterpolatedOpticalFunctionsSet::Kinematics k_be_out;
      ofp.second.transport(k_be_in, k_be_out, true);
 
      a_x -= k_be_out.th_x;
      a_y -= k_be_out.th_y;
      b_x -= k_be_out.x * 1E1;
      b_y -= k_be_out.y * 1E1;  // conversions: cm -> mm
    }
 
    const double z_scoringPlane = ofp.second.getScoringPlaneZ() * 1E1;  // conversion: cm --> mm
 
    if (verbosity_) {
      ssLog << "    proton transported: a_x = " << a_x << " rad, a_y = " << a_y << " rad, b_x = " << b_x
            << " mm, b_y = " << b_y << " mm, z = " << z_scoringPlane << " mm" << std::endl;
    }
 
    // save scoring plane hit
    if (produceScoringPlaneHits_)
      out_tracks.emplace_back(
          rpId.rawId(), b_x, 0., b_y, 0., 0., 0., 0., 0., 0., CTPPSpixelLocalTrackReconstructionInfo::invalid, 0, 0., 0.);
 
    // stop if rec hits are not to be produced
    if (!produceRecHits_)
      continue;
 
    // loop over all sensors in the RP
    for (const auto &detIdInt : geometry.sensorsInRP(rpId)) {
      CTPPSDetId detId(detIdInt);
 
      // determine the track impact point (in global coordinates)
      // !! this assumes that local axes (1, 0, 0) and (0, 1, 0) describe the sensor surface
      const auto &gl_o = geometry.localToGlobal(detId, CTPPSGeometry::Vector(0, 0, 0));
      const auto &gl_a1 = geometry.localToGlobal(detId, CTPPSGeometry::Vector(1, 0, 0)) - gl_o;
      const auto &gl_a2 = geometry.localToGlobal(detId, CTPPSGeometry::Vector(0, 1, 0)) - gl_o;
 
      double gl_o_z = gl_o.z();
      if (detId.subdetId() == CTPPSDetId::sdTimingDiamond)
        gl_o_z = -gl_o_z;  // fix bug in diamond geometry
 
      TMatrixD A(3, 3);
      TVectorD B(3);
      A(0, 0) = a_x;
      A(0, 1) = -gl_a1.x();
      A(0, 2) = -gl_a2.x();
      B(0) = gl_o.x() - b_x;
      A(1, 0) = a_y;
      A(1, 1) = -gl_a1.y();
      A(1, 2) = -gl_a2.y();
      B(1) = gl_o.y() - b_y;
      A(2, 0) = z_sign;
      A(2, 1) = -gl_a1.z();
      A(2, 2) = -gl_a2.z();
      B(2) = gl_o_z - z_scoringPlane;
      TMatrixD Ai(3, 3);
      Ai = A.Invert();
      TVectorD P(3);
      P = Ai * B;
 
      double ze = P(0);
      const CTPPSGeometry::Vector h_glo(a_x * ze + b_x, a_y * ze + b_y, z_sign * ze + z_scoringPlane);
 
      // hit in local coordinates
      auto h_loc = geometry.globalToLocal(detId, h_glo);
 
      if (verbosity_) {
        ssLog << std::endl
              << "    de z = " << P(0) << " mm, p1 = " << P(1) << " mm, p2 = " << P(2) << " mm" << std::endl
              << "    h_glo: x = " << h_glo.x() << " mm, y = " << h_glo.y() << " mm, z = " << h_glo.z() << " mm"
              << std::endl
              << "    h_loc: c1 = " << h_loc.x() << " mm, c2 = " << h_loc.y() << " mm, c3 = " << h_loc.z() << " mm"
              << std::endl;
      }
 
      // strips
      if (detId.subdetId() == CTPPSDetId::sdTrackingStrip) {
        double u = h_loc.x();
        double v = h_loc.y();
 
        if (verbosity_ > 5)
          ssLog << "            u=" << u << ", v=" << v;
 
        // is it within detector?
        if (checkIsHit_ && !RPTopology::IsHit(u, v, insensitiveMarginStrips_)) {
          if (verbosity_ > 5)
            ssLog << " | no hit" << std::endl;
          continue;
        }
 
        // round the measurement
        if (roundToPitch_) {
          double m = stripZeroPosition_ - v;
          signed int strip = (int)floor(m / pitchStrips_ + 0.5);
 
          v = stripZeroPosition_ - pitchStrips_ * strip;
 
          if (verbosity_ > 5)
            ssLog << " | strip=" << strip;
        }
 
        double sigma = pitchStrips_ / sqrt(12.);
 
        if (verbosity_ > 5)
          ssLog << " | m=" << v << ", sigma=" << sigma << std::endl;
 
        edm::DetSet<TotemRPRecHit> &hits = out_strip_hits.find_or_insert(detId);
        hits.push_back(TotemRPRecHit(v, sigma));
 
        edm::DetSet<TotemRPRecHit> &hits_per_particle =
            out_strip_hits_per_particle[in_trk->barcode()].find_or_insert(detId);
        hits_per_particle.push_back(TotemRPRecHit(v, sigma));
      }
 
      // diamonds
      if (detId.subdetId() == CTPPSDetId::sdTimingDiamond) {
        CTPPSDiamondDetId diamondDetId(detIdInt);
 
        const auto *dg = geometry.sensor(detIdInt);
 
        const auto &diamondDimensions = dg->getDiamondDimensions();
        const auto x_half_width = diamondDimensions.xHalfWidth;
        const auto y_half_width = diamondDimensions.yHalfWidth;
        const auto z_half_width = diamondDimensions.zHalfWidth;
 
        const double time_resolution = (diamondDetId.arm() == 0) ? timeResolutionDiamonds45_->Eval(h_glo.x())
                                                                 : timeResolutionDiamonds56_->Eval(h_glo.x());
 
        // check acceptance
        if (h_loc.x() < -x_half_width || h_loc.x() > +x_half_width || h_loc.y() < -y_half_width ||
            h_loc.y() > +y_half_width)
          continue;
 
        // timing information
        const double t0 = time_eff + CLHEP::RandGauss::shoot(rndEngine, 0., time_resolution);
        const double tot = 1.23456;
        const double ch_t_precis = time_resolution;
        const int time_slice = 0;
 
        // build rec hit
        const bool multiHit = false;
 
        CTPPSDiamondRecHit rc(gl_o.x(),
                              2. * x_half_width,
                              gl_o.y(),
                              2. * y_half_width,
                              gl_o_z,
                              2. * z_half_width,
                              t0,
                              tot,
                              ch_t_precis,
                              time_slice,
                              HPTDCErrorFlags(),
                              multiHit);
 
        edm::DetSet<CTPPSDiamondRecHit> &hits = out_diamond_hits.find_or_insert(detId);
        hits.push_back(rc);
 
        edm::DetSet<CTPPSDiamondRecHit> &hits_per_particle =
            out_diamond_hits_per_particle[in_trk->barcode()].find_or_insert(detId);
        hits_per_particle.push_back(rc);
      }
 
      // pixels
      if (detId.subdetId() == CTPPSDetId::sdTrackingPixel) {
        if (verbosity_) {
          CTPPSPixelDetId pixelDetId(detIdInt);
          ssLog << "    pixel plane " << pixelDetId.plane() << ": local hit x = " << h_loc.x()
                << " mm, y = " << h_loc.y() << " mm, z = " << h_loc.z() << " mm" << std::endl;
        }
 
        bool module3By2 = (geometry.sensor(detIdInt)->sensorType() != DDD_CTPPS_PIXELS_SENSOR_TYPE_2x2);
        if (checkIsHit_ && !CTPPSPixelTopology::isPixelHit(h_loc.x(), h_loc.y(), module3By2))
          continue;
 
        if (roundToPitch_) {
          h_loc.SetX(pitchPixelsHor_ * floor(h_loc.x() / pitchPixelsHor_ + 0.5));
          h_loc.SetY(pitchPixelsVer_ * floor(h_loc.y() / pitchPixelsVer_ + 0.5));
        }
 
        if (verbosity_ > 5)
          ssLog << "            hit accepted: m1 = " << h_loc.x() << " mm, m2 = " << h_loc.y() << " mm" << std::endl;
 
        const double sigmaHor = pitchPixelsHor_ / sqrt(12.);
        const double sigmaVer = pitchPixelsVer_ / sqrt(12.);
 
        const LocalPoint lp(h_loc.x(), h_loc.y(), h_loc.z());
        const LocalError le(sigmaHor, 0., sigmaVer);
 
        edm::DetSet<CTPPSPixelRecHit> &hits = out_pixel_hits.find_or_insert(detId);
        hits.push_back(CTPPSPixelRecHit(lp, le));
 
        edm::DetSet<CTPPSPixelRecHit> &hits_per_particle =
            out_pixel_hits_per_particle[in_trk->barcode()].find_or_insert(detId);
        hits_per_particle.push_back(CTPPSPixelRecHit(lp, le));
      }
    }
  }
 
  if (verbosity_)
    edm::LogInfo("CTPPSDirectProtonSimulation") << ssLog.str();

References MaterialEffects_cfi::A, CTPPSDetId::arm(), TtFullHadDaughter::B, checkIsHit_, LHCInfo::crossingAngle(), DDD_CTPPS_PIXELS_SENSOR_TYPE_2x2, empiricalAperture45_, empiricalAperture56_, edm::DetSetVector< T >::find_or_insert(), CTPPSBeamParameters::getBeamMom45(), CTPPSBeamParameters::getBeamMom56(), CTPPSBeamParameters::getHalfXangleX45(), CTPPSBeamParameters::getHalfXangleX56(), hfClusterShapes_cfi::hits, insensitiveMarginStrips_, createfilelist::int, invalid, RPTopology::IsHit(), CTPPSPixelTopology::isPixelHit(), visualization-live-secondInstance_cfg::m, OpticalFunctionsConfig_cfi::opticalFunctions, AlCaHLTBitMon_ParallelJobs::p, pitchPixelsHor_, pitchPixelsVer_, pitchStrips_, CTPPSPixelDetId::plane(), produceHitsRelativeToBeam_, produceRecHits_, produceScoringPlaneHits_, edm::DetSet< T >::push_back(), roundToPitch_, year_2016_postTS2_cff::rpId, CTPPSDetId::sdTimingDiamond, CTPPSDetId::sdTrackingPixel, CTPPSDetId::sdTrackingStrip, mathSSE::sqrt(), digitizers_cfi::strip, stripZeroPosition_, DetId::subdetId(), FrontierCondition_GT_autoExpress_cfi::t0, LHCInterpolatedOpticalFunctionsSet::Kinematics::th_x, LHCInterpolatedOpticalFunctionsSet::Kinematics::th_y, timeResolutionDiamonds45_, timeResolutionDiamonds56_, useEmpiricalApertures_, findQualityFiles::v, verbosity_, LHCInterpolatedOpticalFunctionsSet::Kinematics::x, OpticalFunctionsConfig_cfi::xangle, hybridSuperClusters_cfi::xi, and LHCInterpolatedOpticalFunctionsSet::Kinematics::y.

Referenced by produce().

produce()

void CTPPSDirectProtonSimulation::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 212 of file CTPPSDirectProtonSimulation.cc.

                                                                                       {
  // get input
  edm::Handle<edm::HepMCProduct> hepmc_prod;
  iEvent.getByToken(hepMCToken_, hepmc_prod);
 
  // get conditions
  edm::ESHandle<LHCInfo> hLHCInfo;
  iSetup.get<LHCInfoRcd>().get(lhcInfoLabel_, hLHCInfo);
 
  edm::ESHandle<CTPPSBeamParameters> hBeamParameters;
  iSetup.get<CTPPSBeamParametersRcd>().get(hBeamParameters);
 
  edm::ESHandle<LHCInterpolatedOpticalFunctionsSetCollection> hOpticalFunctions;
  iSetup.get<CTPPSInterpolatedOpticsRcd>().get(opticsLabel_, hOpticalFunctions);
 
  edm::ESHandle<CTPPSGeometry> geometry;
  iSetup.get<VeryForwardMisalignedGeometryRecord>().get(geometry);
 
  // prepare outputs
  std::unique_ptr<edm::DetSetVector<TotemRPRecHit>> pStripRecHits(new edm::DetSetVector<TotemRPRecHit>());
  std::unique_ptr<edm::DetSetVector<CTPPSDiamondRecHit>> pDiamondRecHits(new edm::DetSetVector<CTPPSDiamondRecHit>());
  std::unique_ptr<edm::DetSetVector<CTPPSPixelRecHit>> pPixelRecHits(new edm::DetSetVector<CTPPSPixelRecHit>());
 
  auto pStripRecHitsPerParticle = std::make_unique<std::map<int, edm::DetSetVector<TotemRPRecHit>>>();
  auto pDiamondRecHitsPerParticle = std::make_unique<std::map<int, edm::DetSetVector<CTPPSDiamondRecHit>>>();
  auto pPixelRecHitsPerParticle = std::make_unique<std::map<int, edm::DetSetVector<CTPPSPixelRecHit>>>();
 
  std::unique_ptr<std::vector<CTPPSLocalTrackLite>> pTracks(new std::vector<CTPPSLocalTrackLite>());
 
  // get random engine
  edm::Service<edm::RandomNumberGenerator> rng;
  CLHEP::HepRandomEngine *engine = &rng->getEngine(iEvent.streamID());
 
  // loop over event vertices
  auto evt = hepmc_prod->GetEvent();
  for (auto it_vtx = evt->vertices_begin(); it_vtx != evt->vertices_end(); ++it_vtx) {
    auto vtx = *(it_vtx);
 
    // loop over outgoing particles
    for (auto it_part = vtx->particles_out_const_begin(); it_part != vtx->particles_out_const_end(); ++it_part) {
      auto part = *(it_part);
 
      // accept only stable protons
      if (part->pdg_id() != 2212)
        continue;
 
      if (part->status() != 1 && part->status() < 83)
        continue;
 
      processProton(vtx,
                    part,
                    *geometry,
                    *hLHCInfo,
                    *hBeamParameters,
                    *hOpticalFunctions,
                    engine,
                    *pTracks,
                    *pStripRecHits,
                    *pPixelRecHits,
                    *pDiamondRecHits,
                    *pStripRecHitsPerParticle,
                    *pPixelRecHitsPerParticle,
                    *pDiamondRecHitsPerParticle);
    }
  }
 
  if (produceScoringPlaneHits_)
    iEvent.put(std::move(pTracks));
 
  if (produceRecHits_) {
    iEvent.put(std::move(pStripRecHits));
    iEvent.put(std::move(pPixelRecHits));
    iEvent.put(std::move(pDiamondRecHits));
 
    iEvent.put(std::move(pStripRecHitsPerParticle));
    iEvent.put(std::move(pPixelRecHitsPerParticle));
    iEvent.put(std::move(pDiamondRecHitsPerParticle));
  }

References geometry, edm::EventSetup::get(), get, edm::RandomNumberGenerator::getEngine(), edm::HepMCProduct::GetEvent(), hepMCToken_, iEvent, lhcInfoLabel_, eostools::move(), opticsLabel_, processProton(), produceRecHits_, produceScoringPlaneHits_, and extraflags_cff::vtx.

Member Data Documentation

checkApertures_

bool CTPPSDirectProtonSimulation::checkApertures_

private

simulation parameters

Definition at line 101 of file CTPPSDirectProtonSimulation.cc.

checkIsHit_

bool CTPPSDirectProtonSimulation::checkIsHit_

private

Definition at line 110 of file CTPPSDirectProtonSimulation.cc.

Referenced by processProton().

empiricalAperture45_

std::unique_ptr<TF2> CTPPSDirectProtonSimulation::empiricalAperture45_

private

Definition at line 105 of file CTPPSDirectProtonSimulation.cc.

Referenced by processProton().

empiricalAperture56_

std::unique_ptr<TF2> CTPPSDirectProtonSimulation::empiricalAperture56_

private

Definition at line 106 of file CTPPSDirectProtonSimulation.cc.

Referenced by processProton().

hepMCToken_

edm::EDGetTokenT<edm::HepMCProduct> CTPPSDirectProtonSimulation::hepMCToken_

private

Definition at line 94 of file CTPPSDirectProtonSimulation.cc.

Referenced by produce().

insensitiveMarginStrips_

double CTPPSDirectProtonSimulation::insensitiveMarginStrips_

private

size of insensitive margin at sensor's edge facing the beam, in mm

Definition at line 113 of file CTPPSDirectProtonSimulation.cc.

Referenced by processProton().

lhcInfoLabel_

std::string CTPPSDirectProtonSimulation::lhcInfoLabel_

private

input

Definition at line 91 of file CTPPSDirectProtonSimulation.cc.

Referenced by produce().

opticsLabel_

std::string CTPPSDirectProtonSimulation::opticsLabel_

private

Definition at line 92 of file CTPPSDirectProtonSimulation.cc.

Referenced by produce().

pitchPixelsHor_

double CTPPSDirectProtonSimulation::pitchPixelsHor_

private

Definition at line 115 of file CTPPSDirectProtonSimulation.cc.

Referenced by processProton().

pitchPixelsVer_

double CTPPSDirectProtonSimulation::pitchPixelsVer_

private

Definition at line 116 of file CTPPSDirectProtonSimulation.cc.

Referenced by processProton().

pitchStrips_

double CTPPSDirectProtonSimulation::pitchStrips_

private

strip pitch in mm

Definition at line 112 of file CTPPSDirectProtonSimulation.cc.

Referenced by processProton().

produceHitsRelativeToBeam_

bool CTPPSDirectProtonSimulation::produceHitsRelativeToBeam_

private

Definition at line 108 of file CTPPSDirectProtonSimulation.cc.

Referenced by processProton().

produceRecHits_

bool CTPPSDirectProtonSimulation::produceRecHits_

private

Definition at line 98 of file CTPPSDirectProtonSimulation.cc.

Referenced by CTPPSDirectProtonSimulation(), processProton(), and produce().

produceScoringPlaneHits_

bool CTPPSDirectProtonSimulation::produceScoringPlaneHits_

private

flags what output to be produced

Definition at line 97 of file CTPPSDirectProtonSimulation.cc.

Referenced by CTPPSDirectProtonSimulation(), processProton(), and produce().

roundToPitch_

bool CTPPSDirectProtonSimulation::roundToPitch_

private

Definition at line 109 of file CTPPSDirectProtonSimulation.cc.

Referenced by processProton().

stripZeroPosition_

double CTPPSDirectProtonSimulation::stripZeroPosition_

private

internal variable: v position of strip 0, in mm

Definition at line 126 of file CTPPSDirectProtonSimulation.cc.

Referenced by CTPPSDirectProtonSimulation(), and processProton().

timeResolutionDiamonds45_

std::unique_ptr<TF1> CTPPSDirectProtonSimulation::timeResolutionDiamonds45_

private

Definition at line 118 of file CTPPSDirectProtonSimulation.cc.

Referenced by processProton().

timeResolutionDiamonds56_

std::unique_ptr<TF1> CTPPSDirectProtonSimulation::timeResolutionDiamonds56_

private

x-dependent time resolution of diamonds (per rec hit) in ns

Definition at line 118 of file CTPPSDirectProtonSimulation.cc.

Referenced by processProton().

useEmpiricalApertures_

bool CTPPSDirectProtonSimulation::useEmpiricalApertures_

private

Definition at line 103 of file CTPPSDirectProtonSimulation.cc.

Referenced by processProton().

verbosity_

unsigned int CTPPSDirectProtonSimulation::verbosity_

private

Definition at line 121 of file CTPPSDirectProtonSimulation.cc.

Referenced by processProton().