ZDCDigiStudy.h

Go to the documentation of this file.

// Package:    ZDCDigiStudy
// Class:      ZDCDigiStudy
//
/*
 Description: 
              This code has been developed to be a check for the ZDC sim. In 2009, it was found that the ZDC Simulation was unrealistic and needed repair. The aim of this code is to show the user the input and output of a ZDC MinBias simulation.

 Implementation:
      First a MinBias simulation should be run, it could be pythia,hijin,or hydjet. This will output a .root file which should have information about recoGenParticles, hcalunsuppresseddigis. Use this .root file as the input into the cfg.py which is found in the main directory of this package. This output will be another .root file which is meant to be viewed in a TBrowser.

*/
//
// Original Author: Jaime Gomez (U. of Maryland) with SIGNIFICANT assistance of Dr. Jefferey Temple (U. of Maryland) 
//   
//
//         Created:  Summer 2012



#ifndef SimG4CMS_ZDCDigiStudy_H
#define SimG4CMS_ZDCDigiStudy_H

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"
#include "FWCore/ServiceRegistry/interface/Service.h"


#include "DataFormats/HcalDigi/interface/HcalDigiCollections.h"

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

class ZDCDigiStudy: public edm::EDAnalyzer{
public:

  ZDCDigiStudy(const edm::ParameterSet& ps);
  ~ZDCDigiStudy();

protected:

  void beginJob ();
  void endJob   ();
  void endRun(const edm::Run& run, const edm::EventSetup& c);

  void analyze  (const edm::Event& e, const edm::EventSetup& c);
  int FillHitValHist (int side,int section,int channel,double energy,double time);

private:

  //#   Below all the monitoring elements #
  //# are simply the plots "code names"  #
  //# they will be filled in the .cc file #



  std::string    zdcHits, outFile_;
  bool           verbose_, checkHit_;
  DQMStore       *dbe_;

  edm::EDGetTokenT<ZDCDigiCollection> tok_zdc_;



 MonitorElement* meZdcfCPHAD;
 MonitorElement* meZdcfCPTOT;
 MonitorElement* meZdcfCNHAD;
 MonitorElement* meZdcfCNTOT;
 MonitorElement* meZdcfCPEMvHAD;
 MonitorElement* meZdcfCNEMvHAD;
 MonitorElement* meZdcPEM1fCvsTS;
 MonitorElement* meZdcPEM2fCvsTS;
 MonitorElement* meZdcPEM3fCvsTS;
 MonitorElement* meZdcPEM4fCvsTS;
 MonitorElement* meZdcPEM5fCvsTS;
 MonitorElement* meZdcPHAD1fCvsTS;
 MonitorElement* meZdcPHAD2fCvsTS;
 MonitorElement* meZdcPHAD3fCvsTS;
 MonitorElement* meZdcPHAD4fCvsTS;
 MonitorElement* meZdcNEM1fCvsTS;
 MonitorElement* meZdcNEM2fCvsTS;
 MonitorElement* meZdcNEM3fCvsTS;
 MonitorElement* meZdcNEM4fCvsTS;
 MonitorElement* meZdcNEM5fCvsTS;
 MonitorElement* meZdcNHAD1fCvsTS;
 MonitorElement* meZdcNHAD2fCvsTS;
 MonitorElement* meZdcNHAD3fCvsTS;
 MonitorElement* meZdcNHAD4fCvsTS;



};

#endif