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

Inheritance diagram for MuonAlignmentFromReference:

Public Member Functions
void	initialize (const edm::EventSetup &iSetup, AlignableTracker alignableTracker, AlignableMuon alignableMuon, AlignableExtras extras, AlignmentParameterStore alignmentParameterStore) override
	Call at beginning of job (must be implemented in derived class) More...

	MuonAlignmentFromReference (const edm::ParameterSet &cfg, edm::ConsumesCollector &iC)

void	processMuonResidualsFromTrack (MuonResidualsFromTrack &mrft)

void	run (const edm::EventSetup &iSetup, const EventInfo &eventInfo) override
	Run the algorithm (must be implemented in derived class) More...

void	startNewLoop () override

void	terminate (const edm::EventSetup &iSetup) override
	Call at end of each loop (must be implemented in derived class) More...

	~MuonAlignmentFromReference () 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 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...

Private Member Functions
void	bookNtuple ()

std::string	chamberPrettyNameFromId (unsigned int idx)

void	correctBField ()

void	eraseNotSelectedResiduals ()

void	fiducialCuts ()

void	fillNtuple ()

void	fitAndAlign ()

int	number (std::string s)

bool	numeric (std::string s)

void	parseReference (align::Alignables &reference, const align::Alignables &all_DT_chambers, const align::Alignables &all_CSC_chambers)

void	readTmpFiles ()

void	selectResidualsPeaks ()

void	writeTmpFiles ()

Private Attributes
AlignableNavigator *	m_alignableNavigator

align::Alignables	m_alignables

AlignmentParameterStore *	m_alignmentParameterStore

bool	m_allowTIDTEC

int	m_BFieldCorrection

const MuonResidualsFromTrack::BuilderToken	m_builderToken

bool	m_combineME11

long	m_counter_csc

long	m_counter_cscaligning

long	m_counter_cschits

long	m_counter_cscvalid

long	m_counter_events

long	m_counter_minchambers

long	m_counter_resslopey

long	m_counter_station123

long	m_counter_station123aligning

long	m_counter_station123dt13hits

long	m_counter_station123dt2hits

long	m_counter_station123valid

long	m_counter_station4

long	m_counter_station4aligning

long	m_counter_station4hits

long	m_counter_station4valid

long	m_counter_totchambers

long	m_counter_trackdxy

long	m_counter_trackerchi2

long	m_counter_trackerhits

long	m_counter_trackertidtec

long	m_counter_trackmomentum

long	m_counter_tracks

bool	m_createNtuple

const edm::ESGetToken< CSCGeometry, MuonGeometryRecord >	m_cscGeometryToken

bool	m_debug

const edm::ESGetToken< DetIdAssociator, DetIdAssociatorRecord >	m_DetIdToken

bool	m_doAlignment

bool	m_doCSC

bool	m_doDT

std::map< unsigned int, MuonResidualsTwoBin * >	m_fitterOrder

std::map< Alignable , MuonResidualsTwoBin >	m_fitters

const edm::ESGetToken< GlobalTrackingGeometry, GlobalTrackingGeometryRecord >	m_globTackingToken

std::vector< unsigned int >	m_indexes

const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	m_MagFieldToken

double	m_maxDxy

double	m_maxResSlopeY

double	m_maxTrackerRedChi2

double	m_maxTrackP

double	m_maxTrackPt

std::map< Alignable , Alignable >	m_me11map

int	m_minAlignmentHits

int	m_minCSCHits

int	m_minDT13Hits

int	m_minDT2Hits

int	m_minNCrossedChambers

int	m_minTrackerHits

double	m_minTrackP

double	m_minTrackPt

edm::InputTag	m_muonCollectionTag

double	m_peakNSigma

const edm::ESGetToken< Propagator, TrackingComponentsRecord >	m_propToken

std::vector< std::string >	m_readTemporaryFiles

std::vector< std::string >	m_reference

std::string	m_reportFileName

std::string	m_residualsModel

int	m_strategy

MuonResidualsFitter::MuonAlignmentTreeRow	m_tree_row

TTree *	m_ttree

bool	m_twoBin

std::string	m_useResiduals

bool	m_weightAlignment

std::string	m_writeTemporaryFile

Additional Inherited Members
Public Types inherited from AlignmentAlgorithmBase
typedef std::pair< const Trajectory , const reco::Track >	ConstTrajTrackPair

typedef std::vector< ConstTrajTrackPair >	ConstTrajTrackPairCollection

using	RunNumber = align::RunNumber

using	RunRange = align::RunRange

Detailed Description

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

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

Definition at line 74 of file MuonAlignmentFromReference.cc.

Constructor & Destructor Documentation

◆ MuonAlignmentFromReference()

MuonAlignmentFromReference::MuonAlignmentFromReference	(	const edm::ParameterSet &	cfg,
		edm::ConsumesCollector &	iC
	)

Definition at line 196 of file MuonAlignmentFromReference.cc.

References bookNtuple(), compareTotals::fs, m_counter_csc, m_counter_cscaligning, m_counter_cschits, m_counter_cscvalid, m_counter_events, m_counter_minchambers, m_counter_resslopey, m_counter_station123, m_counter_station123aligning, m_counter_station123dt13hits, m_counter_station123dt2hits, m_counter_station123valid, m_counter_station4, m_counter_station4aligning, m_counter_station4hits, m_counter_station4valid, m_counter_totchambers, m_counter_trackdxy, m_counter_trackerchi2, m_counter_trackerhits, m_counter_trackertidtec, m_counter_trackmomentum, m_counter_tracks, m_createNtuple, m_debug, m_doAlignment, and m_ttree.

     : AlignmentAlgorithmBase(cfg, iC),
       m_cscGeometryToken(iC.esConsumes<edm::Transition::BeginRun>()),
       m_globTackingToken(iC.esConsumes()),
       m_MagFieldToken(iC.esConsumes()),
       m_propToken(iC.esConsumes(edm::ESInputTag("", "SteppingHelixPropagatorAny"))),
       m_DetIdToken(iC.esConsumes(edm::ESInputTag("", "MuonDetIdAssociator"))),
       m_builderToken(iC.esConsumes(MuonResidualsFromTrack::builderESInputTag())),
       m_muonCollectionTag(cfg.getParameter<edm::InputTag>("muonCollectionTag")),
       m_reference(cfg.getParameter<std::vector<std::string> >("reference")),
       m_minTrackPt(cfg.getParameter<double>("minTrackPt")),
       m_maxTrackPt(cfg.getParameter<double>("maxTrackPt")),
       m_minTrackP(cfg.getParameter<double>("minTrackP")),
       m_maxTrackP(cfg.getParameter<double>("maxTrackP")),
       m_maxDxy(cfg.getParameter<double>("maxDxy")),
       m_minTrackerHits(cfg.getParameter<int>("minTrackerHits")),
       m_maxTrackerRedChi2(cfg.getParameter<double>("maxTrackerRedChi2")),
       m_allowTIDTEC(cfg.getParameter<bool>("allowTIDTEC")),
       m_minNCrossedChambers(cfg.getParameter<int>("minNCrossedChambers")),
       m_minDT13Hits(cfg.getParameter<int>("minDT13Hits")),
       m_minDT2Hits(cfg.getParameter<int>("minDT2Hits")),
       m_minCSCHits(cfg.getParameter<int>("minCSCHits")),
       m_writeTemporaryFile(cfg.getParameter<std::string>("writeTemporaryFile")),
       m_readTemporaryFiles(cfg.getParameter<std::vector<std::string> >("readTemporaryFiles")),
       m_doAlignment(cfg.getParameter<bool>("doAlignment")),
       m_strategy(cfg.getParameter<int>("strategy")),
       m_residualsModel(cfg.getParameter<std::string>("residualsModel")),
       m_minAlignmentHits(cfg.getParameter<int>("minAlignmentHits")),
       m_twoBin(cfg.getParameter<bool>("twoBin")),
       m_combineME11(cfg.getParameter<bool>("combineME11")),
       m_weightAlignment(cfg.getParameter<bool>("weightAlignment")),
       m_reportFileName(cfg.getParameter<std::string>("reportFileName")),
       m_maxResSlopeY(cfg.getParameter<double>("maxResSlopeY")),
       m_createNtuple(cfg.getParameter<bool>("createNtuple")),
       m_peakNSigma(cfg.getParameter<double>("peakNSigma")),
       m_BFieldCorrection(cfg.getParameter<int>("bFieldCorrection")),
       m_doDT(cfg.getParameter<bool>("doDT")),
       m_doCSC(cfg.getParameter<bool>("doCSC")),
       m_useResiduals(cfg.getParameter<std::string>("useResiduals")) {
   // alignment requires a TFile to provide plots to check the fit output
   // just filling the residuals lists does not
   // but we don't want to wait until the end of the job to find out that the TFile is missing
   if (m_doAlignment || m_createNtuple) {
     edm::Service<TFileService> fs;
     TFile& tfile = fs->file();
     tfile.ls();
   }
 
   m_ttree = nullptr;
   if (m_createNtuple)
     bookNtuple();
 
   m_counter_events = 0;
   m_counter_tracks = 0;
   m_counter_trackmomentum = 0;
   m_counter_trackdxy = 0;
   m_counter_trackerhits = 0;
   m_counter_trackerchi2 = 0;
   m_counter_trackertidtec = 0;
   m_counter_minchambers = 0;
   m_counter_totchambers = 0;
   m_counter_station123 = 0;
   m_counter_station123valid = 0;
   m_counter_station123dt13hits = 0;
   m_counter_station123dt2hits = 0;
   m_counter_station123aligning = 0;
   m_counter_station4 = 0;
   m_counter_station4valid = 0;
   m_counter_station4hits = 0;
   m_counter_station4aligning = 0;
   m_counter_csc = 0;
   m_counter_cscvalid = 0;
   m_counter_cschits = 0;
   m_counter_cscaligning = 0;
   m_counter_resslopey = 0;
 
   m_debug = false;
 }

◆ ~MuonAlignmentFromReference()

MuonAlignmentFromReference::~MuonAlignmentFromReference ( )

override

Definition at line 275 of file MuonAlignmentFromReference.cc.

References m_alignableNavigator.

275 { delete m_alignableNavigator; }

MuonAlignmentFromReference::m_alignableNavigator

AlignableNavigator * m_alignableNavigator

Definition: MuonAlignmentFromReference.cc:155

Member Function Documentation

◆ bookNtuple()

void MuonAlignmentFromReference::bookNtuple ( )

private

Definition at line 277 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference().

                                             {
   edm::Service<TFileService> fs;
   m_ttree = fs->make<TTree>("mual_ttree", "mual_ttree");
   m_ttree->Branch("is_plus", &m_tree_row.is_plus, "is_plus/O");
   m_ttree->Branch("is_dt", &m_tree_row.is_dt, "is_dt/O");
   m_ttree->Branch("station", &m_tree_row.station, "station/b");
   m_ttree->Branch("ring_wheel", &m_tree_row.ring_wheel, "ring_wheel/B");
   m_ttree->Branch("sector", &m_tree_row.sector, "sector/b");
   m_ttree->Branch("res_x", &m_tree_row.res_x, "res_x/F");
   m_ttree->Branch("res_y", &m_tree_row.res_y, "res_y/F");
   m_ttree->Branch("res_slope_x", &m_tree_row.res_slope_x, "res_slope_x/F");
   m_ttree->Branch("res_slope_y", &m_tree_row.res_slope_y, "res_slope_y/F");
   m_ttree->Branch("pos_x", &m_tree_row.pos_x, "pos_x/F");
   m_ttree->Branch("pos_y", &m_tree_row.pos_y, "pos_y/F");
   m_ttree->Branch("angle_x", &m_tree_row.angle_x, "angle_x/F");
   m_ttree->Branch("angle_y", &m_tree_row.angle_y, "angle_y/F");
   m_ttree->Branch("pz", &m_tree_row.pz, "pz/F");
   m_ttree->Branch("pt", &m_tree_row.pt, "pt/F");
   m_ttree->Branch("q", &m_tree_row.q, "q/B");
   m_ttree->Branch("select", &m_tree_row.select, "select/O");
   //m_ttree->Branch("",&m_tree_row.,"/");
 }

◆ chamberPrettyNameFromId()

std::string MuonAlignmentFromReference::chamberPrettyNameFromId ( unsigned int idx )

private

Definition at line 1756 of file MuonAlignmentFromReference.cc.

References CSCDetId::chamber(), MuonSubdetId::CSC, MuonSubdetId::DT, CSCDetId::endcap(), l1ctLayer2EG_cff::id, heavyIonCSV_trainingSettings::idx, CSCDetId::ring(), DTChamberId::sector(), DTChamberId::station(), CSCDetId::station(), AlCaHLTBitMon_QueryRunRegistry::string, and DTChamberId::wheel().

Referenced by correctBField(), eraseNotSelectedResiduals(), fiducialCuts(), and selectResidualsPeaks().

                                                                               {
   DetId id(idx);
   char cname[40];
   if (id.subdetId() == MuonSubdetId::DT) {
     DTChamberId chamberId(id.rawId());
     sprintf(cname, "MB%+d/%d/%02d", chamberId.wheel(), chamberId.station(), chamberId.sector());
   } else if (id.subdetId() == MuonSubdetId::CSC) {
     CSCDetId chamberId(id.rawId());
     sprintf(cname,
             "ME%s%d/%d/%02d",
             (chamberId.endcap() == 1 ? "+" : "-"),
             chamberId.station(),
             chamberId.ring(),
             chamberId.chamber());
   }
   return std::string(cname);
 }

◆ correctBField()

void MuonAlignmentFromReference::correctBField ( )

private

Definition at line 1658 of file MuonAlignmentFromReference.cc.

References chamberPrettyNameFromId(), MuonResidualsTwoBin::correctBField(), gather_cfg::cout, m_debug, m_fitterOrder, and m_indexes.

Referenced by terminate().

                                                {
   bool m_debug = false;
 
   for (std::vector<unsigned int>::const_iterator index = m_indexes.begin(); index != m_indexes.end(); ++index) {
     if (m_debug)
       std::cout << "correcting B in " << chamberPrettyNameFromId(*index) << std::endl;
     MuonResidualsTwoBin* fitter = m_fitterOrder[*index];
     fitter->correctBField();
   }
 }

◆ eraseNotSelectedResiduals()

void MuonAlignmentFromReference::eraseNotSelectedResiduals ( )

private

Definition at line 1678 of file MuonAlignmentFromReference.cc.

References chamberPrettyNameFromId(), gather_cfg::cout, MuonResidualsTwoBin::eraseNotSelectedResiduals(), m_debug, m_fitterOrder, and m_indexes.

Referenced by terminate().

                                                            {
   for (std::vector<unsigned int>::const_iterator index = m_indexes.begin(); index != m_indexes.end(); ++index) {
     if (m_debug)
       std::cout << "erasing in " << chamberPrettyNameFromId(*index) << std::endl;
     MuonResidualsTwoBin* fitter = m_fitterOrder[*index];
     fitter->eraseNotSelectedResiduals();
   }
 }

◆ fiducialCuts()

void MuonAlignmentFromReference::fiducialCuts ( )

private

Definition at line 1669 of file MuonAlignmentFromReference.cc.

References chamberPrettyNameFromId(), gather_cfg::cout, MuonResidualsTwoBin::fiducialCuts(), m_debug, m_fitterOrder, and m_indexes.

Referenced by terminate().

                                               {
   for (std::vector<unsigned int>::const_iterator index = m_indexes.begin(); index != m_indexes.end(); ++index) {
     if (m_debug)
       std::cout << "applying fiducial cuts in " << chamberPrettyNameFromId(*index) << std::endl;
     MuonResidualsTwoBin* fitter = m_fitterOrder[*index];
     fitter->fiducialCuts();
   }
 }

◆ fillNtuple()

void MuonAlignmentFromReference::fillNtuple ( )

private

Definition at line 1774 of file MuonAlignmentFromReference.cc.

