EcalTestBeam/test/testEcalTBCrystalMap.cpp

0001 #include "Geometry/EcalTestBeam/interface/EcalTBCrystalMap.h"
0002 #include "FWCore/MessageLogger/interface/MessageLogger.h"
0003 #include "FWCore/ParameterSet/interface/FileInPath.h"
0004 #include "CLHEP/Random/RandFlat.h"
0005
0006 #include <string>
0007
0008 int main() {
0009   edm::FileInPath dataFile("Geometry/EcalTestBeam/data/BarrelSM1CrystalCenterElectron120GeV.dat");
0010   EcalTBCrystalMap theTestMap(dataFile.fullPath());
0011
0012   long nCrystal = 1700;
0013
0014   for (int i = 0; i < 10; ++i) {
0015     int thisCrystal = CLHEP::RandFlat::shootInt(nCrystal);
0016     double thisEta = 0.;
0017     double thisPhi = 0.;
0018
0019     theTestMap.findCrystalAngles(thisCrystal, thisEta, thisPhi);
0020
0021     edm::LogVerbatim("EcalTestBeam") << "Crystal number " << thisCrystal << " eta = " << thisEta
0022                                      << " phi = " << thisPhi;
0023
0024     int checkThisCrystal = theTestMap.CrystalIndex(thisEta, thisPhi);
0025
0026     edm::LogVerbatim("EcalTestBeam") << "(eta,phi) = " << thisEta << " , " << thisPhi
0027                                      << " corresponds to crystal n. = " << checkThisCrystal;
0028   }
0029
0030   return 0;
0031 }