CMS 3D CMS Logo

List of all members | Public Member Functions | Private Member Functions | Private Attributes
SurveyInputCSCfromPins Class Reference

#include <SurveyInputCSCfromPins.h>

Inheritance diagram for SurveyInputCSCfromPins:
SurveyInputBase edm::EDAnalyzer edm::EDConsumerBase

Public Member Functions

void analyze (const edm::Event &, const edm::EventSetup &) override
 Read ideal tracker geometry from DB. More...
 
 SurveyInputCSCfromPins (const edm::ParameterSet &)
 
- Public Member Functions inherited from SurveyInputBase
void beginJob () override
 Read data from input. More...
 
 ~SurveyInputBase () override
 
- Public Member Functions inherited from edm::EDAnalyzer
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzer ()
 
SerialTaskQueueglobalLuminosityBlocksQueue ()
 
SerialTaskQueueglobalRunsQueue ()
 
ModuleDescription const & moduleDescription () const
 
std::string workerType () const
 
 ~EDAnalyzer () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBaseoperator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Member Functions

void errors (double a, double b, bool missing1, bool missing2, double &dx_dx, double &dy_dy, double &dz_dz, double &phix_phix, double &phiz_phiz, double &dy_phix)
 
void fillAllRecords (Alignable *ali)
 
void orient (align::LocalVector LC1, align::LocalVector LC2, double a, double b, double &T, double &dx, double &dy, double &dz, double &PhX, double &PhZ)
 

Private Attributes

double m_errorX
 
double m_errorY
 
double m_errorZ
 
double m_missingErrorAngle
 
double m_missingErrorTranslation
 
std::string m_pinPositions
 
std::string m_rootFile
 
double m_stationErrorPhiX
 
double m_stationErrorPhiY
 
double m_stationErrorPhiZ
 
double m_stationErrorX
 
double m_stationErrorY
 
double m_stationErrorZ
 
bool m_verbose
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from SurveyInputBase
static void addComponent (Alignable *)
 
static Alignabledetector ()
 Get alignable detector as read from input. More...
 
- Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &)
 
static bool wantsGlobalLuminosityBlocks ()
 
static bool wantsGlobalRuns ()
 
static bool wantsStreamLuminosityBlocks ()
 
static bool wantsStreamRuns ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
- Protected Attributes inherited from SurveyInputBase
bool theFirstEvent
 

Detailed Description

Class to read ideal tracker from DB.

$Date: Fri Jun 29 09:20:52 CEST 2007

Revision
1.2
Author
Dmitry Yakorev

Definition at line 16 of file SurveyInputCSCfromPins.h.

Constructor & Destructor Documentation

SurveyInputCSCfromPins::SurveyInputCSCfromPins ( const edm::ParameterSet cfg)

Definition at line 25 of file SurveyInputCSCfromPins.cc.

26  : m_pinPositions(cfg.getParameter<std::string>("pinPositions")),
27  m_rootFile(cfg.getParameter<std::string>("rootFile")),
28  m_verbose(cfg.getParameter<bool>("verbose")),
29  m_errorX(cfg.getParameter<double>("errorX")),
30  m_errorY(cfg.getParameter<double>("errorY")),
31  m_errorZ(cfg.getParameter<double>("errorZ")),
32  m_missingErrorTranslation(cfg.getParameter<double>("missingErrorTranslation")),
33  m_missingErrorAngle(cfg.getParameter<double>("missingErrorAngle")),
34  m_stationErrorX(cfg.getParameter<double>("stationErrorX")),
35  m_stationErrorY(cfg.getParameter<double>("stationErrorY")),
36  m_stationErrorZ(cfg.getParameter<double>("stationErrorZ")),
37  m_stationErrorPhiX(cfg.getParameter<double>("stationErrorPhiX")),
38  m_stationErrorPhiY(cfg.getParameter<double>("stationErrorPhiY")),
39  m_stationErrorPhiZ(cfg.getParameter<double>("stationErrorPhiZ")) {}
T getParameter(std::string const &) const

