PhotonFixCMS Class Reference

#include <PhotonFixCMS.h>

Public Member Functions
double	fixedEnergy () const
const PhotonFix &	photonFix () const
	PhotonFixCMS (const reco::Photon &p)
double	sigmaEnergy () const

Static Public Member Functions
static bool	initialise (const edm::EventSetup &iSetup, const std::string &s="Nominal")

Private Attributes
PhotonFix	pf

Detailed Description

Definition at line 49 of file PhotonFixCMS.h.

Constructor & Destructor Documentation

PhotonFixCMS::PhotonFixCMS

(

const reco::Photon &

p

)

Definition at line 26 of file PhotonFixCMS.cc.

                                             :
  pf(p.energy(),p.superCluster()->eta(),p.superCluster()->phi(),p.r9()) {
}

Member Function Documentation

double PhotonFixCMS::fixedEnergy

(

)

const

Definition at line 152 of file PhotonFixCMS.cc.

References PhotonFix::fixedEnergy(), and pf.

                                       {
  
  return pf.fixedEnergy();
}

bool PhotonFixCMS::initialise ( const edm::EventSetup &  iSetup, const std::string &  s = "Nominal"  )

static

Definition at line 30 of file PhotonFixCMS.cc.

References PhotonFix::barrelCGap(), PhotonFix::barrelMGap(), PhotonFix::barrelSGap(), EnergyCorrector::c, gather_cfg::cout, DetId::Ecal, EcalBarrel, EcalEndcap, PhotonFix::endcapCGap(), PhotonFix::endcapCrystal(), PhotonFix::endcapMGap(), PhotonFix::endcapSGap(), eta, PV3DBase< T, PVType, FrameType >::eta(), geometry, edm::EventSetup::get(), PhotonFix::GetaPhi(), CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), CaloGeometry::getSubdetectorGeometry(), i, PhotonFix::initialised(), PhotonFix::initialiseParameters(), j, relval_steps::k, phi, PV3DBase< T, PVType, FrameType >::phi(), EEDetId::validDetId(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                             {
  
  if(PhotonFix::initialised()) return false;
 
  PhotonFix::initialiseParameters(s);
  
  // Get ECAL geometry
  edm::ESHandle<CaloGeometry> geoHandle;
  iSetup.get<CaloGeometryRecord>().get(geoHandle);
  const CaloGeometry& geometry = *geoHandle;
  
  // EB
  const CaloSubdetectorGeometry *barrelGeometry = geometry.getSubdetectorGeometry(DetId::Ecal, EcalBarrel);
  
  double bc[170][360][2];
  for(int iz(0);iz<2;iz++) {
    for(int ie(0);ie<85;ie++) {
      int id = ie+1;
      if (iz==0) id = ie-85; 
      for(int ip(0);ip<360;ip++) {
    EBDetId eb(id,ip+1);
    
    const CaloCellGeometry *cellGeometry = barrelGeometry->getGeometry(eb);
    GlobalPoint crystalPos = cellGeometry->getPosition();
    bc[85*iz+ie][ip][0]=crystalPos.eta();
    bc[85*iz+ie][ip][1]=crystalPos.phi();
      }
    }
  }
  
  for(unsigned i(0);i<169;i++) {
    for(unsigned j(0);j<360;j++) {
      unsigned k((j+1)%360);
     
      double eta = 0.25*(    bc[i][j][0]+bc[i+1][j][0]+
                     bc[i][k][0]+bc[i+1][k][0]);
      double phi = PhotonFix::GetaPhi(PhotonFix::GetaPhi(bc[i][j][1],bc[i+1][j][1]),
                PhotonFix::GetaPhi(bc[i][k][1],bc[i+1][k][1]));

      PhotonFix::barrelCGap(i,j,0,eta);
      PhotonFix::barrelCGap(i,j,1,phi);
      
      if((i%5)==4 && (j%2)==1) {
      PhotonFix::barrelSGap(i/5,j/2,0,eta);
    PhotonFix::barrelSGap(i/5,j/2,1,phi);   
      }
      
      if((j%20)==19) {
      if(i== 19) {PhotonFix::barrelMGap(0,j/20,0,eta); PhotonFix::barrelMGap(0,j/20,1,phi);}
      if(i== 39) {PhotonFix::barrelMGap(1,j/20,0,eta); PhotonFix::barrelMGap(1,j/20,1,phi);}
      if(i== 59) {PhotonFix::barrelMGap(2,j/20,0,eta); PhotonFix::barrelMGap(2,j/20,1,phi);}
      if(i== 84) {PhotonFix::barrelMGap(3,j/20,0,eta); PhotonFix::barrelMGap(3,j/20,1,phi);}
      if(i==109) {PhotonFix::barrelMGap(4,j/20,0,eta); PhotonFix::barrelMGap(4,j/20,1,phi);}
      if(i==129) {PhotonFix::barrelMGap(5,j/20,0,eta); PhotonFix::barrelMGap(5,j/20,1,phi);}
      if(i==149) {PhotonFix::barrelMGap(6,j/20,0,eta); PhotonFix::barrelMGap(6,j/20,1,phi);}
    
      }
    }
  }
  
  // EE
  const CaloSubdetectorGeometry *endcapGeometry = geometry.getSubdetectorGeometry(DetId::Ecal, EcalEndcap);
  
  double ec[2][100][100][2];
  bool valid[100][100];
  int val_count=0;
  for(int iz(0);iz<2;iz++) {
    for(int ix(0);ix<100;ix++) {
      for(int iy(0);iy<100;iy++) {
    valid[ix][iy] = EEDetId::validDetId(ix+1,iy+1,2*iz-1);
if(iz==0) PhotonFix::endcapCrystal(ix,iy,valid[ix][iy]);
    if(valid[ix][iy]) {
      EEDetId ee(ix+1,iy+1,2*iz-1);
    val_count+=1;
      
      const CaloCellGeometry *cellGeometry = endcapGeometry->getGeometry(ee);
      GlobalPoint crystalPos = cellGeometry->getPosition();
      ec[iz][ix][iy][0]=asinh(crystalPos.x()/fabs(crystalPos.z()));
      ec[iz][ix][iy][1]=asinh(crystalPos.y()/fabs(crystalPos.z()));
    }
      }
    }
  }
  std::cout << "GG valid " << val_count << std::endl;
  double c[2];
  for(unsigned iz(0);iz<2;iz++) {
    unsigned nC(0),nS(0);
    for(unsigned i(0);i<99;i++) {
      for(unsigned j(0);j<99;j++) {
    if(valid[i][j  ] && valid[i+1][j  ] && 
       valid[i][j+1] && valid[i+1][j+1]) {
      for(unsigned k(0);k<2;k++) {

      c[k] = 0.25*(ec[iz][i][j][k]+ec[iz][i+1][j][k]+ec[iz][i][j+1][k]+ec[iz][i+1][j+1][k]);

        PhotonFix::endcapCGap(iz,nC,k,c[k]);     
    }
      
      if((i%5)==4 && (j%5)==4) {
        for(unsigned k(0);k<2;k++) {
        PhotonFix::endcapSGap(iz,nS,k,c[k]);     
        }
        nS++;
      }
      nC++;
    }
      }
    }
    std::cout << "Endcap number of crystal, submodule boundaries = "
          << nC << ", " << nS << std::endl;
  }
  
  // Hardcode EE D-module gap to 0,0
    PhotonFix::endcapMGap(0,0,0,0.0);    
    PhotonFix::endcapMGap(0,0,1,0.0);    
    PhotonFix::endcapMGap(1,0,0,0.0);    
    PhotonFix::endcapMGap(1,0,1,0.0);    
  
  return true;
}

const PhotonFix & PhotonFixCMS::photonFix

(

)

const

Definition at line 162 of file PhotonFixCMS.cc.

References pf.

                                               {
  return pf;
}

double PhotonFixCMS::sigmaEnergy

(

)

const

Definition at line 157 of file PhotonFixCMS.cc.

References pf, and PhotonFix::sigmaEnergy().

                                       {
  
  return pf.sigmaEnergy();
}

Member Data Documentation

PhotonFix PhotonFixCMS::pf

private

Definition at line 63 of file PhotonFixCMS.h.

Referenced by fixedEnergy(), photonFix(), and sigmaEnergy().