Referenced by terminate().

                                             {
   // WARNING: does not support two bin option!!!
 
   for (std::vector<unsigned int>::const_iterator index = m_indexes.begin(); index != m_indexes.end(); ++index) {
     DetId detid(*index);
     if (detid.det() != DetId::Muon || !(detid.subdetId() == MuonSubdetId::DT || detid.subdetId() == MuonSubdetId::CSC))
       assert(false);
 
     if (detid.subdetId() == MuonSubdetId::DT) {
       m_tree_row.is_dt = (Bool_t) true;
       DTChamberId id(*index);
       m_tree_row.is_plus = (Bool_t) true;
       m_tree_row.station = (UChar_t)id.station();
       m_tree_row.ring_wheel = (Char_t)id.wheel();
       m_tree_row.sector = (UChar_t)id.sector();
     } else {
       m_tree_row.is_dt = (Bool_t) false;
       CSCDetId id(*index);
       m_tree_row.is_plus = (Bool_t)(id.endcap() == 1);
       m_tree_row.station = (UChar_t)id.station();
       m_tree_row.ring_wheel = (Char_t)id.ring();
       m_tree_row.sector = (UChar_t)id.chamber();
     }
 
     MuonResidualsTwoBin* fitter = m_fitterOrder[*index];
 
     std::vector<double*>::const_iterator residual = fitter->residualsPos_begin();
     std::vector<bool>::const_iterator residual_ok = fitter->residualsPos_ok_begin();
     for (; residual != fitter->residualsPos_end(); ++residual, ++residual_ok) {
       if (fitter->type() == MuonResidualsFitter::k5DOF || fitter->type() == MuonResidualsFitter::k6DOFrphi) {
         m_tree_row.res_x = (Float_t)(*residual)[MuonResiduals5DOFFitter::kResid];
         m_tree_row.res_y = (Float_t)0.;
         m_tree_row.res_slope_x = (Float_t)(*residual)[MuonResiduals5DOFFitter::kResSlope];
         m_tree_row.res_slope_y = (Float_t)0.;
         m_tree_row.pos_x = (Float_t)(*residual)[MuonResiduals5DOFFitter::kPositionX];
         m_tree_row.pos_y = (Float_t)(*residual)[MuonResiduals5DOFFitter::kPositionY];
         m_tree_row.angle_x = (Float_t)(*residual)[MuonResiduals5DOFFitter::kAngleX];
         m_tree_row.angle_y = (Float_t)(*residual)[MuonResiduals5DOFFitter::kAngleY];
         m_tree_row.pz = (Float_t)(*residual)[MuonResiduals5DOFFitter::kPz];
         m_tree_row.pt = (Float_t)(*residual)[MuonResiduals5DOFFitter::kPt];
         m_tree_row.q = (Char_t)(*residual)[MuonResiduals5DOFFitter::kCharge];
         m_tree_row.select = (Bool_t)*residual_ok;
       } else if (fitter->type() == MuonResidualsFitter::k6DOF) {
         m_tree_row.res_x = (Float_t)(*residual)[MuonResiduals6DOFFitter::kResidX];
         m_tree_row.res_y = (Float_t)(*residual)[MuonResiduals6DOFFitter::kResidY];
         m_tree_row.res_slope_x = (Float_t)(*residual)[MuonResiduals6DOFFitter::kResSlopeX];
         m_tree_row.res_slope_y = (Float_t)(*residual)[MuonResiduals6DOFFitter::kResSlopeY];
         m_tree_row.pos_x = (Float_t)(*residual)[MuonResiduals6DOFFitter::kPositionX];
         m_tree_row.pos_y = (Float_t)(*residual)[MuonResiduals6DOFFitter::kPositionY];
         m_tree_row.angle_x = (Float_t)(*residual)[MuonResiduals6DOFFitter::kAngleX];
         m_tree_row.angle_y = (Float_t)(*residual)[MuonResiduals6DOFFitter::kAngleY];
         m_tree_row.pz = (Float_t)(*residual)[MuonResiduals6DOFFitter::kPz];
         m_tree_row.pt = (Float_t)(*residual)[MuonResiduals6DOFFitter::kPt];
         m_tree_row.q = (Char_t)(*residual)[MuonResiduals6DOFFitter::kCharge];
         m_tree_row.select = (Bool_t)*residual_ok;
       } else
         assert(false);
 
       m_ttree->Fill();
     }
   }
 }

◆ fitAndAlign()

void MuonAlignmentFromReference::fitAndAlign ( )

private

Definition at line 835 of file MuonAlignmentFromReference.cc.

