L1TMuonBarrelKalmanStubProcessor.cc

Go to the documentation of this file.

#include "L1Trigger/L1TMuonBarrel/interface/L1TMuonBarrelKalmanStubProcessor.h"
#include <cmath>
#include "CondFormats/L1TObjects/interface/L1MuDTTFParameters.h"
#include "CondFormats/DataRecord/interface/L1MuDTTFParametersRcd.h"
#include "CondFormats/L1TObjects/interface/L1MuDTTFMasks.h"
#include "CondFormats/DataRecord/interface/L1MuDTTFMasksRcd.h"

#include <iostream> 
#include <string> 
#include <sstream> 

L1TMuonBarrelKalmanStubProcessor::L1TMuonBarrelKalmanStubProcessor():
  minPhiQuality_(0),
  minBX_(-3),
  maxBX_(3)
{
 
} 



L1TMuonBarrelKalmanStubProcessor::L1TMuonBarrelKalmanStubProcessor(const edm::ParameterSet& iConfig):
  minPhiQuality_(iConfig.getParameter<int>("minPhiQuality")),
  minBX_(iConfig.getParameter<int>("minBX")),
  maxBX_(iConfig.getParameter<int>("maxBX")),
  eta1_(iConfig.getParameter<std::vector<int> >("cotTheta_1")),
  eta2_(iConfig.getParameter<std::vector<int> >("cotTheta_2")),
  eta3_(iConfig.getParameter<std::vector<int> >("cotTheta_3")),
  disableMasks_(iConfig.getParameter<bool>("disableMasks")),
  verbose_(iConfig.getParameter<int>("verbose"))
{

} 



L1TMuonBarrelKalmanStubProcessor::~L1TMuonBarrelKalmanStubProcessor() {}



bool L1TMuonBarrelKalmanStubProcessor::isGoodPhiStub(const L1MuDTChambPhDigi * stub) {
  if (stub->code()<minPhiQuality_)
    return false;
  return true;
}




L1MuKBMTCombinedStub 
L1TMuonBarrelKalmanStubProcessor::buildStub(const L1MuDTChambPhDigi& phiS,const L1MuDTChambThDigi* etaS) {
  int wheel = phiS.whNum();
  int sector = phiS.scNum();
  int station = phiS.stNum();
  int phi = phiS.phi();
  int phiB = phiS.phiB();
  bool tag = (phiS.Ts2Tag()==1);
  int bx=phiS.bxNum();
  int quality=phiS.code();



  //Now full eta
  int qeta1=0;
  int qeta2=0;
  int eta1=255;
  int eta2=255; 


  bool hasEta=false;
  for (uint i=0;i<7;++i) {
    if (etaS->position(i)==0)
      continue;
    if (!hasEta) {
      hasEta=true;
      eta1=calculateEta(i,etaS->whNum(),etaS->scNum(),etaS->stNum());
      if (etaS->quality(i)==1)
    qeta1=2;
      else
    qeta1=1;
    }
    else {
      eta2=calculateEta(i,etaS->whNum(),etaS->scNum(),etaS->stNum());
      if (etaS->quality(i)==1)
    qeta2=2;
      else
    qeta2=1;
    }
  }
  L1MuKBMTCombinedStub stub(wheel,sector,station,phi,phiB,tag,
                bx,quality,eta1,eta2,qeta1,qeta2);
  
  return stub;

  }




L1MuKBMTCombinedStub 
L1TMuonBarrelKalmanStubProcessor::buildStubNoEta(const L1MuDTChambPhDigi& phiS) {
  int wheel = phiS.whNum();
  int sector = phiS.scNum();
  int station = phiS.stNum();
  int phi = phiS.phi();
  int phiB = phiS.phiB();
  bool tag = (phiS.Ts2Tag()==1);
  int bx=phiS.bxNum();
  int quality=phiS.code();


  //Now full eta
  int qeta1=0;
  int qeta2=0;
  int eta1=7;
  int eta2=7; 
  L1MuKBMTCombinedStub stub(wheel,sector,station,phi,phiB,tag,
                bx,quality,eta1,eta2,qeta1,qeta2);

  return stub;

}