Member Function Documentation

void SurveyInputCSCfromPins::analyze ( const edm::Event ,
const edm::EventSetup iSetup 
)
overridevirtual

Read ideal tracker geometry from DB.

Implements SurveyInputBase.

Definition at line 131 of file SurveyInputCSCfromPins.cc.

References a, SurveyInputBase::addComponent(), AlignableNavigator::alignableFromDetId(), align::AlignableMuon, b, AlignableMuon::CSCEndcaps(), PVValHelper::dx, PVValHelper::dy, PVValHelper::dz, relativeConstraints::error, errors(), edmOneToOneComparison::file1, fillAllRecords(), edm::EventSetup::get(), h, mps_fire::i, cuy::ii, recoMuon::in, m_missingErrorAngle, m_pinPositions, m_rootFile, m_verbose, Alignable::mother(), Alignable::move(), orient(), Alignable::rotateAroundLocalX(), Alignable::rotateAroundLocalZ(), Alignable::setSurvey(), Alignable::surface(), Alignable::survey(), SurveyInputBase::theFirstEvent, AlignableSurface::toGlobal(), AlignableSurface::toLocal(), PV3DBase< T, PVType, FrameType >::x(), globals_cff::x1, globals_cff::x2, PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

131  {
132  if (theFirstEvent) {
133  edm::LogInfo("SurveyInputCSCfromPins") << "***************ENTERING INITIALIZATION******************"
134  << " \n";
135 
136  std::ifstream in;
137  in.open(m_pinPositions.c_str());
138 
139  Double_t x1, y1, z1, x2, y2, z2, a, b, tot = 0.0, maxErr = 0.0, h, s1, dx, dy, dz, PhX, PhZ, T;
140 
141  int ID1, ID2, ID3, ID4, ID5, i = 1, ii = 0;
142 
143  TFile *file1 = new TFile(m_rootFile.c_str(), "recreate");
144  TTree *tree1 = new TTree("tree1", "alignment pins");
145 
146  if (m_verbose) {
147  tree1->Branch("displacement_x_pin1_cm", &x1, "x1/D");
148  tree1->Branch("displacement_y_pin1_cm", &y1, "y1/D");
149  tree1->Branch("displacement_z_pin1_cm", &z1, "z1/D");
150  tree1->Branch("displacement_x_pin2_cm", &x2, "x2/D");
151  tree1->Branch("displacement_y_pin2_cm", &y2, "y2/D");
152  tree1->Branch("displacement_z_pin2_cm", &z2, "z2/D");
153  tree1->Branch("error_vector_length_cm", &h, "h/D");
154  tree1->Branch("stretch_diff_cm", &s1, "s1/D");
155  tree1->Branch("stretch_factor", &T, "T/D");
156  tree1->Branch("chamber_displacement_x_cm", &dx, "dx/D");
157  tree1->Branch("chamber_displacement_y_cm", &dy, "dy/D");
158  tree1->Branch("chamber_displacement_z_cm", &dz, "dz/D");
159  tree1->Branch("chamber_rotation_x_rad", &PhX, "PhX/D");
160  tree1->Branch("chamber_rotation_z_rad", &PhZ, "PhZ/D");
161  }
162 
163  edm::ESHandle<DTGeometry> dtGeometry;
164  edm::ESHandle<CSCGeometry> cscGeometry;
165  iSetup.get<MuonGeometryRecord>().get(dtGeometry);
166  iSetup.get<MuonGeometryRecord>().get(cscGeometry);
167 
168  AlignableMuon *theAlignableMuon = new AlignableMuon(&(*dtGeometry), &(*cscGeometry));
169  AlignableNavigator *theAlignableNavigator = new AlignableNavigator(theAlignableMuon);
170 
171  const auto &theEndcaps = theAlignableMuon->CSCEndcaps();
172 
173  for (const auto &aliiter : theEndcaps) {
174  addComponent(aliiter);
175  }
176 
177  while (in.good()) {
178  bool missing1 = false;
179  bool missing2 = false;
180 
181  in >> ID1 >> ID2 >> ID3 >> ID4 >> ID5 >> x1 >> y1 >> z1 >> x2 >> y2 >> z2 >> a >> b;
182 
183  if (fabs(x1 - 1000.) < 1e5 && fabs(y1 - 1000.) < 1e5 && fabs(z1 - 1000.) < 1e5) {
184  missing1 = true;
185  x1 = x2;
186  y1 = y2;
187  z1 = z2;
188  }
189 
190  if (fabs(x2 - 1000.) < 1e5 && fabs(y2 - 1000.) < 1e5 && fabs(z2 - 1000.) < 1e5) {
191  missing2 = true;
192  x2 = x1;
193  y2 = y1;
194  z2 = z1;
195  }
196 
197  x1 = x1 / 10.0;
198  y1 = y1 / 10.0;
199  z1 = z1 / 10.0;
200  x2 = x2 / 10.0;
201  y2 = y2 / 10.0;
202  z2 = z2 / 10.0;
203 
204  CSCDetId layerID(ID1, ID2, ID3, ID4, 1);
205 
206  // We cannot use chamber ID (when ID5=0), because AlignableNavigator gives the error (aliDet and aliDetUnit are undefined for chambers)
207 
208  Alignable *theAlignable1 = theAlignableNavigator->alignableFromDetId(layerID);
209  Alignable *chamberAli = theAlignable1->mother();
210 
211  LocalVector LC1 = chamberAli->surface().toLocal(GlobalVector(x1, y1, z1));
212  LocalVector LC2 = chamberAli->surface().toLocal(GlobalVector(x2, y2, z2));
213 
214  orient(LC1, LC2, a, b, T, dx, dy, dz, PhX, PhZ);
215 
216  GlobalPoint PG1 = chamberAli->surface().toGlobal(LocalPoint(LC1.x(), LC1.y() + a, LC1.z() + b));
217  chamberAli->surface().toGlobal(LocalPoint(LC2.x(), LC2.y() - a, LC2.z() + b));
218 
219  LocalVector lvector(dx, dy, dz);
220  GlobalVector gvector = (chamberAli->surface()).toGlobal(lvector);
221  chamberAli->move(gvector);
222 
223  chamberAli->rotateAroundLocalX(PhX);
224  chamberAli->rotateAroundLocalZ(PhZ);
225 
226  double dx_dx, dy_dy, dz_dz, phix_phix, phiz_phiz, dy_phix;
227  errors(a, b, missing1, missing2, dx_dx, dy_dy, dz_dz, phix_phix, phiz_phiz, dy_phix);
228  align::ErrorMatrix error = ROOT::Math::SMatrixIdentity();
229  error(0, 0) = dx_dx;
230  error(1, 1) = dy_dy;
231  error(2, 2) = dz_dz;
232  error(3, 3) = phix_phix;
233  error(4, 4) = m_missingErrorAngle;
234  error(5, 5) = phiz_phiz;
235  error(1, 3) = dy_phix;
236  error(3, 1) = dy_phix; // just in case
237 
238  chamberAli->setSurvey(new SurveyDet(chamberAli->surface(), error));
239 
240  if (m_verbose) {
241  edm::LogInfo("SurveyInputCSCfromPins") << " survey information = " << chamberAli->survey() << " \n";
242 
243  LocalPoint LP1n = chamberAli->surface().toLocal(PG1);
244 
245  LocalPoint hiP(LP1n.x(), LP1n.y() - a * T, LP1n.z() - b);
246 
247  h = hiP.mag();
248  s1 = LP1n.y() - a;
249 
250  if (h > maxErr) {
251  maxErr = h;
252 
253  ii = i;
254  }
255 
256  edm::LogInfo("SurveyInputCSCfromPins")
257  << " \n"
258  << "i " << i++ << " " << ID1 << " " << ID2 << " " << ID3 << " " << ID4 << " " << ID5 << " error " << h
259  << " \n"
260  << " x1 " << x1 << " y1 " << y1 << " z1 " << z1 << " x2 " << x2 << " y2 " << y2 << " z2 " << z2 << " \n"
261  << " error " << h << " S1 " << s1 << " \n"
262  << " dx " << dx << " dy " << dy << " dz " << dz << " PhX " << PhX << " PhZ " << PhZ << " \n";
263 
264  tot += h;
265 
266  tree1->Fill();
267  }
268  }
269 
270  in.close();
271 
272  if (m_verbose) {
273  file1->Write();
274  edm::LogInfo("SurveyInputCSCfromPins")
275  << " Total error " << tot << " Max Error " << maxErr << " N " << ii << " \n";
276  }
277 
278  file1->Close();
279 
280  for (const auto &aliiter : theEndcaps) {
281  fillAllRecords(aliiter);
282  }
283 
284  delete theAlignableMuon;
285  delete theAlignableNavigator;
286 
287  edm::LogInfo("SurveyInputCSCfromPins") << "*************END INITIALIZATION***************"
288  << " \n";
289 
290  theFirstEvent = false;
291  }
292 }
void errors(double a, double b, bool missing1, bool missing2, double &dx_dx, double &dy_dy, double &dz_dz, double &phix_phix, double &phiz_phiz, double &dy_phix)
Point3DBase< Scalar, LocalTag > LocalPoint
Definition: Definitions.h:32
FWCore Framework interface EventSetupRecordImplementation h
Helper function to determine trigger accepts.
virtual void rotateAroundLocalZ(Scalar radians)
Rotation around local z-axis.
Definition: Alignable.cc:233
const SurveyDet * survey() const
Return survey info.
Definition: Alignable.h:225
T y() const
Definition: PV3DBase.h:63
virtual void move(const GlobalVector &displacement)=0
Movement with respect to the global reference frame.
static void addComponent(Alignable *)
align::Alignables CSCEndcaps()
align::RotationType toLocal(const align::RotationType &) const
Return in local frame a rotation given in global frame.
void fillAllRecords(Alignable *ali)
T z() const
Definition: PV3DBase.h:64
const AlignableSurface & surface() const
Return the Surface (global position and orientation) of the object.
Definition: Alignable.h:135
ii
Definition: cuy.py:590
void setSurvey(const SurveyDet *)
Set survey info.
Definition: Alignable.cc:340
double b
Definition: hdecay.h:120
double a
Definition: hdecay.h:121
T get() const
Definition: EventSetup.h:71
virtual void rotateAroundLocalX(Scalar radians)
Rotation around local x-axis.
Definition: Alignable.cc:181
align::GlobalPoints toGlobal(const align::LocalPoints &) const
Return in global coord given a set of local points.
math::Error< 6 >::type ErrorMatrix
Definition: Definitions.h:39
long double T
Constructor of the full muon geometry.
Definition: AlignableMuon.h:37
T x() const
Definition: PV3DBase.h:62
void orient(align::LocalVector LC1, align::LocalVector LC2, double a, double b, double &T, double &dx, double &dy, double &dz, double &PhX, double &PhZ)
AlignableDetOrUnitPtr alignableFromDetId(const DetId &detid)
Returns AlignableDetOrUnitPtr corresponding to given DetId.
Alignable * mother() const
Return pointer to container alignable (if any)
Definition: Alignable.h:94
Global3DVector GlobalVector
Definition: GlobalVector.h:10
void SurveyInputCSCfromPins::errors ( double  a,
double  b,
bool  missing1,
bool  missing2,
double &  dx_dx,
double &  dy_dy,
double &  dz_dz,
double &  phix_phix,
double &  phiz_phiz,
double &  dy_phix 
)
private

