DataFormats/EcalDetId/test/testEcalDetId.cpp - master - githubmirror/cmssw

DataFormats/EcalDetId/test/testEcalDetId.cpp
334 строки12 KB
Cms Build
Clang-Format
09 май 2019, 07:14
09 май 2019, 07:1457cea53
/**
   \file
   Test suit for EcalDetId


   \note This test is not exaustive     
*/

#include <Utilities/Testing/interface/CppUnit_testdriver.icpp>
#include <cppunit/extensions/HelperMacros.h>
#include <DataFormats/EcalDetId/interface/EBDetId.h>
#include <DataFormats/EcalDetId/interface/EEDetId.h>
#include <DataFormats/EcalDetId/interface/ESDetId.h>
#include <DataFormats/EcalDetId/interface/EcalTrigTowerDetId.h>
#include <DataFormats/EcalDetId/interface/EcalPnDiodeDetId.h>
#include "FWCore/Utilities/interface/Exception.h"

#include <vector>
#include <algorithm>
#include <iostream>

class testEcalDetId : public CppUnit::TestFixture {
  CPPUNIT_TEST_SUITE(testEcalDetId);

  CPPUNIT_TEST(testEBDetId);
  CPPUNIT_TEST(testEEDetId);
  CPPUNIT_TEST(testESDetId);
  CPPUNIT_TEST(testEcalTrigTowerDetId);
  // CPPUNIT_TEST(testEcalElectronicsId);
  CPPUNIT_TEST(testPnDiodeDetId);
  CPPUNIT_TEST(testDistancePhi);

  CPPUNIT_TEST_SUITE_END();

public:
  void setUp() {}
  void tearDown() {}
  void testEBDetId();
  void testEEDetId();
  void testESDetId();
  void testEcalTrigTowerDetId();
  // void testEcalElectronicsId();
  void testPnDiodeDetId();
  void testDistancePhi();
};

///registration of the test so that the runner can find it
CPPUNIT_TEST_SUITE_REGISTRATION(testEcalDetId);

