CTPPSFastTrackingProducer Class Reference

#include <FastSimulation/CTPPSFastTrackingProducer/plugins/CTPPSFastTrackingProducer.cc>

Inheritance diagram for CTPPSFastTrackingProducer:

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

Private Types
typedef std::vector< CTPPSFastRecHit >	CTPPSFastRecHitContainer

Private Member Functions
void	FastReco (int Direction, H_RecRPObject *station)
void	MatchCellId (int cellId, std::vector< int > vrecCellId, std::vector< double > vrecTof, bool &match, double &recTof)
void	produce (edm::Event &, const edm::EventSetup &) override
void	ProjectToToF (const double x1, const double y1, const double x2, const double y2, double &xt, double &yt)
void	ReadRecHits (edm::Handle< CTPPSFastRecHitContainer > &)
void	ReconstructArm (H_RecRPObject *pps_station, double x1, double y1, double x2, double y2, double &tx, double &ty, double &eloss)
void	Reconstruction ()
bool	SearchTrack (int, int, int Direction, double &xi, double &t, double &partP, double &pt, double &thx, double &thy, double &x0, double &y0, double &xt, double &yt, double &X1d, double &Y1d, double &X2d, double &Y2d)
bool	SetBeamLine ()
void	TrackerStationClear ()
void	TrackerStationStarting ()

Private Attributes
edm::EDGetTokenT< CTPPSFastRecHitContainer >	_recHitToken
std::string	beam1filename
std::string	beam2filename
std::unique_ptr< CTPPSTrkDetector >	det1B
std::unique_ptr< CTPPSTrkDetector >	det1F
std::unique_ptr< CTPPSTrkDetector >	det2B
std::unique_ptr< CTPPSTrkDetector >	det2F
std::unique_ptr< CTPPSToFDetector >	detToF_B
std::unique_ptr< CTPPSToFDetector >	detToF_F
double	fBeamEnergy
double	fBeamMomentum
double	fBeamXRMS_ToF
double	fBeamXRMS_Trk1
double	fBeamXRMS_Trk2
bool	fCrossAngleCorr
double	fCrossingAngleBeam1
double	fCrossingAngleBeam2
double	fImpParcut
double	fTimeSigma
double	fToFCellHeight
std::vector< double >	fToFCellWidth
double	fToFInsertion
int	fToFNCellX
int	fToFNCellY
double	fToFPitchX
double	fToFPitchY
double	fToFXOffset
double	fTrackerHeight
double	fTrackerInsertion
double	fTrackerWidth
double	fTrk1XOffset
double	fTrk2XOffset
double	fz_timing
double	fz_tracker1
double	fz_tracker2
double	lengthctpps
std::unique_ptr< H_BeamLine >	m_beamlineCTPPS1
std::unique_ptr< H_BeamLine >	m_beamlineCTPPS2
bool	m_verbosity
std::unique_ptr< H_RecRPObject >	pps_stationB
std::unique_ptr< H_RecRPObject >	pps_stationF
std::vector< int >	recCellId_B
std::vector< int >	recCellId_F
std::vector< double >	recTof_B
std::vector< double >	recTof_F
std::vector< CTPPSFastTrack >	theCTPPSFastTrack
CTPPSFastTrack	track
std::unique_ptr< CTPPSTrkStation >	TrkStation_B
std::unique_ptr< CTPPSTrkStation >	TrkStation_F

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

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 50 of file CTPPSFastTrackingProducer.cc.

Member Typedef Documentation

CTPPSFastRecHitContainer

typedef std::vector<CTPPSFastRecHit> CTPPSFastTrackingProducer::CTPPSFastRecHitContainer

private

Definition at line 59 of file CTPPSFastTrackingProducer.cc.

Constructor & Destructor Documentation

CTPPSFastTrackingProducer()

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

explicit

