NuclearInteractionEDProducer.h

Go to the documentation of this file.

#ifndef CD_NuclearInteractionEDProducer_H_
#define CD_NuclearInteractionEDProducer_H_
// -*- C++ -*-
//
// Package:    NuclearAssociatonMapEDProducer
// Class:      NuclearInteractionEDProducer
//
//
// Original Author:  Vincent ROBERFROID
//         Created:  Fri Aug 10 12:05:36 CET 2007
// $Id: NuclearInteractionEDProducer.h,v 1.10 2010/01/21 10:53:40 adamwo Exp $
//
//


// system include files
#include <memory>

// user include files
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/ESWatcher.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDProducer.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"

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

#include "DataFormats/TrackReco/interface/TrackFwd.h"

#include "MagneticField/Engine/interface/MagneticField.h"
#include "MagneticField/Records/interface/IdealMagneticFieldRecord.h"
#include "RecoTracker/NuclearSeedGenerator/interface/TrajectoryToSeedMap.h"
#include "TrackingTools/Records/interface/TransientTrackRecord.h"

#include "DataFormats/VertexReco/interface/NuclearInteraction.h"

class NuclearVertexBuilder;
class NuclearLikelihood;

class NuclearInteractionEDProducer : public edm::EDProducer {

public:
      typedef edm::RefVector<TrajectorySeedCollection> TrajectorySeedRefVector;

      explicit NuclearInteractionEDProducer(const edm::ParameterSet&);
      ~NuclearInteractionEDProducer();

   private:
      virtual void beginJob() ;
      virtual void produce(edm::Event&, const edm::EventSetup&);
      virtual void endJob();

      bool isInside( const reco::TrackRef& track, const TrajectorySeedRefVector& seeds);
      void findAdditionalSecondaryTracks( reco::NuclearInteraction& nucl,
                                          const edm::Handle<reco::TrackCollection>& additionalSecTracks);

      // ----------member data ---------------------------
      edm::ParameterSet conf_;
      std::string primaryProducer_;
      std::string seedsProducer_;
      std::string secondaryProducer_;
      std::string additionalSecondaryProducer_;

      std::auto_ptr< NuclearVertexBuilder >  vertexBuilder;
      std::auto_ptr< NuclearLikelihood >     likelihoodCalculator;

  edm::ESWatcher<IdealMagneticFieldRecord> magFieldWatcher_;
  edm::ESWatcher<TransientTrackRecord> transientTrackWatcher_;

};

void print(std::ostringstream& str, const reco::NuclearInteraction& nucl, const std::auto_ptr< NuclearVertexBuilder >& builder);
#endif