PixelToLNKAssociateFromAscii Class Reference

#include <PixelToLNKAssociateFromAscii.h>

Inheritance diagram for PixelToLNKAssociateFromAscii:

Public Types
typedef PixelToFEDAssociate::CablingRocId	CablingRocId
typedef PixelToFEDAssociate::DetectorRocId	DetectorRocId

Public Member Functions
const CablingRocId *	operator() (const DetectorRocId &roc) const override
	LNK id for module. More...
	PixelToLNKAssociateFromAscii (const std::string &fileName, const bool phase1=false)
std::string	version () const override
	version More...
Public Member Functions inherited from PixelToFEDAssociate
virtual int	operator() (const PixelModuleName &) const
	FED id for module. More...
virtual	~PixelToFEDAssociate ()

Private Types
typedef TRange< int >	Range

Private Member Functions
void	addConnections (int fedId, int linkId, std::string module, Range rocDetIds)
void	init (const std::string &fileName)
	initialisatin (read file) More...
Range	readRange (const std::string &) const

Private Attributes
bool	phase1_
std::vector< std::pair< DetectorRocId, CablingRocId > >	theConnection
std::string	theVersion

Detailed Description

Check to which FED pixel module belongs to. The associacions are read from the datafile

Definition at line 20 of file PixelToLNKAssociateFromAscii.h.

Member Typedef Documentation

typedef PixelToFEDAssociate::CablingRocId PixelToLNKAssociateFromAscii::CablingRocId

Definition at line 22 of file PixelToLNKAssociateFromAscii.h.

typedef PixelToFEDAssociate::DetectorRocId PixelToLNKAssociateFromAscii::DetectorRocId

Definition at line 23 of file PixelToLNKAssociateFromAscii.h.

typedef TRange<int> PixelToLNKAssociateFromAscii::Range

private

Definition at line 33 of file PixelToLNKAssociateFromAscii.h.

Constructor & Destructor Documentation

PixelToLNKAssociateFromAscii::PixelToLNKAssociateFromAscii	(	const std::string &	fileName,
		const bool	phase1 = false
	)

Definition at line 14 of file PixelToLNKAssociateFromAscii.cc.

References init, and hcalRecAlgoESProd_cfi::phase.

                                                                                             {
  phase1_ = phase;
  init(fn);
}

Member Function Documentation

void PixelToLNKAssociateFromAscii::addConnections ( int  fedId, int  linkId, std::string  module, Range  rocDetIds  )

private

Definition at line 126 of file PixelToLNKAssociateFromAscii.cc.