Definition at line 136 of file CTPPSFastTrackingProducer.cc.

    : m_verbosity(false), fBeamMomentum(0.), fCrossAngleCorr(false), fCrossingAngleBeam1(0.), fCrossingAngleBeam2(0.) {
  //register your products
  produces<edm::CTPPSFastTrackContainer>("CTPPSFastTrack");
  using namespace edm;
  _recHitToken = consumes<CTPPSFastRecHitContainer>(iConfig.getParameter<edm::InputTag>("recHitTag"));
  m_verbosity = iConfig.getParameter<bool>("Verbosity");
  // User definitons
 
  // Read beam parameters needed for Hector reconstruction
  lengthctpps = iConfig.getParameter<double>("BeamLineLengthCTPPS");
  beam1filename = iConfig.getParameter<string>("Beam1");
  beam2filename = iConfig.getParameter<string>("Beam2");
  fBeamEnergy = iConfig.getParameter<double>("BeamEnergy");  // beam energy in GeV
  fBeamMomentum = sqrt(fBeamEnergy * fBeamEnergy - PPSTools::ProtonMassSQ);
  fCrossingAngleBeam1 = iConfig.getParameter<double>("CrossingAngleBeam1");
  fCrossingAngleBeam2 = iConfig.getParameter<double>("CrossingAngleBeam2");
 
  if (fCrossingAngleBeam1 != 0 || fCrossingAngleBeam2 != 0)
    fCrossAngleCorr = true;
  //Read detectors positions and parameters
 
  fz_tracker1 = iConfig.getParameter<double>("Z_Tracker1");
  fz_tracker2 = iConfig.getParameter<double>("Z_Tracker2");
  fz_timing = iConfig.getParameter<double>("Z_Timing");
  //
  fTrackerWidth = iConfig.getParameter<double>("TrackerWidth");
  fTrackerHeight = iConfig.getParameter<double>("TrackerHeight");
  fTrackerInsertion = iConfig.getParameter<double>("TrackerInsertion");
  fBeamXRMS_Trk1 = iConfig.getParameter<double>("BeamXRMS_Trk1");
  fBeamXRMS_Trk2 = iConfig.getParameter<double>("BeamXRMS_Trk2");
  fTrk1XOffset = iConfig.getParameter<double>("Trk1XOffset");
  fTrk2XOffset = iConfig.getParameter<double>("Trk2XOffset");
  fToFCellWidth = iConfig.getUntrackedParameter<std::vector<double> >("ToFCellWidth");
  fToFCellHeight = iConfig.getParameter<double>("ToFCellHeight");
  fToFPitchX = iConfig.getParameter<double>("ToFPitchX");
  fToFPitchY = iConfig.getParameter<double>("ToFPitchY");
  fToFNCellX = iConfig.getParameter<int>("ToFNCellX");
  fToFNCellY = iConfig.getParameter<int>("ToFNCellY");
  fToFInsertion = iConfig.getParameter<double>("ToFInsertion");
  fBeamXRMS_ToF = iConfig.getParameter<double>("BeamXRMS_ToF");
  fToFXOffset = iConfig.getParameter<double>("ToFXOffset");
  fTimeSigma = iConfig.getParameter<double>("TimeSigma");
  fImpParcut = iConfig.getParameter<double>("ImpParcut");
 
  if (!SetBeamLine()) {
    if (m_verbosity)
      LogDebug("CTPPSFastTrackingProducer") << "CTPPSFastTrackingProducer: WARNING: lengthctpps=  " << lengthctpps;
    return;
  }
 
  // Create a particle to get the beam energy from the beam file
  // Take care: the z inside the station is in meters
  //
  //Tracker Detector Description
  det1F = std::make_unique<CTPPSTrkDetector>(
      fTrackerWidth, fTrackerHeight, fTrackerInsertion * fBeamXRMS_Trk1 + fTrk1XOffset);
  det2F = std::make_unique<CTPPSTrkDetector>(
      fTrackerWidth, fTrackerHeight, fTrackerInsertion * fBeamXRMS_Trk2 + fTrk2XOffset);
  det1B = std::make_unique<CTPPSTrkDetector>(
      fTrackerWidth, fTrackerHeight, fTrackerInsertion * fBeamXRMS_Trk1 + fTrk1XOffset);
  det2B = std::make_unique<CTPPSTrkDetector>(
      fTrackerWidth, fTrackerHeight, fTrackerInsertion * fBeamXRMS_Trk2 + fTrk2XOffset);
 
  //Timing Detector Description
  std::vector<double> vToFCellWidth;
  vToFCellWidth.reserve(8);
  for (int i = 0; i < 8; i++) {
    vToFCellWidth.push_back(fToFCellWidth[i]);
  }
  double pos_tof = fToFInsertion * fBeamXRMS_ToF + fToFXOffset;
  detToF_F = std::make_unique<CTPPSToFDetector>(
      fToFNCellX, fToFNCellY, vToFCellWidth, fToFCellHeight, fToFPitchX, fToFPitchY, pos_tof, fTimeSigma);
  detToF_B = std::make_unique<CTPPSToFDetector>(
      fToFNCellX, fToFNCellY, vToFCellWidth, fToFCellHeight, fToFPitchX, fToFPitchY, pos_tof, fTimeSigma);
  //
}

