PFGeometry.cc

Go to the documentation of this file.

#include "RecoParticleFlow/PFTracking/interface/PFGeometry.h"
#include "DataFormats/GeometrySurface/interface/BoundCylinder.h"
#include "DataFormats/GeometrySurface/interface/BoundDisk.h"
#include "DataFormats/GeometrySurface/interface/SimpleCylinderBounds.h"
#include "DataFormats/GeometrySurface/interface/SimpleDiskBounds.h"

PFGeometry::PFGeometry() {
  if (innerRadius_.empty()) {
    // All distances are in cm
    // BeamPipe
    innerRadius_.push_back(2.5);
    outerRadius_.push_back(2.5);
    innerZ_.push_back(0.);
    outerZ_.push_back(500.);
    // PS1
    innerRadius_.push_back(45.0);
    outerRadius_.push_back(125.0);
    innerZ_.push_back(303.16);
    outerZ_.push_back(303.16);
    // PS2
    innerRadius_.push_back(45.0);
    outerRadius_.push_back(125.0);
    innerZ_.push_back(307.13);
    outerZ_.push_back(307.13);
    // ECALBarrel
    innerRadius_.push_back(129.0);
    outerRadius_.push_back(175.0);
    innerZ_.push_back(0.);
    outerZ_.push_back(304.5);
    // ECALEndcap
    innerRadius_.push_back(31.6);
    outerRadius_.push_back(171.1);
    innerZ_.push_back(317.0);
    outerZ_.push_back(388.0);
    // HCALBarrel
    innerRadius_.push_back(183.0);
    outerRadius_.push_back(285.0);
    innerZ_.push_back(0.);
    outerZ_.push_back(433.2);
    // HCALEndcap
    innerRadius_.push_back(31.6);   // !!! Do not use : Probably wrong !!!
    outerRadius_.push_back(285.0);  // !!! Do not use : Probably wrong !!!
    innerZ_.push_back(388.0);
    outerZ_.push_back(560.0);

    // HO Barrel
    innerRadius_.push_back(387.6);
    outerRadius_.push_back(410.2);
    innerZ_.push_back(0.);
    outerZ_.push_back(700.25);

    // Define reference surfaces
    tanTh_.push_back(innerRadius_[BeamPipe] / outerZ_[BeamPipe]);
    tanTh_.push_back(outerRadius_[PS1] / outerZ_[PS1]);
    tanTh_.push_back(outerRadius_[PS2] / outerZ_[PS2]);
    tanTh_.push_back(innerRadius_[ECALBarrel] / innerZ_[ECALEndcap]);
    tanTh_.push_back(innerRadius_[HCALBarrel] / innerZ_[HCALEndcap]);
    tanTh_.push_back(outerRadius_[HCALBarrel] / outerZ_[HCALEndcap]);
  }
}