Analyzer_minbias.h

Go to the documentation of this file.

#ifndef PhysicsToolsAnalysisAnalyzerMinBias_h
#define PhysicsToolsAnalysisAnalyzerMinBias_h

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

// 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/Framework/interface/EDProducer.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/Common/interface/Ref.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/CaloTowers/interface/CaloTowerCollection.h"
#include "DataFormats/Candidate/interface/Candidate.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/JetReco/interface/Jet.h"
#include "DataFormats/JetReco/interface/CaloJet.h"
#include "DataFormats/JetReco/interface/CaloJetCollection.h"
#include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
#include "Geometry/CaloGeometry/interface/CaloGeometry.h"
#include "Geometry/Records/interface/CaloGeometryRecord.h"
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "DataFormats/CaloTowers/interface/CaloTowerDetId.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/Run.h"
#include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutRecord.h"
//#include "CalibFormats/HcalObjects/interface/HcalDbService.h"
#include "TFile.h"
#include "TH1.h"
#include "TH2.h"
#include "TTree.h"
#include "DataFormats/HcalDetId/interface/HcalDetId.h"
#include <vector>
#include <map>
//#include "CalibCalorimetry/CaloMiscalibTools/interface/CaloMiscalibMapHcal.h"
#include <string>
#include <iostream>
#include <fstream>
#include <sstream>

//
// class declaration
//
namespace cms{
class Analyzer_minbias : public edm::EDAnalyzer {
   public:
      explicit Analyzer_minbias(const edm::ParameterSet&);
      ~Analyzer_minbias();

      virtual void analyze(const edm::Event&, const edm::EventSetup&);
      virtual void beginJob() ;
      virtual void endJob() ;
      virtual void beginRun( const edm::Run& r, const edm::EventSetup& iSetup);
      virtual void endRun( const edm::Run& r, const edm::EventSetup& iSetup);

   private:
  // ----------member data ---------------------------
     std::string fOutputFileName ;
     std::string hcalfile_;
     std::ofstream *myout_hcal;
     
  // names of modules, producing object collections
  
     TFile*      hOutputFile ;
     TTree*      myTree;
     TH1F*       hCalo1[73][43];
     TH1F*       hCalo2[73][43];
     TH1F*       hCalo1mom2[73][43];
     TH1F*       hCalo2mom2[73][43];
     TH1F*       hbheNoiseE;
     TH1F*       hbheSignalE;
     TH1F*       hfNoiseE;
     TH1F*       hfSignalE;
     
     TH2F*       hHBHEsize_vs_run; 
     TH2F*       hHFsize_vs_run;
  // Root tree members
     int nevent_run;   
     int mydet, mysubd, depth, iphi, ieta;
     float phi,eta;
     float mom0_MB,mom1_MB,mom2_MB,mom3_MB,mom4_MB,occup;
     float mom0_Noise,mom1_Noise,mom2_Noise,mom3_Noise,mom4_Noise;
     float mom0_Diff,mom1_Diff,mom2_Diff,mom3_Diff,mom4_Diff;
  
  // Noise subtraction
     
     double meannoise_pl[73][43],meannoise_min[73][43]; 
     double noise_pl[73][43],noise_min[73][43];

  // counters

     double nevent;
     double theMBFillDetMapPl0[5][5][73][43]; 
     double theMBFillDetMapPl1[5][5][73][43]; 
     double theMBFillDetMapPl2[5][5][73][43];
     double theMBFillDetMapPl4[5][5][73][43];

     double theMBFillDetMapMin0[5][5][73][43]; 
     double theMBFillDetMapMin1[5][5][73][43]; 
     double theMBFillDetMapMin2[5][5][73][43];
     double theMBFillDetMapMin4[5][5][73][43];

     double theNSFillDetMapPl0[5][5][73][43]; 
     double theNSFillDetMapPl1[5][5][73][43]; 
     double theNSFillDetMapPl2[5][5][73][43];
     double theNSFillDetMapPl4[5][5][73][43];

     double theNSFillDetMapMin0[5][5][73][43]; 
     double theNSFillDetMapMin1[5][5][73][43]; 
     double theNSFillDetMapMin2[5][5][73][43];
     double theNSFillDetMapMin4[5][5][73][43];

     double theDFFillDetMapPl0[5][5][73][43]; 
     double theDFFillDetMapPl1[5][5][73][43]; 
     double theDFFillDetMapPl2[5][5][73][43];
     double theDFFillDetMapMin0[5][5][73][43]; 
     double theDFFillDetMapMin1[5][5][73][43]; 
     double theDFFillDetMapMin2[5][5][73][43];
     

    edm::EDGetTokenT<HBHERecHitCollection> tok_hbhe_;
    edm::EDGetTokenT<HORecHitCollection> tok_ho_;
    edm::EDGetTokenT<HFRecHitCollection> tok_hf_;
  
    edm::EDGetTokenT<HBHERecHitCollection> tok_hbheNoise_; 
    edm::EDGetTokenT<HORecHitCollection> tok_hoNoise_;
    edm::EDGetTokenT<HFRecHitCollection> tok_hfNoise_;

  // 
  edm::EDGetTokenT<L1GlobalTriggerReadoutRecord> tok_gtRec_;
  edm::EDGetTokenT<HBHERecHitCollection> tok_hbheNorm_;
  
     bool theRecalib;

};
}
#endif