CMS 3D CMS Logo
List of all members | Public Member Functions | Private Attributes
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::RecoProducerFP420 ( const edm::ParameterSet &  conf)

Definition at line 11 of file RecoProducerFP420.cc.

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

11  : conf_(conf) {
12  // Create LHC beam line
13  // double length = param.getParameter<double>("BeamLineLength");
14  // std::cout << " BeamLineLength = " << length << std::endl;
15 
16  length = conf_.getParameter<double>("BeamLineLength"); //m
17  verbosity = conf_.getUntrackedParameter<int>("VerbosityLevel");
18  beam1filename = conf_.getParameter<std::string>("Beam1");
19  beam2filename = conf_.getParameter<std::string>("Beam2");
20 
21  edm::LogInfo("RecoProducerFP420") << "RecoProducerFP420 parameters: \n"
22  << " Verbosity: " << verbosity << "\n"
23  << " length: " << length << "\n";
24  if (verbosity > 1) {
25  std::cout << " RecoProducerFP420: constructor " << std::endl;
26  std::cout << " BeamLineLength: " << length << std::endl;
27  }
28 
29  // edm::FileInPath b1("SimTransport/HectorData/twiss_ip5_b1_v6.5.txt");
30  // edm::FileInPath b2("SimTransport/HectorData/twiss_ip5_b2_v6.5.txt");
31 
32  edm::FileInPath b1(beam1filename.c_str());
33  edm::FileInPath b2(beam2filename.c_str());
34 
35  m_beamline1 = new H_BeamLine(1, length + 0.1); // (direction, length)
36  m_beamline2 = new H_BeamLine(-1, length + 0.1); //
37 
38  try {
39  m_beamline1->fill(b1.fullPath(), 1, "IP5");
40  m_beamline2->fill(b2.fullPath(), -1, "IP5");
41  } catch (const edm::Exception& e) {
42  std::string msg = e.what();
43  msg += " caught in RecoProducerFP420... \nERROR: Could not locate SimTransport/HectorData data files.";
44  edm::LogError("DataNotFound") << msg;
45  }
46 
47  m_beamline1->offsetElements(120, -0.097);
48  m_beamline2->offsetElements(120, +0.097);
49 
50  m_beamline1->calcMatrix();
51  m_beamline2->calcMatrix();
52 
53  edm::LogInfo("RecoProducerFP420") << "==============================\n";
54 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
RecoProducerFP420::beam1filename
std::string beam1filename
Definition: RecoProducerFP420.h:38
RecoProducerFP420::m_beamline2
H_BeamLine * m_beamline2
Definition: RecoProducerFP420.h:33
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:123
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
RecoProducerFP420::conf_
edm::ParameterSet conf_
Definition: RecoProducerFP420.h:23
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
RecoProducerFP420::beam2filename
std::string beam2filename
Definition: RecoProducerFP420.h:39
mps_check.msg
tuple msg
Definition: mps_check.py:285
gather_cfg.cout
cout
Definition: gather_cfg.py:144
RecoProducerFP420::m_beamline1
H_BeamLine * m_beamline1
Definition: RecoProducerFP420.h:32
b2
static constexpr float b2
Definition: L1EGammaCrystalsEmulatorProducer.cc:80
b1
static constexpr float b1
Definition: L1EGammaCrystalsEmulatorProducer.cc:80

◆ ~RecoProducerFP420()

RecoProducerFP420::~RecoProducerFP420 ( )
virtual

Definition at line 56 of file RecoProducerFP420.cc.

56 {}

Member Function Documentation

◆ reconstruct()

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

59  {
60  // ==================================================================================
61  // ==================================================================================
62  std::vector<RecoFP420> rhits;
63  int restracks = 10; // max # tracks
64  rhits.reserve(restracks);
65  rhits.clear();
66  // ==================================================================================
67  // trivial (TM), angle compensation (AM) and position compensation (PM) methods
68  // #define TM 1 #define AM 2 #define PM 3
69  m_tx0 = -100., m_ty0 = -100., m_x0 = -100., m_y0 = -100.;
70  if (direction == 1) {
71  m_rp420_f = new H_RecRPObject(z1_420 * 0.001, z2_420 * 0.001, *m_beamline1); // m
72  if (m_rp420_f) {
73  if (verbosity > 1) {
74  std::cout << " RecoProducerFP420: input coord. in um " << std::endl;
75  std::cout << " x1_420: " << x1_420 << " y1_420: " << y1_420 << std::endl;
76  std::cout << " x2_420: " << x2_420 << " y2_420: " << y2_420 << std::endl;
77  }
78  m_rp420_f->setPositions(x1_420, y1_420, x2_420, y2_420); //input coord. in um
79  m_e = m_rp420_f->getE(AM); // GeV
80  // std::cout << " m_e1: " << m_rp420_f->getE( TM ) << std::endl;
81  // std::cout << " m_e2: " << m_rp420_f->getE( AM ) << std::endl;
82  m_tx0 = m_rp420_f->getTXIP(); // urad
83  m_ty0 = m_rp420_f->getTYIP(); // urad
84  m_x0 = m_rp420_f->getX0(); // um
85  m_y0 = m_rp420_f->getY0(); // um
86  m_q2 = m_rp420_f->getQ2(); // GeV^2
87  } // if ( m_rp420_f
88  } // if ( dire
89  else if (direction == 2) {
90  m_rp420_b = new H_RecRPObject(z1_420 * 0.001, z2_420 * 0.001, *m_beamline2); // m
91  if (m_rp420_b) {
92  m_rp420_b->setPositions(x1_420, y1_420, x2_420, y2_420); // input coord. in um
93  m_e = m_rp420_b->getE(AM); // GeV
94  m_tx0 = m_rp420_b->getTXIP(); // urad
95  m_ty0 = m_rp420_b->getTYIP(); // urad
96  m_x0 = m_rp420_b->getX0(); // um
97  m_y0 = m_rp420_b->getY0(); // um
98  m_q2 = m_rp420_b->getQ2(); // GeV^2
99  } // if ( m_rp420_b
100  } else {
101  return rhits;
102  }
103 
104  // ==============================
105  if (verbosity > 1) {
106  std::cout << " RecoProducerFP420: rhits.push_back " << std::endl;
107  std::cout << " m_e: " << m_e << std::endl;
108  std::cout << " m_x0: " << m_x0 << std::endl;
109  std::cout << " m_y0: " << m_y0 << std::endl;
110  std::cout << " m_tx0: " << m_tx0 << std::endl;
111  std::cout << " m_ty0: " << m_ty0 << std::endl;
112  std::cout << " m_q2: " << m_q2 << std::endl;
113  std::cout << " direction: " << direction << std::endl;
114  }
115  rhits.push_back(RecoFP420(m_e, m_x0, m_y0, m_tx0, m_ty0, m_q2, direction));
116  // ==============================
118  return rhits;
119  //============
120 }
RecoProducerFP420::rhits
std::vector< RecoFP420 > rhits
Definition: RecoProducerFP420.h:24
RecoProducerFP420::m_beamline2
H_BeamLine * m_beamline2
Definition: RecoProducerFP420.h:33
RecoProducerFP420::m_rp420_f
H_RecRPObject * m_rp420_f
Definition: RecoProducerFP420.h:35
RecoProducerFP420::m_rp420_b
H_RecRPObject * m_rp420_b
Definition: RecoProducerFP420.h:36
gather_cfg.cout
cout
Definition: gather_cfg.py:144
RecoProducerFP420::m_beamline1
H_BeamLine * m_beamline1
Definition: RecoProducerFP420.h:32

Member Data Documentation

◆ beam1filename

std::string RecoProducerFP420::beam1filename
private

Definition at line 38 of file RecoProducerFP420.h.

Referenced by RecoProducerFP420().

◆ beam2filename

std::string RecoProducerFP420::beam2filename
private

Definition at line 39 of file RecoProducerFP420.h.

Referenced by RecoProducerFP420().

◆ conf_

edm::ParameterSet RecoProducerFP420::conf_
private

Definition at line 23 of file RecoProducerFP420.h.

Referenced by RecoProducerFP420().

◆ length

double RecoProducerFP420::length
private

Definition at line 28 of file RecoProducerFP420.h.

Referenced by RecoProducerFP420().

◆ m_beamline1

H_BeamLine* RecoProducerFP420::m_beamline1
private

Definition at line 32 of file RecoProducerFP420.h.

Referenced by reconstruct(), and RecoProducerFP420().

◆ m_beamline2

H_BeamLine* RecoProducerFP420::m_beamline2
private

Definition at line 33 of file RecoProducerFP420.h.

Referenced by reconstruct(), and RecoProducerFP420().

◆ m_e

double RecoProducerFP420::m_e
private

Definition at line 29 of file RecoProducerFP420.h.

Referenced by reconstruct().

◆ m_q2

double RecoProducerFP420::m_q2
private

Definition at line 29 of file RecoProducerFP420.h.

Referenced by reconstruct().

◆ m_rp420_b

H_RecRPObject* RecoProducerFP420::m_rp420_b
private

Definition at line 36 of file RecoProducerFP420.h.

Referenced by reconstruct().

◆ m_rp420_f

H_RecRPObject* RecoProducerFP420::m_rp420_f
private

Definition at line 35 of file RecoProducerFP420.h.

Referenced by reconstruct().

◆ m_tx0

double RecoProducerFP420::m_tx0
private

Definition at line 30 of file RecoProducerFP420.h.

Referenced by reconstruct().

◆ m_ty0

double RecoProducerFP420::m_ty0
private

Definition at line 30 of file RecoProducerFP420.h.

Referenced by reconstruct().

◆ m_x0

double RecoProducerFP420::m_x0
private

Definition at line 29 of file RecoProducerFP420.h.

Referenced by reconstruct().

◆ m_y0

double RecoProducerFP420::m_y0
private

Definition at line 29 of file RecoProducerFP420.h.

Referenced by reconstruct().

◆ rhits

std::vector<RecoFP420> RecoProducerFP420::rhits
private

Definition at line 24 of file RecoProducerFP420.h.

Referenced by reconstruct().

◆ verbosity

int RecoProducerFP420::verbosity
private

Definition at line 26 of file RecoProducerFP420.h.

Referenced by reconstruct(), and RecoProducerFP420().

Generated for CMSSW Reference Manual by doxygen 1.8.14