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

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

private

Definition at line 59 of file CTPPSFastTrackingProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 136 of file CTPPSFastTrackingProducer.cc.

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().

    : 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);
  //
}

Member Function Documentation

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

private

Definition at line 435 of file CTPPSFastTrackingProducer.cc.

References funct::cos(), alignCSCRings::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, dqmiolumiharvest::j, PPSTools::LorentzBoost(), MatchCellId(), AlCaHLTBitMon_ParallelJobs::p, phi, Pi, CTPPSTrkDetector::ppsNHits_, PPSTools::ProtonMassSQ, DiDispStaMuonMonitor_cfi::pt, 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, and urad.

Referenced by Reconstruction().

                                                                              {
  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

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.

References mps_fire::i.

Referenced by FastReco().

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

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

overrideprivate

Definition at line 215 of file CTPPSFastTrackingProducer.cc.

References _recHitToken, edm::Event::getByToken(), mps_fire::i, eostools::move(), dqmiodumpmetadata::n, edm::Event::put(), ReadRecHits(), recCellId_B, recCellId_F, HLT_FULL_cff::recHits, Reconstruction(), recTof_B, recTof_F, theCTPPSFastTrack, TrackerStationClear(), and TrackerStationStarting().

                                                                                     {
  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

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.

References fz_timing, fz_tracker1, and fz_tracker2.

Referenced by SearchTrack().

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

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

private

Definition at line 256 of file CTPPSFastTrackingProducer.cc.

References det1B, det1F, det2B, det2F, detToF_B, detToF_F, CTPPSFastRecHit::detUnitId(), CTPPSFastRecHit::entryPoint(), recCellId_B, recCellId_F, recHitDet(), recTof_B, recTof_F, CTPPSFastRecHit::tof(), TrkStation_B, TrkStation_F, x, PV3DBase< T, PVType, FrameType >::x(), y, PV3DBase< T, PVType, FrameType >::y(), z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by produce().

                                                                                        {
  // 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

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.

References relval_parameters_module::energy, edm::isNotFinite(), and mm_to_um.

Referenced by SearchTrack().

                                                                                                                   {
  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();
}

void CTPPSFastTrackingProducer::Reconstruction ( )

private

Definition at line 296 of file CTPPSFastTrackingProducer.cc.

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

Referenced by produce().

                                               {
  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

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.

References funct::cos(), relval_parameters_module::energy, fBeamEnergy, fBeamMomentum, fImpParcut, edm::isNotFinite(), LogDebug, m_verbosity, pps_stationB, pps_stationF, CTPPSTrkDetector::ppsNHits_, CTPPSTrkDetector::ppsX_, CTPPSTrkDetector::ppsY_, ProjectToToF(), PPSTools::ProtonMassSQ, ReconstructArm(), mathSSE::sqrt(), relativeConstraints::station, theta(), TrkStation_B, TrkStation_F, um_to_cm, and urad.

Referenced by FastReco().

                                                         {
  // 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

bool CTPPSFastTrackingProducer::SetBeamLine ( )

private

Definition at line 512 of file CTPPSFastTrackingProducer.cc.

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

Referenced by CTPPSFastTrackingProducer().

                                            {
  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;
}

void CTPPSFastTrackingProducer::TrackerStationClear ( )

private

Definition at line 239 of file CTPPSFastTrackingProducer.cc.

References TrkStation_B, and TrkStation_F.

Referenced by produce().

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

void CTPPSFastTrackingProducer::TrackerStationStarting ( )

private

Definition at line 246 of file CTPPSFastTrackingProducer.cc.

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

Referenced by produce().

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

Member Data Documentation

edm::EDGetTokenT<CTPPSFastRecHitContainer> CTPPSFastTrackingProducer::_recHitToken

private

Definition at line 60 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and produce().

std::string CTPPSFastTrackingProducer::beam1filename

private

Definition at line 98 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and SetBeamLine().

std::string CTPPSFastTrackingProducer::beam2filename

private

Definition at line 99 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and SetBeamLine().

std::unique_ptr<CTPPSTrkDetector> CTPPSFastTrackingProducer::det1B

private

Definition at line 113 of file CTPPSFastTrackingProducer.cc.

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

std::unique_ptr<CTPPSTrkDetector> CTPPSFastTrackingProducer::det1F

private

Definition at line 112 of file CTPPSFastTrackingProducer.cc.

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

std::unique_ptr<CTPPSTrkDetector> CTPPSFastTrackingProducer::det2B

private

Definition at line 115 of file CTPPSFastTrackingProducer.cc.

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

std::unique_ptr<CTPPSTrkDetector> CTPPSFastTrackingProducer::det2F

private

Definition at line 114 of file CTPPSFastTrackingProducer.cc.

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

std::unique_ptr<CTPPSToFDetector> CTPPSFastTrackingProducer::detToF_B

private

Definition at line 117 of file CTPPSFastTrackingProducer.cc.

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

std::unique_ptr<CTPPSToFDetector> CTPPSFastTrackingProducer::detToF_F

private

Definition at line 116 of file CTPPSFastTrackingProducer.cc.

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

double CTPPSFastTrackingProducer::fBeamEnergy

private

Definition at line 104 of file CTPPSFastTrackingProducer.cc.

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

double CTPPSFastTrackingProducer::fBeamMomentum

private

Definition at line 105 of file CTPPSFastTrackingProducer.cc.

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

double CTPPSFastTrackingProducer::fBeamXRMS_ToF

private

Definition at line 131 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

double CTPPSFastTrackingProducer::fBeamXRMS_Trk1

private

Definition at line 127 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer().

double CTPPSFastTrackingProducer::fBeamXRMS_Trk2

private

Definition at line 127 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer().

bool CTPPSFastTrackingProducer::fCrossAngleCorr

private

Definition at line 106 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

double CTPPSFastTrackingProducer::fCrossingAngleBeam1

private

Definition at line 107 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

double CTPPSFastTrackingProducer::fCrossingAngleBeam2

private

Definition at line 108 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

double CTPPSFastTrackingProducer::fImpParcut

private

Definition at line 131 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and SearchTrack().

double CTPPSFastTrackingProducer::fTimeSigma

private

Definition at line 131 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

double CTPPSFastTrackingProducer::fToFCellHeight

private

Definition at line 129 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

std::vector<double> CTPPSFastTrackingProducer::fToFCellWidth

private

Definition at line 128 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

double CTPPSFastTrackingProducer::fToFInsertion

private

Definition at line 131 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

int CTPPSFastTrackingProducer::fToFNCellX

private

Definition at line 130 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

int CTPPSFastTrackingProducer::fToFNCellY

private

Definition at line 130 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

double CTPPSFastTrackingProducer::fToFPitchX

private

Definition at line 129 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

double CTPPSFastTrackingProducer::fToFPitchY

private

Definition at line 129 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

double CTPPSFastTrackingProducer::fToFXOffset

private

Definition at line 131 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

double CTPPSFastTrackingProducer::fTrackerHeight

private

Definition at line 127 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer().

double CTPPSFastTrackingProducer::fTrackerInsertion

private

Definition at line 127 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer().

double CTPPSFastTrackingProducer::fTrackerWidth

private

Definition at line 127 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer().

double CTPPSFastTrackingProducer::fTrk1XOffset

private

Definition at line 127 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer().

double CTPPSFastTrackingProducer::fTrk2XOffset

private

Definition at line 127 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer().

double CTPPSFastTrackingProducer::fz_timing

private

Definition at line 126 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and ProjectToToF().

double CTPPSFastTrackingProducer::fz_tracker1

private

Definition at line 126 of file CTPPSFastTrackingProducer.cc.

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

double CTPPSFastTrackingProducer::fz_tracker2

private

Definition at line 126 of file CTPPSFastTrackingProducer.cc.

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

double CTPPSFastTrackingProducer::lengthctpps

private

Definition at line 102 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and SetBeamLine().

std::unique_ptr<H_BeamLine> CTPPSFastTrackingProducer::m_beamlineCTPPS1

private

Definition at line 93 of file CTPPSFastTrackingProducer.cc.

Referenced by SetBeamLine().

std::unique_ptr<H_BeamLine> CTPPSFastTrackingProducer::m_beamlineCTPPS2

private

Definition at line 94 of file CTPPSFastTrackingProducer.cc.

Referenced by SetBeamLine().

bool CTPPSFastTrackingProducer::m_verbosity

private

Definition at line 103 of file CTPPSFastTrackingProducer.cc.

Referenced by CTPPSFastTrackingProducer(), and SearchTrack().

std::unique_ptr<H_RecRPObject> CTPPSFastTrackingProducer::pps_stationB

private

Definition at line 96 of file CTPPSFastTrackingProducer.cc.

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

std::unique_ptr<H_RecRPObject> CTPPSFastTrackingProducer::pps_stationF

private

Definition at line 95 of file CTPPSFastTrackingProducer.cc.

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

std::vector<int> CTPPSFastTrackingProducer::recCellId_B

private

Definition at line 123 of file CTPPSFastTrackingProducer.cc.

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

std::vector<int> CTPPSFastTrackingProducer::recCellId_F

private

Definition at line 123 of file CTPPSFastTrackingProducer.cc.

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

std::vector<double> CTPPSFastTrackingProducer::recTof_B

private

Definition at line 124 of file CTPPSFastTrackingProducer.cc.

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

std::vector<double> CTPPSFastTrackingProducer::recTof_F

private

Definition at line 124 of file CTPPSFastTrackingProducer.cc.

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

std::vector<CTPPSFastTrack> CTPPSFastTrackingProducer::theCTPPSFastTrack

private

Definition at line 119 of file CTPPSFastTrackingProducer.cc.

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

CTPPSFastTrack CTPPSFastTrackingProducer::track

private

Definition at line 121 of file CTPPSFastTrackingProducer.cc.

Referenced by FastReco().

std::unique_ptr<CTPPSTrkStation> CTPPSFastTrackingProducer::TrkStation_B

private

Definition at line 111 of file CTPPSFastTrackingProducer.cc.

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

std::unique_ptr<CTPPSTrkStation> CTPPSFastTrackingProducer::TrkStation_F

private

Definition at line 110 of file CTPPSFastTrackingProducer.cc.

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