Inheritance diagram for ApeSettingAlgorithm:

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

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

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

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

virtual	~ApeSettingAlgorithm ()
	Destructor. More...

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

	AlignmentAlgorithmBase (const edm::ParameterSet &)
	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 60 of file ApeSettingAlgorithm.cc.

Constructor & Destructor Documentation

ApeSettingAlgorithm::ApeSettingAlgorithm ( const edm::ParameterSet & cfg )

Constructor.

Definition at line 97 of file ApeSettingAlgorithm.cc.

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

                                                                    :
   AlignmentAlgorithmBase(cfg), theConfig(cfg),
   theAlignableNavigator(0)
 {
   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");
   
 }

ApeSettingAlgorithm::~ApeSettingAlgorithm ( )

virtual

Destructor.

Definition at line 114 of file ApeSettingAlgorithm.cc.

References theAlignableNavigator.

 {
   delete theAlignableNavigator;
 }

Member Function Documentation

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 121 of file ApeSettingAlgorithm.cc.

References AlignableNavigator::alignableFromDetId(), edm::FileInPath::fullPath(), edm::ParameterSet::getParameter(), hcalTTPDigis_cfi::id, SequenceTypes::ignore(), readApeFromAscii_, readFullLocalMatrix_, readLocalNotGlobal_, setComposites_, theAlignableNavigator, theConfig, theTracker, mixOne_simraw_on_sim_cfi::tracker, TkRotation< T >::xx(), TkRotation< T >::xy(), TkRotation< T >::xz(), TkRotation< T >::yx(), TkRotation< T >::yy(), TkRotation< T >::yz(), TkRotation< T >::zx(), TkRotation< T >::zy(), and TkRotation< T >::zz().

 {
  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().size()<1) || 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.";
    }
 }

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

overridevirtual

Run the algorithm.

Implements AlignmentAlgorithmBase.

Definition at line 231 of file ApeSettingAlgorithm.cc.

References DEFINE_EDM_PLUGIN.

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

 {
   // nothing to do here?
 }

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

overridevirtual

Call at end of job.

Implements AlignmentAlgorithmBase.

Definition at line 194 of file ApeSettingAlgorithm.cc.

References AlignableNavigator::alignableFromDetId(), AlignableTracker::alignmentErrors(), edm::ParameterSet::getUntrackedParameter(), mps_fire::i, hcalTTPDigis_cfi::id, gen::k, AlignmentErrorsExtended::m_alignError, saveApeToAscii_, saveComposites_, saveLocalNotGlobal_, AlCaHLTBitMon_QueryRunRegistry::string, theAlignableNavigator, theConfig, theTracker, TkRotation< T >::xx(), TkRotation< T >::xy(), TkRotation< T >::xz(), TkRotation< T >::yx(), TkRotation< T >::yy(), TkRotation< T >::yz(), TkRotation< T >::zx(), TkRotation< T >::zy(), and TkRotation< T >::zz().

 {
   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().size()<1) || 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 = 0;
 }