Definition at line 75 of file SurveyInputCSCfromPins.cc.

References PVValHelper::dx, PVValHelper::dy, PVValHelper::dz, mps_fire::i, m_errorX, m_errorY, m_errorZ, m_missingErrorTranslation, and orient().

Referenced by analyze(), and fillAllRecords().

84  {
85  dx_dx = 0.;
86  dy_dy = 0.;
87  dz_dz = 0.;
88  phix_phix = 0.;
89  phiz_phiz = 0.;
90  dy_phix = 0.;
91 
92  const double trials = 10000.; // two significant digits
93 
94  for (int i = 0; i < trials; i++) {
95  LocalVector LC1, LC2;
96  if (missing1) {
97  LC1 = LocalVector(gRandom->Gaus(0., m_missingErrorTranslation),
98  gRandom->Gaus(0., m_missingErrorTranslation),
99  gRandom->Gaus(0., m_missingErrorTranslation));
100  } else {
101  LC1 = LocalVector(gRandom->Gaus(0., m_errorX), gRandom->Gaus(0., m_errorY), gRandom->Gaus(0., m_errorZ));
102  }
103 
104  if (missing2) {
105  LC2 = LocalVector(gRandom->Gaus(0., m_missingErrorTranslation),
106  gRandom->Gaus(0., m_missingErrorTranslation),
107  gRandom->Gaus(0., m_missingErrorTranslation));
108  } else {
109  LC2 = LocalVector(gRandom->Gaus(0., m_errorX), gRandom->Gaus(0., m_errorY), gRandom->Gaus(0., m_errorZ));
110  }
111 
112  double dx, dy, dz, PhX, PhZ, T;
113  orient(LC1, LC2, a, b, T, dx, dy, dz, PhX, PhZ);
114 
115  dx_dx += dx * dx;
116  dy_dy += dy * dy;
117  dz_dz += dz * dz;
118  phix_phix += PhX * PhX;
119  phiz_phiz += PhZ * PhZ;
120  dy_phix += dy * PhX; // the only non-zero off-diagonal element
121  }
122 
123  dx_dx /= trials;
124  dy_dy /= trials;
125  dz_dz /= trials;
126  phix_phix /= trials;
127  phiz_phiz /= trials;
128  dy_phix /= trials;
129 }
Local3DVector LocalVector
Definition: LocalVector.h:12
double b
Definition: hdecay.h:120
double a
Definition: hdecay.h:121
long double T
void orient(align::LocalVector LC1, align::LocalVector LC2, double a, double b, double &T, double &dx, double &dy, double &dz, double &PhX, double &PhZ)
void SurveyInputCSCfromPins::fillAllRecords ( Alignable ali)
private

