CMS 3D CMS Logo

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

#include <Alignment/CSCOverlapsAlignmentAlgorithm/interface/CSCOverlapsAlignmentAlgorithm.h>

Inheritance diagram for CSCOverlapsAlignmentAlgorithm:
AlignmentAlgorithmBase

Public Member Functions

 CSCOverlapsAlignmentAlgorithm (const edm::ParameterSet &iConfig, edm::ConsumesCollector &)
 
void initialize (const edm::EventSetup &iSetup, AlignableTracker *alignableTracker, AlignableMuon *alignableMuon, AlignableExtras *alignableExtras, AlignmentParameterStore *alignmentParameterStore) override
 Call at beginning of job (must be implemented in derived class) More...
 
void run (const edm::EventSetup &iSetup, const EventInfo &eventInfo) override
 Run the algorithm (must be implemented in derived class) More...
 
void terminate (const edm::EventSetup &iSetup) override
 Call at end of each loop (must be implemented in derived class) More...
 
 ~CSCOverlapsAlignmentAlgorithm () override
 
- Public Member Functions inherited from AlignmentAlgorithmBase
virtual bool addCalibrations (const Calibrations &)
 
bool addCalibrations (const CalibrationsOwner &cals)
 
 AlignmentAlgorithmBase (const edm::ParameterSet &, const edm::ConsumesCollector &)
 Constructor. More...
 
virtual void beginLuminosityBlock (const edm::EventSetup &setup)
 called at begin of luminosity block (no lumi block info passed yet) More...
 
virtual void beginRun (const edm::Run &, const edm::EventSetup &, bool changed)
 called at begin of run More...
 
virtual void endLuminosityBlock (const edm::EventSetup &setup)
 called at end of luminosity block (no lumi block info passed yet) More...
 
virtual void endRun (const EndRunInfo &runInfo, const edm::EventSetup &setup)
 called at end of run - order of arguments like in EDProducer etc. More...
 
virtual bool processesEvents ()
 Returns whether algorithm proccesses events in current configuration. More...
 
virtual bool setParametersForRunRange (const RunRange &rr)
 
virtual void startNewLoop ()
 
virtual bool storeAlignments ()
 Returns whether algorithm produced results to be stored. More...
 
virtual bool supportsCalibrations ()
 
virtual void terminate ()
 Called at end of job (must be implemented in derived class) More...
 
virtual ~AlignmentAlgorithmBase ()
 Destructor. More...
 

Public Attributes

bool m_combineME11
 
TH1F * m_drdz
 
bool m_errorFromRMS
 
bool m_fiducial
 
TProfile * m_fiducial_ME11
 
TProfile * m_fiducial_ME12
 
TProfile * m_fiducial_MEx1
 
TProfile * m_fiducial_MEx2
 
TH1F * m_hitsPerChamber
 
bool m_makeHistograms
 
double m_maxdrdz
 
int m_minHitsPerChamber
 
int m_minStationsInTrackRefits
 
int m_minTracksPerOverlap
 
int m_mode
 
TH2F * m_occupancy
 
TH1F * m_offsetResiduals
 
TH1F * m_offsetResiduals_normalized
 
TH1F * m_offsetResiduals_weighted
 
TH2F * m_RPhipos_mem1
 
TH2F * m_RPhipos_mem2
 
TH2F * m_RPhipos_mem3
 
TH2F * m_RPhipos_mem4
 
TH2F * m_RPhipos_mep1
 
TH2F * m_RPhipos_mep2
 
TH2F * m_RPhipos_mep3
 
TH2F * m_RPhipos_mep4
 
TH1F * m_slope
 
TH1F * m_slope_MEm1
 
TH1F * m_slope_MEm2
 
TH1F * m_slope_MEm3
 
TH1F * m_slope_MEm4
 
TH1F * m_slope_MEp1
 
TH1F * m_slope_MEp2
 
TH1F * m_slope_MEp3
 
TH1F * m_slope_MEp4
 
bool m_slopeFromTrackRefit
 
TH1F * m_slopeResiduals
 
TH1F * m_slopeResiduals_normalized
 
TH1F * m_slopeResiduals_weighted
 
double m_truncateOffsetResid
 
double m_truncateSlopeResid
 
bool m_useHitWeights
 
bool m_useTrackWeights
 
TH2F * m_XYpos_mem1
 
TH2F * m_XYpos_mem2
 
TH2F * m_XYpos_mem3
 
TH2F * m_XYpos_mem4
 
TH2F * m_XYpos_mep1
 
TH2F * m_XYpos_mep2
 
TH2F * m_XYpos_mep3
 
TH2F * m_XYpos_mep4
 

Private Attributes

AlignableNavigatorm_alignableNavigator
 
align::Alignables m_alignables
 
AlignmentParameterStorem_alignmentParameterStore
 
const edm::ESGetToken< CSCGeometry, MuonGeometryRecordm_cscGeometryToken
 
bool m_doAlignment
 
std::vector< CSCChamberFitterm_fitters
 
TH1F * m_histP10
 
TH1F * m_histP100
 
TH1F * m_histP1000
 
double m_maxRedChi2
 
double m_minP
 
std::string m_mode_string
 
const Propagatorm_propagatorPointer
 
const edm::ESGetToken< Propagator, TrackingComponentsRecordm_propToken
 
std::map< std::pair< CSCDetId, CSCDetId >, CSCPairResidualsConstraint * > m_quickChamberLookup
 
std::vector< std::string > m_readTemporaryFiles
 
std::string m_reportFileName
 
std::vector< CSCPairResidualsConstraint * > m_residualsConstraints
 
TrackTransformerm_trackTransformer
 
const edm::ESGetToken< TransientTrackBuilder, TransientTrackRecordm_tthbToken
 
std::string m_writeTemporaryFile
 

Additional Inherited Members

- Public Types inherited from AlignmentAlgorithmBase
typedef std::pair< const Trajectory *, const reco::Track * > ConstTrajTrackPair
 
typedef std::vector< ConstTrajTrackPairConstTrajTrackPairCollection
 
using RunNumber = align::RunNumber
 
using RunRange = align::RunRange
 

Detailed Description

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 61 of file CSCOverlapsAlignmentAlgorithm.h.

Constructor & Destructor Documentation

◆ CSCOverlapsAlignmentAlgorithm()

CSCOverlapsAlignmentAlgorithm::CSCOverlapsAlignmentAlgorithm ( const edm::ParameterSet iConfig,
edm::ConsumesCollector iC 
)

Definition at line 3 of file CSCOverlapsAlignmentAlgorithm.cc.

