RecoProducerFP420 Class Reference

#include <RecoProducerFP420.h>

Public Member Functions
std::vector< RecoFP420 >	reconstruct (int, double, double, double, double, double, double)
	RecoProducerFP420 (const edm::ParameterSet &conf)
virtual	~RecoProducerFP420 ()

Private Attributes
std::string	beam1filename
std::string	beam2filename
edm::ParameterSet	conf_
double	length
H_BeamLine *	m_beamline1
H_BeamLine *	m_beamline2
double	m_e
double	m_q2
H_RecRPObject *	m_rp420_b
H_RecRPObject *	m_rp420_f
double	m_tx0
double	m_ty0
double	m_x0
double	m_y0
std::vector< RecoFP420 >	rhits
int	verbosity

Detailed Description

Definition at line 14 of file RecoProducerFP420.h.

Constructor & Destructor Documentation

RecoProducerFP420::RecoProducerFP420

(

const edm::ParameterSet &

conf

)

Definition at line 14 of file RecoProducerFP420.cc.

References b1, b2, beam1filename, beam2filename, conf_, gather_cfg::cout, alignCSCRings::e, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), length, m_beamline1, m_beamline2, mps_check::msg, AlCaHLTBitMon_QueryRunRegistry::string, verbosity, and cms::Exception::what().

                                                                : conf_(conf) {
  // Create LHC beam line
  //  double length  = param.getParameter<double>("BeamLineLength");
  //  std::cout << " BeamLineLength = " << length << std::endl;

  length = conf_.getParameter<double>("BeamLineLength");  //m
  verbosity = conf_.getUntrackedParameter<int>("VerbosityLevel");
  beam1filename = conf_.getParameter<std::string>("Beam1");
  beam2filename = conf_.getParameter<std::string>("Beam2");

  edm::LogInfo("RecoProducerFP420") << "RecoProducerFP420 parameters: \n"
                                    << "   Verbosity: " << verbosity << "\n"
                                    << "   length:    " << length << "\n";
  if (verbosity > 1) {
    std::cout << "  RecoProducerFP420: constructor    " << std::endl;
    std::cout << "   BeamLineLength:    " << length << std::endl;
  }

  //  edm::FileInPath b1("SimTransport/HectorData/twiss_ip5_b1_v6.5.txt");
  //  edm::FileInPath b2("SimTransport/HectorData/twiss_ip5_b2_v6.5.txt");

  edm::FileInPath b1(beam1filename.c_str());
  edm::FileInPath b2(beam2filename.c_str());

  m_beamline1 = new H_BeamLine(1, length + 0.1);   // (direction, length)
  m_beamline2 = new H_BeamLine(-1, length + 0.1);  //

  try {
    m_beamline1->fill(b1.fullPath(), 1, "IP5");
    m_beamline2->fill(b2.fullPath(), -1, "IP5");
  } catch (const edm::Exception& e) {
    std::string msg = e.what();
    msg += " caught in RecoProducerFP420... \nERROR: Could not locate SimTransport/HectorData data files.";
    edm::LogError("DataNotFound") << msg;
  }

  m_beamline1->offsetElements(120, -0.097);
  m_beamline2->offsetElements(120, +0.097);

  m_beamline1->calcMatrix();
  m_beamline2->calcMatrix();

  edm::LogInfo("RecoProducerFP420") << "==============================\n";
}

RecoProducerFP420::~RecoProducerFP420 ( )

virtual

Definition at line 59 of file RecoProducerFP420.cc.

{}

Member Function Documentation

std::vector< RecoFP420 > RecoProducerFP420::reconstruct	(	int	direction,
		double	x1_420,
		double	y1_420,
		double	x2_420,
		double	y2_420,
		double	z1_420,
		double	z2_420
	)

Definition at line 61 of file RecoProducerFP420.cc.

References gather_cfg::cout, m_beamline1, m_beamline2, m_e, m_q2, m_rp420_b, m_rp420_f, m_tx0, m_ty0, m_x0, m_y0, rhits, and verbosity.