void testEcalDetId::testEBDetId() {
  EBDetId smId;

  std::vector<unsigned int> detIds(EBDetId::kSizeForDenseIndexing, 0);

  CPPUNIT_ASSERT(!EBDetId::validDetId(0, 1));
  CPPUNIT_ASSERT(!EBDetId::validDetId(1, 0));
  CPPUNIT_ASSERT(!EBDetId::validDetId(1, -1));
  CPPUNIT_ASSERT(!EBDetId::validDetId(1, EBDetId::MAX_IPHI + 1));
  CPPUNIT_ASSERT(!EBDetId::validDetId(EBDetId::MAX_IETA + 1, 1));
  CPPUNIT_ASSERT(!EBDetId::validDetId(-EBDetId::MAX_IETA - 1, 1));

  for (int ieta = EBDetId::MIN_IETA; ieta <= EBDetId::MAX_IETA; ieta++)
    for (int iphi = EBDetId::MIN_IPHI; iphi <= EBDetId::MAX_IPHI; iphi++) {
      //EBDetId Zside 1
      try {
        {
          EBDetId aPositiveId(ieta, iphi);
          CPPUNIT_ASSERT(EBDetId::validDetId(ieta, iphi));
          CPPUNIT_ASSERT(aPositiveId.ieta() == ieta);
          CPPUNIT_ASSERT(aPositiveId.iphi() == iphi);
          CPPUNIT_ASSERT(aPositiveId.zside() == 1);
          CPPUNIT_ASSERT(aPositiveId.ietaAbs() == ieta);
          int i = 0;
          for (; i != 4; ++i)
            if (aPositiveId.ietaAbs() <= EBDetId::kModuleBoundaries[i])
              break;
          CPPUNIT_ASSERT(aPositiveId.im() == i + 1);
          CPPUNIT_ASSERT(
              !(EBDetId::isNextToEtaBoundary(aPositiveId) ^
                (ieta == 1 || std::binary_search(EBDetId::kModuleBoundaries, EBDetId::kModuleBoundaries + 4, ieta))));

          smId = EBDetId(aPositiveId.ism(), aPositiveId.ic(), EBDetId::SMCRYSTALMODE);
          CPPUNIT_ASSERT(aPositiveId == smId);
          CPPUNIT_ASSERT(aPositiveId.ism() == smId.ism());
          CPPUNIT_ASSERT(aPositiveId.ic() == smId.ic());
          CPPUNIT_ASSERT(smId.ism() >= EBDetId::MIN_SM);
          CPPUNIT_ASSERT(smId.ism() <= EBDetId::MAX_SM);
          CPPUNIT_ASSERT(smId.ic() >= EBDetId::MIN_C);
          CPPUNIT_ASSERT(smId.ic() <= EBDetId::MAX_C);
          CPPUNIT_ASSERT(EBDetId::validHashIndex(aPositiveId.hashedIndex()));
          CPPUNIT_ASSERT(EBDetId::unhashIndex(aPositiveId.hashedIndex()) == aPositiveId);
          detIds.at(aPositiveId.hashedIndex()) = aPositiveId;
        }

        //EBDetId Zside -1
        {
          EBDetId aNegativeId(-1 * ieta, iphi);
          CPPUNIT_ASSERT(EBDetId::validDetId(-ieta, iphi));
          CPPUNIT_ASSERT(aNegativeId.ieta() == -1 * ieta);
          CPPUNIT_ASSERT(aNegativeId.iphi() == iphi);
          CPPUNIT_ASSERT(aNegativeId.zside() == -1);
          CPPUNIT_ASSERT(aNegativeId.ietaAbs() == ieta);

          smId = EBDetId(aNegativeId.ism(), aNegativeId.ic(), EBDetId::SMCRYSTALMODE);
          CPPUNIT_ASSERT(aNegativeId == smId);
          CPPUNIT_ASSERT(aNegativeId.ism() == smId.ism());
          CPPUNIT_ASSERT(aNegativeId.ic() == smId.ic());
          CPPUNIT_ASSERT(smId.ism() >= EBDetId::MIN_SM);
          CPPUNIT_ASSERT(smId.ism() <= EBDetId::MAX_SM);
          CPPUNIT_ASSERT(smId.ic() >= EBDetId::MIN_C);
          CPPUNIT_ASSERT(smId.ic() <= EBDetId::MAX_C);
          CPPUNIT_ASSERT(EBDetId::validHashIndex(aNegativeId.hashedIndex()));
          CPPUNIT_ASSERT(EBDetId::unhashIndex(aNegativeId.hashedIndex()) == aNegativeId);
          detIds.at(aNegativeId.hashedIndex()) = aNegativeId;
        }
      } catch (cms::Exception &e) {
        std::cout << e.what() << " failed for " << ieta << "," << iphi << std::endl;
        bool cmsExceptionCought = false;
        CPPUNIT_ASSERT(cmsExceptionCought);
      } catch (std::exception &e) {
        std::cout << e.what() << " failed for " << ieta << "," << iphi << std::endl;
        bool stdExceptionCought = false;
        CPPUNIT_ASSERT(stdExceptionCought);
      }
    }

  for (int i = 0; i != EBDetId::kSizeForDenseIndexing; ++i) {
    CPPUNIT_ASSERT(detIds[i] != 0);
    CPPUNIT_ASSERT(EBDetId(detIds[i]).hashedIndex() == i);
    CPPUNIT_ASSERT(EBDetId::unhashIndex(i) == detIds[i]);
  }
}

void testEcalDetId::testEEDetId() {
  std::vector<unsigned int> detIds(EEDetId::kSizeForDenseIndexing, 0);

  for (int ix = EEDetId::IX_MIN; ix <= EEDetId::IX_MAX; ix++)
    for (int iy = EEDetId::IY_MIN; iy <= EEDetId::IY_MAX; iy++)
      try {
        bool fastV = EEDetId::fastValidDetId(ix, iy);
        //EEDetId Zside 1
        if (EEDetId::slowValidDetId(ix, iy)) {
          CPPUNIT_ASSERT(EEDetId::validDetId(ix, iy, 1));
          EEDetId aPositiveId(ix, iy, 1);
          CPPUNIT_ASSERT(aPositiveId.ix() == ix);
          CPPUNIT_ASSERT(aPositiveId.iy() == iy);
          CPPUNIT_ASSERT(aPositiveId.zside() == 1);
          CPPUNIT_ASSERT(EEDetId::validHashIndex(aPositiveId.hashedIndex()));
          CPPUNIT_ASSERT(EEDetId::unhashIndex(aPositiveId.hashedIndex()) == aPositiveId);
          detIds.at(aPositiveId.hashedIndex()) = aPositiveId;
        } else {
          CPPUNIT_ASSERT(!fastV);
        }
        //EEDetId Zside -1
        if (EEDetId::slowValidDetId(ix, iy)) {
          CPPUNIT_ASSERT(EEDetId::validDetId(ix, iy, -1));
          EEDetId aNegativeId(ix, iy, -1);
          CPPUNIT_ASSERT(aNegativeId.ix() == ix);
          CPPUNIT_ASSERT(aNegativeId.iy() == iy);
          CPPUNIT_ASSERT(aNegativeId.zside() == -1);
          CPPUNIT_ASSERT(EEDetId::validHashIndex(aNegativeId.hashedIndex()));
          CPPUNIT_ASSERT(EEDetId::unhashIndex(aNegativeId.hashedIndex()) == aNegativeId);
          detIds.at(aNegativeId.hashedIndex()) = aNegativeId;
        }
      } catch (cms::Exception &e) {
        bool cmsExceptionCought = false;
        CPPUNIT_ASSERT(cmsExceptionCought);
      } catch (std::exception &e) {
        bool stdExceptionCought = false;
        CPPUNIT_ASSERT(stdExceptionCought);
      }

  int holes = 0;
  for (int i = 0; i != EEDetId::kSizeForDenseIndexing; ++i) {
    // CPPUNIT_ASSERT(detIds[i]!=0);
    if (detIds[i] == 0) {
      holes++;
      continue;
    }  // there are holes...
    CPPUNIT_ASSERT(EEDetId(detIds[i]).hashedIndex() == i);
    CPPUNIT_ASSERT(EEDetId::unhashIndex(i) == detIds[i]);
  }
  CPPUNIT_ASSERT(holes == 0);
  //FIXME hope a better test...
  //CPPUNIT_ASSERT(holes>EEDetId::kSizeForDenseIndexing/100);
  //CPPUNIT_ASSERT(holes<EEDetId::kSizeForDenseIndexing/10);
}