5  : AlignmentAlgorithmBase(iConfig, iC),
6  m_minHitsPerChamber(iConfig.getParameter<int>("minHitsPerChamber")),
7  m_maxdrdz(iConfig.getParameter<double>("maxdrdz")),
8  m_fiducial(iConfig.getParameter<bool>("fiducial")),
9  m_useHitWeights(iConfig.getParameter<bool>("useHitWeights")),
10  m_slopeFromTrackRefit(iConfig.getParameter<bool>("slopeFromTrackRefit")),
11  m_minStationsInTrackRefits(iConfig.getParameter<int>("minStationsInTrackRefits")),
12  m_truncateSlopeResid(iConfig.getParameter<double>("truncateSlopeResid")),
13  m_truncateOffsetResid(iConfig.getParameter<double>("truncateOffsetResid")),
14  m_combineME11(iConfig.getParameter<bool>("combineME11")),
15  m_useTrackWeights(iConfig.getParameter<bool>("useTrackWeights")),
16  m_errorFromRMS(iConfig.getParameter<bool>("errorFromRMS")),
17  m_minTracksPerOverlap(iConfig.getParameter<int>("minTracksPerOverlap")),
18  m_makeHistograms(iConfig.getParameter<bool>("makeHistograms")),
21  "",
23  ? iConfig.getParameter<edm::ParameterSet>("TrackTransformer").getParameter<std::string>("Propagator")
24  : std::string("")))),
25  m_tthbToken(iC.esConsumes(edm::ESInputTag("", "TransientTrackBuilder"))),
26  m_mode_string(iConfig.getParameter<std::string>("mode")),
27  m_reportFileName(iConfig.getParameter<std::string>("reportFileName")),
28  m_minP(iConfig.getParameter<double>("minP")),
29  m_maxRedChi2(iConfig.getParameter<double>("maxRedChi2")),
30  m_writeTemporaryFile(iConfig.getParameter<std::string>("writeTemporaryFile")),
31  m_readTemporaryFiles(iConfig.getParameter<std::vector<std::string> >("readTemporaryFiles")),
32  m_doAlignment(iConfig.getParameter<bool>("doAlignment")) {
33  if (m_mode_string == std::string("phiy"))
35  else if (m_mode_string == std::string("phipos"))
37  else if (m_mode_string == std::string("phiz"))
39  else if (m_mode_string == std::string("radius"))
41  else
42  throw cms::Exception("BadConfig") << "mode must be one of \"phiy\", \"phipos\", \"phiz\", \"radius\"" << std::endl;
43 
44  std::vector<edm::ParameterSet> fitters = iConfig.getParameter<std::vector<edm::ParameterSet> >("fitters");
45  for (std::vector<edm::ParameterSet>::const_iterator fitter = fitters.begin(); fitter != fitters.end(); ++fitter) {
47  }
48 
49  for (std::vector<CSCPairResidualsConstraint*>::const_iterator residualsConstraint = m_residualsConstraints.begin();
50  residualsConstraint != m_residualsConstraints.end();
51  ++residualsConstraint) {
52  (*residualsConstraint)->configure(this);
53  m_quickChamberLookup[std::pair<CSCDetId, CSCDetId>((*residualsConstraint)->id_i(), (*residualsConstraint)->id_j())] =
54  *residualsConstraint;
55  }
56 
58  m_trackTransformer = new TrackTransformer(iConfig.getParameter<edm::ParameterSet>("TrackTransformer"));
59  } else {
60  m_trackTransformer = nullptr;
61  }
62 
63  m_propagatorPointer = nullptr;
64 
65  if (m_makeHistograms) {
66  edm::Service<TFileService> tFileService;
67  m_histP10 = tFileService->make<TH1F>("P10", "", 100, 0, 10);
68  m_histP100 = tFileService->make<TH1F>("P100", "", 100, 0, 100);
69  m_histP1000 = tFileService->make<TH1F>("P1000", "", 100, 0, 1000);
70 
71  m_hitsPerChamber = tFileService->make<TH1F>("hitsPerChamber", "", 10, -0.5, 9.5);
72 
73  m_fiducial_ME11 = tFileService->make<TProfile>("fiducial_ME11", "", 100, 0.075, 0.100);
74  m_fiducial_ME12 = tFileService->make<TProfile>("fiducial_ME12", "", 100, 0.080, 0.105);
75  m_fiducial_MEx1 = tFileService->make<TProfile>("fiducial_MEx1", "", 100, 0.160, 0.210);
76  m_fiducial_MEx2 = tFileService->make<TProfile>("fiducial_MEx2", "", 100, 0.080, 0.105);
77 
78  m_slope = tFileService->make<TH1F>("slope", "", 100, -0.5, 0.5);
79  m_slope_MEp4 = tFileService->make<TH1F>("slope_MEp4", "", 100, -0.5, 0.5);
80  m_slope_MEp3 = tFileService->make<TH1F>("slope_MEp3", "", 100, -0.5, 0.5);
81  m_slope_MEp2 = tFileService->make<TH1F>("slope_MEp2", "", 100, -0.5, 0.5);
82  m_slope_MEp1 = tFileService->make<TH1F>("slope_MEp1", "", 100, -0.5, 0.5);
83  m_slope_MEm1 = tFileService->make<TH1F>("slope_MEm1", "", 100, -0.5, 0.5);
84  m_slope_MEm2 = tFileService->make<TH1F>("slope_MEm2", "", 100, -0.5, 0.5);
85  m_slope_MEm3 = tFileService->make<TH1F>("slope_MEm3", "", 100, -0.5, 0.5);
86  m_slope_MEm4 = tFileService->make<TH1F>("slope_MEm4", "", 100, -0.5, 0.5);
87 
88  m_slopeResiduals = tFileService->make<TH1F>("slopeResiduals", "mrad", 300, -30., 30.);
89  m_slopeResiduals_weighted = tFileService->make<TH1F>("slopeResiduals_weighted", "mrad", 300, -30., 30.);
90  m_slopeResiduals_normalized = tFileService->make<TH1F>("slopeResiduals_normalized", "", 200, -20., 20.);
91  m_offsetResiduals = tFileService->make<TH1F>("offsetResiduals", "mm", 300, -30., 30.);
92  m_offsetResiduals_weighted = tFileService->make<TH1F>("offsetResiduals_weighted", "mm", 300, -30., 30.);
93  m_offsetResiduals_normalized = tFileService->make<TH1F>("offsetResiduals_normalized", "", 200, -20., 20.);
94 
95  m_drdz = tFileService->make<TH1F>("drdz", "", 100, -0.5, 0.5);
96 
97  m_occupancy = tFileService->make<TH2F>("occupancy", "", 36, 1, 37, 20, 1, 21);
98  for (int i = 1; i <= 36; i++) {
99  std::stringstream pairname;
100  pairname << i << "-";
101  if (i + 1 == 37)
102  pairname << 1;
103  else
104  pairname << (i + 1);
105  m_occupancy->GetXaxis()->SetBinLabel(i, pairname.str().c_str());
106  }
107  m_occupancy->GetYaxis()->SetBinLabel(1, "ME-4/2");
108  m_occupancy->GetYaxis()->SetBinLabel(2, "ME-4/1");
109  m_occupancy->GetYaxis()->SetBinLabel(3, "ME-3/2");
110  m_occupancy->GetYaxis()->SetBinLabel(4, "ME-3/1");
111  m_occupancy->GetYaxis()->SetBinLabel(5, "ME-2/2");
112  m_occupancy->GetYaxis()->SetBinLabel(6, "ME-2/1");
113  m_occupancy->GetYaxis()->SetBinLabel(7, "ME-1/3");
114  m_occupancy->GetYaxis()->SetBinLabel(8, "ME-1/2");
115  if (!m_combineME11) {
116  m_occupancy->GetYaxis()->SetBinLabel(9, "ME-1/1b");
117  m_occupancy->GetYaxis()->SetBinLabel(10, "ME-1/1a");
118  m_occupancy->GetYaxis()->SetBinLabel(11, "ME+1/1a");
119  m_occupancy->GetYaxis()->SetBinLabel(12, "ME+1/1b");
120  } else {
121  m_occupancy->GetYaxis()->SetBinLabel(9, "ME-1/1");
122  m_occupancy->GetYaxis()->SetBinLabel(10, "");
123  m_occupancy->GetYaxis()->SetBinLabel(11, "");
124  m_occupancy->GetYaxis()->SetBinLabel(12, "ME+1/1");
125  }
126  m_occupancy->GetYaxis()->SetBinLabel(13, "ME+1/2");
127  m_occupancy->GetYaxis()->SetBinLabel(14, "ME+1/3");
128  m_occupancy->GetYaxis()->SetBinLabel(15, "ME+2/1");
129  m_occupancy->GetYaxis()->SetBinLabel(16, "ME+2/2");
130  m_occupancy->GetYaxis()->SetBinLabel(17, "ME+3/1");
131  m_occupancy->GetYaxis()->SetBinLabel(18, "ME+3/2");
132  m_occupancy->GetYaxis()->SetBinLabel(19, "ME+4/1");
133  m_occupancy->GetYaxis()->SetBinLabel(20, "ME+4/2");
134 
135  m_XYpos_mep1 = tFileService->make<TH2F>("XYpos_mep1", "Positions: ME+1", 140, -700., 700., 140, -700., 700.);
136  m_XYpos_mep2 = tFileService->make<TH2F>("XYpos_mep2", "Positions: ME+2", 140, -700., 700., 140, -700., 700.);
137  m_XYpos_mep3 = tFileService->make<TH2F>("XYpos_mep3", "Positions: ME+3", 140, -700., 700., 140, -700., 700.);
138  m_XYpos_mep4 = tFileService->make<TH2F>("XYpos_mep4", "Positions: ME+4", 140, -700., 700., 140, -700., 700.);
139  m_XYpos_mem1 = tFileService->make<TH2F>("XYpos_mem1", "Positions: ME-1", 140, -700., 700., 140, -700., 700.);
140  m_XYpos_mem2 = tFileService->make<TH2F>("XYpos_mem2", "Positions: ME-2", 140, -700., 700., 140, -700., 700.);
141  m_XYpos_mem3 = tFileService->make<TH2F>("XYpos_mem3", "Positions: ME-3", 140, -700., 700., 140, -700., 700.);
142  m_XYpos_mem4 = tFileService->make<TH2F>("XYpos_mem4", "Positions: ME-4", 140, -700., 700., 140, -700., 700.);
143  m_RPhipos_mep1 = tFileService->make<TH2F>("RPhipos_mep1", "Positions: ME+1", 144, -M_PI, M_PI, 21, 0., 700.);
144  m_RPhipos_mep2 = tFileService->make<TH2F>("RPhipos_mep2", "Positions: ME+2", 144, -M_PI, M_PI, 21, 0., 700.);
145  m_RPhipos_mep3 = tFileService->make<TH2F>("RPhipos_mep3", "Positions: ME+3", 144, -M_PI, M_PI, 21, 0., 700.);
146  m_RPhipos_mep4 = tFileService->make<TH2F>("RPhipos_mep4", "Positions: ME+4", 144, -M_PI, M_PI, 21, 0., 700.);
147  m_RPhipos_mem1 = tFileService->make<TH2F>("RPhipos_mem1", "Positions: ME-1", 144, -M_PI, M_PI, 21, 0., 700.);
148  m_RPhipos_mem2 = tFileService->make<TH2F>("RPhipos_mem2", "Positions: ME-2", 144, -M_PI, M_PI, 21, 0., 700.);
149  m_RPhipos_mem3 = tFileService->make<TH2F>("RPhipos_mem3", "Positions: ME-3", 144, -M_PI, M_PI, 21, 0., 700.);
150  m_RPhipos_mem4 = tFileService->make<TH2F>("RPhipos_mem4", "Positions: ME-4", 144, -M_PI, M_PI, 21, 0., 700.);
151  } else {
152  m_histP10 = nullptr;
153  m_histP100 = nullptr;
154  m_histP1000 = nullptr;
155  m_hitsPerChamber = nullptr;
156  m_fiducial_ME11 = nullptr;
157  m_fiducial_ME12 = nullptr;
158  m_fiducial_MEx1 = nullptr;
159  m_fiducial_MEx2 = nullptr;
160  m_slope = nullptr;
161  m_slope_MEp4 = nullptr;
162  m_slope_MEp3 = nullptr;
163  m_slope_MEp2 = nullptr;
164  m_slope_MEp1 = nullptr;
165  m_slope_MEm1 = nullptr;
166  m_slope_MEm2 = nullptr;
167  m_slope_MEm3 = nullptr;
168  m_slope_MEm4 = nullptr;
169  m_slopeResiduals = nullptr;
170  m_slopeResiduals_weighted = nullptr;
171  m_slopeResiduals_normalized = nullptr;
172  m_offsetResiduals = nullptr;
173  m_offsetResiduals_weighted = nullptr;
175  m_drdz = nullptr;
176  m_occupancy = nullptr;
177  m_XYpos_mep1 = nullptr;
178  m_XYpos_mep2 = nullptr;
179  m_XYpos_mep3 = nullptr;
180  m_XYpos_mep4 = nullptr;
181  m_XYpos_mem1 = nullptr;
182  m_XYpos_mem2 = nullptr;
183  m_XYpos_mem3 = nullptr;
184  m_XYpos_mem4 = nullptr;
185  m_RPhipos_mep1 = nullptr;
186  m_RPhipos_mep2 = nullptr;
187  m_RPhipos_mep3 = nullptr;
188  m_RPhipos_mep4 = nullptr;
189  m_RPhipos_mem1 = nullptr;
190  m_RPhipos_mem2 = nullptr;
191  m_RPhipos_mem3 = nullptr;
192  m_RPhipos_mem4 = nullptr;
193  }
194 }

