#include <DTMapGenerator.h>

Inheritance diagram for DTMapGenerator:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &setup) override

void	beginJob () override

	DTMapGenerator (const edm::ParameterSet &pset)
	Constructor. More...

void	endJob () override

	~DTMapGenerator () override=default
	Destructor. More...

Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () 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
bool	checkWireExist (const std::set< DTWireId > &wireMap, int wheel, int station, int sector, int sl, int layer, int wire)

Private Attributes
std::string	inputMapName

std::string	outputMapName

int	rosType

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

Class which creates a textual map of the hardware channels in the software detIds

Author: G. Cerminara S. Bolognesi - INFN Torino

Definition at line 17 of file DTMapGenerator.h.

Constructor & Destructor Documentation

◆ DTMapGenerator()

DTMapGenerator::DTMapGenerator ( const edm::ParameterSet & pset )

Constructor.

Definition at line 18 of file DTMapGenerator.cc.

References gather_cfg::cout, and muonDTDigis_cfi::pset.

                                                        {
   // The output file with the map
   outputMapName = pset.getUntrackedParameter<string>("outputFileName", "output.map");
   // The input file with the base map (DDU ROS -> Wheel, Station, Sector)
   inputMapName = pset.getUntrackedParameter<string>("inputFileName", "basemap.txt");
   //The ros type: ROS8 for commissioning, ROS25 otherwise
   rosType = pset.getUntrackedParameter<int>("rosType", 25);
   if (rosType != 8 && rosType != 25) {
     cout << "[DTMapGenerator]: wrong ros type (8 for commissioning, 25 otherwise)" << endl;
     abort();
   }
 }

◆ ~DTMapGenerator()

DTMapGenerator::~DTMapGenerator ( )

overridedefault

Destructor.

Member Function Documentation

◆ analyze()

void DTMapGenerator::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

inlineoverridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 29 of file DTMapGenerator.h.

29 {}

◆ beginJob()

void DTMapGenerator::beginJob ( )

inlineoverridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 27 of file DTMapGenerator.h.

27 {}

◆ checkWireExist()

bool DTMapGenerator::checkWireExist	(	const std::set< DTWireId > &	wireMap,
		int	wheel,
		int	station,
		int	sector,
		int	sl,
		int	layer,
		int	wire
	)

private

Definition at line 164 of file DTMapGenerator.cc.

References gather_cfg::cout, relativeConstraints::station, and makeMuonMisalignmentScenario::wheel.

                                                                                                    {
   DTWireId wireId(wheel, station, sector, sl, layer, wire);
   if (wireMap.find(wireId) == wireMap.end()) {
     cout << "Skipping channel: Wh: " << wheel << ", st: " << station << ", sec: " << sector << ", sl: " << sl
          << ", lay: " << layer << ", wire: " << wire << endl;
     return false;
   }
 
   return true;
 }

◆ endJob()

void DTMapGenerator::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 31 of file DTMapGenerator.cc.

References gather_cfg::cout, MillePedeFileConverter_cfg::fileName, mps_splice::line, relativeConstraints::station, and makeMuonMisalignmentScenario::wheel.

