|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
|
#ifndef Alignment_MuonAlignmentAlgorithms_CSCChamberFitter_H
#define Alignment_MuonAlignmentAlgorithms_CSCChamberFitter_H
/** \class CSCChamberFitter
* $Date: 2010/05/06 05:57:33 $
* $Revision: 1.2 $
* \author J. Pivarski - Texas A&M University <pivarski@physics.tamu.edu>
*/
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/MuonDetId/interface/CSCDetId.h"
#include "DataFormats/CLHEP/interface/AlgebraicObjects.h"
#include "Alignment/CommonAlignment/interface/Alignable.h"
#include "Alignment/CommonAlignment/interface/AlignableNavigator.h"
#include "Alignment/CommonAlignment/interface/AlignableObjectId.h"
#include "Alignment/CommonAlignment/interface/AlignmentParameters.h"
#include "Alignment/CommonAlignmentAlgorithm/interface/AlignmentParameterStore.h"
#include "Alignment/MuonAlignmentAlgorithms/interface/CSCPairConstraint.h"
#include "Alignment/MuonAlignmentAlgorithms/plugins/CSCPairResidualsConstraint.h"
#include "Alignment/MuonAlignmentAlgorithms/plugins/CSCAlignmentCorrections.h"
class CSCChamberFitter {
public:
CSCChamberFitter(const edm::ParameterSet &iConfig, std::vector<CSCPairResidualsConstraint *> &residualsConstraints);
virtual ~CSCChamberFitter() {}
bool fit(std::vector<CSCAlignmentCorrections *> &corrections) const;
void radiusCorrection(AlignableNavigator *alignableNavigator,
AlignmentParameterStore *alignmentParameterStore,
bool combineME11) const;
protected:
int index(std::string alignable) const;
void walk(std::map<int, bool> &touched, int alignable) const;
long alignableId(std::string alignable) const;
bool isFrame(int i) const;
double chi2(const AlgebraicVector &A, double lambda) const;
double lhsVector(int k) const;
double hessian(int k, int l, double lambda) const;
std::string m_name;
std::vector<std::string> m_alignables;
std::vector<int> m_frames;
int m_fixed;
std::vector<CSCPairConstraint *> m_constraints;
};
#endif // Alignment_MuonAlignmentAlgorithms_CSCChamberFitter_H
|