References _recHitToken, beam1filename, beam2filename, det1B, det1F, det2B, det2F, detToF_B, detToF_F, fBeamEnergy, fBeamMomentum, fBeamXRMS_ToF, fBeamXRMS_Trk1, fBeamXRMS_Trk2, fCrossAngleCorr, fCrossingAngleBeam1, fCrossingAngleBeam2, fImpParcut, fTimeSigma, fToFCellHeight, fToFCellWidth, fToFInsertion, fToFNCellX, fToFNCellY, fToFPitchX, fToFPitchY, fToFXOffset, fTrackerHeight, fTrackerInsertion, fTrackerWidth, fTrk1XOffset, fTrk2XOffset, fz_timing, fz_tracker1, fz_tracker2, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), mps_fire::i, lengthctpps, LogDebug, m_verbosity, PPSTools::ProtonMassSQ, SetBeamLine(), and mathSSE::sqrt().

Member Function Documentation

FastReco()

void CTPPSFastTrackingProducer::FastReco ( int  Direction, H_RecRPObject *  station  )

private

Definition at line 435 of file CTPPSFastTrackingProducer.cc.

                                                                              {
  double theta = 0.;
  double xi, t, partP, pt, phi, x0, y0, thx, thy, xt, yt, X1d, Y1d, X2d, Y2d;
  CTPPSTrkDetector* Trk1 = nullptr;
  CTPPSTrkDetector* Trk2 = nullptr;
  double pos_tof = fToFInsertion * fBeamXRMS_ToF + fToFXOffset;
  int cellId = 0;
  std::vector<double> vToFCellWidth;
  vToFCellWidth.reserve(8);
  for (int i = 0; i < 8; i++) {
    vToFCellWidth.push_back(fToFCellWidth[i]);
  }
  CTPPSToFDetector* ToF = new CTPPSToFDetector(
      fToFNCellX, fToFNCellY, vToFCellWidth, fToFCellHeight, fToFPitchX, fToFPitchY, pos_tof, fTimeSigma);
  if (Direction > 0) {
    Trk1 = &(TrkStation_F->first);
    Trk2 = &(TrkStation_F->second);
  } else {
    Trk1 = &(TrkStation_B->first);
    Trk2 = &(TrkStation_B->second);
  }
  // Make a track from EVERY pair of hits combining Tracker1 and Tracker2.
  // The tracks may not be independent as 1 hit may belong to more than 1 track.
  for (int i = 0; i < (int)Trk1->ppsNHits_; i++) {
    for (int j = 0; j < (int)Trk2->ppsNHits_; j++) {
      if (SearchTrack(i, j, Direction, xi, t, partP, pt, thx, thy, x0, y0, xt, yt, X1d, Y1d, X2d, Y2d)) {
        // Check if the hitted timing cell matches the reconstructed track
        cellId = ToF->findCellId(xt, yt);
        double recTof = 0.;
        bool matchCellId = false;
        if (Direction > 0) {
          theta = sqrt(thx * thx + thy * thy) * urad;
          MatchCellId(cellId, recCellId_F, recTof_F, matchCellId, recTof);
        } else if (Direction < 0) {
          theta = TMath::Pi() - sqrt(thx * thx + thy * thy) * urad;
          MatchCellId(cellId, recCellId_B, recTof_B, matchCellId, recTof);
        }
        phi = atan2(thy, thx);  // at this point, thx is already in the cms ref. frame
 
        double px = partP * sin(theta) * cos(phi);
        double py = partP * sin(theta) * sin(phi);
        double pz = partP * cos(theta);
        double e = sqrt(partP * partP + PPSTools::ProtonMassSQ);
        TLorentzVector p(px, py, pz, e);
        // Invert the Lorentz boost made to take into account the crossing angle during simulation
        if (fCrossAngleCorr) {
          PPSTools::LorentzBoost(p, "MC", {fCrossingAngleBeam1, fCrossingAngleBeam2, fBeamMomentum, fBeamEnergy});
        }
        //Getting the Xi and t (squared four momentum transferred) of the reconstructed track
        PPSTools::Get_t_and_xi(const_cast<TLorentzVector*>(&p), t, xi, {fBeamMomentum, fBeamEnergy});
        double pxx = p.Px();
        double pyy = p.Py();
        double pzz = p.Pz();
        math::XYZVector momentum(pxx, pyy, pzz);
        math::XYZPoint vertex(x0, y0, 0);
 
        track.setp(momentum);
        track.setvertex(vertex);
        track.sett(t);
        track.setxi(xi);
        track.setx1(X1d);
        track.sety1(Y1d);
        track.setx2(X2d);
        track.sety2(Y2d);
        if (matchCellId) {
          track.setcellid(cellId);
          track.settof(recTof);
        } else {
          track.setcellid(0);
          track.settof(0.);
        }
        theCTPPSFastTrack.push_back(track);
      }
    }
  }
}  //end FastReco

