CTPPSFastTrackingProducer Class Reference

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

Inheritance diagram for CTPPSFastTrackingProducer:

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

Private Types
typedef std::vector< CTPPSFastRecHit >	CTPPSFastRecHitContainer

Private Member Functions
virtual void	beginEvent (edm::Event &event, const edm::EventSetup &eventSetup)
void	beginStream (edm::StreamID) override
virtual void	endEvent (edm::Event &event, const edm::EventSetup &eventSetup)
void	endStream () override
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
std::map< unsigned int, H_BeamParticle * >	m_beamPart
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<>
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

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 53 of file CTPPSFastTrackingProducer.h.

Member Typedef Documentation

CTPPSFastRecHitContainer

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

private

Definition at line 68 of file CTPPSFastTrackingProducer.h.

Constructor & Destructor Documentation

CTPPSFastTrackingProducer()

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

explicit

Definition at line 34 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::unique_ptr<CTPPSTrkDetector>(
      new CTPPSTrkDetector(fTrackerWidth, fTrackerHeight, fTrackerInsertion * fBeamXRMS_Trk1 + fTrk1XOffset));
  det2F = std::unique_ptr<CTPPSTrkDetector>(
      new CTPPSTrkDetector(fTrackerWidth, fTrackerHeight, fTrackerInsertion * fBeamXRMS_Trk2 + fTrk2XOffset));
  det1B = std::unique_ptr<CTPPSTrkDetector>(
      new CTPPSTrkDetector(fTrackerWidth, fTrackerHeight, fTrackerInsertion * fBeamXRMS_Trk1 + fTrk1XOffset));
  det2B = std::unique_ptr<CTPPSTrkDetector>(
      new CTPPSTrkDetector(fTrackerWidth, fTrackerHeight, fTrackerInsertion * fBeamXRMS_Trk2 + fTrk2XOffset));
 
  //Timing Detector Description
  std::vector<double> vToFCellWidth;
  for (int i = 0; i < 8; i++) {
    vToFCellWidth.push_back(fToFCellWidth[i]);
  }
  double pos_tof = fToFInsertion * fBeamXRMS_ToF + fToFXOffset;
  detToF_F = std::unique_ptr<CTPPSToFDetector>(new CTPPSToFDetector(
      fToFNCellX, fToFNCellY, vToFCellWidth, fToFCellHeight, fToFPitchX, fToFPitchY, pos_tof, fTimeSigma));
  detToF_B = std::unique_ptr<CTPPSToFDetector>(new 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().

~CTPPSFastTrackingProducer()

CTPPSFastTrackingProducer::~CTPPSFastTrackingProducer ( )

override

Definition at line 110 of file CTPPSFastTrackingProducer.cc.

                                                      {
  for (std::map<unsigned int, H_BeamParticle*>::iterator it = m_beamPart.begin(); it != m_beamPart.end(); ++it) {
    delete (*it).second;
  }
}

References m_beamPart.

Member Function Documentation

beginEvent()

void CTPPSFastTrackingProducer::beginEvent ( edm::Event &  event, const edm::EventSetup &  eventSetup  )

privatevirtual

Definition at line 138 of file CTPPSFastTrackingProducer.cc.

                                                                                           {
  TrackerStationStarting();
}

References TrackerStationStarting().

beginStream()

void CTPPSFastTrackingProducer::beginStream ( edm::StreamID  )

overrideprivate

Definition at line 424 of file CTPPSFastTrackingProducer.cc.

{}

endEvent()

void CTPPSFastTrackingProducer::endEvent ( edm::Event &  event, const edm::EventSetup &  eventSetup  )

privatevirtual

Definition at line 142 of file CTPPSFastTrackingProducer.cc.

                                                                                         {
  TrackerStationClear();
}

References TrackerStationClear().

endStream()

void CTPPSFastTrackingProducer::endStream ( )

overrideprivate

Definition at line 428 of file CTPPSFastTrackingProducer.cc.

{}

FastReco()

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

private

Definition at line 346 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;
  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(), OrderedSet::t, theCTPPSFastTrack, theta(), track, TrkStation_B, TrkStation_F, urad, bphysicsOniaDQM_cfi::vertex, and hybridSuperClusters_cfi::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 335 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 116 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 95 of file CTPPSFastTrackingProducer.h.

                                                                                                                {
    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 164 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::unique_ptr<CTPPSTrkStation>(new std::pair<CTPPSTrkDetector, CTPPSTrkDetector>(*det1F, *det2F));
  TrkStation_B = std::unique_ptr<CTPPSTrkStation>(new std::pair<CTPPSTrkDetector, CTPPSTrkDetector>(*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 315 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 (std::isnan(energy) || std::isinf(energy))
    return;
  tx = pps_station->getTXIP();  // change orientation to CMS
  ty = pps_station->getTYIP();
  eloss = pps_station->getE();
}

References HCALHighEnergyHPDFilter_cfi::energy, edm::detail::isnan(), 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 204 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 213 of file CTPPSFastTrackingProducer.cc.

                                                         {
  // Given 1 hit in Tracker1 and 1 hit in Tracker2 try to make a track with Hector
  double theta = 0.;
  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 (std::isnan(eloss) || std::isinf(eloss) || std::isnan(thx) || std::isinf(thx) || std::isnan(thy) ||
      std::isinf(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
  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 = sqrt(pow(partP * thx * urad, 2) + pow(partP * thy * urad, 2));
  if (xi < 0. || xi > 1. || t < 0. || t > 10. || pt <= 0.) {
    xi = 0.;
    t = 0.;
    partP = 0.;
    pt = 0.;
    theta = 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::detail::isnan(), dqmiolumiharvest::j, LogDebug, m_verbosity, funct::pow(), pps_stationB, pps_stationF, CTPPSTrkDetector::ppsNHits_, CTPPSTrkDetector::ppsX_, CTPPSTrkDetector::ppsY_, ProjectToToF(), PPSTools::ProtonMassSQ, DiDispStaMuonMonitor_cfi::pt, ReconstructArm(), mathSSE::sqrt(), relativeConstraints::station, OrderedSet::t, theta(), TrkStation_B, TrkStation_F, um_to_cm, urad, testProducerWithPsetDescEmpty_cfi::x1, testProducerWithPsetDescEmpty_cfi::x2, hybridSuperClusters_cfi::xi, testProducerWithPsetDescEmpty_cfi::y1, and testProducerWithPsetDescEmpty_cfi::y2.

Referenced by FastReco().

SetBeamLine()

bool CTPPSFastTrackingProducer::SetBeamLine ( )

private

Definition at line 429 of file CTPPSFastTrackingProducer.cc.

                                            {
  edm::FileInPath b1(beam1filename.c_str());
  edm::FileInPath b2(beam2filename.c_str());
  if (lengthctpps <= 0)
    return false;
  m_beamlineCTPPS1 = std::unique_ptr<H_BeamLine>(new H_BeamLine(-1, lengthctpps + 0.1));  // (direction, length)
  m_beamlineCTPPS1->fill(b2.fullPath(), 1, "IP5");
  m_beamlineCTPPS2 = std::unique_ptr<H_BeamLine>(new 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::unique_ptr<H_RecRPObject>(new H_RecRPObject(fz_tracker1, fz_tracker2, *m_beamlineCTPPS1));
  pps_stationB = std::unique_ptr<H_RecRPObject>(new 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 147 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 endEvent(), and produce().

TrackerStationStarting()

void CTPPSFastTrackingProducer::TrackerStationStarting ( )

private

Definition at line 154 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 beginEvent(), and produce().

Member Data Documentation

_recHitToken

edm::EDGetTokenT<CTPPSFastRecHitContainer> CTPPSFastTrackingProducer::_recHitToken

private

Definition at line 69 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and produce().

beam1filename

std::string CTPPSFastTrackingProducer::beam1filename

private

Definition at line 108 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and SetBeamLine().

beam2filename

std::string CTPPSFastTrackingProducer::beam2filename

private

Definition at line 109 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and SetBeamLine().

det1B

std::unique_ptr<CTPPSTrkDetector> CTPPSFastTrackingProducer::det1B

private

Definition at line 123 of file CTPPSFastTrackingProducer.h.

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

det1F

std::unique_ptr<CTPPSTrkDetector> CTPPSFastTrackingProducer::det1F

private

Definition at line 122 of file CTPPSFastTrackingProducer.h.

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

det2B

std::unique_ptr<CTPPSTrkDetector> CTPPSFastTrackingProducer::det2B

private

Definition at line 125 of file CTPPSFastTrackingProducer.h.

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

det2F

std::unique_ptr<CTPPSTrkDetector> CTPPSFastTrackingProducer::det2F

private

Definition at line 124 of file CTPPSFastTrackingProducer.h.

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

detToF_B

std::unique_ptr<CTPPSToFDetector> CTPPSFastTrackingProducer::detToF_B

private

Definition at line 127 of file CTPPSFastTrackingProducer.h.

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

detToF_F

std::unique_ptr<CTPPSToFDetector> CTPPSFastTrackingProducer::detToF_F

private

Definition at line 126 of file CTPPSFastTrackingProducer.h.

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

fBeamEnergy

double CTPPSFastTrackingProducer::fBeamEnergy

private

Definition at line 114 of file CTPPSFastTrackingProducer.h.

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

fBeamMomentum

double CTPPSFastTrackingProducer::fBeamMomentum

private

Definition at line 115 of file CTPPSFastTrackingProducer.h.

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

fBeamXRMS_ToF

double CTPPSFastTrackingProducer::fBeamXRMS_ToF

private

Definition at line 141 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fBeamXRMS_Trk1

double CTPPSFastTrackingProducer::fBeamXRMS_Trk1

private

Definition at line 137 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer().

fBeamXRMS_Trk2

double CTPPSFastTrackingProducer::fBeamXRMS_Trk2

private

Definition at line 137 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer().

fCrossAngleCorr

bool CTPPSFastTrackingProducer::fCrossAngleCorr

private

Definition at line 116 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fCrossingAngleBeam1

double CTPPSFastTrackingProducer::fCrossingAngleBeam1

private

Definition at line 117 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fCrossingAngleBeam2

double CTPPSFastTrackingProducer::fCrossingAngleBeam2

private

Definition at line 118 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fImpParcut

double CTPPSFastTrackingProducer::fImpParcut

private

Definition at line 141 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and SearchTrack().

fTimeSigma

double CTPPSFastTrackingProducer::fTimeSigma

private

Definition at line 141 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fToFCellHeight

double CTPPSFastTrackingProducer::fToFCellHeight

private

Definition at line 139 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fToFCellWidth

std::vector<double> CTPPSFastTrackingProducer::fToFCellWidth

private

Definition at line 138 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fToFInsertion

double CTPPSFastTrackingProducer::fToFInsertion

private

Definition at line 141 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fToFNCellX

int CTPPSFastTrackingProducer::fToFNCellX

private

Definition at line 140 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fToFNCellY

int CTPPSFastTrackingProducer::fToFNCellY

private

Definition at line 140 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fToFPitchX

double CTPPSFastTrackingProducer::fToFPitchX

private

Definition at line 139 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fToFPitchY

double CTPPSFastTrackingProducer::fToFPitchY

private

Definition at line 139 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fToFXOffset

double CTPPSFastTrackingProducer::fToFXOffset

private

Definition at line 141 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

fTrackerHeight

double CTPPSFastTrackingProducer::fTrackerHeight

private

Definition at line 137 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer().

fTrackerInsertion

double CTPPSFastTrackingProducer::fTrackerInsertion

private

Definition at line 137 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer().

fTrackerWidth

double CTPPSFastTrackingProducer::fTrackerWidth

private

Definition at line 137 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer().

fTrk1XOffset

double CTPPSFastTrackingProducer::fTrk1XOffset

private

Definition at line 137 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer().

fTrk2XOffset

double CTPPSFastTrackingProducer::fTrk2XOffset

private

Definition at line 137 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer().

fz_timing

double CTPPSFastTrackingProducer::fz_timing

private

Definition at line 136 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and ProjectToToF().

fz_tracker1

double CTPPSFastTrackingProducer::fz_tracker1

private

Definition at line 136 of file CTPPSFastTrackingProducer.h.

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

fz_tracker2

double CTPPSFastTrackingProducer::fz_tracker2

private

Definition at line 136 of file CTPPSFastTrackingProducer.h.

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

lengthctpps

double CTPPSFastTrackingProducer::lengthctpps

private

Definition at line 112 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and SetBeamLine().

m_beamlineCTPPS1

std::unique_ptr<H_BeamLine> CTPPSFastTrackingProducer::m_beamlineCTPPS1

private

Definition at line 103 of file CTPPSFastTrackingProducer.h.

Referenced by SetBeamLine().

m_beamlineCTPPS2

std::unique_ptr<H_BeamLine> CTPPSFastTrackingProducer::m_beamlineCTPPS2

private

Definition at line 104 of file CTPPSFastTrackingProducer.h.

Referenced by SetBeamLine().

m_beamPart

std::map<unsigned int, H_BeamParticle*> CTPPSFastTrackingProducer::m_beamPart

private

Definition at line 102 of file CTPPSFastTrackingProducer.h.

Referenced by ~CTPPSFastTrackingProducer().

m_verbosity

bool CTPPSFastTrackingProducer::m_verbosity

private

Definition at line 113 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and SearchTrack().

pps_stationB

std::unique_ptr<H_RecRPObject> CTPPSFastTrackingProducer::pps_stationB

private

Definition at line 106 of file CTPPSFastTrackingProducer.h.

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

pps_stationF

std::unique_ptr<H_RecRPObject> CTPPSFastTrackingProducer::pps_stationF

private

Definition at line 105 of file CTPPSFastTrackingProducer.h.

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

recCellId_B

std::vector<int> CTPPSFastTrackingProducer::recCellId_B

private

Definition at line 133 of file CTPPSFastTrackingProducer.h.

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

recCellId_F

std::vector<int> CTPPSFastTrackingProducer::recCellId_F

private

Definition at line 133 of file CTPPSFastTrackingProducer.h.

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

recTof_B

std::vector<double> CTPPSFastTrackingProducer::recTof_B

private

Definition at line 134 of file CTPPSFastTrackingProducer.h.

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

recTof_F

std::vector<double> CTPPSFastTrackingProducer::recTof_F

private

Definition at line 134 of file CTPPSFastTrackingProducer.h.

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

theCTPPSFastTrack

std::vector<CTPPSFastTrack> CTPPSFastTrackingProducer::theCTPPSFastTrack

private

Definition at line 129 of file CTPPSFastTrackingProducer.h.

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

track

CTPPSFastTrack CTPPSFastTrackingProducer::track

private

Definition at line 131 of file CTPPSFastTrackingProducer.h.

Referenced by FastReco().

TrkStation_B

std::unique_ptr<CTPPSTrkStation> CTPPSFastTrackingProducer::TrkStation_B

private

Definition at line 121 of file CTPPSFastTrackingProducer.h.

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

TrkStation_F

std::unique_ptr<CTPPSTrkStation> CTPPSFastTrackingProducer::TrkStation_F

private

Definition at line 120 of file CTPPSFastTrackingProducer.h.

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