L1MuKBMTCombinedStubCollection 
L1TMuonBarrelKalmanStubProcessor::makeStubs(const L1MuDTChambPhContainer* phiContainer,const L1MuDTChambThContainer* etaContainer,const L1TMuonBarrelParams& params) {


  //get the masks from th standard BMTF setup!
  //    const L1TMuonBarrelParamsRcd& bmtfParamsRcd = setup.get<L1TMuonBarrelParamsRcd>();
  //  bmtfParamsRcd.get(bmtfParamsHandle);
  //  const L1TMuonBarrelParams& bmtfParams = *bmtfParamsHandle.product();
  //  masks_ =  bmtfParams.l1mudttfmasks;


  //get the masks!
  L1MuDTTFMasks msks = params.l1mudttfmasks;


  L1MuKBMTCombinedStubCollection  out;
  for (int bx=minBX_;bx<=maxBX_;bx++) {
    for (int wheel=-2;wheel<=2;wheel++) {
      for (uint sector=0;sector<12;sector++) {
    for (uint station=1;station<5;station++) {

      //Have to cook up something for the fact that KMTF doesnt use 2 SP at whel=0
      int lwheel1;
      int lwheel2;
      if (wheel<0) {
        lwheel1=wheel-1;
        lwheel2=wheel-1;
      }
      else if (wheel>0) {
        lwheel1=wheel+1;
        lwheel2=wheel+1;
      }
      else {
        lwheel1=-1;
        lwheel2=+1;
      }
        
      bool phiMask=false;
      bool etaMask=false;
      if (station==1) {
        phiMask = msks.get_inrec_chdis_st1(lwheel1, sector) |msks.get_inrec_chdis_st1(lwheel2, sector);  
        etaMask = msks.get_etsoc_chdis_st1(lwheel1, sector) |msks.get_etsoc_chdis_st1(lwheel2, sector);         
      }
      if (station==2) {
        phiMask = msks.get_inrec_chdis_st2(lwheel1, sector) |msks.get_inrec_chdis_st2(lwheel2, sector);  
        etaMask = msks.get_etsoc_chdis_st2(lwheel1, sector) |msks.get_etsoc_chdis_st2(lwheel2, sector);  

      }
      if (station==3) {
        phiMask = msks.get_inrec_chdis_st3(lwheel1, sector) |msks.get_inrec_chdis_st3(lwheel2, sector);  
        etaMask = msks.get_etsoc_chdis_st3(lwheel1, sector) |msks.get_etsoc_chdis_st3(lwheel2, sector);  
      }
      if (station==4) {
        phiMask = msks.get_inrec_chdis_st4(lwheel1, sector) |msks.get_inrec_chdis_st4(lwheel2, sector);  
      }


      if (disableMasks_)
        {
          phiMask=false;
          etaMask=false;
        }


      bool hasEta=false;
      L1MuDTChambThDigi const*  tseta   = etaContainer->chThetaSegm(wheel,station,sector,bx);
      if (tseta && (!etaMask)) {
        hasEta=true;
      }

      //      if (abs(wheel)==2 && station==1)
      //        continue;
        

      L1MuDTChambPhDigi const* high = phiContainer->chPhiSegm1(wheel,station,sector,bx);
      if (high && (!phiMask)) {
        if (high->code()>=minPhiQuality_) {
          const L1MuDTChambPhDigi& stubPhi = *high;
          if (hasEta) {
        out.push_back(buildStub(stubPhi,tseta));
          }
          else {
        out.push_back(buildStubNoEta(stubPhi));
          }
        }
      }

      L1MuDTChambPhDigi const* low = phiContainer->chPhiSegm2(wheel,station,sector,bx-1);
      if (low && ! (phiMask)) {
        if (low->code()>=minPhiQuality_) {
          const L1MuDTChambPhDigi& stubPhi = *low;
          if (hasEta) {
        out.push_back(buildStub(stubPhi,tseta));
          }
          else {
        out.push_back(buildStubNoEta(stubPhi));
          }
        }
      }
    }
      }
    }
  }







  return out;
}