References Exception, l1t::stage2::layer2::fedId, PixelToFEDAssociate::CablingRocId::fedId, PixelBarrelName::isHalfModule(), PixelToFEDAssociate::CablingRocId::linkId, TRange< T >::max(), PixelEndcapName::mI, PixelBarrelName::mI, TRange< T >::min(), PixelBarrelName::mO, PixelEndcapName::mO, PixelToFEDAssociate::DetectorRocId::module, PixelEndcapName::name(), dataset::name, AlCaHLTBitMon_ParallelJobs::p, PixelPannelType::p2x8, PixelPannelType::p3L, PixelPannelType::p3R, PixelPannelType::p4L, PixelPannelType::p4R, PixelEndcapName::pI, PixelBarrelName::pI, PixelEndcapName::pO, PixelBarrelName::pO, relativeConstraints::ring, PixelToFEDAssociate::DetectorRocId::rocDetId, and PixelToFEDAssociate::CablingRocId::rocLinkId.

                                                                                                          {
  string::size_type pos;

  // check for Barrel modules
  pos = module.find("BPix");
  if (pos != string::npos) {
    string module0 = module;
    // strip the trailing spaces
    string::size_type p = module0.find(" ");
    //string::size_type p1 = module0.find_first_of(" ");
    //string::size_type p2 = module0.find_last_not_of(" ");
    //cout<<p<<" "<<p1<<" "<<p2<<endl;
    if (p != string::npos)
      module0 = module0.substr(0, p);
    PixelBarrelName* name = new PixelBarrelName(module0, phase1_);

    // shell
    string strP = module.substr(pos + 6, 2);
    PixelBarrelName::Shell part;
    if (strP == "mO")
      part = PixelBarrelName::mO;
    else if (strP == "mI")
      part = PixelBarrelName::mI;
    else if (strP == "pO")
      part = PixelBarrelName::pO;
    else
      part = PixelBarrelName::pI;

    // // all this can be skipped -----------------------------------
    // module = module.substr(pos+9);
    // // sector
    // pos = module.find("_");
    // if (pos ==  string::npos) throw cms::Exception("problem with sector formatting");
    // // int sector = atoi( module.substr(3,pos-3).c_str());
    // module = module.substr(pos+1);
    // // layer
    // pos = module.find("_");
    // if (pos ==  string::npos) throw cms::Exception("problem with layer formatting");
    // int layer = atoi( module.substr(3,pos-3).c_str());
    // module = module.substr(pos+1);
    // // ladder
    // pos = module.find("_");
    // if (pos ==  string::npos) throw cms::Exception("problem with ladder formatting");
    // int ladder = atoi( module.substr(3,pos-3).c_str());
    // module = module.substr(pos+1);
    // // z-module
    // int zmodule = atoi( module.substr(3,pos-3).c_str());
    // // place modules in connections
    // PixelBarrelName * name0 = new PixelBarrelName(part, layer, zmodule, ladder, phase1_);
    // if(name->name() != module0) cout<<" wrong name "<<fedId<<" "<<linkId<<" "
    //                   <<module0<<" "<<name->name()<<" "<<name0->name()<<endl;
    // if(name->name() != name0->name()) cout<<" wrong name "<<fedId<<" "<<linkId<<" "
    //                   <<module0<<" "<<name->name()<<" "<<name0->name()<<endl;
    // //edm::LogInfo(" module ")<<fedId<<" "<<linkId<<" "<<module0<<" "
    // //                <<name0->name()<<" "<<rocDetIds.max()<<endl;
    // // until here

    int rocLnkId = 0;
    bool loopExecuted = false;
    for (int rocDetId = rocDetIds.min(); rocDetId <= rocDetIds.max(); rocDetId++) {
      loopExecuted = true;
      rocLnkId++;
      DetectorRocId detectorRocId;
      //detectorRocId.module = name0;
      detectorRocId.module = name;
      detectorRocId.rocDetId = rocDetId;

      CablingRocId cablingRocId;
      cablingRocId.fedId = fedId;
      cablingRocId.linkId = linkId;
      cablingRocId.rocLinkId = rocLnkId;
      // fix for type-B modules in barrel
      edm::LogInfo(" roc ") << rocDetId << " " << rocLnkId << " " << name->isHalfModule() << endl;
      if (name->isHalfModule() && (rocDetIds.min() > 7) &&
          (part == PixelBarrelName::mO || part == PixelBarrelName::mI)) {
        //cablingRocId.rocLinkId = 9-rocLnkId;
        // rocDetId=8,...,15
        edm::LogInfo(" special for half modules ");
        cablingRocId.rocLinkId = rocLnkId;      // 1...8    19/11/08 d.k.
        detectorRocId.rocDetId = rocDetId - 8;  // 0...7
      }
      theConnection.push_back(make_pair(detectorRocId, cablingRocId));
    }
    if (!loopExecuted)
      delete name;
  }

  // check for endcap modules
  pos = module.find("FPix");
  if (pos != string::npos) {
    string strH = module.substr(pos + 6, 2);
    PixelEndcapName::HalfCylinder part;
    if (strH == "mO")
      part = PixelEndcapName::mO;
    else if (strH == "mI")
      part = PixelEndcapName::mI;
    else if (strH == "pO")
      part = PixelEndcapName::pO;
    else
      part = PixelEndcapName::pI;
    module = module.substr(pos + 9);

    // disk
    pos = module.find("_");
    if (pos == string::npos)
      throw cms::Exception("problem with disk formatting");
    int disk = atoi(module.substr(1, pos - 1).c_str());
    module = module.substr(pos + 1);

    // blade
    pos = module.find("_");
    if (pos == string::npos)
      throw cms::Exception("problem with blade formatting");
    int blade = atoi(module.substr(3, pos - 3).c_str());
    module = module.substr(pos + 1);

    //pannel
    pos = module.find("_");
    if (pos == string::npos)
      throw cms::Exception("problem with pannel formatting");
    int pannel = atoi(module.substr(3, pos - 3).c_str());
    module = module.substr(pos + 1);

    // plaquete
    //     pos = module.find("_");
    //     if (pos ==  string::npos) throw cms::Exception("problem with plaquette formatting");
    //     int plaq = atoi( module.substr(3,pos-3).c_str());

    int ring = 1;  // preset to 1 so it is ok for pilot blades
    // pannel type
    PixelPannelType::PannelType pannelType;
    if (phase1_) {
      pannelType = PixelPannelType::p2x8;  // only 1 type for phase1

      // this is not really needed, just for testing
      // ring
      pos = module.find("RNG");
      if (pos == string::npos)
        throw cms::Exception("problem with ring  formatting");
      ring = atoi(module.substr(pos + 3, 1).c_str());  //
      //cout<<" ring "<<ring<<" "<<module<<endl;

    } else {  // phase0
      pos = module.find("TYP:");
      if (pos == string::npos)
        throw cms::Exception("problem with pannel type formatting");
      string strT = module.substr(pos + 5, 3);
      string strT4 = module.substr(pos + 5, 4);
      ring = 1;
      if (strT == "P3R")
        pannelType = PixelPannelType::p3R;
      else if (strT == "P3L")
        pannelType = PixelPannelType::p3L;
      else if (strT == "P4R")
        pannelType = PixelPannelType::p4R;
      else if (strT == "P4L")
        pannelType = PixelPannelType::p4L;
      else if (strT4 == "P2X8")
        pannelType = PixelPannelType::p2x8;  // for pilot blades
      else
        throw cms::Exception("problem with pannel type formatting (unrecoginzed word)");
    }

    // Cabling accoring to the panle type
    if (pannelType == PixelPannelType::p4L) {
      //     cout <<"----------- p4L"<<endl;
      int rocLnkId = 0;
      for (int plaq = 1; plaq <= 4; plaq++) {
        Range rocs;
        int firstRoc = 0;
        int step = 0;
        if (plaq == 1) {
          rocs = Range(0, 1);
          firstRoc = 1;
          step = -1;
        }
        if (plaq == 2) {
          rocs = Range(0, 5);
          firstRoc = 0;
          step = +1;
        }
        if (plaq == 3) {
          rocs = Range(0, 7);
          firstRoc = 0;
          step = +1;
        }
        if (plaq == 4) {
          rocs = Range(0, 4);
          firstRoc = 0;
          step = +1;
        }
        PixelEndcapName* name = new PixelEndcapName(part, disk, blade, pannel, plaq);
        for (int iroc = rocs.min(); iroc <= rocs.max(); iroc++) {
          rocLnkId++;
          int rocDetId = firstRoc + step * iroc;

          DetectorRocId detectorRocId;
          //detectorRocId.module = new PixelEndcapName(part,disk,blade,pannel,plaq);
          detectorRocId.module = name;
          detectorRocId.rocDetId = rocDetId;

          CablingRocId cablingRocId;
          cablingRocId.fedId = fedId;
          cablingRocId.linkId = linkId;
          cablingRocId.rocLinkId = rocLnkId;

          theConnection.push_back(make_pair(detectorRocId, cablingRocId));
          //         cout <<"PLAQ:"<<plaq<<" rocDetId: "<<rocDetId<<" rocLnkId:"<<rocLnkId<<endl;
        }
      }
    } else if (pannelType == PixelPannelType::p4R) {
      //     cout <<"----------- p4R"<<endl;
      int rocLnkId = 0;
      for (int plaq = 4; plaq >= 1; plaq--) {
        Range rocs;
        int firstRoc = 0;
        int step = 0;
        if (plaq == 1) {
          rocs = Range(0, 1);
          firstRoc = 1;
          step = -1;
        }
        if (plaq == 2) {
          rocs = Range(0, 5);
          firstRoc = 3;
          step = +1;
        }
        if (plaq == 3) {
          rocs = Range(0, 7);
          firstRoc = 4;
          step = +1;
        }
        if (plaq == 4) {
          rocs = Range(0, 4);
          firstRoc = 0;
          step = +1;
        }
        PixelEndcapName* name = new PixelEndcapName(part, disk, blade, pannel, plaq);
        for (int iroc = rocs.min(); iroc - rocs.max() <= 0; iroc++) {
          rocLnkId++;
          int rocDetId = firstRoc + step * iroc;
          if (rocDetId > rocs.max())
            rocDetId = (rocDetId - 1) % rocs.max();

          DetectorRocId detectorRocId;
          //detectorRocId.module = new PixelEndcapName(part,disk,blade,pannel,plaq);
          detectorRocId.module = name;
          detectorRocId.rocDetId = rocDetId;
          CablingRocId cablingRocId;
          cablingRocId.fedId = fedId;
          cablingRocId.linkId = linkId;
          cablingRocId.rocLinkId = rocLnkId;

          theConnection.push_back(make_pair(detectorRocId, cablingRocId));
          //         cout <<"PLAQ:"<<plaq<<" rocDetId: "<<rocDetId<<" rocLnkId:"<<rocLnkId<<endl;
        }
      }
    } else if (pannelType == PixelPannelType::p3L) {
      //     cout <<"----------- p3L"<<endl;
      int rocLnkId = 0;
      for (int plaq = 1; plaq <= 3; plaq++) {
        Range rocs;
        int firstRoc = 0;
        int step = 0;
        if (plaq == 1) {
          rocs = Range(0, 5);
          firstRoc = 0;
          step = 1;
        }
        if (plaq == 2) {
          rocs = Range(0, 7);
          firstRoc = 0;
          step = 1;
        }
        if (plaq == 3) {
          rocs = Range(0, 9);
          firstRoc = 0;
          step = 1;
        }
        PixelEndcapName* name = new PixelEndcapName(part, disk, blade, pannel, plaq);
        for (int iroc = rocs.min(); iroc <= rocs.max(); iroc++) {
          rocLnkId++;
          int rocDetId = firstRoc + step * iroc;

          DetectorRocId detectorRocId;
          detectorRocId.module = name;
          detectorRocId.rocDetId = rocDetId;

          CablingRocId cablingRocId;
          cablingRocId.fedId = fedId;
          cablingRocId.linkId = linkId;
          cablingRocId.rocLinkId = rocLnkId;

          theConnection.push_back(make_pair(detectorRocId, cablingRocId));
          //         cout <<"PLAQ:"<<plaq<<" rocDetId: "<<rocDetId<<" rocLnkId:"<<rocLnkId<<endl;
        }
      }
    } else if (pannelType == PixelPannelType::p3R) {
      //     cout <<"----------- p3R"<<endl;
      int rocLnkId = 0;
      for (int plaq = 3; plaq >= 1; plaq--) {
        Range rocs;
        int firstRoc = 0;
        int step = 0;
        if (plaq == 1) {
          rocs = Range(0, 5);
          firstRoc = 3;
          step = 1;
        }
        if (plaq == 2) {
          rocs = Range(0, 7);
          firstRoc = 4;
          step = 1;
        }
        if (plaq == 3) {
          rocs = Range(0, 9);
          firstRoc = 5;
          step = 1;
        }
        PixelEndcapName* name = new PixelEndcapName(part, disk, blade, pannel, plaq);
        for (int iroc = rocs.min(); iroc <= rocs.max(); iroc++) {
          rocLnkId++;
          int rocDetId = firstRoc + step * iroc;
          if (rocDetId > rocs.max())
            rocDetId = (rocDetId - 1) % rocs.max();

          DetectorRocId detectorRocId;
          detectorRocId.module = name;
          detectorRocId.rocDetId = rocDetId;

          CablingRocId cablingRocId;
          cablingRocId.fedId = fedId;
          cablingRocId.linkId = linkId;
          cablingRocId.rocLinkId = rocLnkId;

          theConnection.push_back(make_pair(detectorRocId, cablingRocId));
          //         cout <<"PLAQ:"<<plaq<<" rocDetId: "<<rocDetId<<" rocLnkId:"<<rocLnkId<<endl;
        }  // for
      }    // for

    } else if (pannelType == PixelPannelType::p2x8) {  // phase-1 blades
      //       cout <<"----------- p2x8"<<endl;
      int rocLnkId = 0;
      //       Range rocs = Range(0, 15);
      //       for (int rocDetId=rocs.min(); rocDetId <= rocs.max(); rocDetId++) {
      PixelEndcapName* name = new PixelEndcapName(part, disk, blade, pannel, ring, phase1_);
      bool loopExecuted = false;
      for (int rocDetId = rocDetIds.min(); rocDetId <= rocDetIds.max(); rocDetId++) {
        loopExecuted = true;
        rocLnkId++;
        DetectorRocId detectorRocId;
        detectorRocId.module = name;
        detectorRocId.rocDetId = rocDetId;
        CablingRocId cablingRocId;
        cablingRocId.fedId = fedId;
        cablingRocId.linkId = linkId;
        cablingRocId.rocLinkId = rocLnkId;
        theConnection.push_back(make_pair(detectorRocId, cablingRocId));
        edm::LogInfo("PixelToLNKAssociateFromAscii FPix ")
            << " rocDetId: " << rocDetId << " rocLnkId:" << rocLnkId << " fedId = " << fedId << " linkId = " << linkId
            << " name = " << name->name();
        // cout << " rocDetId: " << rocDetId
        //      << " rocLnkId:" << rocLnkId
        //      << " fedId = " << fedId
        //      << " linkId = " << linkId
        //      << " name = " << name0->name()
        //      << endl;
      }  // end for
      if (!loopExecuted) {
        delete name;
      }

    }  // end of type
  }
}

