CTPPSFastTrackingProducer Class Reference

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

Inheritance diagram for CTPPSFastTrackingProducer:

Public Types
typedef CLHEP::HepLorentzVector	LorentzVector
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

Public Member Functions
	CTPPSFastTrackingProducer (const edm::ParameterSet &)
	~CTPPSFastTrackingProducer () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () 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	Get_t_and_xi (const LorentzVector *p, double &t, double &xi)
void	LorentzBoost (LorentzVector &p_out, const string &frame)
void	MatchCellId (int cellId, vector< int > vrecCellId, 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)
void	TrackerStationClear ()
void	TrackerStationStarting ()

Private Attributes
edm::EDGetTokenT< CTPPSFastRecHitContainer >	_recHitToken
string	beam1filename
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	fCrossingAngle
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

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 63 of file CTPPSFastTrackingProducer.h.

Member Typedef Documentation

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

private

Definition at line 79 of file CTPPSFastTrackingProducer.h.

typedef CLHEP::HepLorentzVector CTPPSFastTrackingProducer::LorentzVector

Definition at line 67 of file CTPPSFastTrackingProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 30 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, fCrossingAngle, 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_beamlineCTPPS1, m_beamlineCTPPS2, m_verbosity, pps_stationB, pps_stationF, ProtonMassSQ, and mathSSE::sqrt().

                                                                                  :
    m_verbosity(false), fBeamMomentum(0.), fCrossAngleCorr(false), fCrossingAngle(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 - ProtonMassSQ); 
    fCrossAngleCorr = iConfig.getParameter<bool>("CrossAngleCorr");
    fCrossingAngle = iConfig.getParameter<double>("CrossingAngle");
  
    //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");    



    // reading beamlines
    FileInPath b1(beam1filename.c_str());
    FileInPath b2(beam2filename.c_str());
    //
    if(lengthctpps>0. ) {
        m_beamlineCTPPS1 = std::unique_ptr<H_BeamLine>(new H_BeamLine( -1, lengthctpps + 0.1 )); // (direction, length)
        m_beamlineCTPPS2 = std::unique_ptr<H_BeamLine>(new H_BeamLine( 1, lengthctpps + 0.1 )); //
        m_beamlineCTPPS1->fill( b2.fullPath(), 1, "IP5" );
        m_beamlineCTPPS2->fill( b1.fullPath(), 1, "IP5" );
        m_beamlineCTPPS1->offsetElements( 120, 0.097 );
        m_beamlineCTPPS2->offsetElements( 120, 0.097 );
        m_beamlineCTPPS1->calcMatrix();
        m_beamlineCTPPS2->calcMatrix();
    } else {
        if ( m_verbosity ) LogDebug("CTPPSFastTrackingProducer") << "CTPPSFastTrackingProducer: WARNING: lengthctpps=  " << lengthctpps;
    } 

    // Create a particle to get the beam energy from the beam file
    // Take care: the z inside the station is in meters
    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));
    //
    //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)); 

}

CTPPSFastTrackingProducer::~CTPPSFastTrackingProducer ( )

override

Definition at line 114 of file CTPPSFastTrackingProducer.cc.

References m_beamPart.

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

Member Function Documentation

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

privatevirtual

Definition at line 146 of file CTPPSFastTrackingProducer.cc.

References TrackerStationStarting().

{

    TrackerStationStarting();

}

void CTPPSFastTrackingProducer::beginStream ( edm::StreamID  )

overrideprivate

Definition at line 435 of file CTPPSFastTrackingProducer.cc.

{
}

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

privatevirtual

Definition at line 153 of file CTPPSFastTrackingProducer.cc.

References TrackerStationClear().

{
    TrackerStationClear();

}

void CTPPSFastTrackingProducer::endStream ( )

overrideprivate

Definition at line 441 of file CTPPSFastTrackingProducer.cc.

References DEFINE_FWK_MODULE.

{
}

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

private

Definition at line 348 of file CTPPSFastTrackingProducer.cc.