Referenced by terminate().

                                              {
   bool m_debug = false;
 
   edm::Service<TFileService> tfileService;
   TFileDirectory rootDirectory(tfileService->mkdir("MuonAlignmentFromReference"));
 
   std::ofstream report;
   bool writeReport = (m_reportFileName != std::string(""));
   if (writeReport) {
     report.open(m_reportFileName.c_str());
     report << "nan = None;  NAN = None" << std::endl;
     report << "nan = 0" << std::endl;
     report << "reports = []" << std::endl;
     report << "class ValErr:" << std::endl
            << "    def __init__(self, value, error, antisym):" << std::endl
            << "        self.value, self.error, self.antisym = value, error, antisym" << std::endl
            << "" << std::endl
            << "    def __repr__(self):" << std::endl
            << "        if self.antisym == 0.:" << std::endl
            << "            return \"%g +- %g\" % (self.value, self.error)" << std::endl
            << "        else:" << std::endl
            << "            return \"%g +- %g ~ %g\" % (self.value, self.error, self.antisym)" << std::endl
            << "" << std::endl
            << "class Report:" << std::endl
            << "    def __init__(self, chamberId, postal_address, name):" << std::endl
            << "        self.chamberId, self.postal_address, self.name = chamberId, postal_address, name" << std::endl
            << "        self.status = \"NOFIT\"" << std::endl
            << "        self.fittype = None" << std::endl
            << "" << std::endl
            << "    def add_parameters(self, deltax, deltay, deltaz, deltaphix, deltaphiy, deltaphiz, loglikelihood, "
               "numsegments, sumofweights, redchi2):"
            << std::endl
            << "        self.status = \"PASS\"" << std::endl
            << "        self.deltax, self.deltay, self.deltaz, self.deltaphix, self.deltaphiy, self.deltaphiz = deltax, "
               "deltay, deltaz, deltaphix, deltaphiy, deltaphiz"
            << std::endl
            << "        self.loglikelihood, self.numsegments, self.sumofweights, self.redchi2 = loglikelihood, "
               "numsegments, sumofweights, redchi2"
            << std::endl
            << "" << std::endl
            << "    def add_stats(self, median_x, median_y, median_dxdz, median_dydz, mean30_x, mean30_y, mean20_dxdz, "
               "mean50_dydz, mean15_x, mean15_y, mean10_dxdz, mean25_dydz, wmean30_x, wmean30_y, wmean20_dxdz, "
               "wmean50_dydz, wmean15_x, wmean15_y, wmean10_dxdz, wmean25_dydz, stdev30_x, stdev30_y, stdev20_dxdz, "
               "stdev50_dydz, stdev15_x, stdev15_y, stdev10_dxdz, stdev25_dydz):"
            << std::endl
            << "        self.median_x, self.median_y, self.median_dxdz, self.median_dydz, self.mean30_x, self.mean30_y, "
               "self.mean20_dxdz, self.mean50_dydz, self.mean15_x, self.mean15_y, self.mean10_dxdz, self.mean25_dydz, "
               "self.wmean30_x, self.wmean30_y, self.wmean20_dxdz, self.wmean50_dydz, self.wmean15_x, self.wmean15_y, "
               "self.wmean10_dxdz, self.wmean25_dydz, self.stdev30_x, self.stdev30_y, self.stdev20_dxdz, "
               "self.stdev50_dydz, self.stdev15_x, self.stdev15_y, self.stdev10_dxdz, self.stdev25_dydz = median_x, "
               "median_y, median_dxdz, median_dydz, mean30_x, mean30_y, mean20_dxdz, mean50_dydz, mean15_x, mean15_y, "
               "mean10_dxdz, mean25_dydz, wmean30_x, wmean30_y, wmean20_dxdz, wmean50_dydz, wmean15_x, wmean15_y, "
               "wmean10_dxdz, wmean25_dydz, stdev30_x, stdev30_y, stdev20_dxdz, stdev50_dydz, stdev15_x, stdev15_y, "
               "stdev10_dxdz, stdev25_dydz"
            << std::endl
            << "" << std::endl
            << "    def __repr__(self):" << std::endl
            << "        return \"<Report %s %s %s>\" % (self.postal_address[0], \" \".join(map(str, "
               "self.postal_address[1:])), self.status)"
            << std::endl
            << std::endl;
   }
 
   if (m_debug)
     std::cout << "***** just after report.open" << std::endl;
 
   for (const auto& ali : m_alignables) {
     if (m_debug)
       std::cout << "***** Start loop over alignables" << std::endl;
 
     std::vector<bool> selector = ali->alignmentParameters()->selector();
     bool align_x = selector[0];
     bool align_y = selector[1];
     bool align_z = selector[2];
     bool align_phix = selector[3];
     bool align_phiy = selector[4];
     bool align_phiz = selector[5];
     int numParams = ((align_x ? 1 : 0) + (align_y ? 1 : 0) + (align_z ? 1 : 0) + (align_phix ? 1 : 0) +
                      (align_phiy ? 1 : 0) + (align_phiz ? 1 : 0));
 
     // map from 0-5 to the index of params, above
     std::vector<int> paramIndex;
     int paramIndex_counter = -1;
     if (align_x)
       paramIndex_counter++;
     paramIndex.push_back(paramIndex_counter);
     if (align_y)
       paramIndex_counter++;
     paramIndex.push_back(paramIndex_counter);
     if (align_z)
       paramIndex_counter++;
     paramIndex.push_back(paramIndex_counter);
     if (align_phix)
       paramIndex_counter++;
     paramIndex.push_back(paramIndex_counter);
     if (align_phiy)
       paramIndex_counter++;
     paramIndex.push_back(paramIndex_counter);
     if (align_phiz)
       paramIndex_counter++;
     paramIndex.push_back(paramIndex_counter);
 
     DetId id = ali->geomDetId();
 
     auto thisali = ali;
     if (m_combineME11 && id.subdetId() == MuonSubdetId::CSC) {
       CSCDetId cscid(id.rawId());
       if (cscid.station() == 1 && cscid.ring() == 4)
         thisali = m_me11map[ali];
     }
 
     if (m_debug)
       std::cout << "***** loop over alignables 1" << std::endl;
 
     char cname[40];
     char wheel_label[][2] = {"A", "B", "C", "D", "E"};
 
     if (id.subdetId() == MuonSubdetId::DT) {
       DTChamberId chamberId(id.rawId());
 
       //if ( ! ( (chamberId.station()==1&&chamberId.wheel()==0) || (chamberId.station()==4&&chamberId.wheel()==2) ) ) continue;
 
       sprintf(cname, "MBwh%sst%dsec%02d", wheel_label[chamberId.wheel() + 2], chamberId.station(), chamberId.sector());
       if (writeReport) {
         report << "reports.append(Report(" << id.rawId() << ", (\"DT\", " << chamberId.wheel() << ", "
                << chamberId.station() << ", " << chamberId.sector() << "), \"" << cname << "\"))" << std::endl;
       }
     } else if (id.subdetId() == MuonSubdetId::CSC) {
       CSCDetId chamberId(id.rawId());
       sprintf(cname,
               "ME%s%d%d_%02d",
               (chamberId.endcap() == 1 ? "p" : "m"),
               chamberId.station(),
               chamberId.ring(),
               chamberId.chamber());
 
       //if ( chamberId.chamber()>6 || chamberId.endcap()==2 || ! ( (chamberId.station()==2&&chamberId.ring()==1) || (chamberId.station()==3&&chamberId.ring()==2) ) ) continue;
 
       if (writeReport) {
         report << "reports.append(Report(" << id.rawId() << ", (\"CSC\", " << chamberId.endcap() << ", "
                << chamberId.station() << ", " << chamberId.ring() << ", " << chamberId.chamber() << "), \"" << cname
                << "\"))" << std::endl;
       }
     }
 
     if (m_debug)
       std::cout << "***** loop over alignables 2" << std::endl;
 
     //if(! ( strcmp(cname,"MBwhCst3sec12")==0 || strcmp(cname,"MBwhCst3sec06")==0)) continue;
 
     std::map<Alignable*, MuonResidualsTwoBin*>::const_iterator fitter = m_fitters.find(thisali);
 
     if (m_debug)
       std::cout << "***** loop over alignables 3" << std::endl;
 
     if (fitter != m_fitters.end()) {
       //if (fitter->second->type() != MuonResidualsFitter::k6DOFrphi) continue;
 
       TStopwatch stop_watch;
       stop_watch.Start();
 
       // MINUIT is verbose in std::cout anyway
       if (m_debug)
         std::cout << "============================================================================================="
                   << std::endl;
       if (m_debug)
         std::cout << "Fitting " << cname << std::endl;
 
       if (writeReport) {
         report << "reports[-1].posNum = " << fitter->second->numResidualsPos() << std::endl;
         report << "reports[-1].negNum = " << fitter->second->numResidualsNeg() << std::endl;
       }
 
       if (fitter->second->type() == MuonResidualsFitter::k5DOF) {
         if (!align_x)
           fitter->second->fix(MuonResiduals5DOFFitter::kAlignX);
         if (!align_z)
           fitter->second->fix(MuonResiduals5DOFFitter::kAlignZ);
         if (!align_phix)
           fitter->second->fix(MuonResiduals5DOFFitter::kAlignPhiX);
         if (!align_phiy)
           fitter->second->fix(MuonResiduals5DOFFitter::kAlignPhiY);
         if (!align_phiz)
           fitter->second->fix(MuonResiduals5DOFFitter::kAlignPhiZ);
       } else if (fitter->second->type() == MuonResidualsFitter::k6DOF) {
         if (!align_x)
           fitter->second->fix(MuonResiduals6DOFFitter::kAlignX);
         if (!align_y)
           fitter->second->fix(MuonResiduals6DOFFitter::kAlignY);
         if (!align_z)
           fitter->second->fix(MuonResiduals6DOFFitter::kAlignZ);
         if (!align_phix)
           fitter->second->fix(MuonResiduals6DOFFitter::kAlignPhiX);
         if (!align_phiy)
           fitter->second->fix(MuonResiduals6DOFFitter::kAlignPhiY);
         if (!align_phiz)
           fitter->second->fix(MuonResiduals6DOFFitter::kAlignPhiZ);
       } else if (fitter->second->type() == MuonResidualsFitter::k6DOFrphi) {
         if (!align_x)
           fitter->second->fix(MuonResiduals6DOFrphiFitter::kAlignX);
         if (!align_y)
           fitter->second->fix(MuonResiduals6DOFrphiFitter::kAlignY);
         if (!align_z)
           fitter->second->fix(MuonResiduals6DOFrphiFitter::kAlignZ);
         if (!align_phix)
           fitter->second->fix(MuonResiduals6DOFrphiFitter::kAlignPhiX);
         if (!align_phiy)
           fitter->second->fix(MuonResiduals6DOFrphiFitter::kAlignPhiY);
         if (!align_phiz)
           fitter->second->fix(MuonResiduals6DOFrphiFitter::kAlignPhiZ);
       } else
         assert(false);
 
       if (m_debug)
         std::cout << "***** loop over alignables 4" << std::endl;
 
       AlgebraicVector params(numParams);
       AlgebraicSymMatrix cov(numParams);
 
       if (fitter->second->numsegments() >= m_minAlignmentHits) {
         if (m_debug)
           std::cout << "***** loop over alignables 5" << std::endl;
 
         bool successful_fit = fitter->second->fit(thisali);
 
         if (m_debug)
           std::cout << "***** loop over alignables 6 " << fitter->second->type() << std::endl;
 
         double loglikelihood = fitter->second->loglikelihood();
         double numsegments = fitter->second->numsegments();
         double sumofweights = fitter->second->sumofweights();
         double redchi2 = fitter->second->plot(cname, &rootDirectory, thisali);
 
         if (fitter->second->type() == MuonResidualsFitter::k5DOF) {
           if (m_debug)
             std::cout << "***** loop over alignables k5DOF" << std::endl;
 
           double deltax_value = fitter->second->value(MuonResiduals5DOFFitter::kAlignX);
           double deltax_error = fitter->second->errorerror(MuonResiduals5DOFFitter::kAlignX);
           double deltax_antisym = fitter->second->antisym(MuonResiduals5DOFFitter::kAlignX);
 
           double deltaz_value = fitter->second->value(MuonResiduals5DOFFitter::kAlignZ);
           double deltaz_error = fitter->second->errorerror(MuonResiduals5DOFFitter::kAlignZ);
           double deltaz_antisym = fitter->second->antisym(MuonResiduals5DOFFitter::kAlignZ);
 
           double deltaphix_value = fitter->second->value(MuonResiduals5DOFFitter::kAlignPhiX);
           double deltaphix_error = fitter->second->errorerror(MuonResiduals5DOFFitter::kAlignPhiX);
           double deltaphix_antisym = fitter->second->antisym(MuonResiduals5DOFFitter::kAlignPhiX);
 
           double deltaphiy_value = fitter->second->value(MuonResiduals5DOFFitter::kAlignPhiY);
           double deltaphiy_error = fitter->second->errorerror(MuonResiduals5DOFFitter::kAlignPhiY);
           double deltaphiy_antisym = fitter->second->antisym(MuonResiduals5DOFFitter::kAlignPhiY);
 
           double deltaphiz_value = fitter->second->value(MuonResiduals5DOFFitter::kAlignPhiZ);
           double deltaphiz_error = fitter->second->errorerror(MuonResiduals5DOFFitter::kAlignPhiZ);
           double deltaphiz_antisym = fitter->second->antisym(MuonResiduals5DOFFitter::kAlignPhiZ);
 
           double sigmaresid_value = fitter->second->value(MuonResiduals5DOFFitter::kResidSigma);
           double sigmaresid_error = fitter->second->errorerror(MuonResiduals5DOFFitter::kResidSigma);
           double sigmaresid_antisym = fitter->second->antisym(MuonResiduals5DOFFitter::kResidSigma);
 
           double sigmaresslope_value = fitter->second->value(MuonResiduals5DOFFitter::kResSlopeSigma);
           double sigmaresslope_error = fitter->second->errorerror(MuonResiduals5DOFFitter::kResSlopeSigma);
           double sigmaresslope_antisym = fitter->second->antisym(MuonResiduals5DOFFitter::kResSlopeSigma);
 
           double gammaresid_value, gammaresid_error, gammaresid_antisym, gammaresslope_value, gammaresslope_error,
               gammaresslope_antisym;
           gammaresid_value = gammaresid_error = gammaresid_antisym = gammaresslope_value = gammaresslope_error =
               gammaresslope_antisym = 0.;
 
           if (fitter->second->residualsModel() != MuonResidualsFitter::kPureGaussian &&
               fitter->second->residualsModel() != MuonResidualsFitter::kPureGaussian2D &&
               fitter->second->residualsModel() != MuonResidualsFitter::kGaussPowerTails) {
             gammaresid_value = fitter->second->value(MuonResiduals5DOFFitter::kResidGamma);
             gammaresid_error = fitter->second->errorerror(MuonResiduals5DOFFitter::kResidGamma);
             gammaresid_antisym = fitter->second->antisym(MuonResiduals5DOFFitter::kResidGamma);
 
             gammaresslope_value = fitter->second->value(MuonResiduals5DOFFitter::kResSlopeGamma);
             gammaresslope_error = fitter->second->errorerror(MuonResiduals5DOFFitter::kResSlopeGamma);
             gammaresslope_antisym = fitter->second->antisym(MuonResiduals5DOFFitter::kResSlopeGamma);
           }
 
           if (writeReport) {
             report << "reports[-1].fittype = \"5DOF\"" << std::endl;
             report << "reports[-1].add_parameters(ValErr(" << deltax_value << ", " << deltax_error << ", "
                    << deltax_antisym << "), \\" << std::endl
                    << "                           None, \\" << std::endl
                    << "                           ValErr(" << deltaz_value << ", " << deltaz_error << ", "
                    << deltaz_antisym << "), \\" << std::endl
                    << "                           ValErr(" << deltaphix_value << ", " << deltaphix_error << ", "
                    << deltaphix_antisym << "), \\" << std::endl
                    << "                           ValErr(" << deltaphiy_value << ", " << deltaphiy_error << ", "
                    << deltaphiy_antisym << "), \\" << std::endl
                    << "                           ValErr(" << deltaphiz_value << ", " << deltaphiz_error << ", "
                    << deltaphiz_antisym << "), \\" << std::endl
                    << "                           " << loglikelihood << ", " << numsegments << ", " << sumofweights
                    << ", " << redchi2 << ")" << std::endl;
             report << "reports[-1].sigmaresid = ValErr(" << sigmaresid_value << ", " << sigmaresid_error << ", "
                    << sigmaresid_antisym << ")" << std::endl;
             report << "reports[-1].sigmaresslope = ValErr(" << sigmaresslope_value << ", " << sigmaresslope_error
                    << ", " << sigmaresslope_antisym << ")" << std::endl;
             if (fitter->second->residualsModel() != MuonResidualsFitter::kPureGaussian &&
                 fitter->second->residualsModel() != MuonResidualsFitter::kPureGaussian2D &&
                 fitter->second->residualsModel() != MuonResidualsFitter::kGaussPowerTails) {
               report << "reports[-1].gammaresid = ValErr(" << gammaresid_value << ", " << gammaresid_error << ", "
                      << gammaresid_antisym << ")" << std::endl;
               report << "reports[-1].gammaresslope = ValErr(" << gammaresslope_value << ", " << gammaresslope_error
                      << ", " << gammaresslope_antisym << ")" << std::endl;
             }
 
             report << "reports[-1].add_stats(" << fitter->second->median(MuonResiduals5DOFFitter::kResid) << ", "
                    << "None, " << fitter->second->median(MuonResiduals5DOFFitter::kResSlope) << ", "
                    << "None, " << fitter->second->mean(MuonResiduals5DOFFitter::kResid, 30.) << ", "
                    << "None, " << fitter->second->mean(MuonResiduals5DOFFitter::kResSlope, 20.) << ", "
                    << "None, " << fitter->second->mean(MuonResiduals5DOFFitter::kResid, 15.) << ", "
                    << "None, " << fitter->second->mean(MuonResiduals5DOFFitter::kResSlope, 10.) << ", "
                    << "None, "
                    << fitter->second->wmean(MuonResiduals5DOFFitter::kResid, MuonResiduals5DOFFitter::kRedChi2, 30.)
                    << ", "
                    << "None, "
                    << fitter->second->wmean(MuonResiduals5DOFFitter::kResSlope, MuonResiduals5DOFFitter::kRedChi2, 20.)
                    << ", "
                    << "None, "
                    << fitter->second->wmean(MuonResiduals5DOFFitter::kResid, MuonResiduals5DOFFitter::kRedChi2, 15.)
                    << ", "
                    << "None, "
                    << fitter->second->wmean(MuonResiduals5DOFFitter::kResSlope, MuonResiduals5DOFFitter::kRedChi2, 10.)
                    << ", "
                    << "None, " << fitter->second->stdev(MuonResiduals5DOFFitter::kResid, 30.) << ", "
                    << "None, " << fitter->second->stdev(MuonResiduals5DOFFitter::kResSlope, 20.) << ", "
                    << "None, " << fitter->second->stdev(MuonResiduals5DOFFitter::kResid, 15.) << ", "
                    << "None, " << fitter->second->stdev(MuonResiduals5DOFFitter::kResSlope, 10.) << ", "
                    << "None)" << std::endl;
 
             std::stringstream namesimple_x, namesimple_dxdz, nameweighted_x, nameweighted_dxdz;
             namesimple_x << cname << "_simple_x";
             namesimple_dxdz << cname << "_simple_dxdz";
             nameweighted_x << cname << "_weighted_x";
             nameweighted_dxdz << cname << "_weighted_dxdz";
 
             fitter->second->plotsimple(namesimple_x.str(), &rootDirectory, MuonResiduals5DOFFitter::kResid, 10.);
             fitter->second->plotsimple(
                 namesimple_dxdz.str(), &rootDirectory, MuonResiduals5DOFFitter::kResSlope, 1000.);
 
             fitter->second->plotweighted(nameweighted_x.str(),
                                          &rootDirectory,
                                          MuonResiduals5DOFFitter::kResid,
                                          MuonResiduals5DOFFitter::kRedChi2,
                                          10.);
             fitter->second->plotweighted(nameweighted_dxdz.str(),
                                          &rootDirectory,
                                          MuonResiduals5DOFFitter::kResSlope,
                                          MuonResiduals5DOFFitter::kRedChi2,
                                          1000.);
           }
 
           if (successful_fit) {
             if (align_x)
               params[paramIndex[0]] = deltax_value;
             if (align_z)
               params[paramIndex[2]] = deltaz_value;
             if (align_phix)
               params[paramIndex[3]] = deltaphix_value;
             if (align_phiy)
               params[paramIndex[4]] = deltaphiy_value;
             if (align_phiz)
               params[paramIndex[5]] = deltaphiz_value;
           }
         }  // end if 5DOF
 
         else if (fitter->second->type() == MuonResidualsFitter::k6DOF) {
           if (m_debug)
             std::cout << "***** loop over alignables k6DOF" << std::endl;
 
           double deltax_value = fitter->second->value(MuonResiduals6DOFFitter::kAlignX);
           double deltax_error = fitter->second->errorerror(MuonResiduals6DOFFitter::kAlignX);
           double deltax_antisym = fitter->second->antisym(MuonResiduals6DOFFitter::kAlignX);
 
           double deltay_value = fitter->second->value(MuonResiduals6DOFFitter::kAlignY);
           double deltay_error = fitter->second->errorerror(MuonResiduals6DOFFitter::kAlignY);
           double deltay_antisym = fitter->second->antisym(MuonResiduals6DOFFitter::kAlignY);
 
           double deltaz_value = fitter->second->value(MuonResiduals6DOFFitter::kAlignZ);
           double deltaz_error = fitter->second->errorerror(MuonResiduals6DOFFitter::kAlignZ);
           double deltaz_antisym = fitter->second->antisym(MuonResiduals6DOFFitter::kAlignZ);
 
           double deltaphix_value = fitter->second->value(MuonResiduals6DOFFitter::kAlignPhiX);
           double deltaphix_error = fitter->second->errorerror(MuonResiduals6DOFFitter::kAlignPhiX);
           double deltaphix_antisym = fitter->second->antisym(MuonResiduals6DOFFitter::kAlignPhiX);
 
           double deltaphiy_value = fitter->second->value(MuonResiduals6DOFFitter::kAlignPhiY);
           double deltaphiy_error = fitter->second->errorerror(MuonResiduals6DOFFitter::kAlignPhiY);
           double deltaphiy_antisym = fitter->second->antisym(MuonResiduals6DOFFitter::kAlignPhiY);
 
           double deltaphiz_value = fitter->second->value(MuonResiduals6DOFFitter::kAlignPhiZ);
           double deltaphiz_error = fitter->second->errorerror(MuonResiduals6DOFFitter::kAlignPhiZ);
           double deltaphiz_antisym = fitter->second->antisym(MuonResiduals6DOFFitter::kAlignPhiZ);
 
           double sigmax_value = fitter->second->value(MuonResiduals6DOFFitter::kResidXSigma);
           double sigmax_error = fitter->second->errorerror(MuonResiduals6DOFFitter::kResidXSigma);
           double sigmax_antisym = fitter->second->antisym(MuonResiduals6DOFFitter::kResidXSigma);
 
           double sigmay_value = fitter->second->value(MuonResiduals6DOFFitter::kResidYSigma);
           double sigmay_error = fitter->second->errorerror(MuonResiduals6DOFFitter::kResidYSigma);
           double sigmay_antisym = fitter->second->antisym(MuonResiduals6DOFFitter::kResidYSigma);
 
           double sigmadxdz_value = fitter->second->value(MuonResiduals6DOFFitter::kResSlopeXSigma);
           double sigmadxdz_error = fitter->second->errorerror(MuonResiduals6DOFFitter::kResSlopeXSigma);
           double sigmadxdz_antisym = fitter->second->antisym(MuonResiduals6DOFFitter::kResSlopeXSigma);
 
           double sigmadydz_value = fitter->second->value(MuonResiduals6DOFFitter::kResSlopeYSigma);
           double sigmadydz_error = fitter->second->errorerror(MuonResiduals6DOFFitter::kResSlopeYSigma);
           double sigmadydz_antisym = fitter->second->antisym(MuonResiduals6DOFFitter::kResSlopeYSigma);
 
           double gammax_value, gammax_error, gammax_antisym, gammay_value, gammay_error, gammay_antisym,
               gammadxdz_value, gammadxdz_error, gammadxdz_antisym, gammadydz_value, gammadydz_error, gammadydz_antisym;
           gammax_value = gammax_error = gammax_antisym = gammay_value = gammay_error = gammay_antisym =
               gammadxdz_value = gammadxdz_error = gammadxdz_antisym = gammadydz_value = gammadydz_error =
                   gammadydz_antisym = 0.;
           if (fitter->second->residualsModel() != MuonResidualsFitter::kPureGaussian &&
               fitter->second->residualsModel() != MuonResidualsFitter::kPureGaussian2D &&
               fitter->second->residualsModel() != MuonResidualsFitter::kGaussPowerTails) {
             gammax_value = fitter->second->value(MuonResiduals6DOFFitter::kResidXGamma);
             gammax_error = fitter->second->errorerror(MuonResiduals6DOFFitter::kResidXGamma);
             gammax_antisym = fitter->second->antisym(MuonResiduals6DOFFitter::kResidXGamma);
 
             gammay_value = fitter->second->value(MuonResiduals6DOFFitter::kResidYGamma);
             gammay_error = fitter->second->errorerror(MuonResiduals6DOFFitter::kResidYGamma);
             gammay_antisym = fitter->second->antisym(MuonResiduals6DOFFitter::kResidYGamma);
 
             gammadxdz_value = fitter->second->value(MuonResiduals6DOFFitter::kResSlopeXGamma);
             gammadxdz_error = fitter->second->errorerror(MuonResiduals6DOFFitter::kResSlopeXGamma);
             gammadxdz_antisym = fitter->second->antisym(MuonResiduals6DOFFitter::kResSlopeXGamma);
 
             gammadydz_value = fitter->second->value(MuonResiduals6DOFFitter::kResSlopeYGamma);
             gammadydz_error = fitter->second->errorerror(MuonResiduals6DOFFitter::kResSlopeYGamma);
             gammadydz_antisym = fitter->second->antisym(MuonResiduals6DOFFitter::kResSlopeYGamma);
           }
 
           if (writeReport) {
             report << "reports[-1].fittype = \"6DOF\"" << std::endl;
             report << "reports[-1].add_parameters(ValErr(" << deltax_value << ", " << deltax_error << ", "
                    << deltax_antisym << "), \\" << std::endl
                    << "                           ValErr(" << deltay_value << ", " << deltay_error << ", "
                    << deltay_antisym << "), \\" << std::endl
                    << "                           ValErr(" << deltaz_value << ", " << deltaz_error << ", "
                    << deltaz_antisym << "), \\" << std::endl
                    << "                           ValErr(" << deltaphix_value << ", " << deltaphix_error << ", "
                    << deltaphix_antisym << "), \\" << std::endl
                    << "                           ValErr(" << deltaphiy_value << ", " << deltaphiy_error << ", "
                    << deltaphiy_antisym << "), \\" << std::endl
                    << "                           ValErr(" << deltaphiz_value << ", " << deltaphiz_error << ", "
                    << deltaphiz_antisym << "), \\" << std::endl
                    << "                           " << loglikelihood << ", " << numsegments << ", " << sumofweights
                    << ", " << redchi2 << ")" << std::endl;
             report << "reports[-1].sigmax = ValErr(" << sigmax_value << ", " << sigmax_error << ", " << sigmax_antisym
                    << ")" << std::endl;
             report << "reports[-1].sigmay = ValErr(" << sigmay_value << ", " << sigmay_error << ", " << sigmay_antisym
                    << ")" << std::endl;
             report << "reports[-1].sigmadxdz = ValErr(" << sigmadxdz_value << ", " << sigmadxdz_error << ", "
                    << sigmadxdz_antisym << ")" << std::endl;
             report << "reports[-1].sigmadydz = ValErr(" << sigmadydz_value << ", " << sigmadydz_error << ", "
                    << sigmadydz_antisym << ")" << std::endl;
             if (fitter->second->residualsModel() != MuonResidualsFitter::kPureGaussian &&
                 fitter->second->residualsModel() != MuonResidualsFitter::kPureGaussian2D &&
                 fitter->second->residualsModel() != MuonResidualsFitter::kGaussPowerTails) {
               report << "reports[-1].gammax = ValErr(" << gammax_value << ", " << gammax_error << ", " << gammax_antisym
                      << ")" << std::endl;
               report << "reports[-1].gammay = ValErr(" << gammay_value << ", " << gammay_error << ", " << gammay_antisym
                      << ")" << std::endl;
               report << "reports[-1].gammadxdz = ValErr(" << gammadxdz_value << ", " << gammadxdz_error << ", "
                      << gammadxdz_antisym << ")" << std::endl;
               report << "reports[-1].gammadydz = ValErr(" << gammadydz_value << ", " << gammadydz_error << ", "
                      << gammadydz_antisym << ")" << std::endl;
             }
 
             report << "reports[-1].add_stats(" << fitter->second->median(MuonResiduals6DOFFitter::kResidX) << ", "
                    << fitter->second->median(MuonResiduals6DOFFitter::kResidY) << ", "
                    << fitter->second->median(MuonResiduals6DOFFitter::kResSlopeX) << ", "
                    << fitter->second->median(MuonResiduals6DOFFitter::kResSlopeY) << ", "
                    << fitter->second->mean(MuonResiduals6DOFFitter::kResidX, 30.) << ", "
                    << fitter->second->mean(MuonResiduals6DOFFitter::kResidY, 30.) << ", "
                    << fitter->second->mean(MuonResiduals6DOFFitter::kResSlopeX, 20.) << ", "
                    << fitter->second->mean(MuonResiduals6DOFFitter::kResSlopeY, 50.) << ", "
                    << fitter->second->mean(MuonResiduals6DOFFitter::kResidX, 15.) << ", "
                    << fitter->second->mean(MuonResiduals6DOFFitter::kResidY, 15.) << ", "
                    << fitter->second->mean(MuonResiduals6DOFFitter::kResSlopeX, 10.) << ", "
                    << fitter->second->mean(MuonResiduals6DOFFitter::kResSlopeY, 25.) << ", "
                    << fitter->second->wmean(MuonResiduals6DOFFitter::kResidX, MuonResiduals6DOFFitter::kRedChi2, 30.)
                    << ", "
                    << fitter->second->wmean(MuonResiduals6DOFFitter::kResidY, MuonResiduals6DOFFitter::kRedChi2, 30.)
                    << ", "
                    << fitter->second->wmean(MuonResiduals6DOFFitter::kResSlopeX, MuonResiduals6DOFFitter::kRedChi2, 20.)
                    << ", "
                    << fitter->second->wmean(MuonResiduals6DOFFitter::kResSlopeY, MuonResiduals6DOFFitter::kRedChi2, 50.)
                    << ", "
                    << fitter->second->wmean(MuonResiduals6DOFFitter::kResidX, MuonResiduals6DOFFitter::kRedChi2, 15.)
                    << ", "
                    << fitter->second->wmean(MuonResiduals6DOFFitter::kResidY, MuonResiduals6DOFFitter::kRedChi2, 15.)
                    << ", "
                    << fitter->second->wmean(MuonResiduals6DOFFitter::kResSlopeX, MuonResiduals6DOFFitter::kRedChi2, 10.)
                    << ", "
                    << fitter->second->wmean(MuonResiduals6DOFFitter::kResSlopeY, MuonResiduals6DOFFitter::kRedChi2, 25.)
                    << ", " << fitter->second->stdev(MuonResiduals6DOFFitter::kResidX, 30.) << ", "
                    << fitter->second->stdev(MuonResiduals6DOFFitter::kResidY, 30.) << ", "
                    << fitter->second->stdev(MuonResiduals6DOFFitter::kResSlopeX, 20.) << ", "
                    << fitter->second->stdev(MuonResiduals6DOFFitter::kResSlopeY, 50.) << ", "
                    << fitter->second->stdev(MuonResiduals6DOFFitter::kResidX, 15.) << ", "
                    << fitter->second->stdev(MuonResiduals6DOFFitter::kResidY, 15.) << ", "
                    << fitter->second->stdev(MuonResiduals6DOFFitter::kResSlopeX, 10.) << ", "
                    << fitter->second->stdev(MuonResiduals6DOFFitter::kResSlopeY, 25.) << ")" << std::endl;
 
             std::stringstream namesimple_x, namesimple_y, namesimple_dxdz, namesimple_dydz, nameweighted_x,
                 nameweighted_y, nameweighted_dxdz, nameweighted_dydz;
             namesimple_x << cname << "_simple_x";
             namesimple_y << cname << "_simple_y";
             namesimple_dxdz << cname << "_simple_dxdz";
             namesimple_dydz << cname << "_simple_dydz";
             nameweighted_x << cname << "_weighted_x";
             nameweighted_y << cname << "_weighted_y";
             nameweighted_dxdz << cname << "_weighted_dxdz";
             nameweighted_dydz << cname << "_weighted_dydz";
 
             fitter->second->plotsimple(namesimple_x.str(), &rootDirectory, MuonResiduals6DOFFitter::kResidX, 10.);
             fitter->second->plotsimple(namesimple_y.str(), &rootDirectory, MuonResiduals6DOFFitter::kResidY, 10.);
             fitter->second->plotsimple(
                 namesimple_dxdz.str(), &rootDirectory, MuonResiduals6DOFFitter::kResSlopeX, 1000.);
             fitter->second->plotsimple(
                 namesimple_dydz.str(), &rootDirectory, MuonResiduals6DOFFitter::kResSlopeY, 1000.);
 
             fitter->second->plotweighted(nameweighted_x.str(),
                                          &rootDirectory,
                                          MuonResiduals6DOFFitter::kResidX,
                                          MuonResiduals6DOFFitter::kRedChi2,
                                          10.);
             fitter->second->plotweighted(nameweighted_y.str(),
                                          &rootDirectory,
                                          MuonResiduals6DOFFitter::kResidY,
                                          MuonResiduals6DOFFitter::kRedChi2,
                                          10.);
             fitter->second->plotweighted(nameweighted_dxdz.str(),
                                          &rootDirectory,
                                          MuonResiduals6DOFFitter::kResSlopeX,
                                          MuonResiduals6DOFFitter::kRedChi2,
                                          1000.);
             fitter->second->plotweighted(nameweighted_dydz.str(),
                                          &rootDirectory,
                                          MuonResiduals6DOFFitter::kResSlopeY,
                                          MuonResiduals6DOFFitter::kRedChi2,
                                          1000.);
           }
 
           if (successful_fit) {
             if (align_x)
               params[paramIndex[0]] = deltax_value;
             if (align_y)
               params[paramIndex[1]] = deltay_value;
             if (align_z)
               params[paramIndex[2]] = deltaz_value;
             if (align_phix)
               params[paramIndex[3]] = deltaphix_value;
             if (align_phiy)
               params[paramIndex[4]] = deltaphiy_value;
             if (align_phiz)
               params[paramIndex[5]] = deltaphiz_value;
           }
         }  // end if 6DOF
 
         else if (fitter->second->type() == MuonResidualsFitter::k6DOFrphi) {
           if (m_debug)
             std::cout << "***** loop over alignables k6DOFrphi" << std::endl;
 
           double deltax_value = fitter->second->value(MuonResiduals6DOFrphiFitter::kAlignX);
           double deltax_error = fitter->second->errorerror(MuonResiduals6DOFrphiFitter::kAlignX);
           double deltax_antisym = fitter->second->antisym(MuonResiduals6DOFrphiFitter::kAlignX);
 
           double deltay_value = fitter->second->value(MuonResiduals6DOFrphiFitter::kAlignY);
           double deltay_error = fitter->second->errorerror(MuonResiduals6DOFrphiFitter::kAlignY);
           double deltay_antisym = fitter->second->antisym(MuonResiduals6DOFrphiFitter::kAlignY);
 
           double deltaz_value = fitter->second->value(MuonResiduals6DOFrphiFitter::kAlignZ);
           double deltaz_error = fitter->second->errorerror(MuonResiduals6DOFrphiFitter::kAlignZ);
           double deltaz_antisym = fitter->second->antisym(MuonResiduals6DOFrphiFitter::kAlignZ);
 
           double deltaphix_value = fitter->second->value(MuonResiduals6DOFrphiFitter::kAlignPhiX);
           double deltaphix_error = fitter->second->errorerror(MuonResiduals6DOFrphiFitter::kAlignPhiX);
           double deltaphix_antisym = fitter->second->antisym(MuonResiduals6DOFrphiFitter::kAlignPhiX);
 
           double deltaphiy_value = fitter->second->value(MuonResiduals6DOFrphiFitter::kAlignPhiY);
           double deltaphiy_error = fitter->second->errorerror(MuonResiduals6DOFrphiFitter::kAlignPhiY);
           double deltaphiy_antisym = fitter->second->antisym(MuonResiduals6DOFrphiFitter::kAlignPhiY);
 
           double deltaphiz_value = fitter->second->value(MuonResiduals6DOFrphiFitter::kAlignPhiZ);
           double deltaphiz_error = fitter->second->errorerror(MuonResiduals6DOFrphiFitter::kAlignPhiZ);
           double deltaphiz_antisym = fitter->second->antisym(MuonResiduals6DOFrphiFitter::kAlignPhiZ);
 
           double sigmaresid_value = fitter->second->value(MuonResiduals6DOFrphiFitter::kResidSigma);
           double sigmaresid_error = fitter->second->errorerror(MuonResiduals6DOFrphiFitter::kResidSigma);
           double sigmaresid_antisym = fitter->second->antisym(MuonResiduals6DOFrphiFitter::kResidSigma);
 
           double sigmaresslope_value = fitter->second->value(MuonResiduals6DOFrphiFitter::kResSlopeSigma);
           double sigmaresslope_error = fitter->second->errorerror(MuonResiduals6DOFrphiFitter::kResSlopeSigma);
           double sigmaresslope_antisym = fitter->second->antisym(MuonResiduals6DOFrphiFitter::kResSlopeSigma);
 
           double gammaresid_value, gammaresid_error, gammaresid_antisym, gammaresslope_value, gammaresslope_error,
               gammaresslope_antisym;
           gammaresid_value = gammaresid_error = gammaresid_antisym = gammaresslope_value = gammaresslope_error =
               gammaresslope_antisym = 0.;
           if (fitter->second->residualsModel() != MuonResidualsFitter::kPureGaussian &&
               fitter->second->residualsModel() != MuonResidualsFitter::kPureGaussian2D &&
               fitter->second->residualsModel() != MuonResidualsFitter::kGaussPowerTails) {
             gammaresid_value = fitter->second->value(MuonResiduals6DOFrphiFitter::kResidGamma);
             gammaresid_error = fitter->second->errorerror(MuonResiduals6DOFrphiFitter::kResidGamma);
             gammaresid_antisym = fitter->second->antisym(MuonResiduals6DOFrphiFitter::kResidGamma);
 
             gammaresslope_value = fitter->second->value(MuonResiduals6DOFrphiFitter::kResSlopeGamma);
             gammaresslope_error = fitter->second->errorerror(MuonResiduals6DOFrphiFitter::kResSlopeGamma);
             gammaresslope_antisym = fitter->second->antisym(MuonResiduals6DOFrphiFitter::kResSlopeGamma);
           }
 
           if (writeReport) {
             report << "reports[-1].fittype = \"6DOFrphi\"" << std::endl;
             report << "reports[-1].add_parameters(ValErr(" << deltax_value << ", " << deltax_error << ", "
                    << deltax_antisym << "), \\" << std::endl
                    << "                           ValErr(" << deltay_value << ", " << deltay_error << ", "
                    << deltay_antisym << "), \\" << std::endl
                    << "                           ValErr(" << deltaz_value << ", " << deltaz_error << ", "
                    << deltaz_antisym << "), \\" << std::endl
                    << "                           ValErr(" << deltaphix_value << ", " << deltaphix_error << ", "
                    << deltaphix_antisym << "), \\" << std::endl
                    << "                           ValErr(" << deltaphiy_value << ", " << deltaphiy_error << ", "
                    << deltaphiy_antisym << "), \\" << std::endl
                    << "                           ValErr(" << deltaphiz_value << ", " << deltaphiz_error << ", "
                    << deltaphiz_antisym << "), \\" << std::endl
                    << "                           " << loglikelihood << ", " << numsegments << ", " << sumofweights
                    << ", " << redchi2 << ")" << std::endl;
             report << "reports[-1].sigmaresid = ValErr(" << sigmaresid_value << ", " << sigmaresid_error << ", "
                    << sigmaresid_antisym << ")" << std::endl;
             report << "reports[-1].sigmaresslope = ValErr(" << sigmaresslope_value << ", " << sigmaresslope_error
                    << ", " << sigmaresslope_antisym << ")" << std::endl;
             if (fitter->second->residualsModel() != MuonResidualsFitter::kPureGaussian &&
                 fitter->second->residualsModel() != MuonResidualsFitter::kPureGaussian2D &&
                 fitter->second->residualsModel() != MuonResidualsFitter::kGaussPowerTails) {
               report << "reports[-1].gammaresid = ValErr(" << gammaresid_value << ", " << gammaresid_error << ", "
                      << gammaresid_antisym << ")" << std::endl;
               report << "reports[-1].gammaresslope = ValErr(" << gammaresslope_value << ", " << gammaresslope_error
                      << ", " << gammaresslope_antisym << ")" << std::endl;
             }
 
             report << "reports[-1].add_stats(" << fitter->second->median(MuonResiduals6DOFrphiFitter::kResid) << ", "
                    << "None, " << fitter->second->median(MuonResiduals6DOFrphiFitter::kResSlope) << ", "
                    << "None, " << fitter->second->mean(MuonResiduals6DOFrphiFitter::kResid, 30.) << ", "
                    << "None, " << fitter->second->mean(MuonResiduals6DOFrphiFitter::kResSlope, 20.) << ", "
                    << "None, " << fitter->second->mean(MuonResiduals6DOFrphiFitter::kResid, 15.) << ", "
                    << "None, " << fitter->second->mean(MuonResiduals6DOFrphiFitter::kResSlope, 10.) << ", "
                    << "None, "
                    << fitter->second->wmean(
                           MuonResiduals6DOFrphiFitter::kResid, MuonResiduals6DOFrphiFitter::kRedChi2, 30.)
                    << ", "
                    << "None, "
                    << fitter->second->wmean(
                           MuonResiduals6DOFrphiFitter::kResSlope, MuonResiduals6DOFrphiFitter::kRedChi2, 20.)
                    << ", "
                    << "None, "
                    << fitter->second->wmean(
                           MuonResiduals6DOFrphiFitter::kResid, MuonResiduals6DOFrphiFitter::kRedChi2, 15.)
                    << ", "
                    << "None, "
                    << fitter->second->wmean(
                           MuonResiduals6DOFrphiFitter::kResSlope, MuonResiduals6DOFrphiFitter::kRedChi2, 10.)
                    << ", "
                    << "None, " << fitter->second->stdev(MuonResiduals6DOFrphiFitter::kResid, 30.) << ", "
                    << "None, " << fitter->second->stdev(MuonResiduals6DOFrphiFitter::kResSlope, 20.) << ", "
                    << "None, " << fitter->second->stdev(MuonResiduals6DOFrphiFitter::kResid, 15.) << ", "
                    << "None, " << fitter->second->stdev(MuonResiduals6DOFrphiFitter::kResSlope, 10.) << ", "
                    << "None)" << std::endl;
 
             std::stringstream namesimple_x, namesimple_dxdz, nameweighted_x, nameweighted_dxdz;
             namesimple_x << cname << "_simple_x";
             namesimple_dxdz << cname << "_simple_dxdz";
             nameweighted_x << cname << "_weighted_x";
             nameweighted_dxdz << cname << "_weighted_dxdz";
 
             fitter->second->plotsimple(namesimple_x.str(), &rootDirectory, MuonResiduals6DOFrphiFitter::kResid, 10.);
             fitter->second->plotsimple(
                 namesimple_dxdz.str(), &rootDirectory, MuonResiduals6DOFrphiFitter::kResSlope, 1000.);
 
             fitter->second->plotweighted(nameweighted_x.str(),
                                          &rootDirectory,
                                          MuonResiduals6DOFrphiFitter::kResid,
                                          MuonResiduals6DOFrphiFitter::kRedChi2,
                                          10.);
             fitter->second->plotweighted(nameweighted_dxdz.str(),
                                          &rootDirectory,
                                          MuonResiduals6DOFrphiFitter::kResSlope,
                                          MuonResiduals6DOFrphiFitter::kRedChi2,
                                          1000.);
           }
 
           if (successful_fit) {
             if (align_x)
               params[paramIndex[0]] = deltax_value;
             if (align_y)
               params[paramIndex[1]] = deltay_value;
             if (align_z)
               params[paramIndex[2]] = deltaz_value;
             if (align_phix)
               params[paramIndex[3]] = deltaphix_value;
             if (align_phiy)
               params[paramIndex[4]] = deltaphiy_value;
             if (align_phiz)
               params[paramIndex[5]] = deltaphiz_value;
           }
         }  // end if 6DOFrphi
 
         if (successful_fit) {
           align::Alignables oneortwo;
           oneortwo.push_back(ali);
           if (thisali != ali)
             oneortwo.push_back(thisali);
           m_alignmentParameterStore->setAlignmentPositionError(oneortwo, 0., 0.);
         } else {
           if (m_debug)
             std::cout << "MINUIT fit failed!" << std::endl;
           if (writeReport) {
             report << "reports[-1].status = \"MINUITFAIL\"" << std::endl;
           }
 
           for (int i = 0; i < numParams; i++)
             cov[i][i] = 1000.;
 
           align::Alignables oneortwo;
           oneortwo.push_back(ali);
           if (thisali != ali)
             oneortwo.push_back(thisali);
           m_alignmentParameterStore->setAlignmentPositionError(oneortwo, 1000., 0.);
         }
       } else {  // too few hits
         if (m_debug)
           std::cout << "Too few hits!" << std::endl;
         if (writeReport) {
           report << "reports[-1].status = \"TOOFEWHITS\"" << std::endl;
         }
 
         for (int i = 0; i < numParams; i++)
           cov[i][i] = 1000.;
 
         align::Alignables oneortwo;
         oneortwo.push_back(ali);
         if (thisali != ali)
           oneortwo.push_back(thisali);
         m_alignmentParameterStore->setAlignmentPositionError(oneortwo, 1000., 0.);
       }
 
       AlignmentParameters* parnew = ali->alignmentParameters()->cloneFromSelected(params, cov);
       ali->setAlignmentParameters(parnew);
       m_alignmentParameterStore->applyParameters(ali);
       ali->alignmentParameters()->setValid(true);
 
       if (m_debug)
         std::cout << cname << " fittime= " << stop_watch.CpuTime() << " sec" << std::endl;
     }  // end we have a fitter for this alignable
 
     if (writeReport)
       report << std::endl;
 
   }  // end loop over alignables
 
   if (writeReport)
     report.close();
 }

