L1TMuonEndCapTrackProducer Class Reference

#include <L1TMuonEndCapTrackProducer.h>

Inheritance diagram for L1TMuonEndCapTrackProducer:

Public Member Functions
void	beginJob ()
void	endJob ()
	L1TMuonEndCapTrackProducer (const PSet &)
void	produce (edm::Event &, const edm::EventSetup &)
	~L1TMuonEndCapTrackProducer ()
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Public Attributes
const float	ptscale [33]

Private Attributes
edm::EDGetTokenT < CSCCorrelatedLCTDigiCollection >	inputTokenCSC
PrimitiveConverter	primConv_
l1t::EmtfPtAssignment	ptAssignment_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
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 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

Definition at line 49 of file L1TMuonEndCapTrackProducer.h.

Constructor & Destructor Documentation

L1TMuonEndCapTrackProducer::L1TMuonEndCapTrackProducer

(

const PSet &

p

)

Definition at line 40 of file L1TMuonEndCapTrackProducer.cc.

References edm::ParameterSet::getParameter().

                                                                    {

  inputTokenCSC = consumes<CSCCorrelatedLCTDigiCollection>(p.getParameter<edm::InputTag>("CSCInput"));
  
  produces<l1t::RegionalMuonCandBxCollection >("EMTF");
  produces< l1t::EMTFTrackCollection >("EMTF");
  produces< l1t::EMTFHitCollection >("EMTF");  
  produces< l1t::EMTFTrackExtraCollection >("EMTF");
  produces< l1t::EMTFHitExtraCollection >("EMTF");  

}

L1TMuonEndCapTrackProducer::~L1TMuonEndCapTrackProducer ( )

inline

Definition at line 52 of file L1TMuonEndCapTrackProducer.h.

{}

Member Function Documentation

void L1TMuonEndCapTrackProducer::beginJob ( void  )

virtual

Reimplemented from edm::EDProducer.

Definition at line 542 of file L1TMuonEndCapTrackProducer.cc.

{
}

void L1TMuonEndCapTrackProducer::endJob ( void  )

virtual

Reimplemented from edm::EDProducer.

Definition at line 545 of file L1TMuonEndCapTrackProducer.cc.

{
}

void L1TMuonEndCapTrackProducer::produce ( edm::Event &  ev, const edm::EventSetup &  es  )

virtual

Collect Muons from all sectors ////

Make Internal track if //////// //// tracks are found //////////

Implements edm::EDProducer.

Definition at line 53 of file L1TMuonEndCapTrackProducer.cc.