References funct::cos(), MillePedeFileConverter_cfg::e, fBeamXRMS_ToF, fCrossAngleCorr, CTPPSToFDetector::findCellId(), fTimeSigma, fToFCellHeight, fToFCellWidth, fToFInsertion, fToFNCellX, fToFNCellY, fToFPitchX, fToFPitchY, fToFXOffset, Get_t_and_xi(), mps_fire::i, createfilelist::int, LorentzBoost(), MatchCellId(), AlCaHLTBitMon_ParallelJobs::p, phi, pi, CTPPSTrkDetector::ppsNHits_, ProtonMassSQ, EnergyCorrector::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(), lumiQTWidget::t, theCTPPSFastTrack, theta(), track, TrkStation_B, TrkStation_F, urad, and hybridSuperClusters_cfi::xi.

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;
    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 = CLHEP::pi - sqrt(thx*thx+thy*thy)*urad;
                    MatchCellId(cellId, recCellId_B, recTof_B, matchCellId, recTof); 
                }
                phi   = (Direction>0)?-atan2(thy,-thx):atan2(thy,thx); // defined according to the positive direction


                double px = partP*sin(theta)*cos(phi);
                double py = partP*sin(theta)*sin(phi);
                double pz = partP*cos(theta);
                double  e = sqrt(partP*partP+ProtonMassSQ);
                LorentzVector p(px,py,pz,e);
                // Invert the Lorentz boost made to take into account the crossing angle during simulation
                if (fCrossAngleCorr) LorentzBoost(p,"MC");
                //Getting the Xi and t (squared four momentum transferred) of the reconstructed track
                Get_t_and_xi(const_cast<LorentzVector*>(&p),t,xi);
                double pxx = p.px(); double pyy = p.py(); double pzz = p.pz(); //double ee = p.E(); 
                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::Get_t_and_xi ( const LorentzVector *  p, double &  t, double &  xi  )

private

Definition at line 421 of file CTPPSFastTrackingProducer.cc.

References funct::cos(), fBeamEnergy, fBeamMomentum, pi, ProtonMassSQ, mathSSE::sqrt(), and theta().

Referenced by FastReco().

                                                                                             {
    t = 0.;
    xi = -1.;
    if (!proton) return;
    double mom = sqrt((proton->px())*(proton->px())+(proton->py())*(proton->py())+(proton->pz())*(proton->pz()));
    if (mom>fBeamMomentum) mom=fBeamMomentum;
    double energy = proton->e();
    double theta  = (proton->pz()>0)?proton->theta():CLHEP::pi-proton->theta();
    t      = -2.*(ProtonMassSQ-fBeamEnergy*energy+fBeamMomentum*mom*cos(theta));
    xi     = (1.0-energy/fBeamEnergy);
}

void CTPPSFastTrackingProducer::LorentzBoost ( LorentzVector &  p_out, const string &  frame  )

private

Definition at line 299 of file CTPPSFastTrackingProducer.cc.

References funct::cos(), fCrossingAngle, p4, funct::sin(), and funct::tan().

Referenced by FastReco().

{
    // Use a matrix
    double microrad = 1.e-6;
    TMatrixD tmpboost(4,4);
    double alpha_ = 0.;
    double phi_  = fCrossingAngle*microrad;
    if (p_out.pz()<0) phi_*=-1;
    tmpboost(0,0) = 1./cos(phi_);
    tmpboost(0,1) = - cos(alpha_)*sin(phi_);
    tmpboost(0,2) = - tan(phi_)*sin(phi_);
    tmpboost(0,3) = - sin(alpha_)*sin(phi_);
    tmpboost(1,0) = - cos(alpha_)*tan(phi_);
    tmpboost(1,1) = 1.;
    tmpboost(1,2) = cos(alpha_)*tan(phi_);
    tmpboost(1,3) = 0.;
    tmpboost(2,0) = 0.;
    tmpboost(2,1) = - cos(alpha_)*sin(phi_);
    tmpboost(2,2) = cos(phi_);
    tmpboost(2,3) = - sin(alpha_)*sin(phi_);
    tmpboost(3,0) = - sin(alpha_)*tan(phi_);
    tmpboost(3,1) = 0.;
    tmpboost(3,2) = sin(alpha_)*tan(phi_);
    tmpboost(3,3) = 1.;

    if(frame=="LAB") tmpboost.Invert();

    TMatrixD p4(4,1);
    p4(0,0) = p_out.e();
    p4(1,0) = p_out.px();
    p4(2,0) = p_out.py();
    p4(3,0) = p_out.pz();
    TMatrixD p4lab(4,1);
    p4lab = tmpboost * p4;
    p_out.setPx(p4lab(1,0));
    p_out.setPy(p4lab(2,0));
    p_out.setPz(p4lab(3,0));
    p_out.setE(p4lab(0,0));
}

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

private