int L1TMuonBarrelKalmanStubProcessor::calculateEta(uint i, int wheel,uint sector,uint station) {
  int eta=0;
  if (wheel>0) {
    eta=7*wheel+3-i;
      }
  else if (wheel<0) {
    eta=7*wheel+i-3;
  }
  else {
    if (sector==0 || sector==3 ||sector==4 ||sector==7 ||sector==8 ||sector==11)
      eta=i-3;
    else
      eta=3-i;
  }


  if (station==1)
    eta=-eta1_[eta+17];
  else if (station==2)
    eta=-eta2_[eta+17];
  else 
    eta=-eta3_[eta+17];



  return eta;



}



void  L1TMuonBarrelKalmanStubProcessor::makeInputPattern(const L1MuDTChambPhContainer* phiContainer,const L1MuDTChambThContainer* etaContainer,int sector) {

  int previousSector;
  int nextSector;

  if (sector==11)
    nextSector=0;
  else
    nextSector=sector+1;

  if (sector==0)
    previousSector=11;
  else
    previousSector=sector-1;


  ostringstream os;
  os<<"I "<<sector<<" ";


  bool hasStub=false;

  //start by previous sector
  for (int wheel=-2;wheel<3;++wheel) {
    const L1MuDTChambPhDigi* seg1 = phiContainer->chPhiSegm1(wheel,1,previousSector,0);
    if (seg1 && seg1->phi()>111) {
      os<<seg1->phi()-2144<<" "<<seg1->phiB()<<" "<<seg1->code()<<" 1 4 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }
    const L1MuDTChambPhDigi* seg2 = phiContainer->chPhiSegm1(wheel,2,previousSector,0);
    if (seg2 && seg2->phi()>111) {
      os<<seg2->phi()-2144<<" "<<seg2->phiB()<<" "<<seg2->code()<<" 1 4 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }
    const L1MuDTChambPhDigi* seg3 = phiContainer->chPhiSegm1(wheel,3,previousSector,0);
    if (seg3 && seg3->phi()>111) {
      os<<seg3->phi()-2144<<" "<<seg3->phiB()<<" "<<seg3->code()<<" 1 4 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }
    const L1MuDTChambPhDigi* seg4 = phiContainer->chPhiSegm1(wheel,4,previousSector,0);
    if (seg4 && seg4->phi()>111) {
      os<<seg4->phi()-2144<<" "<<seg4->phiB()<<" "<<seg4->code()<<" 1 4 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }

    const L1MuDTChambThDigi*  eta1   = etaContainer->chThetaSegm(wheel,1,previousSector,0);
    if (eta1) {
      int etaPos=eta1->position(0)+(eta1->position(1)<<1) +(eta1->position(2)<<2) +(eta1->position(3)<<3) +(eta1->position(4)<<4)+(eta1->position(5)<<5)+(eta1->position(6)<<6);
      os<< etaPos<<" ";
    }
    else {
      os<< " 0 ";
    }
    const L1MuDTChambThDigi*  eta2   = etaContainer->chThetaSegm(wheel,2,previousSector,0);
    if (eta2) {
      int etaPos=eta2->position(0)+(eta2->position(1)<<1) +(eta2->position(2)<<2) +(eta2->position(3)<<3) +(eta2->position(4)<<4)+(eta2->position(5)<<5)+(eta2->position(6)<<6);
      os<< etaPos<<" ";
    }
    else {
      os<< " 0 ";
    }

    const L1MuDTChambThDigi*  eta3   = etaContainer->chThetaSegm(wheel,3,previousSector,0);
    if (eta3) {
      int etaPos=eta3->position(0)+(eta3->position(1)<<1) +(eta3->position(2)<<2) +(eta3->position(3)<<3) +(eta3->position(4)<<4)+(eta3->position(5)<<5)+(eta3->position(6)<<6);
      os<< etaPos<<" ";
    }
    else {
      os<< " 0 ";
    }
    const L1MuDTChambPhDigi* seg5 = phiContainer->chPhiSegm2(wheel,1,previousSector,-1);
    if (seg5 && seg5->phi()>111) {
      os<<seg5->phi()-2144<<" "<<seg5->phiB()<<" "<<seg5->code()<<" 1 5 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }

    const L1MuDTChambPhDigi* seg6 = phiContainer->chPhiSegm2(wheel,2,previousSector,-1);
    if (seg6 && seg6->phi()>111) {
      os<<seg6->phi()-2144<<" "<<seg6->phiB()<<" "<<seg6->code()<<" 1 5 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }
    const L1MuDTChambPhDigi* seg7 = phiContainer->chPhiSegm2(wheel,3,previousSector,-1);
    if (seg7 && seg7->phi()>111) {
      os<<seg7->phi()-2144<<" "<<seg7->phiB()<<" "<<seg7->code()<<" 1 5 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }
    const L1MuDTChambPhDigi* seg8 = phiContainer->chPhiSegm2(wheel,4,previousSector,1);
    if (seg8 && seg8->phi()>111) {
      os<<seg8->phi()-2144<<" "<<seg8->phiB()<<" "<<seg8->code()<<" 1 5 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }
    if (eta1) {
      int etaPos=eta1->quality(0)+(eta1->quality(1)<<1) +(eta1->quality(2)<<2) +(eta1->quality(3)<<3) +(eta1->quality(4)<<4)+(eta1->quality(5)<<5)+(eta1->quality(6)<<6);
      os<< etaPos<<" ";
    }
    else {
      os<< " 0 ";
    }
    if (eta2) {
      int etaPos=eta2->quality(0)+(eta2->quality(1)<<1) +(eta2->quality(2)<<2) +(eta2->quality(3)<<3) +(eta2->quality(4)<<4)+(eta2->quality(5)<<5)+(eta2->quality(6)<<6);
      os<< etaPos<<" ";
    }
    else {
      os<< " 0 ";
    }
    if (eta3) {
      int etaPos=eta3->quality(0)+(eta3->quality(1)<<1) +(eta3->quality(2)<<2) +(eta3->quality(3)<<3) +(eta3->quality(4)<<4)+(eta3->quality(5)<<5)+(eta3->quality(6)<<6);
      os<< etaPos<<" ";
    }
    else {
      os<< " 0 ";
    }

    const L1MuDTChambPhDigi* seg9 = phiContainer->chPhiSegm1(wheel,1,sector,0);
    if (seg9 ) {
      os<<seg9->phi()<<" "<<seg9->phiB()<<" "<<seg9->code()<<" 1 0 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }
    const L1MuDTChambPhDigi* seg10 = phiContainer->chPhiSegm1(wheel,2,sector,0);
    if (seg10 ) {
      hasStub=true;
      os<<seg10->phi()<<" "<<seg10->phiB()<<" "<<seg10->code()<<" 1 0 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }
    const L1MuDTChambPhDigi* seg11 = phiContainer->chPhiSegm1(wheel,3,sector,0);
    if (seg11 ) {
      hasStub=true;
      os<<seg11->phi()<<" "<<seg11->phiB()<<" "<<seg11->code()<<" 1 0 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }
    const L1MuDTChambPhDigi* seg12 = phiContainer->chPhiSegm1(wheel,4,sector,0);
    if (seg12 ) {
      hasStub=true;
      os<<seg12->phi()<<" "<<seg12->phiB()<<" "<<seg12->code()<<" 1 0 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }
    const L1MuDTChambThDigi*  eta4   = etaContainer->chThetaSegm(wheel,1,sector,0);
    if (eta4) {
      int etaPos=eta4->position(0)+(eta4->position(1)<<1) +(eta4->position(2)<<2) +(eta4->position(3)<<3) +(eta4->position(4)<<4)+(eta4->position(5)<<5)+(eta4->position(6)<<6);
      os<< etaPos<<" ";
    }
    else {
      os<< " 0 ";
    }

    const L1MuDTChambThDigi*  eta5   = etaContainer->chThetaSegm(wheel,2,sector,0);
    if (eta5) {
      int etaPos=eta5->position(0)+(eta5->position(1)<<1) +(eta5->position(2)<<2) +(eta5->position(3)<<3) +(eta5->position(4)<<4)+(eta5->position(5)<<5)+(eta5->position(6)<<6);
      os<< etaPos<<" ";
    }
    else {
      os<< " 0 ";
    }
    const L1MuDTChambThDigi*  eta6   = etaContainer->chThetaSegm(wheel,3,sector,0);
    if (eta6) {
      int etaPos=eta6->position(0)+(eta6->position(1)<<1) +(eta6->position(2)<<2) +(eta6->position(3)<<3) +(eta6->position(4)<<4)+(eta6->position(5)<<5)+(eta6->position(6)<<6);
      os<< etaPos<<" ";
    }
    else {
      os<< " 0 ";
    }

    const L1MuDTChambPhDigi* seg13 = phiContainer->chPhiSegm2(wheel,1,sector,-1);
    if (seg13 ) {
      os<<seg13->phi()<<" "<<seg13->phiB()<<" "<<seg13->code()<<" 1 1 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }

    const L1MuDTChambPhDigi* seg14 = phiContainer->chPhiSegm2(wheel,2,sector,-1);
    if (seg14 ) {
      hasStub=true;

      os<<seg14->phi()<<" "<<seg14->phiB()<<" "<<seg14->code()<<" 1 1 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }

    const L1MuDTChambPhDigi* seg15 = phiContainer->chPhiSegm2(wheel,3,sector,-1);
    if (seg15 ) {
      hasStub=true;

      os<<seg15->phi()<<" "<<seg15->phiB()<<" "<<seg15->code()<<" 1 1 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }
    const L1MuDTChambPhDigi* seg16 = phiContainer->chPhiSegm2(wheel,4,sector,-1);
    if (seg16 ) {
      hasStub=true;

      os<<seg16->phi()<<" "<<seg16->phiB()<<" "<<seg16->code()<<" 1 1 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }


    if (eta4) {
      int etaPos=eta4->quality(0)+(eta4->quality(1)<<1) +(eta4->quality(2)<<2) +(eta4->quality(3)<<3) +(eta4->quality(4)<<4)+(eta4->quality(5)<<5)+(eta4->quality(6)<<6);
      os<< etaPos<<" ";
    }
    else {
      os<< " 0 ";
    }

    if (eta5) {
      int etaPos=eta5->quality(0)+(eta5->quality(1)<<1) +(eta5->quality(2)<<2) +(eta5->quality(3)<<3) +(eta5->quality(4)<<4)+(eta5->quality(5)<<5)+(eta5->quality(6)<<6);
      os<< etaPos<<" ";
    }
    else {
      os<< " 0 ";
    }
    if (eta6) {
      int etaPos=eta6->quality(0)+(eta6->quality(1)<<1) +(eta6->quality(2)<<2) +(eta6->quality(3)<<3) +(eta6->quality(4)<<4)+(eta6->quality(5)<<5)+(eta6->quality(6)<<6);
      os<< etaPos<<" ";
    }
    else {
      os<< " 0 ";
    }
    const L1MuDTChambPhDigi* seg17 = phiContainer->chPhiSegm1(wheel,1,nextSector,0);
    if (seg17 && seg17->phi()<-112) {
      os<<seg17->phi()+2144<<" "<<seg17->phiB()<<" "<<seg17->code()<<" 1 2 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }
    const L1MuDTChambPhDigi* seg18 = phiContainer->chPhiSegm1(wheel,2,nextSector,0);
    if (seg18 && seg18->phi()<-112) {
      os<<seg18->phi()+2144<<" "<<seg18->phiB()<<" "<<seg18->code()<<" 1 2 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }
    const L1MuDTChambPhDigi* seg19 = phiContainer->chPhiSegm1(wheel,3,nextSector,0);
    if (seg19 && seg19->phi()<-112) {
      os<<seg19->phi()+2144<<" "<<seg19->phiB()<<" "<<seg19->code()<<" 1 2 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }
    const L1MuDTChambPhDigi* seg20 = phiContainer->chPhiSegm1(wheel,4,nextSector,0);
    if (seg20 && seg20->phi()<-112) {
      os<<seg20->phi()+2144<<" "<<seg20->phiB()<<" "<<seg20->code()<<" 1 2 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }
    const L1MuDTChambThDigi*  eta7   = etaContainer->chThetaSegm(wheel,1,nextSector,0);
    if (eta7) {
      int etaPos=eta7->position(0)+(eta7->position(1)<<1) +(eta7->position(2)<<2) +(eta7->position(3)<<3) +(eta7->position(4)<<4)+(eta7->position(5)<<5)+(eta7->position(6)<<6);
      os<< etaPos<<" ";
    }
    else {
      os<< " 0 ";
    }
    const L1MuDTChambThDigi*  eta8   = etaContainer->chThetaSegm(wheel,2,nextSector,0);
    if (eta8) {
      int etaPos=eta8->position(0)+(eta8->position(1)<<1) +(eta8->position(2)<<2) +(eta8->position(3)<<3) +(eta8->position(4)<<4)+(eta8->position(5)<<5)+(eta8->position(6)<<6);
      os<< etaPos<<" ";
    }
    else {
      os<< " 0 ";
    }
    const L1MuDTChambThDigi*  eta9   = etaContainer->chThetaSegm(wheel,3,nextSector,0);
    if (eta9) {
      int etaPos=eta9->position(0)+(eta9->position(1)<<1) +(eta9->position(2)<<2) +(eta9->position(3)<<3) +(eta9->position(4)<<4)+(eta9->position(5)<<5)+(eta9->position(6)<<6);
      os<< etaPos<<" ";
    }
    else {
      os<< " 0 ";
    }
    const L1MuDTChambPhDigi* seg21 = phiContainer->chPhiSegm2(wheel,1,nextSector,-1);
    if (seg21 && seg21->phi()<-112) {
      os<<seg21->phi()+2144<<" "<<seg21->phiB()<<" "<<seg21->code()<<" 1 3 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }

    const L1MuDTChambPhDigi* seg22 = phiContainer->chPhiSegm2(wheel,2,nextSector,-1);
    if (seg22 && seg22->phi()<-112) {
      os<<seg22->phi()+2144<<" "<<seg22->phiB()<<" "<<seg22->code()<<" 1 3 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }
    const L1MuDTChambPhDigi* seg23 = phiContainer->chPhiSegm2(wheel,3,nextSector,-1);
    if (seg23 && seg23->phi()<-112) {
      os<<seg23->phi()+2144<<" "<<seg23->phiB()<<" "<<seg23->code()<<" 1 2 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }
     const L1MuDTChambPhDigi* seg24 = phiContainer->chPhiSegm2(wheel,4,nextSector,-1);
    if (seg24 && seg24->phi()<-112) {
      os<<seg24->phi()+2144<<" "<<seg24->phiB()<<" "<<seg24->code()<<" 1 2 "  ;
    }
    else {
      os<<"-2048 0 0 0 15 "; 
    }
    if (eta7) {
      int etaPos=eta7->quality(0)+(eta7->quality(1)<<1) +(eta7->quality(2)<<2) +(eta7->quality(3)<<3) +(eta7->quality(4)<<4)+(eta7->quality(5)<<5)+(eta7->quality(6)<<6);
      os<< etaPos<<" ";
    }
    else {
      os<< " 0 ";
    }
    if (eta8) {
      int etaPos=eta8->quality(0)+(eta8->quality(1)<<1) +(eta8->quality(2)<<2) +(eta8->quality(3)<<3) +(eta8->quality(4)<<4)+(eta8->quality(5)<<5)+(eta8->quality(6)<<6);
      os<< etaPos<<" ";
    }
    else {
      os<< " 0 ";
    }
    if (eta9) {
      int etaPos=eta9->quality(0)+(eta9->quality(1)<<1) +(eta9->quality(2)<<2) +(eta9->quality(3)<<3) +(eta9->quality(4)<<4)+(eta9->quality(5)<<5)+(eta9->quality(6)<<6);
      os<< etaPos<<" ";
    }
    else {
      os<< " 0 ";
    }

  }
  if (hasStub) {
    std::cout <<os.str()<<std::endl;
  }
}