References funct::cos(), MillePedeFileConverter_cfg::e, fBeamEnergy, fBeamMomentum, fBeamXRMS_ToF, fCrossAngleCorr, fCrossingAngleBeam1, fCrossingAngleBeam2, CTPPSToFDetector::findCellId(), fTimeSigma, fToFCellHeight, fToFCellWidth, fToFInsertion, fToFNCellX, fToFNCellY, fToFPitchX, fToFPitchY, fToFXOffset, PPSTools::Get_t_and_xi(), mps_fire::i, createfilelist::int, dqmiolumiharvest::j, PPSTools::LorentzBoost(), MatchCellId(), AlCaHLTBitMon_ParallelJobs::p, phi, Pi, CTPPSTrkDetector::ppsNHits_, PPSTools::ProtonMassSQ, DiDispStaMuonMonitor_cfi::pt, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, recCellId_B, recCellId_F, recTof_B, recTof_F, SearchTrack(), CTPPSFastTrack::setcellid(), CTPPSFastTrack::setp(), CTPPSFastTrack::sett(), CTPPSFastTrack::settof(), CTPPSFastTrack::setvertex(), CTPPSFastTrack::setx1(), CTPPSFastTrack::setx2(), CTPPSFastTrack::setxi(), CTPPSFastTrack::sety1(), CTPPSFastTrack::sety2(), funct::sin(), mathSSE::sqrt(), submitPVValidationJobs::t, theCTPPSFastTrack, theta(), track, TrkStation_B, TrkStation_F, urad, bphysicsOniaDQM_cfi::vertex, and protons_cff::xi.

Referenced by Reconstruction().

MatchCellId()

void CTPPSFastTrackingProducer::MatchCellId ( int  cellId, std::vector< int >  vrecCellId, std::vector< double >  vrecTof, bool &  match, double &  recTof  )

private

Definition at line 424 of file CTPPSFastTrackingProducer.cc.

                                                                                                   {
  for (unsigned int i = 0; i < vrecCellId.size(); i++) {
    if (cellId == vrecCellId.at(i)) {
      match = true;
      recTof = vrecTof.at(i);
      continue;
    }
  }
}

References mps_fire::i, and match().

Referenced by FastReco().

produce()

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

overrideprivate

Definition at line 215 of file CTPPSFastTrackingProducer.cc.

                                                                                     {
  using namespace edm;
  TrackerStationStarting();
  Handle<CTPPSFastRecHitContainer> recHits;
  iEvent.getByToken(_recHitToken, recHits);
  recCellId_F.clear();
  recCellId_B.clear();
  recTof_F.clear();
  recTof_B.clear();
  ReadRecHits(recHits);
  Reconstruction();
  TrackerStationClear();
 
  std::unique_ptr<CTPPSFastTrackContainer> output_tracks(new CTPPSFastTrackContainer);
  int n = 0;
  for (std::vector<CTPPSFastTrack>::const_iterator i = theCTPPSFastTrack.begin(); i != theCTPPSFastTrack.end(); i++) {
    output_tracks->push_back(*i);
    n += 1;
  }
 
  iEvent.put(std::move(output_tracks), "CTPPSFastTrack");
}  //end

References _recHitToken, mps_fire::i, iEvent, eostools::move(), dqmiodumpmetadata::n, ReadRecHits(), recCellId_B, recCellId_F, FastTrackerRecHitMaskProducer_cfi::recHits, Reconstruction(), recTof_B, recTof_F, theCTPPSFastTrack, TrackerStationClear(), and TrackerStationStarting().

ProjectToToF()

void CTPPSFastTrackingProducer::ProjectToToF ( const double  x1, const double  y1, const double  x2, const double  y2, double &  xt, double &  yt  )

inlineprivate

Definition at line 86 of file CTPPSFastTrackingProducer.cc.

                                                                                                                {
    xt = ((fz_timing - fz_tracker2) * (x2 - x1) / (fz_tracker2 - fz_tracker1)) + x2;
    yt = ((fz_timing - fz_tracker2) * (y2 - y1) / (fz_tracker2 - fz_tracker1)) + y2;
  };

References fz_timing, fz_tracker1, fz_tracker2, testProducerWithPsetDescEmpty_cfi::x1, testProducerWithPsetDescEmpty_cfi::x2, testProducerWithPsetDescEmpty_cfi::y1, and testProducerWithPsetDescEmpty_cfi::y2.

Referenced by SearchTrack().

ReadRecHits()

void CTPPSFastTrackingProducer::ReadRecHits ( edm::Handle< CTPPSFastRecHitContainer > &  recHits )

private