void testEcalDetId::testESDetId() {
  for (int istrip = ESDetId::ISTRIP_MIN; istrip <= ESDetId::ISTRIP_MAX; istrip++)
    for (int ix = ESDetId::IX_MIN; ix <= ESDetId::IX_MAX; ix++)
      for (int iy = ESDetId::IY_MIN; iy <= ESDetId::IY_MAX; iy++) {
        try {
          //ESDetId Plane 0 Zside 1
          {
            ESDetId aPositiveId(istrip, ix, iy, 1, 1);
            CPPUNIT_ASSERT(aPositiveId.strip() == istrip);
            CPPUNIT_ASSERT(aPositiveId.six() == ix);
            CPPUNIT_ASSERT(aPositiveId.siy() == iy);
            CPPUNIT_ASSERT(aPositiveId.plane() == 1);
            CPPUNIT_ASSERT(aPositiveId.zside() == 1);
          }
          //ESDetId Plane 1 Zside 1
          {
            ESDetId aPositiveId(istrip, ix, iy, 2, 1);
            CPPUNIT_ASSERT(aPositiveId.strip() == istrip);
            CPPUNIT_ASSERT(aPositiveId.six() == ix);
            CPPUNIT_ASSERT(aPositiveId.siy() == iy);
            CPPUNIT_ASSERT(aPositiveId.plane() == 2);
            CPPUNIT_ASSERT(aPositiveId.zside() == 1);
          }
          //ESDetId Plane 0 Zside -1
          {
            ESDetId aNegativeId(istrip, ix, iy, 1, -1);
            CPPUNIT_ASSERT(aNegativeId.strip() == istrip);
            CPPUNIT_ASSERT(aNegativeId.six() == ix);
            CPPUNIT_ASSERT(aNegativeId.siy() == iy);
            CPPUNIT_ASSERT(aNegativeId.plane() == 1);
            CPPUNIT_ASSERT(aNegativeId.zside() == -1);
          }
          //ESDetId Plane 1 Zside -1
          {
            ESDetId aNegativeId(istrip, ix, iy, 2, -1);
            CPPUNIT_ASSERT(aNegativeId.strip() == istrip);
            CPPUNIT_ASSERT(aNegativeId.six() == ix);
            CPPUNIT_ASSERT(aNegativeId.siy() == iy);
            CPPUNIT_ASSERT(aNegativeId.plane() == 2);
            CPPUNIT_ASSERT(aNegativeId.zside() == -1);
          }
        } catch (cms::Exception &e) {
        }
      }
}

