MBLTCollection.h

Go to the documentation of this file.

#ifndef MBLTCollection_h_
#define MBLTCollection_h_
//
// Class: L1ITMu::
//
// Info: This track represents a DT(1 station) plus eventual RPC station(s)
//       based track seed, from which a full multi-station track can be
//       built.
//
// Author:
// Carlo Battilana - Giuseppe Codispoti
//

#include <iostream>

#include "L1Trigger/L1TMuon/interface/deprecate/MuonTriggerPrimitive.h"

#include "DataFormats/Common/interface/RefToBase.h"
#include "DataFormats/MuonDetId/interface/DTChamberId.h"
#include "DataFormats/Common/interface/Ref.h"
#include "DataFormats/Common/interface/Ptr.h"


using namespace L1TMuon;

namespace L1TwinMux {

   typedef edm::Ref<TriggerPrimitiveCollection> TriggerPrimitiveRef;
   typedef std::vector<TriggerPrimitiveRef> TriggerPrimitiveList;

  class MBLTCollection {

  public:

    enum subsystem_offset{ kDT, kRPCb, kCSC, kRPCf };

    struct primitiveAssociation {
      std::vector< size_t > rpcIn;
      std::vector< size_t > rpcOut;
    };

    enum bxMatch { NOMATCH, INMATCH, OUTMATCH, FULLMATCH };

    MBLTCollection() :_wheel(0),_sector(0),_station(0) {};

    MBLTCollection( const DTChamberId & dtId );
    ~MBLTCollection() {};

    void addStub( const TriggerPrimitiveRef & stub );


    const  TriggerPrimitiveList & getDtSegments() const {
      return _dtAssociatedStubs;
    }

    const  TriggerPrimitiveList & getRpcInner() const {
      return _rpcInAssociatedStubs;
    }

    const  TriggerPrimitiveList & getRpcOuter() const {
      return _rpcOutAssociatedStubs;
    }

    inline int wheel() const { return _wheel; }

    inline int sector() const { return _sector; }

    inline int station() const { return _station; }

    inline DTChamberId detId() const { return DTChamberId( _wheel, _station, _sector ); }

     TriggerPrimitiveList getRpcInAssociatedStubs( size_t dtIndex ) const;

     TriggerPrimitiveList getRpcOutAssociatedStubs( size_t dtIndex ) const;

    // build association map among dt and rpc primitives
    void associate( double );

    // look for common rpc hists among 2 dt primitives
    bxMatch haveCommonRpc( size_t dt1, size_t dt2 ) const;

      TriggerPrimitiveList getRpcInUnassociatedStubs() const;

      TriggerPrimitiveList getRpcOutUnassociatedStubs() const;

     std::vector< std::pair<  TriggerPrimitiveList,  TriggerPrimitiveList > >
       getUnassociatedRpcClusters( double minRpcPhi ) const;

   private :
     void
       getUnassociatedRpcClusters( const std::vector< size_t > & rpcUnass,
                                 const  TriggerPrimitiveList & rpcList,
                                 double minRpcPhi,
                                 std::vector< std::vector <size_t> > & clusters ) const;
     //
     bool areCloseClusters( std::vector< size_t > & cluster1,
                          std::vector< size_t > & cluster2,
                          const  TriggerPrimitiveList & rpcList1,
                          const  TriggerPrimitiveList & rpcList2,
                          double minRpcPhi ) const;
     //
     size_t reduceRpcClusters( std::vector< std::vector <size_t> > & tmpClusters,
                             const  TriggerPrimitiveList & rpcList,
                             double minRpcPhi ) const;
  private :

     TriggerPrimitiveList _dtAssociatedStubs;

     TriggerPrimitiveList _rpcInAssociatedStubs;

     TriggerPrimitiveList _rpcOutAssociatedStubs;

    int _wheel, _sector, _station;

    // association map among dt and rpc primitives
    std::vector< primitiveAssociation > _dtMapAss;
    // association map among dt and rpc primitives
    primitiveAssociation _rpcMapUnass;

  };
}

#endif