Definition at line 256 of file CTPPSFastTrackingProducer.cc.

                                                                                        {
  // DetId codification for PSimHit taken from CTPPSPixel- It will be replaced by CTPPSDetId
  // 2014314496 -> Tracker1 zPositive
  // 2014838784 -> Tracker2 zPositive
  // 2046820352 -> Timing   zPositive
  // 2031091712 -> Tracker1 zNegative
  // 2031616000 -> Tracker2 zNegative
  // 2063597568 -> Timing   zNegative
 
  for (unsigned int irecHits = 0; irecHits < recHits->size(); ++irecHits) {
    const CTPPSFastRecHit* recHitDet = &(*recHits)[irecHits];
    unsigned int detlayerId = recHitDet->detUnitId();
    double x = recHitDet->entryPoint().x();
    double y = recHitDet->entryPoint().y();
    double z = recHitDet->entryPoint().z();
    float tof = recHitDet->tof();
    if (detlayerId == 2014314496)
      det1F->AddHit(detlayerId, x, y, z);
    else if (detlayerId == 2014838784)
      det2F->AddHit(detlayerId, x, y, z);
    else if (detlayerId == 2031091712)
      det1B->AddHit(detlayerId, x, y, z);
    else if (detlayerId == 2031616000)
      det2B->AddHit(detlayerId, x, y, z);
    else if (detlayerId == 2046820352) {
      detToF_F->AddHit(x, y, tof);
      recCellId_F.push_back(detToF_F->findCellId(x, y));
      recTof_F.push_back(tof);
    } else if (detlayerId == 2063597568) {
      detToF_B->AddHit(x, y, tof);
      recCellId_B.push_back(detToF_B->findCellId(x, y));
      recTof_B.push_back(tof);
    }
 
  }  //LOOP TRK
  //creating Stations
  TrkStation_F = std::make_unique<CTPPSTrkStation>(*det1F, *det2F);
  TrkStation_B = std::make_unique<CTPPSTrkStation>(*det1B, *det2B);
}  // end function

References det1B, det1F, det2B, det2F, detToF_B, detToF_F, recCellId_B, recCellId_F, recHitDet(), FastTrackerRecHitMaskProducer_cfi::recHits, recTof_B, recTof_F, TrkStation_B, TrkStation_F, x, y, and z.

Referenced by produce().

ReconstructArm()

void CTPPSFastTrackingProducer::ReconstructArm ( H_RecRPObject *  pps_station, double  x1, double  y1, double  x2, double  y2, double &  tx, double &  ty, double &  eloss  )

private

Definition at line 404 of file CTPPSFastTrackingProducer.cc.

                                                                                                                   {
  tx = 0.;
  ty = 0.;
  eloss = 0.;
  if (!pps_station)
    return;
  x1 *= mm_to_um;
  x2 *= mm_to_um;
  y1 *= mm_to_um;
  y2 *= mm_to_um;
  pps_station->setPositions(x1, y1, x2, y2);
  double energy = pps_station->getE(AM);  // dummy call needed to calculate some Hector internal parameter
  if (edm::isNotFinite(energy))
    return;
  tx = pps_station->getTXIP();  // change orientation to CMS
  ty = pps_station->getTYIP();
  eloss = pps_station->getE();
}

References HCALHighEnergyHPDFilter_cfi::energy, edm::isNotFinite(), mm_to_um, testProducerWithPsetDescEmpty_cfi::x1, testProducerWithPsetDescEmpty_cfi::x2, testProducerWithPsetDescEmpty_cfi::y1, and testProducerWithPsetDescEmpty_cfi::y2.

Referenced by SearchTrack().

Reconstruction()

void CTPPSFastTrackingProducer::Reconstruction ( )

private

Definition at line 296 of file CTPPSFastTrackingProducer.cc.

                                               {
  theCTPPSFastTrack.clear();
  int Direction;
  Direction = 1;  //cms positive Z / forward
  FastReco(Direction, &*pps_stationF);
  Direction = -1;  //cms negative Z / backward
  FastReco(Direction, &*pps_stationB);
}  //end Reconstruction

References FastReco(), pps_stationB, pps_stationF, and theCTPPSFastTrack.

Referenced by produce().

SearchTrack()

bool CTPPSFastTrackingProducer::SearchTrack ( int  i, int  j, int  Direction, double &  xi, double &  t, double &  partP, double &  pt, double &  thx, double &  thy, double &  x0, double &  y0, double &  xt, double &  yt, double &  X1d, double &  Y1d, double &  X2d, double &  Y2d  )

private

