PixelToLNKAssociateFromAscii Class Reference

#include <PixelToLNKAssociateFromAscii.h>

Inheritance diagram for PixelToLNKAssociateFromAscii:

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

Public Member Functions
virtual const CablingRocId *	operator() (const DetectorRocId &roc) const
	LNK id for module. More...
	PixelToLNKAssociateFromAscii (const std::string &fileName, const bool phase1=false)
virtual std::string	version () const
	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 23 of file PixelToLNKAssociateFromAscii.h.

typedef PixelToFEDAssociate::DetectorRocId PixelToLNKAssociateFromAscii::DetectorRocId

Definition at line 24 of file PixelToLNKAssociateFromAscii.h.

typedef TRange<int> PixelToLNKAssociateFromAscii::Range

private

Definition at line 34 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 StGenEvtProducer_cfi::init.

                                                                                              {
  phase1_ = phase;
  init(fn);

}

Member Function Documentation

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

private

Definition at line 135 of file PixelToLNKAssociateFromAscii.cc.

References Exception, l1t::stage2::layer2::fedId, PixelToFEDAssociate::CablingRocId::fedId, PixelBarrelName::isHalfModule(), PixelToFEDAssociate::CablingRocId::linkId, TRange< T >::max(), PixelBarrelName::mI, PixelEndcapName::mI, TRange< T >::min(), PixelEndcapName::mO, PixelBarrelName::mO, PixelToFEDAssociate::DetectorRocId::module, python.rootplot.argparse::module, PixelEndcapName::name(), dataset::name, AlCaHLTBitMon_ParallelJobs::p, PixelPannelType::p2x8, PixelPannelType::p3L, PixelPannelType::p3R, PixelPannelType::p4L, PixelPannelType::p4R, PixelEndcapName::pI, PixelBarrelName::pI, PixelBarrelName::pO, PixelEndcapName::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; 
     for (int rocDetId=rocDetIds.min(); rocDetId <= rocDetIds.max(); rocDetId++) {
       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));
     } 
  }

  // 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_); 
       for (int rocDetId=rocDetIds.min(); rocDetId <= rocDetIds.max(); rocDetId++) {
     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 
       
     } // end of type

  }
}

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

private

initialisatin (read file)

Definition at line 36 of file PixelToLNKAssociateFromAscii.cc.

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

                                                               {

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

  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

virtual

LNK id for module.

Reimplemented from PixelToFEDAssociate.

Definition at line 23 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 0;
}

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

private

Definition at line 431 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

virtual

version

Implements PixelToFEDAssociate.

Definition at line 19 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 42 of file PixelToLNKAssociateFromAscii.h.

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

private

Definition at line 40 of file PixelToLNKAssociateFromAscii.h.

std::string PixelToLNKAssociateFromAscii::theVersion

private

Definition at line 39 of file PixelToLNKAssociateFromAscii.h.