CSCSkim.h

Go to the documentation of this file.

#ifndef DPGAnalysisSkims_CSCSkim_H
#define DPGAnalysisSkims_CSCSkim_H

#include <memory>
#include <iostream>
#include <vector>
#include <map>
#include <string>
#include <iomanip>
#include <fstream>

#include "FWCore/Framework/interface/ConsumesCollector.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/EDFilter.h"
#include "FWCore/Framework/interface/Event.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "DataFormats/MuonDetId/interface/CSCDetId.h"
#include "DataFormats/CSCDigi/interface/CSCWireDigi.h"
#include "DataFormats/CSCDigi/interface/CSCWireDigiCollection.h"
#include "DataFormats/CSCDigi/interface/CSCStripDigi.h"
#include "DataFormats/CSCDigi/interface/CSCStripDigiCollection.h"
#include <DataFormats/CSCRecHit/interface/CSCRecHit2D.h>
#include "DataFormats/CSCRecHit/interface/CSCRecHit2DCollection.h"
#include <DataFormats/CSCRecHit/interface/CSCSegmentCollection.h>

#include "DataFormats/DTRecHit/interface/DTRecSegment4D.h"
#include "DataFormats/DTRecHit/interface/DTRecSegment4DCollection.h"

#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"
#include "Geometry/Records/interface/GlobalTrackingGeometryRecord.h"
#include "Geometry/CommonDetUnit/interface/GlobalTrackingGeometry.h"

#include "Geometry/DTGeometry/interface/DTGeometry.h"
#include "Geometry/CSCGeometry/interface/CSCGeometry.h"
#include "Geometry/CSCGeometry/interface/CSCChamber.h"
#include "Geometry/CSCGeometry/interface/CSCLayer.h"
#include "Geometry/CSCGeometry/interface/CSCLayerGeometry.h"
#include "Geometry/Records/interface/MuonGeometryRecord.h"

#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "DataFormats/GeometryVector/interface/GlobalVector.h"
#include "DataFormats/GeometryVector/interface/LocalPoint.h"
#include "DataFormats/GeometryVector/interface/LocalVector.h"
#include "DataFormats/CLHEP/interface/AlgebraicObjects.h"

// #include "SimDataFormats/TrackingHit/interface/PSimHitContainer.h"

#include "TVector3.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TFile.h"
#include "TString.h"
#include "TTree.h"


class CSCSkim : public edm::EDFilter {
 public:
  // Constructor
  explicit CSCSkim(const edm::ParameterSet& pset);
   
  // Destructor
  ~CSCSkim() override;
     
  // Analysis routines
  void beginJob() override;
  bool filter(edm::Event& event, const edm::EventSetup& eventSetup) override;
  void endJob() override;
       
protected:
       
private: 

  // main skimming routine
  bool doCSCSkimming(edm::Handle<CSCRecHit2DCollection> cscRecHits, 
                     edm::Handle<CSCSegmentCollection> cscSegments);

  // extra skimming routine for alignment studies
  bool doOverlapSkimming(edm::Handle<CSCSegmentCollection> cscSegments);

  // skimming routine for messy events
  bool doMessyEventSkimming(edm::Handle<CSCRecHit2DCollection> cscRecHits, 
                            edm::Handle<CSCSegmentCollection> cscSegments);

  // select events with DIGIs in a certain chamber
  bool doCertainChamberSelection(edm::Handle<CSCWireDigiCollection> wires,
                 edm::Handle<CSCStripDigiCollection> strips);

  // select events which might probe the DT-CSC overlap region
  bool doDTOverlap(edm::Handle<CSCSegmentCollection> cscSegments);

  // select muons which go through inner part of stations 1,2,3,4
  bool doHaloLike(edm::Handle<CSCSegmentCollection> cscSegments);

  // select events with long stand-alone (cosmic) muons
  bool doLongSATrack(edm::Handle<reco::TrackCollection> saTracks);

  // select events suitable for magnetic field studies - they have a track in the tracker
  bool doBFieldStudySelection(edm::Handle<reco::TrackCollection> saTracks, edm::Handle<reco::TrackCollection> Tracks, edm::Handle<reco::MuonCollection> gMuons);

  // some useful functions
  int    chamberSerial(int kE, int kS, int kR, int kCh);

  // counters
  int nEventsAnalyzed;
  int nEventsSelected;
  int nEventsChambersBothSides;
  int nEventsOverlappingChambers;
  int nEventsMessy;
  int nEventsCertainChamber;
  int nEventsDTOverlap;
  int nEventsHaloLike;
  int nEventsLongSATrack;
  int nEventsForBFieldStudies;

  // run and event number
  int iRun;
  int iEvent;

  // The root file for the histograms.
  TFile *theHistogramFile;

  // file names
  std::string outputFileName;
  std::string histogramFileName;

  // token names
  edm::EDGetTokenT<CSCWireDigiCollection> wds_token;
  edm::EDGetTokenT<CSCStripDigiCollection> sds_token;
  edm::EDGetTokenT<CSCWireDigiCollection> wdr_token;
  edm::EDGetTokenT<CSCStripDigiCollection> sdr_token;

  edm::EDGetTokenT<CSCRecHit2DCollection> rh_token;
  edm::EDGetTokenT<CSCSegmentCollection> seg_token;
  edm::EDGetTokenT<reco::TrackCollection> sam_token;
  edm::EDGetTokenT<reco::TrackCollection> trk_token;
  edm::EDGetTokenT<reco::MuonCollection> glm_token;

  // parameters for the selection
  bool isSimulation;
  int typeOfSkim;
  int nLayersWithHitsMinimum;
  int minimumHitChambers;
  int minimumSegments;
  bool demandChambersBothSides;
  bool makeHistograms;
  bool makeHistogramsForMessyEvents;
  int whichEndcap;
  int whichStation;
  int whichRing;
  int whichChamber;
  // parameters for B-field study skim
  float pMin;
  float zLengthMin;
  int nCSCHitsMin;
  float zInnerMax;
  int nTrHitsMin;
  float zLengthTrMin;
  float rExtMax;
  float redChiSqMax;
  int nValidHitsMin;


  // histograms for skimming module
  TH1F *hxnRecHits;
  TH1F *hxnSegments;
  TH1F *hxnHitChambers;
  TH1F *hxnRecHitsSel;
  TH1F *mevnRecHits0;
  TH1F *mevnChambers0;
  TH1F *mevnSegments0;
  TH1F *mevnRecHits1;
  TH1F *mevnChambers1;
  TH1F *mevnSegments1;

  TH1F *xxP, *xxnValidHits, *xxnCSCHits, *xxredChiSq, *xxnTrackerHits;

};
#endif