void testEcalDetId::testEcalTrigTowerDetId() {
  for (int ieta = EcalTrigTowerDetId::MIN_I; ieta <= EcalTrigTowerDetId::MAX_I; ieta++)
    for (int iphi = EcalTrigTowerDetId::MIN_I; iphi <= EcalTrigTowerDetId::MAX_I; iphi++) {
      try {
        //EcalTrigTowerDetId Zside 1
        {
          EcalTrigTowerDetId aPositiveId(1, EcalBarrel, ieta, iphi);
          CPPUNIT_ASSERT(aPositiveId.subDet() == EcalBarrel);
          CPPUNIT_ASSERT(aPositiveId.ieta() == ieta);
          CPPUNIT_ASSERT(aPositiveId.iphi() == iphi);
          CPPUNIT_ASSERT(aPositiveId.zside() == 1);
          CPPUNIT_ASSERT(aPositiveId.ietaAbs() == ieta);
        }
        //EcalTrigTowerDetId Zside -1
        {
          EcalTrigTowerDetId aNegativeId(-1, EcalBarrel, ieta, iphi);
          CPPUNIT_ASSERT(aNegativeId.subDet() == EcalBarrel);
          CPPUNIT_ASSERT(aNegativeId.ieta() == -1 * ieta);
          CPPUNIT_ASSERT(aNegativeId.iphi() == iphi);
          CPPUNIT_ASSERT(aNegativeId.zside() == -1);
          CPPUNIT_ASSERT(aNegativeId.ietaAbs() == ieta);
        }
      } catch (cms::Exception &e) {
      }
      try {
        //EcalTrigTowerDetId Zside 1
        {
          EcalTrigTowerDetId aPositiveId(1, EcalEndcap, ieta, iphi);
          CPPUNIT_ASSERT(aPositiveId.ieta() == ieta);
          CPPUNIT_ASSERT(aPositiveId.subDet() == EcalEndcap);
          CPPUNIT_ASSERT(aPositiveId.iphi() == iphi);
          CPPUNIT_ASSERT(aPositiveId.zside() == 1);
        }
        //EcalTrigTowerDetId Zside -1
        {
          EcalTrigTowerDetId aNegativeId(-1, EcalEndcap, ieta, iphi);
          CPPUNIT_ASSERT(aNegativeId.subDet() == EcalEndcap);
          CPPUNIT_ASSERT(aNegativeId.ieta() == -1 * ieta);
          CPPUNIT_ASSERT(aNegativeId.iphi() == iphi);
          CPPUNIT_ASSERT(aNegativeId.zside() == -1);
        }
      } catch (cms::Exception &e) {
      }
    }
}

/*
void testEcalDetId::testEcalElectronicsId() {
  for (int idcc=EcalElectronicsId::MIN_DCCID;idcc<=EcalElectronicsId::MAX_DCCID;idcc++)
    for (int itower=EcalElectronicsId::MIN_TOWERID;itower<EcalElectronicsId::MAX_TOWERID;itower++)
      for (int ichannel=EcalElectronicsId::MIN_CHANNELID;ichannel<EcalElectronicsId::MAX_CHANNELID;ichannel++)
	{
	  try
	    {
	      EcalElectronicsId electronicsId(idcc,itower,ichannel);
	      CPPUNIT_ASSERT(electronicsId.dccId()==idcc);
	      CPPUNIT_ASSERT(electronicsId.towerId()==itower);
	      CPPUNIT_ASSERT(electronicsId.channelId()==ichannel);
	    }
	  catch ( cms::Exception &e ) 
	    { 
	    }
	}
}
*/

void testEcalDetId::testPnDiodeDetId() {
  for (int iEcalSubDet = EcalBarrel; iEcalSubDet <= EcalEndcap; iEcalSubDet++)
    for (int iDCC = 1; iDCC <= 36; iDCC++)
      for (int iPn = 1; iPn <= 10; iPn++) {
        try {
          //EcalTrigTowerDetId Zside 1
          {
            EcalPnDiodeDetId aId(iEcalSubDet, iDCC, iPn);
            CPPUNIT_ASSERT(aId.iEcalSubDetectorId() == iEcalSubDet);
            CPPUNIT_ASSERT(aId.iDCCId() == iDCC);
            CPPUNIT_ASSERT(aId.iPnId() == iPn);
          }
        } catch (cms::Exception &e) {
        }
      }
}

void testEcalDetId::testDistancePhi() {
  EBDetId a1(1, 360);
  EBDetId b1(1, 1);
  CPPUNIT_ASSERT(EBDetId::distancePhi(a1, b1) == 1);

  EBDetId a2(1, 1);
  EBDetId b2(1, 360);
  CPPUNIT_ASSERT(EBDetId::distancePhi(a2, b2) == 1);

  EBDetId a3(1, 175);
  EBDetId b3(1, 185);
  CPPUNIT_ASSERT(EBDetId::distancePhi(a3, b3) == 10);

  EBDetId a4(1, 350);
  EBDetId b4(1, 3);
  CPPUNIT_ASSERT(EBDetId::distancePhi(a4, b4) == 13);
}