L1AnalysisDTTF.cc

Go to the documentation of this file.

#include <cmath>
#include "L1Trigger/L1TNtuples/interface/L1AnalysisDTTF.h"
#include <cmath>


L1Analysis::L1AnalysisDTTF::L1AnalysisDTTF()
{        
}


L1Analysis::L1AnalysisDTTF::~L1AnalysisDTTF()
{

}

void L1Analysis::L1AnalysisDTTF::SetDTPH(const edm::Handle<L1MuDTChambPhContainer > L1MuDTChambPhContainer, unsigned int maxDTPH)
{
   
  const L1MuDTChambPhContainer::Phi_Container *PhContainer = L1MuDTChambPhContainer->getContainer();
    
   dttf_.phSize = PhContainer->size();
   int iphtr=0;
   double pig=acos(-1);
   for( L1MuDTChambPhContainer::Phi_Container::const_iterator 
       DTPhDigiItr =  PhContainer->begin() ;
       DTPhDigiItr != PhContainer->end() ;
       ++DTPhDigiItr ) 
    {        
      if((unsigned int) iphtr>maxDTPH-1) continue;
      dttf_.phBx.push_back     (  DTPhDigiItr->bxNum() - DTPhDigiItr->Ts2Tag()+1 );
      dttf_.phWh.push_back     (  DTPhDigiItr->whNum() );
      dttf_.phSe.push_back     (  DTPhDigiItr->scNum() );
      dttf_.phSt.push_back     (  DTPhDigiItr->stNum() );
      dttf_.phAng.push_back    (  DTPhDigiItr->phi()   );
      dttf_.phBandAng.push_back(  DTPhDigiItr->phiB()  );
      double locphi = (((double)DTPhDigiItr->phi())/4096.0);
      double newphi = locphi+(DTPhDigiItr->scNum()*(pig/6.));
      if (newphi>2*pig) newphi-=2*pig;
      if (newphi<0    ) newphi+=2*pig;
      dttf_.phGlobPhi.push_back(  newphi               );
      dttf_.phCode.push_back   (  DTPhDigiItr->code()  );
      dttf_.phX.push_back      (  DTPhDigiItr->scNum() );
      dttf_.phY.push_back      (  DTPhDigiItr->stNum()+4*(DTPhDigiItr->whNum()+2) );
      
      iphtr++;
    } 
  
}


void L1Analysis::L1AnalysisDTTF::SetDTTH(const edm::Handle<L1MuDTChambThContainer > L1MuDTChambThContainer, unsigned int maxDTTH)
{
   
   const L1MuDTChambThContainer::The_Container* ThContainer = L1MuDTChambThContainer->getContainer();

   int ithtr=0;
   dttf_.thSize = ThContainer->size();
   dttf_.thTheta.ResizeTo(dttf_.thSize,7);
   dttf_.thCode.ResizeTo(dttf_.thSize,7);  


   for( L1MuDTChambThContainer::The_Container::const_iterator 
     DTThDigiItr =  ThContainer->begin() ;
       DTThDigiItr != ThContainer->end() ;
       ++DTThDigiItr ) 
     {  
     
      if((unsigned int) ithtr>maxDTTH-1) continue;
      dttf_.thBx.push_back( DTThDigiItr->bxNum() + 1 );
      dttf_.thWh.push_back( DTThDigiItr->whNum() );
      dttf_.thSe.push_back( DTThDigiItr->scNum() );
      dttf_.thSt.push_back( DTThDigiItr->stNum() );
      dttf_.thX.push_back( DTThDigiItr->stNum()+4*(DTThDigiItr->whNum()+2) );
//    int xpos = iwh*4+ist+1; ????
      dttf_.thY.push_back( DTThDigiItr->scNum() );

      for (unsigned int j = 0; j < 7; j++)
    {
         dttf_.thTheta[ithtr][j] = DTThDigiItr->position(j);
         dttf_.thCode[ithtr][j]  = DTThDigiItr->code(j);
    }
      
      ithtr++;
     
    }
}


void L1Analysis::L1AnalysisDTTF::SetDTTR(const edm::Handle<L1MuDTTrackContainer > L1MuDTTrackContainer, unsigned int maxDTTR)
{
  double pig=acos(-1);
  const L1MuDTTrackContainer::TrackContainer *tr =  L1MuDTTrackContainer->getContainer();
  int idttr = 0 ;
  dttf_.trSize = tr->size();
  for ( L1MuDTTrackContainer::TrackContainer::const_iterator i 
      = tr->begin(); i != tr->end(); ++i ) {
        if((unsigned int)idttr>maxDTTR-1) continue; 
    dttf_.trBx.push_back(  i->bx()+1 );  
    dttf_.trTag.push_back(  i->TrkTag() );  
    dttf_.trQual.push_back(  i->quality_packed() ); 
    dttf_.trPtPck.push_back(  i->pt_packed() );
    dttf_.trPtVal.push_back(  i->ptValue() );
    dttf_.trPhiPck.push_back(  i->phi_packed() ); 
    dttf_.trPhiVal.push_back(  i->phiValue() ); 
    dttf_.trEtaPck.push_back(  i->eta_packed() );
    dttf_.trEtaVal.push_back(  i->etaValue() );

        int phi_local = i->phi_packed();//range: 0 < phi_local < 31 
        if(phi_local > 15) phi_local -= 32; //range: -16 < phi_local < 15
        double phi_global = (phi_local*(pig/72.))+((pig/6.)*i->scNum());// + 12*i->scNum(); //range: -16 < phi_global < 147
        if(phi_global < 0) phi_global+=2*pig; //range: 0 < phi_global < 147
        if(phi_global > 2*pig) phi_global-=2*pig; //range: 0 < phi_global < 143
    dttf_.trPhiGlob.push_back(  phi_global );
    dttf_.trChPck.push_back(  i->charge_packed() ); 
    dttf_.trWh.push_back(  i->whNum() );
    dttf_.trSc.push_back(  i->scNum() );
    unsigned int addresses=0;
    for (int KK = 0; KK<4; KK++){
      addresses|=((unsigned int)(i->stNum(KK+1))<<(KK*4));
    }
    dttf_.trAddress.push_back(addresses);
    idttr++;
  }
 
}