References Exception, alignBH_cfg::fitters, edm::ParameterSet::getParameter(), mps_fire::i, CSCPairResidualsConstraint::kModePhiPos, CSCPairResidualsConstraint::kModePhiy, CSCPairResidualsConstraint::kModePhiz, CSCPairResidualsConstraint::kModeRadius, m_combineME11, m_drdz, m_fiducial_ME11, m_fiducial_ME12, m_fiducial_MEx1, m_fiducial_MEx2, m_fitters, m_histP10, m_histP100, m_histP1000, m_hitsPerChamber, m_makeHistograms, m_mode, m_mode_string, m_occupancy, m_offsetResiduals, m_offsetResiduals_normalized, m_offsetResiduals_weighted, M_PI, m_propagatorPointer, m_quickChamberLookup, m_residualsConstraints, m_RPhipos_mem1, m_RPhipos_mem2, m_RPhipos_mem3, m_RPhipos_mem4, m_RPhipos_mep1, m_RPhipos_mep2, m_RPhipos_mep3, m_RPhipos_mep4, m_slope, m_slope_MEm1, m_slope_MEm2, m_slope_MEm3, m_slope_MEm4, m_slope_MEp1, m_slope_MEp2, m_slope_MEp3, m_slope_MEp4, m_slopeFromTrackRefit, m_slopeResiduals, m_slopeResiduals_normalized, m_slopeResiduals_weighted, m_trackTransformer, m_XYpos_mem1, m_XYpos_mem2, m_XYpos_mem3, m_XYpos_mem4, m_XYpos_mep1, m_XYpos_mep2, m_XYpos_mep3, m_XYpos_mep4, TFileService::make(), AlCaHLTBitMon_QueryRunRegistry::string, and HLT_FULL_cff::TrackTransformer.

◆ ~CSCOverlapsAlignmentAlgorithm()

CSCOverlapsAlignmentAlgorithm::~CSCOverlapsAlignmentAlgorithm ( )
override

Definition at line 196 of file CSCOverlapsAlignmentAlgorithm.cc.

196 {}

Member Function Documentation

◆ initialize()

void CSCOverlapsAlignmentAlgorithm::initialize ( const edm::EventSetup setup,
AlignableTracker tracker,
AlignableMuon muon,
AlignableExtras extras,
AlignmentParameterStore store 
)
overridevirtual

Call at beginning of job (must be implemented in derived class)

Implements AlignmentAlgorithmBase.

Definition at line 198 of file CSCOverlapsAlignmentAlgorithm.cc.

202  {
203  m_alignmentParameterStore = alignmentParameterStore;
205 
206  if (alignableTracker == nullptr)
207  m_alignableNavigator = new AlignableNavigator(alignableMuon);
208  else
209  m_alignableNavigator = new AlignableNavigator(alignableTracker, alignableMuon);
210 
211  for (const auto& alignable : m_alignables) {
212  DetId id = alignable->geomDetId();
213  if (id.det() != DetId::Muon || id.subdetId() != MuonSubdetId::CSC || CSCDetId(id.rawId()).layer() != 0) {
214  throw cms::Exception("BadConfig") << "Only CSC chambers may be alignable" << std::endl;
215  }
216 
217  std::vector<bool> selector = alignable->alignmentParameters()->selector();
218  for (std::vector<bool>::const_iterator i = selector.begin(); i != selector.end(); ++i) {
219  if (!(*i))
220  throw cms::Exception("BadConfig") << "All selector strings should be \"111111\"" << std::endl;
221  }
222  }
223 
224  const CSCGeometry* cscGeometry = &iSetup.getData(m_cscGeometryToken);
225 
226  for (std::vector<CSCPairResidualsConstraint*>::const_iterator residualsConstraint = m_residualsConstraints.begin();
227  residualsConstraint != m_residualsConstraints.end();
228  ++residualsConstraint) {
229  (*residualsConstraint)->setZplane(cscGeometry);
230  }
231 
232  if (!m_readTemporaryFiles.empty()) {
233  std::vector<std::ifstream*> input;
234  for (std::vector<std::string>::const_iterator fileName = m_readTemporaryFiles.begin();
236  ++fileName) {
237  input.push_back(new std::ifstream(fileName->c_str()));
238  }
239 
240  for (std::vector<CSCPairResidualsConstraint*>::const_iterator residualsConstraint = m_residualsConstraints.begin();
241  residualsConstraint != m_residualsConstraints.end();
242  ++residualsConstraint) {
243  (*residualsConstraint)->read(input, m_readTemporaryFiles);
244  }
245 
246  for (std::vector<std::ifstream*>::const_iterator file = input.begin(); file != input.end(); ++file) {
247  delete (*file);
248  }
249  }
250 }