Definition at line 294 of file SurveyInputCSCfromPins.cc.

References a, funct::abs(), b, Alignable::components(), DEFINE_FWK_MODULE, MillePedeFileConverter_cfg::e, relativeConstraints::error, errors(), Alignable::geomDetId(), m_missingErrorAngle, m_missingErrorTranslation, m_stationErrorPhiX, m_stationErrorPhiY, m_stationErrorPhiZ, m_stationErrorX, m_stationErrorY, m_stationErrorZ, Alignable::setSurvey(), Alignable::surface(), and Alignable::survey().

Referenced by analyze().

294  {
295  if (ali->survey() == nullptr) {
296  AlignableCSCChamber *ali_AlignableCSCChamber = dynamic_cast<AlignableCSCChamber *>(ali);
297  AlignableCSCStation *ali_AlignableCSCStation = dynamic_cast<AlignableCSCStation *>(ali);
298 
299  if (ali_AlignableCSCChamber != nullptr) {
300  CSCDetId detid(ali->geomDetId());
301  if (abs(detid.station()) == 1 && (detid.ring() == 1 || detid.ring() == 4)) {
302  align::ErrorMatrix error = ROOT::Math::SMatrixIdentity();
306  error(3, 3) = m_missingErrorAngle;
307  error(4, 4) = m_missingErrorAngle;
308  error(5, 5) = m_missingErrorAngle;
309 
310  ali->setSurvey(new SurveyDet(ali->surface(), error));
311  } else {
312  double a = 100.;
313  double b = -9.4034;
314  if (abs(detid.station()) == 1 && detid.ring() == 2)
315  a = -90.260;
316  else if (abs(detid.station()) == 1 && detid.ring() == 3)
317  a = -85.205;
318  else if (abs(detid.station()) == 2 && detid.ring() == 1)
319  a = -97.855;
320  else if (abs(detid.station()) == 2 && detid.ring() == 2)
321  a = -164.555;
322  else if (abs(detid.station()) == 3 && detid.ring() == 1)
323  a = -87.870;
324  else if (abs(detid.station()) == 3 && detid.ring() == 2)
325  a = -164.555;
326  else if (abs(detid.station()) == 4 && detid.ring() == 1)
327  a = -77.890;
328 
329  double dx_dx, dy_dy, dz_dz, phix_phix, phiz_phiz, dy_phix;
330  errors(a, b, true, true, dx_dx, dy_dy, dz_dz, phix_phix, phiz_phiz, dy_phix);
331  align::ErrorMatrix error = ROOT::Math::SMatrixIdentity();
332  error(0, 0) = dx_dx;
333  error(1, 1) = dy_dy;
334  error(2, 2) = dz_dz;
335  error(3, 3) = phix_phix;
336  error(4, 4) = m_missingErrorAngle;
337  error(5, 5) = phiz_phiz;
338  error(1, 3) = dy_phix;
339  error(3, 1) = dy_phix; // just in case
340 
341  ali->setSurvey(new SurveyDet(ali->surface(), error));
342  }
343  }
344 
345  else if (ali_AlignableCSCStation != nullptr) {
346  align::ErrorMatrix error = ROOT::Math::SMatrixIdentity();
347  error(0, 0) = m_stationErrorX;
348  error(1, 1) = m_stationErrorY;
349  error(2, 2) = m_stationErrorZ;
350  error(3, 3) = m_stationErrorPhiX;
351  error(4, 4) = m_stationErrorPhiY;
352  error(5, 5) = m_stationErrorPhiZ;
353 
354  ali->setSurvey(new SurveyDet(ali->surface(), error));
355  }
356 
357  else {
358  align::ErrorMatrix error = ROOT::Math::SMatrixIdentity();
359  ali->setSurvey(new SurveyDet(ali->surface(), error * (1e-10)));
360  }
361  }
362 
363  for (const auto &iter : ali->components()) {
364  fillAllRecords(iter);
365  }
366 }
void errors(double a, double b, bool missing1, bool missing2, double &dx_dx, double &dy_dy, double &dz_dz, double &phix_phix, double &phiz_phiz, double &dy_phix)
const SurveyDet * survey() const
Return survey info.
Definition: Alignable.h:225
virtual const Alignables & components() const =0
Return vector of all direct components.
void fillAllRecords(Alignable *ali)
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
const AlignableSurface & surface() const
Return the Surface (global position and orientation) of the object.
Definition: Alignable.h:135
void setSurvey(const SurveyDet *)
Set survey info.
Definition: Alignable.cc:340
double b
Definition: hdecay.h:120
double a
Definition: hdecay.h:121
math::Error< 6 >::type ErrorMatrix
Definition: Definitions.h:39
const DetId & geomDetId() const
Definition: Alignable.h:186
A muon CSC Chamber( an AlignableDet )
void SurveyInputCSCfromPins::orient ( align::LocalVector  LC1,
align::LocalVector  LC2,
double  a,
double  b,
double &  T,
double &  dx,
double &  dy,
double &  dz,
double &  PhX,
double &  PhZ 
)
private