Definition at line 338 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 122 of file CTPPSFastTrackingProducer.cc.

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

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), and HTMLExport.HTMLExportStatic::export().

{
    using namespace edm;
    //using namespace std;  
    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 91 of file CTPPSFastTrackingProducer.h.

References fz_timing, fz_tracker1, fz_tracker2, and globals_cff::x2.

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 181 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::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

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 280 of file CTPPSFastTrackingProducer.cc.

References edm::detail::isnan(), and mm_to_um.

Referenced by SearchTrack().

{
    tx=0.;
    ty=0.;
    eloss=0.;
    if (!pps_station) return;
    // Change the orientation and units according to Hector and LHC coordinates
    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();
}

void CTPPSFastTrackingProducer::Reconstruction ( )

private

Definition at line 212 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 222 of file CTPPSFastTrackingProducer.cc.

References funct::cos(), gather_cfg::cout, fBeamEnergy, fBeamMomentum, fImpParcut, edm::detail::isnan(), m_verbosity, funct::pow(), pps_stationB, pps_stationF, CTPPSTrkDetector::ppsNHits_, CTPPSTrkDetector::ppsX_, CTPPSTrkDetector::ppsY_, ProjectToToF(), ProtonMassSQ, ReconstructArm(), mathSSE::sqrt(), relativeConstraints::station, theta(), TrkStation_B, TrkStation_F, um_to_cm, urad, globals_cff::x1, globals_cff::x2, and hybridSuperClusters_cfi::xi.

Referenced by FastReco().

{
    // 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,x1,y1,x2,y2,thx,thy,eloss);
    // 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 (-thx<-100||-thx>300) return false;
    if (thy<-200||thy>200) return false;
    //
    if ( m_verbosity ) std::cout << "thx " << thx << " thy " << thy << " eloss " << eloss << std::endl;

    // Get the start point of the reconstructed track near the origin made by Hector
    x0 = -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-ProtonMassSQ);
    t     = -2.*(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);
    X1d = x1;
    Y1d = y1;
    X2d = x2;
    Y2d = y2;
    return true;
}//end  SearchTrack

void CTPPSFastTrackingProducer::TrackerStationClear ( )

private

Definition at line 160 of file CTPPSFastTrackingProducer.cc.

References TrkStation_B, and TrkStation_F.

Referenced by endEvent(), and produce().

                                                   {

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


}

void CTPPSFastTrackingProducer::TrackerStationStarting ( )

private

Definition at line 168 of file CTPPSFastTrackingProducer.cc.

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

Referenced by beginEvent(), and 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 80 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and produce().

string CTPPSFastTrackingProducer::beam1filename

private

Definition at line 104 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer().

string CTPPSFastTrackingProducer::beam2filename

private

Definition at line 105 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer().

std::unique_ptr<CTPPSTrkDetector> CTPPSFastTrackingProducer::det1B

private

Definition at line 118 of file CTPPSFastTrackingProducer.h.

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

std::unique_ptr<CTPPSTrkDetector> CTPPSFastTrackingProducer::det1F

private

Definition at line 117 of file CTPPSFastTrackingProducer.h.

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

std::unique_ptr<CTPPSTrkDetector> CTPPSFastTrackingProducer::det2B

private

Definition at line 120 of file CTPPSFastTrackingProducer.h.

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

std::unique_ptr<CTPPSTrkDetector> CTPPSFastTrackingProducer::det2F

private

Definition at line 119 of file CTPPSFastTrackingProducer.h.

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

std::unique_ptr<CTPPSToFDetector> CTPPSFastTrackingProducer::detToF_B

private

Definition at line 122 of file CTPPSFastTrackingProducer.h.

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

std::unique_ptr<CTPPSToFDetector> CTPPSFastTrackingProducer::detToF_F

private

Definition at line 121 of file CTPPSFastTrackingProducer.h.

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

double CTPPSFastTrackingProducer::fBeamEnergy

private

Definition at line 110 of file CTPPSFastTrackingProducer.h.

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

double CTPPSFastTrackingProducer::fBeamMomentum

private

Definition at line 111 of file CTPPSFastTrackingProducer.h.

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

double CTPPSFastTrackingProducer::fBeamXRMS_ToF

private

Definition at line 136 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

double CTPPSFastTrackingProducer::fBeamXRMS_Trk1

private

Definition at line 132 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer().

double CTPPSFastTrackingProducer::fBeamXRMS_Trk2

private

Definition at line 132 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer().