References AlignmentParameterStore::alignables(), MuonSubdetId::CSC, Exception, geometryDiff::file, MillePedeFileConverter_cfg::fileName, edm::EventSetup::getData(), mps_fire::i, input, phase1PixelTopology::layer, m_alignableNavigator, m_alignables, m_alignmentParameterStore, m_cscGeometryToken, m_readTemporaryFiles, m_residualsConstraints, and DetId::Muon.

◆ run()

void CSCOverlapsAlignmentAlgorithm::run ( const edm::EventSetup setup,
const EventInfo eventInfo 
)
overridevirtual

Run the algorithm (must be implemented in derived class)

Implements AlignmentAlgorithmBase.

Definition at line 252 of file CSCOverlapsAlignmentAlgorithm.cc.

252  {
254  if (m_slopeFromTrackRefit) {
255  iSetup.getHandle(m_propToken);
256  if (m_propagatorPointer != &*propagator) {
258 
259  for (std::vector<CSCPairResidualsConstraint*>::const_iterator residualsConstraint =
260  m_residualsConstraints.begin();
261  residualsConstraint != m_residualsConstraints.end();
262  ++residualsConstraint) {
263  (*residualsConstraint)->setPropagator(m_propagatorPointer);
264  }
265  }
266  }
267 
268  const TransientTrackBuilder* transientTrackBuilder = &iSetup.getData(m_tthbToken);
269 
270  if (m_trackTransformer != nullptr)
272 
273  const ConstTrajTrackPairCollection& trajtracks = eventInfo.trajTrackPairs();
274  for (ConstTrajTrackPairCollection::const_iterator trajtrack = trajtracks.begin(); trajtrack != trajtracks.end();
275  ++trajtrack) {
276  const Trajectory* traj = (*trajtrack).first;
277  const reco::Track* track = (*trajtrack).second;
278 
279  if (m_makeHistograms) {
280  m_histP10->Fill(track->p());
281  m_histP100->Fill(track->p());
282  m_histP1000->Fill(track->p());
283  }
284  if (track->p() >= m_minP) {
285  std::vector<TrajectoryMeasurement> measurements = traj->measurements();
286  reco::TransientTrack transientTrack = transientTrackBuilder->build(track);
287 
288  std::map<int, std::map<CSCDetId, bool> > stationsToChambers;
289  for (std::vector<TrajectoryMeasurement>::const_iterator measurement = measurements.begin();
290  measurement != measurements.end();
291  ++measurement) {
292  DetId id = measurement->recHit()->geographicalId();
293  if (id.det() == DetId::Muon && id.subdetId() == MuonSubdetId::CSC) {
294  CSCDetId cscid(id.rawId());
295  CSCDetId chamberId(cscid.endcap(), cscid.station(), cscid.ring(), cscid.chamber(), 0);
296  if (m_combineME11 && cscid.station() == 1 && cscid.ring() == 4)
297  chamberId = CSCDetId(cscid.endcap(), 1, 1, cscid.chamber(), 0);
298  int station = (cscid.endcap() == 1 ? 1 : -1) * cscid.station();
299 
300  if (stationsToChambers.find(station) == stationsToChambers.end())
301  stationsToChambers[station] = std::map<CSCDetId, bool>();
302  stationsToChambers[station][chamberId] = true;
303 
304  if (m_makeHistograms) {
305  GlobalPoint pos = measurement->recHit()->globalPosition();
306  if (cscid.endcap() == 1 && cscid.station() == 1) {
307  m_XYpos_mep1->Fill(pos.x(), pos.y());
308  m_RPhipos_mep1->Fill(pos.phi(), pos.perp());
309  }
310  if (cscid.endcap() == 1 && cscid.station() == 2) {
311  m_XYpos_mep2->Fill(pos.x(), pos.y());
312  m_RPhipos_mep2->Fill(pos.phi(), pos.perp());
313  }
314  if (cscid.endcap() == 1 && cscid.station() == 3) {
315  m_XYpos_mep3->Fill(pos.x(), pos.y());
316  m_RPhipos_mep3->Fill(pos.phi(), pos.perp());
317  }
318  if (cscid.endcap() == 1 && cscid.station() == 4) {
319  m_XYpos_mep4->Fill(pos.x(), pos.y());
320  m_RPhipos_mep4->Fill(pos.phi(), pos.perp());
321  }
322  if (cscid.endcap() == 2 && cscid.station() == 1) {
323  m_XYpos_mem1->Fill(pos.x(), pos.y());
324  m_RPhipos_mem1->Fill(pos.phi(), pos.perp());
325  }
326  if (cscid.endcap() == 2 && cscid.station() == 2) {
327  m_XYpos_mem2->Fill(pos.x(), pos.y());
328  m_RPhipos_mem2->Fill(pos.phi(), pos.perp());
329  }
330  if (cscid.endcap() == 2 && cscid.station() == 3) {
331  m_XYpos_mem3->Fill(pos.x(), pos.y());
332  m_RPhipos_mem3->Fill(pos.phi(), pos.perp());
333  }
334  if (cscid.endcap() == 2 && cscid.station() == 4) {
335  m_XYpos_mem4->Fill(pos.x(), pos.y());
336  m_RPhipos_mem4->Fill(pos.phi(), pos.perp());
337  }
338  }
339  }
340  }
341 
342  std::map<CSCPairResidualsConstraint*, bool> residualsConstraints;
343  for (std::map<int, std::map<CSCDetId, bool> >::const_iterator iter = stationsToChambers.begin();
344  iter != stationsToChambers.end();
345  ++iter) {
346  for (std::map<CSCDetId, bool>::const_iterator one = iter->second.begin(); one != iter->second.end(); ++one) {
347  for (std::map<CSCDetId, bool>::const_iterator two = one; two != iter->second.end(); ++two) {
348  if (one != two) {
349  std::map<std::pair<CSCDetId, CSCDetId>, CSCPairResidualsConstraint*>::const_iterator quick;
350 
351  quick = m_quickChamberLookup.find(std::pair<CSCDetId, CSCDetId>(one->first, two->first));
352  if (quick != m_quickChamberLookup.end())
353  residualsConstraints[quick->second] = true;
354 
355  quick = m_quickChamberLookup.find(std::pair<CSCDetId, CSCDetId>(two->first, one->first));
356  if (quick != m_quickChamberLookup.end())
357  residualsConstraints[quick->second] = true;
358  }
359  }
360  }
361  }
362 
363  for (std::map<CSCPairResidualsConstraint*, bool>::const_iterator residualsConstraint =
364  residualsConstraints.begin();
365  residualsConstraint != residualsConstraints.end();
366  ++residualsConstraint) {
367  residualsConstraint->first->addTrack(measurements, transientTrack, m_trackTransformer);
368  }
369  }
370  }
371 }

References TransientTrackBuilder::build(), CSCDetId::chamber(), MuonSubdetId::CSC, CSCDetId::endcap(), dumpTauVariables_cfi::eventInfo, edm::EventSetup::getData(), edm::EventSetup::getHandle(), m_combineME11, m_histP10, m_histP100, m_histP1000, m_makeHistograms, m_minP, m_propagatorPointer, m_propToken, m_quickChamberLookup, m_residualsConstraints, m_RPhipos_mem1, m_RPhipos_mem2, m_RPhipos_mem3, m_RPhipos_mem4, m_RPhipos_mep1, m_RPhipos_mep2, m_RPhipos_mep3, m_RPhipos_mep4, m_slopeFromTrackRefit, m_trackTransformer, m_tthbToken, m_XYpos_mem1, m_XYpos_mem2, m_XYpos_mem3, m_XYpos_mem4, m_XYpos_mep1, m_XYpos_mep2, m_XYpos_mep3, m_XYpos_mep4, genParticles_cff::map, Trajectory::measurements(), DetId::Muon, SiPixelPI::one, TrackCandidateProducer_cfi::propagator, CSCDetId::ring(), TrackTransformer::setServices(), relativeConstraints::station, CSCDetId::station(), HLT_FULL_cff::track, and SiPixelPI::two.

Referenced by DTWorkflow.DTWorkflow::all(), Types.EventID::cppID(), Types.LuminosityBlockID::cppID(), and o2olib.O2OTool::execute().

◆ terminate()

void CSCOverlapsAlignmentAlgorithm::terminate ( const edm::EventSetup iSetup)
overridevirtual

Call at end of each loop (must be implemented in derived class)