◆ initialize()

void MuonAlignmentFromReference::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 308 of file MuonAlignmentFromReference.cc.

                                                                                               {
   if (alignableMuon == nullptr)
     throw cms::Exception("MuonAlignmentFromReference") << "doMuon must be set to True" << std::endl;
 
   m_alignableNavigator = new AlignableNavigator(alignableMuon);
   m_alignmentParameterStore = alignmentParameterStore;
   m_alignables = m_alignmentParameterStore->alignables();
 
   int residualsModel;
   if (m_residualsModel == std::string("pureGaussian"))
     residualsModel = MuonResidualsFitter::kPureGaussian;
   else if (m_residualsModel == std::string("pureGaussian2D"))
     residualsModel = MuonResidualsFitter::kPureGaussian2D;
   else if (m_residualsModel == std::string("powerLawTails"))
     residualsModel = MuonResidualsFitter::kPowerLawTails;
   else if (m_residualsModel == std::string("ROOTVoigt"))
     residualsModel = MuonResidualsFitter::kROOTVoigt;
   else if (m_residualsModel == std::string("GaussPowerTails"))
     residualsModel = MuonResidualsFitter::kGaussPowerTails;
   else
     throw cms::Exception("MuonAlignmentFromReference")
         << "unrecognized residualsModel: \"" << m_residualsModel << "\"" << std::endl;
 
   int useResiduals;
   if (m_useResiduals == std::string("1111"))
     useResiduals = MuonResidualsFitter::k1111;
   else if (m_useResiduals == std::string("1110"))
     useResiduals = MuonResidualsFitter::k1110;
   else if (m_useResiduals == std::string("1100"))
     useResiduals = MuonResidualsFitter::k1100;
   else if (m_useResiduals == std::string("1000"))
     useResiduals = MuonResidualsFitter::k1000;
   else if (m_useResiduals == std::string("1010"))
     useResiduals = MuonResidualsFitter::k1010;
   else if (m_useResiduals == std::string("0010"))
     useResiduals = MuonResidualsFitter::k0010;
   else
     throw cms::Exception("MuonAlignmentFromReference")
         << "unrecognized useResiduals: \"" << m_useResiduals << "\"" << std::endl;
 
   const CSCGeometry* cscGeometry = &iSetup.getData(m_cscGeometryToken);
 
   // set up the MuonResidualsFitters (which also collect residuals for fitting)
   m_me11map.clear();
   m_fitters.clear();
   m_indexes.clear();
   m_fitterOrder.clear();
 
   for (const auto& ali : m_alignables) {
     bool made_fitter = false;
 
     // fitters for DT
     if (ali->alignableObjectId() == align::AlignableDTChamber) {
       DTChamberId id(ali->geomDetId().rawId());
 
       if (id.station() == 4) {
         m_fitters[ali] = new MuonResidualsTwoBin(
             m_twoBin,
             new MuonResiduals5DOFFitter(residualsModel, m_minAlignmentHits, useResiduals, m_weightAlignment),
             new MuonResiduals5DOFFitter(residualsModel, m_minAlignmentHits, useResiduals, m_weightAlignment));
         made_fitter = true;
       } else {
         m_fitters[ali] = new MuonResidualsTwoBin(
             m_twoBin,
             new MuonResiduals6DOFFitter(residualsModel, m_minAlignmentHits, useResiduals, m_weightAlignment),
             new MuonResiduals6DOFFitter(residualsModel, m_minAlignmentHits, useResiduals, m_weightAlignment));
         made_fitter = true;
       }
     }
 
     // fitters for CSC
     else if (ali->alignableObjectId() == align::AlignableCSCChamber) {
       auto thisali = ali;
       CSCDetId id(ali->geomDetId().rawId());
 
       // take care of ME1/1a
       if (m_combineME11 && id.station() == 1 && id.ring() == 4) {
         CSCDetId pairid(id.endcap(), 1, 1, id.chamber());
 
         for (const auto& ali2 : m_alignables) {
           if (ali2->alignableObjectId() == align::AlignableCSCChamber && ali2->geomDetId().rawId() == pairid.rawId()) {
             thisali = ali2;
             break;
           }
         }
         m_me11map[ali] = thisali;  // points from each ME1/4 chamber to the corresponding ME1/1 chamber
       }
 
       if (thisali == ali)  // don't make fitters for ME1/4; they get taken care of in ME1/1
       {
         m_fitters[ali] = new MuonResidualsTwoBin(
             m_twoBin,
             new MuonResiduals6DOFrphiFitter(
                 residualsModel, m_minAlignmentHits, useResiduals, cscGeometry, m_weightAlignment),
             new MuonResiduals6DOFrphiFitter(
                 residualsModel, m_minAlignmentHits, useResiduals, cscGeometry, m_weightAlignment));
         made_fitter = true;
       }
     }
 
     else {
       throw cms::Exception("MuonAlignmentFromReference")
           << "only DTChambers and CSCChambers can be aligned with this module" << std::endl;
     }
 
     if (made_fitter) {
       m_fitters[ali]->setStrategy(m_strategy);
 
       int index = ali->geomDetId().rawId();
       m_indexes.push_back(index);
       m_fitterOrder[index] = m_fitters[ali];
     }
   }  // end loop over chambers chosen for alignment
 
   // cannonical order of fitters in the file
   std::sort(m_indexes.begin(), m_indexes.end());
 
   // de-weight all chambers but the reference
   const auto& all_DT_chambers = alignableMuon->DTChambers();
   const auto& all_CSC_chambers = alignableMuon->CSCChambers();
   align::Alignables reference;
   if (!m_reference.empty())
     parseReference(reference, all_DT_chambers, all_CSC_chambers);
 
   alignmentParameterStore->setAlignmentPositionError(all_DT_chambers, 100000000., 0.);
   alignmentParameterStore->setAlignmentPositionError(all_CSC_chambers, 100000000., 0.);
   alignmentParameterStore->setAlignmentPositionError(reference, 0., 0.);
 }

