DTMapGenerator Class Reference

#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
	Destructor. More...
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
	~EDAnalyzer () 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
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::vector< ModuleDescription const * > &modules, 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
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::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
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<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)

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

(

const edm::ParameterSet &

pset

)

Constructor.

Definition at line 18 of file DTMapGenerator.cc.

References gather_cfg::cout, and edm::ParameterSet::getUntrackedParameter().

                                                       {
  // 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 ( )

override

Destructor.

Definition at line 31 of file DTMapGenerator.cc.

{}

Member Function Documentation

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

inlineoverride

Definition at line 29 of file DTMapGenerator.h.

References checkWireExist(), endJob(), relativeConstraints::station, and makeMuonMisalignmentScenario::wheel.

{}

void DTMapGenerator::beginJob ( void  )

inlineoverridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 27 of file DTMapGenerator.h.

{}

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

private

Definition at line 167 of file DTMapGenerator.cc.

References gather_cfg::cout.

Referenced by analyze().

                                                                                                   {
  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;
}

void DTMapGenerator::endJob ( void  )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 33 of file DTMapGenerator.cc.

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

Referenced by analyze(), and 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;
            }
          }
          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

std::string DTMapGenerator::inputMapName

private

Definition at line 43 of file DTMapGenerator.h.

std::string DTMapGenerator::outputMapName

private

Definition at line 41 of file DTMapGenerator.h.

int DTMapGenerator::rosType

private

Definition at line 45 of file DTMapGenerator.h.