/data/refman/pasoursint/CMSSW_5_3_8_patch3/src/Validation/RecoTrack/interface/MTVHistoProducerAlgo.h

Go to the documentation of this file.

#ifndef Validation_RecoTrack_MTVHistoProducerAlgo_h
#define Validation_RecoTrack_MTVHistoProducerAlgo_h

/* \author B.Mangano, UCSD
 *
 * Base class which defines the interface of a generic HistoProducerAlogs 
 * to be used within the MultiTrackValidator module. 
 * The concrete algorithms will be plugged into the MTV to produce all 
 * the validation plots that the user wants.
 */
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "SimDataFormats/TrackingAnalysis/interface/TrackingParticle.h"
#include "DataFormats/TrackReco/interface/Track.h"

#include "DQMServices/Core/interface/MonitorElement.h"
#include "DQMServices/Core/interface/DQMStore.h"

#include "DataFormats/TrackReco/interface/DeDxData.h"
#include "DataFormats/Common/interface/ValueMap.h"

#include <TH1F.h>
#include <TH2F.h>

class MTVHistoProducerAlgo{
 public:
  
 MTVHistoProducerAlgo(const edm::ParameterSet& pset) : pset_(pset){};
  virtual ~MTVHistoProducerAlgo() {}
  // to be implemented in the concrete classes
  virtual void initialize()=0; 

  void setDQMStore(DQMStore* dbe) {dbe_ = dbe;}

  virtual void bookSimHistos()=0;
  virtual void bookRecoHistos()=0;
  virtual void bookRecoHistosForStandaloneRunning()=0;

  virtual void fill_generic_simTrack_histos(int counter,ParticleBase::Vector,ParticleBase::Point vertex, int bx)=0;

  virtual void fill_recoAssociated_simTrack_histos(int count,
                                                   const TrackingParticle& tp,
                                                   ParticleBase::Vector momentumTP,ParticleBase::Point vertexTP,
                                                   double dxy, double dz, int nSimHits,
                                                   const reco::Track* track,
                           int numVertices, double vertz)=0;

  virtual void fill_generic_recoTrack_histos(int count,
                                             const reco::Track& track,
                                             math::XYZPoint bsPosition,
                                             bool isMatched,
                                             bool isSigMatched,
                                             bool isChargeMatched,
                             int numAssocRecoTracks,
                                             int numVertices, 
                                             int tpbunchcrossing,
                                             int nSimHits,
                                             double sharedFraction)=0;

  virtual void fill_dedx_recoTrack_histos(int count, edm::RefToBase<reco::Track>& trackref, std::vector< edm::ValueMap<reco::DeDxData> > v_dEdx)=0;
  //  virtual void fill_dedx_recoTrack_histos(reco::TrackRef trackref, std::vector< edm::ValueMap<reco::DeDxData> > v_dEdx)=0;

  virtual void fill_simAssociated_recoTrack_histos(int count,
                                                   const reco::Track& track)=0;

  virtual void fill_trackBased_histos(int count, 
                                      int assTracks, 
                                      int numRecoTracks,
                                      int numSimTracks)=0;

  virtual void fill_ResoAndPull_recoTrack_histos(int count,
                                                 ParticleBase::Vector momentumTP,
                                                 ParticleBase::Point vertexTP,
                                                 int chargeTP,
                                                 const reco::Track& track,
                                                 math::XYZPoint bsPosition)=0;

  virtual void finalHistoFits(int counter)=0;


  virtual void fillHistosFromVectors(int counter)=0;
  virtual void fillProfileHistosFromVectors(int counter)=0;


 protected: 
  //protected functions 

  virtual double getEta(double eta)=0;
  
  virtual double getPt(double pt)=0;

  void doProfileX(TH2 * th2, MonitorElement* me);

  void doProfileX(MonitorElement * th2m, MonitorElement* me) {
    doProfileX(th2m->getTH2F(), me);
  }
  
  void fillPlotFromVector(MonitorElement* h, std::vector<int>& vec);
  
  void fillPlotFromVectors(MonitorElement* h, 
                           std::vector<int>& numerator, 
                           std::vector<int>& denominator,
                           std::string type);

  void BinLogX(TH1*h);

  DQMStore* dbe_;
 private:
  //private data members
  const edm::ParameterSet& pset_;


};

#endif