Implements AlignmentAlgorithmBase.

Definition at line 373 of file CSCOverlapsAlignmentAlgorithm.cc.

373  {
374  // write residuals partial fits to temporary files for collection
375  if (m_writeTemporaryFile != std::string("")) {
376  std::ofstream output(m_writeTemporaryFile.c_str());
377  for (std::vector<CSCPairResidualsConstraint*>::const_iterator residualsConstraint = m_residualsConstraints.begin();
378  residualsConstraint != m_residualsConstraints.end();
379  ++residualsConstraint) {
380  (*residualsConstraint)->write(output);
381  }
382  }
383 
384  // write report for alignment results
385  if (m_doAlignment) {
386  std::ofstream report;
387  bool writeReport = (m_reportFileName != std::string(""));
388  if (writeReport) {
389  report.open(m_reportFileName.c_str());
390  report << "cscReports = []" << std::endl
391  << std::endl
392  << "class CSCChamberCorrection:" << std::endl
393  << " def __init__(self, name, detid, value):" << std::endl
394  << " self.name, self.detid, self.value = name, detid, value" << std::endl
395  << std::endl
396  << "class CSCErrorMode:" << std::endl
397  << " def __init__(self, error):" << std::endl
398  << " self.error = error" << std::endl
399  << " self.terms = {}" << std::endl
400  << " self.detids = {}" << std::endl
401  << " def addTerm(self, name, detid, coefficient):" << std::endl
402  << " self.terms[name] = coefficient" << std::endl
403  << " self.detids[name] = detid" << std::endl
404  << std::endl
405  << "class CSCConstraintResidual:" << std::endl
406  << " def __init__(self, i, j, before, uncert, residual, pull):" << std::endl
407  << " self.i, self.j, self.before, self.error, self.residual, self.pull = i, j, before, uncert, "
408  "residual, pull"
409  << std::endl
410  << std::endl
411  << "class CSCFitterReport:" << std::endl
412  << " def __init__(self, name, oldchi2, newchi2):" << std::endl
413  << " self.name, self.oldchi2, self.newchi2 = name, oldchi2, newchi2" << std::endl
414  << " self.chamberCorrections = []" << std::endl
415  << " self.errorModes = []" << std::endl
416  << " self.constraintResiduals = []" << std::endl
417  << std::endl
418  << " def addChamberCorrection(self, name, detid, value):" << std::endl
419  << " self.chamberCorrections.append(CSCChamberCorrection(name, detid, value))" << std::endl
420  << std::endl
421  << " def addErrorMode(self, error):" << std::endl
422  << " self.errorModes.append(CSCErrorMode(error))" << std::endl
423  << std::endl
424  << " def addErrorModeTerm(self, name, detid, coefficient):" << std::endl
425  << " self.errorModes[-1].addTerm(name, detid, coefficient)" << std::endl
426  << std::endl
427  << " def addCSCConstraintResidual(self, i, j, before, uncert, residual, pull):" << std::endl
428  << " self.constraintResiduals.append(CSCConstraintResidual(i, j, before, uncert, residual, pull))"
429  << std::endl
430  << std::endl
431  << "import re" << std::endl
432  << "def nameToKey(name):" << std::endl
433  << " match = re.match(\"ME([\\+\\-])([1-4])/([1-4])/([0-9]{2})\", name)" << std::endl
434  << " if match is None: return None" << std::endl
435  << " endcap, station, ring, chamber = match.groups()" << std::endl
436  << " if endcap == \"+\": endcap = 1" << std::endl
437  << " else: endcap = 2" << std::endl
438  << " station = int(station)" << std::endl
439  << " ring = int(ring)" << std::endl
440  << " chamber = int(chamber)" << std::endl
441  << " return endcap, station, ring, chamber" << std::endl
442  << std::endl;
443  }
444 
445  for (std::vector<CSCChamberFitter>::const_iterator fitter = m_fitters.begin(); fitter != m_fitters.end();
446  ++fitter) {
449 
450  } else {
451  std::vector<CSCAlignmentCorrections*> corrections;
452  fitter->fit(corrections);
453 
454  // corrections only exist if the fit was successful
455  for (std::vector<CSCAlignmentCorrections*>::iterator correction = corrections.begin();
456  correction != corrections.end();
457  ++correction) {
458  (*correction)->applyAlignment(m_alignableNavigator, m_alignmentParameterStore, m_mode, m_combineME11);
459  if (m_makeHistograms)
460  (*correction)->plot();
461  if (writeReport)
462  (*correction)->report(report);
463  }
464  }
465  }
466  }
467 }

References pfMETCorrectionType0_cfi::correction, CSCPairResidualsConstraint::kModeRadius, m_alignableNavigator, m_alignmentParameterStore, m_combineME11, m_doAlignment, m_fitters, m_makeHistograms, m_mode, m_reportFileName, m_residualsConstraints, m_writeTemporaryFile, convertSQLitetoXML_cfg::output, edmIntegrityCheck::report, and AlCaHLTBitMon_QueryRunRegistry::string.

Member Data Documentation

◆ m_alignableNavigator

AlignableNavigator* CSCOverlapsAlignmentAlgorithm::m_alignableNavigator
private

Definition at line 150 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by initialize(), and terminate().

◆ m_alignables

align::Alignables CSCOverlapsAlignmentAlgorithm::m_alignables
private

Definition at line 149 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by initialize().

◆ m_alignmentParameterStore

AlignmentParameterStore* CSCOverlapsAlignmentAlgorithm::m_alignmentParameterStore
private

Definition at line 148 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by initialize(), and terminate().

◆ m_combineME11

bool CSCOverlapsAlignmentAlgorithm::m_combineME11

◆ m_cscGeometryToken

const edm::ESGetToken<CSCGeometry, MuonGeometryRecord> CSCOverlapsAlignmentAlgorithm::m_cscGeometryToken
private

Definition at line 136 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by initialize().

◆ m_doAlignment

bool CSCOverlapsAlignmentAlgorithm::m_doAlignment
private

Definition at line 146 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by terminate().

◆ m_drdz

TH1F* CSCOverlapsAlignmentAlgorithm::m_drdz

◆ m_errorFromRMS

bool CSCOverlapsAlignmentAlgorithm::m_errorFromRMS

Definition at line 131 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCPairResidualsConstraint::error().

◆ m_fiducial

bool CSCOverlapsAlignmentAlgorithm::m_fiducial

◆ m_fiducial_ME11

TProfile* CSCOverlapsAlignmentAlgorithm::m_fiducial_ME11

◆ m_fiducial_ME12

TProfile* CSCOverlapsAlignmentAlgorithm::m_fiducial_ME12

◆ m_fiducial_MEx1

TProfile* CSCOverlapsAlignmentAlgorithm::m_fiducial_MEx1

◆ m_fiducial_MEx2

TProfile* CSCOverlapsAlignmentAlgorithm::m_fiducial_MEx2

◆ m_fitters

std::vector<CSCChamberFitter> CSCOverlapsAlignmentAlgorithm::m_fitters
private

Definition at line 151 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm(), and terminate().

◆ m_histP10

TH1F* CSCOverlapsAlignmentAlgorithm::m_histP10
private

Definition at line 158 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm(), and run().

◆ m_histP100

TH1F* CSCOverlapsAlignmentAlgorithm::m_histP100
private

Definition at line 159 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm(), and run().

◆ m_histP1000

TH1F* CSCOverlapsAlignmentAlgorithm::m_histP1000
private

Definition at line 160 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm(), and run().

◆ m_hitsPerChamber

TH1F* CSCOverlapsAlignmentAlgorithm::m_hitsPerChamber

◆ m_makeHistograms

bool CSCOverlapsAlignmentAlgorithm::m_makeHistograms

◆ m_maxdrdz

double CSCOverlapsAlignmentAlgorithm::m_maxdrdz

◆ m_maxRedChi2

double CSCOverlapsAlignmentAlgorithm::m_maxRedChi2
private

Definition at line 143 of file CSCOverlapsAlignmentAlgorithm.h.

◆ m_minHitsPerChamber

int CSCOverlapsAlignmentAlgorithm::m_minHitsPerChamber

◆ m_minP

double CSCOverlapsAlignmentAlgorithm::m_minP
private

Definition at line 142 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by run().

◆ m_minStationsInTrackRefits

int CSCOverlapsAlignmentAlgorithm::m_minStationsInTrackRefits