Definition at line 305 of file CTPPSFastTrackingProducer.cc.

                                                         {
  // Given 1 hit in Tracker1 and 1 hit in Tracker2 try to make a track with Hector
  xi = 0;
  t = 0;
  partP = 0;
  pt = 0;
  x0 = 0.;
  y0 = 0.;
  xt = 0.;
  yt = 0.;
  X1d = 0.;
  Y1d = 0.;
  X2d = 0.;
  Y2d = 0.;
  CTPPSTrkDetector* det1 = nullptr;
  CTPPSTrkDetector* det2 = nullptr;
  H_RecRPObject* station = nullptr;
  // Separate in forward and backward stations according to direction
  if (Direction > 0) {
    det1 = &(TrkStation_F->first);
    det2 = &(TrkStation_F->second);
    station = &*pps_stationF;
  } else {
    det1 = &(TrkStation_B->first);
    det2 = &(TrkStation_B->second);
    station = &*pps_stationB;
  }
  if (det1->ppsNHits_ <= i || det2->ppsNHits_ <= j)
    return false;
  //
  double x1 = det1->ppsX_.at(i);
  double y1 = det1->ppsY_.at(i);
  double x2 = det2->ppsX_.at(j);
  double y2 = det2->ppsY_.at(j);
  double eloss;
 
  //thx and thy are returned in microrad
  ReconstructArm(
      station, Direction * x1, y1, Direction * x2, y2, thx, thy, eloss);  // Pass the hits in the LHC ref. frame
  thx *= -Direction;                                                      //  invert to the CMS ref frame
 
  // Protect for unphysical results
  if (edm::isNotFinite(eloss) || edm::isNotFinite(thx) || edm::isNotFinite(thy))
    return false;
  //
 
  if (m_verbosity)
    LogDebug("CTPPSFastTrackingProducer::SearchTrack:") << "thx " << thx << " thy " << thy << " eloss " << eloss;
 
  // Get the start point of the reconstructed track near the origin made by Hector in the CMS ref. frame
  x0 = -Direction * station->getX0() * um_to_cm;
  y0 = station->getY0() * um_to_cm;
  double ImpPar = sqrt(x0 * x0 + y0 * y0);
  if (ImpPar > fImpParcut)
    return false;
  if (eloss < 0. || eloss > fBeamEnergy)
    return false;
  //
  // Calculate the reconstructed track parameters
  double theta = sqrt(thx * thx + thy * thy) * urad;
  xi = eloss / fBeamEnergy;
  double energy = fBeamEnergy * (1. - xi);
  partP = sqrt(energy * energy - PPSTools::ProtonMassSQ);
  t = -2. * (PPSTools::ProtonMassSQ - fBeamEnergy * energy + fBeamMomentum * partP * cos(theta));
  pt = partP * theta;
  if (xi < 0. || xi > 1. || t < 0. || t > 10. || pt <= 0.) {
    xi = 0.;
    t = 0.;
    partP = 0.;
    pt = 0.;
    x0 = 0.;
    y0 = 0.;
    return false;  // unphysical values
  }
  //Try to include the timing detector in the track
  ProjectToToF(x1, y1, x2, y2, xt, yt);  // the projections is done in the CMS ref frame
  X1d = x1;
  Y1d = y1;
  X2d = x2;
  Y2d = y2;
  return true;
}  //end  SearchTrack

References funct::cos(), HCALHighEnergyHPDFilter_cfi::energy, fBeamEnergy, fBeamMomentum, fImpParcut, mps_fire::i, edm::isNotFinite(), dqmiolumiharvest::j, LogDebug, m_verbosity, pps_stationB, pps_stationF, CTPPSTrkDetector::ppsNHits_, CTPPSTrkDetector::ppsX_, CTPPSTrkDetector::ppsY_, ProjectToToF(), PPSTools::ProtonMassSQ, DiDispStaMuonMonitor_cfi::pt, ReconstructArm(), mathSSE::sqrt(), relativeConstraints::station, submitPVValidationJobs::t, theta(), TrkStation_B, TrkStation_F, um_to_cm, urad, testProducerWithPsetDescEmpty_cfi::x1, testProducerWithPsetDescEmpty_cfi::x2, protons_cff::xi, testProducerWithPsetDescEmpty_cfi::y1, and testProducerWithPsetDescEmpty_cfi::y2.

Referenced by FastReco().

SetBeamLine()

bool CTPPSFastTrackingProducer::SetBeamLine ( )

private

Definition at line 512 of file CTPPSFastTrackingProducer.cc.

                                            {
  edm::FileInPath b1(beam1filename.c_str());
  edm::FileInPath b2(beam2filename.c_str());
  if (lengthctpps <= 0)
    return false;
  m_beamlineCTPPS1 = std::make_unique<H_BeamLine>(-1, lengthctpps + 0.1);  // (direction, length)
  m_beamlineCTPPS1->fill(b2.fullPath(), 1, "IP5");
  m_beamlineCTPPS2 = std::make_unique<H_BeamLine>(1, lengthctpps + 0.1);  //
  m_beamlineCTPPS2->fill(b1.fullPath(), 1, "IP5");
  m_beamlineCTPPS1->offsetElements(120, 0.097);
  m_beamlineCTPPS2->offsetElements(120, -0.097);
  pps_stationF = std::make_unique<H_RecRPObject>(fz_tracker1, fz_tracker2, *m_beamlineCTPPS1);
  pps_stationB = std::make_unique<H_RecRPObject>(fz_tracker1, fz_tracker2, *m_beamlineCTPPS2);
  return true;
}

