Inheritance diagram for ApeSettingAlgorithm:

Public Member Functions
	ApeSettingAlgorithm (const edm::ParameterSet &cfg, const edm::ConsumesCollector &iC)
	Constructor. More...

void	initialize (const edm::EventSetup &setup, AlignableTracker tracker, AlignableMuon muon, AlignableExtras extras, AlignmentParameterStore store) override
	Call at beginning of job. More...

void	run (const edm::EventSetup &setup, const EventInfo &eventInfo) override
	Run the algorithm. More...

void	terminate (const edm::EventSetup &iSetup) override
	Call at end of job. More...

	~ApeSettingAlgorithm () override
	Destructor. More...

Public Member Functions inherited from AlignmentAlgorithmBase
virtual bool	addCalibrations (const Calibrations &)

bool	addCalibrations (const CalibrationsOwner &cals)

	AlignmentAlgorithmBase (const edm::ParameterSet &, const edm::ConsumesCollector &)
	Constructor. More...

virtual void	beginLuminosityBlock (const edm::EventSetup &setup)
	called at begin of luminosity block (no lumi block info passed yet) More...

virtual void	beginRun (const edm::Run &, const edm::EventSetup &, bool changed)
	called at begin of run More...

virtual void	endLuminosityBlock (const edm::EventSetup &setup)
	called at end of luminosity block (no lumi block info passed yet) More...

virtual void	endRun (const EndRunInfo &runInfo, const edm::EventSetup &setup)
	called at end of run - order of arguments like in EDProducer etc. More...

virtual bool	processesEvents ()
	Returns whether algorithm proccesses events in current configuration. More...

virtual bool	setParametersForRunRange (const RunRange &rr)

virtual void	startNewLoop ()

virtual bool	storeAlignments ()
	Returns whether algorithm produced results to be stored. More...

virtual bool	supportsCalibrations ()

virtual void	terminate ()
	Called at end of job (must be implemented in derived class) More...

virtual	~AlignmentAlgorithmBase ()
	Destructor. More...

Private Attributes
bool	readApeFromAscii_

bool	readFullLocalMatrix_

bool	readLocalNotGlobal_

bool	saveApeToAscii_

bool	saveComposites_

bool	saveLocalNotGlobal_

bool	setComposites_

AlignableNavigator *	theAlignableNavigator

edm::ParameterSet	theConfig

AlignableTracker *	theTracker

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

Definition at line 61 of file ApeSettingAlgorithm.cc.

Constructor & Destructor Documentation

◆ ApeSettingAlgorithm()

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

Constructor.

Definition at line 98 of file ApeSettingAlgorithm.cc.

     : AlignmentAlgorithmBase(cfg, iC), theConfig(cfg), theAlignableNavigator(nullptr) {
   edm::LogInfo("Alignment") << "@SUB=ApeSettingAlgorithm"
                             << "Start.";
   saveApeToAscii_ = theConfig.getUntrackedParameter<bool>("saveApeToASCII");
   saveComposites_ = theConfig.getUntrackedParameter<bool>("saveComposites");
   saveLocalNotGlobal_ = theConfig.getUntrackedParameter<bool>("saveLocalNotGlobal");
   readApeFromAscii_ = theConfig.getParameter<bool>("readApeFromASCII");
   readLocalNotGlobal_ = theConfig.getParameter<bool>("readLocalNotGlobal");
   readFullLocalMatrix_ = theConfig.getParameter<bool>("readFullLocalMatrix");
   setComposites_ = theConfig.getParameter<bool>("setComposites");
 }

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), readApeFromAscii_, readFullLocalMatrix_, readLocalNotGlobal_, saveApeToAscii_, saveComposites_, saveLocalNotGlobal_, setComposites_, and theConfig.

◆ ~ApeSettingAlgorithm()

ApeSettingAlgorithm::~ApeSettingAlgorithm ( )

override

Destructor.

Definition at line 113 of file ApeSettingAlgorithm.cc.

113 { delete theAlignableNavigator; }

References theAlignableNavigator.

Member Function Documentation

◆ initialize()

void ApeSettingAlgorithm::initialize	(	const edm::EventSetup &	setup,
		AlignableTracker *	tracker,
		AlignableMuon *	muon,
		AlignableExtras *	extras,
		AlignmentParameterStore *	store
	)

overridevirtual

Call at beginning of job.

Implements AlignmentAlgorithmBase.