◆ m_minTracksPerOverlap

int CSCOverlapsAlignmentAlgorithm::m_minTracksPerOverlap

Definition at line 132 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCPairResidualsConstraint::valid().

◆ m_mode

int CSCOverlapsAlignmentAlgorithm::m_mode

◆ m_mode_string

std::string CSCOverlapsAlignmentAlgorithm::m_mode_string
private

Definition at line 140 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm().

◆ m_occupancy

TH2F* CSCOverlapsAlignmentAlgorithm::m_occupancy

◆ m_offsetResiduals

TH1F* CSCOverlapsAlignmentAlgorithm::m_offsetResiduals

◆ m_offsetResiduals_normalized

TH1F* CSCOverlapsAlignmentAlgorithm::m_offsetResiduals_normalized

◆ m_offsetResiduals_weighted

TH1F* CSCOverlapsAlignmentAlgorithm::m_offsetResiduals_weighted

◆ m_propagatorPointer

const Propagator* CSCOverlapsAlignmentAlgorithm::m_propagatorPointer
private

Definition at line 156 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm(), and run().

◆ m_propToken

const edm::ESGetToken<Propagator, TrackingComponentsRecord> CSCOverlapsAlignmentAlgorithm::m_propToken
private

Definition at line 137 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by run().

◆ m_quickChamberLookup

std::map<std::pair<CSCDetId, CSCDetId>, CSCPairResidualsConstraint *> CSCOverlapsAlignmentAlgorithm::m_quickChamberLookup
private

Definition at line 153 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm(), and run().

◆ m_readTemporaryFiles

std::vector<std::string> CSCOverlapsAlignmentAlgorithm::m_readTemporaryFiles
private

Definition at line 145 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by initialize().

◆ m_reportFileName

std::string CSCOverlapsAlignmentAlgorithm::m_reportFileName
private

Definition at line 141 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by terminate().

◆ m_residualsConstraints

std::vector<CSCPairResidualsConstraint *> CSCOverlapsAlignmentAlgorithm::m_residualsConstraints
private

◆ m_RPhipos_mem1

TH2F* CSCOverlapsAlignmentAlgorithm::m_RPhipos_mem1

Definition at line 115 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm(), and run().

◆ m_RPhipos_mem2

TH2F* CSCOverlapsAlignmentAlgorithm::m_RPhipos_mem2

Definition at line 116 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm(), and run().

◆ m_RPhipos_mem3

TH2F* CSCOverlapsAlignmentAlgorithm::m_RPhipos_mem3

Definition at line 117 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm(), and run().

◆ m_RPhipos_mem4

TH2F* CSCOverlapsAlignmentAlgorithm::m_RPhipos_mem4

Definition at line 118 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm(), and run().

◆ m_RPhipos_mep1

TH2F* CSCOverlapsAlignmentAlgorithm::m_RPhipos_mep1

Definition at line 111 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm(), and run().

◆ m_RPhipos_mep2

TH2F* CSCOverlapsAlignmentAlgorithm::m_RPhipos_mep2

Definition at line 112 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm(), and run().

◆ m_RPhipos_mep3

TH2F* CSCOverlapsAlignmentAlgorithm::m_RPhipos_mep3

Definition at line 113 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm(), and run().

◆ m_RPhipos_mep4

TH2F* CSCOverlapsAlignmentAlgorithm::m_RPhipos_mep4

Definition at line 114 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm(), and run().

◆ m_slope

TH1F* CSCOverlapsAlignmentAlgorithm::m_slope

◆ m_slope_MEm1

TH1F* CSCOverlapsAlignmentAlgorithm::m_slope_MEm1

◆ m_slope_MEm2

TH1F* CSCOverlapsAlignmentAlgorithm::m_slope_MEm2

◆ m_slope_MEm3

TH1F* CSCOverlapsAlignmentAlgorithm::m_slope_MEm3

◆ m_slope_MEm4

TH1F* CSCOverlapsAlignmentAlgorithm::m_slope_MEm4

◆ m_slope_MEp1

TH1F* CSCOverlapsAlignmentAlgorithm::m_slope_MEp1

◆ m_slope_MEp2

TH1F* CSCOverlapsAlignmentAlgorithm::m_slope_MEp2

◆ m_slope_MEp3

TH1F* CSCOverlapsAlignmentAlgorithm::m_slope_MEp3

◆ m_slope_MEp4

TH1F* CSCOverlapsAlignmentAlgorithm::m_slope_MEp4

◆ m_slopeFromTrackRefit

bool CSCOverlapsAlignmentAlgorithm::m_slopeFromTrackRefit

◆ m_slopeResiduals

TH1F* CSCOverlapsAlignmentAlgorithm::m_slopeResiduals

◆ m_slopeResiduals_normalized

TH1F* CSCOverlapsAlignmentAlgorithm::m_slopeResiduals_normalized

◆ m_slopeResiduals_weighted

TH1F* CSCOverlapsAlignmentAlgorithm::m_slopeResiduals_weighted

◆ m_trackTransformer

TrackTransformer* CSCOverlapsAlignmentAlgorithm::m_trackTransformer
private

Definition at line 155 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm(), and run().

◆ m_truncateOffsetResid

double CSCOverlapsAlignmentAlgorithm::m_truncateOffsetResid

◆ m_truncateSlopeResid

double CSCOverlapsAlignmentAlgorithm::m_truncateSlopeResid

◆ m_tthbToken

const edm::ESGetToken<TransientTrackBuilder, TransientTrackRecord> CSCOverlapsAlignmentAlgorithm::m_tthbToken
private

Definition at line 138 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by run().

◆ m_useHitWeights

bool CSCOverlapsAlignmentAlgorithm::m_useHitWeights

◆ m_useTrackWeights

bool CSCOverlapsAlignmentAlgorithm::m_useTrackWeights

◆ m_writeTemporaryFile

std::string CSCOverlapsAlignmentAlgorithm::m_writeTemporaryFile
private

Definition at line 144 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by terminate().

◆ m_XYpos_mem1

TH2F* CSCOverlapsAlignmentAlgorithm::m_XYpos_mem1

Definition at line 107 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm(), and run().

◆ m_XYpos_mem2

TH2F* CSCOverlapsAlignmentAlgorithm::m_XYpos_mem2

Definition at line 108 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm(), and run().

◆ m_XYpos_mem3

TH2F* CSCOverlapsAlignmentAlgorithm::m_XYpos_mem3

Definition at line 109 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm(), and run().

◆ m_XYpos_mem4

TH2F* CSCOverlapsAlignmentAlgorithm::m_XYpos_mem4

Definition at line 110 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm(), and run().

◆ m_XYpos_mep1

TH2F* CSCOverlapsAlignmentAlgorithm::m_XYpos_mep1

Definition at line 103 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm(), and run().

◆ m_XYpos_mep2

TH2F* CSCOverlapsAlignmentAlgorithm::m_XYpos_mep2

Definition at line 104 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm(), and run().

◆ m_XYpos_mep3

TH2F* CSCOverlapsAlignmentAlgorithm::m_XYpos_mep3

Definition at line 105 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm(), and run().

◆ m_XYpos_mep4

TH2F* CSCOverlapsAlignmentAlgorithm::m_XYpos_mep4

Definition at line 106 of file CSCOverlapsAlignmentAlgorithm.h.

Referenced by CSCOverlapsAlignmentAlgorithm(), and run().

