Migration of L1 tracker primitives from SLHCDEV

DataFormats/L1TrackTrigger/BuildFile.xml

@@ -0,0 +1,9 @@
+<use   name="DataFormats/Common"/>
+<use   name="DataFormats/DetId"/>
+<use   name="DataFormats/GeometryCommonDetAlgo"/>
+<use   name="DataFormats/GeometryVector"/>
+<use   name="DataFormats/Phase2TrackerDigi"/>
+<export>
+  <lib   name="1"/>
+</export>
+

DataFormats/L1TrackTrigger/interface/TTCluster.h

@@ -0,0 +1,185 @@
+/*! \class   TTCluster
+ *  \brief   Class to store the L1 Track Trigger clusters
+ *  \details After moving from SimDataFormats to DataFormats,
+ *           the template structure of the class was maintained
+ *           in order to accomodate any types other than Phase2TrackerDigis
+ *           in case there is such a need in the future.
+ *
+ *  \author Nicola Pozzobon
+ *  \author Emmanuele Salvati
+ *  \date   2013, Jul 12
+ *
+ */
+
+#ifndef L1_TRACK_TRIGGER_CLUSTER_FORMAT_H
+#define L1_TRACK_TRIGGER_CLUSTER_FORMAT_H
+
+#include "DataFormats/Common/interface/Ref.h"
+#include "DataFormats/Common/interface/Ptr.h"
+#include "DataFormats/Common/interface/DetSet.h"
+#include "DataFormats/Common/interface/DetSetVector.h"
+#include "DataFormats/DetId/interface/DetId.h"
+#include "DataFormats/Phase2TrackerDigi/interface/Phase2TrackerDigi.h"
+#include "DataFormats/GeometryCommonDetAlgo/interface/MeasurementPoint.h"
+#include "DataFormats/GeometryVector/interface/GlobalPoint.h" /// NOTE: this is needed even if it seems not
+
+template< typename T >
+class TTCluster
+{
+  public:
+    /// Constructors
+    TTCluster();
+    TTCluster( std::vector< T > aHits,
+               DetId aDetId,
+               unsigned int aStackMember,
+               bool storeLocal );
+
+    /// Destructor
+    ~TTCluster();
+
+    /// Data members:   getABC( ... )
+    /// Helper methods: findABC( ... )
+
+    /// Hits in the Cluster
+    std::vector< T > getHits() const                   { return theHits; }
+    void             setHits( std::vector< T > aHits ) { theHits = aHits; }
+
+    /// Detector element
+    DetId        getDetId() const         { return theDetId; }
+    void         setDetId( DetId aDetId ) { theDetId = aDetId; }
+    unsigned int getStackMember() const                      { return theStackMember; }
+    void         setStackMember( unsigned int aStackMember ) { theStackMember = aStackMember; }
+
+    /// Rows and columns to get rid of Digi collection
+    std::vector< int > findRows() const;
+    std::vector< int > findCols() const;
+    void setCoordinates( std::vector< int > a, std::vector< int > b ) { theRows = a; theCols = b; }
+    std::vector< int > getRows() const { return theRows; }
+    std::vector< int > getCols() const { return theCols; }
+
+    /// Cluster width
+    unsigned int findWidth() const;
+
+    /// Single hit coordinates
+    /// Average cluster coordinates
+    MeasurementPoint findHitLocalCoordinates( unsigned int hitIdx ) const;
+    MeasurementPoint findAverageLocalCoordinates() const;
+
+    /// Information
+    std::string print( unsigned int i = 0 ) const;
+
+  private:
+    /// Data members
+    std::vector< T >  theHits;
+    DetId             theDetId;
+    unsigned int      theStackMember;
+
+    std::vector< int > theRows;
+    std::vector< int > theCols;
+
+}; /// Close class
+
+/*! \brief   Implementation of methods
+ *  \details Here, in the header file, the methods which do not depend
+ *           on the specific type <T> that can fit the template.
+ *           Other methods, with type-specific features, are implemented
+ *           in the source file.
+ */
+
+/// Default Constructor
+/// NOTE: to be used with setSomething(...) methods
+template< typename T >
+TTCluster< T >::TTCluster()
+{
+  /// Set default data members
+  theHits.clear();
+  theDetId = 0;
+  theStackMember = 0;
+
+  theRows.clear();
+  theCols.clear();
+}
+
+/// Another Constructor
+template< typename T >
+TTCluster< T >::TTCluster( std::vector< T > aHits,
+                           DetId aDetId,
+                           unsigned int aStackMember,
+                           bool storeLocal )
+{
+  /// Set data members
+  this->setHits( aHits );
+  this->setDetId( aDetId );
+  this->setStackMember( aStackMember );    
+
+  theRows.clear();
+  theCols.clear();
+  if ( storeLocal )
+  {
+    this->setCoordinates( this->findRows(), this->findCols() );
+  }
+}
+
+/// Destructor
+template< typename T >
+TTCluster< T >::~TTCluster(){}
+
+/// Cluster width
+template< >
+unsigned int TTCluster< edm::Ref< edm::DetSetVector< Phase2TrackerDigi >, Phase2TrackerDigi > >::findWidth() const;
+
+/// Single hit coordinates
+/// Average cluster coordinates
+template< >
+MeasurementPoint TTCluster< edm::Ref< edm::DetSetVector< Phase2TrackerDigi >, Phase2TrackerDigi > >::findHitLocalCoordinates( unsigned int hitIdx ) const;
+
+template< >
+MeasurementPoint TTCluster< edm::Ref< edm::DetSetVector< Phase2TrackerDigi >, Phase2TrackerDigi > >::findAverageLocalCoordinates() const;
+
+/// Operations with coordinates stored locally
+template< typename T > 
+std::vector< int > TTCluster< T >::findRows() const
+{
+  std::vector< int > temp;
+  return temp;
+}
+
+template< typename T > 
+std::vector< int > TTCluster< T >::findCols() const
+{
+  std::vector< int > temp;
+  return temp;
+}
+
+template< >
+std::vector< int > TTCluster< edm::Ref< edm::DetSetVector< Phase2TrackerDigi >, Phase2TrackerDigi > >::findRows() const;
+
+template< >
+std::vector< int > TTCluster< edm::Ref< edm::DetSetVector< Phase2TrackerDigi >, Phase2TrackerDigi > >::findCols() const;
+
+/// Information
+template< typename T >
+std::string TTCluster< T >::print( unsigned int i ) const
+{
+  std::string padding("");
+  for ( unsigned int j = 0; j != i; ++j )
+  {
+    padding+="\t";
+  }
+
+  std::stringstream output;
+  output<<padding<<"TTCluster:\n";
+  padding+='\t';
+  output << padding << "DetId: " << theDetId.rawId() << '\n';
+  output << padding << "member: " << theStackMember << ", cluster size: " << theHits.size() << '\n';
+  return output.str();
+}
+
+template< typename T >
+std::ostream& operator << ( std::ostream& os, const TTCluster< T >& aTTCluster )
+{
+  return ( os << aTTCluster.print() );
+}
+
+#endif
+