◆ number()

int MuonAlignmentFromReference::number ( std::string s )

private

Definition at line 302 of file MuonAlignmentFromReference.cc.

References cms::cuda::assert(), numeric(), and alignCSCRings::s.

Referenced by parseReference().

                                                   {
   if (!numeric(s))
     assert(false);
   return atoi(s.c_str());
 }

◆ numeric()

bool MuonAlignmentFromReference::numeric ( std::string s )

private

Definition at line 300 of file MuonAlignmentFromReference.cc.

References alignCSCRings::s.

Referenced by number(), and parseReference().

300 { return s.length() == 1 && std::isdigit(s[0]); }

alignCSCRings.s

s

Definition: alignCSCRings.py:92

◆ parseReference()

void MuonAlignmentFromReference::parseReference	(	align::Alignables &	reference,
		const align::Alignables &	all_DT_chambers,
		const align::Alignables &	all_CSC_chambers
	)

private

Definition at line 1837 of file MuonAlignmentFromReference.cc.

References funct::abs(), Reference_intrackfit_cff::barrel, relativeConstraints::chamber, makeMuonMisalignmentScenario::endcap, Exception, l1ctLayer2EG_cff::id, m_reference, Skims_PA_cff::name, number(), numeric(), reco::HaloData::plus, RecoTauValidation_cfi::reference, relativeConstraints::ring, relativeConstraints::station, AlCaHLTBitMon_QueryRunRegistry::string, dtNoiseDBValidation_cfg::trial, and makeMuonMisalignmentScenario::wheel.

Referenced by initialize().

                                                                                          {
   std::map<Alignable*, bool> already_seen;
 
   for (std::vector<std::string>::const_iterator name = m_reference.begin(); name != m_reference.end(); ++name) {
     bool parsing_error = false;
 
     bool barrel = (name->substr(0, 2) == std::string("MB"));
     bool endcap = (name->substr(0, 2) == std::string("ME"));
     if (!barrel && !endcap)
       parsing_error = true;
 
     if (!parsing_error && barrel) {
       int index = 2;
       if (name->substr(index, 1) == std::string(" "))
         index++;
 
       bool plus = true;
       if (name->substr(index, 1) == std::string("+")) {
         plus = true;
         index++;
       } else if (name->substr(index, 1) == std::string("-")) {
         plus = false;
         index++;
       } else if (numeric(name->substr(index, 1))) {
       } else
         parsing_error = true;
 
       int wheel = 0;
       bool wheel_digit = false;
       while (!parsing_error && numeric(name->substr(index, 1))) {
         wheel *= 10;
         wheel += number(name->substr(index, 1));
         wheel_digit = true;
         index++;
       }
       if (!plus)
         wheel *= -1;
       if (!wheel_digit)
         parsing_error = true;
 
       if (name->substr(index, 1) != std::string(" "))
         parsing_error = true;
       index++;
 
       int station = 0;
       bool station_digit = false;
       while (!parsing_error && numeric(name->substr(index, 1))) {
         station *= 10;
         station += number(name->substr(index, 1));
         station_digit = true;
         index++;
       }
       if (!station_digit)
         parsing_error = true;
 
       if (name->substr(index, 1) != std::string(" "))
         parsing_error = true;
       index++;
 
       int sector = 0;
       bool sector_digit = false;
       while (!parsing_error && numeric(name->substr(index, 1))) {
         sector *= 10;
         sector += number(name->substr(index, 1));
         sector_digit = true;
         index++;
       }
       if (!sector_digit)
         parsing_error = true;
 
       if (!parsing_error) {
         bool no_such_chamber = false;
 
         if (wheel < -2 || wheel > 2)
           no_such_chamber = true;
         if (station < 1 || station > 4)
           no_such_chamber = true;
         if (station == 4 && (sector < 1 || sector > 14))
           no_such_chamber = true;
         if (station < 4 && (sector < 1 || sector > 12))
           no_such_chamber = true;
 
         if (no_such_chamber)
           throw cms::Exception("MuonAlignmentFromReference")
               << "reference chamber doesn't exist: " << (*name) << std::endl;
 
         DTChamberId id(wheel, station, sector);
         for (const auto& ali : all_DT_chambers) {
           if (ali->geomDetId().rawId() == id.rawId()) {
             std::map<Alignable*, bool>::const_iterator trial = already_seen.find(ali);
             if (trial == already_seen.end()) {
               reference.push_back(ali);
               already_seen[ali] = true;
             }
           }
         }
       }  // if (!parsing_error)
     }
 
     if (!parsing_error && endcap) {
       int index = 2;
       if (name->substr(index, 1) == std::string(" "))
         index++;
 
       bool plus = true;
       if (name->substr(index, 1) == std::string("+")) {
         plus = true;
         index++;
       } else if (name->substr(index, 1) == std::string("-")) {
         plus = false;
         index++;
       } else if (numeric(name->substr(index, 1))) {
       } else
         parsing_error = true;
 
       int station = 0;
       bool station_digit = false;
       while (!parsing_error && numeric(name->substr(index, 1))) {
         station *= 10;
         station += number(name->substr(index, 1));
         station_digit = true;
         index++;
       }
       if (!plus)
         station *= -1;
       if (!station_digit)
         parsing_error = true;
 
       if (name->substr(index, 1) != std::string("/"))
         parsing_error = true;
       index++;
 
       int ring = 0;
       bool ring_digit = false;
       while (!parsing_error && numeric(name->substr(index, 1))) {
         ring *= 10;
         ring += number(name->substr(index, 1));
         ring_digit = true;
         index++;
       }
       if (!ring_digit)
         parsing_error = true;
 
       if (name->substr(index, 1) != std::string(" "))
         parsing_error = true;
       index++;
 
       int chamber = 0;
       bool chamber_digit = false;
       while (!parsing_error && numeric(name->substr(index, 1))) {
         chamber *= 10;
         chamber += number(name->substr(index, 1));
         chamber_digit = true;
         index++;
       }
       if (!chamber_digit)
         parsing_error = true;
 
       if (!parsing_error) {
         bool no_such_chamber = false;
 
         int endcap = (station > 0 ? 1 : 2);
         station = abs(station);
         if (station < 1 || station > 4)
           no_such_chamber = true;
         if (station == 1 && (ring < 1 || ring > 4))
           no_such_chamber = true;
         if (station > 1 && (ring < 1 || ring > 2))
           no_such_chamber = true;
         if (station == 1 && (chamber < 1 || chamber > 36))
           no_such_chamber = true;
         if (station > 1 && ring == 1 && (chamber < 1 || chamber > 18))
           no_such_chamber = true;
         if (station > 1 && ring == 2 && (chamber < 1 || chamber > 36))
           no_such_chamber = true;
 
         if (no_such_chamber)
           throw cms::Exception("MuonAlignmentFromReference")
               << "reference chamber doesn't exist: " << (*name) << std::endl;
 
         CSCDetId id(endcap, station, ring, chamber);
         for (const auto& ali : all_CSC_chambers) {
           if (ali->geomDetId().rawId() == id.rawId()) {
             std::map<Alignable*, bool>::const_iterator trial = already_seen.find(ali);
             if (trial == already_seen.end()) {
               reference.push_back(ali);
               already_seen[ali] = true;
             }
           }
         }
       }  // if (!parsing_error)
     }    // endcap
 
     if (parsing_error)
       throw cms::Exception("MuonAlignmentFromReference")
           << "reference chamber name is malformed: " << (*name) << std::endl;
   }
 }