void PixelToLNKAssociateFromAscii::init ( const std::string &  fileName )

private

initialisatin (read file)

Definition at line 32 of file PixelToLNKAssociateFromAscii.cc.

References cppFunctionSkipper::exception, l1t::stage2::layer2::fedId, FrontierConditions_GlobalTag_cff::file, mps_splice::line, and str.

                                                              {
  edm::LogInfo(" init, input file: ") << cfg_name;

  std::ifstream file(cfg_name.c_str());
  if (!file) {
    edm::LogError(" ** PixelToLNKAssociateFromAscii,init ** ") << " cant open data file: " << cfg_name;
    return;
  } else {
    edm::LogInfo("PixelToLNKAssociateFromAscii, read data from: ") << cfg_name;
  }

  string line;
  int fedId = -1;
  int linkId = -1;

  try {
    while (getline(file, line)) {
      //
      // treat # lines
      //
      string::size_type pos = line.find("#");
      if (pos != string::npos)
        line = line.erase(pos);

      string::size_type posF = line.find("FED:");
      string::size_type posL = line.find("LNK:");
      string::size_type posM = line.find("MOD:");
      string::size_type posR = line.find("ROC:");

      edm::LogInfo("") << " read line: " << line;

      //
      // treat version lines, reset date
      //
      if (line.compare(0, 3, "VER") == 0) {
        edm::LogInfo("version: ") << line;
        theVersion = line;
      }

      //
      // fed id line
      //
      else if (posF != string::npos) {
        line = line.substr(posF + 4);
        fedId = atoi(line.c_str());
      }

      //
      // link id linke
      //
      else if (posL != string::npos) {
        string srtL = line.substr(posL + 4);
        linkId = atoi(srtL.c_str());
      }

      //
      // module description
      //
      if (posM != string::npos) {
        if (posR != string::npos) {
          string strM = line.substr(posM + 4, posR - posM - 5);
          string::size_type pos = strM.find(" ");
          if (pos != string::npos)
            strM = strM.substr(pos + 1);
          string strR = line.substr(posR + 4);
          Range range = readRange(strR);
          //cout<<" range find "<<strR<<" "<<strR.size()<<" "<<range.min()<<" "<<range.max()<<endl;
          addConnections(fedId, linkId, strM, range);
        } else {
          string strM = line.substr(posM + 4);
          string::size_type pos = strM.find(" ");
          if (pos != string::npos)
            strM = strM.substr(pos + 1);
          addConnections(fedId, linkId, strM, Range(0, 0));
        }
      }
    }
  } catch (exception& err) {
    edm::LogError("**PixelToLNKAssociateFromAscii**  exception") << err.what();
  }

  //
  // for debug
  //
  std::ostringstream str;
  str << " **PixelToLNKAssociateFromAscii ** CONNECTIONS: " << endl;
  typedef vector<pair<DetectorRocId, CablingRocId> >::const_iterator ICON;
  for (ICON ic = theConnection.begin(); ic != theConnection.end(); ic++) {
    str << (*ic).first.module->name() << ", rocDetId=" << (*ic).first.rocDetId << ", fedId=" << ic->second.fedId
        << ", linkId=" << ic->second.linkId << ", rocLinkId=" << ic->second.rocLinkId << endl;
  }
  edm::LogInfo("PixelToLNKAssociateFromAscii") << str.str();
}