bool CTPPSFastTrackingProducer::fCrossAngleCorr

private

Definition at line 112 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

double CTPPSFastTrackingProducer::fCrossingAngle

private

Definition at line 113 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and LorentzBoost().

double CTPPSFastTrackingProducer::fImpParcut

private

Definition at line 136 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and SearchTrack().

double CTPPSFastTrackingProducer::fTimeSigma

private

Definition at line 136 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

double CTPPSFastTrackingProducer::fToFCellHeight

private

Definition at line 134 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

std::vector<double> CTPPSFastTrackingProducer::fToFCellWidth

private

Definition at line 133 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

double CTPPSFastTrackingProducer::fToFInsertion

private

Definition at line 136 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

int CTPPSFastTrackingProducer::fToFNCellX

private

Definition at line 135 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

int CTPPSFastTrackingProducer::fToFNCellY

private

Definition at line 135 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

double CTPPSFastTrackingProducer::fToFPitchX

private

Definition at line 134 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

double CTPPSFastTrackingProducer::fToFPitchY

private

Definition at line 134 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

double CTPPSFastTrackingProducer::fToFXOffset

private

Definition at line 136 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and FastReco().

double CTPPSFastTrackingProducer::fTrackerHeight

private

Definition at line 132 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer().

double CTPPSFastTrackingProducer::fTrackerInsertion

private

Definition at line 132 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer().

double CTPPSFastTrackingProducer::fTrackerWidth

private

Definition at line 132 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer().

double CTPPSFastTrackingProducer::fTrk1XOffset

private

Definition at line 132 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer().

double CTPPSFastTrackingProducer::fTrk2XOffset

private

Definition at line 132 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer().

double CTPPSFastTrackingProducer::fz_timing

private

Definition at line 131 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and ProjectToToF().

double CTPPSFastTrackingProducer::fz_tracker1

private

Definition at line 131 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and ProjectToToF().

double CTPPSFastTrackingProducer::fz_tracker2

private

Definition at line 131 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and ProjectToToF().

double CTPPSFastTrackingProducer::lengthctpps

private

Definition at line 108 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer().

std::unique_ptr<H_BeamLine> CTPPSFastTrackingProducer::m_beamlineCTPPS1

private

Definition at line 98 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer().

std::unique_ptr<H_BeamLine> CTPPSFastTrackingProducer::m_beamlineCTPPS2

private

Definition at line 99 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer().

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

private

Definition at line 94 of file CTPPSFastTrackingProducer.h.

Referenced by ~CTPPSFastTrackingProducer().

bool CTPPSFastTrackingProducer::m_verbosity

private

Definition at line 109 of file CTPPSFastTrackingProducer.h.

Referenced by CTPPSFastTrackingProducer(), and SearchTrack().

std::unique_ptr<H_RecRPObject> CTPPSFastTrackingProducer::pps_stationB

private

Definition at line 102 of file CTPPSFastTrackingProducer.h.

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

std::unique_ptr<H_RecRPObject> CTPPSFastTrackingProducer::pps_stationF

private

Definition at line 101 of file CTPPSFastTrackingProducer.h.

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

std::vector<int> CTPPSFastTrackingProducer::recCellId_B

private

Definition at line 128 of file CTPPSFastTrackingProducer.h.

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

std::vector<int> CTPPSFastTrackingProducer::recCellId_F

private

Definition at line 128 of file CTPPSFastTrackingProducer.h.

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

std::vector<double> CTPPSFastTrackingProducer::recTof_B

private

Definition at line 129 of file CTPPSFastTrackingProducer.h.

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

std::vector<double> CTPPSFastTrackingProducer::recTof_F

private

Definition at line 129 of file CTPPSFastTrackingProducer.h.

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

std::vector<CTPPSFastTrack> CTPPSFastTrackingProducer::theCTPPSFastTrack

private

Definition at line 124 of file CTPPSFastTrackingProducer.h.

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

CTPPSFastTrack CTPPSFastTrackingProducer::track

private

Definition at line 126 of file CTPPSFastTrackingProducer.h.

Referenced by FastReco().

std::unique_ptr<CTPPSTrkStation> CTPPSFastTrackingProducer::TrkStation_B

private

Definition at line 116 of file CTPPSFastTrackingProducer.h.

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

std::unique_ptr<CTPPSTrkStation> CTPPSFastTrackingProducer::TrkStation_F

private

Definition at line 115 of file CTPPSFastTrackingProducer.h.

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