Definition at line 117 of file ApeSettingAlgorithm.cc.

                                                                      {
   theAlignableNavigator = new AlignableNavigator(tracker, muon);
   theTracker = tracker;
  
   if (readApeFromAscii_) {
     std::ifstream apeReadFile(
         theConfig.getParameter<edm::FileInPath>("apeASCIIReadFile").fullPath().c_str());  //requires <fstream>
     if (!apeReadFile.good()) {
       edm::LogInfo("Alignment") << "@SUB=initialize"
                                 << "Problem opening APE file: skipping"
                                 << theConfig.getParameter<edm::FileInPath>("apeASCIIReadFile").fullPath();
       return;
     }
     std::set<int> apeList;  //To avoid duplicates
     while (!apeReadFile.eof()) {
       int apeId = 0;
       double x11, x21, x22, x31, x32, x33, ignore;
       if (!readLocalNotGlobal_ || readFullLocalMatrix_) {
         apeReadFile >> apeId >> x11 >> x21 >> x22 >> x31 >> x32 >> x33 >> ignore >> ignore >> ignore >> ignore >>
             ignore >> ignore >> ignore >> ignore >> ignore >> ignore >> ignore >> ignore >> ignore >> ignore >>
             ignore >> std::ws;
       } else {
         apeReadFile >> apeId >> x11 >> x22 >> x33 >> std::ws;
       }
       //idr What sanity checks do we need to put here?
       if (apeId != 0)  //read appears valid?
       {
         if (apeList.find(apeId) == apeList.end())  //Not previously done
         {
           DetId id(apeId);
           AlignableDetOrUnitPtr alidet(theAlignableNavigator->alignableFromDetId(id));  //NULL if none
           if (alidet) {
             if ((alidet->components().empty()) || setComposites_)  //the problem with glued dets...
             {
               GlobalErrorExtended globErr;
               if (readLocalNotGlobal_) {
                 AlgebraicSymMatrix33 as;
                 if (readFullLocalMatrix_) {
                   as[0][0] = x11;
                   as[1][0] = x21;
                   as[1][1] = x22;
                   as[2][0] = x31;
                   as[2][1] = x32;
                   as[2][2] = x33;
                 } else {
                   as[0][0] = x11 * x11;
                   as[1][1] = x22 * x22;
                   as[2][2] = x33 * x33;
                 }  //local cov.
                 align::RotationType rt = alidet->globalRotation();
                 AlgebraicMatrix33 am;
                 am[0][0] = rt.xx();
                 am[0][1] = rt.xy();
                 am[0][2] = rt.xz();
                 am[1][0] = rt.yx();
                 am[1][1] = rt.yy();
                 am[1][2] = rt.yz();
                 am[2][0] = rt.zx();
                 am[2][1] = rt.zy();
                 am[2][2] = rt.zz();
                 globErr = GlobalErrorExtended(ROOT::Math::SimilarityT(am, as));
               } else {
                 if (readFullLocalMatrix_)
                   globErr =
                       GlobalErrorExtended(x11, x21, x31, 0, 0, 0, x22, x32, 0, 0, 0, x33, 0, 0, 0, 0, 0, 0, 0, 0, 0);
                 else {
                   globErr = GlobalErrorExtended(
                       x11 * x11, 0, 0, 0, 0, 0, x22 * x22, 0, 0, 0, 0, x33 * x33, 0, 0, 0, 0, 0, 0, 0, 0, 0);
                 }
               }
               alidet->setAlignmentPositionError(globErr, false);  // do not propagate down!
               apeList.insert(apeId);                              //Flag it's been set
             } else {
               edm::LogInfo("Alignment") << "@SUB=initialize"
                                         << "Not Setting APE for Composite DetId " << apeId;
             }
           }
         } else {
           edm::LogInfo("Alignment") << "@SUB=initialize"
                                     << "Skipping duplicate APE for DetId " << apeId;
         }
       }
     }
     apeReadFile.close();
     edm::LogInfo("Alignment") << "@SUB=initialize"
                               << "Set " << apeList.size() << " APE values.";
   }
 }

References AlignableNavigator::alignableFromDetId(), contentValuesFiles::fullPath, edm::FileInPath::fullPath(), edm::ParameterSet::getParameter(), triggerObjects_cff::id, SequenceTypes::ignore(), readApeFromAscii_, readFullLocalMatrix_, readLocalNotGlobal_, hcal_runs::rt, setComposites_, theAlignableNavigator, theConfig, theTracker, PbPb_ZMuSkimMuonDPG_cff::tracker, and gpuVertexFinder::ws.

◆ run()

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

overridevirtual

Run the algorithm.

Implements AlignmentAlgorithmBase.

Definition at line 255 of file ApeSettingAlgorithm.cc.

                                                                                     {
   // nothing to do here?
 }

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

◆ terminate()

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

overridevirtual

Call at end of job.

Implements AlignmentAlgorithmBase.

Definition at line 212 of file ApeSettingAlgorithm.cc.

                                                                {
   if (saveApeToAscii_) {
     AlignmentErrorsExtended *aliErr = theTracker->alignmentErrors();
     int theSize = aliErr->m_alignError.size();
     std::ofstream apeSaveFile(
         theConfig.getUntrackedParameter<std::string>("apeASCIISaveFile").c_str());  //requires <fstream>
     for (int i = 0; i < theSize; ++i) {
       int id = aliErr->m_alignError[i].rawId();
       AlignableDetOrUnitPtr alidet(theAlignableNavigator->alignableFromDetId(DetId(id)));  //NULL if none
       if (alidet && ((alidet->components().empty()) || saveComposites_)) {
         apeSaveFile << id;
         CLHEP::HepSymMatrix sm = aliErr->m_alignError[i].matrix();
         if (saveLocalNotGlobal_) {
           align::RotationType rt = alidet->globalRotation();
           AlgebraicMatrix am(3, 3);
           am[0][0] = rt.xx();
           am[0][1] = rt.xy();
           am[0][2] = rt.xz();
           am[1][0] = rt.yx();
           am[1][1] = rt.yy();
           am[1][2] = rt.yz();
           am[2][0] = rt.zx();
           am[2][1] = rt.zy();
           am[2][2] = rt.zz();
           sm = sm.similarity(am);  //symmetric matrix
         }                          //transform to local
         for (int j = 0; j < sm.num_row(); ++j)
           for (int k = 0; k <= j; ++k)
             apeSaveFile << "  " << sm[j][k];  //always write full matrix
  
         apeSaveFile << std::endl;
       }
     }
     delete aliErr;
     apeSaveFile.close();
   }
   // clean up at end:  // FIXME: should we delete here or in destructor?
   delete theAlignableNavigator;
   theAlignableNavigator = nullptr;
 }