TrackParameterAnalyzer.h

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    TrackParameterAnalyzer
// Class:      TrackParameterAnalyzer
// 
//
// Original Author:  Wolfram Erdmann
//         Created:  Fri Jun  2 10:54:05 CEST 2006
// $Id: TrackParameterAnalyzer.h,v 1.12 2009/10/30 19:55:01 drell Exp $
//
//


// system include files
#include <memory>
#include <string>
#include <vector>

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

// Hep MC stuff from CLHEP, add  <use name=clhep> to the buildfile
#include "HepMC/GenEvent.h"
#include "HepMC/GenVertex.h"
#include "HepMC/GenParticle.h"
#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"
#include "CLHEP/Vector/LorentzVector.h"

// simulated vertices,..., add <use name=SimDataFormats/Vertex> and <../Track>
#include <SimDataFormats/Vertex/interface/SimVertex.h>
#include <SimDataFormats/Vertex/interface/SimVertexContainer.h>
#include <SimDataFormats/Track/interface/SimTrack.h>
#include <SimDataFormats/Track/interface/SimTrackContainer.h>

// vertex stuff
#include <DataFormats/VertexReco/interface/Vertex.h>
#include "RecoVertex/VertexPrimitives/interface/TransientVertex.h"
// perigee 
#include <TrackingTools/TrajectoryParametrization/interface/PerigeeTrajectoryParameters.h>
//#include <TrackingTools/TrajectoryState/interface/TrajectoryStateClosestToPoint.h>

// Root
#include <TH1.h>
#include <TH2.h>
#include <TFile.h>


// class declaration
//
typedef reco::TrackBase::ParameterVector ParameterVector;

class TrackParameterAnalyzer : public edm::EDAnalyzer {
   public:
      explicit TrackParameterAnalyzer(const edm::ParameterSet&);
      ~TrackParameterAnalyzer();

      virtual void analyze(const edm::Event&, const edm::EventSetup&);
      virtual void beginJob();
      virtual void endJob();

   private:
       bool match(const ParameterVector  &a, const ParameterVector &b);
      // ----------member data ---------------------------
      std::string recoTrackProducer_;
      // root file to store histograms
      std::string outputFile_; // output file
      TFile*  rootFile_;
      TH1*   h1_pull0_; 
      TH1*   h1_pull1_;
      TH1*   h1_pull2_;
      TH1*   h1_pull3_;
      TH1*   h1_pull4_;
      TH1*   h1_res0_; 
      TH1*   h1_res1_;
      TH1*   h1_res2_;
      TH1*   h1_res3_;
      TH1*   h1_res4_;
      TH1*   h1_Beff_;
      TH2*   h2_dvsphi_;
      TH1*   h1_par0_; 
      TH1*   h1_par1_;
      TH1*   h1_par2_;
      TH1*   h1_par3_;
      TH1*   h1_par4_;

      edm::InputTag simG4_;
      bool verbose_;
      double simUnit_;               
};