CMS 3D CMS Logo
List of all members | Public Member Functions | Private Member Functions | Private Attributes
LaserAlignment Class Reference

#include <LaserAlignment.h>

Inheritance diagram for LaserAlignment:
edm::one::EDProducer< edm::EndRunProducer > edm::one::EDProducerBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

void beginJob () override
 
void endJob () override
 
void endRunProduce (edm::Run &, const edm::EventSetup &) override
 
 LaserAlignment (edm::ParameterSet const &theConf)
 
void produce (edm::Event &, edm::EventSetup const &) override
 
void testRoutine (void)
 for debugging & testing only, will disappear.. More...
 
 ~LaserAlignment () override
 
- Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer >
 EDProducer ()=default
 
 EDProducer (const EDProducer &)=delete
 
SerialTaskQueue * globalLuminosityBlocksQueue () final
 
SerialTaskQueue * globalRunsQueue () final
 
bool hasAbilityToProduceInBeginLumis () const final
 
bool hasAbilityToProduceInBeginProcessBlocks () const final
 
bool hasAbilityToProduceInBeginRuns () const final
 
bool hasAbilityToProduceInEndLumis () const final
 
bool hasAbilityToProduceInEndProcessBlocks () const final
 
bool hasAbilityToProduceInEndRuns () const final
 
const EDProducer & operator= (const EDProducer &)=delete
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
bool wantsInputProcessBlocks () const final
 
bool wantsProcessBlocks () const final
 
- Public Member Functions inherited from edm::one::EDProducerBase
 EDProducerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const
 
bool wantsStreamRuns () const
 
 ~EDProducerBase () override
 
- Public Member Functions inherited from edm::ProducerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
std::vector< edm::ProductResolverIndex > const & indiciesForPutProducts (BranchType iBranchType) const
 
 ProducerBase ()
 
std::vector< edm::ProductResolverIndex > const & putTokenIndexToProductResolverIndex () const
 
std::vector< bool > const & recordProvenanceList () const
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription const &)> registrationCallback () const
 used by the fwk to register list of products More...
 
void resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
 
TypeLabelList const & typeLabelList () const
 used by the fwk to register the list of products of this module More...
 
 ~ProducerBase () noexcept(false) override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
std::vector< ESProxyIndex > const & esGetTokenIndicesVector (edm::Transition iTrans) const
 
