#ifndef DTANALYZER_H
#define DTANALYZER_H

/** \class DTAnalyzer
 *
 * Description:
 *  
 *  detailed description
 *
 * \author : Stefano Lacaprara - INFN Padova <stefano.lacaprara@pd.infn.it>
 * $date   : 20/11/2006 16:51:04 CET $
 *
 * Modification:
 *
 */

/* Base Class Headers */
#include "FWCore/Framework/interface/EDAnalyzer.h"
namespace edm {
  class ParameterSet;
  class Event;
  class EventSetup;
}

/* Collaborating Class Declarations */
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/ESHandle.h"
class TFile;
class TH1F;
class TH2F;
class DTLayerId;
class DTSuperLayerId;
class DTChamberId;
#include "Geometry/DTGeometry/interface/DTGeometry.h"
#include "DataFormats/DTRecHit/interface/DTRecSegment4DCollection.h"
class DTTTrigBaseSync;
class Propagator;

/* C++ Headers */
#include <iosfwd>
#include <bitset>

/* ====================================================================== */

/* Class DTAnalyzer Interface */

class DTAnalyzer : public edm::EDAnalyzer {

  public:

/* Constructor */ 
    DTAnalyzer(const edm::ParameterSet& pset) ;

/* Destructor */ 
    ~DTAnalyzer() ;

/* Operations */ 
    void analyze(const edm::Event & event, const edm::EventSetup& eventSetup);

  private:
    void analyzeDTHits(const edm::Event & event, const edm::EventSetup& eventSetup);
    void analyzeDTSegments(const edm::Event & event, const edm::EventSetup& eventSetup);
    void analyzeSATrack(const edm::Event & event, const edm::EventSetup& eventSetup);
    void analyzeTrigger(const edm::Event & event, const edm::EventSetup& eventSetup);
    void LabelTriggerMatrix();

    const DTRecSegment4D& getBestSegment(const DTRecSegment4DCollection::range& segs) const;
    const DTRecSegment4D* getBestSegment(const DTRecSegment4D* s1,
                                         const DTRecSegment4D* s2) const;
    LocalPoint interpolate(const DTRecSegment4D& seg1,
                           const DTRecSegment4D& seg3,
                           const DTChamberId& MB2) const;
    LocalPoint extrapolate(const DTRecSegment4D& seg1,
                           const DTChamberId& MB2) const;
    
    TH1F* histo(const std::string& name) const;
    TH2F* histo2d(const std::string& name) const;

    void createTH1F(const std::string& name,
                    const std::string& title,
                    const std::string& suffix,
                    int nbin, const double& binMin, const double& binMax) const;

    void createTH2F(const std::string& name,
                    const std::string& title,
                    const std::string& suffix,
                    int nBinX,
                    const double& binXMin,
                    const double& binXMax,
                    int nBinY,
                    const double& binYMin,
                    const double& binYMax) const ;

    enum LCTType { DT, CSC, RPC_W1, RPC_W2 };
    bool getLCT(LCTType) const;
    bool selectEvent() const ;

    std::string toString(const DTLayerId& id) const;
    std::string toString(const DTSuperLayerId& id) const;
    std::string toString(const DTChamberId& id) const;
    template<class T> std::string hName(const std::string& s, const T& id) const;
  private:
    bool LCT_RPC, LCT_DT, LCT_CSC;
    bool debug;
    int _ev; 
    std::string theRootFileName;
    TFile* theFile;
    //static std::string theAlgoName;
    int sects[12];
    std::string theDTLocalTriggerLabel;
    std::string theRecHits4DLabel;
    std::string theRecHits2DLabel;     
    std::string theRecHits1DLabel;     
    std::string theSTAMuonLabel;
    std::string thePropagatorName;
 
    mutable Propagator* thePropagator;

    bool mc;

    bool doHits;
    bool doSegs;
    bool doSA;
    bool doTrig;
    
    std::bitset<6> LCT;

    DTTTrigBaseSync *theSync;

    edm::ESHandle<DTGeometry> dtGeom;
    edm::Handle<DTRecSegment4DCollection> segs;

  protected:

};
#endif // DTANALYZER_H