DataFormats/L1TrackTrigger/interface/TTStub.h

@@ -0,0 +1,171 @@
+/*! \class   TTStub
+ *  \brief   Class to store the L1 Track Trigger stubs
+ *  \details After moving from SimDataFormats to DataFormats,
+ *           the template structure of the class was maintained
+ *           in order to accomodate any types other than Phase2TrackerDigis
+ *           in case there is such a need in the future.
+ *
+ *  \author Andrew W. Rose
+ *  \author Nicola Pozzobon
+ *  \date   2013, Jul 12
+ *
+ */
+
+#ifndef L1_TRACK_TRIGGER_STUB_FORMAT_H
+#define L1_TRACK_TRIGGER_STUB_FORMAT_H
+
+#include "DataFormats/GeometryVector/interface/GlobalVector.h"
+#include "DataFormats/L1TrackTrigger/interface/TTCluster.h"
+#include "DataFormats/Common/interface/DetSetVectorNew.h"
+
+template< typename T >
+class TTStub
+{
+  public:
+    /// Constructors
+    TTStub();
+    TTStub( DetId aDetId );
+
+    /// Destructor
+    ~TTStub();
+
+    /// Data members:   getABC( ... )
+    /// Helper methods: findABC( ... )
+
+    /// Clusters composing the Stub
+    const edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > >&         getClusterRef( unsigned int hitIdentifier ) const;
+    void addClusterRef( edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > > aTTCluster );
+
+    /// Detector element
+    DetId getDetId() const         { return theDetId; }
+    void  setDetId( DetId aDetId ) { theDetId = aDetId; }
+
+    /// Trigger information
+    double getTriggerDisplacement() const;              /// In FULL-STRIP units! (hence, not implemented herein)
+    void   setTriggerDisplacement( int aDisplacement ); /// In HALF-STRIP units!
+    double getTriggerOffset() const;         /// In FULL-STRIP units! (hence, not implemented herein)
+    void   setTriggerOffset( int anOffset ); /// In HALF-STRIP units!
+
+    /// CBC3-style trigger information
+    /// for sake of simplicity, these methods are
+    /// slightly out of the getABC(...)/findABC(...) rule
+    double getTriggerPosition() const; /// In FULL-STRIP units!
+    double getTriggerBend() const;     /// In FULL-STRIP units!
+
+    /// Information
+    std::string print( unsigned int i = 0 ) const;
+
+  private:
+    /// Data members
+    DetId theDetId;
+    edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > > theClusterRef0;
+    edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > > theClusterRef1;
+    int theDisplacement;
+    int theOffset;
+
+}; /// Close class
+
+/*! \brief   Implementation of methods
+ *  \details Here, in the header file, the methods which do not depend
+ *           on the specific type <T> that can fit the template.
+ *           Other methods, with type-specific features, are implemented
+ *           in the source file.
+ */
+
+/// Default Constructor
+template< typename T >
+TTStub< T >::TTStub()
+{
+  /// Set default data members
+  theDetId = 0;
+  theDisplacement = 999999;
+  theOffset = 0;
+}
+
+/// Another Constructor
+template< typename T >
+TTStub< T >::TTStub( DetId aDetId )
+{
+  /// Set data members
+  this->setDetId( aDetId );
+
+  /// Set default data members
+  theDisplacement = 999999;
+  theOffset = 0;
+}
+
+/// Destructor
+template< typename T >
+TTStub< T >::~TTStub(){}
+
+/// Get the Reference to a Cluster
+template< typename T >
+const edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > >& TTStub< T >::getClusterRef( unsigned int hitIdentifier ) const
+{
+  return hitIdentifier==0 ? theClusterRef0 : theClusterRef1;
+}
+
+template< typename T >
+void TTStub< T >::addClusterRef( edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > > aTTCluster )
+{
+  if(aTTCluster->getStackMember() == 0) theClusterRef0 = aTTCluster;
+  else if (aTTCluster->getStackMember() == 1) theClusterRef1 = aTTCluster;
+}
+
+/// Trigger info
+template< typename T >
+double TTStub< T >::getTriggerDisplacement() const { return 0.5*theDisplacement; }
+
+template< typename T >
+void TTStub< T >::setTriggerDisplacement( int aDisplacement ) { theDisplacement = aDisplacement; }
+
+template< typename T >
+double TTStub< T >::getTriggerOffset() const { return 0.5*theOffset; }
+
+template< typename T >
+void TTStub< T >::setTriggerOffset( int anOffset ) { theOffset = anOffset; }
+
+/// CBC3-style trigger info
+template< typename T >
+double TTStub< T >::getTriggerPosition() const
+{
+  return this->getClusterRef(0)->findAverageLocalCoordinates().x();
+}
+
+template< typename T >
+double TTStub< T >::getTriggerBend() const
+{
+  if ( theDisplacement == 999999 )
+    return theDisplacement;
+
+  return 0.5*( theDisplacement - theOffset );
+}
+
+/// Information
+template< typename T >
+std::string TTStub< T >::print( unsigned int i ) const
+{
+  std::string padding("");
+  for ( unsigned int j = 0; j != i; ++j )
+  {
+    padding+="\t";
+  }
+
+  std::stringstream output;
+  output<<padding<<"TTStub:\n";
+  padding+='\t';
+  output << padding << "DetId: " << theDetId.rawId() << ", position: " << this->getTriggerPosition();
+  output << ", bend: " << this->getTriggerBend() << '\n';
+  unsigned int iClu = 0;
+  output << padding << "cluster 0: address: " << theClusterRef0.get();
+  output << ", cluster size: " << theClusterRef0->getHits().size() << '\n';
+  output << padding << "cluster 1: address: " << theClusterRef1.get();
+  output << ", cluster size: " << theClusterRef1->getHits().size() << '\n';
+  return output.str();
+}
+
+template< typename T >
+std::ostream& operator << ( std::ostream& os, const TTStub< T >& aTTStub ) { return ( os << aTTStub.print() ); }
+
+#endif
+

DataFormats/L1TrackTrigger/interface/TTTypes.h

@@ -0,0 +1,24 @@
+/*! \brief   Definition of all the relevant data types
+ *
+ *  \author Andrew W. Rose
+ *  \author Nicola Pozzobon
+ *  \date   2013, Jul 12
+ *
+ */
+
+#ifndef L1_TRACK_TRIGGER_TYPES_H
+#define L1_TRACK_TRIGGER_TYPES_H
+
+/// Standard CMS Formats
+#include "DataFormats/Common/interface/Ref.h"
+#include "DataFormats/Common/interface/Ptr.h"
+#include "DataFormats/Common/interface/DetSetVector.h"
+#include "DataFormats/Phase2TrackerDigi/interface/Phase2TrackerDigi.h"
+#include "DataFormats/L1TrackTrigger/interface/TTCluster.h"
+#include "DataFormats/L1TrackTrigger/interface/TTStub.h"
+
+/// The reference types
+typedef edm::Ref< edm::DetSetVector< Phase2TrackerDigi >, Phase2TrackerDigi > Ref_Phase2TrackerDigi_;
+
+#endif
+