Definition at line 41 of file SurveyInputCSCfromPins.cc.

References a, b, PV3DBase< T, PVType, FrameType >::mag(), SQR, mathSSE::sqrt(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by analyze(), and errors().

50  {
51  double cosPhX, sinPhX, cosPhZ, sinPhZ;
52 
53  LocalPoint LP1(LC1.x(), LC1.y() + a, LC1.z() + b);
54  LocalPoint LP2(LC2.x(), LC2.y() - a, LC2.z() + b);
55 
56  LocalPoint P((LP1.x() - LP2.x()) / (2. * a), (LP1.y() - LP2.y()) / (2. * a), (LP1.z() - LP2.z()) / (2. * a));
57  LocalPoint Pp((LP1.x() + LP2.x()) / (2.), (LP1.y() + LP2.y()) / (2.), (LP1.z() + LP2.z()) / (2.));
58 
59  T = P.mag();
60 
61  sinPhX = P.z() / T;
62  cosPhX = sqrt(1 - SQR(sinPhX));
63  cosPhZ = P.y() / (T * cosPhX);
64  sinPhZ = -P.x() / (T * cosPhZ);
65 
66  PhX = atan2(sinPhX, cosPhX);
67 
68  PhZ = atan2(sinPhZ, cosPhZ);
69 
70  dx = Pp.x() - sinPhZ * sinPhX * b;
71  dy = Pp.y() + cosPhZ * sinPhX * b;
72  dz = Pp.z() - cosPhX * b;
73 }
T y() const
Definition: PV3DBase.h:63
T mag() const
Definition: PV3DBase.h:67
T sqrt(T t)
Definition: SSEVec.h:18
T z() const
Definition: PV3DBase.h:64
double b
Definition: hdecay.h:120
std::pair< OmniClusterRef, TrackingParticleRef > P
double a
Definition: hdecay.h:121
long double T
T x() const
Definition: PV3DBase.h:62
#define SQR(x)

Member Data Documentation

double SurveyInputCSCfromPins::m_errorX
private

Definition at line 50 of file SurveyInputCSCfromPins.h.

Referenced by errors().

double SurveyInputCSCfromPins::m_errorY
private

Definition at line 50 of file SurveyInputCSCfromPins.h.

Referenced by errors().

double SurveyInputCSCfromPins::m_errorZ
private

Definition at line 50 of file SurveyInputCSCfromPins.h.

Referenced by errors().

double SurveyInputCSCfromPins::m_missingErrorAngle
private

Definition at line 51 of file SurveyInputCSCfromPins.h.

Referenced by analyze(), and fillAllRecords().

double SurveyInputCSCfromPins::m_missingErrorTranslation
private

Definition at line 51 of file SurveyInputCSCfromPins.h.

Referenced by errors(), and fillAllRecords().

std::string SurveyInputCSCfromPins::m_pinPositions
private

Definition at line 47 of file SurveyInputCSCfromPins.h.

Referenced by analyze().

std::string SurveyInputCSCfromPins::m_rootFile
private

Definition at line 48 of file SurveyInputCSCfromPins.h.

Referenced by analyze().

double SurveyInputCSCfromPins::m_stationErrorPhiX
private

Definition at line 52 of file SurveyInputCSCfromPins.h.

Referenced by fillAllRecords().

double SurveyInputCSCfromPins::m_stationErrorPhiY
private

Definition at line 52 of file SurveyInputCSCfromPins.h.

Referenced by fillAllRecords().

double SurveyInputCSCfromPins::m_stationErrorPhiZ
private

Definition at line 52 of file SurveyInputCSCfromPins.h.

Referenced by fillAllRecords().

double SurveyInputCSCfromPins::m_stationErrorX
private

Definition at line 52 of file SurveyInputCSCfromPins.h.

Referenced by fillAllRecords().

double SurveyInputCSCfromPins::m_stationErrorY
private

Definition at line 52 of file SurveyInputCSCfromPins.h.

Referenced by fillAllRecords().

double SurveyInputCSCfromPins::m_stationErrorZ
private

Definition at line 52 of file SurveyInputCSCfromPins.h.

Referenced by fillAllRecords().

bool SurveyInputCSCfromPins::m_verbose
private

Definition at line 49 of file SurveyInputCSCfromPins.h.

Referenced by analyze().