◆ processMuonResidualsFromTrack()

void MuonAlignmentFromReference::processMuonResidualsFromTrack ( MuonResidualsFromTrack & mrft )

Definition at line 529 of file MuonAlignmentFromReference.cc.

Referenced by run().

                                                                                            {
   // std::cout << "minTrackerHits: " << mrft.trackerNumHits() << std::endl;
   if (mrft.trackerNumHits() >= m_minTrackerHits) {
     m_counter_trackerhits++;
     // std::cout << "mrft.normalizedChi2(): " << mrft.normalizedChi2() << std::endl;
 
     if (mrft.normalizedChi2() < m_maxTrackerRedChi2) {
       m_counter_trackerchi2++;
       if (m_allowTIDTEC || !mrft.contains_TIDTEC()) {
         m_counter_trackertidtec++;
 
         std::vector<DetId> chamberIds = mrft.chamberIds();
 
         if ((int)chamberIds.size() >= m_minNCrossedChambers) {
           m_counter_minchambers++;
 
           char charge = (mrft.getTrack()->charge() > 0 ? 1 : -1);
 
           for (std::vector<DetId>::const_iterator chamberId = chamberIds.begin(); chamberId != chamberIds.end();
                ++chamberId) {
             if (chamberId->det() != DetId::Muon)
               continue;
             m_counter_totchambers++;
 
             // DT station 1,2,3
             if (m_doDT && chamberId->subdetId() == MuonSubdetId::DT && DTChamberId(chamberId->rawId()).station() != 4) {
               MuonChamberResidual* dt13 = mrft.chamberResidual(*chamberId, MuonChamberResidual::kDT13);
               MuonChamberResidual* dt2 = mrft.chamberResidual(*chamberId, MuonChamberResidual::kDT2);
 
               m_counter_station123++;
               if (dt13 != nullptr && dt2 != nullptr) {
                 m_counter_station123valid++;
                 if (dt13->numHits() >= m_minDT13Hits) {
                   m_counter_station123dt13hits++;
                   if (dt2->numHits() >= m_minDT2Hits) {
                     m_counter_station123dt2hits++;
                     std::map<Alignable*, MuonResidualsTwoBin*>::const_iterator fitter =
                         m_fitters.find(dt13->chamberAlignable());
                     if (fitter != m_fitters.end()) {
                       m_counter_station123aligning++;
                       if (fabs(dt2->resslope()) < m_maxResSlopeY && (dt2->chi2() / double(dt2->ndof())) < 2.0) {
                         m_counter_resslopey++;
                         double* residdata = new double[MuonResiduals6DOFFitter::kNData];
                         residdata[MuonResiduals6DOFFitter::kResidX] = dt13->residual();
                         residdata[MuonResiduals6DOFFitter::kResidY] = dt2->residual();
                         residdata[MuonResiduals6DOFFitter::kResSlopeX] = dt13->resslope();
                         residdata[MuonResiduals6DOFFitter::kResSlopeY] = dt2->resslope();
                         residdata[MuonResiduals6DOFFitter::kPositionX] = dt13->trackx();
                         residdata[MuonResiduals6DOFFitter::kPositionY] = dt13->tracky();
                         residdata[MuonResiduals6DOFFitter::kAngleX] = dt13->trackdxdz();
                         residdata[MuonResiduals6DOFFitter::kAngleY] = dt13->trackdydz();
                         residdata[MuonResiduals6DOFFitter::kRedChi2] =
                             (dt13->chi2() + dt2->chi2()) / double(dt13->ndof() + dt2->ndof());
                         residdata[MuonResiduals6DOFFitter::kPz] = mrft.getTrack()->pz();
                         residdata[MuonResiduals6DOFFitter::kPt] = mrft.getTrack()->pt();
                         residdata[MuonResiduals6DOFFitter::kCharge] = mrft.getTrack()->charge();
                         residdata[MuonResiduals6DOFFitter::kStation] = DTChamberId(chamberId->rawId()).station();
                         residdata[MuonResiduals6DOFFitter::kWheel] = DTChamberId(chamberId->rawId()).wheel();
                         residdata[MuonResiduals6DOFFitter::kSector] = DTChamberId(chamberId->rawId()).sector();
                         residdata[MuonResiduals6DOFFitter::kChambW] = dt13->ChambW();
                         residdata[MuonResiduals6DOFFitter::kChambl] = dt13->Chambl();
 
                         if (m_debug) {
                           std::cout << "processMuonResidualsFromTrack 6DOF dt13->residual()  " << dt13->residual()
                                     << std::endl;
                           std::cout << "                                   dt2->residual()   " << dt2->residual()
                                     << std::endl;
                           std::cout << "                                   dt13->resslope()  " << dt13->resslope()
                                     << std::endl;
                           std::cout << "                                   dt2->resslope()   " << dt2->resslope()
                                     << std::endl;
                           std::cout << "                                   dt13->trackx()    " << dt13->trackx()
                                     << std::endl;
                           std::cout << "                                   dt13->tracky()    " << dt13->tracky()
                                     << std::endl;
                           std::cout << "                                   dt13->trackdxdz() " << dt13->trackdxdz()
                                     << std::endl;
                           std::cout << "                                   dt13->trackdydz() " << dt13->trackdydz()
                                     << std::endl;
                         }
 
                         fitter->second->fill(charge, residdata);
                         // the MuonResidualsFitter will delete the array when it is destroyed
                       }
                     }
                   }
                 }
               }
             }
 
             // DT 4th station
             else if (m_doDT && chamberId->subdetId() == MuonSubdetId::DT &&
                      DTChamberId(chamberId->rawId()).station() == 4) {
               MuonChamberResidual* dt13 = mrft.chamberResidual(*chamberId, MuonChamberResidual::kDT13);
 
               m_counter_station4++;
               if (dt13 != nullptr) {
                 m_counter_station4valid++;
                 if (dt13->numHits() >= m_minDT13Hits) {
                   m_counter_station4hits++;
 
                   std::map<Alignable*, MuonResidualsTwoBin*>::const_iterator fitter =
                       m_fitters.find(dt13->chamberAlignable());
                   if (fitter != m_fitters.end()) {
                     m_counter_station4aligning++;
 
                     double* residdata = new double[MuonResiduals5DOFFitter::kNData];
                     residdata[MuonResiduals5DOFFitter::kResid] = dt13->residual();
                     residdata[MuonResiduals5DOFFitter::kResSlope] = dt13->resslope();
                     residdata[MuonResiduals5DOFFitter::kPositionX] = dt13->trackx();
                     residdata[MuonResiduals5DOFFitter::kPositionY] = dt13->tracky();
                     residdata[MuonResiduals5DOFFitter::kAngleX] = dt13->trackdxdz();
                     residdata[MuonResiduals5DOFFitter::kAngleY] = dt13->trackdydz();
                     residdata[MuonResiduals5DOFFitter::kRedChi2] = dt13->chi2() / double(dt13->ndof());
                     residdata[MuonResiduals5DOFFitter::kPz] = mrft.getTrack()->pz();
                     residdata[MuonResiduals5DOFFitter::kPt] = mrft.getTrack()->pt();
                     residdata[MuonResiduals5DOFFitter::kCharge] = mrft.getTrack()->charge();
                     residdata[MuonResiduals5DOFFitter::kStation] = DTChamberId(chamberId->rawId()).station();
                     residdata[MuonResiduals5DOFFitter::kWheel] = DTChamberId(chamberId->rawId()).wheel();
                     residdata[MuonResiduals5DOFFitter::kSector] = DTChamberId(chamberId->rawId()).sector();
                     residdata[MuonResiduals5DOFFitter::kChambW] = dt13->ChambW();
                     residdata[MuonResiduals5DOFFitter::kChambl] = dt13->Chambl();
 
                     if (m_debug) {
                       std::cout << "processMuonResidualsFromTrack 5DOF dt13->residual()  " << dt13->residual()
                                 << std::endl;
                       std::cout << "                                   dt13->resslope()  " << dt13->resslope()
                                 << std::endl;
                       std::cout << "                                   dt13->trackx()    " << dt13->trackx()
                                 << std::endl;
                       std::cout << "                                   dt13->tracky()    " << dt13->tracky()
                                 << std::endl;
                       std::cout << "                                   dt13->trackdxdz() " << dt13->trackdxdz()
                                 << std::endl;
                       std::cout << "                                   dt13->trackdydz() " << dt13->trackdydz()
                                 << std::endl;
                     }
 
                     fitter->second->fill(charge, residdata);
                     // the MuonResidualsFitter will delete the array when it is destroyed
                   }
                 }
               }
             }  // end DT 4th station
 
             // CSC
             else if (m_doCSC && chamberId->subdetId() == MuonSubdetId::CSC) {
               MuonChamberResidual* csc = mrft.chamberResidual(*chamberId, MuonChamberResidual::kCSC);
               m_counter_csc++;
               if (csc != nullptr) {
                 m_counter_cscvalid++;
                 if (csc->numHits() >= m_minCSCHits) {
                   m_counter_cschits++;
                   Alignable* ali = csc->chamberAlignable();
 
                   CSCDetId id(ali->geomDetId().rawId());
                   if (m_combineME11 && id.station() == 1 && id.ring() == 4)
                     ali = m_me11map[ali];
 
                   std::map<Alignable*, MuonResidualsTwoBin*>::const_iterator fitter = m_fitters.find(ali);
                   if (fitter != m_fitters.end()) {
                     m_counter_cscaligning++;
                     double* residdata = new double[MuonResiduals6DOFrphiFitter::kNData];
                     residdata[MuonResiduals6DOFrphiFitter::kResid] = csc->residual();
                     residdata[MuonResiduals6DOFrphiFitter::kResSlope] = csc->resslope();
                     residdata[MuonResiduals6DOFrphiFitter::kPositionX] = csc->trackx();
                     residdata[MuonResiduals6DOFrphiFitter::kPositionY] = csc->tracky();
                     residdata[MuonResiduals6DOFrphiFitter::kAngleX] = csc->trackdxdz();
                     residdata[MuonResiduals6DOFrphiFitter::kAngleY] = csc->trackdydz();
                     residdata[MuonResiduals6DOFrphiFitter::kRedChi2] = csc->chi2() / double(csc->ndof());
                     residdata[MuonResiduals6DOFrphiFitter::kPz] = mrft.getTrack()->pz();
                     residdata[MuonResiduals6DOFrphiFitter::kPt] = mrft.getTrack()->pt();
                     residdata[MuonResiduals6DOFrphiFitter::kCharge] = mrft.getTrack()->charge();
 
                     if (m_debug) {
                       std::cout << "processMuonResidualsFromTrack 6DOFrphi csc->residual()  " << csc->residual()
                                 << std::endl;
                       std::cout << "                                       csc->resslope()  " << csc->resslope()
                                 << std::endl;
                       std::cout << "                                       csc->trackx()    " << csc->trackx()
                                 << std::endl;
                       std::cout << "                                       csc->tracky()    " << csc->tracky()
                                 << std::endl;
                       std::cout << "                                       csc->trackdxdz() " << csc->trackdxdz()
                                 << std::endl;
                       std::cout << "                                       csc->trackdydz() " << csc->trackdydz()
                                 << std::endl;
                     }
 
                     fitter->second->fill(charge, residdata);
                     // the MuonResidualsFitter will delete the array when it is destroyed
                   }
                 }
               }
             }  // end CSC
 
             else if (m_doDT && m_doCSC)
               assert(false);
 
           }  // end loop over chamberIds
         }    // # crossed muon chambers ok
       }      // endcap tracker ok
     }        // chi2 ok
   }          // trackerNumHits ok
 }

◆ readTmpFiles()

void MuonAlignmentFromReference::readTmpFiles ( )

private

Definition at line 1607 of file MuonAlignmentFromReference.cc.

References Exception, geometryDiff::file, MillePedeFileConverter_cfg::fileName, mps_fire::i, m_fitterOrder, m_indexes, m_readTemporaryFiles, MuonResidualsTwoBin::read(), and findQualityFiles::size.

Referenced by terminate().

                                               {
   for (std::vector<std::string>::const_iterator fileName = m_readTemporaryFiles.begin();
        fileName != m_readTemporaryFiles.end();
        ++fileName) {
     FILE* file;
     int size;
     file = fopen((*fileName).c_str(), "r");
     if (file == nullptr)
       throw cms::Exception("MuonAlignmentFromReference")
           << "file \"" << *fileName << "\" can't be opened (doesn't exist?)" << std::endl;
 
     fread(&size, sizeof(int), 1, file);
     if (int(m_indexes.size()) != size)
       throw cms::Exception("MuonAlignmentFromReference")
           << "file \"" << *fileName << "\" has " << size << " fitters, but this job has " << m_indexes.size()
           << " fitters (probably corresponds to the wrong alignment job)" << std::endl;
 
     int i = 0;
     for (std::vector<unsigned int>::const_iterator index = m_indexes.begin(); index != m_indexes.end(); ++index, ++i) {
       MuonResidualsTwoBin* fitter = m_fitterOrder[*index];
       unsigned int index_toread;
       fread(&index_toread, sizeof(unsigned int), 1, file);
       if (*index != index_toread)
         throw cms::Exception("MuonAlignmentFromReference")
             << "file \"" << *fileName << "\" has index " << index_toread << " at position " << i
             << ", but this job is expecting " << *index << " (probably corresponds to the wrong alignment job)"
             << std::endl;
       fitter->read(file, i);
     }
 
     fclose(file);
   }
 }

◆ run()

void MuonAlignmentFromReference::run	(	const edm::EventSetup &	setup,
		const EventInfo &	eventInfo
	)

overridevirtual

Run the algorithm (must be implemented in derived class)

Implements AlignmentAlgorithmBase.

Definition at line 441 of file MuonAlignmentFromReference.cc.

