#include <FamosManager.h>

Public Member Functions
CalorimetryManager *	calorimetryManager () const
	The calorimeter. More...

	FamosManager (edm::ParameterSet const &p)
	Constructor. More...

void	reconstruct (const HepMC::GenEvent evt, const reco::GenParticleCollection particles, const HepMC::GenEvent pu, const TrackerTopology tTopo, RandomEngineAndDistribution const *)
	The real thing is done here. More...

void	reconstruct (const reco::GenParticleCollection particles, const TrackerTopology tTopo, RandomEngineAndDistribution const *)

void	setupGeometryAndField (edm::Run const &run, const edm::EventSetup &es)
	Get information from the Event Setup. More...

FSimEvent *	simEvent () const
	The generated event. More...

TrajectoryManager *	trackerManager () const
	The tracker. More...

	~FamosManager ()
	Destructor. More...

Private Attributes
int	iEvent

bool	m_Alignment

bool	m_Calorimetry

int	m_pRunNumber

bool	m_pUseMagneticField

int	m_pVerbose

bool	m_Tracking

CalorimetryManager *	myCalorimetry

PileUpSimulator *	myPileUpSimulator

FSimEvent *	mySimEvent

TrajectoryManager *	myTrajectoryManager

double	weight_

Detailed Description

Definition at line 30 of file FamosManager.h.

Constructor & Destructor Documentation

FamosManager::FamosManager ( edm::ParameterSet const & p )

Constructor.

Initialize the TrajectoryManager

Definition at line 47 of file FamosManager.cc.

References edm::ParameterSet::getParameter(), m_Calorimetry, myCalorimetry, myPileUpSimulator, mySimEvent, and myTrajectoryManager.

     : iEvent(0),
       myPileUpSimulator(0),
       myCalorimetry(0),
       m_pUseMagneticField(p.getParameter<bool>("UseMagneticField")),
       m_Tracking(p.getParameter<bool>("SimulateTracking")),
       m_Calorimetry(p.getParameter<bool>("SimulateCalorimetry")),
       m_Alignment(p.getParameter<bool>("ApplyAlignment")),
       m_pRunNumber(p.getUntrackedParameter<int>("RunNumber",1)),
       m_pVerbose(p.getUntrackedParameter<int>("Verbosity",1))
 {
   // Initialize the FSimEvent
   mySimEvent = 
     new FSimEvent(p.getParameter<edm::ParameterSet>("ParticleFilter"));
 
   myTrajectoryManager = 
     new TrajectoryManager(mySimEvent,
               p.getParameter<edm::ParameterSet>("MaterialEffects"),
               p.getParameter<edm::ParameterSet>("TrackerSimHits"),
               p.getParameter<edm::ParameterSet>("ActivateDecays"));
 
   // Initialize PileUp Producer (if requested)
   myPileUpSimulator = new PileUpSimulator(mySimEvent);
 
   // Initialize Calorimetry Fast Simulation (if requested)
   if ( m_Calorimetry) 
     myCalorimetry = 
       new CalorimetryManager(mySimEvent,
                  p.getParameter<edm::ParameterSet>("Calorimetry"),               
                  p.getParameter<edm::ParameterSet>("MaterialEffectsForMuonsInECAL"),
                  p.getParameter<edm::ParameterSet>("MaterialEffectsForMuonsInHCAL"),
                              p.getParameter<edm::ParameterSet>("GFlash"));
 }

FamosManager::~FamosManager ( )

Destructor.

Definition at line 82 of file FamosManager.cc.

References myCalorimetry, myPileUpSimulator, mySimEvent, and myTrajectoryManager.

 { 
   if ( mySimEvent ) delete mySimEvent; 
   if ( myTrajectoryManager ) delete myTrajectoryManager; 
   if ( myPileUpSimulator ) delete myPileUpSimulator;
   if ( myCalorimetry) delete myCalorimetry;
 }

Member Function Documentation

CalorimetryManager* FamosManager::calorimetryManager ( ) const

inline

The calorimeter.

Definition at line 66 of file FamosManager.h.

References myCalorimetry.

Referenced by FamosProducer::produce().

66 {return myCalorimetry;}

FamosManager::myCalorimetry

CalorimetryManager * myCalorimetry

Definition: FamosManager.h:76

void FamosManager::reconstruct	(	const HepMC::GenEvent *	evt,
		const reco::GenParticleCollection *	particles,
		const HepMC::GenEvent *	pu,
		const TrackerTopology *	tTopo,
		RandomEngineAndDistribution const *	random
	)

The real thing is done here.

Definition at line 155 of file FamosManager.cc.

References FSimEvent::fill(), iEvent, m_pRunNumber, myCalorimetry, myPileUpSimulator, mySimEvent, myTrajectoryManager, FSimEvent::nGenParts(), FSimEvent::nTracks(), FSimEvent::nVertices(), PileUpSimulator::produce(), CalorimetryManager::reconstruct(), TrajectoryManager::reconstruct(), and FSimEvent::weight().