Referenced by FP420RecoMain::run().

                                                                                                             {
  // ==================================================================================
  // ==================================================================================
  std::vector<RecoFP420> rhits;
  int restracks = 10;  // max # tracks
  rhits.reserve(restracks);
  rhits.clear();
  // ==================================================================================
  // trivial (TM), angle compensation (AM) and position compensation (PM) methods
  // #define TM 1    #define AM 2   #define PM 3
  m_tx0 = -100., m_ty0 = -100., m_x0 = -100., m_y0 = -100.;
  if (direction == 1) {
    m_rp420_f = new H_RecRPObject(z1_420 * 0.001, z2_420 * 0.001, *m_beamline1);  // m
    if (m_rp420_f) {
      if (verbosity > 1) {
        std::cout << "  RecoProducerFP420: input coord. in um   " << std::endl;
        std::cout << "   x1_420:    " << x1_420 << "   y1_420:    " << y1_420 << std::endl;
        std::cout << "   x2_420:    " << x2_420 << "   y2_420:    " << y2_420 << std::endl;
      }
      m_rp420_f->setPositions(x1_420, y1_420, x2_420, y2_420);  //input coord. in um
      m_e = m_rp420_f->getE(AM);                                // GeV
      //  std::cout << "   m_e1:    " << m_rp420_f->getE( TM ) << std::endl;
      //  std::cout << "   m_e2:    " << m_rp420_f->getE( AM ) << std::endl;
      m_tx0 = m_rp420_f->getTXIP();  // urad
      m_ty0 = m_rp420_f->getTYIP();  // urad
      m_x0 = m_rp420_f->getX0();     // um
      m_y0 = m_rp420_f->getY0();     // um
      m_q2 = m_rp420_f->getQ2();     // GeV^2
    }                                // if ( m_rp420_f
  }                                  // if ( dire
  else if (direction == 2) {
    m_rp420_b = new H_RecRPObject(z1_420 * 0.001, z2_420 * 0.001, *m_beamline2);  // m
    if (m_rp420_b) {
      m_rp420_b->setPositions(x1_420, y1_420, x2_420, y2_420);  // input coord. in um
      m_e = m_rp420_b->getE(AM);                                // GeV
      m_tx0 = m_rp420_b->getTXIP();                             // urad
      m_ty0 = m_rp420_b->getTYIP();                             // urad
      m_x0 = m_rp420_b->getX0();                                // um
      m_y0 = m_rp420_b->getY0();                                // um
      m_q2 = m_rp420_b->getQ2();                                // GeV^2
    }                                                           // if ( m_rp420_b
  } else {
    return rhits;
  }

  // ==============================
  if (verbosity > 1) {
    std::cout << "  RecoProducerFP420: rhits.push_back    " << std::endl;
    std::cout << "   m_e:    " << m_e << std::endl;
    std::cout << "   m_x0:    " << m_x0 << std::endl;
    std::cout << "   m_y0:    " << m_y0 << std::endl;
    std::cout << "   m_tx0:    " << m_tx0 << std::endl;
    std::cout << "   m_ty0:    " << m_ty0 << std::endl;
    std::cout << "   m_q2:    " << m_q2 << std::endl;
    std::cout << "   direction:    " << direction << std::endl;
  }
  rhits.push_back(RecoFP420(m_e, m_x0, m_y0, m_tx0, m_ty0, m_q2, direction));
  // ==============================
  return rhits;
  //============
}

Member Data Documentation

std::string RecoProducerFP420::beam1filename

private

Definition at line 38 of file RecoProducerFP420.h.

Referenced by RecoProducerFP420().

std::string RecoProducerFP420::beam2filename

private

Definition at line 39 of file RecoProducerFP420.h.

Referenced by RecoProducerFP420().

edm::ParameterSet RecoProducerFP420::conf_

private

Definition at line 23 of file RecoProducerFP420.h.

Referenced by RecoProducerFP420().

double RecoProducerFP420::length

private

Definition at line 28 of file RecoProducerFP420.h.

Referenced by RecoProducerFP420().

H_BeamLine* RecoProducerFP420::m_beamline1

private

Definition at line 32 of file RecoProducerFP420.h.

Referenced by reconstruct(), and RecoProducerFP420().

H_BeamLine* RecoProducerFP420::m_beamline2

private

Definition at line 33 of file RecoProducerFP420.h.

Referenced by reconstruct(), and RecoProducerFP420().

double RecoProducerFP420::m_e

private

Definition at line 29 of file RecoProducerFP420.h.

Referenced by reconstruct().

double RecoProducerFP420::m_q2

private

Definition at line 29 of file RecoProducerFP420.h.

Referenced by reconstruct().

H_RecRPObject* RecoProducerFP420::m_rp420_b

private

Definition at line 36 of file RecoProducerFP420.h.

Referenced by reconstruct().

H_RecRPObject* RecoProducerFP420::m_rp420_f

private

Definition at line 35 of file RecoProducerFP420.h.

Referenced by reconstruct().

double RecoProducerFP420::m_tx0

private

Definition at line 30 of file RecoProducerFP420.h.

Referenced by reconstruct().

double RecoProducerFP420::m_ty0

private

Definition at line 30 of file RecoProducerFP420.h.

Referenced by reconstruct().

double RecoProducerFP420::m_x0

private

Definition at line 29 of file RecoProducerFP420.h.

Referenced by reconstruct().

double RecoProducerFP420::m_y0

private

Definition at line 29 of file RecoProducerFP420.h.

Referenced by reconstruct().

std::vector<RecoFP420> RecoProducerFP420::rhits

private

Definition at line 24 of file RecoProducerFP420.h.

Referenced by reconstruct().

int RecoProducerFP420::verbosity

private

Definition at line 26 of file RecoProducerFP420.h.

Referenced by reconstruct(), and RecoProducerFP420().