Referenced by o2olib.O2ORunMgr::executeJob().

                             {
   cout << "DTMapGenerator: Output Map: " << outputMapName << " ROS Type: " << rosType << endl;
 
   // Read the existing wires
   ifstream existingChannels("/afs/cern.ch/cms/Physics/muon/CMSSW/DT/channelsMaps/existing_channels.txt");
 
   set<DTWireId> wireMap;  //FIXME:MAYBE YOU NEED THE > and == operators to use set?
 
   // Read the map between DDU - ROS and Chambers
   string lineMap;
   while (getline(existingChannels, lineMap)) {
     if (lineMap.empty() || lineMap[0] == '#')
       continue;  // Skip comments and empty lines
     stringstream linestr;
     linestr << lineMap;
     int wheelEx, stationEx, sectorEx, slEx, layerEx, wireEx;
     linestr >> wheelEx >> sectorEx >> stationEx >> slEx >> layerEx >> wireEx;
     DTWireId wireIdEx(wheelEx, stationEx, sectorEx, slEx, layerEx, wireEx);
     wireMap.insert(wireIdEx);
   }
   cout << "Map size: " << wireMap.size() << endl;
 
   // The map between DDU - ROS and Chambers
   ifstream skeletonMap(inputMapName.c_str());
 
   // The output map in the CMSSW format
   ofstream outputMap(outputMapName.c_str());
 
   // Read the map between DDU - ROS and Chambers
   string line;
   while (getline(skeletonMap, line)) {
     if (line.empty() || line[0] == '#')
       continue;  // Skip comments and empty lines
     stringstream linestr;
     linestr << line;
     int ddu, ros, wheel, station, sector;
     linestr >> ddu >> ros >> wheel >> station >> sector;
     cout << "DDU: " << ddu << endl
          << "ROS: " << ros << endl
          << "Connected to chamber in Wh: " << wheel << " St: " << station << " Sec: " << sector << endl;
 
     int previousROB = -1;
     int robCounter = -1;
     // The chamber map in ORCA commissioning format
     string fileName;
     stringstream nameTmp;
     nameTmp << "/afs/cern.ch/cms/Physics/muon/CMSSW/DT/channelsMaps/templates/MB" << station << "_" << sector << ".map";
     nameTmp >> fileName;
     ifstream chamberMap(fileName.c_str());
 
     string lineChamberMap;
     while (getline(chamberMap, lineChamberMap)) {
       if (lineChamberMap.empty() || lineChamberMap[0] == '#')
         continue;  // Skip comments and empty lines
       stringstream chamberMapStr;
       chamberMapStr << lineChamberMap;
 
       int rob, tdc, tdcChannel, sl, layer, wire;
       int unusedRos, unusedChamberCode;
       int outRob = -1;
       chamberMapStr >> unusedRos >> rob >> tdc >> tdcChannel >> unusedChamberCode >> sl >> layer >> wire;
 
       // Check if the channel really exists
       if (!checkWireExist(wireMap, wheel, station, sector, sl, layer, wire))
         continue;
 
       if (rob > previousROB) {
         previousROB = rob;
         robCounter++;
       } else if (rob < previousROB) {
         cout << "Error: ROB number is not uniformly increasing!" << endl;
         abort();
       }
       // Set the ROB id within the ros
       if (rosType == 25) {
         if (station == 1) {  //MB1
           outRob = robCounter;
         } else if (station == 2) {  //MB2
           outRob = robCounter + 6;
         } else if (station == 3) {  //MB3
           if (robCounter < 3)
             outRob = robCounter + 12;
           else if (robCounter == 3)
             outRob = 24;
           else if (robCounter > 3)
             outRob = robCounter + 11;
         } else if (station == 4) {  //MB4
           if (sector == 14) {
             if (robCounter == 3) {
               continue;
             }
             outRob = robCounter + 18;
           } else if (sector == 10) {
             if (robCounter == 3) {
               continue;
             } else if (robCounter == 0) {
               outRob = 21;
             } else {
               outRob = robCounter + 21;
             }
           } else if (sector == 4) {
             if (robCounter == 3 || robCounter == 4) {
               continue;
             }
             outRob = robCounter + 18;
           } else if (sector == 13) {
             if (robCounter == 3 || robCounter == 4) {
               continue;
             } else if (robCounter == 0) {
               outRob = 21;
             } else {
               outRob = robCounter + 21;
             }
           } else if (sector == 11 || sector == 9) {
             outRob = robCounter + 18;
             if (robCounter == 3) {
               continue;
             }
           }
           //else if(sector==12 || sector == 8 || sector == 7 || sector == 6 || sector == 5 ||  sector == 3 || sector == 2 ||sector == 1 ){
           else {
             outRob = robCounter + 18;
           }
         }
       } else {
         outRob = rob;
       }
       outputMap << ddu << " " << ros << " " << outRob << " " << tdc << " " << tdcChannel << " " << wheel << " "
                 << station << " " << sector << " " << sl << " " << layer << " " << wire << endl;
     }
   }
 }

Member Data Documentation

◆ inputMapName

std::string DTMapGenerator::inputMapName

private

Definition at line 43 of file DTMapGenerator.h.

◆ outputMapName

std::string DTMapGenerator::outputMapName

private

Definition at line 41 of file DTMapGenerator.h.

◆ rosType

int DTMapGenerator::rosType

private

Definition at line 45 of file DTMapGenerator.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ DTMapGenerator()

◆ ~DTMapGenerator()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ checkWireExist()

◆ endJob()

Member Data Documentation

◆ inputMapName

◆ outputMapName

◆ rosType