References gather_cfg::cout, dumpTauVariables_cfi::eventInfo, edm::EventSetup::getData(), edm::EventSetup::getHandle(), edm::InputTag::label(), m_alignableNavigator, m_builderToken, m_counter_events, m_counter_trackdxy, m_counter_trackmomentum, m_counter_tracks, m_debug, m_DetIdToken, m_globTackingToken, m_MagFieldToken, m_maxDxy, m_maxTrackP, m_maxTrackPt, m_minTrackP, m_muonCollectionTag, m_propToken, HLT_2022v15_cff::magneticField, HLT_2022v15_cff::muonDetIdAssociator, processMuonResidualsFromTrack(), DiDispStaMuonMonitor_cfi::pt, and HLT_2022v15_cff::track.

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

                                                                                             {
   if (m_debug)
     std::cout << "****** EVENT START *******" << std::endl;
   m_counter_events++;
 
   const GlobalTrackingGeometry* globalGeometry = &iSetup.getData(m_globTackingToken);
   const MagneticField* magneticField = &iSetup.getData(m_MagFieldToken);
   const Propagator* prop = &iSetup.getData(m_propToken);
   const DetIdAssociator* muonDetIdAssociator = &iSetup.getData(m_DetIdToken);
   auto builder = iSetup.getHandle(m_builderToken);
 
   if (m_muonCollectionTag.label().empty())  // use trajectories
   {
     if (m_debug)
       std::cout << "JUST BEFORE LOOP OVER trajTrackPairs" << std::endl;
     // const ConstTrajTrackPairCollection &trajtracks = eventInfo.trajTrackPairs_; // trajTrackPairs_ now private
     const ConstTrajTrackPairCollection& trajtracks = eventInfo.trajTrackPairs();
 
     for (ConstTrajTrackPairCollection::const_iterator trajtrack = trajtracks.begin(); trajtrack != trajtracks.end();
          ++trajtrack) {
       m_counter_tracks++;
 
       const Trajectory* traj = (*trajtrack).first;
       const reco::Track* track = (*trajtrack).second;
 
       if (m_minTrackPt < track->pt() && track->pt() < m_maxTrackPt && m_minTrackP < track->p() &&
           track->p() < m_maxTrackP) {
         m_counter_trackmomentum++;
 
         if (fabs(track->dxy(eventInfo.beamSpot().position())) < m_maxDxy) {
           m_counter_trackdxy++;
           if (m_debug)
             std::cout << "JUST BEFORE muonResidualsFromTrack" << std::endl;
           MuonResidualsFromTrack muonResidualsFromTrack(builder,
                                                         magneticField,
                                                         globalGeometry,
                                                         muonDetIdAssociator,
                                                         prop,
                                                         traj,
                                                         track,
                                                         m_alignableNavigator,
                                                         1000.);
           if (m_debug)
             std::cout << "JUST AFTER muonResidualsFromTrack" << std::endl;
 
           if (m_debug)
             std::cout << "JUST BEFORE PROCESS" << std::endl;
           processMuonResidualsFromTrack(muonResidualsFromTrack);
           if (m_debug)
             std::cout << "JUST AFTER PROCESS" << std::endl;
         }
       }  // end if track p is within range
     }    // end if track pT is within range
     if (m_debug)
       std::cout << "JUST AFTER LOOP OVER trajTrackPairs" << std::endl;
 
   } else  // use muons
   {
     /*
            for (reco::MuonCollection::const_iterator muon = eventInfo.muonCollection_->begin();  muon != eventInfo.muonCollection_->end();  ++muon)
            {
            if ( !(muon->isTrackerMuon() && muon->innerTrack().isNonnull() ) ) continue;
 
            m_counter_tracks++;
 
            if (m_minTrackPt < muon->pt()  &&  muon->pt() < m_maxTrackPt && m_minTrackP < muon->p()  &&  muon->p() < m_maxTrackP)
            {
            m_counter_trackmomentum++;
 
            if (fabs(muon->innerTrack()->dxy(eventInfo.beamSpot_.position())) < m_maxDxy)
            {
            m_counter_trackdxy++;
 
         //std::cout<<"    *** will make MuonResidualsFromTrack ***"<<std::endl;
         MuonResidualsFromTrack muonResidualsFromTrack(globalGeometry, &(*muon), m_alignableNavigator, 100.);
         //std::cout<<"    *** have made MuonResidualsFromTrack ***"<<std::endl;
 
         //std::cout<<" trk eta="<<muon->eta()<<" ndof="<<muon->innerTrack()->ndof()<<" nchi2="<<muon->innerTrack()->normalizedChi2()
         //         <<" muresnchi2="<<muonResidualsFromTrack.normalizedChi2()<<" muresnhits="<<muonResidualsFromTrack.trackerNumHits()<<std::endl;
 
         processMuonResidualsFromTrack(muonResidualsFromTrack);
         } // end if track p is within range
         } // end if track pT is within range
         } // end loop over tracks
         */
   }
 }

◆ selectResidualsPeaks()

void MuonAlignmentFromReference::selectResidualsPeaks ( )

private

Definition at line 1687 of file MuonAlignmentFromReference.cc.

Referenced by terminate().

                                                       {
   // should not be called with negative peakNSigma
   assert(m_peakNSigma > 0.);
 
   for (std::vector<unsigned int>::const_iterator index = m_indexes.begin(); index != m_indexes.end(); ++index) {
     MuonResidualsTwoBin* fitter = m_fitterOrder[*index];
 
     int nvar = 2;
     int vars_index[10] = {0, 1};
     if (fitter->type() == MuonResidualsFitter::k5DOF) {
       if (fitter->useRes() == MuonResidualsFitter::k1111 || fitter->useRes() == MuonResidualsFitter::k1110 ||
           fitter->useRes() == MuonResidualsFitter::k1010) {
         nvar = 2;
         vars_index[0] = MuonResiduals5DOFFitter::kResid;
         vars_index[1] = MuonResiduals5DOFFitter::kResSlope;
       } else if (fitter->useRes() == MuonResidualsFitter::k1100) {
         nvar = 1;
         vars_index[0] = MuonResiduals5DOFFitter::kResid;
       } else if (fitter->useRes() == MuonResidualsFitter::k0010) {
         nvar = 1;
         vars_index[0] = MuonResiduals5DOFFitter::kResSlope;
       }
     } else if (fitter->type() == MuonResidualsFitter::k6DOF) {
       if (fitter->useRes() == MuonResidualsFitter::k1111) {
         nvar = 4;
         vars_index[0] = MuonResiduals6DOFFitter::kResidX;
         vars_index[1] = MuonResiduals6DOFFitter::kResidY;
         vars_index[2] = MuonResiduals6DOFFitter::kResSlopeX;
         vars_index[3] = MuonResiduals6DOFFitter::kResSlopeY;
       } else if (fitter->useRes() == MuonResidualsFitter::k1110) {
         nvar = 3;
         vars_index[0] = MuonResiduals6DOFFitter::kResidX;
         vars_index[1] = MuonResiduals6DOFFitter::kResidY;
         vars_index[2] = MuonResiduals6DOFFitter::kResSlopeX;
       } else if (fitter->useRes() == MuonResidualsFitter::k1010) {
         nvar = 2;
         vars_index[0] = MuonResiduals6DOFFitter::kResidX;
         vars_index[2] = MuonResiduals6DOFFitter::kResSlopeX;
       } else if (fitter->useRes() == MuonResidualsFitter::k1100) {
         nvar = 2;
         vars_index[0] = MuonResiduals6DOFFitter::kResidX;
         vars_index[1] = MuonResiduals6DOFFitter::kResidY;
       } else if (fitter->useRes() == MuonResidualsFitter::k0010) {
         nvar = 1;
         vars_index[0] = MuonResiduals6DOFFitter::kResSlopeX;
       }
     } else if (fitter->type() == MuonResidualsFitter::k6DOFrphi) {
       if (fitter->useRes() == MuonResidualsFitter::k1111 || fitter->useRes() == MuonResidualsFitter::k1110 ||
           fitter->useRes() == MuonResidualsFitter::k1010) {
         nvar = 2;
         vars_index[0] = MuonResiduals6DOFrphiFitter::kResid;
         vars_index[1] = MuonResiduals6DOFrphiFitter::kResSlope;
       } else if (fitter->useRes() == MuonResidualsFitter::k1100) {
         nvar = 1;
         vars_index[0] = MuonResiduals6DOFrphiFitter::kResid;
       } else if (fitter->useRes() == MuonResidualsFitter::k0010) {
         nvar = 1;
         vars_index[0] = MuonResiduals6DOFrphiFitter::kResSlope;
       }
     } else
       assert(false);
 
     if (m_debug)
       std::cout << "selecting in " << chamberPrettyNameFromId(*index) << std::endl;
 
     fitter->selectPeakResiduals(m_peakNSigma, nvar, vars_index);
   }
 }

◆ startNewLoop()

void MuonAlignmentFromReference::startNewLoop ( void )

inlineoverridevirtual

Called at start of loop, default implementation is dummy for non-iterative algorithms

Reimplemented from AlignmentAlgorithmBase.

Definition at line 85 of file MuonAlignmentFromReference.cc.

85 {};

◆ terminate()

void MuonAlignmentFromReference::terminate ( const edm::EventSetup & iSetup )

overridevirtual

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

Implements AlignmentAlgorithmBase.

Definition at line 735 of file MuonAlignmentFromReference.cc.

References correctBField(), gather_cfg::cout, eraseNotSelectedResiduals(), fiducialCuts(), fillNtuple(), fitAndAlign(), m_BFieldCorrection, m_counter_csc, m_counter_cscaligning, m_counter_cschits, m_counter_cscvalid, m_counter_events, m_counter_minchambers, m_counter_resslopey, m_counter_station123, m_counter_station123aligning, m_counter_station123dt13hits, m_counter_station123dt2hits, m_counter_station123valid, m_counter_station4, m_counter_station4aligning, m_counter_station4hits, m_counter_station4valid, m_counter_totchambers, m_counter_trackdxy, m_counter_trackerchi2, m_counter_trackerhits, m_counter_trackertidtec, m_counter_trackmomentum, m_counter_tracks, m_createNtuple, m_debug, m_doAlignment, m_doCSC, m_peakNSigma, m_readTemporaryFiles, m_writeTemporaryFile, readTmpFiles(), selectResidualsPeaks(), AlCaHLTBitMon_QueryRunRegistry::string, and writeTmpFiles().

                                                                       {
   bool m_debug = false;
 
   // one-time print-out
   std::cout << "Counters:" << std::endl
             << "COUNT{ events: " << m_counter_events << " }" << std::endl
             << "COUNT{  tracks: " << m_counter_tracks << " }" << std::endl
             << "COUNT{   trackppt: " << m_counter_trackmomentum << " }" << std::endl
             << "COUNT{    trackdxy: " << m_counter_trackdxy << " }" << std::endl
             << "COUNT{     trackerhits: " << m_counter_trackerhits << " }" << std::endl
             << "COUNT{      trackerchi2: " << m_counter_trackerchi2 << " }" << std::endl
             << "COUNT{       trackertidtec: " << m_counter_trackertidtec << " }" << std::endl
             << "COUNT{        minnchambers: " << m_counter_minchambers << " }" << std::endl
             << "COUNT{          totchambers: " << m_counter_totchambers << " }" << std::endl
             << "COUNT{            station123:             " << m_counter_station123 << " }" << std::endl
             << "COUNT{             station123valid:       " << m_counter_station123valid << " }" << std::endl
             << "COUNT{              station123dt13hits:   " << m_counter_station123dt13hits << " }" << std::endl
             << "COUNT{               station123dt2hits:   " << m_counter_station123dt2hits << " }" << std::endl
             << "COUNT{                station123aligning: " << m_counter_station123aligning << " }" << std::endl
             << "COUNT{                 resslopey: " << m_counter_resslopey << " }" << std::endl
             << "COUNT{            station4:            " << m_counter_station4 << " }" << std::endl
             << "COUNT{             station4valid:      " << m_counter_station4valid << " }" << std::endl
             << "COUNT{              station4hits:      " << m_counter_station4hits << " }" << std::endl
             << "COUNT{               station4aligning: " << m_counter_station4aligning << " }" << std::endl
             << "COUNT{            csc:            " << m_counter_csc << " }" << std::endl
             << "COUNT{             cscvalid:      " << m_counter_cscvalid << " }" << std::endl
             << "COUNT{              cschits:      " << m_counter_cschits << " }" << std::endl
             << "COUNT{               cscaligning: " << m_counter_cscaligning << " }" << std::endl
             << "That's all!" << std::endl;
 
   TStopwatch stop_watch;
 
   // collect temporary files
   if (!m_readTemporaryFiles.empty()) {
     stop_watch.Start();
     readTmpFiles();
     if (m_debug)
       std::cout << "readTmpFiles took " << stop_watch.CpuTime() << " sec" << std::endl;
     stop_watch.Stop();
   }
 
   // select residuals peaks and discard tails if peakNSigma>0 (only while doing alignment)
   if (m_peakNSigma > 0. && m_doAlignment) {
     stop_watch.Start();
     selectResidualsPeaks();
     if (m_debug)
       std::cout << "selectResidualsPeaks took " << stop_watch.CpuTime() << " sec" << std::endl;
     stop_watch.Stop();
   }
 
   if (m_BFieldCorrection > 0 && m_doAlignment) {
     stop_watch.Start();
     correctBField();
     if (m_debug)
       std::cout << "correctBField took " << stop_watch.CpuTime() << " sec" << std::endl;
     stop_watch.Stop();
   }
 
   if (m_doAlignment && !m_doCSC)  // for now apply fiducial cuts to DT only
   {
     stop_watch.Start();
     fiducialCuts();
     if (m_debug)
       std::cout << "fiducialCuts took " << stop_watch.CpuTime() << " sec" << std::endl;
     stop_watch.Stop();
   }
 
   // optionally, create an nutuple for easy debugging
   if (m_createNtuple) {
     stop_watch.Start();
     fillNtuple();
     if (m_debug)
       std::cout << "fillNtuple took " << stop_watch.CpuTime() << " sec" << std::endl;
     stop_watch.Stop();
   }
 
   if (m_doAlignment) {
     stop_watch.Start();
     eraseNotSelectedResiduals();
     if (m_debug)
       std::cout << "eraseNotSelectedResiduals took " << stop_watch.CpuTime() << " sec" << std::endl;
     stop_watch.Stop();
   }
 
   // fit and align (time-consuming, so the user can turn it off if in a residuals-gathering job)
   if (m_doAlignment) {
     stop_watch.Start();
     fitAndAlign();
     if (m_debug)
       std::cout << "fitAndAlign took " << stop_watch.CpuTime() << " sec" << std::endl;
     stop_watch.Stop();
   }
 
   // write out the pseudontuples for a later job to collect
   if (m_writeTemporaryFile != std::string(""))
     writeTmpFiles();
   if (m_debug)
     std::cout << "end: MuonAlignmentFromReference::terminate()" << std::endl;
 }

◆ writeTmpFiles()

void MuonAlignmentFromReference::writeTmpFiles ( )

private

Definition at line 1641 of file MuonAlignmentFromReference.cc.

References geometryDiff::file, mps_fire::i, m_fitterOrder, m_indexes, m_writeTemporaryFile, findQualityFiles::size, and MuonResidualsTwoBin::write().

Referenced by terminate().

                                                {
   FILE* file;
   file = fopen(m_writeTemporaryFile.c_str(), "w");
   int size = m_indexes.size();
   fwrite(&size, sizeof(int), 1, file);
 
   int i = 0;
   for (std::vector<unsigned int>::const_iterator index = m_indexes.begin(); index != m_indexes.end(); ++index, ++i) {
     MuonResidualsTwoBin* fitter = m_fitterOrder[*index];
     unsigned int index_towrite = *index;
     fwrite(&index_towrite, sizeof(unsigned int), 1, file);
     fitter->write(file, i);
   }
 
   fclose(file);
 }

Member Data Documentation

◆ m_alignableNavigator

AlignableNavigator* MuonAlignmentFromReference::m_alignableNavigator

private

Definition at line 155 of file MuonAlignmentFromReference.cc.

Referenced by initialize(), run(), and ~MuonAlignmentFromReference().

◆ m_alignables

align::Alignables MuonAlignmentFromReference::m_alignables

private

Definition at line 157 of file MuonAlignmentFromReference.cc.

Referenced by fitAndAlign(), and initialize().

◆ m_alignmentParameterStore

AlignmentParameterStore* MuonAlignmentFromReference::m_alignmentParameterStore