References b1, b2, beam1filename, beam2filename, fz_tracker1, fz_tracker2, lengthctpps, m_beamlineCTPPS1, m_beamlineCTPPS2, pps_stationB, and pps_stationF.

Referenced by CTPPSFastTrackingProducer().

TrackerStationClear()

void CTPPSFastTrackingProducer::TrackerStationClear ( )

private

Definition at line 239 of file CTPPSFastTrackingProducer.cc.

                                                    {
  TrkStation_F->first.clear();
  TrkStation_F->second.clear();
  TrkStation_B->first.clear();
  TrkStation_B->second.clear();
}

References TrkStation_B, and TrkStation_F.

Referenced by produce().

TrackerStationStarting()

void CTPPSFastTrackingProducer::TrackerStationStarting ( )

private

Definition at line 246 of file CTPPSFastTrackingProducer.cc.

                                                       {
  det1F->clear();
  det1B->clear();
  det2F->clear();
  det2B->clear();
  detToF_F->clear();
  detToF_B->clear();
}

References det1B, det1F, det2B, det2F, detToF_B, and detToF_F.

Referenced by produce().

Member Data Documentation

_recHitToken

edm::EDGetTokenT<CTPPSFastRecHitContainer> CTPPSFastTrackingProducer::_recHitToken

private

Definition at line 60 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and produce().

beam1filename

std::string CTPPSFastTrackingProducer::beam1filename

private

Definition at line 98 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and SetBeamLine().

beam2filename

std::string CTPPSFastTrackingProducer::beam2filename

private

Definition at line 99 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and SetBeamLine().

det1B

std::unique_ptr<CTPPSTrkDetector> CTPPSFastTrackingProducer::det1B

private

Definition at line 113 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), ReadRecHits(), and TrackerStationStarting().

det1F

std::unique_ptr<CTPPSTrkDetector> CTPPSFastTrackingProducer::det1F

private

Definition at line 112 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), ReadRecHits(), and TrackerStationStarting().

det2B

std::unique_ptr<CTPPSTrkDetector> CTPPSFastTrackingProducer::det2B

private

Definition at line 115 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), ReadRecHits(), and TrackerStationStarting().

det2F

std::unique_ptr<CTPPSTrkDetector> CTPPSFastTrackingProducer::det2F

private

Definition at line 114 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), ReadRecHits(), and TrackerStationStarting().

detToF_B

std::unique_ptr<CTPPSToFDetector> CTPPSFastTrackingProducer::detToF_B

private

Definition at line 117 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), ReadRecHits(), and TrackerStationStarting().

detToF_F

std::unique_ptr<CTPPSToFDetector> CTPPSFastTrackingProducer::detToF_F

private

Definition at line 116 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), ReadRecHits(), and TrackerStationStarting().

fBeamEnergy

double CTPPSFastTrackingProducer::fBeamEnergy

private

Definition at line 104 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), FastReco(), and SearchTrack().

fBeamMomentum

double CTPPSFastTrackingProducer::fBeamMomentum

private

Definition at line 105 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), FastReco(), and SearchTrack().

fBeamXRMS_ToF

double CTPPSFastTrackingProducer::fBeamXRMS_ToF

private

Definition at line 131 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fBeamXRMS_Trk1

double CTPPSFastTrackingProducer::fBeamXRMS_Trk1

private

Definition at line 127 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer().

fBeamXRMS_Trk2

double CTPPSFastTrackingProducer::fBeamXRMS_Trk2

private

Definition at line 127 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer().

fCrossAngleCorr

bool CTPPSFastTrackingProducer::fCrossAngleCorr

private

Definition at line 106 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fCrossingAngleBeam1

double CTPPSFastTrackingProducer::fCrossingAngleBeam1

private

Definition at line 107 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fCrossingAngleBeam2

double CTPPSFastTrackingProducer::fCrossingAngleBeam2

private

Definition at line 108 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fImpParcut

double CTPPSFastTrackingProducer::fImpParcut

private

Definition at line 131 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and SearchTrack().

fTimeSigma

double CTPPSFastTrackingProducer::fTimeSigma

private

Definition at line 131 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fToFCellHeight

double CTPPSFastTrackingProducer::fToFCellHeight

private

Definition at line 129 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fToFCellWidth

std::vector<double> CTPPSFastTrackingProducer::fToFCellWidth

private

Definition at line 128 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fToFInsertion

double CTPPSFastTrackingProducer::fToFInsertion

private

Definition at line 131 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fToFNCellX

int CTPPSFastTrackingProducer::fToFNCellX

private

Definition at line 130 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fToFNCellY

