|
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
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
|
#ifndef Alignment_MuonAlignmentAlgorithms_CSCAlignmentCorrections_H
#define Alignment_MuonAlignmentAlgorithms_CSCAlignmentCorrections_H
/** \class CSCAlignmentCorrections
* $Date: 2010/05/12 21:57:51 $
* $Revision: 1.2 $
* \author J. Pivarski - Texas A&M University <pivarski@physics.tamu.edu>
*/
#include <fstream>
#include "TH1F.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/MuonDetId/interface/CSCDetId.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"
class CSCAlignmentCorrections {
public:
CSCAlignmentCorrections(std::string fitterName, double oldchi2, double newchi2)
: m_fitterName(fitterName), m_oldchi2(oldchi2), m_newchi2(newchi2) {}
virtual ~CSCAlignmentCorrections() {}
void insertCorrection(std::string name, CSCDetId id, double value) {
m_name.push_back(name);
m_id.push_back(id);
m_value.push_back(value);
};
void insertMode(const std::vector<double>& coefficient,
const std::vector<std::string>& modename,
const std::vector<long>& modeid,
double error) {
m_coefficient.push_back(coefficient);
m_modename.push_back(modename);
m_modeid.push_back(modeid);
m_error.push_back(error);
};
void insertResidual(std::string i, std::string j, double before, double uncert, double residual, double pull) {
m_i.push_back(i);
m_j.push_back(j);
m_before.push_back(before);
m_uncert.push_back(uncert);
m_residual.push_back(residual);
m_pull.push_back(pull);
};
void applyAlignment(AlignableNavigator* alignableNavigator,
AlignmentParameterStore* alignmentParameterStore,
int mode,
bool combineME11);
void plot();
void report(std::ofstream& report);
protected:
std::string m_fitterName;
double m_oldchi2, m_newchi2;
// there's one of these for each chamber
std::vector<std::string> m_name;
std::vector<CSCDetId> m_id;
std::vector<double> m_value;
// there's one of these for each error mode
std::vector<std::vector<double> > m_coefficient;
std::vector<std::vector<std::string> > m_modename;
std::vector<std::vector<long> > m_modeid;
std::vector<double> m_error;
// there's one of these for each constraint
std::vector<std::string> m_i;
std::vector<std::string> m_j;
std::vector<double> m_before;
std::vector<double> m_uncert;
std::vector<double> m_residual;
std::vector<double> m_pull;
std::vector<TH1F*> th1f_modes;
};
#endif // Alignment_MuonAlignmentAlgorithms_CSCAlignmentCorrections_H
|