std::vector< ESRecordIndex > const & esGetTokenRecordIndicesVector (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Member Functions

void ApplyATMaskingCorrections (LASGlobalData< LASCoordinateSet > &, LASGlobalData< LASCoordinateSet > &, LASBarrelAlignmentParameterSet &)
 same for alignment tube modules More...
 
void ApplyEndcapMaskingCorrections (LASGlobalData< LASCoordinateSet > &, LASGlobalData< LASCoordinateSet > &, LASEndcapAlignmentParameterSet &)
 apply endcap correction to masked modules in TEC More...
 
void CalculateNominalCoordinates (void)
 fills a LASGlobalData<LASCoordinateSet> with nominal module positions More...
 
double ConvertAngle (double)
 convert an angle in the [-pi,pi] range to the [0,2*pi] range More...
 
void DumpHitmaps (LASGlobalData< int > &)
 for debugging only, will disappear More...
 
void DumpPosFileSet (LASGlobalData< LASCoordinateSet > &)
 for debugging only, will disappear More...
 
void DumpStripFileSet (LASGlobalData< std::pair< float, float > > &)
 for debugging only, will disappear More...
 
void fillDataProfiles (edm::Event const &, edm::EventSetup const &)
 fill profiles from SiStrip(Raw)Digi container More...
 
void fillDetectorId (void)
 fill hard coded detIds More...
 
void fillPedestalProfiles (edm::ESHandle< SiStripPedestals > &)
 fill pedestals from dbase More...
 
double getTEC2TECNominalBeamOffset (unsigned int, unsigned int, unsigned int)
 returns the nominal beam position (strips) in TEC (AT) for the profileJudge More...
 
double getTIBTOBNominalBeamOffset (unsigned int, unsigned int, unsigned int)
 returns the nominal beam position (strips) in TOB for the profileJudge More...
 
bool isATBeam (void)
 
bool isTECBeam (void)
 decide whether TEC or AT beams have fired More...
 

Private Attributes

LASGlobalData< LASModuleProfile > collectedDataProfiles
 
LASGlobalData< LASModuleProfile > currentDataProfiles
 data profiles for the current event More...
 
LASGlobalData< unsigned int > detectorId
 
bool enableJudgeZeroFilter
 config switch More...
 
bool firstEvent_
 
edm::ESHandle< GeometricDet > gD
 tracker geometry; More...
 
const edm::ESGetToken< GeometricDet, IdealGeometryRecord > geomDetToken_
 
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > geomToken_
 
const edm::ESGetToken< Alignments, GlobalPositionRcd > gprToken_
 
LASGlobalData< int > isAcceptedProfile
 
LASProfileJudge judge
 
unsigned int judgeOverdriveThreshold
 config parameters for the LASProfileJudge More...
 
bool misalignedByRefGeometry
 config switch More...
 
LASGlobalLoop moduleLoop
 
LASGlobalData< LASCoordinateSet > nominalCoordinates
 
LASGlobalData< int > numberOfAcceptedProfiles
 
double peakFinderThreshold
 config parameter More...
 
LASGlobalData< LASModuleProfile > pedestalProfiles
 
const edm::ESGetToken< PTrackerAdditionalParametersPerDet, PTrackerAdditionalParametersPerDetRcd > ptitpToken_
 
const edm::ESGetToken< PTrackerParameters, PTrackerParametersRcd > ptpToken_
 
TDirectory * singleModulesDir
 
const edm::ESGetToken< SiStripPedestals, SiStripPedestalsRcd > stripPedestalsToken_
 
LASGlobalData< TH1D * > summedHistograms
 
std::vector< unsigned int > tecDoubleHitDetId
 
AlignableTracker * theAlignableTracker
 
std::string theAlignRecordName
 
bool theApplyBeamKinkCorrections
 config switch More...
 
int theCompression
 config parameter (histograms file compression level) More...
 
std::vector< edm::ParameterSet > theDigiProducersList
 
bool theDoPedestalSubtraction
 config switch More...
 
std::string theErrorRecordName
 
int theEvents
 counter for the total number of events processed More...
 
TFile * theFile
 Tree stuff. More...
 
std::string theFileName
 config parameter (histograms file output name) More...
 
edm::ESHandle< Alignments > theGlobalPositionRcd
 
LASConstants theLasConstants
 
std::vector< unsigned int > theMaskAtModules
 
std::vector< unsigned int > theMaskTecModules
 config parameters More...
 
LASGlobalData< std::string > theProfileNames
 
bool theSaveHistograms
 config switch More...
 
bool theSetNominalStrips
 config switch More...
 
bool theStoreToDB
 config switch More...
 
edm::ESHandle< TrackerGeometry > theTrackerGeometry
 
bool theUseMinuitAlgorithm
 config switch More...
 
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > topoToken_
 Tokens for ESconsumes. More...
 
bool updateFromInputGeometry
 config switch More...
 

Additional Inherited Members

- Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase ModuleType
 
- Public Types inherited from edm::ProducerBase
using ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex > >
 
typedef ProductRegistryHelper::TypeLabelList TypeLabelList
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- Protected Member Functions inherited from edm::ProducerBase
template<Transition Tr = Transition::Event>
auto produces (std::string instanceName) noexcept
 declare what type of product will make and with which optional label More...
 
template<Transition B>
BranchAliasSetter produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
 
template<BranchType B>
BranchAliasSetter produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
 
BranchAliasSetter produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
 
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType > produces (std::string instanceName)
 
template<class ProductType >
BranchAliasSetterT< ProductType > produces ()
 
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType > produces (std::string instanceName)
 
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType > produces ()
 
template<class ProductType >
BranchAliasSetterT< ProductType > produces (std::string instanceName)
 
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType > produces ()
 
template<Transition Tr = Transition::Event>
auto produces () noexcept
 
ProducesCollector producesCollector ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B > consumes (edm::InputTag tag) noexcept
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<Transition Tr = Transition::Event>
constexpr auto esConsumes ()
 
template<Transition Tr = Transition::Event>
auto esConsumes (ESInputTag tag)
 
template<Transition Tr = Transition::Event>
ESGetTokenGeneric esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
 Used with EventSetupRecord::doGet. More...
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
void resetItemsToGetFrom (BranchType iType)
 

Detailed Description

Main reconstruction module for the Laser Alignment System

Date
2012/12/26 20:34:44
Revision
1.32
Author
Maarten Thomas
Jan Olzem

Definition at line 73 of file LaserAlignment.h.

Constructor & Destructor Documentation

◆ LaserAlignment()

LaserAlignment::LaserAlignment ( edm::ParameterSet const &  theConf)
explicit

Definition at line 21 of file LaserAlignment.cc.

References SiStripOfflineCRack_cfg::alias, gather_cfg::cout, enableJudgeZeroFilter, LASProfileJudge::EnableZeroFilter(), edm::ParameterSet::getParameter(), judge, judgeOverdriveThreshold, misalignedByRefGeometry, peakFinderThreshold, LASProfileJudge::SetOverdriveThreshold(), AlCaHLTBitMon_QueryRunRegistry::string, theApplyBeamKinkCorrections, theCompression, theDoPedestalSubtraction, theFileName, theMaskAtModules, theMaskTecModules, theSaveHistograms, theSetNominalStrips, theStoreToDB, theUseMinuitAlgorithm, and updateFromInputGeometry.

22  : topoToken_(esConsumes()),
23  geomToken_(esConsumes()),
24  geomDetToken_(esConsumes()),
25  ptpToken_(esConsumes()),
26  ptitpToken_(esConsumes()),
27  gprToken_(esConsumes()),
28  stripPedestalsToken_(esConsumes()),
29  theEvents(0),
30  theDoPedestalSubtraction(theConf.getUntrackedParameter<bool>("SubtractPedestals", true)),
31  theUseMinuitAlgorithm(theConf.getUntrackedParameter<bool>("RunMinuitAlignmentTubeAlgorithm", false)),
32  theApplyBeamKinkCorrections(theConf.getUntrackedParameter<bool>("ApplyBeamKinkCorrections", true)),
33  peakFinderThreshold(theConf.getUntrackedParameter<double>("PeakFinderThreshold", 10.)),
34  enableJudgeZeroFilter(theConf.getUntrackedParameter<bool>("EnableJudgeZeroFilter", true)),
35  judgeOverdriveThreshold(theConf.getUntrackedParameter<unsigned int>("JudgeOverdriveThreshold", 220)),
36  updateFromInputGeometry(theConf.getUntrackedParameter<bool>("UpdateFromInputGeometry", false)),
37  misalignedByRefGeometry(theConf.getUntrackedParameter<bool>("MisalignedByRefGeometry", false)),
38  theStoreToDB(theConf.getUntrackedParameter<bool>("SaveToDbase", false)),
39  theDigiProducersList(theConf.getParameter<std::vector<edm::ParameterSet> >("DigiProducersList")),
40  theSaveHistograms(theConf.getUntrackedParameter<bool>("SaveHistograms", false)),
41  theCompression(theConf.getUntrackedParameter<int>("ROOTFileCompression", 1)),
42  theFileName(theConf.getUntrackedParameter<std::string>("ROOTFileName", "test.root")),
43  theMaskTecModules(theConf.getUntrackedParameter<std::vector<unsigned int> >("MaskTECModules")),
44  theMaskAtModules(theConf.getUntrackedParameter<std::vector<unsigned int> >("MaskATModules")),
45  theSetNominalStrips(theConf.getUntrackedParameter<bool>("ForceFitterToNominalStrips", false)),
46  theLasConstants(theConf.getUntrackedParameter<std::vector<edm::ParameterSet> >("LaserAlignmentConstants")),
47  theFile(),
48  theAlignableTracker(),
49  theAlignRecordName("TrackerAlignmentRcd"),
50  theErrorRecordName("TrackerAlignmentErrorExtendedRcd"),
51  firstEvent_(true) {
52  std::cout << std::endl;
53  std::cout << "=============================================================="
54  << "\n=== LaserAlignment module configuration ==="
55  << "\n"
56  << "\n Write histograms to file = " << (theSaveHistograms ? "true" : "false")
57  << "\n Histogram file name = " << theFileName
58  << "\n Histogram file compression = " << theCompression
59  << "\n Subtract pedestals = " << (theDoPedestalSubtraction ? "true" : "false")
60  << "\n Run Minuit AT algorithm = " << (theUseMinuitAlgorithm ? "true" : "false")
61  << "\n Apply beam kink corrections = " << (theApplyBeamKinkCorrections ? "true" : "false")
62  << "\n Peak Finder Threshold = " << peakFinderThreshold
63  << "\n EnableJudgeZeroFilter = " << (enableJudgeZeroFilter ? "true" : "false")
64  << "\n JudgeOverdriveThreshold = " << judgeOverdriveThreshold
65  << "\n Update from input geometry = " << (updateFromInputGeometry ? "true" : "false")
66  << "\n Misalignment from ref geometry = " << (misalignedByRefGeometry ? "true" : "false")
67  << "\n Number of TEC modules masked = " << theMaskTecModules.size() << " (s. below list if > 0)"
68  << "\n Number of AT modules masked = " << theMaskAtModules.size() << " (s. below list if > 0)"
69  << "\n Store to database = " << (theStoreToDB ? "true" : "false")
70  << "\n ----------------------------------------------- ----------"
71  << (theSetNominalStrips ? "\n Set strips to nominal = true" : "\n")
72  << "\n=============================================================" << std::endl;
73 
74  // tell about masked modules
75  if (!theMaskTecModules.empty()) {
76  std::cout << " ===============================================================================================\n"
77  << std::flush;
78  std::cout << " The following " << theMaskTecModules.size()
79  << " TEC modules have been masked out and will not be considered by the TEC algorithm:\n " << std::flush;
80  for (std::vector<unsigned int>::iterator moduleIt = theMaskTecModules.begin(); moduleIt != theMaskTecModules.end();
81  ++moduleIt) {
82  std::cout << *moduleIt << (moduleIt != --theMaskTecModules.end() ? ", " : "") << std::flush;
83  }
84  std::cout << std::endl << std::flush;
85  std::cout << " ===============================================================================================\n\n"
86  << std::flush;
87  }
88  if (!theMaskAtModules.empty()) {
89  std::cout << " ===============================================================================================\n"
90  << std::flush;
91  std::cout << " The following " << theMaskAtModules.size()
92  << " AT modules have been masked out and will not be considered by the AT algorithm:\n " << std::flush;
93  for (std::vector<unsigned int>::iterator moduleIt = theMaskAtModules.begin(); moduleIt != theMaskAtModules.end();
94  ++moduleIt) {
95  std::cout << *moduleIt << (moduleIt != --theMaskAtModules.end() ? ", " : "") << std::flush;
96  }
97  std::cout << std::endl << std::flush;
98  std::cout << " ===============================================================================================\n\n"
99  << std::flush;
100  }
101 
102  // alias for the Branches in the root files
103  std::string alias(theConf.getParameter<std::string>("@module_label"));
104 
105  // declare the product to produce
106  produces<TkLasBeamCollection, edm::Transition::EndRun>("tkLaserBeams").setBranchAlias(alias + "TkLasBeamCollection");
107 
108  // switch judge's zero filter depending on cfg
109  judge.EnableZeroFilter(enableJudgeZeroFilter);
110 
111  // set the upper threshold for zero suppressed data
112  judge.SetOverdriveThreshold(judgeOverdriveThreshold);
113 }
LaserAlignment::judgeOverdriveThreshold
unsigned int judgeOverdriveThreshold
config parameters for the LASProfileJudge
Definition: LaserAlignment.h:158
LaserAlignment::ptpToken_
const edm::ESGetToken< PTrackerParameters, PTrackerParametersRcd > ptpToken_
Definition: LaserAlignment.h:134
LaserAlignment::theSaveHistograms
bool theSaveHistograms
config switch
Definition: LaserAlignment.h:173
LaserAlignment::geomDetToken_
const edm::ESGetToken< GeometricDet, IdealGeometryRecord > geomDetToken_
Definition: LaserAlignment.h:133
LaserAlignment::stripPedestalsToken_
const edm::ESGetToken< SiStripPedestals, SiStripPedestalsRcd > stripPedestalsToken_
Definition: LaserAlignment.h:137
LaserAlignment::misalignedByRefGeometry
bool misalignedByRefGeometry
config switch
Definition: LaserAlignment.h:164
LaserAlignment::theAlignRecordName
std::string theAlignRecordName
Definition: LaserAlignment.h:242
LaserAlignment::ptitpToken_
const edm::ESGetToken< PTrackerAdditionalParametersPerDet, PTrackerAdditionalParametersPerDetRcd > ptitpToken_
Definition: LaserAlignment.h:135
LaserAlignment::peakFinderThreshold
double peakFinderThreshold
config parameter
Definition: LaserAlignment.h:152
SiStripOfflineCRack_cfg.alias
alias
Definition: SiStripOfflineCRack_cfg.py:128
LaserAlignment::theDoPedestalSubtraction
bool theDoPedestalSubtraction
config switch
Definition: LaserAlignment.h:143
LaserAlignment::theMaskTecModules
std::vector< unsigned int > theMaskTecModules
config parameters
Definition: LaserAlignment.h:182
LaserAlignment::theUseMinuitAlgorithm
bool theUseMinuitAlgorithm
config switch
Definition: LaserAlignment.h:146
LaserAlignment::theFile
TFile * theFile
Tree stuff.
Definition: LaserAlignment.h:232
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
LaserAlignment::theDigiProducersList
std::vector< edm::ParameterSet > theDigiProducersList
Definition: LaserAlignment.h:170
LaserAlignment::geomToken_
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > geomToken_
Definition: LaserAlignment.h:132
LaserAlignment::theFileName
std::string theFileName
config parameter (histograms file output name)
Definition: LaserAlignment.h:179
LaserAlignment::topoToken_
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > topoToken_
Tokens for ESconsumes.
Definition: LaserAlignment.h:131
LaserAlignment::theApplyBeamKinkCorrections
bool theApplyBeamKinkCorrections
config switch
Definition: LaserAlignment.h:149
LaserAlignment::theErrorRecordName
std::string theErrorRecordName
Definition: LaserAlignment.h:242
LaserAlignment::theStoreToDB
bool theStoreToDB
config switch
Definition: LaserAlignment.h:167
LaserAlignment::theAlignableTracker
AlignableTracker * theAlignableTracker
Definition: LaserAlignment.h:240
LaserAlignment::theEvents
int theEvents
counter for the total number of events processed
Definition: LaserAlignment.h:140
LaserAlignment::theLasConstants
LASConstants theLasConstants
Definition: LaserAlignment.h:193
LASProfileJudge::SetOverdriveThreshold
void SetOverdriveThreshold(unsigned int)
Definition: LASProfileJudge.cc:90
LaserAlignment::theMaskAtModules
std::vector< unsigned int > theMaskAtModules
Definition: LaserAlignment.h:183
LaserAlignment::theCompression
int theCompression
config parameter (histograms file compression level)
Definition: LaserAlignment.h:176
LASProfileJudge::EnableZeroFilter
void EnableZeroFilter(bool)
Definition: LASProfileJudge.cc:79
LaserAlignment::gprToken_
const edm::ESGetToken< Alignments, GlobalPositionRcd > gprToken_
Definition: LaserAlignment.h:136
gather_cfg.cout
cout
Definition: gather_cfg.py:144
LaserAlignment::enableJudgeZeroFilter
bool enableJudgeZeroFilter
config switch
Definition: LaserAlignment.h:155
LaserAlignment::judge
LASProfileJudge judge
Definition: LaserAlignment.h:190
LaserAlignment::updateFromInputGeometry
bool updateFromInputGeometry
config switch
Definition: LaserAlignment.h:161
LaserAlignment::theSetNominalStrips
bool theSetNominalStrips
config switch
Definition: LaserAlignment.h:186

◆ ~LaserAlignment()

LaserAlignment::~LaserAlignment ( )
override

Definition at line 118 of file LaserAlignment.cc.

References theAlignableTracker, theFile, and theSaveHistograms.

118  {
119  if (theSaveHistograms)
120  theFile->Write();
121  if (theFile) {
122  delete theFile;
123  }
124  if (theAlignableTracker) {
125  delete theAlignableTracker;
126  }
127 }
LaserAlignment::theSaveHistograms
bool theSaveHistograms
config switch
Definition: LaserAlignment.h:173
LaserAlignment::theFile
TFile * theFile
Tree stuff.
Definition: LaserAlignment.h:232
LaserAlignment::theAlignableTracker
AlignableTracker * theAlignableTracker
Definition: LaserAlignment.h:240

Member Function Documentation

◆ ApplyATMaskingCorrections()

void LaserAlignment::ApplyATMaskingCorrections ( LASGlobalData< LASCoordinateSet > &  measuredCoordinates,
LASGlobalData< LASCoordinateSet > &  nominalCoordinates,
LASBarrelAlignmentParameterSet &  atParameters 
)
private

same for alignment tube modules

loop the list of alignment tube modules to be masked and apply the corrections from the "barrelParameters" to them

Definition at line 1625 of file LaserAlignment.cc.

References EcalCondDBWriter_cfi::beam, detectorId, LASCoordinateSet::GetPhi(), LASAlignmentTubeAlgorithm::GetTEC2TECAlignmentParameterCorrection(), LASGlobalData< T >::GetTEC2TECEntry(), LASAlignmentTubeAlgorithm::GetTIBTOBAlignmentParameterCorrection(), LASGlobalData< T >::GetTIBTOBEntry(), moduleLoop, nominalCoordinates, LASCoordinateSet::SetPhi(), LASGlobalLoop::TEC2TECLoop(), theMaskAtModules, and LASGlobalLoop::TIBTOBLoop().

Referenced by endRunProduce().

1627  {
1628  // loop the list of modules to be masked
1629  for (std::vector<unsigned int>::iterator moduleIt = theMaskAtModules.begin(); moduleIt != theMaskAtModules.end();
1630  ++moduleIt) {
1631  // loop variables
1632  LASGlobalLoop moduleLoop;
1633  int det, beam, disk, pos;
1634 
1635  // this will calculate the corrections from the alignment parameters
1636  LASAlignmentTubeAlgorithm atAlgorithm;
1637 
1638  // find the location of the respective module in the container with these loops:
1639 
1640  // first TIB+TOB
1641  det = 2;
1642  beam = 0;
1643  pos = 0;
1644  do {
1645  // here we got it
1646  if (detectorId.GetTIBTOBEntry(det, beam, pos) == *moduleIt) {
1647  // the nominal phi value for this module
1648  const double nominalPhi = nominalCoordinates.GetTIBTOBEntry(det, beam, pos).GetPhi();
1649 
1650  // the offset from the alignment parameters
1651  const double phiCorrection =
1652  atAlgorithm.GetTIBTOBAlignmentParameterCorrection(det, beam, pos, nominalCoordinates, atParameters);
1653 
1654  // apply the corrections
1655  measuredCoordinates.GetTIBTOBEntry(det, beam, pos).SetPhi(nominalPhi - phiCorrection);
1656  }
1657 
1658  } while (moduleLoop.TIBTOBLoop(det, beam, pos));
1659 
1660  // then TEC(AT)
1661  det = 0;
1662  beam = 0;
1663  disk = 0;
1664  do {
1665  // here we got it
1666  if (detectorId.GetTEC2TECEntry(det, beam, disk) == *moduleIt) {
1667  // the nominal phi value for this module
1668  const double nominalPhi = nominalCoordinates.GetTEC2TECEntry(det, beam, disk).GetPhi();
1669 
1670  // the offset from the alignment parameters
1671  const double phiCorrection =
1672  atAlgorithm.GetTEC2TECAlignmentParameterCorrection(det, beam, disk, nominalCoordinates, atParameters);
1673 
1674  // apply the corrections
1675  measuredCoordinates.GetTEC2TECEntry(det, beam, disk).SetPhi(nominalPhi - phiCorrection);
1676  }
1677 
1678  } while (moduleLoop.TEC2TECLoop(det, beam, disk));
1679  }
1680 }
LASCoordinateSet::SetPhi
void SetPhi(double aPhi)
Definition: LASCoordinateSet.h:27
LASGlobalLoop::TEC2TECLoop
bool TEC2TECLoop(int &, int &, int &) const
Definition: LASGlobalLoop.cc:74
EcalCondDBWriter_cfi.beam
beam
Definition: EcalCondDBWriter_cfi.py:45
LASGlobalData::GetTIBTOBEntry
T & GetTIBTOBEntry(int subdetector, int beam, int tibTobPosition)
Definition: LASGlobalData.h:107
LASAlignmentTubeAlgorithm::GetTEC2TECAlignmentParameterCorrection
double GetTEC2TECAlignmentParameterCorrection(int, int, int, LASGlobalData< LASCoordinateSet > &, LASBarrelAlignmentParameterSet &)
Definition: LASAlignmentTubeAlgorithm.cc:609
LASGlobalLoop::TIBTOBLoop
bool TIBTOBLoop(int &, int &, int &) const
Definition: LASGlobalLoop.cc:46
LASCoordinateSet::GetPhi
double GetPhi(void) const
Definition: LASCoordinateSet.h:18
LASGlobalData::GetTEC2TECEntry
T & GetTEC2TECEntry(int subdetector, int beam, int tecDisk)
Definition: LASGlobalData.h:129
LaserAlignment::nominalCoordinates
LASGlobalData< LASCoordinateSet > nominalCoordinates
Definition: LaserAlignment.h:225
LaserAlignment::theMaskAtModules
std::vector< unsigned int > theMaskAtModules
Definition: LaserAlignment.h:183
LaserAlignment::detectorId
LASGlobalData< unsigned int > detectorId
Definition: LaserAlignment.h:196
LASAlignmentTubeAlgorithm::GetTIBTOBAlignmentParameterCorrection
double GetTIBTOBAlignmentParameterCorrection(int, int, int, LASGlobalData< LASCoordinateSet > &, LASBarrelAlignmentParameterSet &)
Definition: LASAlignmentTubeAlgorithm.cc:490
LaserAlignment::moduleLoop
LASGlobalLoop moduleLoop
Definition: LaserAlignment.h:229

◆ ApplyEndcapMaskingCorrections()

void LaserAlignment::ApplyEndcapMaskingCorrections ( LASGlobalData< LASCoordinateSet > &  measuredCoordinates,
LASGlobalData< LASCoordinateSet > &  nominalCoordinates,
LASEndcapAlignmentParameterSet &  endcapParameters 
)
private

apply endcap correction to masked modules in TEC

loop the list of endcap modules to be masked and apply the corrections from the "endcapParameters" to them

Definition at line 1585 of file LaserAlignment.cc.

References EcalCondDBWriter_cfi::beam, detectorId, LASEndcapAlgorithm::GetAlignmentParameterCorrection(), LASCoordinateSet::GetPhi(), LASGlobalData< T >::GetTECEntry(), moduleLoop, nominalCoordinates, relativeConstraints::ring, LASCoordinateSet::SetPhi(), LASGlobalLoop::TECLoop(), and theMaskTecModules.

Referenced by endRunProduce().

1587  {
1588  // loop the list of modules to be masked
1589  for (std::vector<unsigned int>::iterator moduleIt = theMaskTecModules.begin(); moduleIt != theMaskTecModules.end();
1590  ++moduleIt) {
1591  // loop variables
1592  LASGlobalLoop moduleLoop;
1593  int det, ring, beam, disk;
1594 
1595  // this will calculate the corrections from the alignment parameters
1596  LASEndcapAlgorithm endcapAlgorithm;
1597 
1598  // find the location of the respective module in the container with this loop
1599  det = 0;
1600  ring = 0;
1601  beam = 0;
1602  disk = 0;
1603  do {
1604  // here we got it
1605  if (detectorId.GetTECEntry(det, ring, beam, disk) == *moduleIt) {
1606  // the nominal phi value for this module
1607  const double nominalPhi = nominalCoordinates.GetTECEntry(det, ring, beam, disk).GetPhi();
1608 
1609  // the offset from the alignment parameters
1610  const double phiCorrection = endcapAlgorithm.GetAlignmentParameterCorrection(
1611  det, ring, beam, disk, nominalCoordinates, endcapParameters);
1612 
1613  // apply the corrections
1614  measuredCoordinates.GetTECEntry(det, ring, beam, disk).SetPhi(nominalPhi - phiCorrection);
1615  }
1616 
1617  } while (moduleLoop.TECLoop(det, ring, beam, disk));
1618  }
1619 }
LASCoordinateSet::SetPhi
void SetPhi(double aPhi)
Definition: LASCoordinateSet.h:27
LaserAlignment::theMaskTecModules
std::vector< unsigned int > theMaskTecModules
config parameters
Definition: LaserAlignment.h:182
EcalCondDBWriter_cfi.beam
beam
Definition: EcalCondDBWriter_cfi.py:45
LASGlobalLoop::TECLoop
bool TECLoop(int &, int &, int &, int &) const
Definition: LASGlobalLoop.cc:13
LASCoordinateSet::GetPhi
double GetPhi(void) const
Definition: LASCoordinateSet.h:18
LaserAlignment::nominalCoordinates
LASGlobalData< LASCoordinateSet > nominalCoordinates
Definition: LaserAlignment.h:225
LASGlobalData::GetTECEntry
T & GetTECEntry(int subdetector, int tecRing, int beam, int tecDisk)
Definition: LASGlobalData.h:84
LASEndcapAlgorithm::GetAlignmentParameterCorrection
double GetAlignmentParameterCorrection(int, int, int, int, LASGlobalData< LASCoordinateSet > &, LASEndcapAlignmentParameterSet &)
Definition: LASEndcapAlgorithm.cc:312
LaserAlignment::detectorId
LASGlobalData< unsigned int > detectorId
Definition: LaserAlignment.h:196
LaserAlignment::moduleLoop
LASGlobalLoop moduleLoop
Definition: LaserAlignment.h:229

◆ beginJob()

void LaserAlignment::beginJob ( )
overridevirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 132 of file LaserAlignment.cc.

References EcalCondDBWriter_cfi::beam, collectedDataProfiles, currentDataProfiles, Exception, fillDetectorId(), firstEvent_, LASGlobalData< T >::GetTEC2TECEntry(), LASGlobalData< T >::GetTECEntry(), LASGlobalData< T >::GetTIBTOBEntry(), isAcceptedProfile, moduleLoop, numberOfAcceptedProfiles, pedestalProfiles, relativeConstraints::ring, LASModuleProfile::SetAllValuesTo(), LASGlobalData< T >::SetTEC2TECEntry(), LASGlobalData< T >::SetTECEntry(), LASGlobalData< T >::SetTIBTOBEntry(), singleModulesDir, summedHistograms, LASGlobalLoop::TEC2TECLoop(), LASGlobalLoop::TECLoop(), theCompression, theFile, theFileName, theProfileNames, theSaveHistograms, and LASGlobalLoop::TIBTOBLoop().

132  {
133  // write sumed histograms to file (if selected in cfg)
134  if (theSaveHistograms) {
135  // creating a new file
136  theFile = new TFile(theFileName.c_str(), "RECREATE", "CMS ROOT file");
137 
138  // initialize the histograms
139  if (theFile) {
140  theFile->SetCompressionLevel(theCompression);
141  singleModulesDir = theFile->mkdir("single modules");
142  } else
143  throw cms::Exception(" [LaserAlignment::beginJob]")
144  << " ** ERROR: could not open file:" << theFileName.c_str() << " for writing." << std::endl;
145  }
146 
147  // detector id maps (hard coded)
148  fillDetectorId();
149 
150  // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
151  // PROFILE, HISTOGRAM & FITFUNCTION INITIALIZATION
152  // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
153 
154  // object used to build various strings for names and labels
155  std::stringstream nameBuilder;
156 
157  // loop variables for use with LASGlobalLoop object
158  int det, ring, beam, disk, pos;
159 
160  // loop TEC modules
161  det = 0;
162  ring = 0;
163  beam = 0;
164  disk = 0;
165  do { // loop using LASGlobalLoop functionality
166  // init the profiles
167  pedestalProfiles.GetTECEntry(det, ring, beam, disk).SetAllValuesTo(0.);
168  currentDataProfiles.GetTECEntry(det, ring, beam, disk).SetAllValuesTo(0.);
169  collectedDataProfiles.GetTECEntry(det, ring, beam, disk).SetAllValuesTo(0.);
170 
171  // init the hit maps
172  isAcceptedProfile.SetTECEntry(det, ring, beam, disk, 0);
173  numberOfAcceptedProfiles.SetTECEntry(det, ring, beam, disk, 0);
174 
175  // create strings for histo names
176  nameBuilder.clear();
177  nameBuilder.str("");
178  nameBuilder << "TEC";
179  if (det == 0)
180  nameBuilder << "+";
181  else
182  nameBuilder << "-";
183  nameBuilder << "_Ring";
184  if (ring == 0)
185  nameBuilder << "4";
186  else
187  nameBuilder << "6";
188  nameBuilder << "_Beam" << beam;
189  nameBuilder << "_Disk" << disk;
190  theProfileNames.SetTECEntry(det, ring, beam, disk, nameBuilder.str());
191 
192  // init the histograms
193  if (theSaveHistograms) {
194  nameBuilder << "_Histo";
195  summedHistograms.SetTECEntry(
196  det, ring, beam, disk, new TH1D(nameBuilder.str().c_str(), nameBuilder.str().c_str(), 512, 0, 512));
197  summedHistograms.GetTECEntry(det, ring, beam, disk)->SetDirectory(singleModulesDir);
198  }
199 
200  } while (moduleLoop.TECLoop(det, ring, beam, disk));
201 
202  // TIB & TOB section
203  det = 2;
204  beam = 0;
205  pos = 0;
206  do { // loop using LASGlobalLoop functionality
207  // init the profiles
208  pedestalProfiles.GetTIBTOBEntry(det, beam, pos).SetAllValuesTo(0.);
209  currentDataProfiles.GetTIBTOBEntry(det, beam, pos).SetAllValuesTo(0.);
210  collectedDataProfiles.GetTIBTOBEntry(det, beam, pos).SetAllValuesTo(0.);
211 
212  // init the hit maps
213  isAcceptedProfile.SetTIBTOBEntry(det, beam, pos, 0);
214  numberOfAcceptedProfiles.SetTIBTOBEntry(det, beam, pos, 0);
215 
216  // create strings for histo names
217  nameBuilder.clear();
218  nameBuilder.str("");
219  if (det == 2)
220  nameBuilder << "TIB";
221  else
222  nameBuilder << "TOB";
223  nameBuilder << "_Beam" << beam;
224  nameBuilder << "_Zpos" << pos;
225 
226  theProfileNames.SetTIBTOBEntry(det, beam, pos, nameBuilder.str());
227 
228  // init the histograms
229  if (theSaveHistograms) {
230  nameBuilder << "_Histo";
231  summedHistograms.SetTIBTOBEntry(
232  det, beam, pos, new TH1D(nameBuilder.str().c_str(), nameBuilder.str().c_str(), 512, 0, 512));
233  summedHistograms.GetTIBTOBEntry(det, beam, pos)->SetDirectory(singleModulesDir);
234  }
235 
236  } while (moduleLoop.TIBTOBLoop(det, beam, pos));
237 
238  // TEC2TEC AT section
239  det = 0;
240  beam = 0;
241  disk = 0;
242  do { // loop using LASGlobalLoop functionality
243  // init the profiles
244  pedestalProfiles.GetTEC2TECEntry(det, beam, disk).SetAllValuesTo(0.);
245  currentDataProfiles.GetTEC2TECEntry(det, beam, disk).SetAllValuesTo(0.);
246  collectedDataProfiles.GetTEC2TECEntry(det, beam, disk).SetAllValuesTo(0.);
247 
248  // init the hit maps
249  isAcceptedProfile.SetTEC2TECEntry(det, beam, disk, 0);
250  numberOfAcceptedProfiles.SetTEC2TECEntry(det, beam, disk, 0);
251 
252  // create strings for histo names
253  nameBuilder.clear();
254  nameBuilder.str("");
255  nameBuilder << "TEC(AT)";
256  if (det == 0)
257  nameBuilder << "+";
258  else
259  nameBuilder << "-";
260  nameBuilder << "_Beam" << beam;
261  nameBuilder << "_Disk" << disk;
262  theProfileNames.SetTEC2TECEntry(det, beam, disk, nameBuilder.str());
263 
264  // init the histograms
265  if (theSaveHistograms) {
266  nameBuilder << "_Histo";
267  summedHistograms.SetTEC2TECEntry(
268  det, beam, disk, new TH1D(nameBuilder.str().c_str(), nameBuilder.str().c_str(), 512, 0, 512));
269  summedHistograms.GetTEC2TECEntry(det, beam, disk)->SetDirectory(singleModulesDir);
270  }
271 
272  } while (moduleLoop.TEC2TECLoop(det, beam, disk));
273 
274  firstEvent_ = true;
275 }
LASGlobalData::SetTEC2TECEntry
void SetTEC2TECEntry(int subdetector, int beam, int tecDisk, T)
Definition: LASGlobalData.h:195
LaserAlignment::theSaveHistograms
bool theSaveHistograms
config switch
Definition: LaserAlignment.h:173
LASGlobalLoop::TEC2TECLoop
bool TEC2TECLoop(int &, int &, int &) const
Definition: LASGlobalLoop.cc:74
LaserAlignment::numberOfAcceptedProfiles
LASGlobalData< int > numberOfAcceptedProfiles
Definition: LaserAlignment.h:214
LASModuleProfile::SetAllValuesTo
void SetAllValuesTo(const double &)
Definition: LASModuleProfile.cc:76
LaserAlignment::currentDataProfiles
LASGlobalData< LASModuleProfile > currentDataProfiles
data profiles for the current event
Definition: LaserAlignment.h:205
LaserAlignment::theFile
TFile * theFile
Tree stuff.
Definition: LaserAlignment.h:232
LaserAlignment::isAcceptedProfile
LASGlobalData< int > isAcceptedProfile
Definition: LaserAlignment.h:218
EcalCondDBWriter_cfi.beam
beam
Definition: EcalCondDBWriter_cfi.py:45
LaserAlignment::summedHistograms
LASGlobalData< TH1D * > summedHistograms
Definition: LaserAlignment.h:222
LaserAlignment::singleModulesDir
TDirectory * singleModulesDir
Definition: LaserAlignment.h:233
LASGlobalLoop::TECLoop
bool TECLoop(int &, int &, int &, int &) const
Definition: LASGlobalLoop.cc:13
LaserAlignment::theFileName
std::string theFileName
config parameter (histograms file output name)
Definition: LaserAlignment.h:179
LASGlobalData::GetTIBTOBEntry
T & GetTIBTOBEntry(int subdetector, int beam, int tibTobPosition)
Definition: LASGlobalData.h:107
LASGlobalData::SetTIBTOBEntry
void SetTIBTOBEntry(int subdetector, int beam, int tibTobPosition, T)
Definition: LASGlobalData.h:173
LASGlobalLoop::TIBTOBLoop
bool TIBTOBLoop(int &, int &, int &) const
Definition: LASGlobalLoop.cc:46
LaserAlignment::pedestalProfiles
LASGlobalData< LASModuleProfile > pedestalProfiles
Definition: LaserAlignment.h:202
LASGlobalData::GetTEC2TECEntry
T & GetTEC2TECEntry(int subdetector, int beam, int tecDisk)
Definition: LASGlobalData.h:129
LASGlobalData::GetTECEntry
T & GetTECEntry(int subdetector, int tecRing, int beam, int tecDisk)
Definition: LASGlobalData.h:84
LaserAlignment::theProfileNames
LASGlobalData< std::string > theProfileNames
Definition: LaserAlignment.h:211
LaserAlignment::collectedDataProfiles
LASGlobalData< LASModuleProfile > collectedDataProfiles
Definition: LaserAlignment.h:208
LaserAlignment::fillDetectorId
void fillDetectorId(void)
fill hard coded detIds
Definition: LaserAlignmentFillDetId.cc:7
LaserAlignment::theCompression
int theCompression
config parameter (histograms file compression level)
Definition: LaserAlignment.h:176
LASGlobalData::SetTECEntry
void SetTECEntry(int subdetector, int tecRing, int beam, int tecDisk, T)
Definition: LASGlobalData.h:150
LaserAlignment::moduleLoop
LASGlobalLoop moduleLoop
Definition: LaserAlignment.h:229

◆ CalculateNominalCoordinates()

void LaserAlignment::CalculateNominalCoordinates ( void  )
private

fills a LASGlobalData<LASCoordinateSet> with nominal module positions

Definition at line 1332 of file LaserAlignment.cc.

References EcalCondDBWriter_cfi::beam, moduleLoop, nominalCoordinates, relativeConstraints::ring, LASGlobalData< T >::SetTEC2TECEntry(), LASGlobalData< T >::SetTECEntry(), LASGlobalData< T >::SetTIBTOBEntry(), LASGlobalLoop::TEC2TECLoop(), LASGlobalLoop::TECLoop(), and LASGlobalLoop::TIBTOBLoop().

Referenced by endRunProduce().

1332  {
1333  //
1334  // hard coded data yet...
1335  //
1336 
1337  // nominal phi values of tec beam / alignment tube hits (parameter is beam 0-7)
1338  const double tecPhiPositions[8] = {
1339  0.392699, 1.178097, 1.963495, 2.748894, 3.534292, 4.319690, 5.105088, 5.890486}; // new values calculated by maple
1340  const double atPhiPositions[8] = {
1341  0.392699, 1.289799, 1.851794, 2.748894, 3.645995, 4.319690, 5.216791, 5.778784}; // new values calculated by maple
1342 
1343  // nominal r values (mm) of hits
1344  const double tobRPosition = 600.;
1345  const double tibRPosition = 514.;
1346  const double tecRPosition[2] = {564., 840.}; // ring 4,6
1347 
1348  // nominal z values (mm) of hits in barrel (parameter is pos 0-6)
1349  const double tobZPosition[6] = {1040., 580., 220., -140., -500., -860.};
1350  const double tibZPosition[6] = {620., 380., 180., -100., -340., -540.};
1351 
1352  // nominal z values (mm) of hits in tec (parameter is disk 0-8); FOR TEC-: (* -1.)
1353  const double tecZPosition[9] = {1322.5, 1462.5, 1602.5, 1742.5, 1882.5, 2057.5, 2247.5, 2452.5, 2667.5};
1354 
1355  //
1356  // now we fill these into the nominalCoordinates container;
1357  // errors are zero for nominal values..
1358  //
1359 
1360  // loop object and its variables
1361  LASGlobalLoop moduleLoop;
1362  int det, ring, beam, disk, pos;
1363 
1364  // TEC+- section
1365  det = 0;
1366  ring = 0, beam = 0;
1367  disk = 0;
1368  do {
1369  if (det == 0) { // this is TEC+
1370  nominalCoordinates.SetTECEntry(
1371  det,
1372  ring,
1373  beam,
1374  disk,
1375  LASCoordinateSet(tecPhiPositions[beam], 0., tecRPosition[ring], 0., tecZPosition[disk], 0.));
1376  } else { // now TEC-
1377  nominalCoordinates.SetTECEntry(
1378  det,
1379  ring,
1380  beam,
1381  disk,
1382  LASCoordinateSet(
1383  tecPhiPositions[beam], 0., tecRPosition[ring], 0., -1. * tecZPosition[disk], 0.)); // just * -1.
1384  }
1385 
1386  } while (moduleLoop.TECLoop(det, ring, beam, disk));
1387 
1388  // TIB & TOB section
1389  det = 2;
1390  beam = 0;
1391  pos = 0;
1392  do {
1393  if (det == 2) { // this is TIB
1394  nominalCoordinates.SetTIBTOBEntry(
1395  det, beam, pos, LASCoordinateSet(atPhiPositions[beam], 0., tibRPosition, 0., tibZPosition[pos], 0.));
1396  } else { // now TOB
1397  nominalCoordinates.SetTIBTOBEntry(
1398  det, beam, pos, LASCoordinateSet(atPhiPositions[beam], 0., tobRPosition, 0., tobZPosition[pos], 0.));
1399  }
1400 
1401  } while (moduleLoop.TIBTOBLoop(det, beam, pos));
1402 
1403  // TEC2TEC AT section
1404  det = 0;
1405  beam = 0;
1406  disk = 0;
1407  do {
1408  if (det == 0) { // this is TEC+, ring4 only
1409  nominalCoordinates.SetTEC2TECEntry(
1410  det, beam, disk, LASCoordinateSet(atPhiPositions[beam], 0., tecRPosition[0], 0., tecZPosition[disk], 0.));
1411  } else { // now TEC-
1412  nominalCoordinates.SetTEC2TECEntry(
1413  det,
1414  beam,
1415  disk,
1416  LASCoordinateSet(atPhiPositions[beam], 0., tecRPosition[0], 0., -1. * tecZPosition[disk], 0.)); // just * -1.
1417  }
1418 
1419  } while (moduleLoop.TEC2TECLoop(det, beam, disk));
1420 }
LASGlobalData::SetTEC2TECEntry
void SetTEC2TECEntry(int subdetector, int beam, int tecDisk, T)
Definition: LASGlobalData.h:195
LASGlobalLoop::TEC2TECLoop
bool TEC2TECLoop(int &, int &, int &) const
Definition: LASGlobalLoop.cc:74
EcalCondDBWriter_cfi.beam
beam
Definition: EcalCondDBWriter_cfi.py:45
LASGlobalLoop::TECLoop
bool TECLoop(int &, int &, int &, int &) const
Definition: LASGlobalLoop.cc:13
LASGlobalData::SetTIBTOBEntry
void SetTIBTOBEntry(int subdetector, int beam, int tibTobPosition, T)
Definition: LASGlobalData.h:173
LASGlobalLoop::TIBTOBLoop
bool TIBTOBLoop(int &, int &, int &) const
Definition: LASGlobalLoop.cc:46
LaserAlignment::nominalCoordinates
LASGlobalData< LASCoordinateSet > nominalCoordinates
Definition: LaserAlignment.h:225
LASGlobalData::SetTECEntry
void SetTECEntry(int subdetector, int tecRing, int beam, int tecDisk, T)
Definition: LASGlobalData.h:150
LaserAlignment::moduleLoop
LASGlobalLoop moduleLoop
Definition: LaserAlignment.h:229

◆ ConvertAngle()

double LaserAlignment::ConvertAngle ( double  angle)
private

convert an angle in the [-pi,pi] range to the [0,2*pi] range

convert an angle in the [-pi,pi] range to the [0,2*pi] range

Definition at line 1426 of file LaserAlignment.cc.

References angle(), Exception, and M_PI.

Referenced by endRunProduce().

1426  {
1427  if (angle < -1. * M_PI || angle > M_PI) {
1428  throw cms::Exception(" [LaserAlignment::ConvertAngle] ")
1429  << "** ERROR: Called with illegal input angle: " << angle << "." << std::endl;
1430  }
1431 
1432  if (angle >= 0.)
1433  return angle;
1434  else
1435  return (angle + 2. * M_PI);
1436 }
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:50
angle
T angle(T x1, T y1, T z1, T x2, T y2, T z2)
Definition: angle.h:11

◆ DumpHitmaps()

void LaserAlignment::DumpHitmaps ( LASGlobalData< int > &  numberOfAcceptedProfiles)
private

for debugging only, will disappear

Definition at line 1531 of file LaserAlignment.cc.

References EcalCondDBWriter_cfi::beam, gather_cfg::cout, LASGlobalData< T >::GetTECEntry(), and numberOfAcceptedProfiles.

Referenced by endRunProduce().

1531  {
1532  std::cout << " [LaserAlignment::DumpHitmaps] -- Dumping hitmap for TEC+:" << std::endl;
1533  std::cout << " [LaserAlignment::DumpHitmaps] -- Ring4:" << std::endl;
1534  std::cout << " disk0 disk1 disk2 disk3 disk4 disk5 disk6 disk7 disk8" << std::endl;
1535 
1536  for (int beam = 0; beam < 8; ++beam) {
1537  std::cout << " beam" << beam << ":";
1538  for (int disk = 0; disk < 9; ++disk) {
1539  std::cout << "\t" << numberOfAcceptedProfiles.GetTECEntry(0, 0, beam, disk);
1540  }
1541  std::cout << std::endl;
1542  }
1543 
1544  std::cout << " [LaserAlignment::DumpHitmaps] -- Ring6:" << std::endl;
1545  std::cout << " disk0 disk1 disk2 disk3 disk4 disk5 disk6 disk7 disk8" << std::endl;
1546 
1547  for (int beam = 0; beam < 8; ++beam) {
1548  std::cout << " beam" << beam << ":";
1549  for (int disk = 0; disk < 9; ++disk) {
1550  std::cout << "\t" << numberOfAcceptedProfiles.GetTECEntry(0, 1, beam, disk);
1551  }
1552  std::cout << std::endl;
1553  }
1554 
1555  std::cout << " [LaserAlignment::DumpHitmaps] -- Dumping hitmap for TEC-:" << std::endl;
1556  std::cout << " [LaserAlignment::DumpHitmaps] -- Ring4:" << std::endl;
1557  std::cout << " disk0 disk1 disk2 disk3 disk4 disk5 disk6 disk7 disk8" << std::endl;
1558 
1559  for (int beam = 0; beam < 8; ++beam) {
1560  std::cout << " beam" << beam << ":";
1561  for (int disk = 0; disk < 9; ++disk) {
1562  std::cout << "\t" << numberOfAcceptedProfiles.GetTECEntry(1, 0, beam, disk);
1563  }
1564  std::cout << std::endl;
1565  }
1566 
1567  std::cout << " [LaserAlignment::DumpHitmaps] -- Ring6:" << std::endl;
1568  std::cout << " disk0 disk1 disk2 disk3 disk4 disk5 disk6 disk7 disk8" << std::endl;
1569 
1570  for (int beam = 0; beam < 8; ++beam) {
1571  std::cout << " beam" << beam << ":";
1572  for (int disk = 0; disk < 9; ++disk) {
1573  std::cout << "\t" << numberOfAcceptedProfiles.GetTECEntry(1, 1, beam, disk);
1574  }
1575  std::cout << std::endl;
1576  }
1577 
1578  std::cout << " [LaserAlignment::DumpHitmaps] -- End of dump." << std::endl << std::endl;
1579 }
LaserAlignment::numberOfAcceptedProfiles
LASGlobalData< int > numberOfAcceptedProfiles
Definition: LaserAlignment.h:214
EcalCondDBWriter_cfi.beam
beam
Definition: EcalCondDBWriter_cfi.py:45
LASGlobalData::GetTECEntry
T & GetTECEntry(int subdetector, int tecRing, int beam, int tecDisk)
Definition: LASGlobalData.h:84
gather_cfg.cout
cout
Definition: gather_cfg.py:144

◆ DumpPosFileSet()

void LaserAlignment::DumpPosFileSet ( LASGlobalData< LASCoordinateSet > &  coordinates)
private

for debugging only, will disappear

debug only, will disappear

Definition at line 1441 of file LaserAlignment.cc.

References EcalCondDBWriter_cfi::beam, gather_cfg::cout, LASCoordinateSet::GetPhi(), LASCoordinateSet::GetPhiError(), LASGlobalData< T >::GetTEC2TECEntry(), LASGlobalData< T >::GetTECEntry(), LASGlobalData< T >::GetTIBTOBEntry(), heppy_loop::loop, and relativeConstraints::ring.

1441  {
1442  LASGlobalLoop loop;
1443  int det, ring, beam, disk, pos;
1444 
1445  std::cout << std::endl << " [LaserAlignment::DumpPosFileSet] -- Dump: " << std::endl;
1446 
1447  // TEC INTERNAL
1448  det = 0;
1449  ring = 0;
1450  beam = 0;
1451  disk = 0;
1452  do {
1453  std::cout << "POS " << det << "\t" << beam << "\t" << disk << "\t" << ring << "\t"
1454  << coordinates.GetTECEntry(det, ring, beam, disk).GetPhi() << "\t"
1455  << coordinates.GetTECEntry(det, ring, beam, disk).GetPhiError() << std::endl;
1456  } while (loop.TECLoop(det, ring, beam, disk));
1457 
1458  // TIBTOB
1459  det = 2;
1460  beam = 0;
1461  pos = 0;
1462  do {
1463  std::cout << "POS " << det << "\t" << beam << "\t" << pos << "\t"
1464  << "-1"
1465  << "\t" << coordinates.GetTIBTOBEntry(det, beam, pos).GetPhi() << "\t"
1466  << coordinates.GetTIBTOBEntry(det, beam, pos).GetPhiError() << std::endl;
1467  } while (loop.TIBTOBLoop(det, beam, pos));
1468 
1469  // TEC2TEC
1470  det = 0;
1471  beam = 0;
1472  disk = 0;
1473  do {
1474  std::cout << "POS " << det << "\t" << beam << "\t" << disk << "\t"
1475  << "-1"
1476  << "\t" << coordinates.GetTEC2TECEntry(det, beam, disk).GetPhi() << "\t"
1477  << coordinates.GetTEC2TECEntry(det, beam, disk).GetPhiError() << std::endl;
1478  } while (loop.TEC2TECLoop(det, beam, disk));
1479 
1480  std::cout << std::endl << " [LaserAlignment::DumpPosFileSet] -- End dump: " << std::endl;
1481 }
heppy_loop.loop
loop
Definition: heppy_loop.py:28
LASCoordinateSet::GetPhiError
double GetPhiError(void) const
Definition: LASCoordinateSet.h:19
EcalCondDBWriter_cfi.beam
beam
Definition: EcalCondDBWriter_cfi.py:45
LASGlobalData::GetTIBTOBEntry
T & GetTIBTOBEntry(int subdetector, int beam, int tibTobPosition)
Definition: LASGlobalData.h:107
LASCoordinateSet::GetPhi
double GetPhi(void) const
Definition: LASCoordinateSet.h:18
LASGlobalData::GetTEC2TECEntry
T & GetTEC2TECEntry(int subdetector, int beam, int tecDisk)
Definition: LASGlobalData.h:129
LASGlobalData::GetTECEntry
T & GetTECEntry(int subdetector, int tecRing, int beam, int tecDisk)
Definition: LASGlobalData.h:84
gather_cfg.cout
cout
Definition: gather_cfg.py:144

◆ DumpStripFileSet()

void LaserAlignment::DumpStripFileSet ( LASGlobalData< std::pair< float, float > > &  measuredStripPositions)
private

for debugging only, will disappear

Definition at line 1486 of file LaserAlignment.cc.

References EcalCondDBWriter_cfi::beam, gather_cfg::cout, heppy_loop::loop, and relativeConstraints::ring.

1486  {
1487  LASGlobalLoop loop;
1488  int det, ring, beam, disk, pos;
1489 
1490  std::cout << std::endl << " [LaserAlignment::DumpStripFileSet] -- Dump: " << std::endl;
1491 
1492  // TEC INTERNAL
1493  det = 0;
1494  ring = 0;
1495  beam = 0;
1496  disk = 0;
1497  do {
1498  std::cout << "STRIP " << det << "\t" << beam << "\t" << disk << "\t" << ring << "\t"
1499  << measuredStripPositions.GetTECEntry(det, ring, beam, disk).first << "\t"
1500  << measuredStripPositions.GetTECEntry(det, ring, beam, disk).second << std::endl;
1501  } while (loop.TECLoop(det, ring, beam, disk));
1502 
1503  // TIBTOB
1504  det = 2;
1505  beam = 0;
1506  pos = 0;
1507  do {
1508  std::cout << "STRIP " << det << "\t" << beam << "\t" << pos << "\t"
1509  << "-1"
1510  << "\t" << measuredStripPositions.GetTIBTOBEntry(det, beam, pos).first << "\t"
1511  << measuredStripPositions.GetTIBTOBEntry(det, beam, pos).second << std::endl;
1512  } while (loop.TIBTOBLoop(det, beam, pos));
1513 
1514  // TEC2TEC
1515  det = 0;
1516  beam = 0;
1517  disk = 0;
1518  do {
1519  std::cout << "STRIP " << det << "\t" << beam << "\t" << disk << "\t"
1520  << "-1"
1521  << "\t" << measuredStripPositions.GetTEC2TECEntry(det, beam, disk).first << "\t"
1522  << measuredStripPositions.GetTEC2TECEntry(det, beam, disk).second << std::endl;
1523  } while (loop.TEC2TECLoop(det, beam, disk));
1524 
1525  std::cout << std::endl << " [LaserAlignment::DumpStripFileSet] -- End dump: " << std::endl;
1526 }
heppy_loop.loop
loop
Definition: heppy_loop.py:28
EcalCondDBWriter_cfi.beam
beam
Definition: EcalCondDBWriter_cfi.py:45
LASGlobalData::GetTIBTOBEntry
T & GetTIBTOBEntry(int subdetector, int beam, int tibTobPosition)
Definition: LASGlobalData.h:107
LASGlobalData::GetTEC2TECEntry
T & GetTEC2TECEntry(int subdetector, int beam, int tecDisk)
Definition: LASGlobalData.h:129
LASGlobalData::GetTECEntry
T & GetTECEntry(int subdetector, int tecRing, int beam, int tecDisk)
Definition: LASGlobalData.h:84
gather_cfg.cout
cout
Definition: gather_cfg.py:144

◆ endJob()

void LaserAlignment::endJob ( )
overridevirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 935 of file LaserAlignment.cc.

Referenced by o2olib.O2ORunMgr::executeJob().

935 {}

◆ endRunProduce()

void LaserAlignment::endRunProduce ( edm::Run &  theRun,
const edm::EventSetup &  theSetup 
)
override

laser hit section for trackbased interface

due to the peculiar order of beams in TkLasBeamCollection, we cannot use the LASGlobalLoop object here

Definition at line 489 of file LaserAlignment.cc.

References AlignableTracker::alignmentErrors(), AlignableTracker::alignments(), ApplyATMaskingCorrections(), LASGeometryUpdater::ApplyBeamKinkCorrections(), ApplyEndcapMaskingCorrections(), PV3DBase< T, PVType, FrameType >::barePhi(), EcalCondDBWriter_cfi::beam, cond::service::PoolDBOutputService::beginOfTime(), CalculateNominalCoordinates(), LASAlignmentTubeAlgorithm::CalculateParameters(), LASBarrelAlgorithm::CalculateParameters(), LASEndcapAlgorithm::CalculateParameters(), collectedDataProfiles, ConvertAngle(), gather_cfg::cout, detectorId, DumpHitmaps(), LASGeometryUpdater::EndcapUpdate(), LASPeakFinder::FindPeakIn(), LASGlobalData< T >::GetTEC2TECEntry(), getTEC2TECNominalBeamOffset(), LASGlobalData< T >::GetTECEntry(), LASGlobalData< T >::GetTIBTOBEntry(), getTIBTOBNominalBeamOffset(), TrackerGeometry::idToDet(), edm::Service< T >::isAvailable(), align_cfg::iteration, misalignedByRefGeometry, moduleLoop, eostools::move(), nominalCoordinates, numberOfAcceptedProfiles, peakFinderThreshold, LASBarrelAlignmentParameterSet::Print(), LASEndcapAlignmentParameterSet::Print(), TkLasBeam::push_back(), edm::Run::put(), relativeConstraints::ring, LASPeakFinder::SetAmplitudeThreshold(), LASGeometryUpdater::SetMisalignmentFromRefGeometry(), LASCoordinateSet::SetPhi(), LASCoordinateSet::SetPhiError(), LASGeometryUpdater::SetReverseDirection(), LASGlobalData< T >::SetTEC2TECEntry(), LASGlobalData< T >::SetTECEntry(), LASGlobalData< T >::SetTIBTOBEntry(), summedHistograms, LASGlobalLoop::TEC2TECLoop(), LASGlobalLoop::TECLoop(), theAlignableTracker, theAlignRecordName, theApplyBeamKinkCorrections, theErrorRecordName, theEvents, theLasConstants, theMaskAtModules, theMaskTecModules, theSetNominalStrips, theStoreToDB, theTrackerGeometry, theUseMinuitAlgorithm, LASGlobalLoop::TIBTOBLoop(), LASGeometryUpdater::TrackerUpdate(), updateFromInputGeometry, and cond::service::PoolDBOutputService::writeOneIOV().

489  {
490  std::cout << " [LaserAlignment::endRun] -- Total number of events processed: " << theEvents << std::endl;
491 
492  // for debugging only..
493  DumpHitmaps(numberOfAcceptedProfiles);
494 
495  // index variables for the LASGlobalLoop objects
496  int det, ring, beam, disk, pos;
497 
498  // measured positions container for the algorithms
499  LASGlobalData<LASCoordinateSet> measuredCoordinates;
500 
501  // fitted peak positions in units of strips (pair for value,error)
502  LASGlobalData<std::pair<float, float> > measuredStripPositions;
503 
504  // the peak finder, a pair (pos/posErr in units of strips) for its results, and the success confirmation
505  LASPeakFinder peakFinder;
506  peakFinder.SetAmplitudeThreshold(peakFinderThreshold);
507  std::pair<double, double> peakFinderResults;
508  bool isGoodFit;
509 
510  // tracker geom. object for calculating the global beam positions
511  const TrackerGeometry& theTracker(*theTrackerGeometry);
512 
513  // fill LASGlobalData<LASCoordinateSet> nominalCoordinates
514  CalculateNominalCoordinates();
515 
516  // for determining the phi errors
517  // ErrorFrameTransformer errorTransformer; // later...
518 
519  // do the fits for TEC+- internal
520  det = 0;
521  ring = 0;
522  beam = 0;
523  disk = 0;
524  do {
525  // do the fit
526  isGoodFit = peakFinder.FindPeakIn(collectedDataProfiles.GetTECEntry(det, ring, beam, disk),
527  peakFinderResults,
528  summedHistograms.GetTECEntry(det, ring, beam, disk),
529  0); // offset is 0 for TEC
530 
531  // now we have the measured positions in units of strips.
532  if (!isGoodFit)
533  std::cout << " [LaserAlignment::endRun] ** WARNING: Fit failed for TEC det: " << det << ", ring: " << ring
534  << ", beam: " << beam << ", disk: " << disk << " (id: " << detectorId.GetTECEntry(det, ring, beam, disk)
535  << ")." << std::endl;
536 
537  // <- here we will later implement the kink corrections
538 
539  // access the tracker geometry for this module
540  const DetId theDetId(detectorId.GetTECEntry(det, ring, beam, disk));
541  const StripGeomDetUnit* const theStripDet = dynamic_cast<const StripGeomDetUnit*>(theTracker.idToDet(theDetId));
542 
543  if (theStripDet) {
544  // first, set the measured coordinates to their nominal values
545  measuredCoordinates.SetTECEntry(det, ring, beam, disk, nominalCoordinates.GetTECEntry(det, ring, beam, disk));
546 
547  if (isGoodFit) { // convert strip position to global phi and replace the nominal phi value/error
548 
549  measuredStripPositions.GetTECEntry(det, ring, beam, disk) = peakFinderResults;
550  const float positionInStrips =
551  theSetNominalStrips
552  ? 256.
553  : peakFinderResults.first; // implementation of "ForceFitterToNominalStrips" config parameter
554  const GlobalPoint& globalPoint =
555  theStripDet->surface().toGlobal(theStripDet->specificTopology().localPosition(positionInStrips));
556  measuredCoordinates.GetTECEntry(det, ring, beam, disk).SetPhi(ConvertAngle(globalPoint.barePhi()));
557 
558  // const GlobalError& globalError = errorTransformer.transform( theStripDet->specificTopology().localError( peakFinderResults.first, pow( peakFinderResults.second, 2 ) ), theStripDet->surface() );
559  // measuredCoordinates.GetTECEntry( det, ring, beam, disk ).SetPhiError( globalError.phierr( globalPoint ) );
560  measuredCoordinates.GetTECEntry(det, ring, beam, disk).SetPhiError(0.00046); // PRELIMINARY ESTIMATE
561 
562  } else { // keep nominal position (middle-of-module) but set a giant phi error so that the module can be ignored by the alignment algorithm
563  measuredStripPositions.GetTECEntry(det, ring, beam, disk) = std::pair<float, float>(256., 1000.);
564  const GlobalPoint& globalPoint =
565  theStripDet->surface().toGlobal(theStripDet->specificTopology().localPosition(256.));
566  measuredCoordinates.GetTECEntry(det, ring, beam, disk).SetPhi(ConvertAngle(globalPoint.barePhi()));
567  measuredCoordinates.GetTECEntry(det, ring, beam, disk).SetPhiError(1000.);
568  }
569  }
570 
571  } while (moduleLoop.TECLoop(det, ring, beam, disk));
572 
573  // do the fits for TIB/TOB
574  det = 2;
575  beam = 0;
576  pos = 0;
577  do {
578  // do the fit
579  isGoodFit = peakFinder.FindPeakIn(collectedDataProfiles.GetTIBTOBEntry(det, beam, pos),
580  peakFinderResults,
581  summedHistograms.GetTIBTOBEntry(det, beam, pos),
582  getTIBTOBNominalBeamOffset(det, beam, pos));
583 
584  // now we have the measured positions in units of strips.
585  if (!isGoodFit)
586  std::cout << " [LaserAlignment::endJob] ** WARNING: Fit failed for TIB/TOB det: " << det << ", beam: " << beam
587  << ", pos: " << pos << " (id: " << detectorId.GetTIBTOBEntry(det, beam, pos) << ")." << std::endl;
588 
589  // <- here we will later implement the kink corrections
590 
591  // access the tracker geometry for this module
592  const DetId theDetId(detectorId.GetTIBTOBEntry(det, beam, pos));
593  const StripGeomDetUnit* const theStripDet = dynamic_cast<const StripGeomDetUnit*>(theTracker.idToDet(theDetId));
594 
595  if (theStripDet) {
596  // first, set the measured coordinates to their nominal values
597  measuredCoordinates.SetTIBTOBEntry(det, beam, pos, nominalCoordinates.GetTIBTOBEntry(det, beam, pos));
598 
599  if (isGoodFit) { // convert strip position to global phi and replace the nominal phi value/error
600  measuredStripPositions.GetTIBTOBEntry(det, beam, pos) = peakFinderResults;
601  const float positionInStrips =
602  theSetNominalStrips
603  ? 256. + getTIBTOBNominalBeamOffset(det, beam, pos)
604  : peakFinderResults.first; // implementation of "ForceFitterToNominalStrips" config parameter
605  const GlobalPoint& globalPoint =
606  theStripDet->surface().toGlobal(theStripDet->specificTopology().localPosition(positionInStrips));
607  measuredCoordinates.GetTIBTOBEntry(det, beam, pos).SetPhi(ConvertAngle(globalPoint.barePhi()));
608  measuredCoordinates.GetTIBTOBEntry(det, beam, pos).SetPhiError(0.00028); // PRELIMINARY ESTIMATE
609  } else { // keep nominal position but set a giant phi error so that the module can be ignored by the alignment algorithm
610  measuredStripPositions.GetTIBTOBEntry(det, beam, pos) =
611  std::pair<float, float>(256. + getTIBTOBNominalBeamOffset(det, beam, pos), 1000.);
612  const GlobalPoint& globalPoint = theStripDet->surface().toGlobal(
613  theStripDet->specificTopology().localPosition(256. + getTIBTOBNominalBeamOffset(det, beam, pos)));
614  measuredCoordinates.GetTIBTOBEntry(det, beam, pos).SetPhi(ConvertAngle(globalPoint.barePhi()));
615  measuredCoordinates.GetTIBTOBEntry(det, beam, pos).SetPhiError(1000.);
616  }
617  }
618 
619  } while (moduleLoop.TIBTOBLoop(det, beam, pos));
620 
621  // do the fits for TEC AT
622  det = 0;
623  beam = 0;
624  disk = 0;
625  do {
626  // do the fit
627  isGoodFit = peakFinder.FindPeakIn(collectedDataProfiles.GetTEC2TECEntry(det, beam, disk),
628  peakFinderResults,
629  summedHistograms.GetTEC2TECEntry(det, beam, disk),
630  getTEC2TECNominalBeamOffset(det, beam, disk));
631  // now we have the positions in units of strips.
632  if (!isGoodFit)
633  std::cout << " [LaserAlignment::endRun] ** WARNING: Fit failed for TEC2TEC det: " << det << ", beam: " << beam
634  << ", disk: " << disk << " (id: " << detectorId.GetTEC2TECEntry(det, beam, disk) << ")." << std::endl;
635 
636  // <- here we will later implement the kink corrections
637 
638  // access the tracker geometry for this module
639  const DetId theDetId(detectorId.GetTEC2TECEntry(det, beam, disk));
640  const StripGeomDetUnit* const theStripDet = dynamic_cast<const StripGeomDetUnit*>(theTracker.idToDet(theDetId));
641 
642  if (theStripDet) {
643  // first, set the measured coordinates to their nominal values
644  measuredCoordinates.SetTEC2TECEntry(det, beam, disk, nominalCoordinates.GetTEC2TECEntry(det, beam, disk));
645 
646  if (isGoodFit) { // convert strip position to global phi and replace the nominal phi value/error
647  measuredStripPositions.GetTEC2TECEntry(det, beam, disk) = peakFinderResults;
648  const float positionInStrips =
649  theSetNominalStrips
650  ? 256. + getTEC2TECNominalBeamOffset(det, beam, disk)
651  : peakFinderResults.first; // implementation of "ForceFitterToNominalStrips" config parameter
652  const GlobalPoint& globalPoint =
653  theStripDet->surface().toGlobal(theStripDet->specificTopology().localPosition(positionInStrips));
654  measuredCoordinates.GetTEC2TECEntry(det, beam, disk).SetPhi(ConvertAngle(globalPoint.barePhi()));
655  measuredCoordinates.GetTEC2TECEntry(det, beam, disk).SetPhiError(0.00047); // PRELIMINARY ESTIMATE
656  } else { // keep nominal position but set a giant phi error so that the module can be ignored by the alignment algorithm
657  measuredStripPositions.GetTEC2TECEntry(det, beam, disk) =
658  std::pair<float, float>(256. + getTEC2TECNominalBeamOffset(det, beam, disk), 1000.);
659  const GlobalPoint& globalPoint = theStripDet->surface().toGlobal(
660  theStripDet->specificTopology().localPosition(256. + getTEC2TECNominalBeamOffset(det, beam, disk)));
661  measuredCoordinates.GetTEC2TECEntry(det, beam, disk).SetPhi(ConvertAngle(globalPoint.barePhi()));
662  measuredCoordinates.GetTEC2TECEntry(det, beam, disk).SetPhiError(1000.);
663  }
664  }
665 
666  } while (moduleLoop.TEC2TECLoop(det, beam, disk));
667 
668  // see what we got (for debugging)
669  // DumpStripFileSet( measuredStripPositions );
670  // DumpPosFileSet( measuredCoordinates );
671 
672  // CALCULATE PARAMETERS AND UPDATE DB OBJECT
673  // for beam kink corrections, reconstructing the geometry and updating the db object
674  LASGeometryUpdater geometryUpdater(nominalCoordinates, theLasConstants);
675 
676  // apply all beam corrections
677  if (theApplyBeamKinkCorrections)
678  geometryUpdater.ApplyBeamKinkCorrections(measuredCoordinates);
679 
680  // if we start with input geometry instead of IDEAL,
681  // reverse the adjustments in the AlignableTracker object
682  if (updateFromInputGeometry)
683  geometryUpdater.SetReverseDirection(true);
684 
685  // if we have "virtual" misalignment which is introduced via the reference geometry,
686  // tell the LASGeometryUpdater to reverse x & y adjustments
687  if (misalignedByRefGeometry)
688  geometryUpdater.SetMisalignmentFromRefGeometry(true);
689 
690  // run the endcap algorithm
691  LASEndcapAlgorithm endcapAlgorithm;
692  LASEndcapAlignmentParameterSet endcapParameters;
693 
694  // this basically sets all the endcap modules to be masked
695  // to their nominal positions (since endcapParameters is overall zero)
696  if (!theMaskTecModules.empty()) {
697  ApplyEndcapMaskingCorrections(measuredCoordinates, nominalCoordinates, endcapParameters);
698  }
699 
700  // run the algorithm
701  endcapParameters = endcapAlgorithm.CalculateParameters(measuredCoordinates, nominalCoordinates);
702 
703  //
704  // loop to mask out events
705  // DESCRIPTION:
706  //
707 
708  // do this only if there are modules to be masked..
709  if (!theMaskTecModules.empty()) {
710  const unsigned int nIterations = 30;
711  for (unsigned int iteration = 0; iteration < nIterations; ++iteration) {
712  // set the endcap modules to be masked to their positions
713  // according to the reconstructed parameters
714  ApplyEndcapMaskingCorrections(measuredCoordinates, nominalCoordinates, endcapParameters);
715 
716  // modifications applied, so re-run the algorithm
717  endcapParameters = endcapAlgorithm.CalculateParameters(measuredCoordinates, nominalCoordinates);
718  }
719  }
720 
721  // these are now final, so:
722  endcapParameters.Print();
723 
724  // do a pre-alignment of the endcaps (TEC2TEC only)
725  // so that the alignment tube algorithms finds orderly disks
726  geometryUpdater.EndcapUpdate(endcapParameters, measuredCoordinates);
727 
728  // the alignment tube algorithms, choose from config
729  LASBarrelAlignmentParameterSet alignmentTubeParameters;
730  // the MINUIT-BASED alignment tube algorithm
731  LASBarrelAlgorithm barrelAlgorithm;
732  // the ANALYTICAL alignment tube algorithm
733  LASAlignmentTubeAlgorithm alignmentTubeAlgorithm;
734 
735  // this basically sets all the modules to be masked
736  // to their nominal positions (since alignmentTubeParameters is overall zero)
737  if (!theMaskAtModules.empty()) {
738  ApplyATMaskingCorrections(measuredCoordinates, nominalCoordinates, alignmentTubeParameters);
739  }
740 
741  if (theUseMinuitAlgorithm) {
742  // run the MINUIT-BASED alignment tube algorithm
743  alignmentTubeParameters = barrelAlgorithm.CalculateParameters(measuredCoordinates, nominalCoordinates);
744  } else {
745  // the ANALYTICAL alignment tube algorithm
746  alignmentTubeParameters = alignmentTubeAlgorithm.CalculateParameters(measuredCoordinates, nominalCoordinates);
747  }
748 
749  //
750  // loop to mask out events
751  // DESCRIPTION:
752  //
753 
754  // do this only if there are modules to be masked..
755  if (!theMaskAtModules.empty()) {
756  const unsigned int nIterations = 30;
757  for (unsigned int iteration = 0; iteration < nIterations; ++iteration) {
758  // set the AT modules to be masked to their positions
759  // according to the reconstructed parameters
760  ApplyATMaskingCorrections(measuredCoordinates, nominalCoordinates, alignmentTubeParameters);
761 
762  // modifications applied, so re-run the algorithm
763  if (theUseMinuitAlgorithm) {
764  alignmentTubeParameters = barrelAlgorithm.CalculateParameters(measuredCoordinates, nominalCoordinates);
765  } else {
766  alignmentTubeParameters = alignmentTubeAlgorithm.CalculateParameters(measuredCoordinates, nominalCoordinates);
767  }
768  }
769  }
770 
771  // these are now final, so:
772  alignmentTubeParameters.Print();
773 
774  // combine the results and update the db object
775  geometryUpdater.TrackerUpdate(endcapParameters, alignmentTubeParameters, *theAlignableTracker);
776 
781 
782  // the collection container
783  auto laserBeams = std::make_unique<TkLasBeamCollection>();
784 
785  // first for the endcap internal beams
786  for (det = 0; det < 2; ++det) {
787  for (ring = 0; ring < 2; ++ring) {
788  for (beam = 0; beam < 8; ++beam) {
789  // the beam and its identifier (see TkLasTrackBasedInterface TWiki)
790  TkLasBeam currentBeam(100 * det + 10 * beam + ring);
791 
792  // order the hits in the beam by increasing z
793  const int firstDisk = det == 0 ? 0 : 8;
794  const int lastDisk = det == 0 ? 8 : 0;
795 
796  // count upwards or downwards
797  for (disk = firstDisk; det == 0 ? disk <= lastDisk : disk >= lastDisk; det == 0 ? ++disk : --disk) {
798  // detId for the SiStripLaserRecHit2D
799  const SiStripDetId theDetId(detectorId.GetTECEntry(det, ring, beam, disk));
800 
801  // need this to calculate the localPosition and its error
802  const StripGeomDetUnit* const theStripDet =
803  dynamic_cast<const StripGeomDetUnit*>(theTracker.idToDet(theDetId));
804 
805  // the hit container
806  const SiStripLaserRecHit2D currentHit(theStripDet->specificTopology().localPosition(
807  measuredStripPositions.GetTECEntry(det, ring, beam, disk).first),
808  theStripDet->specificTopology().localError(
809  measuredStripPositions.GetTECEntry(det, ring, beam, disk).first,
810  measuredStripPositions.GetTECEntry(det, ring, beam, disk).second),
811  theDetId);
812 
813  currentBeam.push_back(currentHit);
814  }
815 
816  laserBeams->push_back(currentBeam);
817  }
818  }
819  }
820 
821  // then, following the convention in TkLasTrackBasedInterface TWiki, the alignment tube beams;
822  // they comprise hits in TIBTOB & TEC2TEC
823 
824  for (beam = 0; beam < 8; ++beam) {
825  // the beam and its identifier (see TkLasTrackBasedInterface TWiki)
826  TkLasBeam currentBeam(100 * 2 /*beamGroup=AT=2*/ + 10 * beam + 0 /*ring=0*/);
827 
828  // first: tec-
829  det = 1;
830  for (disk = 4; disk >= 0; --disk) {
831  // detId for the SiStripLaserRecHit2D
832  const SiStripDetId theDetId(detectorId.GetTEC2TECEntry(det, beam, disk));
833 
834  // need this to calculate the localPosition and its error
835  const StripGeomDetUnit* const theStripDet = dynamic_cast<const StripGeomDetUnit*>(theTracker.idToDet(theDetId));
836 
837  // the hit container
838  const SiStripLaserRecHit2D currentHit(
839  theStripDet->specificTopology().localPosition(measuredStripPositions.GetTEC2TECEntry(det, beam, disk).first),
840  theStripDet->specificTopology().localError(measuredStripPositions.GetTEC2TECEntry(det, beam, disk).first,
841  measuredStripPositions.GetTEC2TECEntry(det, beam, disk).second),
842  theDetId);
843 
844  currentBeam.push_back(currentHit);
845  }
846 
847  // now TIB and TOB in one go
848  for (det = 2; det < 4; ++det) {
849  for (pos = 5; pos >= 0; --pos) { // stupidly, pos is defined from +z to -z in LASGlobalLoop
850 
851  // detId for the SiStripLaserRecHit2D
852  const SiStripDetId theDetId(detectorId.GetTIBTOBEntry(det, beam, pos));
853 
854  // need this to calculate the localPosition and its error
855  const StripGeomDetUnit* const theStripDet = dynamic_cast<const StripGeomDetUnit*>(theTracker.idToDet(theDetId));
856 
857  // the hit container
858  const SiStripLaserRecHit2D currentHit(
859  theStripDet->specificTopology().localPosition(measuredStripPositions.GetTIBTOBEntry(det, beam, pos).first),
860  theStripDet->specificTopology().localError(measuredStripPositions.GetTIBTOBEntry(det, beam, pos).first,
861  measuredStripPositions.GetTIBTOBEntry(det, beam, pos).second),
862  theDetId);
863 
864  currentBeam.push_back(currentHit);
865  }
866  }
867 
868  // then: tec+
869  det = 0;
870  for (disk = 0; disk < 5; ++disk) {
871  // detId for the SiStripLaserRecHit2D
872  const SiStripDetId theDetId(detectorId.GetTEC2TECEntry(det, beam, disk));
873 
874  // need this to calculate the localPosition and its error
875  const StripGeomDetUnit* const theStripDet = dynamic_cast<const StripGeomDetUnit*>(theTracker.idToDet(theDetId));
876 
877  // the hit container
878  const SiStripLaserRecHit2D currentHit(
879  theStripDet->specificTopology().localPosition(measuredStripPositions.GetTEC2TECEntry(det, beam, disk).first),
880  theStripDet->specificTopology().localError(measuredStripPositions.GetTEC2TECEntry(det, beam, disk).first,
881  measuredStripPositions.GetTEC2TECEntry(det, beam, disk).second),
882  theDetId);
883 
884  currentBeam.push_back(currentHit);
885  }
886 
887  // save this beam to the beamCollection
888  laserBeams->push_back(currentBeam);
889 
890  } // (close beam loop)
891 
892  // now attach the collection to the run
893  theRun.put(std::move(laserBeams), "tkLaserBeams");
894 
895  // store the estimated alignment parameters into the DB
896  // first get them
897  Alignments alignments = *(theAlignableTracker->alignments());
898  AlignmentErrorsExtended alignmentErrors = *(theAlignableTracker->alignmentErrors());
899 
900  if (theStoreToDB) {
901  std::cout << " [LaserAlignment::endRun] -- Storing the calculated alignment parameters to the DataBase:"
902  << std::endl;
903 
904  // Call service
905  edm::Service<cond::service::PoolDBOutputService> poolDbService;
906  if (!poolDbService.isAvailable()) // Die if not available
907  throw cms::Exception("NotAvailable") << "PoolDBOutputService not available";
908 
909  // Store
910 
911  // if ( poolDbService->isNewTagRequest(theAlignRecordName) ) {
912  // poolDbService->createOneIOV<Alignments>( alignments, poolDbService->currentTime(), theAlignRecordName );
913  // }
914  // else {
915  // poolDbService->appendOneIOV<Alignments>( alignments, poolDbService->currentTime(), theAlignRecordName );
916  // }
917  poolDbService->writeOneIOV<Alignments>(alignments, poolDbService->beginOfTime(), theAlignRecordName);
918 
919  // if ( poolDbService->isNewTagRequest(theErrorRecordName) ) {
920  // poolDbService->createOneIOV<AlignmentErrorsExtended>( alignmentErrors, poolDbService->currentTime(), poolDbService->endOfTime(), theErrorRecordName );
921  // }
922  // else {
923  // poolDbService->appendOneIOV<AlignmentErrorsExtended>( alignmentErrors, poolDbService->currentTime(), theErrorRecordName );
924  // }
925  poolDbService->writeOneIOV<AlignmentErrorsExtended>(
926  alignmentErrors, poolDbService->beginOfTime(), theErrorRecordName);
927 
928  std::cout << " [LaserAlignment::endRun] -- Storing done." << std::endl;
929  }
930 }
LASGlobalData::SetTEC2TECEntry
void SetTEC2TECEntry(int subdetector, int beam, int tecDisk, T)
Definition: LASGlobalData.h:195
LASCoordinateSet::SetPhi
void SetPhi(double aPhi)
Definition: LASCoordinateSet.h:27
align_cfg.iteration
iteration
Definition: align_cfg.py:5
LASCoordinateSet::SetPhiError
void SetPhiError(double aPhiError)
Definition: LASCoordinateSet.h:32
LaserAlignment::ApplyEndcapMaskingCorrections
void ApplyEndcapMaskingCorrections(LASGlobalData< LASCoordinateSet > &, LASGlobalData< LASCoordinateSet > &, LASEndcapAlignmentParameterSet &)
apply endcap correction to masked modules in TEC
Definition: LaserAlignment.cc:1585
AlignableTracker::alignments
Alignments * alignments() const override
Return alignments, sorted by DetId.
Definition: AlignableTracker.cc:42
LaserAlignment::misalignedByRefGeometry
bool misalignedByRefGeometry
config switch
Definition: LaserAlignment.h:164
LaserAlignment::theAlignRecordName
std::string theAlignRecordName
Definition: LaserAlignment.h:242
LaserAlignment::getTIBTOBNominalBeamOffset
double getTIBTOBNominalBeamOffset(unsigned int, unsigned int, unsigned int)
returns the nominal beam position (strips) in TOB for the profileJudge
Definition: LaserAlignment.cc:1279
LaserAlignment::theTrackerGeometry
edm::ESHandle< TrackerGeometry > theTrackerGeometry
Definition: LaserAlignment.h:237
LaserAlignment::peakFinderThreshold
double peakFinderThreshold
config parameter
Definition: LaserAlignment.h:152
LASGlobalLoop::TEC2TECLoop
bool TEC2TECLoop(int &, int &, int &) const
Definition: LASGlobalLoop.cc:74
LaserAlignment::DumpHitmaps
void DumpHitmaps(LASGlobalData< int > &)
for debugging only, will disappear
Definition: LaserAlignment.cc:1531
LaserAlignment::numberOfAcceptedProfiles
LASGlobalData< int > numberOfAcceptedProfiles
Definition: LaserAlignment.h:214
LaserAlignment::theMaskTecModules
std::vector< unsigned int > theMaskTecModules
config parameters
Definition: LaserAlignment.h:182
LaserAlignment::theUseMinuitAlgorithm
bool theUseMinuitAlgorithm
config switch
Definition: LaserAlignment.h:146
EcalCondDBWriter_cfi.beam
beam
Definition: EcalCondDBWriter_cfi.py:45
PV3DBase::barePhi
T barePhi() const
Definition: PV3DBase.h:65
LaserAlignment::summedHistograms
LASGlobalData< TH1D * > summedHistograms
Definition: LaserAlignment.h:222
LASPeakFinder::FindPeakIn
bool FindPeakIn(const LASModuleProfile &, std::pair< double, double > &, TH1D *, const double)
Definition: LASPeakFinder.cc:16
LaserAlignment::ConvertAngle
double ConvertAngle(double)
convert an angle in the [-pi,pi] range to the [0,2*pi] range
Definition: LaserAlignment.cc:1426
LASGlobalLoop::TECLoop
bool TECLoop(int &, int &, int &, int &) const
Definition: LASGlobalLoop.cc:13
LASBarrelAlgorithm::CalculateParameters
LASBarrelAlignmentParameterSet CalculateParameters(LASGlobalData< LASCoordinateSet > &, LASGlobalData< LASCoordinateSet > &)
Definition: LASBarrelAlgorithm.cc:20
cond::service::PoolDBOutputService::writeOneIOV
Hash writeOneIOV(const T &payload, Time_t time, const std::string &recordName)
Definition: PoolDBOutputService.h:62
LASPeakFinder::SetAmplitudeThreshold
void SetAmplitudeThreshold(double)
Definition: LASPeakFinder.cc:106
LaserAlignment::theApplyBeamKinkCorrections
bool theApplyBeamKinkCorrections
config switch
Definition: LaserAlignment.h:149
LASGlobalData::GetTIBTOBEntry
T & GetTIBTOBEntry(int subdetector, int beam, int tibTobPosition)
Definition: LASGlobalData.h:107
LASGlobalData::SetTIBTOBEntry
void SetTIBTOBEntry(int subdetector, int beam, int tibTobPosition, T)
Definition: LASGlobalData.h:173
LaserAlignment::CalculateNominalCoordinates
void CalculateNominalCoordinates(void)
fills a LASGlobalData<LASCoordinateSet> with nominal module positions
Definition: LaserAlignment.cc:1332
LASEndcapAlgorithm::CalculateParameters
LASEndcapAlignmentParameterSet CalculateParameters(LASGlobalData< LASCoordinateSet > &, LASGlobalData< LASCoordinateSet > &)
Definition: LASEndcapAlgorithm.cc:15
LaserAlignment::theErrorRecordName
std::string theErrorRecordName
Definition: LaserAlignment.h:242
LASGlobalLoop::TIBTOBLoop
bool TIBTOBLoop(int &, int &, int &) const
Definition: LASGlobalLoop.cc:46
LaserAlignment::theStoreToDB
bool theStoreToDB
config switch
Definition: LaserAlignment.h:167
SiStripDetId
Detector identifier class for the strip tracker.
Definition: SiStripDetId.h:18
DetId
Definition: DetId.h:17
LASGlobalData::GetTEC2TECEntry
T & GetTEC2TECEntry(int subdetector, int beam, int tecDisk)
Definition: LASGlobalData.h:129
LaserAlignment::nominalCoordinates
LASGlobalData< LASCoordinateSet > nominalCoordinates
Definition: LaserAlignment.h:225
LaserAlignment::theAlignableTracker
AlignableTracker * theAlignableTracker
Definition: LaserAlignment.h:240
LaserAlignment::theEvents
int theEvents
counter for the total number of events processed
Definition: LaserAlignment.h:140
LASAlignmentTubeAlgorithm::CalculateParameters
LASBarrelAlignmentParameterSet CalculateParameters(LASGlobalData< LASCoordinateSet > &, LASGlobalData< LASCoordinateSet > &)
Definition: LASAlignmentTubeAlgorithm.cc:12
LaserAlignment::theLasConstants
LASConstants theLasConstants
Definition: LaserAlignment.h:193
LaserAlignment::getTEC2TECNominalBeamOffset
double getTEC2TECNominalBeamOffset(unsigned int, unsigned int, unsigned int)
returns the nominal beam position (strips) in TEC (AT) for the profileJudge
Definition: LaserAlignment.cc:1314
edm::Run::put
void put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Run.h:109
LASGlobalData::GetTECEntry
T & GetTECEntry(int subdetector, int tecRing, int beam, int tecDisk)
Definition: LASGlobalData.h:84
LaserAlignment::theMaskAtModules
std::vector< unsigned int > theMaskAtModules
Definition: LaserAlignment.h:183
AlignableTracker::alignmentErrors
AlignmentErrorsExtended * alignmentErrors() const override
Return alignment errors, sorted by DetId.
Definition: AlignableTracker.cc:59
LaserAlignment::ApplyATMaskingCorrections
void ApplyATMaskingCorrections(LASGlobalData< LASCoordinateSet > &, LASGlobalData< LASCoordinateSet > &, LASBarrelAlignmentParameterSet &)
same for alignment tube modules
Definition: LaserAlignment.cc:1625
LaserAlignment::collectedDataProfiles
LASGlobalData< LASModuleProfile > collectedDataProfiles
Definition: LaserAlignment.h:208
LaserAlignment::detectorId
LASGlobalData< unsigned int > detectorId
Definition: LaserAlignment.h:196
edm::Service::isAvailable
bool isAvailable() const
Definition: Service.h:40
gather_cfg.cout
cout
Definition: gather_cfg.py:144
LASGlobalData::SetTECEntry
void SetTECEntry(int subdetector, int tecRing, int beam, int tecDisk, T)
Definition: LASGlobalData.h:150
eostools.move
def move(src, dest)
Definition: eostools.py:511
LaserAlignment::moduleLoop
LASGlobalLoop moduleLoop
Definition: LaserAlignment.h:229
LaserAlignment::updateFromInputGeometry
bool updateFromInputGeometry
config switch
Definition: LaserAlignment.h:161
LaserAlignment::theSetNominalStrips
bool theSetNominalStrips
config switch
Definition: LaserAlignment.h:186

◆ fillDataProfiles()

void LaserAlignment::fillDataProfiles ( edm::Event const &  theEvent,
edm::EventSetup const &  theSetup 
)
private

fill profiles from SiStrip(Raw)Digi container

fills the module profiles (LASGlobalLoop<LASModuleProfile> currentDataProfiles) from the event digi containers, distinguishing between SiStripDigi or SiStripRawDigi.

Definition at line 941 of file LaserAlignment.cc.

References SiStripDigi::adc(), SiStripRawDigi::adc(), EcalCondDBWriter_cfi::beam, edm::DetSetVector< T >::begin(), currentDataProfiles, detectorId, DataMixerDataOnData_cff::digiLabel, EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::digiProducer, HLT_2022v12_cff::distance, edm::DetSetVector< T >::end(), Exception, edm::Event::getByLabel(), LASGlobalData< T >::GetTEC2TECEntry(), LASGlobalData< T >::GetTECEntry(), LASGlobalData< T >::GetTIBTOBEntry(), moduleLoop, relativeConstraints::ring, LASModuleProfile::SetAllValuesTo(), LASModuleProfile::SetValue(), AlCaHLTBitMon_QueryRunRegistry::string, SiStripDigi::strip(), LASGlobalLoop::TEC2TECLoop(), LASGlobalLoop::TECLoop(), theDigiProducersList, and LASGlobalLoop::TIBTOBLoop().

Referenced by produce().

941  {
942  // two handles for the two different kinds of digis
943  edm::Handle<edm::DetSetVector<SiStripRawDigi> > theStripRawDigis;
944  edm::Handle<edm::DetSetVector<SiStripDigi> > theStripDigis;
945 
946  bool isRawDigi = false;
947 
948  // indices for the LASGlobalLoop object
949  int det = 0, ring = 0, beam = 0, disk = 0, pos = 0;
950 
951  // query config set and loop over all PSets in the VPSet
952  for (std::vector<edm::ParameterSet>::iterator itDigiProducersList = theDigiProducersList.begin();
953  itDigiProducersList != theDigiProducersList.end();
954  ++itDigiProducersList) {
955  std::string digiProducer = itDigiProducersList->getParameter<std::string>("DigiProducer");
956  std::string digiLabel = itDigiProducersList->getParameter<std::string>("DigiLabel");
957  std::string digiType = itDigiProducersList->getParameter<std::string>("DigiType");
958 
959  // now branch according to digi type (raw or processed);
960  // first we go for raw digis => SiStripRawDigi
961  if (digiType == "Raw") {
962  theEvent.getByLabel(digiProducer, digiLabel, theStripRawDigis);
963  isRawDigi = true;
964  } else if (digiType == "Processed") {
965  theEvent.getByLabel(digiProducer, digiLabel, theStripDigis);
966  isRawDigi = false;
967  } else {
968  throw cms::Exception(" [LaserAlignment::fillDataProfiles]")
969  << " ** ERROR: Invalid digi type: \"" << digiType << "\" specified in configuration." << std::endl;
970  }
971 
972  // loop TEC internal modules
973  det = 0;
974  ring = 0;
975  beam = 0;
976  disk = 0;
977  do {
978  // first clear the profile
979  currentDataProfiles.GetTECEntry(det, ring, beam, disk).SetAllValuesTo(0.);
980 
981  // retrieve the raw id of that module
982  const int detRawId = detectorId.GetTECEntry(det, ring, beam, disk);
983 
984  if (isRawDigi) { // we have raw SiStripRawDigis
985 
986  // search the digis for the raw id
987  edm::DetSetVector<SiStripRawDigi>::const_iterator detSetIter = theStripRawDigis->find(detRawId);
988  if (detSetIter == theStripRawDigis->end()) {
989  throw cms::Exception("[Laser Alignment::fillDataProfiles]")
990  << " ** ERROR: No raw DetSet found for det: " << detRawId << "." << std::endl;
991  }
992 
993  // fill the digis to the profiles
994  edm::DetSet<SiStripRawDigi>::const_iterator digiRangeIterator = detSetIter->data.begin(); // for the loop
995  edm::DetSet<SiStripRawDigi>::const_iterator digiRangeStart = digiRangeIterator; // save starting positions
996 
997  // loop all digis
998  for (; digiRangeIterator != detSetIter->data.end(); ++digiRangeIterator) {
999  const SiStripRawDigi& digi = *digiRangeIterator;
1000  const int channel = distance(digiRangeStart, digiRangeIterator);
1001  if (channel >= 0 && channel < 512)
1002  currentDataProfiles.GetTECEntry(det, ring, beam, disk).SetValue(channel, digi.adc());
1003  else
1004  throw cms::Exception("[Laser Alignment::fillDataProfiles]")
1005  << " ** ERROR: raw digi channel: " << channel << " out of range for det: " << detRawId << "."
1006  << std::endl;
1007  }
1008 
1009  }
1010 
1011  else { // we have zero suppressed SiStripDigis
1012 
1013  // search the digis for the raw id
1014  edm::DetSetVector<SiStripDigi>::const_iterator detSetIter = theStripDigis->find(detRawId);
1015 
1016  // processed DetSets may be missing, just skip
1017  if (detSetIter == theStripDigis->end())
1018  continue;
1019 
1020  // fill the digis to the profiles
1021  edm::DetSet<SiStripDigi>::const_iterator digiRangeIterator = detSetIter->data.begin(); // for the loop
1022 
1023  for (; digiRangeIterator != detSetIter->data.end(); ++digiRangeIterator) {
1024  const SiStripDigi& digi = *digiRangeIterator;
1025  if (digi.strip() < 512)
1026  currentDataProfiles.GetTECEntry(det, ring, beam, disk).SetValue(digi.strip(), digi.adc());
1027  else
1028  throw cms::Exception("[Laser Alignment::fillDataProfiles]")
1029  << " ** ERROR: digi strip: " << digi.strip() << " out of range for det: " << detRawId << "."
1030  << std::endl;
1031  }
1032  }
1033 
1034  } while (moduleLoop.TECLoop(det, ring, beam, disk));
1035 
1036  // loop TIBTOB modules
1037  det = 2;
1038  beam = 0;
1039  pos = 0;
1040  do {
1041  // first clear the profile
1042  currentDataProfiles.GetTIBTOBEntry(det, beam, pos).SetAllValuesTo(0.);
1043 
1044  // retrieve the raw id of that module
1045  const int detRawId = detectorId.GetTIBTOBEntry(det, beam, pos);
1046 
1047  if (isRawDigi) { // we have raw SiStripRawDigis
1048 
1049  // search the digis for the raw id
1050  edm::DetSetVector<SiStripRawDigi>::const_iterator detSetIter = theStripRawDigis->find(detRawId);
1051  if (detSetIter == theStripRawDigis->end()) {
1052  throw cms::Exception("[Laser Alignment::fillDataProfiles]")
1053  << " ** ERROR: No raw DetSet found for det: " << detRawId << "." << std::endl;
1054  }
1055 
1056  // fill the digis to the profiles
1057  edm::DetSet<SiStripRawDigi>::const_iterator digiRangeIterator = detSetIter->data.begin(); // for the loop
1058  edm::DetSet<SiStripRawDigi>::const_iterator digiRangeStart = digiRangeIterator; // save starting positions
1059 
1060  // loop all digis
1061  for (; digiRangeIterator != detSetIter->data.end(); ++digiRangeIterator) {
1062  const SiStripRawDigi& digi = *digiRangeIterator;
1063  const int channel = distance(digiRangeStart, digiRangeIterator);
1064  if (channel >= 0 && channel < 512)
1065  currentDataProfiles.GetTIBTOBEntry(det, beam, pos).SetValue(channel, digi.adc());
1066  else
1067  throw cms::Exception("[Laser Alignment::fillDataProfiles]")
1068  << " ** ERROR: raw digi channel: " << channel << " out of range for det: " << detRawId << "."
1069  << std::endl;
1070  }
1071 
1072  }
1073 
1074  else { // we have zero suppressed SiStripDigis
1075 
1076  // search the digis for the raw id
1077  edm::DetSetVector<SiStripDigi>::const_iterator detSetIter = theStripDigis->find(detRawId);
1078 
1079  // processed DetSets may be missing, just skip
1080  if (detSetIter == theStripDigis->end())
1081  continue;
1082 
1083  // fill the digis to the profiles
1084  edm::DetSet<SiStripDigi>::const_iterator digiRangeIterator = detSetIter->data.begin(); // for the loop
1085 
1086  for (; digiRangeIterator != detSetIter->data.end(); ++digiRangeIterator) {
1087  const SiStripDigi& digi = *digiRangeIterator;
1088  if (digi.strip() < 512)
1089  currentDataProfiles.GetTIBTOBEntry(det, beam, pos).SetValue(digi.strip(), digi.adc());
1090  else
1091  throw cms::Exception("[Laser Alignment::fillDataProfiles]")
1092  << " ** ERROR: digi strip: " << digi.strip() << " out of range for det: " << detRawId << "."
1093  << std::endl;
1094  }
1095  }
1096 
1097  } while (moduleLoop.TIBTOBLoop(det, beam, pos));
1098 
1099  // loop TEC AT modules
1100  det = 0;
1101  beam = 0;
1102  disk = 0;
1103  do {
1104  // first clear the profile
1105  currentDataProfiles.GetTEC2TECEntry(det, beam, disk).SetAllValuesTo(0.);
1106 
1107  // retrieve the raw id of that module
1108  const int detRawId = detectorId.GetTEC2TECEntry(det, beam, disk);
1109 
1110  if (isRawDigi) { // we have raw SiStripRawDigis
1111 
1112  // search the digis for the raw id
1113  edm::DetSetVector<SiStripRawDigi>::const_iterator detSetIter = theStripRawDigis->find(detRawId);
1114  if (detSetIter == theStripRawDigis->end()) {
1115  throw cms::Exception("[Laser Alignment::fillDataProfiles]")
1116  << " ** ERROR: No raw DetSet found for det: " << detRawId << "." << std::endl;
1117  }
1118 
1119  // fill the digis to the profiles
1120  edm::DetSet<SiStripRawDigi>::const_iterator digiRangeIterator = detSetIter->data.begin(); // for the loop
1121  edm::DetSet<SiStripRawDigi>::const_iterator digiRangeStart = digiRangeIterator; // save starting positions
1122 
1123  // loop all digis
1124  for (; digiRangeIterator != detSetIter->data.end(); ++digiRangeIterator) {
1125  const SiStripRawDigi& digi = *digiRangeIterator;
1126  const int channel = distance(digiRangeStart, digiRangeIterator);
1127  if (channel >= 0 && channel < 512)
1128  currentDataProfiles.GetTEC2TECEntry(det, beam, disk).SetValue(channel, digi.adc());
1129  else
1130  throw cms::Exception("[Laser Alignment::fillDataProfiles]")
1131  << " ** ERROR: raw digi channel: " << channel << " out of range for det: " << detRawId << "."
1132  << std::endl;
1133  }
1134 
1135  }
1136 
1137  else { // we have zero suppressed SiStripDigis
1138 
1139  // search the digis for the raw id
1140  edm::DetSetVector<SiStripDigi>::const_iterator detSetIter = theStripDigis->find(detRawId);
1141 
1142  // processed DetSets may be missing, just skip
1143  if (detSetIter == theStripDigis->end())
1144  continue;
1145 
1146  // fill the digis to the profiles
1147  edm::DetSet<SiStripDigi>::const_iterator digiRangeIterator = detSetIter->data.begin(); // for the loop
1148 
1149  for (; digiRangeIterator != detSetIter->data.end(); ++digiRangeIterator) {
1150  const SiStripDigi& digi = *digiRangeIterator;
1151  if (digi.strip() < 512)
1152  currentDataProfiles.GetTEC2TECEntry(det, beam, disk).SetValue(digi.strip(), digi.adc());
1153  else
1154  throw cms::Exception("[Laser Alignment::fillDataProfiles]")
1155  << " ** ERROR: digi strip: " << digi.strip() << " out of range for det: " << detRawId << "."
1156  << std::endl;
1157  }
1158  }
1159 
1160  } while (moduleLoop.TEC2TECLoop(det, beam, disk));
1161 
1162  } // theDigiProducersList loop
1163 }
LASGlobalLoop::TEC2TECLoop
bool TEC2TECLoop(int &, int &, int &) const
Definition: LASGlobalLoop.cc:74
LASModuleProfile::SetAllValuesTo
void SetAllValuesTo(const double &)
Definition: LASModuleProfile.cc:76
LaserAlignment::currentDataProfiles
LASGlobalData< LASModuleProfile > currentDataProfiles
data profiles for the current event
Definition: LaserAlignment.h:205
SiStripRawDigi::adc
uint16_t adc() const
Definition: SiStripRawDigi.h:22
DataMixerDataOnData_cff.digiLabel
digiLabel
Definition: DataMixerDataOnData_cff.py:46
EcalCondDBWriter_cfi.beam
beam
Definition: EcalCondDBWriter_cfi.py:45
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
LaserAlignment::theDigiProducersList
std::vector< edm::ParameterSet > theDigiProducersList
Definition: LaserAlignment.h:170
SiStripDigi::strip
const uint16_t & strip() const
Definition: SiStripDigi.h:33
SiStripDigi::adc
const uint16_t & adc() const
Definition: SiStripDigi.h:34
EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0.digiProducer
digiProducer
Definition: EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0.py:50
LASGlobalLoop::TECLoop
bool TECLoop(int &, int &, int &, int &) const
Definition: LASGlobalLoop.cc:13
LASGlobalData::GetTIBTOBEntry
T & GetTIBTOBEntry(int subdetector, int beam, int tibTobPosition)
Definition: LASGlobalData.h:107
SiStripDigi
A Digi for the silicon strip detector, containing both strip and adc information, and suitable for st...
Definition: SiStripDigi.h:12
edm::DetSetVector::end
iterator end()
Return the off-the-end iterator.
Definition: DetSetVector.h:325
LASGlobalLoop::TIBTOBLoop
bool TIBTOBLoop(int &, int &, int &) const
Definition: LASGlobalLoop.cc:46
LASModuleProfile::SetValue
void SetValue(unsigned int theStripNumber, const double &theValue)
Definition: LASModuleProfile.h:20
LASGlobalData::GetTEC2TECEntry
T & GetTEC2TECEntry(int subdetector, int beam, int tecDisk)
Definition: LASGlobalData.h:129
LASGlobalData::GetTECEntry
T & GetTECEntry(int subdetector, int tecRing, int beam, int tecDisk)
Definition: LASGlobalData.h:84
LaserAlignment::detectorId
LASGlobalData< unsigned int > detectorId
Definition: LaserAlignment.h:196
edm::DetSetVector::begin
iterator begin()
Return an iterator to the first DetSet.
Definition: DetSetVector.h:314
edm::DetSet::const_iterator
collection_type::const_iterator const_iterator
Definition: DetSet.h:31
edm::DetSetVector::const_iterator
collection_type::const_iterator const_iterator
Definition: DetSetVector.h:102
SiStripRawDigi
A Digi for the silicon strip detector, containing only adc information, and suitable for storing raw ...
Definition: SiStripRawDigi.h:15
HLT_2022v12_cff.distance
distance
Definition: HLT_2022v12_cff.py:5871
LaserAlignment::moduleLoop
LASGlobalLoop moduleLoop
Definition: LaserAlignment.h:229

◆ fillDetectorId()

void LaserAlignment::fillDetectorId ( void  )
private

fill hard coded detIds

all the detector ids for the LAS modules hard-coded

Definition at line 7 of file LaserAlignmentFillDetId.cc.

References detectorId, LASGlobalData< T >::SetTEC2TECEntry(), LASGlobalData< T >::SetTECEntry(), LASGlobalData< T >::SetTIBTOBEntry(), and tecDoubleHitDetId.

Referenced by beginJob().

7  {
8  // these are the detids of the TEC modules hit
9  // by the AT as well as the TEC beams
10  tecDoubleHitDetId.push_back(470307208);
11  tecDoubleHitDetId.push_back(470323592);
12  tecDoubleHitDetId.push_back(470339976);
13  tecDoubleHitDetId.push_back(470356360);
14  tecDoubleHitDetId.push_back(470372744);
15  tecDoubleHitDetId.push_back(470307976);
16  tecDoubleHitDetId.push_back(470324360);
17  tecDoubleHitDetId.push_back(470340744);
18  tecDoubleHitDetId.push_back(470357128);
19  tecDoubleHitDetId.push_back(470373512);
20  tecDoubleHitDetId.push_back(470308488);
21  tecDoubleHitDetId.push_back(470324872);
22  tecDoubleHitDetId.push_back(470341256);
23  tecDoubleHitDetId.push_back(470357640);
24  tecDoubleHitDetId.push_back(470374024);
25  tecDoubleHitDetId.push_back(470045064);
26  tecDoubleHitDetId.push_back(470061448);
27  tecDoubleHitDetId.push_back(470077832);
28  tecDoubleHitDetId.push_back(470094216);
29  tecDoubleHitDetId.push_back(470110600);
30  tecDoubleHitDetId.push_back(470045832);
31  tecDoubleHitDetId.push_back(470062216);
32  tecDoubleHitDetId.push_back(470078600);
33  tecDoubleHitDetId.push_back(470094984);
34  tecDoubleHitDetId.push_back(470111368);
35  tecDoubleHitDetId.push_back(470046344);
36  tecDoubleHitDetId.push_back(470062728);
37  tecDoubleHitDetId.push_back(470079112);
38  tecDoubleHitDetId.push_back(470095496);
39  tecDoubleHitDetId.push_back(470111880);
40 
41  // now all the modules (above included)
42 
43  // TEC+
44  detectorId.SetTECEntry(0, 0, 0, 0, 470307208);
45  detectorId.SetTECEntry(0, 0, 0, 1, 470323592);
46  detectorId.SetTECEntry(0, 0, 0, 2, 470339976);
47  detectorId.SetTECEntry(0, 0, 0, 3, 470356360);
48  detectorId.SetTECEntry(0, 0, 0, 4, 470372744);
49  detectorId.SetTECEntry(0, 0, 0, 5, 470389128);
50  detectorId.SetTECEntry(0, 0, 0, 6, 470405512);
51  detectorId.SetTECEntry(0, 0, 0, 7, 470421896);
52  detectorId.SetTECEntry(0, 0, 0, 8, 470438280);
53  detectorId.SetTECEntry(0, 0, 1, 0, 470307464);
54  detectorId.SetTECEntry(0, 0, 1, 1, 470323848);
55  detectorId.SetTECEntry(0, 0, 1, 2, 470340232);
56  detectorId.SetTECEntry(0, 0, 1, 3, 470356616);
57  detectorId.SetTECEntry(0, 0, 1, 4, 470373000);
58  detectorId.SetTECEntry(0, 0, 1, 5, 470389384);
59  detectorId.SetTECEntry(0, 0, 1, 6, 470405768);
60  detectorId.SetTECEntry(0, 0, 1, 7, 470422152);
61  detectorId.SetTECEntry(0, 0, 1, 8, 470438536);
62  detectorId.SetTECEntry(0, 0, 2, 0, 470307720);
63  detectorId.SetTECEntry(0, 0, 2, 1, 470324104);
64  detectorId.SetTECEntry(0, 0, 2, 2, 470340488);
65  detectorId.SetTECEntry(0, 0, 2, 3, 470356872);
66  detectorId.SetTECEntry(0, 0, 2, 4, 470373256);
67  detectorId.SetTECEntry(0, 0, 2, 5, 470389640);
68  detectorId.SetTECEntry(0, 0, 2, 6, 470406024);
69  detectorId.SetTECEntry(0, 0, 2, 7, 470422408);
70  detectorId.SetTECEntry(0, 0, 2, 8, 470438792);
71  detectorId.SetTECEntry(0, 0, 3, 0, 470307976);
72  detectorId.SetTECEntry(0, 0, 3, 1, 470324360);
73  detectorId.SetTECEntry(0, 0, 3, 2, 470340744);
74  detectorId.SetTECEntry(0, 0, 3, 3, 470357128);
75  detectorId.SetTECEntry(0, 0, 3, 4, 470373512);
76  detectorId.SetTECEntry(0, 0, 3, 5, 470389896);
77  detectorId.SetTECEntry(0, 0, 3, 6, 470406280);
78  detectorId.SetTECEntry(0, 0, 3, 7, 470422664);
79  detectorId.SetTECEntry(0, 0, 3, 8, 470439048);
80  detectorId.SetTECEntry(0, 0, 4, 0, 470308232);
81  detectorId.SetTECEntry(0, 0, 4, 1, 470324616);
82  detectorId.SetTECEntry(0, 0, 4, 2, 470341000);
83  detectorId.SetTECEntry(0, 0, 4, 3, 470357384);
84  detectorId.SetTECEntry(0, 0, 4, 4, 470373768);
85  detectorId.SetTECEntry(0, 0, 4, 5, 470390152);
86  detectorId.SetTECEntry(0, 0, 4, 6, 470406536);
87  detectorId.SetTECEntry(0, 0, 4, 7, 470422920);
88  detectorId.SetTECEntry(0, 0, 4, 8, 470439304);
89  detectorId.SetTECEntry(0, 0, 5, 0, 470308488);
90  detectorId.SetTECEntry(0, 0, 5, 1, 470324872);
91  detectorId.SetTECEntry(0, 0, 5, 2, 470341256);
92  detectorId.SetTECEntry(0, 0, 5, 3, 470357640);
93  detectorId.SetTECEntry(0, 0, 5, 4, 470374024);
94  detectorId.SetTECEntry(0, 0, 5, 5, 470390408);
95  detectorId.SetTECEntry(0, 0, 5, 6, 470406792);
96  detectorId.SetTECEntry(0, 0, 5, 7, 470423176);
97  detectorId.SetTECEntry(0, 0, 5, 8, 470439560);
98  detectorId.SetTECEntry(0, 0, 6, 0, 470308744);
99  detectorId.SetTECEntry(0, 0, 6, 1, 470325128);
100  detectorId.SetTECEntry(0, 0, 6, 2, 470341512);
101  detectorId.SetTECEntry(0, 0, 6, 3, 470357896);
102  detectorId.SetTECEntry(0, 0, 6, 4, 470374280);
103  detectorId.SetTECEntry(0, 0, 6, 5, 470390664);
104  detectorId.SetTECEntry(0, 0, 6, 6, 470407048);
105  detectorId.SetTECEntry(0, 0, 6, 7, 470423432);
106  detectorId.SetTECEntry(0, 0, 6, 8, 470439816);
107  detectorId.SetTECEntry(0, 0, 7, 0, 470309000);
108  detectorId.SetTECEntry(0, 0, 7, 1, 470325384);
109  detectorId.SetTECEntry(0, 0, 7, 2, 470341768);
110  detectorId.SetTECEntry(0, 0, 7, 3, 470358152);
111  detectorId.SetTECEntry(0, 0, 7, 4, 470374536);
112  detectorId.SetTECEntry(0, 0, 7, 5, 470390920);
113  detectorId.SetTECEntry(0, 0, 7, 6, 470407304);
114  detectorId.SetTECEntry(0, 0, 7, 7, 470423688);
115  detectorId.SetTECEntry(0, 0, 7, 8, 470440072);
116  detectorId.SetTECEntry(0, 1, 0, 0, 470307272);
117  detectorId.SetTECEntry(0, 1, 0, 1, 470323656);
118  detectorId.SetTECEntry(0, 1, 0, 2, 470340040);
119  detectorId.SetTECEntry(0, 1, 0, 3, 470356424);
120  detectorId.SetTECEntry(0, 1, 0, 4, 470372808);
121  detectorId.SetTECEntry(0, 1, 0, 5, 470389192);
122  detectorId.SetTECEntry(0, 1, 0, 6, 470405576);
123  detectorId.SetTECEntry(0, 1, 0, 7, 470421960);
124  detectorId.SetTECEntry(0, 1, 0, 8, 470438344);
125  detectorId.SetTECEntry(0, 1, 1, 0, 470307528);
126  detectorId.SetTECEntry(0, 1, 1, 1, 470323912);
127  detectorId.SetTECEntry(0, 1, 1, 2, 470340296);
128  detectorId.SetTECEntry(0, 1, 1, 3, 470356680);
129  detectorId.SetTECEntry(0, 1, 1, 4, 470373064);
130  detectorId.SetTECEntry(0, 1, 1, 5, 470389448);
131  detectorId.SetTECEntry(0, 1, 1, 6, 470405832);
132  detectorId.SetTECEntry(0, 1, 1, 7, 470422216);
133  detectorId.SetTECEntry(0, 1, 1, 8, 470438600);
134  detectorId.SetTECEntry(0, 1, 2, 0, 470307784);
135  detectorId.SetTECEntry(0, 1, 2, 1, 470324168);
136  detectorId.SetTECEntry(0, 1, 2, 2, 470340552);
137  detectorId.SetTECEntry(0, 1, 2, 3, 470356936);
138  detectorId.SetTECEntry(0, 1, 2, 4, 470373320);
139  detectorId.SetTECEntry(0, 1, 2, 5, 470389704);
140  detectorId.SetTECEntry(0, 1, 2, 6, 470406088);
141  detectorId.SetTECEntry(0, 1, 2, 7, 470422472);
142  detectorId.SetTECEntry(0, 1, 2, 8, 470438856);
143  detectorId.SetTECEntry(0, 1, 3, 0, 470308040);
144  detectorId.SetTECEntry(0, 1, 3, 1, 470324424);
145  detectorId.SetTECEntry(0, 1, 3, 2, 470340808);
146  detectorId.SetTECEntry(0, 1, 3, 3, 470357192);
147  detectorId.SetTECEntry(0, 1, 3, 4, 470373576);
148  detectorId.SetTECEntry(0, 1, 3, 5, 470389960);
149  detectorId.SetTECEntry(0, 1, 3, 6, 470406344);
150  detectorId.SetTECEntry(0, 1, 3, 7, 470422728);
151  detectorId.SetTECEntry(0, 1, 3, 8, 470439112);
152  detectorId.SetTECEntry(0, 1, 4, 0, 470308296);
153  detectorId.SetTECEntry(0, 1, 4, 1, 470324680);
154  detectorId.SetTECEntry(0, 1, 4, 2, 470341064);
155  detectorId.SetTECEntry(0, 1, 4, 3, 470357448);
156  detectorId.SetTECEntry(0, 1, 4, 4, 470373832);
157  detectorId.SetTECEntry(0, 1, 4, 5, 470390216);
158  detectorId.SetTECEntry(0, 1, 4, 6, 470406600);
159  detectorId.SetTECEntry(0, 1, 4, 7, 470422984);
160  detectorId.SetTECEntry(0, 1, 4, 8, 470439368);
161  detectorId.SetTECEntry(0, 1, 5, 0, 470308552);
162  detectorId.SetTECEntry(0, 1, 5, 1, 470324936);
163  detectorId.SetTECEntry(0, 1, 5, 2, 470341320);
164  detectorId.SetTECEntry(0, 1, 5, 3, 470357704);
165  detectorId.SetTECEntry(0, 1, 5, 4, 470374088);
166  detectorId.SetTECEntry(0, 1, 5, 5, 470390472);
167  detectorId.SetTECEntry(0, 1, 5, 6, 470406856);
168  detectorId.SetTECEntry(0, 1, 5, 7, 470423240);
169  detectorId.SetTECEntry(0, 1, 5, 8, 470439624);
170  detectorId.SetTECEntry(0, 1, 6, 0, 470308808);
171  detectorId.SetTECEntry(0, 1, 6, 1, 470325192);
172  detectorId.SetTECEntry(0, 1, 6, 2, 470341576);
173  detectorId.SetTECEntry(0, 1, 6, 3, 470357960);
174  detectorId.SetTECEntry(0, 1, 6, 4, 470374344);
175  detectorId.SetTECEntry(0, 1, 6, 5, 470390728);
176  detectorId.SetTECEntry(0, 1, 6, 6, 470407112);
177  detectorId.SetTECEntry(0, 1, 6, 7, 470423496);
178  detectorId.SetTECEntry(0, 1, 6, 8, 470439880);
179  detectorId.SetTECEntry(0, 1, 7, 0, 470309064);
180  detectorId.SetTECEntry(0, 1, 7, 1, 470325448);
181  detectorId.SetTECEntry(0, 1, 7, 2, 470341832);
182  detectorId.SetTECEntry(0, 1, 7, 3, 470358216);
183  detectorId.SetTECEntry(0, 1, 7, 4, 470374600);
184  detectorId.SetTECEntry(0, 1, 7, 5, 470390984);
185  detectorId.SetTECEntry(0, 1, 7, 6, 470407368);
186  detectorId.SetTECEntry(0, 1, 7, 7, 470423752);
187  detectorId.SetTECEntry(0, 1, 7, 8, 470440136);
188 
189  // TEC-
190  detectorId.SetTECEntry(1, 0, 0, 0, 470045064);
191  detectorId.SetTECEntry(1, 0, 0, 1, 470061448);
192  detectorId.SetTECEntry(1, 0, 0, 2, 470077832);
193  detectorId.SetTECEntry(1, 0, 0, 3, 470094216);
194  detectorId.SetTECEntry(1, 0, 0, 4, 470110600);
195  detectorId.SetTECEntry(1, 0, 0, 5, 470126984);
196  detectorId.SetTECEntry(1, 0, 0, 6, 470143368);
197  detectorId.SetTECEntry(1, 0, 0, 7, 470159752);
198  detectorId.SetTECEntry(1, 0, 0, 8, 470176136);
199  detectorId.SetTECEntry(1, 0, 1, 0, 470045320);
200  detectorId.SetTECEntry(1, 0, 1, 1, 470061704);
201  detectorId.SetTECEntry(1, 0, 1, 2, 470078088);
202  detectorId.SetTECEntry(1, 0, 1, 3, 470094472);
203  detectorId.SetTECEntry(1, 0, 1, 4, 470110856);
204  detectorId.SetTECEntry(1, 0, 1, 5, 470127240);
205  detectorId.SetTECEntry(1, 0, 1, 6, 470143624);
206  detectorId.SetTECEntry(1, 0, 1, 7, 470160008);
207  detectorId.SetTECEntry(1, 0, 1, 8, 470176392);
208  detectorId.SetTECEntry(1, 0, 2, 0, 470045576);
209  detectorId.SetTECEntry(1, 0, 2, 1, 470061960);
210  detectorId.SetTECEntry(1, 0, 2, 2, 470078344);
211  detectorId.SetTECEntry(1, 0, 2, 3, 470094728);
212  detectorId.SetTECEntry(1, 0, 2, 4, 470111112);
213  detectorId.SetTECEntry(1, 0, 2, 5, 470127496);
214  detectorId.SetTECEntry(1, 0, 2, 6, 470143880);
215  detectorId.SetTECEntry(1, 0, 2, 7, 470160264);
216  detectorId.SetTECEntry(1, 0, 2, 8, 470176648);
217  detectorId.SetTECEntry(1, 0, 3, 0, 470045832);
218  detectorId.SetTECEntry(1, 0, 3, 1, 470062216);
219  detectorId.SetTECEntry(1, 0, 3, 2, 470078600);
220  detectorId.SetTECEntry(1, 0, 3, 3, 470094984);
221  detectorId.SetTECEntry(1, 0, 3, 4, 470111368);
222  detectorId.SetTECEntry(1, 0, 3, 5, 470127752);
223  detectorId.SetTECEntry(1, 0, 3, 6, 470144136);
224  detectorId.SetTECEntry(1, 0, 3, 7, 470160520);
225  detectorId.SetTECEntry(1, 0, 3, 8, 470176904);
226  detectorId.SetTECEntry(1, 0, 4, 0, 470046088);
227  detectorId.SetTECEntry(1, 0, 4, 1, 470062472);
228  detectorId.SetTECEntry(1, 0, 4, 2, 470078856);
229  detectorId.SetTECEntry(1, 0, 4, 3, 470095240);
230  detectorId.SetTECEntry(1, 0, 4, 4, 470111624);
231  detectorId.SetTECEntry(1, 0, 4, 5, 470128008);
232  detectorId.SetTECEntry(1, 0, 4, 6, 470144392);
233  detectorId.SetTECEntry(1, 0, 4, 7, 470160776);
234  detectorId.SetTECEntry(1, 0, 4, 8, 470177160);
235  detectorId.SetTECEntry(1, 0, 5, 0, 470046344);
236  detectorId.SetTECEntry(1, 0, 5, 1, 470062728);
237  detectorId.SetTECEntry(1, 0, 5, 2, 470079112);
238  detectorId.SetTECEntry(1, 0, 5, 3, 470095496);
239  detectorId.SetTECEntry(1, 0, 5, 4, 470111880);
240  detectorId.SetTECEntry(1, 0, 5, 5, 470128264);
241  detectorId.SetTECEntry(1, 0, 5, 6, 470144648);
242  detectorId.SetTECEntry(1, 0, 5, 7, 470161032);
243  detectorId.SetTECEntry(1, 0, 5, 8, 470177416);
244  detectorId.SetTECEntry(1, 0, 6, 0, 470046600);
245  detectorId.SetTECEntry(1, 0, 6, 1, 470062984);
246  detectorId.SetTECEntry(1, 0, 6, 2, 470079368);
247  detectorId.SetTECEntry(1, 0, 6, 3, 470095752);
248  detectorId.SetTECEntry(1, 0, 6, 4, 470112136);
249  detectorId.SetTECEntry(1, 0, 6, 5, 470128520);
250  detectorId.SetTECEntry(1, 0, 6, 6, 470144904);
251  detectorId.SetTECEntry(1, 0, 6, 7, 470161288);
252  detectorId.SetTECEntry(1, 0, 6, 8, 470177672);
253  detectorId.SetTECEntry(1, 0, 7, 0, 470046856);
254  detectorId.SetTECEntry(1, 0, 7, 1, 470063240);
255  detectorId.SetTECEntry(1, 0, 7, 2, 470079624);
256  detectorId.SetTECEntry(1, 0, 7, 3, 470096008);
257  detectorId.SetTECEntry(1, 0, 7, 4, 470112392);
258  detectorId.SetTECEntry(1, 0, 7, 5, 470128776);
259  detectorId.SetTECEntry(1, 0, 7, 6, 470145160);
260  detectorId.SetTECEntry(1, 0, 7, 7, 470161544);
261  detectorId.SetTECEntry(1, 0, 7, 8, 470177928);
262  detectorId.SetTECEntry(1, 1, 0, 0, 470045128);
263  detectorId.SetTECEntry(1, 1, 0, 1, 470061512);
264  detectorId.SetTECEntry(1, 1, 0, 2, 470077896);
265  detectorId.SetTECEntry(1, 1, 0, 3, 470094280);
266  detectorId.SetTECEntry(1, 1, 0, 4, 470110664);
267  detectorId.SetTECEntry(1, 1, 0, 5, 470127048);
268  detectorId.SetTECEntry(1, 1, 0, 6, 470143432);
269  detectorId.SetTECEntry(1, 1, 0, 7, 470159816);
270  detectorId.SetTECEntry(1, 1, 0, 8, 470176200);
271  detectorId.SetTECEntry(1, 1, 1, 0, 470045384);
272  detectorId.SetTECEntry(1, 1, 1, 1, 470061768);
273  detectorId.SetTECEntry(1, 1, 1, 2, 470078152);
274  detectorId.SetTECEntry(1, 1, 1, 3, 470094536);
275  detectorId.SetTECEntry(1, 1, 1, 4, 470110920);
276  detectorId.SetTECEntry(1, 1, 1, 5, 470127304);
277  detectorId.SetTECEntry(1, 1, 1, 6, 470143688);
278  detectorId.SetTECEntry(1, 1, 1, 7, 470160072);
279  detectorId.SetTECEntry(1, 1, 1, 8, 470176456);
280  detectorId.SetTECEntry(1, 1, 2, 0, 470045640);
281  detectorId.SetTECEntry(1, 1, 2, 1, 470062024);
282  detectorId.SetTECEntry(1, 1, 2, 2, 470078408);
283  detectorId.SetTECEntry(1, 1, 2, 3, 470094792);
284  detectorId.SetTECEntry(1, 1, 2, 4, 470111176);
285  detectorId.SetTECEntry(1, 1, 2, 5, 470127560);
286  detectorId.SetTECEntry(1, 1, 2, 6, 470143944);
287  detectorId.SetTECEntry(1, 1, 2, 7, 470160328);
288  detectorId.SetTECEntry(1, 1, 2, 8, 470176712);
289  detectorId.SetTECEntry(1, 1, 3, 0, 470045896);
290  detectorId.SetTECEntry(1, 1, 3, 1, 470062280);
291  detectorId.SetTECEntry(1, 1, 3, 2, 470078664);
292  detectorId.SetTECEntry(1, 1, 3, 3, 470095048);
293  detectorId.SetTECEntry(1, 1, 3, 4, 470111432);
294  detectorId.SetTECEntry(1, 1, 3, 5, 470127816);
295  detectorId.SetTECEntry(1, 1, 3, 6, 470144200);
296  detectorId.SetTECEntry(1, 1, 3, 7, 470160584);
297  detectorId.SetTECEntry(1, 1, 3, 8, 470176968);
298  detectorId.SetTECEntry(1, 1, 4, 0, 470046152);
299  detectorId.SetTECEntry(1, 1, 4, 1, 470062536);
300  detectorId.SetTECEntry(1, 1, 4, 2, 470078920);
301  detectorId.SetTECEntry(1, 1, 4, 3, 470095304);
302  detectorId.SetTECEntry(1, 1, 4, 4, 470111688);
303  detectorId.SetTECEntry(1, 1, 4, 5, 470128072);
304  detectorId.SetTECEntry(1, 1, 4, 6, 470144456);
305  detectorId.SetTECEntry(1, 1, 4, 7, 470160840);
306  detectorId.SetTECEntry(1, 1, 4, 8, 470177224);
307  detectorId.SetTECEntry(1, 1, 5, 0, 470046408);
308  detectorId.SetTECEntry(1, 1, 5, 1, 470062792);
309  detectorId.SetTECEntry(1, 1, 5, 2, 470079176);
310  detectorId.SetTECEntry(1, 1, 5, 3, 470095560);
311  detectorId.SetTECEntry(1, 1, 5, 4, 470111944);
312  detectorId.SetTECEntry(1, 1, 5, 5, 470128328);
313  detectorId.SetTECEntry(1, 1, 5, 6, 470144712);
314  detectorId.SetTECEntry(1, 1, 5, 7, 470161096);
315  detectorId.SetTECEntry(1, 1, 5, 8, 470177480);
316  detectorId.SetTECEntry(1, 1, 6, 0, 470046664);
317  detectorId.SetTECEntry(1, 1, 6, 1, 470063048);
318  detectorId.SetTECEntry(1, 1, 6, 2, 470079432);
319  detectorId.SetTECEntry(1, 1, 6, 3, 470095816);
320  detectorId.SetTECEntry(1, 1, 6, 4, 470112200);
321  detectorId.SetTECEntry(1, 1, 6, 5, 470128584);
322  detectorId.SetTECEntry(1, 1, 6, 6, 470144968);
323  detectorId.SetTECEntry(1, 1, 6, 7, 470161352);
324  detectorId.SetTECEntry(1, 1, 6, 8, 470177736);
325  detectorId.SetTECEntry(1, 1, 7, 0, 470046920);
326  detectorId.SetTECEntry(1, 1, 7, 1, 470063304);
327  detectorId.SetTECEntry(1, 1, 7, 2, 470079688);
328  detectorId.SetTECEntry(1, 1, 7, 3, 470096072);
329  detectorId.SetTECEntry(1, 1, 7, 4, 470112456);
330  detectorId.SetTECEntry(1, 1, 7, 5, 470128840);
331  detectorId.SetTECEntry(1, 1, 7, 6, 470145224);
332  detectorId.SetTECEntry(1, 1, 7, 7, 470161608);
333  detectorId.SetTECEntry(1, 1, 7, 8, 470177992);
334 
335  // TIB
336  detectorId.SetTIBTOBEntry(2, 0, 0, 369174604);
337  detectorId.SetTIBTOBEntry(2, 0, 1, 369174600);
338  detectorId.SetTIBTOBEntry(2, 0, 2, 369174596);
339  detectorId.SetTIBTOBEntry(2, 0, 3, 369170500);
340  detectorId.SetTIBTOBEntry(2, 0, 4, 369170504);
341  detectorId.SetTIBTOBEntry(2, 0, 5, 369170508);
342  detectorId.SetTIBTOBEntry(2, 1, 0, 369174732);
343  detectorId.SetTIBTOBEntry(2, 1, 1, 369174728);
344  detectorId.SetTIBTOBEntry(2, 1, 2, 369174724);
345  detectorId.SetTIBTOBEntry(2, 1, 3, 369170628);
346  detectorId.SetTIBTOBEntry(2, 1, 4, 369170632);
347  detectorId.SetTIBTOBEntry(2, 1, 5, 369170636);
348  detectorId.SetTIBTOBEntry(2, 2, 0, 369174812);
349  detectorId.SetTIBTOBEntry(2, 2, 1, 369174808);
350  detectorId.SetTIBTOBEntry(2, 2, 2, 369174804);
351  detectorId.SetTIBTOBEntry(2, 2, 3, 369170708);
352  detectorId.SetTIBTOBEntry(2, 2, 4, 369170712);
353  detectorId.SetTIBTOBEntry(2, 2, 5, 369170716);
354  detectorId.SetTIBTOBEntry(2, 3, 0, 369174940);
355  detectorId.SetTIBTOBEntry(2, 3, 1, 369174936);
356  detectorId.SetTIBTOBEntry(2, 3, 2, 369174932);
357  detectorId.SetTIBTOBEntry(2, 3, 3, 369170836);
358  detectorId.SetTIBTOBEntry(2, 3, 4, 369170840);
359  detectorId.SetTIBTOBEntry(2, 3, 5, 369170844);
360  detectorId.SetTIBTOBEntry(2, 4, 0, 369175068);
361  detectorId.SetTIBTOBEntry(2, 4, 1, 369175064);
362  detectorId.SetTIBTOBEntry(2, 4, 2, 369175060);
363  detectorId.SetTIBTOBEntry(2, 4, 3, 369170964);
364  detectorId.SetTIBTOBEntry(2, 4, 4, 369170968);
365  detectorId.SetTIBTOBEntry(2, 4, 5, 369170972);
366  detectorId.SetTIBTOBEntry(2, 5, 0, 369175164);
367  detectorId.SetTIBTOBEntry(2, 5, 1, 369175160);
368  detectorId.SetTIBTOBEntry(2, 5, 2, 369175156);
369  detectorId.SetTIBTOBEntry(2, 5, 3, 369171060);
370  detectorId.SetTIBTOBEntry(2, 5, 4, 369171064);
371  detectorId.SetTIBTOBEntry(2, 5, 5, 369171068);
372  detectorId.SetTIBTOBEntry(2, 6, 0, 369175292);
373  detectorId.SetTIBTOBEntry(2, 6, 1, 369175288);
374  detectorId.SetTIBTOBEntry(2, 6, 2, 369175284);
375  detectorId.SetTIBTOBEntry(2, 6, 3, 369171188);
376  detectorId.SetTIBTOBEntry(2, 6, 4, 369171192);
377  detectorId.SetTIBTOBEntry(2, 6, 5, 369171196);
378  detectorId.SetTIBTOBEntry(2, 7, 0, 369175372);
379  detectorId.SetTIBTOBEntry(2, 7, 1, 369175368);
380  detectorId.SetTIBTOBEntry(2, 7, 2, 369175364);
381  detectorId.SetTIBTOBEntry(2, 7, 3, 369171268);
382  detectorId.SetTIBTOBEntry(2, 7, 4, 369171272);
383  detectorId.SetTIBTOBEntry(2, 7, 5, 369171276);
384 
385  // TOB
386  detectorId.SetTIBTOBEntry(3, 0, 0, 436232314);
387  detectorId.SetTIBTOBEntry(3, 0, 1, 436232306);
388  detectorId.SetTIBTOBEntry(3, 0, 2, 436232298);
389  detectorId.SetTIBTOBEntry(3, 0, 3, 436228198);
390  detectorId.SetTIBTOBEntry(3, 0, 4, 436228206);
391  detectorId.SetTIBTOBEntry(3, 0, 5, 436228214);
392  detectorId.SetTIBTOBEntry(3, 1, 0, 436232506);
393  detectorId.SetTIBTOBEntry(3, 1, 1, 436232498);
394  detectorId.SetTIBTOBEntry(3, 1, 2, 436232490);
395  detectorId.SetTIBTOBEntry(3, 1, 3, 436228390);
396  detectorId.SetTIBTOBEntry(3, 1, 4, 436228398);
397  detectorId.SetTIBTOBEntry(3, 1, 5, 436228406);
398  detectorId.SetTIBTOBEntry(3, 2, 0, 436232634);
399  detectorId.SetTIBTOBEntry(3, 2, 1, 436232626);
400  detectorId.SetTIBTOBEntry(3, 2, 2, 436232618);
401  detectorId.SetTIBTOBEntry(3, 2, 3, 436228518);
402  detectorId.SetTIBTOBEntry(3, 2, 4, 436228526);
403  detectorId.SetTIBTOBEntry(3, 2, 5, 436228534);
404  detectorId.SetTIBTOBEntry(3, 3, 0, 436232826);
405  detectorId.SetTIBTOBEntry(3, 3, 1, 436232818);
406  detectorId.SetTIBTOBEntry(3, 3, 2, 436232810);
407  detectorId.SetTIBTOBEntry(3, 3, 3, 436228710);
408  detectorId.SetTIBTOBEntry(3, 3, 4, 436228718);
409  detectorId.SetTIBTOBEntry(3, 3, 5, 436228726);
410  detectorId.SetTIBTOBEntry(3, 4, 0, 436233018);
411  detectorId.SetTIBTOBEntry(3, 4, 1, 436233010);
412  detectorId.SetTIBTOBEntry(3, 4, 2, 436233002);
413  detectorId.SetTIBTOBEntry(3, 4, 3, 436228902);
414  detectorId.SetTIBTOBEntry(3, 4, 4, 436228910);
415  detectorId.SetTIBTOBEntry(3, 4, 5, 436228918);
416  detectorId.SetTIBTOBEntry(3, 5, 0, 436233146);
417  detectorId.SetTIBTOBEntry(3, 5, 1, 436233138);
418  detectorId.SetTIBTOBEntry(3, 5, 2, 436233130);
419  detectorId.SetTIBTOBEntry(3, 5, 3, 436229030);
420  detectorId.SetTIBTOBEntry(3, 5, 4, 436229038);
421  detectorId.SetTIBTOBEntry(3, 5, 5, 436229046);
422  detectorId.SetTIBTOBEntry(3, 6, 0, 436233338);
423  detectorId.SetTIBTOBEntry(3, 6, 1, 436233330);
424  detectorId.SetTIBTOBEntry(3, 6, 2, 436233322);
425  detectorId.SetTIBTOBEntry(3, 6, 3, 436229222);
426  detectorId.SetTIBTOBEntry(3, 6, 4, 436229230);
427  detectorId.SetTIBTOBEntry(3, 6, 5, 436229238);
428  detectorId.SetTIBTOBEntry(3, 7, 0, 436233466);
429  detectorId.SetTIBTOBEntry(3, 7, 1, 436233458);
430  detectorId.SetTIBTOBEntry(3, 7, 2, 436233450);
431  detectorId.SetTIBTOBEntry(3, 7, 3, 436229350);
432  detectorId.SetTIBTOBEntry(3, 7, 4, 436229358);
433  detectorId.SetTIBTOBEntry(3, 7, 5, 436229366);
434 
435  // TEC+ AT
436  detectorId.SetTEC2TECEntry(0, 0, 0, 470307208);
437  detectorId.SetTEC2TECEntry(0, 0, 1, 470323592);
438  detectorId.SetTEC2TECEntry(0, 0, 2, 470339976);
439  detectorId.SetTEC2TECEntry(0, 0, 3, 470356360);
440  detectorId.SetTEC2TECEntry(0, 0, 4, 470372744);
441  detectorId.SetTEC2TECEntry(0, 1, 0, 470307468);
442  detectorId.SetTEC2TECEntry(0, 1, 1, 470323852);
443  detectorId.SetTEC2TECEntry(0, 1, 2, 470340236);
444  detectorId.SetTEC2TECEntry(0, 1, 3, 470356620);
445  detectorId.SetTEC2TECEntry(0, 1, 4, 470373004);
446  detectorId.SetTEC2TECEntry(0, 2, 0, 470307716);
447  detectorId.SetTEC2TECEntry(0, 2, 1, 470324100);
448  detectorId.SetTEC2TECEntry(0, 2, 2, 470340484);
449  detectorId.SetTEC2TECEntry(0, 2, 3, 470356868);
450  detectorId.SetTEC2TECEntry(0, 2, 4, 470373252);
451  detectorId.SetTEC2TECEntry(0, 3, 0, 470307976);
452  detectorId.SetTEC2TECEntry(0, 3, 1, 470324360);
453  detectorId.SetTEC2TECEntry(0, 3, 2, 470340744);
454  detectorId.SetTEC2TECEntry(0, 3, 3, 470357128);
455  detectorId.SetTEC2TECEntry(0, 3, 4, 470373512);
456  detectorId.SetTEC2TECEntry(0, 4, 0, 470308236);
457  detectorId.SetTEC2TECEntry(0, 4, 1, 470324620);
458  detectorId.SetTEC2TECEntry(0, 4, 2, 470341004);
459  detectorId.SetTEC2TECEntry(0, 4, 3, 470357388);
460  detectorId.SetTEC2TECEntry(0, 4, 4, 470373772);
461  detectorId.SetTEC2TECEntry(0, 5, 0, 470308488);
462  detectorId.SetTEC2TECEntry(0, 5, 1, 470324872);
463  detectorId.SetTEC2TECEntry(0, 5, 2, 470341256);
464  detectorId.SetTEC2TECEntry(0, 5, 3, 470357640);
465  detectorId.SetTEC2TECEntry(0, 5, 4, 470374024);
466  detectorId.SetTEC2TECEntry(0, 6, 0, 470308748);
467  detectorId.SetTEC2TECEntry(0, 6, 1, 470325132);
468  detectorId.SetTEC2TECEntry(0, 6, 2, 470341516);
469  detectorId.SetTEC2TECEntry(0, 6, 3, 470357900);
470  detectorId.SetTEC2TECEntry(0, 6, 4, 470374284);
471  detectorId.SetTEC2TECEntry(0, 7, 0, 470308996);
472  detectorId.SetTEC2TECEntry(0, 7, 1, 470325380);
473  detectorId.SetTEC2TECEntry(0, 7, 2, 470341764);
474  detectorId.SetTEC2TECEntry(0, 7, 3, 470358148);
475  detectorId.SetTEC2TECEntry(0, 7, 4, 470374532);
476 
477  // TEC- AT
478  detectorId.SetTEC2TECEntry(1, 0, 0, 470045064);
479  detectorId.SetTEC2TECEntry(1, 0, 1, 470061448);
480  detectorId.SetTEC2TECEntry(1, 0, 2, 470077832);
481  detectorId.SetTEC2TECEntry(1, 0, 3, 470094216);
482  detectorId.SetTEC2TECEntry(1, 0, 4, 470110600);
483  detectorId.SetTEC2TECEntry(1, 1, 0, 470045316);
484  detectorId.SetTEC2TECEntry(1, 1, 1, 470061700);
485  detectorId.SetTEC2TECEntry(1, 1, 2, 470078084);
486  detectorId.SetTEC2TECEntry(1, 1, 3, 470094468);
487  detectorId.SetTEC2TECEntry(1, 1, 4, 470110852);
488  detectorId.SetTEC2TECEntry(1, 2, 0, 470045580);
489  detectorId.SetTEC2TECEntry(1, 2, 1, 470061964);
490  detectorId.SetTEC2TECEntry(1, 2, 2, 470078348);
491  detectorId.SetTEC2TECEntry(1, 2, 3, 470094732);
492  detectorId.SetTEC2TECEntry(1, 2, 4, 470111116);
493  detectorId.SetTEC2TECEntry(1, 3, 0, 470045832);
494  detectorId.SetTEC2TECEntry(1, 3, 1, 470062216);
495  detectorId.SetTEC2TECEntry(1, 3, 2, 470078600);
496  detectorId.SetTEC2TECEntry(1, 3, 3, 470094984);
497  detectorId.SetTEC2TECEntry(1, 3, 4, 470111368);
498  detectorId.SetTEC2TECEntry(1, 4, 0, 470046084);
499  detectorId.SetTEC2TECEntry(1, 4, 1, 470062468);
500  detectorId.SetTEC2TECEntry(1, 4, 2, 470078852);
501  detectorId.SetTEC2TECEntry(1, 4, 3, 470095236);
502  detectorId.SetTEC2TECEntry(1, 4, 4, 470111620);
503  detectorId.SetTEC2TECEntry(1, 5, 0, 470046344);
504  detectorId.SetTEC2TECEntry(1, 5, 1, 470062728);
505  detectorId.SetTEC2TECEntry(1, 5, 2, 470079112);
506  detectorId.SetTEC2TECEntry(1, 5, 3, 470095496);
507  detectorId.SetTEC2TECEntry(1, 5, 4, 470111880);
508  detectorId.SetTEC2TECEntry(1, 6, 0, 470046596);
509  detectorId.SetTEC2TECEntry(1, 6, 1, 470062980);
510  detectorId.SetTEC2TECEntry(1, 6, 2, 470079364);
511  detectorId.SetTEC2TECEntry(1, 6, 3, 470095748);
512  detectorId.SetTEC2TECEntry(1, 6, 4, 470112132);
513  detectorId.SetTEC2TECEntry(1, 7, 0, 470046860);
514  detectorId.SetTEC2TECEntry(1, 7, 1, 470063244);
515  detectorId.SetTEC2TECEntry(1, 7, 2, 470079628);
516  detectorId.SetTEC2TECEntry(1, 7, 3, 470096012);
517  detectorId.SetTEC2TECEntry(1, 7, 4, 470112396);
518 }
LASGlobalData::SetTEC2TECEntry
void SetTEC2TECEntry(int subdetector, int beam, int tecDisk, T)
Definition: LASGlobalData.h:195
LaserAlignment::tecDoubleHitDetId
std::vector< unsigned int > tecDoubleHitDetId
Definition: LaserAlignment.h:199
LASGlobalData::SetTIBTOBEntry
void SetTIBTOBEntry(int subdetector, int beam, int tibTobPosition, T)
Definition: LASGlobalData.h:173
LaserAlignment::detectorId
LASGlobalData< unsigned int > detectorId
Definition: LaserAlignment.h:196
LASGlobalData::SetTECEntry
void SetTECEntry(int subdetector, int tecRing, int beam, int tecDisk, T)
Definition: LASGlobalData.h:150

◆ fillPedestalProfiles()

void LaserAlignment::fillPedestalProfiles ( edm::ESHandle< SiStripPedestals > &  pedestalsHandle)
private

fill pedestals from dbase

This function fills the pedestal profiles (LASGlobalData<LASModuleProfiles> pedestalProfiles) from the ESHandle (from file or DB)

Argument: readily connected SiStripPedestals object (get() alredy called) The functionality inside the loops is basically taken from: CommonTools/SiStripZeroSuppression/src/SiStripPedestalsSubtractor.cc

Definition at line 1173 of file LaserAlignment.cc.

References EcalCondDBWriter_cfi::beam, detectorId, SiStripPedestals::getPed(), SiStripPedestals::getRange(), LASGlobalData< T >::GetTEC2TECEntry(), LASGlobalData< T >::GetTECEntry(), LASGlobalData< T >::GetTIBTOBEntry(), createfilelist::int, moduleLoop, pedestalProfiles, relativeConstraints::ring, LASModuleProfile::SetValue(), digitizers_cfi::strip, LASGlobalLoop::TEC2TECLoop(), LASGlobalLoop::TECLoop(), and LASGlobalLoop::TIBTOBLoop().

Referenced by produce().

1173  {
1174  int det, ring, beam, disk, pos;
1175 
1176  // loop TEC modules (yet without AT)
1177  det = 0;
1178  ring = 0;
1179  beam = 0;
1180  disk = 0;
1181  do { // loop using LASGlobalLoop functionality
1182  SiStripPedestals::Range pedRange = pedestalsHandle->getRange(detectorId.GetTECEntry(det, ring, beam, disk));
1183  for (int strip = 0; strip < 512; ++strip) {
1184  int thePedestal = int(pedestalsHandle->getPed(strip, pedRange));
1185  if (thePedestal > 895)
1186  thePedestal -= 1024;
1187  pedestalProfiles.GetTECEntry(det, ring, beam, disk).SetValue(strip, thePedestal);
1188  }
1189  } while (moduleLoop.TECLoop(det, ring, beam, disk));
1190 
1191  // TIB & TOB section
1192  det = 2;
1193  beam = 0;
1194  pos = 0;
1195  do { // loop using LASGlobalLoop functionality
1196  SiStripPedestals::Range pedRange = pedestalsHandle->getRange(detectorId.GetTIBTOBEntry(det, beam, pos));
1197  for (int strip = 0; strip < 512; ++strip) {
1198  int thePedestal = int(pedestalsHandle->getPed(strip, pedRange));
1199  if (thePedestal > 895)
1200  thePedestal -= 1024;
1201  pedestalProfiles.GetTIBTOBEntry(det, beam, pos).SetValue(strip, thePedestal);
1202  }
1203  } while (moduleLoop.TIBTOBLoop(det, beam, pos));
1204 
1205  // TEC2TEC AT section
1206  det = 0;
1207  beam = 0;
1208  disk = 0;
1209  do { // loop using LASGlobalLoop functionality
1210  SiStripPedestals::Range pedRange = pedestalsHandle->getRange(detectorId.GetTEC2TECEntry(det, beam, disk));
1211  for (int strip = 0; strip < 512; ++strip) {
1212  int thePedestal = int(pedestalsHandle->getPed(strip, pedRange));
1213  if (thePedestal > 895)
1214  thePedestal -= 1024;
1215  pedestalProfiles.GetTEC2TECEntry(det, beam, disk).SetValue(strip, thePedestal);
1216  }
1217  } while (moduleLoop.TEC2TECLoop(det, beam, disk));
1218 }
SiStripPedestals::getRange
const Range getRange(const uint32_t &detID) const
Definition: SiStripPedestals.cc:28
LASGlobalLoop::TEC2TECLoop
bool TEC2TECLoop(int &, int &, int &) const
Definition: LASGlobalLoop.cc:74
SiStripPedestals::Range
std::pair< ContainerIterator, ContainerIterator > Range
Definition: SiStripPedestals.h:51
EcalCondDBWriter_cfi.beam
beam
Definition: EcalCondDBWriter_cfi.py:45
SiStripPedestals::getPed
float getPed(const uint16_t &strip, const Range &range) const
Definition: SiStripPedestals.cc:51
LASGlobalLoop::TECLoop
bool TECLoop(int &, int &, int &, int &) const
Definition: LASGlobalLoop.cc:13
LASGlobalData::GetTIBTOBEntry
T & GetTIBTOBEntry(int subdetector, int beam, int tibTobPosition)
Definition: LASGlobalData.h:107
LASGlobalLoop::TIBTOBLoop
bool TIBTOBLoop(int &, int &, int &) const
Definition: LASGlobalLoop.cc:46
LASModuleProfile::SetValue
void SetValue(unsigned int theStripNumber, const double &theValue)
Definition: LASModuleProfile.h:20
LaserAlignment::pedestalProfiles
LASGlobalData< LASModuleProfile > pedestalProfiles
Definition: LaserAlignment.h:202
LASGlobalData::GetTEC2TECEntry
T & GetTEC2TECEntry(int subdetector, int beam, int tecDisk)
Definition: LASGlobalData.h:129
digitizers_cfi.strip
strip
Definition: digitizers_cfi.py:19
LASGlobalData::GetTECEntry
T & GetTECEntry(int subdetector, int tecRing, int beam, int tecDisk)
Definition: LASGlobalData.h:84
LaserAlignment::detectorId
LASGlobalData< unsigned int > detectorId
Definition: LaserAlignment.h:196
LaserAlignment::moduleLoop
LASGlobalLoop moduleLoop
Definition: LaserAlignment.h:229

◆ getTEC2TECNominalBeamOffset()

double LaserAlignment::getTEC2TECNominalBeamOffset ( unsigned int  det,
unsigned int  beam,
unsigned int  disk 
)
private

returns the nominal beam position (strips) in TEC (AT) for the profileJudge

not all TEC-AT modules are hit in the center; this func returns the nominal beam offset locally on a module (in strips) for the ProfileJudge and the LASPeakFinder in strips. (offset = middle of module - nominal position)

the hard coded numbers will later be supplied by a special geometry class..

Definition at line 1314 of file LaserAlignment.cc.

References EcalCondDBWriter_cfi::beam, and Exception.

Referenced by endRunProduce().

1314  {
1315  if (det > 1 || beam > 7 || disk > 5) {
1316  throw cms::Exception("[LaserAlignment::getTEC2TECNominalBeamOffset]")
1317  << " ERROR ** Called with nonexisting parameter set: det " << det << " beam " << beam << " disk " << disk << "."
1318  << std::endl;
1319  }
1320 
1321  const double nominalOffsets[8] = {0., 2.220, -2.221, 0., 2.214, 0., 2.214, -2.217};
1322 
1323  if (det == 0)
1324  return -1. * nominalOffsets[beam];
1325  else
1326  return nominalOffsets[beam];
1327 }
EcalCondDBWriter_cfi.beam
beam
Definition: EcalCondDBWriter_cfi.py:45

◆ getTIBTOBNominalBeamOffset()

double LaserAlignment::getTIBTOBNominalBeamOffset ( unsigned int  det,
unsigned int  beam,
unsigned int  pos 
)
private

returns the nominal beam position (strips) in TOB for the profileJudge

not all TIB & TOB modules are hit in the center; this func returns the nominal beam offset locally on a module (in strips) for the ProfileJudge and the LASPeakFinder in strips. (offset = middle of module - nominal position)

the hard coded numbers will later be supplied by a special geometry class..

Definition at line 1279 of file LaserAlignment.cc.

References EcalCondDBWriter_cfi::beam, Exception, and or.

Referenced by endRunProduce(), and produce().

1279  {
1280  if (det < 2 || det > 3 || beam > 7 || pos > 5) {
1281  throw cms::Exception("[LaserAlignment::getTIBTOBNominalBeamOffset]")
1282  << " ERROR ** Called with nonexisting parameter set: det " << det << " beam " << beam << " pos " << pos << "."
1283  << std::endl;
1284  }
1285 
1286  const double nominalOffsetsTIB[8] = {
1287  0.00035, 2.10687, -2.10827, -0.00173446, 2.10072, -0.00135114, 2.10105, -2.10401};
1288 
1289  // in tob, modules have alternating orientations along the rods.
1290  // this is described by the following pattern.
1291  // (even more confusing, this pattern is inversed for beams 0, 5, 6, 7)
1292  const int orientationPattern[6] = {-1, 1, 1, -1, -1, 1};
1293  const double nominalOffsetsTOB[8] = {0.00217408, 1.58678, 117.733, 119.321, 120.906, 119.328, 117.743, 1.58947};
1294 
1295  if (det == 2)
1296  return (-1. * nominalOffsetsTIB[beam]);
1297 
1298  else {
1299  if (beam == 0 or beam > 4)
1300  return (nominalOffsetsTOB[beam] * orientationPattern[pos]);
1301  else
1302  return (-1. * nominalOffsetsTOB[beam] * orientationPattern[pos]);
1303  }
1304 }
EcalCondDBWriter_cfi.beam
beam
Definition: EcalCondDBWriter_cfi.py:45
or
The Signals That Services Can Subscribe To This is based on ActivityRegistry and is current per Services can connect to the signals distributed by the ActivityRegistry in order to monitor the activity of the application Each possible callback has some defined which we here list in angle e< void, edm::EventID const &, edm::Timestamp const & > We also list in braces which AR_WATCH_USING_METHOD_ is used for those or
Definition: Activities.doc:12

◆ isATBeam()

bool LaserAlignment::isATBeam ( void  )
private

count useable profiles in TIBTOB, operates on LASGlobalData<bool> LaserAlignment::isAcceptedProfile to allow for more elaborate patterns in the future

Definition at line 1252 of file LaserAlignment.cc.

References EcalCondDBWriter_cfi::beam, gather_cfg::cout, LASGlobalData< T >::GetTIBTOBEntry(), isAcceptedProfile, LogDebug, moduleLoop, and LASGlobalLoop::TIBTOBLoop().

Referenced by produce().

1252  {
1253  int numberOfProfiles = 0;
1254 
1255  int det = 2;
1256  int beam = 0;
1257  int pos = 0; // search all TIB/TOB for signals
1258  do {
1259  if (isAcceptedProfile.GetTIBTOBEntry(det, beam, pos) == 1)
1260  numberOfProfiles++;
1261  } while (moduleLoop.TIBTOBLoop(det, beam, pos));
1262 
1263  LogDebug("[LaserAlignment::isATBeam]") << " Found: " << numberOfProfiles << "hits." << std::endl;
1264  std::cout << " [LaserAlignment::isATBeam] -- Found: " << numberOfProfiles << " hits." << std::endl;
1265 
1266  if (numberOfProfiles > 10)
1267  return (true);
1268  return (false);
1269 }
LaserAlignment::isAcceptedProfile
LASGlobalData< int > isAcceptedProfile
Definition: LaserAlignment.h:218
EcalCondDBWriter_cfi.beam
beam
Definition: EcalCondDBWriter_cfi.py:45
LASGlobalData::GetTIBTOBEntry
T & GetTIBTOBEntry(int subdetector, int beam, int tibTobPosition)
Definition: LASGlobalData.h:107
LASGlobalLoop::TIBTOBLoop
bool TIBTOBLoop(int &, int &, int &) const
Definition: LASGlobalLoop.cc:46
gather_cfg.cout
cout
Definition: gather_cfg.py:144
LaserAlignment::moduleLoop
LASGlobalLoop moduleLoop
Definition: LaserAlignment.h:229
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:233

◆ isTECBeam()

bool LaserAlignment::isTECBeam ( void  )
private

decide whether TEC or AT beams have fired

count useable profiles in TEC, operates on LASGlobalData<int> LaserAlignment::isAcceptedProfile to allow for more elaborate patterns in the future

Definition at line 1225 of file LaserAlignment.cc.

References EcalCondDBWriter_cfi::beam, gather_cfg::cout, LASGlobalData< T >::GetTECEntry(), isAcceptedProfile, LogDebug, and relativeConstraints::ring.

Referenced by produce().

1225  {
1226  int numberOfProfiles = 0;
1227 
1228  int ring = 1; // search all ring6 modules for signals
1229  for (int det = 0; det < 2; ++det) {
1230  for (int beam = 0; beam < 8; ++beam) {
1231  for (int disk = 0; disk < 9; ++disk) {
1232  if (isAcceptedProfile.GetTECEntry(det, ring, beam, disk) == 1)
1233  numberOfProfiles++;
1234  }
1235  }
1236  }
1237 
1238  LogDebug("[LaserAlignment::isTECBeam]") << " Found: " << numberOfProfiles << "hits." << std::endl;
1239  std::cout << " [LaserAlignment::isTECBeam] -- Found: " << numberOfProfiles << " hits." << std::endl;
1240 
1241  if (numberOfProfiles > 10)
1242  return (true);
1243  return (false);
1244 }
LaserAlignment::isAcceptedProfile
LASGlobalData< int > isAcceptedProfile
Definition: LaserAlignment.h:218
EcalCondDBWriter_cfi.beam
beam
Definition: EcalCondDBWriter_cfi.py:45
LASGlobalData::GetTECEntry
T & GetTECEntry(int subdetector, int tecRing, int beam, int tecDisk)
Definition: LASGlobalData.h:84
gather_cfg.cout
cout
Definition: gather_cfg.py:144
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:233

◆ produce()

void LaserAlignment::produce ( edm::Event &  theEvent,
edm::EventSetup const &  theSetup 
)
overridevirtual

Implements edm::one::EDProducerBase.

Definition at line 280 of file LaserAlignment.cc.

References EcalCondDBWriter_cfi::beam, TrackerGeomBuilderFromGeometricDet::build(), collectedDataProfiles, gather_cfg::cout, currentDataProfiles, detectorId, edm::EventID::event(), fillDataProfiles(), fillPedestalProfiles(), spr::find(), firstEvent_, gD, geomDetToken_, geomToken_, edm::EventSetup::getData(), edm::EventSetup::getHandle(), LASGlobalData< T >::GetTEC2TECEntry(), LASGlobalData< T >::GetTECEntry(), LASGlobalData< T >::GetTIBTOBEntry(), getTIBTOBNominalBeamOffset(), gprToken_, edm::EventBase::id(), isAcceptedProfile, isATBeam(), LASProfileJudge::IsSignalIn(), isTECBeam(), judge, LASProfileJudge::JudgeProfile(), LogDebug, moduleLoop, numberOfAcceptedProfiles, pedestalProfiles, ptitpToken_, ptpToken_, relativeConstraints::ring, edm::EventID::run(), LASGlobalData< T >::SetTECEntry(), LASGlobalData< T >::SetTIBTOBEntry(), stripPedestalsToken_, LASGlobalLoop::TEC2TECLoop(), tecDoubleHitDetId, LASGlobalLoop::TECLoop(), theAlignableTracker, theDoPedestalSubtraction, theEvents, theGlobalPositionRcd, theProfileNames, theTrackerGeometry, LASGlobalLoop::TIBTOBLoop(), topoToken_, and updateFromInputGeometry.

280  {
281  if (firstEvent_) {
282  //Retrieve tracker topology from geometry
283  const TrackerTopology* const tTopo = &theSetup.getData(topoToken_);
284 
285  // access the tracker
286  gD = theSetup.getHandle(geomDetToken_);
287  theTrackerGeometry = theSetup.getHandle(geomToken_);
288 
289  // access pedestals (from db..) if desired
290  edm::ESHandle<SiStripPedestals> pedestalsHandle;
291  if (theDoPedestalSubtraction) {
292  pedestalsHandle = theSetup.getHandle(stripPedestalsToken_);
293  fillPedestalProfiles(pedestalsHandle);
294  }
295 
296  // global positions
297  theGlobalPositionRcd = &theSetup.getData(gprToken_);
298 
299  // select the reference geometry
300  if (!updateFromInputGeometry) {
301  // the AlignableTracker object is initialized with the ideal geometry
302  const GeometricDet* theGeometricDet = &theSetup.getData(geomDetToken_);
303  const PTrackerParameters* ptp = &theSetup.getData(ptpToken_);
304  const PTrackerAdditionalParametersPerDet* ptitp = &theSetup.getData(ptitpToken_);
305 
306  TrackerGeomBuilderFromGeometricDet trackerBuilder;
307  TrackerGeometry* theRefTracker = trackerBuilder.build(&*theGeometricDet, &*ptitp, *ptp, tTopo);
308 
309  theAlignableTracker = new AlignableTracker(&(*theRefTracker), tTopo);
310  } else {
311  // the AlignableTracker object is initialized with the input geometry from DB
312  theAlignableTracker = new AlignableTracker(&(*theTrackerGeometry), tTopo);
313  }
314 
315  firstEvent_ = false;
316  }
317 
318  LogDebug("LaserAlignment") << "==========================================================="
319  << "\n Private analysis of event #" << theEvent.id().event() << " in run #"
320  << theEvent.id().run();
321 
322  // do the Tracker Statistics to retrieve the current profiles
323  fillDataProfiles(theEvent, theSetup);
324 
325  // index variables for the LASGlobalLoop object
326  int det, ring, beam, disk, pos;
327 
328  //
329  // first pre-loop on selected entries to find out
330  // whether the TEC or the AT beams have fired
331  // (pedestal profiles are left empty if false in cfg)
332  //
333 
334  // TEC+- (only ring 6)
335  ring = 1;
336  for (det = 0; det < 2; ++det) {
337  for (beam = 0; beam < 8; ++beam) {
338  for (disk = 0; disk < 9; ++disk) {
339  if (judge.IsSignalIn(currentDataProfiles.GetTECEntry(det, ring, beam, disk) -
340  pedestalProfiles.GetTECEntry(det, ring, beam, disk),
341  0)) {
342  isAcceptedProfile.SetTECEntry(det, ring, beam, disk, 1);
343  } else { // assume no initialization
344  isAcceptedProfile.SetTECEntry(det, ring, beam, disk, 0);
345  }
346  }
347  }
348  }
349 
350  // TIBTOB
351  det = 2;
352  beam = 0;
353  pos = 0;
354  do {
355  // add current event's data and subtract pedestals
356  if (judge.IsSignalIn(
357  currentDataProfiles.GetTIBTOBEntry(det, beam, pos) - pedestalProfiles.GetTIBTOBEntry(det, beam, pos),
358  getTIBTOBNominalBeamOffset(det, beam, pos))) {
359  isAcceptedProfile.SetTIBTOBEntry(det, beam, pos, 1);
360  } else { // dto.
361  isAcceptedProfile.SetTIBTOBEntry(det, beam, pos, 0);
362  }
363 
364  } while (moduleLoop.TIBTOBLoop(det, beam, pos));
365 
366  // now come the beam finders
367  bool isTECMode = isTECBeam();
368  // LogDebug( " [LaserAlignment::produce]" ) << "LaserAlignment::isTECBeam declares this event " << ( isTECMode ? "" : "NOT " ) << "a TEC event." << std::endl;
369  std::cout << " [LaserAlignment::produce] -- LaserAlignment::isTECBeam declares this event "
370  << (isTECMode ? "" : "NOT ") << "a TEC event." << std::endl;
371 
372  bool isATMode = isATBeam();
373  // LogDebug( " [LaserAlignment::produce]" ) << "LaserAlignment::isATBeam declares this event " << ( isATMode ? "" : "NOT " ) << "an AT event." << std::endl;
374  std::cout << " [LaserAlignment::produce] -- LaserAlignment::isATBeam declares this event " << (isATMode ? "" : "NOT ")
375  << "an AT event." << std::endl;
376 
377  //
378  // now pass the pedestal subtracted profiles to the judge
379  // if they're accepted, add them on the collectedDataProfiles
380  // (pedestal profiles are left empty if false in cfg)
381  //
382 
383  // loop TEC+- modules
384  det = 0;
385  ring = 0;
386  beam = 0;
387  disk = 0;
388  do {
389  LogDebug("[LaserAlignment::produce]")
390  << "Profile is: " << theProfileNames.GetTECEntry(det, ring, beam, disk) << "." << std::endl;
391 
392  // this now depends on the AT/TEC mode, is this a doubly hit module? -> look for it in vector<int> tecDoubleHitDetId
393  // (ring == 0 not necessary but makes it a little faster)
394  if (ring == 0 &&
395  find(tecDoubleHitDetId.begin(), tecDoubleHitDetId.end(), detectorId.GetTECEntry(det, ring, beam, disk)) !=
396  tecDoubleHitDetId.end()) {
397  if (isTECMode) { // add profile to TEC collection
398  // add current event's data and subtract pedestals
399  if (judge.JudgeProfile(currentDataProfiles.GetTECEntry(det, ring, beam, disk) -
400  pedestalProfiles.GetTECEntry(det, ring, beam, disk),
401  0)) {
402  collectedDataProfiles.GetTECEntry(det, ring, beam, disk) +=
403  currentDataProfiles.GetTECEntry(det, ring, beam, disk) -
404  pedestalProfiles.GetTECEntry(det, ring, beam, disk);
405  numberOfAcceptedProfiles.GetTECEntry(det, ring, beam, disk)++;
406  }
407  }
408  }
409 
410  else { // not a doubly hit module, don't care about the mode
411  // add current event's data and subtract pedestals
412  if (judge.JudgeProfile(currentDataProfiles.GetTECEntry(det, ring, beam, disk) -
413  pedestalProfiles.GetTECEntry(det, ring, beam, disk),
414  0)) {
415  collectedDataProfiles.GetTECEntry(det, ring, beam, disk) +=
416  currentDataProfiles.GetTECEntry(det, ring, beam, disk) -
417  pedestalProfiles.GetTECEntry(det, ring, beam, disk);
418  numberOfAcceptedProfiles.GetTECEntry(det, ring, beam, disk)++;
419  }
420  }
421 
422  } while (moduleLoop.TECLoop(det, ring, beam, disk));
423 
424  // loop TIB/TOB modules
425  det = 2;
426  beam = 0;
427  pos = 0;
428  do {
429  LogDebug("[LaserAlignment::produce]")
430  << "Profile is: " << theProfileNames.GetTIBTOBEntry(det, beam, pos) << "." << std::endl;
431 
432  // add current event's data and subtract pedestals
433  if (judge.JudgeProfile(
434  currentDataProfiles.GetTIBTOBEntry(det, beam, pos) - pedestalProfiles.GetTIBTOBEntry(det, beam, pos),
435  getTIBTOBNominalBeamOffset(det, beam, pos))) {
436  collectedDataProfiles.GetTIBTOBEntry(det, beam, pos) +=
437  currentDataProfiles.GetTIBTOBEntry(det, beam, pos) - pedestalProfiles.GetTIBTOBEntry(det, beam, pos);
438  numberOfAcceptedProfiles.GetTIBTOBEntry(det, beam, pos)++;
439  }
440 
441  } while (moduleLoop.TIBTOBLoop(det, beam, pos));
442 
443  // loop TEC2TEC modules
444  det = 0;
445  beam = 0;
446  disk = 0;
447  do {
448  LogDebug("[LaserAlignment::produce]")
449  << "Profile is: " << theProfileNames.GetTEC2TECEntry(det, beam, disk) << "." << std::endl;
450 
451  // this again depends on the AT/TEC mode, is this a doubly hit module?
452  // (ring == 0 not necessary but makes it a little faster)
453  if (ring == 0 &&
454  find(tecDoubleHitDetId.begin(), tecDoubleHitDetId.end(), detectorId.GetTECEntry(det, ring, beam, disk)) !=
455  tecDoubleHitDetId.end()) {
456  if (isATMode) { // add profile to TEC2TEC collection
457  // add current event's data and subtract pedestals
458  if (judge.JudgeProfile(currentDataProfiles.GetTEC2TECEntry(det, beam, disk) -
459  pedestalProfiles.GetTEC2TECEntry(det, beam, disk),
460  0)) {
461  collectedDataProfiles.GetTEC2TECEntry(det, beam, disk) +=
462  currentDataProfiles.GetTEC2TECEntry(det, beam, disk) - pedestalProfiles.GetTEC2TECEntry(det, beam, disk);
463  numberOfAcceptedProfiles.GetTEC2TECEntry(det, beam, disk)++;
464  }
465  }
466 
467  }
468 
469  else { // not a doubly hit module, don't care about the mode
470  // add current event's data and subtract pedestals
471  if (judge.JudgeProfile(
472  currentDataProfiles.GetTEC2TECEntry(det, beam, disk) - pedestalProfiles.GetTEC2TECEntry(det, beam, disk),
473  0)) {
474  collectedDataProfiles.GetTEC2TECEntry(det, beam, disk) +=
475  currentDataProfiles.GetTEC2TECEntry(det, beam, disk) - pedestalProfiles.GetTEC2TECEntry(det, beam, disk);
476  numberOfAcceptedProfiles.GetTEC2TECEntry(det, beam, disk)++;
477  }
478  }
479 
480  } while (moduleLoop.TEC2TECLoop(det, beam, disk));
481 
482  // total event number counter
483  theEvents++;
484 }
LASProfileJudge::IsSignalIn
bool IsSignalIn(const LASModuleProfile &, double)
Definition: LASProfileJudge.cc:24
LaserAlignment::isATBeam
bool isATBeam(void)
Definition: LaserAlignment.cc:1252
LaserAlignment::fillDataProfiles
void fillDataProfiles(edm::Event const &, edm::EventSetup const &)
fill profiles from SiStrip(Raw)Digi container
Definition: LaserAlignment.cc:941
LaserAlignment::ptpToken_
const edm::ESGetToken< PTrackerParameters, PTrackerParametersRcd > ptpToken_
Definition: LaserAlignment.h:134
LaserAlignment::geomDetToken_
const edm::ESGetToken< GeometricDet, IdealGeometryRecord > geomDetToken_
Definition: LaserAlignment.h:133
LaserAlignment::stripPedestalsToken_
const edm::ESGetToken< SiStripPedestals, SiStripPedestalsRcd > stripPedestalsToken_
Definition: LaserAlignment.h:137
LaserAlignment::ptitpToken_
const edm::ESGetToken< PTrackerAdditionalParametersPerDet, PTrackerAdditionalParametersPerDetRcd > ptitpToken_
Definition: LaserAlignment.h:135
LaserAlignment::isTECBeam
bool isTECBeam(void)
decide whether TEC or AT beams have fired
Definition: LaserAlignment.cc:1225
TrackerGeomBuilderFromGeometricDet::build
TrackerGeometry * build(const GeometricDet *gd, const PTrackerAdditionalParametersPerDet *ptitp, const PTrackerParameters &ptp, const TrackerTopology *tTopo)
Definition: TrackerGeomBuilderFromGeometricDet.cc:44
LaserAlignment::getTIBTOBNominalBeamOffset
double getTIBTOBNominalBeamOffset(unsigned int, unsigned int, unsigned int)
returns the nominal beam position (strips) in TOB for the profileJudge
Definition: LaserAlignment.cc:1279
LaserAlignment::theTrackerGeometry
edm::ESHandle< TrackerGeometry > theTrackerGeometry
Definition: LaserAlignment.h:237
LASGlobalLoop::TEC2TECLoop
bool TEC2TECLoop(int &, int &, int &) const
Definition: LASGlobalLoop.cc:74
LaserAlignment::numberOfAcceptedProfiles
LASGlobalData< int > numberOfAcceptedProfiles
Definition: LaserAlignment.h:214
LaserAlignment::theDoPedestalSubtraction
bool theDoPedestalSubtraction
config switch
Definition: LaserAlignment.h:143
LaserAlignment::currentDataProfiles
LASGlobalData< LASModuleProfile > currentDataProfiles
data profiles for the current event
Definition: LaserAlignment.h:205
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:19
LaserAlignment::isAcceptedProfile
LASGlobalData< int > isAcceptedProfile
Definition: LaserAlignment.h:218
EcalCondDBWriter_cfi.beam
beam
Definition: EcalCondDBWriter_cfi.py:45
LaserAlignment::fillPedestalProfiles
void fillPedestalProfiles(edm::ESHandle< SiStripPedestals > &)
fill pedestals from dbase
Definition: LaserAlignment.cc:1173
LaserAlignment::geomToken_
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > geomToken_
Definition: LaserAlignment.h:132
LASGlobalLoop::TECLoop
bool TECLoop(int &, int &, int &, int &) const
Definition: LASGlobalLoop.cc:13
LaserAlignment::theGlobalPositionRcd
edm::ESHandle< Alignments > theGlobalPositionRcd
Definition: LaserAlignment.h:238
LaserAlignment::topoToken_
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > topoToken_
Tokens for ESconsumes.
Definition: LaserAlignment.h:131
LaserAlignment::tecDoubleHitDetId
std::vector< unsigned int > tecDoubleHitDetId
Definition: LaserAlignment.h:199
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:59
LASGlobalData::GetTIBTOBEntry
T & GetTIBTOBEntry(int subdetector, int beam, int tibTobPosition)
Definition: LASGlobalData.h:107
LASProfileJudge::JudgeProfile
bool JudgeProfile(const LASModuleProfile &, double)
Definition: LASProfileJudge.cc:45
LASGlobalData::SetTIBTOBEntry
void SetTIBTOBEntry(int subdetector, int beam, int tibTobPosition, T)
Definition: LASGlobalData.h:173
LASGlobalLoop::TIBTOBLoop
bool TIBTOBLoop(int &, int &, int &) const
Definition: LASGlobalLoop.cc:46
edm::EventID::run
RunNumber_t run() const
Definition: EventID.h:38
LaserAlignment::pedestalProfiles
LASGlobalData< LASModuleProfile > pedestalProfiles
Definition: LaserAlignment.h:202
LASGlobalData::GetTEC2TECEntry
T & GetTEC2TECEntry(int subdetector, int beam, int tecDisk)
Definition: LASGlobalData.h:129
LaserAlignment::theAlignableTracker
AlignableTracker * theAlignableTracker
Definition: LaserAlignment.h:240
LaserAlignment::theEvents
int theEvents
counter for the total number of events processed
Definition: LaserAlignment.h:140
LASGlobalData::GetTECEntry
T & GetTECEntry(int subdetector, int tecRing, int beam, int tecDisk)
Definition: LASGlobalData.h:84
LaserAlignment::theProfileNames
LASGlobalData< std::string > theProfileNames
Definition: LaserAlignment.h:211
LaserAlignment::collectedDataProfiles
LASGlobalData< LASModuleProfile > collectedDataProfiles
Definition: LaserAlignment.h:208
LaserAlignment::detectorId
LASGlobalData< unsigned int > detectorId
Definition: LaserAlignment.h:196
LaserAlignment::gprToken_
const edm::ESGetToken< Alignments, GlobalPositionRcd > gprToken_
Definition: LaserAlignment.h:136
gather_cfg.cout
cout
Definition: gather_cfg.py:144
LASGlobalData::SetTECEntry
void SetTECEntry(int subdetector, int tecRing, int beam, int tecDisk, T)
Definition: LASGlobalData.h:150
LaserAlignment::gD
edm::ESHandle< GeometricDet > gD
tracker geometry;
Definition: LaserAlignment.h:236
edm::EventID::event
EventNumber_t event() const
Definition: EventID.h:40
LaserAlignment::moduleLoop
LASGlobalLoop moduleLoop
Definition: LaserAlignment.h:229
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:233
LaserAlignment::judge
LASProfileJudge judge
Definition: LaserAlignment.h:190
LaserAlignment::updateFromInputGeometry
bool updateFromInputGeometry
config switch
Definition: LaserAlignment.h:161

◆ testRoutine()

void LaserAlignment::testRoutine ( void  )

for debugging & testing only, will disappear..

this function is for debugging and testing only and will disappear..

Definition at line 1686 of file LaserAlignment.cc.

References EcalCondDBWriter_cfi::beam, gather_cfg::cout, detectorId, LASGlobalData< T >::GetTEC2TECEntry(), LASGlobalData< T >::GetTECEntry(), LASGlobalData< T >::GetTIBTOBEntry(), runTauDisplay::gp, TrackerGeometry::idToDet(), moduleLoop, CosmicsPD_Skims::radius, relativeConstraints::ring, LASGlobalLoop::TEC2TECLoop(), LASGlobalLoop::TECLoop(), theTrackerGeometry, theZ, and LASGlobalLoop::TIBTOBLoop().

1686  {
1687  // tracker geom. object for calculating the global beam positions
1688  const TrackerGeometry& theTracker(*theTrackerGeometry);
1689 
1690  const double atPhiPositions[8] = {0.392699, 1.289799, 1.851794, 2.748894, 3.645995, 4.319690, 5.216791, 5.778784};
1691  const double tecPhiPositions[8] = {0.392699, 1.178097, 1.963495, 2.748894, 3.534292, 4.319690, 5.105088, 5.890486};
1692  const double zPositions[9] = {125.0, 139.0, 153.0, 167.0, 181.0, 198.5, 217.5, 238.0, 259.5};
1693  const double zPositionsTIB[6] = {62.0, 38.0, 18.0, -10.0, -34.0, -54.0};
1694  const double zPositionsTOB[6] = {104.0, 58.0, 22.0, -14.0, -50.0, -86.0};
1695 
1696  int det, beam, disk, pos, ring;
1697 
1698  // loop TEC+- internal
1699  det = 0;
1700  ring = 0;
1701  beam = 0;
1702  disk = 0;
1703  do {
1704  const double radius = ring ? 84.0 : 56.4;
1705 
1706  // access the tracker geometry for this module
1707  const DetId theDetId(detectorId.GetTECEntry(det, ring, beam, disk));
1708  const StripGeomDetUnit* const theStripDet = dynamic_cast<const StripGeomDetUnit*>(theTracker.idToDet(theDetId));
1709 
1710  if (theStripDet) {
1711  const GlobalPoint gp(GlobalPoint::Cylindrical(radius, tecPhiPositions[beam], zPositions[disk]));
1712 
1713  const LocalPoint lp(theStripDet->surface().toLocal(gp));
1714  std::cout << "__TEC: " << 256. - theStripDet->specificTopology().strip(lp)
1715  << std::endl;
1716  }
1717 
1718  } while (moduleLoop.TECLoop(det, ring, beam, disk));
1719 
1720  // loop TIBTOB
1721  det = 2;
1722  beam = 0;
1723  pos = 0;
1724  do {
1725  const double radius =
1726  (det == 2 ? 51.4 : 58.4);
1727  const double theZ = (det == 2 ? zPositionsTIB[pos] : zPositionsTOB[pos]);
1728 
1729  // access the tracker geometry for this module
1730  const DetId theDetId(detectorId.GetTIBTOBEntry(det, beam, pos));
1731  const StripGeomDetUnit* const theStripDet = dynamic_cast<const StripGeomDetUnit*>(theTracker.idToDet(theDetId));
1732 
1733  if (theStripDet) {
1734  const GlobalPoint gp(GlobalPoint::Cylindrical(radius, atPhiPositions[beam], theZ));
1735 
1736  const LocalPoint lp(theStripDet->surface().toLocal(gp));
1737  std::cout << "__TIBTOB det " << det << " beam " << beam << " pos " << pos << " "
1738  << 256. - theStripDet->specificTopology().strip(lp);
1739  std::cout << " " << theStripDet->position().perp() << std::endl;
1740  }
1741 
1742  } while (moduleLoop.TIBTOBLoop(det, beam, pos));
1743 
1744  // loop TEC2TEC
1745  det = 0;
1746  beam = 0;
1747  disk = 0;
1748  do {
1749  const double radius = 56.4;
1750 
1751  // access the tracker geometry for this module
1752  const DetId theDetId(detectorId.GetTEC2TECEntry(det, beam, disk));
1753  const StripGeomDetUnit* const theStripDet = dynamic_cast<const StripGeomDetUnit*>(theTracker.idToDet(theDetId));
1754 
1755  if (theStripDet) {
1756  const GlobalPoint gp(GlobalPoint::Cylindrical(radius, atPhiPositions[beam], zPositions[disk]));
1757 
1758  const LocalPoint lp(theStripDet->surface().toLocal(gp));
1759  std::cout << "__TEC2TEC det " << det << " beam " << beam << " disk " << disk << " "
1760  << 256. - theStripDet->specificTopology().strip(lp) << std::endl;
1761  }
1762 
1763  } while (moduleLoop.TEC2TECLoop(det, beam, disk));
1764 }
theZ
T theZ
Definition: Basic3DVectorLD.h:243
LaserAlignment::theTrackerGeometry
edm::ESHandle< TrackerGeometry > theTrackerGeometry
Definition: LaserAlignment.h:237
LASGlobalLoop::TEC2TECLoop
bool TEC2TECLoop(int &, int &, int &) const
Definition: LASGlobalLoop.cc:74
EcalCondDBWriter_cfi.beam
beam
Definition: EcalCondDBWriter_cfi.py:45
CosmicsPD_Skims.radius
radius
Definition: CosmicsPD_Skims.py:135
LASGlobalLoop::TECLoop
bool TECLoop(int &, int &, int &, int &) const
Definition: LASGlobalLoop.cc:13
LASGlobalData::GetTIBTOBEntry
T & GetTIBTOBEntry(int subdetector, int beam, int tibTobPosition)
Definition: LASGlobalData.h:107
LASGlobalLoop::TIBTOBLoop
bool TIBTOBLoop(int &, int &, int &) const
Definition: LASGlobalLoop.cc:46
DetId
Definition: DetId.h:17
LASGlobalData::GetTEC2TECEntry
T & GetTEC2TECEntry(int subdetector, int beam, int tecDisk)
Definition: LASGlobalData.h:129
LASGlobalData::GetTECEntry
T & GetTECEntry(int subdetector, int tecRing, int beam, int tecDisk)
Definition: LASGlobalData.h:84
LaserAlignment::detectorId
LASGlobalData< unsigned int > detectorId
Definition: LaserAlignment.h:196
gather_cfg.cout
cout
Definition: gather_cfg.py:144
LaserAlignment::moduleLoop
LASGlobalLoop moduleLoop
Definition: LaserAlignment.h:229

Member Data Documentation

◆ collectedDataProfiles

LASGlobalData<LASModuleProfile> LaserAlignment::collectedDataProfiles
private

Definition at line 208 of file LaserAlignment.h.

Referenced by beginJob(), endRunProduce(), and produce().

◆ currentDataProfiles

LASGlobalData<LASModuleProfile> LaserAlignment::currentDataProfiles
private

data profiles for the current event

Definition at line 205 of file LaserAlignment.h.

Referenced by beginJob(), fillDataProfiles(), and produce().

◆ detectorId

LASGlobalData<unsigned int> LaserAlignment::detectorId
private

Definition at line 196 of file LaserAlignment.h.

Referenced by ApplyATMaskingCorrections(), ApplyEndcapMaskingCorrections(), endRunProduce(), fillDataProfiles(), fillDetectorId(), fillPedestalProfiles(), produce(), and testRoutine().

◆ enableJudgeZeroFilter

bool LaserAlignment::enableJudgeZeroFilter
private

config switch

Definition at line 155 of file LaserAlignment.h.

Referenced by LaserAlignment().

◆ firstEvent_

bool LaserAlignment::firstEvent_
private

Definition at line 244 of file LaserAlignment.h.

Referenced by beginJob(), and produce().

◆ gD

edm::ESHandle<GeometricDet> LaserAlignment::gD
private

tracker geometry;

Definition at line 236 of file LaserAlignment.h.

Referenced by produce().

◆ geomDetToken_

const edm::ESGetToken<GeometricDet, IdealGeometryRecord> LaserAlignment::geomDetToken_
private

Definition at line 133 of file LaserAlignment.h.

Referenced by produce().

◆ geomToken_

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> LaserAlignment::geomToken_
private

Definition at line 132 of file LaserAlignment.h.

Referenced by produce().

◆ gprToken_

const edm::ESGetToken<Alignments, GlobalPositionRcd> LaserAlignment::gprToken_
private

Definition at line 136 of file LaserAlignment.h.

Referenced by produce().

◆ isAcceptedProfile

LASGlobalData<int> LaserAlignment::isAcceptedProfile
private

Definition at line 218 of file LaserAlignment.h.

Referenced by beginJob(), isATBeam(), isTECBeam(), and produce().

◆ judge

LASProfileJudge LaserAlignment::judge
private

Definition at line 190 of file LaserAlignment.h.

Referenced by LaserAlignment(), and produce().

◆ judgeOverdriveThreshold

unsigned int LaserAlignment::judgeOverdriveThreshold
private

config parameters for the LASProfileJudge

Definition at line 158 of file LaserAlignment.h.

Referenced by LaserAlignment().

◆ misalignedByRefGeometry

bool LaserAlignment::misalignedByRefGeometry
private

config switch

Definition at line 164 of file LaserAlignment.h.

Referenced by endRunProduce(), and LaserAlignment().

◆ moduleLoop

LASGlobalLoop LaserAlignment::moduleLoop
private

Definition at line 229 of file LaserAlignment.h.

Referenced by ApplyATMaskingCorrections(), ApplyEndcapMaskingCorrections(), beginJob(), CalculateNominalCoordinates(), endRunProduce(), fillDataProfiles(), fillPedestalProfiles(), isATBeam(), produce(), and testRoutine().

◆ nominalCoordinates

LASGlobalData<LASCoordinateSet> LaserAlignment::nominalCoordinates
private

Definition at line 225 of file LaserAlignment.h.

Referenced by ApplyATMaskingCorrections(), ApplyEndcapMaskingCorrections(), CalculateNominalCoordinates(), and endRunProduce().

◆ numberOfAcceptedProfiles

LASGlobalData<int> LaserAlignment::numberOfAcceptedProfiles
private

Definition at line 214 of file LaserAlignment.h.

Referenced by beginJob(), DumpHitmaps(), endRunProduce(), and produce().

◆ peakFinderThreshold

double LaserAlignment::peakFinderThreshold
private

config parameter

Definition at line 152 of file LaserAlignment.h.

Referenced by endRunProduce(), and LaserAlignment().

◆ pedestalProfiles

LASGlobalData<LASModuleProfile> LaserAlignment::pedestalProfiles
private

Definition at line 202 of file LaserAlignment.h.

Referenced by beginJob(), fillPedestalProfiles(), and produce().

◆ ptitpToken_

const edm::ESGetToken<PTrackerAdditionalParametersPerDet, PTrackerAdditionalParametersPerDetRcd> LaserAlignment::ptitpToken_
private

Definition at line 135 of file LaserAlignment.h.

Referenced by produce().

◆ ptpToken_

const edm::ESGetToken<PTrackerParameters, PTrackerParametersRcd> LaserAlignment::ptpToken_
private

Definition at line 134 of file LaserAlignment.h.

Referenced by produce().

◆ singleModulesDir

TDirectory* LaserAlignment::singleModulesDir
private

Definition at line 233 of file LaserAlignment.h.

Referenced by beginJob().

◆ stripPedestalsToken_

const edm::ESGetToken<SiStripPedestals, SiStripPedestalsRcd> LaserAlignment::stripPedestalsToken_
private

Definition at line 137 of file LaserAlignment.h.

Referenced by produce().

◆ summedHistograms

LASGlobalData<TH1D*> LaserAlignment::summedHistograms
private

Definition at line 222 of file LaserAlignment.h.

Referenced by beginJob(), and endRunProduce().

◆ tecDoubleHitDetId

std::vector<unsigned int> LaserAlignment::tecDoubleHitDetId
private

Definition at line 199 of file LaserAlignment.h.

Referenced by fillDetectorId(), and produce().

◆ theAlignableTracker

AlignableTracker* LaserAlignment::theAlignableTracker
private

Definition at line 240 of file LaserAlignment.h.

Referenced by endRunProduce(), produce(), and ~LaserAlignment().

◆ theAlignRecordName

std::string LaserAlignment::theAlignRecordName
private

Definition at line 242 of file LaserAlignment.h.

Referenced by endRunProduce().

◆ theApplyBeamKinkCorrections

bool LaserAlignment::theApplyBeamKinkCorrections
private

config switch

Definition at line 149 of file LaserAlignment.h.

Referenced by endRunProduce(), and LaserAlignment().

◆ theCompression

int LaserAlignment::theCompression
private

config parameter (histograms file compression level)

Definition at line 176 of file LaserAlignment.h.

Referenced by beginJob(), and LaserAlignment().

◆ theDigiProducersList

std::vector<edm::ParameterSet> LaserAlignment::theDigiProducersList
private

Definition at line 170 of file LaserAlignment.h.

Referenced by fillDataProfiles().

◆ theDoPedestalSubtraction

bool LaserAlignment::theDoPedestalSubtraction
private

config switch

Definition at line 143 of file LaserAlignment.h.

Referenced by LaserAlignment(), and produce().

◆ theErrorRecordName

std::string LaserAlignment::theErrorRecordName
private

Definition at line 242 of file LaserAlignment.h.

Referenced by endRunProduce().

◆ theEvents

int LaserAlignment::theEvents
private

counter for the total number of events processed

Definition at line 140 of file LaserAlignment.h.

Referenced by endRunProduce(), and produce().

◆ theFile

TFile* LaserAlignment::theFile
private

Tree stuff.

Definition at line 232 of file LaserAlignment.h.

Referenced by beginJob(), and ~LaserAlignment().

◆ theFileName

std::string LaserAlignment::theFileName
private

config parameter (histograms file output name)

Definition at line 179 of file LaserAlignment.h.

Referenced by beginJob(), and LaserAlignment().

◆ theGlobalPositionRcd

edm::ESHandle<Alignments> LaserAlignment::theGlobalPositionRcd
private

Definition at line 238 of file LaserAlignment.h.

Referenced by produce().

◆ theLasConstants

LASConstants LaserAlignment::theLasConstants
private

Definition at line 193 of file LaserAlignment.h.

Referenced by endRunProduce().

◆ theMaskAtModules

std::vector<unsigned int> LaserAlignment::theMaskAtModules
private

Definition at line 183 of file LaserAlignment.h.

Referenced by ApplyATMaskingCorrections(), endRunProduce(), and LaserAlignment().

◆ theMaskTecModules

std::vector<unsigned int> LaserAlignment::theMaskTecModules
private

config parameters

Definition at line 182 of file LaserAlignment.h.

Referenced by ApplyEndcapMaskingCorrections(), endRunProduce(), and LaserAlignment().

◆ theProfileNames

LASGlobalData<std::string> LaserAlignment::theProfileNames
private

Definition at line 211 of file LaserAlignment.h.

Referenced by beginJob(), and produce().

◆ theSaveHistograms

bool LaserAlignment::theSaveHistograms
private

config switch

Definition at line 173 of file LaserAlignment.h.

Referenced by beginJob(), LaserAlignment(), and ~LaserAlignment().

◆ theSetNominalStrips

bool LaserAlignment::theSetNominalStrips
private

config switch

Definition at line 186 of file LaserAlignment.h.

Referenced by endRunProduce(), and LaserAlignment().

◆ theStoreToDB

bool LaserAlignment::theStoreToDB
private

config switch

Definition at line 167 of file LaserAlignment.h.

Referenced by endRunProduce(), and LaserAlignment().

◆ theTrackerGeometry

edm::ESHandle<TrackerGeometry> LaserAlignment::theTrackerGeometry
private

Definition at line 237 of file LaserAlignment.h.

Referenced by endRunProduce(), produce(), and testRoutine().

◆ theUseMinuitAlgorithm

bool LaserAlignment::theUseMinuitAlgorithm
private

config switch

Definition at line 146 of file LaserAlignment.h.

Referenced by endRunProduce(), and LaserAlignment().

◆ topoToken_

const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> LaserAlignment::topoToken_
private

Tokens for ESconsumes.

Definition at line 131 of file LaserAlignment.h.

Referenced by produce().

◆ updateFromInputGeometry

bool LaserAlignment::updateFromInputGeometry
private

config switch

Definition at line 161 of file LaserAlignment.h.

Referenced by endRunProduce(), LaserAlignment(), and produce().

Generated for CMSSW Reference Manual by doxygen 1.8.14