#include <EcalTBMCInfoProducer.h>

Inheritance diagram for EcalTBMCInfoProducer:

Public Member Functions
	EcalTBMCInfoProducer (const edm::ParameterSet &ps)
	Constructor. More...

void	produce (edm::Event &event, const edm::EventSetup &eventSetup)
	Produce digis out of raw data. More...

virtual	~EcalTBMCInfoProducer ()
	Destructor. More...

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndex	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
double	beamEta

double	beamPhi

double	beamTheta

double	beamXoff

double	beamYoff

int	crysNumber

CLHEP::RandFlat *	flatDistribution_

ROOT::Math::Rotation3D *	fromCMStoTB

std::string	GenVtxLabel

double	partXhodo

double	partYhodo

EcalTBCrystalMap *	theTestMap

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

typedef WorkerT< EDProducer >	WorkerType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDProducer
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 34 of file EcalTBMCInfoProducer.h.

Constructor & Destructor Documentation

EcalTBMCInfoProducer::EcalTBMCInfoProducer ( const edm::ParameterSet & ps )

Constructor.

Definition at line 19 of file EcalTBMCInfoProducer.cc.

References beamEta, beamPhi, beamTheta, beamXoff, beamYoff, funct::cos(), crysNumber, reco_application_tbsim_DetSim-Digi_cfg::CrystalMapFile, HLTFastRecoForTau_cff::deltaEta, SiPixelRawToDigiRegional_cfi::deltaPhi, eta(), edm::hlt::Exception, create_public_lumi_plots::exp, EcalTBCrystalMap::findCrystalAngles(), flatDistribution_, fromCMStoTB, GenVtxLabel, edm::RandomNumberGenerator::getEngine(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), edm::Service< T >::isAvailable(), EcalTBCrystalMap::NCRYSTAL, phi, funct::sin(), theTestMap, and create_public_lumi_plots::xy.

                                                                     : flatDistribution_(0) {
   
   produces<PEcalTBInfo>();
 
   edm::FileInPath CrystalMapFile = ps.getParameter<edm::FileInPath>("CrystalMapFile");
   GenVtxLabel = ps.getUntrackedParameter<string>("moduleLabelVtx","source");
   double fMinEta = ps.getParameter<double>("MinEta");
   double fMaxEta = ps.getParameter<double>("MaxEta");
   double fMinPhi = ps.getParameter<double>("MinPhi");
   double fMaxPhi = ps.getParameter<double>("MaxPhi");
   beamEta = (fMaxEta+fMinEta)/2.;
   beamPhi = (fMaxPhi+fMinPhi)/2.;
   beamTheta = 2.0*atan(exp(-beamEta));
   beamXoff = ps.getParameter<double>("BeamMeanX");
   beamYoff = ps.getParameter<double>("BeamMeanX");
    
   string fullMapName = CrystalMapFile.fullPath();
   theTestMap = new EcalTBCrystalMap(fullMapName);
   crysNumber = 0;
 
   double deltaEta = 999.;
   double deltaPhi = 999.;
   for ( int cryIndex = 1; cryIndex <= EcalTBCrystalMap::NCRYSTAL; ++cryIndex) {
     double eta = 0;
     double phi = 0.;
     theTestMap->findCrystalAngles(cryIndex, eta, phi);
     if ( fabs(beamEta - eta) < deltaEta && fabs(beamPhi - phi) < deltaPhi ) {
       deltaEta = fabs(beamEta - eta);
       deltaPhi = fabs(beamPhi - phi);
       crysNumber = cryIndex;
     }
     else if (fabs(beamEta - eta)<deltaEta && fabs(beamPhi - phi)>deltaPhi ) {
       if ( fabs(beamPhi - phi) < 0.017 ) {
         deltaEta = fabs(beamEta - eta);
         deltaPhi = fabs(beamPhi - phi);
         crysNumber = cryIndex;
       }
     }
     else if (fabs(beamEta - eta)>deltaEta && fabs(beamPhi - phi)<deltaPhi ) {
       if ( fabs(beamEta - eta) < 0.017 ) {
         deltaEta = fabs(beamEta - eta);
         deltaPhi = fabs(beamPhi - phi);
         crysNumber = cryIndex;
       }
     }
   }
 
   edm::LogInfo("EcalTBInfo") << "Initialize TB MC ECAL info producer with parameters: \n"
                              << "Crystal map file:  " << CrystalMapFile << "\n"
                              << "Beam average eta = " << beamEta << "\n"
                              << "Beam average phi = " << beamPhi << "\n"
                              << "Corresponding to crystal number = " << crysNumber << "\n"
                              << "Beam X offset =    " << beamXoff << "\n"
                              << "Beam Y offset =    " << beamYoff;
 
   // rotation matrix to move from the CMS reference frame to the test beam one
 
   double xx = -cos(beamTheta)*cos(beamPhi);
   double xy = -cos(beamTheta)*sin(beamPhi);
   double xz = sin(beamTheta);
   
   double yx = sin(beamPhi);
   double yy = -cos(beamPhi);
   double yz = 0.;
   
   double zx = sin(beamTheta)*cos(beamPhi);
   double zy = sin(beamTheta)*sin(beamPhi);
   double zz = cos(beamTheta);
 
   fromCMStoTB = new ROOT::Math::Rotation3D(xx, xy, xz, yx, yy, yz, zx, zy, zz);
 
   // random number
   edm::Service<edm::RandomNumberGenerator> rng;
    if ( ! rng.isAvailable()) {
      throw cms::Exception("Configuration")
        << "EcalTBMCInfoProducer requires the RandomNumberGeneratorService\n"
           "which is not present in the configuration file.  You must add the service\n"
           "in the configuration file or remove the modules that require it.";
    }
    CLHEP::HepRandomEngine& engine = rng->getEngine();
    flatDistribution_ = new CLHEP::RandFlat(engine);
 
 }

