CMS 3D CMS Logo

EcalTBMCInfoProducer Class Reference

#include <SimG4CMS/EcalTestBeam/interface/EcalTBMCInfoProducer.h>

Inheritance diagram for EcalTBMCInfoProducer:

edm::EDProducer edm::ProducerBase edm::ProductRegistryHelper

List of all members.

Public Member Functions

 EcalTBMCInfoProducer (const edm::ParameterSet &ps)
 Constructor.
void produce (edm::Event &event, const edm::EventSetup &eventSetup)
 Produce digis out of raw data.
virtual ~EcalTBMCInfoProducer ()
 Destructor.

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

Detailed Description

Definition at line 35 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, deltaPhi(), eta, Exception, funct::exp(), EcalTBCrystalMap::findCrystalAngles(), flatDistribution_, fromCMStoTB, GenVtxLabel, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), edm::Service< T >::isAvailable(), EcalTBCrystalMap::NCRYSTAL, phi, funct::sin(), and theTestMap.

00019                                                                     : flatDistribution_(0) {
00020   
00021   produces<PEcalTBInfo>();
00022 
00023   edm::FileInPath CrystalMapFile = ps.getParameter<edm::FileInPath>("CrystalMapFile");
00024   GenVtxLabel = ps.getUntrackedParameter<string>("moduleLabelVtx","source");
00025   double fMinEta = ps.getUntrackedParameter<double>("MinEta");
00026   double fMaxEta = ps.getUntrackedParameter<double>("MaxEta");
00027   double fMinPhi = ps.getUntrackedParameter<double>("MinPhi");
00028   double fMaxPhi = ps.getUntrackedParameter<double>("MaxPhi");
00029   beamEta = (fMaxEta+fMinEta)/2.;
00030   beamPhi = (fMaxPhi+fMinPhi)/2.;
00031   beamTheta = 2.0*atan(exp(-beamEta));
00032   beamXoff = ps.getUntrackedParameter<double>("BeamMeanX",0.0);
00033   beamYoff = ps.getUntrackedParameter<double>("BeamMeanX",0.0);
00034    
00035   string fullMapName = CrystalMapFile.fullPath();
00036   theTestMap = new EcalTBCrystalMap(fullMapName);
00037   crysNumber = 0;
00038 
00039   double deltaEta = 999.;
00040   double deltaPhi = 999.;
00041   for ( int cryIndex = 1; cryIndex <= EcalTBCrystalMap::NCRYSTAL; ++cryIndex) {
00042     double eta = 0;
00043     double phi = 0.;
00044     theTestMap->findCrystalAngles(cryIndex, eta, phi);
00045     if ( fabs(beamEta - eta) < deltaEta && fabs(beamPhi - phi) < deltaPhi ) {
00046       deltaEta = fabs(beamEta - eta);
00047       deltaPhi = fabs(beamPhi - phi);
00048       crysNumber = cryIndex;
00049     }
00050     else if (fabs(beamEta - eta)<deltaEta && fabs(beamPhi - phi)>deltaPhi ) {
00051       if ( fabs(beamPhi - phi) < 0.017 ) {
00052         deltaEta = fabs(beamEta - eta);
00053         deltaPhi = fabs(beamPhi - phi);
00054         crysNumber = cryIndex;
00055       }
00056     }
00057     else if (fabs(beamEta - eta)>deltaEta && fabs(beamPhi - phi)<deltaPhi ) {
00058       if ( fabs(beamEta - eta) < 0.017 ) {
00059         deltaEta = fabs(beamEta - eta);
00060         deltaPhi = fabs(beamPhi - phi);
00061         crysNumber = cryIndex;
00062       }
00063     }
00064   }
00065 
00066   edm::LogInfo("EcalTBInfo") << "Initialize TB MC ECAL info producer with parameters: \n"
00067                              << "Crystal map file:  " << CrystalMapFile << "\n"
00068                              << "Beam average eta = " << beamEta << "\n"
00069                              << "Beam average phi = " << beamPhi << "\n"
00070                              << "Corresponding to crystal number = " << crysNumber << "\n"
00071                              << "Beam X offset =    " << beamXoff << "\n"
00072                              << "Beam Y offset =    " << beamYoff;
00073 
00074   // rotation matrix to move from the CMS reference frame to the test beam one
00075 
00076   double xx = -cos(beamTheta)*cos(beamPhi);
00077   double xy = -cos(beamTheta)*sin(beamPhi);
00078   double xz = sin(beamTheta);
00079   
00080   double yx = sin(beamPhi);
00081   double yy = -cos(beamPhi);
00082   double yz = 0.;
00083   
00084   double zx = sin(beamTheta)*cos(beamPhi);
00085   double zy = sin(beamTheta)*sin(beamPhi);
00086   double zz = cos(beamTheta);
00087 
00088   fromCMStoTB = new ROOT::Math::Rotation3D(xx, xy, xz, yx, yy, yz, zx, zy, zz);
00089 
00090   // random number
00091   edm::Service<edm::RandomNumberGenerator> rng;
00092    if ( ! rng.isAvailable()) {
00093      throw cms::Exception("Configuration")
00094        << "EcalTBMCInfoProducer requires the RandomNumberGeneratorService\n"
00095           "which is not present in the configuration file.  You must add the service\n"
00096           "in the configuration file or remove the modules that require it.";
00097    }
00098    CLHEP::HepRandomEngine& engine = rng->getEngine();
00099    flatDistribution_ = new CLHEP::RandFlat(engine);
00100 
00101 }