private

Definition at line 156 of file MuonAlignmentFromReference.cc.

Referenced by fitAndAlign(), and initialize().

◆ m_allowTIDTEC

bool MuonAlignmentFromReference::m_allowTIDTEC

private

Definition at line 131 of file MuonAlignmentFromReference.cc.

Referenced by processMuonResidualsFromTrack().

◆ m_BFieldCorrection

int MuonAlignmentFromReference::m_BFieldCorrection

private

Definition at line 149 of file MuonAlignmentFromReference.cc.

Referenced by terminate().

◆ m_builderToken

const MuonResidualsFromTrack::BuilderToken MuonAlignmentFromReference::m_builderToken

private

Definition at line 119 of file MuonAlignmentFromReference.cc.

Referenced by run().

◆ m_combineME11

bool MuonAlignmentFromReference::m_combineME11

private

Definition at line 143 of file MuonAlignmentFromReference.cc.

Referenced by fitAndAlign(), initialize(), and processMuonResidualsFromTrack().

◆ m_counter_csc

long MuonAlignmentFromReference::m_counter_csc

private

Definition at line 182 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), processMuonResidualsFromTrack(), and terminate().

◆ m_counter_cscaligning

long MuonAlignmentFromReference::m_counter_cscaligning

private

Definition at line 185 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), processMuonResidualsFromTrack(), and terminate().

◆ m_counter_cschits

long MuonAlignmentFromReference::m_counter_cschits

private

Definition at line 184 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), processMuonResidualsFromTrack(), and terminate().

◆ m_counter_cscvalid

long MuonAlignmentFromReference::m_counter_cscvalid

private

Definition at line 183 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), processMuonResidualsFromTrack(), and terminate().

◆ m_counter_events

long MuonAlignmentFromReference::m_counter_events

private

Definition at line 164 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), run(), and terminate().

◆ m_counter_minchambers

long MuonAlignmentFromReference::m_counter_minchambers

private

Definition at line 171 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), processMuonResidualsFromTrack(), and terminate().

◆ m_counter_resslopey

long MuonAlignmentFromReference::m_counter_resslopey

private

Definition at line 186 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), processMuonResidualsFromTrack(), and terminate().

◆ m_counter_station123

long MuonAlignmentFromReference::m_counter_station123

private

Definition at line 173 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), processMuonResidualsFromTrack(), and terminate().

◆ m_counter_station123aligning

long MuonAlignmentFromReference::m_counter_station123aligning

private

Definition at line 177 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), processMuonResidualsFromTrack(), and terminate().

◆ m_counter_station123dt13hits

long MuonAlignmentFromReference::m_counter_station123dt13hits

private

Definition at line 175 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), processMuonResidualsFromTrack(), and terminate().

◆ m_counter_station123dt2hits

long MuonAlignmentFromReference::m_counter_station123dt2hits

private

Definition at line 176 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), processMuonResidualsFromTrack(), and terminate().

◆ m_counter_station123valid

long MuonAlignmentFromReference::m_counter_station123valid

private

Definition at line 174 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), processMuonResidualsFromTrack(), and terminate().

◆ m_counter_station4

long MuonAlignmentFromReference::m_counter_station4

private

Definition at line 178 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), processMuonResidualsFromTrack(), and terminate().

◆ m_counter_station4aligning

long MuonAlignmentFromReference::m_counter_station4aligning

private

Definition at line 181 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), processMuonResidualsFromTrack(), and terminate().

◆ m_counter_station4hits

long MuonAlignmentFromReference::m_counter_station4hits

private

Definition at line 180 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), processMuonResidualsFromTrack(), and terminate().

◆ m_counter_station4valid

long MuonAlignmentFromReference::m_counter_station4valid

private

Definition at line 179 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), processMuonResidualsFromTrack(), and terminate().

◆ m_counter_totchambers

long MuonAlignmentFromReference::m_counter_totchambers

private

Definition at line 172 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), processMuonResidualsFromTrack(), and terminate().

◆ m_counter_trackdxy

long MuonAlignmentFromReference::m_counter_trackdxy

private

Definition at line 167 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), run(), and terminate().

◆ m_counter_trackerchi2

long MuonAlignmentFromReference::m_counter_trackerchi2

private

Definition at line 169 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), processMuonResidualsFromTrack(), and terminate().

◆ m_counter_trackerhits

long MuonAlignmentFromReference::m_counter_trackerhits

private

Definition at line 168 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), processMuonResidualsFromTrack(), and terminate().

◆ m_counter_trackertidtec

long MuonAlignmentFromReference::m_counter_trackertidtec

private

Definition at line 170 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), processMuonResidualsFromTrack(), and terminate().

◆ m_counter_trackmomentum

long MuonAlignmentFromReference::m_counter_trackmomentum

private

Definition at line 166 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), run(), and terminate().

◆ m_counter_tracks

long MuonAlignmentFromReference::m_counter_tracks

private

Definition at line 165 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), run(), and terminate().

◆ m_createNtuple

bool MuonAlignmentFromReference::m_createNtuple

private

Definition at line 147 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), and terminate().

◆ m_cscGeometryToken

const edm::ESGetToken<CSCGeometry, MuonGeometryRecord> MuonAlignmentFromReference::m_cscGeometryToken

private

Definition at line 114 of file MuonAlignmentFromReference.cc.

Referenced by initialize().

◆ m_debug

bool MuonAlignmentFromReference::m_debug

private

Definition at line 193 of file MuonAlignmentFromReference.cc.

Referenced by correctBField(), eraseNotSelectedResiduals(), fiducialCuts(), fitAndAlign(), MuonAlignmentFromReference(), processMuonResidualsFromTrack(), run(), selectResidualsPeaks(), and terminate().

◆ m_DetIdToken

const edm::ESGetToken<DetIdAssociator, DetIdAssociatorRecord> MuonAlignmentFromReference::m_DetIdToken

private

Definition at line 118 of file MuonAlignmentFromReference.cc.

Referenced by run().

◆ m_doAlignment

bool MuonAlignmentFromReference::m_doAlignment

private

Definition at line 138 of file MuonAlignmentFromReference.cc.

Referenced by MuonAlignmentFromReference(), and terminate().

◆ m_doCSC

bool MuonAlignmentFromReference::m_doCSC

private

Definition at line 151 of file MuonAlignmentFromReference.cc.

Referenced by processMuonResidualsFromTrack(), and terminate().

◆ m_doDT

bool MuonAlignmentFromReference::m_doDT

private

Definition at line 150 of file MuonAlignmentFromReference.cc.

Referenced by processMuonResidualsFromTrack().

◆ m_fitterOrder

std::map<unsigned int, MuonResidualsTwoBin*> MuonAlignmentFromReference::m_fitterOrder

private

Definition at line 161 of file MuonAlignmentFromReference.cc.

Referenced by correctBField(), eraseNotSelectedResiduals(), fiducialCuts(), fillNtuple(), initialize(), readTmpFiles(), selectResidualsPeaks(), and writeTmpFiles().

◆ m_fitters

std::map<Alignable*, MuonResidualsTwoBin*> MuonAlignmentFromReference::m_fitters

private

Definition at line 159 of file MuonAlignmentFromReference.cc.

Referenced by fitAndAlign(), initialize(), and processMuonResidualsFromTrack().

◆ m_globTackingToken

const edm::ESGetToken<GlobalTrackingGeometry, GlobalTrackingGeometryRecord> MuonAlignmentFromReference::m_globTackingToken

private

Definition at line 115 of file MuonAlignmentFromReference.cc.

Referenced by run().

◆ m_indexes

std::vector<unsigned int> MuonAlignmentFromReference::m_indexes

private

Definition at line 160 of file MuonAlignmentFromReference.cc.

Referenced by correctBField(), eraseNotSelectedResiduals(), fiducialCuts(), fillNtuple(), initialize(), readTmpFiles(), selectResidualsPeaks(), and writeTmpFiles().

◆ m_MagFieldToken

const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> MuonAlignmentFromReference::m_MagFieldToken

private

Definition at line 116 of file MuonAlignmentFromReference.cc.

Referenced by run().

◆ m_maxDxy

double MuonAlignmentFromReference::m_maxDxy

private

Definition at line 128 of file MuonAlignmentFromReference.cc.

Referenced by run().

◆ m_maxResSlopeY

double MuonAlignmentFromReference::m_maxResSlopeY

private

Definition at line 146 of file MuonAlignmentFromReference.cc.

Referenced by processMuonResidualsFromTrack().

◆ m_maxTrackerRedChi2

double MuonAlignmentFromReference::m_maxTrackerRedChi2

private

Definition at line 130 of file MuonAlignmentFromReference.cc.

Referenced by processMuonResidualsFromTrack().

◆ m_maxTrackP

double MuonAlignmentFromReference::m_maxTrackP

private

Definition at line 127 of file MuonAlignmentFromReference.cc.

Referenced by run().

◆ m_maxTrackPt

double MuonAlignmentFromReference::m_maxTrackPt

private

Definition at line 125 of file MuonAlignmentFromReference.cc.

Referenced by run().

◆ m_me11map

std::map<Alignable*, Alignable*> MuonAlignmentFromReference::m_me11map

private

Definition at line 158 of file MuonAlignmentFromReference.cc.

Referenced by fitAndAlign(), initialize(), and processMuonResidualsFromTrack().

◆ m_minAlignmentHits

int MuonAlignmentFromReference::m_minAlignmentHits

private

Definition at line 141 of file MuonAlignmentFromReference.cc.

Referenced by fitAndAlign(), and initialize().

◆ m_minCSCHits

int MuonAlignmentFromReference::m_minCSCHits

private

Definition at line 135 of file MuonAlignmentFromReference.cc.

Referenced by processMuonResidualsFromTrack().

◆ m_minDT13Hits

int MuonAlignmentFromReference::m_minDT13Hits

private

Definition at line 133 of file MuonAlignmentFromReference.cc.

Referenced by processMuonResidualsFromTrack().

◆ m_minDT2Hits

int MuonAlignmentFromReference::m_minDT2Hits

private

Definition at line 134 of file MuonAlignmentFromReference.cc.

Referenced by processMuonResidualsFromTrack().

◆ m_minNCrossedChambers

int MuonAlignmentFromReference::m_minNCrossedChambers

private

Definition at line 132 of file MuonAlignmentFromReference.cc.

Referenced by processMuonResidualsFromTrack().

◆ m_minTrackerHits

int MuonAlignmentFromReference::m_minTrackerHits

private

Definition at line 129 of file MuonAlignmentFromReference.cc.

Referenced by processMuonResidualsFromTrack().

◆ m_minTrackP

double MuonAlignmentFromReference::m_minTrackP

private

Definition at line 126 of file MuonAlignmentFromReference.cc.

Referenced by run().

◆ m_minTrackPt

double MuonAlignmentFromReference::m_minTrackPt

private

Definition at line 124 of file MuonAlignmentFromReference.cc.

◆ m_muonCollectionTag

edm::InputTag MuonAlignmentFromReference::m_muonCollectionTag

private

Definition at line 122 of file MuonAlignmentFromReference.cc.

Referenced by run().

◆ m_peakNSigma

double MuonAlignmentFromReference::m_peakNSigma

private

Definition at line 148 of file MuonAlignmentFromReference.cc.

Referenced by selectResidualsPeaks(), and terminate().

◆ m_propToken

const edm::ESGetToken<Propagator, TrackingComponentsRecord> MuonAlignmentFromReference::m_propToken

private

Definition at line 117 of file MuonAlignmentFromReference.cc.

Referenced by run().

◆ m_readTemporaryFiles

std::vector<std::string> MuonAlignmentFromReference::m_readTemporaryFiles

private

Definition at line 137 of file MuonAlignmentFromReference.cc.

Referenced by readTmpFiles(), and terminate().

◆ m_reference

std::vector<std::string> MuonAlignmentFromReference::m_reference

private

Definition at line 123 of file MuonAlignmentFromReference.cc.

Referenced by initialize(), and parseReference().

◆ m_reportFileName

std::string MuonAlignmentFromReference::m_reportFileName

private

Definition at line 145 of file MuonAlignmentFromReference.cc.

Referenced by fitAndAlign().

◆ m_residualsModel

std::string MuonAlignmentFromReference::m_residualsModel

private

Definition at line 140 of file MuonAlignmentFromReference.cc.

Referenced by initialize().

◆ m_strategy

int MuonAlignmentFromReference::m_strategy

private

Definition at line 139 of file MuonAlignmentFromReference.cc.

Referenced by initialize().

◆ m_tree_row

MuonResidualsFitter::MuonAlignmentTreeRow MuonAlignmentFromReference::m_tree_row

private

Definition at line 191 of file MuonAlignmentFromReference.cc.

Referenced by bookNtuple(), and fillNtuple().

◆ m_ttree

TTree* MuonAlignmentFromReference::m_ttree

private

Definition at line 190 of file MuonAlignmentFromReference.cc.

Referenced by bookNtuple(), fillNtuple(), and MuonAlignmentFromReference().

◆ m_twoBin

bool MuonAlignmentFromReference::m_twoBin

private

Definition at line 142 of file MuonAlignmentFromReference.cc.

Referenced by initialize().

◆ m_useResiduals

std::string MuonAlignmentFromReference::m_useResiduals

private

Definition at line 152 of file MuonAlignmentFromReference.cc.

Referenced by initialize().

◆ m_weightAlignment

bool MuonAlignmentFromReference::m_weightAlignment

private

Definition at line 144 of file MuonAlignmentFromReference.cc.

Referenced by initialize().

◆ m_writeTemporaryFile

std::string MuonAlignmentFromReference::m_writeTemporaryFile

private

Definition at line 136 of file MuonAlignmentFromReference.cc.

Referenced by terminate(), and writeTmpFiles().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ MuonAlignmentFromReference()

◆ ~MuonAlignmentFromReference()

Member Function Documentation

◆ bookNtuple()

◆ chamberPrettyNameFromId()

◆ correctBField()

◆ eraseNotSelectedResiduals()

◆ fiducialCuts()

◆ fillNtuple()

◆ fitAndAlign()

◆ initialize()

◆ number()

◆ numeric()

◆ parseReference()

◆ processMuonResidualsFromTrack()

◆ readTmpFiles()

◆ run()

◆ selectResidualsPeaks()

◆ startNewLoop()

◆ terminate()

◆ writeTmpFiles()

Member Data Documentation

◆ m_alignableNavigator

◆ m_alignables

◆ m_alignmentParameterStore

◆ m_allowTIDTEC

◆ m_BFieldCorrection

◆ m_builderToken

◆ m_combineME11

◆ m_counter_csc

◆ m_counter_cscaligning

◆ m_counter_cschits

◆ m_counter_cscvalid

◆ m_counter_events

◆ m_counter_minchambers

◆ m_counter_resslopey

◆ m_counter_station123

◆ m_counter_station123aligning

◆ m_counter_station123dt13hits

◆ m_counter_station123dt2hits

◆ m_counter_station123valid

◆ m_counter_station4

◆ m_counter_station4aligning

◆ m_counter_station4hits

◆ m_counter_station4valid

◆ m_counter_totchambers

◆ m_counter_trackdxy

◆ m_counter_trackerchi2

◆ m_counter_trackerhits

◆ m_counter_trackertidtec

◆ m_counter_trackmomentum

◆ m_counter_tracks

◆ m_createNtuple

◆ m_cscGeometryToken

◆ m_debug

◆ m_DetIdToken

◆ m_doAlignment

◆ m_doCSC

◆ m_doDT

◆ m_fitterOrder

◆ m_fitters

◆ m_globTackingToken

◆ m_indexes

◆ m_MagFieldToken

◆ m_maxDxy

◆ m_maxResSlopeY

◆ m_maxTrackerRedChi2

◆ m_maxTrackP

◆ m_maxTrackPt

◆ m_me11map

◆ m_minAlignmentHits

◆ m_minCSCHits

◆ m_minDT13Hits

◆ m_minDT2Hits