/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_14/src/SimG4CMS/EcalTestBeam/interface/EcalTBMCInfoProducer.h

Go to the documentation of this file.

#ifndef EcalTBMCInfoProducer_H
#define EcalTBMCInfoProducer_H
/*
 * \file EcalTBMCInfoProducer.h
 *
 * $Id: EcalTBMCInfoProducer.h,v 1.9 2011/05/20 17:17:34 wmtan Exp $
 *
*/


#include "FWCore/Framework/interface/EDProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"
#include "SimDataFormats/EcalTestBeam/interface/PEcalTBInfo.h"
#include "Geometry/EcalTestBeam/interface/EcalTBCrystalMap.h"

#include "Math/GenVector/Rotation3D.h"

#include <iostream>
#include <fstream>
#include <vector>

namespace CLHEP {
  class RandFlat;
}

class EcalTBMCInfoProducer: public edm::EDProducer{
  
 public:
  
  EcalTBMCInfoProducer(const edm::ParameterSet& ps);
  
  virtual ~EcalTBMCInfoProducer();
  
  void produce(edm::Event & event, const edm::EventSetup& eventSetup);
  
  // BeginJob
  //void beginJob();
  
  // EndJob
  //void endJob(void);
  

private:

  double beamEta;
  double beamPhi;
  double beamTheta;

  int crysNumber;

  double beamXoff;
  double beamYoff;

  double partXhodo;
  double partYhodo;

  EcalTBCrystalMap * theTestMap;

  ROOT::Math::Rotation3D * fromCMStoTB;

  std::string GenVtxLabel;

  CLHEP::RandFlat *flatDistribution_;
 

};

#endif