int CTPPSFastTrackingProducer::fToFNCellY

private

Definition at line 130 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fToFPitchX

double CTPPSFastTrackingProducer::fToFPitchX

private

Definition at line 129 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fToFPitchY

double CTPPSFastTrackingProducer::fToFPitchY

private

Definition at line 129 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fToFXOffset

double CTPPSFastTrackingProducer::fToFXOffset

private

Definition at line 131 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fTrackerHeight

double CTPPSFastTrackingProducer::fTrackerHeight

private

Definition at line 127 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer().

fTrackerInsertion

double CTPPSFastTrackingProducer::fTrackerInsertion

private

Definition at line 127 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer().

fTrackerWidth

double CTPPSFastTrackingProducer::fTrackerWidth

private

Definition at line 127 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer().

fTrk1XOffset

double CTPPSFastTrackingProducer::fTrk1XOffset

private

Definition at line 127 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer().

fTrk2XOffset

double CTPPSFastTrackingProducer::fTrk2XOffset

private

Definition at line 127 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer().

fz_timing

double CTPPSFastTrackingProducer::fz_timing

private

Definition at line 126 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and ProjectToToF().

fz_tracker1

double CTPPSFastTrackingProducer::fz_tracker1

private

Definition at line 126 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), ProjectToToF(), and SetBeamLine().

fz_tracker2

double CTPPSFastTrackingProducer::fz_tracker2

private

Definition at line 126 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), ProjectToToF(), and SetBeamLine().

lengthctpps

double CTPPSFastTrackingProducer::lengthctpps

private

Definition at line 102 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and SetBeamLine().

m_beamlineCTPPS1

std::unique_ptr<H_BeamLine> CTPPSFastTrackingProducer::m_beamlineCTPPS1

private

Definition at line 93 of file CTPPSFastTrackingProducer.cc.

Referenced by SetBeamLine().

m_beamlineCTPPS2

std::unique_ptr<H_BeamLine> CTPPSFastTrackingProducer::m_beamlineCTPPS2

private

Definition at line 94 of file CTPPSFastTrackingProducer.cc.

Referenced by SetBeamLine().

m_verbosity

bool CTPPSFastTrackingProducer::m_verbosity

private

Definition at line 103 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and SearchTrack().

pps_stationB

std::unique_ptr<H_RecRPObject> CTPPSFastTrackingProducer::pps_stationB

private

Definition at line 96 of file CTPPSFastTrackingProducer.cc.

Referenced by Reconstruction(), SearchTrack(), and SetBeamLine().

pps_stationF

std::unique_ptr<H_RecRPObject> CTPPSFastTrackingProducer::pps_stationF

private

Definition at line 95 of file CTPPSFastTrackingProducer.cc.

Referenced by Reconstruction(), SearchTrack(), and SetBeamLine().

recCellId_B

std::vector<int> CTPPSFastTrackingProducer::recCellId_B

private

Definition at line 123 of file CTPPSFastTrackingProducer.cc.

Referenced by FastReco(), produce(), and ReadRecHits().

recCellId_F

std::vector<int> CTPPSFastTrackingProducer::recCellId_F

private

Definition at line 123 of file CTPPSFastTrackingProducer.cc.

Referenced by FastReco(), produce(), and ReadRecHits().

recTof_B

std::vector<double> CTPPSFastTrackingProducer::recTof_B

private

Definition at line 124 of file CTPPSFastTrackingProducer.cc.

Referenced by FastReco(), produce(), and ReadRecHits().

recTof_F

std::vector<double> CTPPSFastTrackingProducer::recTof_F

private

Definition at line 124 of file CTPPSFastTrackingProducer.cc.

Referenced by FastReco(), produce(), and ReadRecHits().

theCTPPSFastTrack

std::vector<CTPPSFastTrack> CTPPSFastTrackingProducer::theCTPPSFastTrack

private

Definition at line 119 of file CTPPSFastTrackingProducer.cc.

Referenced by FastReco(), produce(), and Reconstruction().

track

CTPPSFastTrack CTPPSFastTrackingProducer::track

private

Definition at line 121 of file CTPPSFastTrackingProducer.cc.

Referenced by FastReco().

TrkStation_B

std::unique_ptr<CTPPSTrkStation> CTPPSFastTrackingProducer::TrkStation_B

private

Definition at line 111 of file CTPPSFastTrackingProducer.cc.

Referenced by FastReco(), ReadRecHits(), SearchTrack(), and TrackerStationClear().

TrkStation_F

std::unique_ptr<CTPPSTrkStation> CTPPSFastTrackingProducer::TrkStation_F

private

Definition at line 110 of file CTPPSFastTrackingProducer.cc.

Referenced by FastReco(), ReadRecHits(), SearchTrack(), and TrackerStationClear().