References funct::A, L1TMuon::InternalTrack::addStub(), BestTracks(), l1t::calc_eta_from_theta_rad(), l1t::calc_phi_glob_deg(), l1t::calc_phi_glob_deg_hit(), l1t::calc_phi_glob_rad(), l1t::calc_phi_glob_rad_hit(), l1t::calc_phi_loc_deg(), l1t::calc_phi_loc_rad(), l1t::calc_sector_from_index(), l1t::calc_sector_GMT(), l1t::calc_theta_deg_from_int(), l1t::calc_theta_rad_from_int(), CalcDeltas(), l1t::EMTFHit::Chamber(), RecoTauCleanerPlugins::charge, l1t::EMTFTrackExtra::CreateEMTFTrack(), DeleteDuplicatePatterns(), L1TMuon::InternalTrack::deltas, l1t::EMTFHit::Endcap(), l1t::EMTFTrack::Endcap(), eta, l1t::EMTFTrackExtra::First_BX(), edm::Event::getByToken(), getCharge(), GroupBX(), h, l1t::RegionalMuonCand::hwEta(), l1t::RegionalMuonCand::hwPhi(), l1t::RegionalMuonCand::hwPt(), l1t::RegionalMuonCand::hwQual(), i, l1t::EMTFHitExtra::ImportCSCCorrelatedLCTDigi(), l1t::EMTFHit::ImportCSCDetId(), l1t::EMTFTrack::ImportPtLUT(), j, MakeRegionalCand(), Reference_intrackfit_cff::ME13, alignBH_cfg::mode, l1t::EMTFTrack::Mode(), NUM_SECTORS, GenerateHcalLaserBadRunList::out, Patterns(), L1TMuon::InternalTrack::phi, phi(), l1t::EMTFTrack::Phi_loc_deg(), l1t::EMTFTrack::Phi_loc_int(), l1t::EMTFTrack::Phi_loc_rad(), PhiMatching(), L1TMuon::InternalTrack::phis, Geom::pi(), L1TMuon::InternalTrack::pt, l1t::EMTFTrack::Pt_LUT_addr(), l1t::EMTFTrackExtra::push_HitExtra(), l1t::EMTFTrackExtra::push_HitExtraIndex(), edm::Event::put(), l1t::EMTFTrackExtra::Rank(), L1TMuon::InternalTrack::rank, relativeConstraints::ring, l1t::EMTFHit::Ring(), edm::second(), l1t::EMTFHit::Sector(), l1t::EMTFTrack::Sector(), l1t::EMTFHit::Sector_index(), l1t::EMTFTrack::set_bx(), l1t::EMTFTrack::set_charge(), l1t::EMTFTrack::set_charge_GMT(), l1t::EMTFTrack::set_endcap(), l1t::EMTFTrack::set_eta(), l1t::EMTFTrack::set_eta_GMT(), l1t::EMTFTrackExtra::set_first_bx(), l1t::EMTFTrackExtra::set_isGMT(), l1t::EMTFTrack::set_mode(), l1t::EMTFHit::set_neighbor(), l1t::EMTFTrack::set_phi_glob_deg(), l1t::EMTFTrack::set_phi_glob_rad(), l1t::EMTFTrack::set_phi_GMT(), l1t::EMTFTrack::set_phi_loc_deg(), l1t::EMTFTrack::set_phi_loc_int(), l1t::EMTFTrack::set_phi_loc_rad(), l1t::EMTFTrackExtra::set_phis(), l1t::EMTFTrack::set_pt(), l1t::EMTFTrack::set_pt_GMT(), l1t::EMTFTrack::set_pt_LUT_addr(), l1t::EMTFTrackExtra::set_pt_XML(), l1t::EMTFTrack::set_quality(), l1t::EMTFTrackExtra::set_rank(), l1t::EMTFHit::set_ring(), l1t::EMTFTrackExtra::set_second_bx(), l1t::EMTFTrack::set_sector(), l1t::EMTFTrack::set_sector_GMT(), l1t::EMTFTrack::set_sector_index(), l1t::EMTFTrackExtra::set_straightness(), l1t::EMTFTrackExtra::set_theta_deg(), l1t::EMTFTrackExtra::set_theta_int(), l1t::EMTFTrackExtra::set_theta_rad(), l1t::EMTFTrackExtra::set_thetas(), L1TMuon::InternalTrack::setType(), SortSect(), relativeConstraints::station, l1t::EMTFHit::Station(), l1t::EMTFHit::Strip(), L1TMuon::InternalTrack::theta, theta(), L1TMuon::InternalTrack::thetas, and Zones().

                                            {

  //bool verbose = false;

  //fprintf (write,"12345\n"); //<-- part of printing text file to send verilog code, not needed if George's package is included


  //std::auto_ptr<L1TMuon::InternalTrackCollection> FoundTracks (new L1TMuon::InternalTrackCollection);
  std::auto_ptr<l1t::RegionalMuonCandBxCollection > OutputCands (new l1t::RegionalMuonCandBxCollection);
  std::auto_ptr<l1t::EMTFTrackCollection> OutTracks (new l1t::EMTFTrackCollection);
  std::auto_ptr<l1t::EMTFHitCollection> OutHits (new l1t::EMTFHitCollection);
  std::auto_ptr<l1t::EMTFTrackExtraCollection> OutputTracks (new l1t::EMTFTrackExtraCollection);
  std::auto_ptr<l1t::EMTFHitExtraCollection> OutputHits (new l1t::EMTFHitExtraCollection);

  std::vector<BTrack> PTracks[NUM_SECTORS];

  std::vector<TriggerPrimitive> tester;
  //std::vector<InternalTrack> FoundTracks;
  
  
  
  edm::Handle<CSCCorrelatedLCTDigiCollection> MDC;
  ev.getByToken(inputTokenCSC, MDC);
  std::vector<TriggerPrimitive> out;
  
  auto chamber = MDC->begin();
  auto chend  = MDC->end();
  for( ; chamber != chend; ++chamber ) {
    auto digi = (*chamber).second.first;
    auto dend = (*chamber).second.second;
    for( ; digi != dend; ++digi ) {
      out.push_back(TriggerPrimitive((*chamber).first,*digi));
      l1t::EMTFHitExtra thisHit;
      thisHit.ImportCSCDetId( (*chamber).first );
      thisHit.ImportCSCCorrelatedLCTDigi( *digi );
      if (thisHit.Station() == 1 && thisHit.Ring() == 1 && thisHit.Strip() > 127) thisHit.set_ring(4);
      thisHit.set_neighbor(0);
      OutputHits->push_back( thisHit );
      if ( ((thisHit.Ring() != 1 || thisHit.Station() == 1) && (thisHit.Chamber() % 6 == 2)) ||
       ((thisHit.Ring() == 1 && thisHit.Station()  > 1) && (thisHit.Chamber() % 3 == 1)) ) {
    l1t::EMTFHitExtra neighborHit = thisHit;
    neighborHit.set_neighbor(1);
    OutputHits->push_back( neighborHit );
      }
    }
  }
  


 // auto tpsrc = _tpinputs.cbegin();
  //auto tpend = _tpinputs.cend();
 // for( ; tpsrc != tpend; ++tpsrc ) {
   // edm::Handle<TriggerPrimitiveCollection> tps;
   // ev.getByLabel(*tpsrc,tps);
    auto tp = out.cbegin();
    auto tpend = out.cend();

    for( ; tp != tpend; ++tp ) {
      if(tp->subsystem() == 1)
      {
        //TriggerPrimitiveRef tpref(out,tp - out.cbegin());

        tester.push_back(*tp);

        //std::cout<<"\ntrigger prim found station:"<<tp->detId<CSCDetId>().station()<<std::endl;
      }

     }
   //}
  std::vector<ConvertedHit> CHits[NUM_SECTORS];
  MatchingOutput MO[NUM_SECTORS];

for(int SectIndex=0;SectIndex<NUM_SECTORS;SectIndex++){//perform TF on all 12 sectors




  
    std::vector<ConvertedHit> ConvHits = primConv_.convert(tester,SectIndex);
    CHits[SectIndex] = ConvHits;

    // Fill OutputHits with ConvertedHit information
    for (uint iCHit = 0; iCHit < ConvHits.size(); iCHit++) {
      // bool isMatched = false;

      for (uint iHit = 0; iHit < OutputHits->size(); iHit++) {
        if ( ConvHits.at(iCHit).Station()   == OutputHits->at(iHit).Station() &&
         ( ConvHits.at(iCHit).Id()      == OutputHits->at(iHit).CSC_ID()  ||
           ConvHits.at(iCHit).Id()      == ( (OutputHits->at(iHit).Ring() != 4) // Account for either ME1/1a 
                           ? OutputHits->at(iHit).CSC_ID()    // CSC ID numbering convention
                           : OutputHits->at(iHit).CSC_ID() + 9 ) ) &&
         ConvHits.at(iCHit).Wire()       == OutputHits->at(iHit).Wire()    &&
         ConvHits.at(iCHit).Strip()      == OutputHits->at(iHit).Strip()   &&
         ConvHits.at(iCHit).BX() - 6     == OutputHits->at(iHit).BX()      &&
         ConvHits.at(iCHit).IsNeighbor() == OutputHits->at(iHit).Neighbor() ) {
          // isMatched = true;
          OutputHits->at(iHit).set_neighbor    ( ConvHits.at(iCHit).IsNeighbor());
          OutputHits->at(iHit).set_sector_index( ConvHits.at(iCHit).SectorIndex() );
          OutputHits->at(iHit).set_zone_hit    ( ConvHits.at(iCHit).Zhit()   );
          OutputHits->at(iHit).set_phi_hit     ( ConvHits.at(iCHit).Ph_hit() );
          OutputHits->at(iHit).set_phi_z_val   ( ConvHits.at(iCHit).Phzvl()  );
          OutputHits->at(iHit).set_phi_loc_int ( ConvHits.at(iCHit).Phi()    );
          OutputHits->at(iHit).set_theta_int   ( ConvHits.at(iCHit).Theta()  );

          OutputHits->at(iHit).SetZoneContribution ( ConvHits.at(iCHit).ZoneContribution() );
          OutputHits->at(iHit).set_phi_loc_deg  ( l1t::calc_phi_loc_deg( OutputHits->at(iHit).Phi_loc_int() ) );
          OutputHits->at(iHit).set_phi_loc_rad  ( l1t::calc_phi_loc_rad( OutputHits->at(iHit).Phi_loc_int() ) );
          OutputHits->at(iHit).set_phi_glob_deg ( l1t::calc_phi_glob_deg_hit( OutputHits->at(iHit).Phi_loc_deg(), OutputHits->at(iHit).Sector_index() ) );
          OutputHits->at(iHit).set_phi_glob_rad ( l1t::calc_phi_glob_rad_hit( OutputHits->at(iHit).Phi_loc_rad(), OutputHits->at(iHit).Sector_index() ) );
          OutputHits->at(iHit).set_theta_deg    ( l1t::calc_theta_deg_from_int( OutputHits->at(iHit).Theta_int() ) );
          OutputHits->at(iHit).set_theta_rad    ( l1t::calc_theta_rad_from_int( OutputHits->at(iHit).Theta_int() ) );
          OutputHits->at(iHit).set_eta( l1t::calc_eta_from_theta_rad( OutputHits->at(iHit).Theta_rad() ) * OutputHits->at(iHit).Endcap() );

          OutHits->push_back( OutputHits->at(iHit).CreateEMTFHit() );
        }
      } // End loop: for (uint iHit = 0; iHit < OutputHits->size(); iHit++)

      // if (isMatched == false) {
      //   std::cout << "***********************************************" << std::endl;
      //   std::cout << "Unmatched ConvHit in event " << ev.id().event() << ", SectIndex " << SectIndex << std::endl;
      //   std::cout << "ConvHit: station = " << ConvHits.at(iCHit).Station() << ", CSC ID = " << ConvHits.at(iCHit).Id()
      //          << ", sector index = " << ConvHits.at(iCHit).SectorIndex() << ", subsector = " << ConvHits.at(iCHit).Sub()
      //          << ", wire = " << ConvHits.at(iCHit).Wire() << ", strip = " << ConvHits.at(iCHit).Strip()
      //          << ", BX = " << ConvHits.at(iCHit).BX() << ", neighbor = " << ConvHits.at(iCHit).IsNeighbor() << std::endl;
        
      //   for (uint iHit = 0; iHit < OutputHits->size(); iHit++) {
      //     std::cout << "EMTFHitExtra: station = " << OutputHits->at(iHit).Station() << ", CSC ID = " << OutputHits->at(iHit).CSC_ID()
      //        << ", sector index = " << OutputHits->at(iHit).Sector_index() << ", subsector = " << OutputHits->at(iHit).Subsector() 
      //        << ", wire = " << OutputHits->at(iHit).Wire() << ", strip = " << OutputHits->at(iHit).Strip()
      //        << ", BX = " << OutputHits->at(iHit).BX() << ", neighbor = " << OutputHits->at(iHit).Neighbor()
      //        << ", chamber = " << OutputHits->at(iHit).Chamber() << ", ring = " << OutputHits->at(iHit).Ring() 
      //        << ", endcap = " << OutputHits->at(iHit).Endcap() << ", sector = " << OutputHits->at(iHit).Sector() << std::endl;
      //   }
      // }

    } // End loop: for (uint iCHit = 0; iCHit < ConvHits.size(); iCHit++)
    

    //for(std::vector<ConvertedHit>::iterator h = ConvHits.begin();h != ConvHits.end();h++){

        //if((h->Id()) > 9){h->SetId(h->Id() - 9);h->SetStrip(h->Strip() + 128);}
        //fprintf (write,"0 1   1   %d  %d\n",h->Sub(),h->Station());
        //fprintf (write,"1 %d  %d  %d\n",h->Quality(),h->Pattern(),h->Wire());
        //fprintf (write,"%d    0   %d\n",h->Id(),h->Strip());
    //}




 std::vector<std::vector<ConvertedHit>> GroupedHits = GroupBX(ConvHits);



 std::vector<ZonesOutput> Zout = Zones(GroupedHits);




  std::vector<PatternOutput> Pout = Patterns(Zout);

  PatternOutput Test = DeleteDuplicatePatterns(Pout);

  //PrintQuality(Test.detected);



  SortingOutput Sout = SortSect(Test);



  MatchingOutput Mout = PhiMatching(Sout);
  MO[SectIndex] = Mout;


 std::vector<std::vector<DeltaOutput>> Dout = CalcDeltas(Mout);



  std::vector<BTrack> Bout = BestTracks(Dout);
   PTracks[SectIndex] = Bout;

   
 } // End loop: for(int SectIndex=0;SectIndex<NUM_SECTORS;SectIndex++)



 std::vector<BTrack> PTemp[NUM_SECTORS];
 std::vector<BTrack> AllTracks;
 for (int i=0; i<NUM_SECTORS; i++) PTemp[i] = PTracks[i];


    for(int j=0;j<36;j++){


            if(PTemp[j/3][j%3].phi)//no track
                AllTracks.push_back(PTemp[j/3][j%3]);

        

    }

  
  std::vector<l1t::RegionalMuonCand> tester1;
  std::vector<std::pair<int,l1t::RegionalMuonCand>> holder;

  for(unsigned int fbest=0;fbest<AllTracks.size();fbest++){

    if(AllTracks[fbest].phi){

        InternalTrack tempTrack;
        tempTrack.setType(2);
        tempTrack.phi = AllTracks[fbest].phi;
        tempTrack.theta = AllTracks[fbest].theta;
        tempTrack.rank = AllTracks[fbest].winner.Rank();
        tempTrack.deltas = AllTracks[fbest].deltas;
        std::vector<int> ps, ts;

        l1t::EMTFTrackExtra thisTrack;
        thisTrack.set_phi_loc_int ( AllTracks[fbest].phi           );
        thisTrack.set_theta_int   ( AllTracks[fbest].theta         );
        thisTrack.set_rank        ( AllTracks[fbest].winner.Rank() );
        // thisTrack.set_deltas        ( AllTracks[fbest].deltas        );
        int tempStraightness = 0;
        int tempRank = thisTrack.Rank();
        if (tempRank & 64)
          tempStraightness |= 4;
        if (tempRank & 16)
          tempStraightness |= 2;
        if (tempRank & 4)
          tempStraightness |= 1;
        thisTrack.set_straightness ( tempStraightness );


        int sector = -1;
        bool ME13 = false;
        int me1address = 0, me2address = 0, CombAddress = 0, mode_uncorr = 0;
        int ebx = 20, sebx = 20;
        int phis[4] = {-99,-99,-99,-99};

        for(std::vector<ConvertedHit>::iterator A = AllTracks[fbest].AHits.begin();A != AllTracks[fbest].AHits.end();A++){

            if(A->Phi() != -999){
              
                    l1t::EMTFHitExtra thisHit;
                    // thisHit.ImportCSCDetId( A->TP().detId<CSCDetId>() );

                for (uint iHit = 0; iHit < OutputHits->size(); iHit++) {
                  if ( A->TP().detId<CSCDetId>().station() == OutputHits->at(iHit).Station() and
                       A->TP().getCSCData().cscID          == OutputHits->at(iHit).CSC_ID()  and
                       A->Wire()                           == OutputHits->at(iHit).Wire()    and
                       A->Strip()                          == OutputHits->at(iHit).Strip()   and
                       A->TP().getCSCData().bx - 6         == OutputHits->at(iHit).BX() ) {
                    thisHit = OutputHits->at(iHit);
                    thisTrack.push_HitExtraIndex(iHit);
                    thisTrack.push_HitExtra(thisHit); // Done before theta windows are applied ... how can we do it after? - AWB 29.04.16
                    break;
                  }
                }
                thisTrack.set_endcap       ( thisHit.Endcap() );
                thisTrack.set_sector_index ( thisHit.Sector_index() );
                thisTrack.set_sector       ( l1t::calc_sector_from_index( thisHit.Sector_index() ) );
                thisTrack.set_sector_GMT   ( l1t::calc_sector_GMT( thisHit.Sector() ) );

                int station = A->TP().detId<CSCDetId>().station();
                int id = A->TP().getCSCData().cscID;
                int trknm = A->TP().getCSCData().trknmb;
                
                phis[station-1] = A->Phi();
                
                
                if(A->TP().getCSCData().bx < ebx){
                    sebx = ebx;
                    ebx = A->TP().getCSCData().bx;
                }
                else if(A->TP().getCSCData().bx < sebx){
                    sebx = A->TP().getCSCData().bx;
                }

                tempTrack.addStub(A->TP());
                ps.push_back(A->Phi());
                ts.push_back(A->Theta());
                
                sector = A->SectorIndex();//(A->TP().detId<CSCDetId>().endcap() -1)*6 + A->TP().detId<CSCDetId>().triggerSector() - 1;
                //std::cout<<"Q: "<<A->Quality()<<", keywire: "<<A->Wire()<<", strip: "<<A->Strip()<<std::endl;

                switch(station){
                    case 1: mode_uncorr |= 8;break;
                    case 2: mode_uncorr |= 4;break;
                    case 3: mode_uncorr |= 2;break;
                    case 4: mode_uncorr |= 1;break;
                    default: mode_uncorr |= 0;
                }


                if(A->TP().detId<CSCDetId>().station() == 1 && A->TP().detId<CSCDetId>().ring() == 3)
                    ME13 = true;

                if(station == 1 && id > 3 && id < 7){

                    int sub = 2;
                    if(A->TP().detId<CSCDetId>().chamber()%6 > 2)
                        sub = 1;

                    me1address = id;
                    me1address -= 3;
                    me1address += 3*(sub - 1);
                    me1address = id<<1;//
                    me1address |= trknm-1;

                }

                if(station == 2 && id > 3){

                    me2address = id;
                    me2address -= 3;
                    me2address = me2address<<1;
                    me2address |= trknm-1;

                }


            }

        }
        
        int mode = 0;
        if(tempTrack.rank & 32)
            mode |= 8;
        if(tempTrack.rank & 8)
            mode |= 4;
        if(tempTrack.rank & 2)
            mode |= 2;
        if(tempTrack.rank & 1)
            mode |= 1;

        tempTrack.phis = ps;
        tempTrack.thetas = ts;

        // // Before Mulhearn cleanup, May 11
        // unsigned long xmlpt_address = 0;
        // float xmlpt = CalculatePt(tempTrack, es, mode, &xmlpt_address);
        
        // After Mulhearn cleanup, May 11
        unsigned long xmlpt_address = ptAssignment_.calculateAddress(tempTrack, es, mode);
        float xmlpt = ptAssignment_.calculatePt(xmlpt_address);

        tempTrack.pt = xmlpt*1.4;
        //FoundTracks->push_back(tempTrack);

        CombAddress = (me2address<<4) | me1address;

        int charge = getCharge(phis[0],phis[1],phis[2],phis[3],mode);

        l1t::RegionalMuonCand outCand = MakeRegionalCand(xmlpt*1.4,AllTracks[fbest].phi,AllTracks[fbest].theta,
                                 charge,mode,CombAddress,sector);

        float theta_angle = l1t::calc_theta_rad_from_int( AllTracks[fbest].theta ); 
        float eta = l1t::calc_eta_from_theta_rad( theta_angle );

        thisTrack.set_phi_loc_deg  ( l1t::calc_phi_loc_deg( thisTrack.Phi_loc_int() ) );
        thisTrack.set_phi_loc_rad  ( l1t::calc_phi_loc_rad( thisTrack.Phi_loc_int() ) );
        thisTrack.set_phi_glob_deg ( l1t::calc_phi_glob_deg( thisTrack.Phi_loc_deg(), thisTrack.Sector() ) );
        thisTrack.set_phi_glob_rad ( l1t::calc_phi_glob_rad( thisTrack.Phi_loc_rad(), thisTrack.Sector() ) );
        thisTrack.set_quality    ( outCand.hwQual());
        thisTrack.set_mode       ( mode            ); 
        thisTrack.set_first_bx   ( ebx - 6         ); 
        thisTrack.set_second_bx  ( sebx - 6        ); 
        thisTrack.set_bx         ( thisTrack.First_BX() );
        thisTrack.set_phis       ( ps              );
        thisTrack.set_thetas     ( ts              );
        thisTrack.set_pt         ( xmlpt*1.4       );
        thisTrack.set_pt_XML     ( xmlpt           );
        thisTrack.set_pt_LUT_addr( xmlpt_address   );
        thisTrack.set_charge     ( (charge == 1) ? -1 : 1 ); // uGMT uses opposite of physical charge (to match pdgID)
        thisTrack.set_charge_GMT ( charge          );
        thisTrack.set_theta_rad  ( theta_angle     );
        thisTrack.set_theta_deg  ( theta_angle * 180/Geom::pi() );
        thisTrack.set_eta        ( eta  * thisTrack.Endcap() );
        thisTrack.set_pt_GMT     ( outCand.hwPt()  ); 
        thisTrack.set_phi_GMT    ( outCand.hwPhi() );
        thisTrack.set_eta_GMT    ( outCand.hwEta() );

        thisTrack.ImportPtLUT    ( thisTrack.Mode(), thisTrack.Pt_LUT_addr() );

        // thisTrack.phi_loc_rad(); // Need to implement - AWB 04.04.16
        // thisTrack.phi_glob_rad(); // Need to implement - AWB 04.04.16

        // // Optimal emulator configuration - AWB 29.03.16
        // std::pair<int,l1t::RegionalMuonCand> outPair(sebx,outCand);

        // Actual setting in firmware - AWB 12.04.16
        std::pair<int,l1t::RegionalMuonCand> outPair(ebx,outCand);

        if(!ME13 && fabs(eta) > 1.1) {
          // // Extra debugging output - AWB 29.03.16
          // std::cout << "Input: eBX = " << ebx << ", seBX = " << sebx << ", pt = " << xmlpt*1.4 
          //        << ", phi = " << AllTracks[fbest].phi << ", eta = " << eta 
          //        << ", theta = " << AllTracks[fbest].theta << ", sign = " << 1 
          //        << ", quality = " << mode << ", trackaddress = " << 1 
          //        << ", sector = " << sector << std::endl;
          // std::cout << "Output: BX = " << ebx << ", hwPt = " << outCand.hwPt() << ", hwPhi = " << outCand.hwPhi() 
          //        << ", hwEta = " << outCand.hwEta() << ", hwSign = " << outCand.hwSign() 
          //        << ", hwQual = " << outCand.hwQual() << ", link = " << outCand.link()
          //        << ", processor = " << outCand.processor() 
          //        << ", trackFinderType = " << outCand.trackFinderType() << std::endl;
            holder.push_back(outPair);
            thisTrack.set_isGMT( 1 );
        }
        OutputTracks->push_back( thisTrack );
        OutTracks->push_back( thisTrack.CreateEMTFTrack() );
    }
  }
  
OutputCands->setBXRange(-2,2);

for(int sect=0;sect<12;sect++){

    for(unsigned int h=0;h<holder.size();h++){
    
        int bx = holder[h].first - 6;
        int sector = holder[h].second.processor();
        if(holder[h].second.trackFinderType() == 3)
            sector += 6;
    
        if(sector == sect){
            OutputCands->push_back(bx,holder[h].second);
        }
        
    }
}


//ev.put( FoundTracks, "DataITC");
ev.put( OutputCands, "EMTF");
 ev.put( OutputHits, "EMTF"); 
 ev.put( OutputTracks, "EMTF");
  //std::cout<<"End Upgraded Track Finder Prducer:::::::::::::::::::::::::::\n:::::::::::::::::::::::::::::::::::::::::::::::::\n\n";

}//analyzer

Member Data Documentation

edm::EDGetTokenT<CSCCorrelatedLCTDigiCollection> L1TMuonEndCapTrackProducer::inputTokenCSC

private

Definition at line 87 of file L1TMuonEndCapTrackProducer.h.

PrimitiveConverter L1TMuonEndCapTrackProducer::primConv_

private

Definition at line 85 of file L1TMuonEndCapTrackProducer.h.

l1t::EmtfPtAssignment L1TMuonEndCapTrackProducer::ptAssignment_

private

Definition at line 88 of file L1TMuonEndCapTrackProducer.h.

const float L1TMuonEndCapTrackProducer::ptscale[33]

Initial value:

= { 
    -1.,   0.0,   1.5,   2.0,   2.5,   3.0,   3.5,   4.0,
    4.5,   5.0,   6.0,   7.0,   8.0,  10.0,  12.0,  14.0,  
    16.0,  18.0,  20.0,  25.0,  30.0,  35.0,  40.0,  45.0, 
    50.0,  60.0,  70.0,  80.0,  90.0, 100.0, 120.0, 140.0, 1.E6 }

Definition at line 77 of file L1TMuonEndCapTrackProducer.h.