ValidationMisalignedTracker.h

Go to the documentation of this file.

#ifndef ValidationMisalignedTracker_h
#define ValidationMisalignedTracker_h

// system include files
#include <memory>

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

// 
#include "FWCore/Framework/interface/ESHandle.h" 
#include "SimTracker/Records/interface/TrackAssociatorRecord.h"
#include "SimTracker/TrackAssociation/interface/TrackAssociatorBase.h"
#include "SimDataFormats/TrackingAnalysis/interface/TrackingParticleFwd.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"
#include "SimDataFormats/TrackingAnalysis/interface/TrackingVertex.h"
#include "DataFormats/GeometryVector/interface/GlobalVector.h" 
#include "MagneticField/Engine/interface/MagneticField.h"

#include "TTree.h"
#include "TFile.h"

#include <TROOT.h>
#include <TSystem.h>
#include <TCanvas.h>
#include <cmath>
#include "TStyle.h"

//
// class decleration
//

class ValidationMisalignedTracker : public edm::EDAnalyzer {
public:
  explicit ValidationMisalignedTracker(const edm::ParameterSet&);
  ~ValidationMisalignedTracker();
  
  
private:
  virtual void analyze(const edm::Event&, const edm::EventSetup&);
  virtual void endJob() ;
  
  // ----------member data ---------------------------
  
  std::string simobject,trackassociator;
  bool selection_eff,selection_fake,ZmassSelection_;
  std::string rootfile_;
  
  bool skip;
  int eventCount_;
  TTree* tree_eff;
  TTree* tree_fake;
  TFile* file_;
  
  
  int irun, ievt;
 
  float mzmu,recmzmu,ptzmu,recptzmu,etazmu,thetazmu,phizmu,recetazmu,recthetazmu,recphizmu;
  float recenezmu, enezmu, pLzmu, recpLzmu,yzmu,recyzmu,mxptmu,recmxptmu,minptmu,recminptmu;
  int countpart[2],countpartrec[2];
  int flag,flagrec,count,countrec;
  // int countsim;
  float ene[2][2],p[2][2],px[2][2],py[2][2],pz[2][2],ptmu[2][2];
  float recene[2][2],recp[2][2],recpx[2][2],recpy[2][2],recpz[2][2],recptmu[2][2];
  
  int  trackType;
  float pt, eta, cottheta, theta, costheta, phi, d0, z0; 
  int nhit;
  float recpt, receta, rectheta, reccottheta, recphi, recd0, recz0;
  int nAssoc, recnhit;
  float recchiq;
  float reseta, respt, resd0, resz0, resphi,rescottheta,eff;
  
  float fakemzmu,fakerecmzmu,fakeptzmu,fakerecptzmu,fakeetazmu,fakethetazmu,fakephizmu,fakerecetazmu,fakerecthetazmu,fakerecphizmu;
  float fakerecenezmu, fakeenezmu, fakepLzmu, fakerecpLzmu,fakeyzmu,fakerecyzmu,fakemxptmu,fakerecmxptmu,fakeminptmu,fakerecminptmu;
  int fakecountpart[2],fakecountpartrec[2],fakeflag,fakeflagrec,fakecount,fakecountsim,fakecountrec;
  float fakeene[2][2],fakep[2][2],fakepx[2][2],fakepy[2][2],fakepz[2][2],fakeptmu[2][2];
  float fakerecene[2][2],fakerecp[2][2],fakerecpx[2][2],fakerecpy[2][2],fakerecpz[2][2],fakerecptmu[2][2];
  
  int  faketrackType;
  float fakept, fakeeta, fakecottheta, faketheta, fakecostheta, fakephi, faked0, fakez0; 
  int fakenhit;
  float fakerecpt, fakereceta, fakerectheta, fakereccottheta, fakerecphi, fakerecd0, fakerecz0;
  int fakenAssoc, fakerecnhit;
  float fakerecchiq;
  float fakereseta, fakerespt, fakeresd0, fakeresz0, fakeresphi,fakerescottheta,fake;
  
  double chi2tmp;
  float fractiontmp;
  bool onlyDiag;
  edm::ESHandle<MagneticField> theMF;
  std::vector<std::string> associators;
  std::vector<const TrackAssociatorBase*> associatore;
 
  std::vector<edm::InputTag> label;
  edm::InputTag label_tp_effic;
  edm::InputTag label_tp_fake;

  GlobalVector magField;
  std::vector<float> ptused;

  edm::ESWatcher<TrackAssociatorRecord> watchTrackAssociatorRecord_;
};

#endif