Referenced by FamosProducer::produce().

 {
   //  myGenEvent = evt;
   
   if (evt != 0 || particles != 0) {
     iEvent++;
     edm::EventID id(m_pRunNumber,1U,iEvent);
 
 
     // Fill the event from the original generated event
     if (evt ) 
       mySimEvent->fill(*evt,id);
     
     else 
       mySimEvent->fill(*particles,id);
     
     //    mySimEvent->printMCTruth(*evt);
     /*
       mySimEvent->print();
       std::cout << "----------------------------------------" << std::endl;
     */
     
     // Get the pileup events and add the particles to the main event
     myPileUpSimulator->produce(pu);
     /*
       mySimEvent->print();
     std::cout << "----------------------------------------" << std::endl;
     */
     
     // And propagate the particles through the detector
     myTrajectoryManager->reconstruct(tTopo, random);
     /*
       mySimEvent->print();
       std::cout << "=========================================" 
       << std::endl
       << std::endl;
     */
 
     if ( myCalorimetry ) myCalorimetry->reconstruct(random);
   }
 
   // Should be moved to LogInfo
   edm::LogInfo("FamosManager")  << " saved : Event  " << iEvent
                 << " of weight " << mySimEvent->weight()
                 << " with " << mySimEvent->nTracks()
                 << " tracks and " << mySimEvent->nVertices()
                 << " vertices, generated by "
                 << mySimEvent->nGenParts() << " particles " << std::endl;
 }

void FamosManager::reconstruct	(	const reco::GenParticleCollection *	particles,
		const TrackerTopology *	tTopo,
		RandomEngineAndDistribution const *	random
	)

Definition at line 209 of file FamosManager.cc.

References FSimEvent::fill(), iEvent, m_pRunNumber, myCalorimetry, mySimEvent, myTrajectoryManager, CalorimetryManager::reconstruct(), and TrajectoryManager::reconstruct().

                                                                                                                                                  {
   iEvent++;
   edm::EventID id(m_pRunNumber,1U,iEvent);
   mySimEvent->fill(*particles, id);
   myTrajectoryManager->reconstruct(tTopo, random);
   if ( myCalorimetry ) myCalorimetry->reconstruct(random);
 }

void FamosManager::setupGeometryAndField	(	edm::Run const &	run,
		const edm::EventSetup &	es
	)

Get information from the Event Setup.

Definition at line 91 of file FamosManager.cc.

References g, edm::EventSetup::get(), CalorimetryManager::getCalorimeter(), edm::EventSetup::getData(), CalorimetryManager::getHFShowerLibrary(), FastHFShowerLibrary::initHFShowerLibrary(), CaloGeometryHelper::initialize(), FBaseSimEvent::initializePdt(), TrajectoryManager::initializeRecoGeometry(), TrajectoryManager::initializeTrackerGeometry(), m_Alignment, m_pRunNumber, m_pUseMagneticField, m_Tracking, myCalorimetry, mySimEvent, myTrajectoryManager, edm::RunBase::run(), Calorimeter::setupGeometry(), Calorimeter::setupTopology(), AlCaHLTBitMon_QueryRunRegistry::string, and patCandidatesForDimuonsSequences_cff::tracker.

Referenced by FamosProducer::beginRun().

 {
   // Particle data table (from Pythia)
   edm::ESHandle < HepPDT::ParticleDataTable > pdt;
   es.getData(pdt);
   mySimEvent->initializePdt(&(*pdt));
   
   // Initialize the full (misaligned) tracker geometry 
   // (only if tracking is requested)
   std::string misAligned = m_Alignment ? "MisAligned" : "";
   // 1) By default, the aligned geometry is chosen (m_Alignment = false)
   // 2) By default, the misaligned geometry is aligned
   edm::ESHandle<TrackerGeometry> tracker;
   es.get<TrackerDigiGeometryRecord>().get(misAligned,tracker);
   if (m_Tracking)  myTrajectoryManager->initializeTrackerGeometry(&(*tracker)); 
 
   // Initialize the tracker misaligned reco geometry (always needed)
   // By default, the misaligned geometry is aligned
   edm::ESHandle<GeometricSearchTracker>       theGeomSearchTracker;
   es.get<TrackerRecoGeometryRecord>().get(misAligned, theGeomSearchTracker );
 
   // Initialize the misaligned tracker interaction geometry 
   edm::ESHandle<TrackerInteractionGeometry>  theTrackerInteractionGeometry;
   es.get<TrackerInteractionGeometryRecord>().get(misAligned, theTrackerInteractionGeometry );
 
   // Initialize the magnetic field
   double bField000 = 0.;
   if (m_pUseMagneticField) {
     edm::ESHandle<MagneticFieldMap> theMagneticFieldMap;
     es.get<MagneticFieldMapRecord>().get(misAligned, theMagneticFieldMap);
     const GlobalPoint g(0.,0.,0.);
     bField000 = theMagneticFieldMap->inTeslaZ(g);
     myTrajectoryManager->initializeRecoGeometry(&(*theGeomSearchTracker),
                         &(*theTrackerInteractionGeometry),
                         &(*theMagneticFieldMap));
   } else { 
     myTrajectoryManager->initializeRecoGeometry(&(*theGeomSearchTracker),
                         &(*theTrackerInteractionGeometry),
                         0);
     bField000 = 4.0;
   }
   // The following should be on LogInfo
   //std::cout << "B-field(T) at (0,0,0)(cm): " << bField000 << std::endl;      
     
   //  Initialize the calorimeter geometry
   if ( myCalorimetry ) {
     edm::ESHandle<CaloGeometry> pG;
     es.get<CaloGeometryRecord>().get(pG);   
     myCalorimetry->getCalorimeter()->setupGeometry(*pG);
 
     edm::ESHandle<CaloTopology> theCaloTopology;
     es.get<CaloTopologyRecord>().get(theCaloTopology);     
     myCalorimetry->getCalorimeter()->setupTopology(*theCaloTopology);
     myCalorimetry->getCalorimeter()->initialize(bField000);
 
     myCalorimetry->getHFShowerLibrary()->initHFShowerLibrary(es);
   }
 
   m_pRunNumber = run.run();
 
 }