const PixelToLNKAssociateFromAscii::CablingRocId * PixelToLNKAssociateFromAscii::operator() ( const DetectorRocId &  roc ) const

overridevirtual

LNK id for module.

Reimplemented from PixelToFEDAssociate.

Definition at line 20 of file PixelToLNKAssociateFromAscii.cc.

References PixelToFEDAssociate::DetectorRocId::module, and PixelToFEDAssociate::DetectorRocId::rocDetId.

                                                                {
  typedef std::vector<std::pair<DetectorRocId, CablingRocId> >::const_iterator IM;
  for (IM im = theConnection.begin(); im != theConnection.end(); im++) {
    if ((*(im->first.module) == *roc.module) && (im->first.rocDetId == roc.rocDetId)) {
      return &(im->second);
    }
  }
  return nullptr;
}

PixelToLNKAssociateFromAscii::Range PixelToLNKAssociateFromAscii::readRange ( const std::string &  ) const

private

Definition at line 502 of file PixelToLNKAssociateFromAscii.cc.

References AlCaHLTBitMon_ParallelJobs::p.

                                                                                               {
  //cout<<l<<" in range "<<l.size()<<endl;
  string l1, l2;
  int i1 = -1, i2 = -1;
  int len = l.size();
  //for(int i=0; i<len;i++) {
  // cout<<i<<" "<<l[i]<<endl;
  //}
  string::size_type p = l.find(",");
  if (p != string::npos) {
    //cout<<p<<" "<<len<<endl;
    l1 = l.substr(0, p - 1 + 1);
    l2 = l.substr(p + 1, len - 1 - p);
    i1 = stoi(l1);
    i2 = stoi(l2);
    //cout<<l1<<" "<<l2<<" "<<i1<<" "<<i2<<endl;
  }

  return Range(i1, i2);

  // this method is very stupid it relies on a space being present after the last number!
  // exchange with string opertaions (above)
  //bool first = true;
  //int num1 = -1;
  //int num2 = -1;
  // const char * line = l.c_str();
  // int i=0;
  // while (line) {
  //   i++;
  //   char * evp = 0;
  //   int num = strtol(line, &evp, 10);
  //   //cout<<i<<" "<<num<<" "<<evp<<" "<<line<<endl;
  //   //{ stringstream s; s<<"read from line: "; s<<num; LogTrace("") << s.str(); }
  //   if (evp != line) {
  //     line = evp +1;
  //     //cout<<i<<" "<<num<<" "<<evp<<" "<<line<<endl;
  //     if (first) { num1 = num; first = false;}
  //     num2 = num;
  //     //cout<<" not first "<<num2<<endl;
  //   } else line = 0;
  // }
  // if (first) {
  //   string s = "** PixelToLNKAssociateFromAscii, read data, cant intrpret: " ;
  //   edm::LogInfo(s) << endl
  //            << l << endl
  //             <<"=====> send exception " << endl;
  //   s += l;
  //   throw cms::Exception(s);
  // }
  //if(i1!=num1) cout<<" something wrong with min range "<<i1<<" "<<num1<<endl;
  //if(!phase1_ && (i2!=num2)) cout<<" something wrong with max range "<<i2<<" "<<num2<<endl;
  //cout<<" min max "<<num1<<" "<<num2<<endl;
  //return Range(num1,num2);
}

std::string PixelToLNKAssociateFromAscii::version ( ) const

overridevirtual

version

Implements PixelToFEDAssociate.

Definition at line 18 of file PixelToLNKAssociateFromAscii.cc.

Referenced by validation.Sample::datasetpattern(), and validation.Sample::filename().

{ return theVersion; }

Member Data Documentation

bool PixelToLNKAssociateFromAscii::phase1_

private

Definition at line 41 of file PixelToLNKAssociateFromAscii.h.

std::vector<std::pair<DetectorRocId, CablingRocId> > PixelToLNKAssociateFromAscii::theConnection

private

Definition at line 39 of file PixelToLNKAssociateFromAscii.h.

std::string PixelToLNKAssociateFromAscii::theVersion

private

Definition at line 38 of file PixelToLNKAssociateFromAscii.h.