EcalTBMCInfoProducer::~EcalTBMCInfoProducer ( )

virtual

Destructor.

Definition at line 103 of file EcalTBMCInfoProducer.cc.

References flatDistribution_, and theTestMap.

                                             {
 
   delete flatDistribution_;
   delete theTestMap;
   
 }

Member Function Documentation

void EcalTBMCInfoProducer::produce	(	edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)

virtual

Produce digis out of raw data.

Implements edm::EDProducer.

Definition at line 110 of file EcalTBMCInfoProducer.cc.

References beamEta, beamPhi, beamXoff, beamYoff, crysNumber, flatDistribution_, GenVtxLabel, LogDebug, partXhodo, and partYhodo.

 {
   auto_ptr<PEcalTBInfo> product(new PEcalTBInfo());
 
   // Fill the run information
 
   product->setCrystal(crysNumber);
 
   product->setBeamDirection(beamEta, beamPhi);
   product->setBeamOffset(beamXoff, beamYoff);
 
   // Compute the event x,y vertex coordinates in the beam reference system
   // e.g. in the place orthogonal to the beam average direction
 
   partXhodo = partYhodo = 0.;
 
   edm::Handle<edm::HepMCProduct> GenEvt;
   event.getByLabel(GenVtxLabel,GenEvt);
 
   const HepMC::GenEvent* Evt = GenEvt->GetEvent() ;
   HepMC::GenEvent::vertex_const_iterator Vtx = Evt->vertices_begin();
 
   math::XYZPoint eventCMSVertex((*Vtx)->position().x(),
                                 (*Vtx)->position().y(),
                                 (*Vtx)->position().z());
   
   LogDebug("EcalTBInfo") << "Generated vertex position = " 
                          << eventCMSVertex.x() << " " 
                          << eventCMSVertex.y() << " " 
                          << eventCMSVertex.z();
 
   math::XYZPoint eventTBVertex = (*fromCMStoTB)*eventCMSVertex;
 
   LogDebug("EcalTBInfo") << "Rotated vertex position   = "
                          << eventTBVertex.x() << " " 
                          << eventTBVertex.y() << " " 
                          << eventTBVertex.z();
 
   partXhodo = eventTBVertex.x();
   partYhodo = eventTBVertex.y();
 
   product->setBeamPosition(partXhodo, partYhodo);
 
   // Asynchronous phase shift
   double thisPhaseShift = flatDistribution_->fire();
 
   product->setPhaseShift(thisPhaseShift);
   LogDebug("EcalTBInfo") << "Asynchronous Phaseshift = " << thisPhaseShift;
 
   // store the object in the framework event
 
   event.put(product);
 }