EcalTBMCInfoProducer::~EcalTBMCInfoProducer (  )  [virtual]

Destructor.

Definition at line 103 of file EcalTBMCInfoProducer.cc.

References flatDistribution_, and theTestMap.

00103                                             {
00104 
00105   delete flatDistribution_;
00106   delete theTestMap;
00107   
00108 }

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.

00111 {
00112   auto_ptr<PEcalTBInfo> product(new PEcalTBInfo());
00113 
00114   // Fill the run information
00115 
00116   product->setCrystal(crysNumber);
00117 
00118   product->setBeamDirection(beamEta, beamPhi);
00119   product->setBeamOffset(beamXoff, beamYoff);
00120 
00121   // Compute the event x,y vertex coordinates in the beam reference system
00122   // e.g. in the place orthogonal to the beam average direction
00123 
00124   partXhodo = partYhodo = 0.;
00125 
00126   edm::Handle<edm::HepMCProduct> GenEvt;
00127   event.getByLabel(GenVtxLabel,GenEvt);
00128 
00129   const HepMC::GenEvent* Evt = GenEvt->GetEvent() ;
00130   HepMC::GenEvent::vertex_const_iterator Vtx = Evt->vertices_begin();
00131 
00132   math::XYZPoint eventCMSVertex((*Vtx)->position().x(),
00133                                 (*Vtx)->position().y(),
00134                                 (*Vtx)->position().z());
00135   
00136   LogDebug("EcalTBInfo") << "Generated vertex position = " 
00137                          << eventCMSVertex.x() << " " 
00138                          << eventCMSVertex.y() << " " 
00139                          << eventCMSVertex.z();
00140 
00141   math::XYZPoint eventTBVertex = (*fromCMStoTB)*eventCMSVertex;
00142 
00143   LogDebug("EcalTBInfo") << "Rotated vertex position   = "
00144                          << eventTBVertex.x() << " " 
00145                          << eventTBVertex.y() << " " 
00146                          << eventTBVertex.z();
00147 
00148   partXhodo = eventTBVertex.x();
00149   partYhodo = eventTBVertex.y();
00150 
00151   product->setBeamPosition(partXhodo, partYhodo);
00152 
00153   // Asynchronous phase shift
00154   double thisPhaseShift = flatDistribution_->fire();
00155 
00156   product->setPhaseShift(thisPhaseShift);
00157   LogDebug("EcalTBInfo") << "Asynchronous Phaseshift = " << thisPhaseShift;
00158 
00159   // store the object in the framework event
00160 
00161   event.put(product);
00162 }

Member Data Documentation

double EcalTBMCInfoProducer::beamEta [private]

Definition at line 57 of file EcalTBMCInfoProducer.h.

Referenced by EcalTBMCInfoProducer(), and produce().

double EcalTBMCInfoProducer::beamPhi [private]

Definition at line 58 of file EcalTBMCInfoProducer.h.

Referenced by EcalTBMCInfoProducer(), and produce().

double EcalTBMCInfoProducer::beamTheta [private]

Definition at line 59 of file EcalTBMCInfoProducer.h.

Referenced by EcalTBMCInfoProducer().

double EcalTBMCInfoProducer::beamXoff [private]

Definition at line 63 of file EcalTBMCInfoProducer.h.

Referenced by EcalTBMCInfoProducer(), and produce().

double EcalTBMCInfoProducer::beamYoff [private]

Definition at line 64 of file EcalTBMCInfoProducer.h.

Referenced by EcalTBMCInfoProducer(), and produce().

int EcalTBMCInfoProducer::crysNumber [private]

Definition at line 61 of file EcalTBMCInfoProducer.h.

Referenced by EcalTBMCInfoProducer(), and produce().

CLHEP::RandFlat* EcalTBMCInfoProducer::flatDistribution_ [private]

Definition at line 75 of file EcalTBMCInfoProducer.h.

Referenced by EcalTBMCInfoProducer(), produce(), and ~EcalTBMCInfoProducer().

ROOT::Math::Rotation3D* EcalTBMCInfoProducer::fromCMStoTB [private]

Definition at line 71 of file EcalTBMCInfoProducer.h.

Referenced by EcalTBMCInfoProducer().

std::string EcalTBMCInfoProducer::GenVtxLabel [private]

Definition at line 73 of file EcalTBMCInfoProducer.h.

Referenced by EcalTBMCInfoProducer(), and produce().

double EcalTBMCInfoProducer::partXhodo [private]

Definition at line 66 of file EcalTBMCInfoProducer.h.

Referenced by produce().

double EcalTBMCInfoProducer::partYhodo [private]

Definition at line 67 of file EcalTBMCInfoProducer.h.

Referenced by produce().

EcalTBCrystalMap* EcalTBMCInfoProducer::theTestMap [private]

Definition at line 69 of file EcalTBMCInfoProducer.h.

Referenced by EcalTBMCInfoProducer(), and ~EcalTBMCInfoProducer().

The documentation for this class was generated from the following files:
Generated on Tue Jun 9 18:19:57 2009 for CMSSW by  doxygen 1.5.4