CSCOverlapsAlignmentAlgorithm::m_slope_MEm4
TH1F * m_slope_MEm4
Definition: CSCOverlapsAlignmentAlgorithm.h:91
alignBH_cfg.fitters
fitters
Definition: alignBH_cfg.py:41
CSCOverlapsAlignmentAlgorithm::m_alignables
align::Alignables m_alignables
Definition: CSCOverlapsAlignmentAlgorithm.h:149
CSCOverlapsAlignmentAlgorithm::m_useTrackWeights
bool m_useTrackWeights
Definition: CSCOverlapsAlignmentAlgorithm.h:130
MuonSubdetId::CSC
static constexpr int CSC
Definition: MuonSubdetId.h:12
CSCOverlapsAlignmentAlgorithm::m_XYpos_mep2
TH2F * m_XYpos_mep2
Definition: CSCOverlapsAlignmentAlgorithm.h:104
CSCOverlapsAlignmentAlgorithm::m_offsetResiduals
TH1F * m_offsetResiduals
Definition: CSCOverlapsAlignmentAlgorithm.h:96
mps_fire.i
i
Definition: mps_fire.py:428
SiPixelPI::one
Definition: SiPixelPayloadInspectorHelper.h:39
edm::ESInputTag
Definition: ESInputTag.h:87
AlignmentAlgorithmBase::ConstTrajTrackPairCollection
std::vector< ConstTrajTrackPair > ConstTrajTrackPairCollection
Definition: AlignmentAlgorithmBase.h:68
HLT_FULL_cff.track
track
Definition: HLT_FULL_cff.py:11724
input
static const std::string input
Definition: EdmProvDump.cc:48
CSCOverlapsAlignmentAlgorithm::m_trackTransformer
TrackTransformer * m_trackTransformer
Definition: CSCOverlapsAlignmentAlgorithm.h:155
CSCOverlapsAlignmentAlgorithm::m_histP10
TH1F * m_histP10
Definition: CSCOverlapsAlignmentAlgorithm.h:158
CSCOverlapsAlignmentAlgorithm::m_reportFileName
std::string m_reportFileName
Definition: CSCOverlapsAlignmentAlgorithm.h:141
AlignmentParameterStore::alignables
const align::Alignables & alignables(void) const
get all alignables
Definition: AlignmentParameterStore.h:47
CSCOverlapsAlignmentAlgorithm::m_fitters
std::vector< CSCChamberFitter > m_fitters
Definition: CSCOverlapsAlignmentAlgorithm.h:151
CSCOverlapsAlignmentAlgorithm::m_fiducial_ME12
TProfile * m_fiducial_ME12
Definition: CSCOverlapsAlignmentAlgorithm.h:79
CSCOverlapsAlignmentAlgorithm::m_fiducial_ME11
TProfile * m_fiducial_ME11
Definition: CSCOverlapsAlignmentAlgorithm.h:78
TrackTransformer::setServices
void setServices(const edm::EventSetup &) override
set the services needed by the TrackTransformer
Definition: TrackTransformer.cc:79
convertSQLitetoXML_cfg.output
output
Definition: convertSQLitetoXML_cfg.py:72
relativeConstraints.station
station
Definition: relativeConstraints.py:67
CSCOverlapsAlignmentAlgorithm::m_maxdrdz
double m_maxdrdz
Definition: CSCOverlapsAlignmentAlgorithm.h:122
CSCOverlapsAlignmentAlgorithm::m_readTemporaryFiles
std::vector< std::string > m_readTemporaryFiles
Definition: CSCOverlapsAlignmentAlgorithm.h:145
CSCOverlapsAlignmentAlgorithm::m_propToken
const edm::ESGetToken< Propagator, TrackingComponentsRecord > m_propToken
Definition: CSCOverlapsAlignmentAlgorithm.h:137
CSCOverlapsAlignmentAlgorithm::m_truncateSlopeResid
double m_truncateSlopeResid
Definition: CSCOverlapsAlignmentAlgorithm.h:127
pos
Definition: PixelAliasList.h:18
CSCOverlapsAlignmentAlgorithm::m_RPhipos_mep1
TH2F * m_RPhipos_mep1
Definition: CSCOverlapsAlignmentAlgorithm.h:111
CSCOverlapsAlignmentAlgorithm::m_histP1000
TH1F * m_histP1000
Definition: CSCOverlapsAlignmentAlgorithm.h:160
CSCOverlapsAlignmentAlgorithm::m_occupancy
TH2F * m_occupancy
Definition: CSCOverlapsAlignmentAlgorithm.h:102
CSCOverlapsAlignmentAlgorithm::m_makeHistograms
bool m_makeHistograms
Definition: CSCOverlapsAlignmentAlgorithm.h:133
CSCOverlapsAlignmentAlgorithm::m_quickChamberLookup
std::map< std::pair< CSCDetId, CSCDetId >, CSCPairResidualsConstraint * > m_quickChamberLookup
Definition: CSCOverlapsAlignmentAlgorithm.h:153
MillePedeFileConverter_cfg.fileName
fileName
Definition: MillePedeFileConverter_cfg.py:32
CSCOverlapsAlignmentAlgorithm::m_XYpos_mep3
TH2F * m_XYpos_mep3
Definition: CSCOverlapsAlignmentAlgorithm.h:105
CSCOverlapsAlignmentAlgorithm::m_writeTemporaryFile
std::string m_writeTemporaryFile
Definition: CSCOverlapsAlignmentAlgorithm.h:144
edm::ConsumesCollector::esConsumes
auto esConsumes()
Definition: ConsumesCollector.h:97
CSCOverlapsAlignmentAlgorithm::m_minStationsInTrackRefits
int m_minStationsInTrackRefits
Definition: CSCOverlapsAlignmentAlgorithm.h:126
CSCOverlapsAlignmentAlgorithm::m_slope_MEm1
TH1F * m_slope_MEm1
Definition: CSCOverlapsAlignmentAlgorithm.h:88
CSCOverlapsAlignmentAlgorithm::m_truncateOffsetResid
double m_truncateOffsetResid
Definition: CSCOverlapsAlignmentAlgorithm.h:128
CSCOverlapsAlignmentAlgorithm::m_alignmentParameterStore
AlignmentParameterStore * m_alignmentParameterStore
Definition: CSCOverlapsAlignmentAlgorithm.h:148
CSCGeometry
Definition: CSCGeometry.h:24
DetId
Definition: DetId.h:17
CSCOverlapsAlignmentAlgorithm::m_fiducial_MEx2
TProfile * m_fiducial_MEx2
Definition: CSCOverlapsAlignmentAlgorithm.h:81
CSCOverlapsAlignmentAlgorithm::m_slopeResiduals
TH1F * m_slopeResiduals
Definition: CSCOverlapsAlignmentAlgorithm.h:93
CSCOverlapsAlignmentAlgorithm::m_drdz
TH1F * m_drdz
Definition: CSCOverlapsAlignmentAlgorithm.h:100
CSCOverlapsAlignmentAlgorithm::m_cscGeometryToken
const edm::ESGetToken< CSCGeometry, MuonGeometryRecord > m_cscGeometryToken
Definition: CSCOverlapsAlignmentAlgorithm.h:136
CSCOverlapsAlignmentAlgorithm::m_slope_MEp3
TH1F * m_slope_MEp3
Definition: CSCOverlapsAlignmentAlgorithm.h:85
TrackCandidateProducer_cfi.propagator
propagator
Definition: TrackCandidateProducer_cfi.py:17
CSCOverlapsAlignmentAlgorithm::m_slopeResiduals_weighted
TH1F * m_slopeResiduals_weighted
Definition: CSCOverlapsAlignmentAlgorithm.h:94
CSCOverlapsAlignmentAlgorithm::m_mode_string
std::string m_mode_string
Definition: CSCOverlapsAlignmentAlgorithm.h:140
CSCOverlapsAlignmentAlgorithm::m_RPhipos_mem4
TH2F * m_RPhipos_mem4
Definition: CSCOverlapsAlignmentAlgorithm.h:118
CSCOverlapsAlignmentAlgorithm::m_minP
double m_minP
Definition: CSCOverlapsAlignmentAlgorithm.h:142
reco::Track
Definition: Track.h:27
edm::ESHandle< Propagator >
CSCOverlapsAlignmentAlgorithm::m_propagatorPointer
const Propagator * m_propagatorPointer
Definition: CSCOverlapsAlignmentAlgorithm.h:156
CSCOverlapsAlignmentAlgorithm::m_offsetResiduals_weighted
TH1F * m_offsetResiduals_weighted
Definition: CSCOverlapsAlignmentAlgorithm.h:97
CSCPairResidualsConstraint::kModePhiy
Definition: CSCPairResidualsConstraint.h:50
CSCPairResidualsConstraint
Definition: CSCPairResidualsConstraint.h:27
Point3DBase< float, GlobalTag >
AlignableNavigator
Definition: AlignableNavigator.h:25
HLT_FULL_cff.TrackTransformer
TrackTransformer
Definition: HLT_FULL_cff.py:9725
phase1PixelTopology::layer
constexpr std::array< uint8_t, layerIndexSize > layer
Definition: phase1PixelTopology.h:99
CSCOverlapsAlignmentAlgorithm::m_tthbToken
const edm::ESGetToken< TransientTrackBuilder, TransientTrackRecord > m_tthbToken
Definition: CSCOverlapsAlignmentAlgorithm.h:138
CSCOverlapsAlignmentAlgorithm::m_histP100
TH1F * m_histP100
Definition: CSCOverlapsAlignmentAlgorithm.h:159
geometryDiff.file
file
Definition: geometryDiff.py:13
CSCOverlapsAlignmentAlgorithm::m_offsetResiduals_normalized
TH1F * m_offsetResiduals_normalized
Definition: CSCOverlapsAlignmentAlgorithm.h:98
CSCOverlapsAlignmentAlgorithm::m_slope_MEp1
TH1F * m_slope_MEp1
Definition: CSCOverlapsAlignmentAlgorithm.h:87
CSCOverlapsAlignmentAlgorithm::m_minHitsPerChamber
int m_minHitsPerChamber
Definition: CSCOverlapsAlignmentAlgorithm.h:121
edm::ParameterSet
Definition: ParameterSet.h:47
CSCOverlapsAlignmentAlgorithm::m_XYpos_mep4
TH2F * m_XYpos_mep4
Definition: CSCOverlapsAlignmentAlgorithm.h:106
CSCOverlapsAlignmentAlgorithm::m_RPhipos_mep4
TH2F * m_RPhipos_mep4
Definition: CSCOverlapsAlignmentAlgorithm.h:114
CSCOverlapsAlignmentAlgorithm::m_RPhipos_mem3
TH2F * m_RPhipos_mem3
Definition: CSCOverlapsAlignmentAlgorithm.h:117
CSCPairResidualsConstraint::kModeRadius
Definition: CSCPairResidualsConstraint.h:50
CSCPairResidualsConstraint::kModePhiz
Definition: CSCPairResidualsConstraint.h:50
pfMETCorrectionType0_cfi.correction
correction
Definition: pfMETCorrectionType0_cfi.py:39
CSCOverlapsAlignmentAlgorithm::m_slope_MEm3
TH1F * m_slope_MEm3
Definition: CSCOverlapsAlignmentAlgorithm.h:90
CSCOverlapsAlignmentAlgorithm::m_errorFromRMS
bool m_errorFromRMS
Definition: CSCOverlapsAlignmentAlgorithm.h:131
CSCDetId
Definition: CSCDetId.h:26
CSCOverlapsAlignmentAlgorithm::m_mode
int m_mode
Definition: CSCOverlapsAlignmentAlgorithm.h:120
edm::Service< TFileService >
CSCPairResidualsConstraint::kModePhiPos
Definition: CSCPairResidualsConstraint.h:50
CSCOverlapsAlignmentAlgorithm::m_useHitWeights
bool m_useHitWeights
Definition: CSCOverlapsAlignmentAlgorithm.h:124
CSCOverlapsAlignmentAlgorithm::m_XYpos_mep1
TH2F * m_XYpos_mep1
Definition: CSCOverlapsAlignmentAlgorithm.h:103
CSCOverlapsAlignmentAlgorithm::m_fiducial_MEx1
TProfile * m_fiducial_MEx1
Definition: CSCOverlapsAlignmentAlgorithm.h:80
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:49
TransientTrackBuilder
Definition: TransientTrackBuilder.h:16
CSCOverlapsAlignmentAlgorithm::m_residualsConstraints
std::vector< CSCPairResidualsConstraint * > m_residualsConstraints
Definition: CSCOverlapsAlignmentAlgorithm.h:152
TFileService::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileService.h:64
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
edmIntegrityCheck.report
report
Definition: edmIntegrityCheck.py:348
Trajectory::measurements
DataContainer const & measurements() const
Definition: Trajectory.h:178
CSCOverlapsAlignmentAlgorithm::m_fiducial
bool m_fiducial
Definition: CSCOverlapsAlignmentAlgorithm.h:123
CSCOverlapsAlignmentAlgorithm::m_XYpos_mem1
TH2F * m_XYpos_mem1
Definition: CSCOverlapsAlignmentAlgorithm.h:107
CSCOverlapsAlignmentAlgorithm::m_doAlignment
bool m_doAlignment
Definition: CSCOverlapsAlignmentAlgorithm.h:146
dumpTauVariables_cfi.eventInfo
eventInfo
add run, event number and lumi section
Definition: dumpTauVariables_cfi.py:12
SiPixelPI::two
Definition: SiPixelPayloadInspectorHelper.h:39
reco::TransientTrack
Definition: TransientTrack.h:19
CSCChamberFitter
Definition: CSCChamberFitter.h:24
CSCOverlapsAlignmentAlgorithm::m_RPhipos_mep3
TH2F * m_RPhipos_mep3
Definition: CSCOverlapsAlignmentAlgorithm.h:113
CSCOverlapsAlignmentAlgorithm::m_slope_MEp4
TH1F * m_slope_MEp4
Definition: CSCOverlapsAlignmentAlgorithm.h:84
Trajectory
Definition: Trajectory.h:38
AlignmentAlgorithmBase::AlignmentAlgorithmBase
AlignmentAlgorithmBase(const edm::ParameterSet &, const edm::ConsumesCollector &)
Constructor.
Definition: AlignmentAlgorithmBase.h:115
edm::Transition::BeginRun
Exception
Definition: hltDiff.cc:245
CSCOverlapsAlignmentAlgorithm::m_slopeResiduals_normalized
TH1F * m_slopeResiduals_normalized
Definition: CSCOverlapsAlignmentAlgorithm.h:95
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
CSCOverlapsAlignmentAlgorithm::m_XYpos_mem3
TH2F * m_XYpos_mem3
Definition: CSCOverlapsAlignmentAlgorithm.h:109
CSCOverlapsAlignmentAlgorithm::m_hitsPerChamber
TH1F * m_hitsPerChamber
Definition: CSCOverlapsAlignmentAlgorithm.h:76
CSCOverlapsAlignmentAlgorithm::m_XYpos_mem2
TH2F * m_XYpos_mem2
Definition: CSCOverlapsAlignmentAlgorithm.h:108
TransientTrackBuilder::build
reco::TransientTrack build(const reco::Track *p) const
Definition: TransientTrackBuilder.cc:20
CSCOverlapsAlignmentAlgorithm::m_RPhipos_mem1
TH2F * m_RPhipos_mem1
Definition: CSCOverlapsAlignmentAlgorithm.h:115
CSCOverlapsAlignmentAlgorithm::m_slope
TH1F * m_slope
Definition: CSCOverlapsAlignmentAlgorithm.h:83
cms::Exception
Definition: Exception.h:70
DetId::Muon
Definition: DetId.h:26
CSCOverlapsAlignmentAlgorithm::m_slope_MEp2
TH1F * m_slope_MEp2
Definition: CSCOverlapsAlignmentAlgorithm.h:86
genParticles_cff.map
map
Definition: genParticles_cff.py:11
CSCOverlapsAlignmentAlgorithm::m_maxRedChi2
double m_maxRedChi2
Definition: CSCOverlapsAlignmentAlgorithm.h:143
CSCOverlapsAlignmentAlgorithm::m_slope_MEm2
TH1F * m_slope_MEm2
Definition: CSCOverlapsAlignmentAlgorithm.h:89
CSCOverlapsAlignmentAlgorithm::m_XYpos_mem4
TH2F * m_XYpos_mem4
Definition: CSCOverlapsAlignmentAlgorithm.h:110
CSCOverlapsAlignmentAlgorithm::m_alignableNavigator
AlignableNavigator * m_alignableNavigator
Definition: CSCOverlapsAlignmentAlgorithm.h:150
CSCOverlapsAlignmentAlgorithm::m_RPhipos_mep2
TH2F * m_RPhipos_mep2
Definition: CSCOverlapsAlignmentAlgorithm.h:112
CSCOverlapsAlignmentAlgorithm::m_RPhipos_mem2
TH2F * m_RPhipos_mem2
Definition: CSCOverlapsAlignmentAlgorithm.h:116
CSCOverlapsAlignmentAlgorithm::m_slopeFromTrackRefit
bool m_slopeFromTrackRefit
Definition: CSCOverlapsAlignmentAlgorithm.h:125
CSCOverlapsAlignmentAlgorithm::m_minTracksPerOverlap
int m_minTracksPerOverlap
Definition: CSCOverlapsAlignmentAlgorithm.h:132
CSCOverlapsAlignmentAlgorithm::m_combineME11
bool m_combineME11
Definition: CSCOverlapsAlignmentAlgorithm.h:129