cscdqm::Detector Class Reference

Detector geometry and addressing related imformation and routines. More...

#include <CSCDQM_Detector.h>

Public Member Functions
const bool	AddressFromString (const std::string &str_address, Address &adr) const
	Construct address from string. More...
const float	Area (const unsigned int station) const
	Calculate station area in eta/phi space. More...
const float	Area (const Address &adr) const
	Calculate address area in eta/phi space. More...
	Detector (const unsigned int p_partitions_x=0, const unsigned int p_partitions_y=0)
	Constructor. More...
unsigned int	GlobalChamberIndex (unsigned int side, unsigned int station, unsigned int ring, unsigned int chamber) const
	Return global chamber index on his geometric location. More...
const bool	NextAddress (unsigned int &i, const Address *&adr, const Address &mask) const
	Address iterator by mask. More...
const bool	NextAddressBox (unsigned int &i, const AddressBox *&box, const Address &mask) const
	Address box iterator by mask. More...
const bool	NextAddressBoxByPartition (unsigned int &i, const unsigned int px, const unsigned int py, AddressBox *&box)
	Address box iterator by partition. More...
const unsigned int	NumberOfChamberCFEBs (const unsigned int station, const unsigned int ring) const
	Returns the number of CFEBs per Chamber on given Station/Ring. More...
const unsigned int	NumberOfChamberHVs (const unsigned int station, const unsigned int ring) const
	Returns the number of HVs per Chamber on given Station/Ring. More...
const unsigned int	NumberOfChambers (const unsigned int station, const unsigned int ring) const
	Returns the number of chambers for the given station and ring. More...
const unsigned int	NumberOfRings (const unsigned int station) const
	Returns the number of rings for the given station. More...
void	PrintAddress (const Address &adr) const
	Prints address for debugging. More...

Private Member Functions
const float	Eta (const float r, const float z) const
const float	EtaToX (const float eta) const
	Transform eta coordinate to local canvas coordinate. More...
const float	PhiMaxCFEB (const int station, const int ring, const int chamber, const int cfeb) const
	Get Max phi boundary for particular CFEB. More...
const float	PhiMinCFEB (const int station, const int ring, const int chamber, const int cfeb) const
	Get Min phi boundary for particular CFEB. More...
const float	PhiToY (const float phi) const
	Transform phi coordinate to local canvas coordinate. More...
const float	RMaxHV (const int station, const int ring, const int n_hv) const
	Get R max parameter (used in address eta/phi calculation) More...
const float	RMinHV (const int station, const int ring, const int n_hv) const
	Get R min parameter (used in address eta/phi calculation) More...
const float	Z (const int station, const int ring) const
	Get Z parameter (used in address eta/phi calculation) More...

Private Attributes
AddressBox	boxes [9540]
PartitionMap	partitions
unsigned int	partitions_x
unsigned int	partitions_y
float	station_area [4]

Detailed Description

Detector geometry and addressing related imformation and routines.

Definition at line 169 of file CSCDQM_Detector.h.

Constructor & Destructor Documentation

Detector()

cscdqm::Detector::Detector	(	const unsigned int	p_partitions_x = 0,
		const unsigned int	p_partitions_y = 0
	)

Constructor.

Parameters

p_partition_x	Number of efficiency partitions on X axis
p_partition_y	Number of efficiency partitions on Y axis

Returns

Creating real eta/phi boxes for available addresses

Address box calculated successfully. Now lets cache its partition elements for performace.
Cached the most frequently used areas

Definition at line 33 of file CSCDQM_Detector.cc.

References cscdqm::AddressBox::adr, Area(), boxes, reco::ceil(), cscdqm::AddressMask::cfeb, cscdqm::Address::cfeb, cscdqm::AddressMask::chamber, cscdqm::Address::chamber, Eta(), egammaIdentification::eta_max, egammaIdentification::eta_min, EtaToX(), cscdqm::AddressMask::hv, cscdqm::Address::hv, mps_fire::i, createfilelist::int, cscdqm::AddressMask::layer, cscdqm::Address::mask, N_SIDES, N_STATIONS, NumberOfChamberCFEBs(), NumberOfChamberHVs(), NumberOfChambers(), NumberOfRings(), PARTITION_INDEX, PARTITION_STEP_X, PARTITION_STEP_Y, partitions, partitions_x, partitions_y, PhiMaxCFEB(), PhiMinCFEB(), PhiToY(), cscdqm::AddressMask::ring, cscdqm::Address::ring, RMaxHV(), RMinHV(), cscdqm::AddressMask::side, cscdqm::Address::side, Validation_hcalonly_cfi::sign, cscdqm::AddressMask::station, cscdqm::Address::station, station_area, findQualityFiles::v, x, cscdqm::AddressBox::xmax, cscdqm::AddressBox::xmin, testProducerWithPsetDescEmpty_cfi::y1, testProducerWithPsetDescEmpty_cfi::y2, cscdqm::AddressBox::ymax, cscdqm::AddressBox::ymin, and Z().

                                                                                         {
    partitions_x = p_partitions_x;
    partitions_y = p_partitions_y;

    unsigned int i = 0;
    Address adr;

    adr.mask.layer = false;
    adr.mask.side = adr.mask.station = adr.mask.ring = adr.mask.chamber = adr.mask.cfeb = adr.mask.hv = true;

    for (adr.side = 1; adr.side <= N_SIDES; adr.side++) {
      float sign = +1.0;
      if (adr.side == 2)
        sign = -1.0;
      for (adr.station = 1; adr.station <= N_STATIONS; adr.station++) {
        for (adr.ring = 1; adr.ring <= NumberOfRings(adr.station); adr.ring++) {
          for (adr.chamber = 1; adr.chamber <= NumberOfChambers(adr.station, adr.ring); adr.chamber++) {
            for (adr.cfeb = 1; adr.cfeb <= NumberOfChamberCFEBs(adr.station, adr.ring); adr.cfeb++) {
              for (adr.hv = 1; adr.hv <= NumberOfChamberHVs(adr.station, adr.ring); adr.hv++) {
                float z = Z(adr.station, adr.ring);
                float r_min = RMinHV(adr.station, adr.ring, adr.hv);
                float r_max = RMaxHV(adr.station, adr.ring, adr.hv);
                float eta_min = sign * Eta(r_min, z);
                float eta_max = sign * Eta(r_max, z);
                float x_min = EtaToX(eta_min);
                float x_max = EtaToX(eta_max);
                float phi_min = 0;
                float phi_max = 0;

                if (adr.station == 1 && adr.ring == 1 && adr.hv == 1) {
                  phi_min = PhiMinCFEB(adr.station, adr.ring, adr.chamber, 1);
                  phi_max = PhiMaxCFEB(adr.station, adr.ring, adr.chamber, NumberOfChamberCFEBs(adr.station, adr.ring));
                } else {
                  phi_min = PhiMinCFEB(adr.station, adr.ring, adr.chamber, adr.cfeb);
                  phi_max = PhiMaxCFEB(adr.station, adr.ring, adr.chamber, adr.cfeb);
                }

                float y_min = PhiToY(phi_min);
                float y_max = PhiToY(phi_max);

                boxes[i].adr = adr;

                float xboxmin = (x_min < x_max ? x_min : x_max);
                float xboxmax = (x_max > x_min ? x_max : x_min);
                float yboxmin = (y_min < y_max ? y_min : y_max);
                float yboxmax = (y_max > y_min ? y_max : y_min);

                boxes[i].xmin = xboxmin;
                boxes[i].xmax = xboxmax;
                boxes[i].ymin = yboxmin;
                boxes[i].ymax = yboxmax;

                unsigned int x1 = int(floor(xboxmin / PARTITION_STEP_X)) + int(partitions_x / 2);
                unsigned int x2 = int(ceil(xboxmax / PARTITION_STEP_X)) + int(partitions_x / 2);
                unsigned int y1 = int(floor(yboxmin / PARTITION_STEP_Y));
                unsigned int y2 = int(ceil(yboxmax / PARTITION_STEP_Y));

                for (unsigned int x = x1; x < x2; x++) {
                  for (unsigned int y = y1; y < y2; y++) {
                    unsigned int index = PARTITION_INDEX(x, y);
                    PartitionMapIterator iter = partitions.find(index);
                    if (iter == partitions.end()) {
                      std::vector<unsigned int> v;
                      partitions.insert(std::make_pair(index, v));
                    }
                    partitions[index].push_back(i);
                  }
                }

                i++;
              }
            }
          }
        }
      }
    }

    adr.mask.side = adr.mask.ring = adr.mask.chamber = adr.mask.layer = adr.mask.cfeb = adr.mask.hv = false;
    adr.mask.station = true;
    adr.station = 1;
    station_area[0] = Area(adr);
    adr.station = 2;
    station_area[1] = Area(adr);
    adr.station = 3;
    station_area[2] = Area(adr);
    adr.station = 4;
    station_area[3] = Area(adr);
  }

Member Function Documentation

AddressFromString()

const bool cscdqm::Detector::AddressFromString	(	const std::string &	str_address,
		Address &	adr
	)		const

Construct address from string.

Parameters

str_address	Address in string
adr	Address to return

Returns

true if address was successfully created, false - otherwise

Definition at line 658 of file CSCDQM_Detector.cc.

References ADDR_SIZE, cscdqm::AddressMask::cfeb, cscdqm::Address::cfeb, cscdqm::AddressMask::chamber, cscdqm::Address::chamber, TauDecayModes::dec, cscdqm::AddressMask::hv, cscdqm::Address::hv, cscdqm::AddressMask::layer, cscdqm::Address::layer, ALPAKA_ACCELERATOR_NAMESPACE::pixelClustering::pixelStatus::mask, cscdqm::Address::mask, EgammaValidation_cff::num, cscdqm::AddressMask::ring, cscdqm::Address::ring, cscdqm::AddressMask::side, cscdqm::Address::side, cscdqm::Utility::splitString(), cscdqm::AddressMask::station, cscdqm::Address::station, AlCaHLTBitMon_QueryRunRegistry::string, unpackBuffers-CaloStage2::token, and cscdqm::Utility::trimString().

Referenced by cscdqm::Summary::setMaskedHWElements().

                                                                                         {
    std::vector<std::string> tokens;
    Utility::splitString(str_address, ",", tokens);

    if (tokens.size() != ADDR_SIZE)
      return false;

    for (unsigned int r = 0; r < ADDR_SIZE; r++) {
      std::string token = tokens.at(r);
      Utility::trimString(token);
      bool mask = false;
      unsigned int num = 0;

      if (token != "*") {
        if (stringToNumber<unsigned int>(num, token, std::dec)) {
          mask = true;
        } else {
          return false;
        }
      }

      switch (r) {
        case 0:
          adr.mask.side = mask;
          adr.side = num;
          break;
        case 1:
          adr.mask.station = mask;
          adr.station = num;
          break;
        case 2:
          adr.mask.ring = mask;
          adr.ring = num;
          break;
        case 3:
          adr.mask.chamber = mask;
          adr.chamber = num;
          break;
        case 4:
          adr.mask.layer = mask;
          adr.layer = num;
          break;
        case 5:
          adr.mask.cfeb = mask;
          adr.cfeb = num;
          break;
        case 6:
          adr.mask.hv = mask;
          adr.hv = num;
      }
    }

    return true;
  }

Area() [1/2]

const float cscdqm::Detector::Area

(

const unsigned int

station

)

const

Calculate station area in eta/phi space.

Parameters

station

Station number

Returns

Area that is being covered by station

Definition at line 132 of file CSCDQM_Detector.cc.

References N_STATIONS, relativeConstraints::station, and station_area.

Referenced by Detector(), cscdqm::Summary::GetEfficiencyArea(), and cscdqm::Summary::GetReportingArea().

                                                             {
    if (station > 0 && station <= N_STATIONS) {
      return station_area[station - 1];
    }
    return 0;
  }

Area() [2/2]

const float cscdqm::Detector::Area

(

const Address &

adr

)

const

Calculate address area in eta/phi space.

Parameters

adr

Address

Returns

Area that is being covered by address

Definition at line 182 of file CSCDQM_Detector.cc.

References a, boxes, mps_fire::i, N_ELEMENTS, TrackerOfflineValidation_Dqm_cff::xmax, TrackerOfflineValidation_Dqm_cff::xmin, L1TOccupancyClient_cfi::ymax, and L1TOccupancyClient_cfi::ymin.

                                                     {
    float a = 0;
    for (unsigned int i = 0; i < N_ELEMENTS; i++) {
      if (boxes[i].adr == adr) {
        a += fabs((boxes[i].xmax - boxes[i].xmin) * (boxes[i].ymax - boxes[i].ymin));
      }
    }
    return a;
  }

Eta()

const float cscdqm::Detector::Eta ( const float  r, const float  z  ) const

private

Definition at line 391 of file CSCDQM_Detector.cc.

References dqm-mbProfile::log, and mathSSE::sqrt().

Referenced by Detector().

                                                              {
    if (r > 0.0 || z > 0.0) {
      float sin_theta = r / sqrt(r * r + z * z);
      float cos_theta = z / sqrt(r * r + z * z);
      return -log(sin_theta / (cos_theta + 1));
    }
    if (r == 0.0)
      return FLT_MAX;
    return 0.0;
  }

EtaToX()

const float cscdqm::Detector::EtaToX ( const float  eta ) const

private

Transform eta coordinate to local canvas coordinate.

Parameters

eta	Eta coordinate

Returns

local canvas coordinate

Definition at line 407 of file CSCDQM_Detector.cc.

References a, b, PVValHelper::eta, egammaIdentification::eta_max, and egammaIdentification::eta_min.

Referenced by Detector().

                                                    {
    float x_min = -2.5;
    float x_max = 2.5;
    float eta_min = -2.5;
    float eta_max = 2.5;
    float a = (x_max - x_min) / (eta_max - eta_min);
    float b = (eta_max * x_min - eta_min * x_max) / (eta_max - eta_min);
    return a * eta + b;
  }

GlobalChamberIndex()

unsigned int cscdqm::Detector::GlobalChamberIndex	(	unsigned int	side,
		unsigned int	station,
		unsigned int	ring,
		unsigned int	chamber
	)		const

Return global chamber index on his geometric location.

Parameters

side	Side (1,2)
station	Station
ring	Ring\
chamber	Chamber position

Returns

Global chamber index starting 1. If chamber is not existing - returns 0

Definition at line 147 of file CSCDQM_Detector.cc.

References cscdqm::AddressMask::cfeb, cscdqm::Address::cfeb, relativeConstraints::chamber, cscdqm::AddressMask::chamber, cscdqm::Address::chamber, cscdqm::AddressMask::hv, cscdqm::Address::hv, mps_fire::i, cscdqm::AddressMask::layer, cscdqm::Address::layer, cscdqm::Address::mask, N_SIDES, N_STATIONS, NumberOfChambers(), NumberOfRings(), relativeConstraints::ring, cscdqm::AddressMask::ring, cscdqm::Address::ring, cscdqm::AddressMask::side, cscdqm::Address::side, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::side(), relativeConstraints::station, cscdqm::AddressMask::station, and cscdqm::Address::station.

Referenced by cscdqm::EventProcessor::processCSC().

                                                                        {
    Address adr, iadr;
    adr.mask.side = adr.mask.station = adr.mask.ring = adr.mask.chamber = true;
    adr.mask.layer = adr.mask.cfeb = adr.mask.hv = false;
    adr.layer = adr.cfeb = adr.hv = 0;
    adr.side = side;
    adr.station = station;
    adr.ring = ring;
    adr.chamber = chamber;
    iadr = adr;

    unsigned int i = 1;
    for (iadr.side = 1; iadr.side <= N_SIDES; iadr.side++) {
      for (iadr.station = 1; iadr.station <= N_STATIONS; iadr.station++) {
        for (iadr.ring = 1; iadr.ring <= NumberOfRings(iadr.station); iadr.ring++) {
          for (iadr.chamber = 1; iadr.chamber <= NumberOfChambers(iadr.station, iadr.ring); iadr.chamber++) {
            if (iadr == adr) {
              return i;
            }
            i += 1;
          }
        }
      }
    }
    return 0;
  }

NextAddress()

const bool cscdqm::Detector::NextAddress	(	unsigned int &	i,
		const Address *&	adr,
		const Address &	mask
	)		const

Address iterator by mask.

Parameters

i	Iterator
adr	Address to return
mask	for addresses

Returns

true if address was found and filled in, false - otherwise

Definition at line 337 of file CSCDQM_Detector.cc.

References cscdqm::AddressBox::adr, boxes, mps_fire::i, ALPAKA_ACCELERATOR_NAMESPACE::pixelClustering::pixelStatus::mask, and N_ELEMENTS.

                                                                                                  {
    for (; i < N_ELEMENTS; i++) {
      if (boxes[i].adr == mask) {
        adr = &boxes[i].adr;
        i++;
        return true;
      }
    }
    return false;
  }

NextAddressBox()

const bool cscdqm::Detector::NextAddressBox	(	unsigned int &	i,
		const AddressBox *&	box,
		const Address &	mask
	)		const

Address box iterator by mask.

Parameters

i	Iterator
adr	AddressBox to return
mask	for addresses

Returns

true if address box was found and filled in, false - otherwise

Definition at line 355 of file CSCDQM_Detector.cc.

References boxes, mps_fire::i, ALPAKA_ACCELERATOR_NAMESPACE::pixelClustering::pixelStatus::mask, and N_ELEMENTS.

Referenced by cscdqm::Summary::Write().

                                                                                                        {
    for (; i < N_ELEMENTS; i++) {
      if (boxes[i].adr == mask) {
        box = &boxes[i];
        i++;
        return true;
      }
    }
    return false;
  }

NextAddressBoxByPartition()

const bool cscdqm::Detector::NextAddressBoxByPartition	(	unsigned int &	i,
		const unsigned int	px,
		const unsigned int	py,
		AddressBox *&	box
	)

Address box iterator by partition.

Parameters

i	Iterator
px	Partition x index
py	Partition y index
box	AddressBox to return

Returns

true if address box was found and filled in, false - otherwise

Definition at line 374 of file CSCDQM_Detector.cc.

References boxes, mps_fire::i, PARTITION_INDEX, partitions, multPhiCorr_741_25nsDY_cfi::px, and multPhiCorr_741_25nsDY_cfi::py.

Referenced by cscdqm::Summary::IsPhysicsReady().

                                                                   {
    unsigned int index = PARTITION_INDEX(px, py);

    PartitionMapIterator iter = partitions.find(index);
    if (iter != partitions.end()) {
      if (i < partitions[index].size()) {
        box = &boxes[partitions[index].at(i)];
        i++;
        return true;
      }
    }
    return false;
  }

NumberOfChamberCFEBs()

const unsigned int cscdqm::Detector::NumberOfChamberCFEBs	(	const unsigned int	station,
		const unsigned int	ring
	)		const

Returns the number of CFEBs per Chamber on given Station/Ring.

Parameters

station	Station number (1...4)
ring	Ring number (1...3)

Returns

Number of CFEBs per Chamber

Definition at line 243 of file CSCDQM_Detector.cc.

References relativeConstraints::ring, and relativeConstraints::station.

Referenced by Detector(), cscdqm::Summary::GetEfficiencyHW(), cscdqm::Summary::GetReportingArea(), cscdqm::Summary::GetValue(), PhiMaxCFEB(), PhiMinCFEB(), and cscdqm::Summary::SetValue().

                                                                                                             {
    if (station == 1 && ring == 1)
      return 4;
    if (station == 1 && ring == 2)
      return 5;
    if (station == 1 && ring == 3)
      return 4;
    if (station == 2 && ring == 1)
      return 5;
    if (station == 2 && ring == 2)
      return 5;
    if (station == 3 && ring == 1)
      return 5;
    if (station == 3 && ring == 2)
      return 5;
    if (station == 4 && ring == 1)
      return 5;
    if (station == 4 && ring == 2)
      return 5;
    return 0;
  }

NumberOfChamberHVs()

const unsigned int cscdqm::Detector::NumberOfChamberHVs	(	const unsigned int	station,
		const unsigned int	ring
	)		const

Returns the number of HVs per Chamber on given Station/Ring.

Parameters

station	Station number (1...4)
ring	Ring number (1...3)

Returns

Number of HVs per Chamber

Definition at line 271 of file CSCDQM_Detector.cc.

References relativeConstraints::ring, and relativeConstraints::station.

Referenced by Detector(), cscdqm::Summary::GetEfficiencyHW(), cscdqm::Summary::GetReportingArea(), cscdqm::Summary::GetValue(), and cscdqm::Summary::SetValue().

                                                                                                           {
    if (station == 1 && ring == 1)
      return 2;
    if (station == 1 && ring == 2)
      return 3;
    if (station == 1 && ring == 3)
      return 3;
    if (station == 2 && ring == 1)
      return 3;
    if (station == 2 && ring == 2)
      return 5;
    if (station == 3 && ring == 1)
      return 3;
    if (station == 3 && ring == 2)
      return 5;
    if (station == 4 && ring == 1)
      return 3;
    if (station == 4 && ring == 2)
      return 5;
    return 0;
  }

NumberOfChambers()

const unsigned int cscdqm::Detector::NumberOfChambers	(	const unsigned int	station,
		const unsigned int	ring
	)		const

Returns the number of chambers for the given station and ring.

Parameters

station	Station number (1...4)
ring	Ring number (1...3)

Returns

number of chambers

Definition at line 215 of file CSCDQM_Detector.cc.

References relativeConstraints::ring, and relativeConstraints::station.

Referenced by Detector(), cscdqm::Summary::GetEfficiencyHW(), cscdqm::Summary::GetReportingArea(), cscdqm::Summary::GetValue(), GlobalChamberIndex(), PhiMaxCFEB(), PhiMinCFEB(), cscdqm::Summary::SetValue(), cscdqm::EventProcessor::updateEfficiencyHistos(), and cscdqm::Summary::WriteChamberState().

                                                                                                         {
    if (station == 1 && ring == 1)
      return 36;
    if (station == 1 && ring == 2)
      return 36;
    if (station == 1 && ring == 3)
      return 36;
    if (station == 2 && ring == 1)
      return 18;
    if (station == 2 && ring == 2)
      return 36;
    if (station == 3 && ring == 1)
      return 18;
    if (station == 3 && ring == 2)
      return 36;
    if (station == 4 && ring == 1)
      return 18;
    if (station == 4 && ring == 2)
      return 36;
    return 0;
  }

NumberOfRings()

const unsigned int cscdqm::Detector::NumberOfRings

(

const unsigned int

station

)

const

Returns the number of rings for the given station.

Parameters

station

Station number (1, 2, 3, 4)

Returns

number of rings for the given station

Definition at line 197 of file CSCDQM_Detector.cc.

References relativeConstraints::station.

Referenced by Detector(), cscdqm::Summary::GetEfficiencyHW(), cscdqm::Summary::GetReportingArea(), cscdqm::Summary::GetValue(), GlobalChamberIndex(), cscdqm::Summary::SetValue(), cscdqm::EventProcessor::updateEfficiencyHistos(), and cscdqm::Summary::WriteChamberState().

                                                                             {
    if (station == 1)
      return 3;
    if (station == 2)
      return 2;
    if (station == 3)
      return 2;
    if (station == 4)
      return 2;
    return 0;
  }

PhiMaxCFEB()

const float cscdqm::Detector::PhiMaxCFEB ( const int  station, const int  ring, const int  chamber, const int  cfeb  ) const

private

Get Max phi boundary for particular CFEB.

Parameters

station	Station number
ring	Ring number
chamber	Chamber number
cfeb	CFEB number

Returns

Max phi CFEB boundary

Definition at line 641 of file CSCDQM_Detector.cc.

References relativeConstraints::chamber, nano_mu_digi_cff::float, NumberOfChamberCFEBs(), NumberOfChambers(), relativeConstraints::ring, and relativeConstraints::station.

Referenced by Detector().

                                                                                                             {
    float phi_max_cfeb;

    int n_cfeb = NumberOfChamberCFEBs(station, ring);
    int n_chambers = NumberOfChambers(station, ring);

    phi_max_cfeb = 0.0 + 2.0 * 3.14159 / (float)n_chambers * ((float)(chamber - 1) + (float)(cfeb) / (float)n_cfeb);

    return phi_max_cfeb;
  }

PhiMinCFEB()

const float cscdqm::Detector::PhiMinCFEB ( const int  station, const int  ring, const int  chamber, const int  cfeb  ) const

private

Get Min phi boundary for particular CFEB.

Parameters

station	Station number
ring	Ring number
chamber	Chamber number
cfeb	CFEB number

Returns

Min phi CFEB boundary

Definition at line 621 of file CSCDQM_Detector.cc.

References relativeConstraints::chamber, nano_mu_digi_cff::float, NumberOfChamberCFEBs(), NumberOfChambers(), relativeConstraints::ring, and relativeConstraints::station.

Referenced by Detector().

                                                                                                             {
    float phi_min_cfeb;

    int n_cfeb = NumberOfChamberCFEBs(station, ring);
    int n_chambers = NumberOfChambers(station, ring);

    phi_min_cfeb =
        0.0 + 2.0 * 3.14159 / ((float)(n_chambers)) * ((float)(chamber - 1) + (float)(cfeb - 1) / (float)(n_cfeb));

    return phi_min_cfeb;
  }

PhiToY()

const float cscdqm::Detector::PhiToY ( const float  phi ) const

private

Transform phi coordinate to local canvas coordinate.

Parameters

phi	Phi coordinate

Returns

local canvas coordinate

Definition at line 422 of file CSCDQM_Detector.cc.

References a, and b.

Referenced by Detector().

                                                    {
    float y_min = 0.0;
    float y_max = 2.0 * 3.14159;
    float phi_min = 0.0;
    float phi_max = 2.0 * 3.14159;
    float a = (y_max - y_min) / (phi_max - phi_min);
    float b = (phi_max * y_min - phi_min * y_max) / (phi_max - phi_min);
    return a * phi + b;
  }

PrintAddress()

void cscdqm::Detector::PrintAddress

(

const Address &

adr

)

const

Prints address for debugging.

Parameters

adr	Address to print

Returns

Definition at line 298 of file CSCDQM_Detector.cc.

References cscdqm::AddressMask::cfeb, cscdqm::Address::cfeb, cscdqm::AddressMask::chamber, cscdqm::Address::chamber, gather_cfg::cout, cscdqm::AddressMask::hv, cscdqm::Address::hv, cscdqm::AddressMask::layer, cscdqm::Address::layer, cscdqm::Address::mask, cscdqm::AddressMask::ring, cscdqm::Address::ring, cscdqm::AddressMask::side, cscdqm::Address::side, cscdqm::AddressMask::station, and cscdqm::Address::station.

                                                      {
    std::cout << "Side (" << std::boolalpha << adr.mask.side << ")";
    if (adr.mask.side)
      std::cout << " = " << adr.side;

    std::cout << ", Station (" << std::boolalpha << adr.mask.station << ")";
    if (adr.mask.station)
      std::cout << " = " << adr.station;

    std::cout << ", Ring (" << std::boolalpha << adr.mask.ring << ")";
    if (adr.mask.ring)
      std::cout << " = " << adr.ring;

    std::cout << ", Chamber (" << std::boolalpha << adr.mask.chamber << ")";
    if (adr.mask.chamber)
      std::cout << " = " << adr.chamber;

    std::cout << ", Layer (" << std::boolalpha << adr.mask.layer << ")";
    if (adr.mask.layer)
      std::cout << " = " << adr.layer;

    std::cout << ", CFEB (" << std::boolalpha << adr.mask.cfeb << ")";
    if (adr.mask.cfeb)
      std::cout << " = " << adr.cfeb;

    std::cout << ", HV (" << std::boolalpha << adr.mask.hv << ")";
    if (adr.mask.hv)
      std::cout << " = " << adr.hv;

    std::cout << std::endl;
  }

RMaxHV()

const float cscdqm::Detector::RMaxHV ( const int  station, const int  ring, const int  n_hv  ) const

private

Get R max parameter (used in address eta/phi calculation)

Parameters

station	Station Id
ring	Ring Id
n_hv	HV number

Returns

R max value

Definition at line 542 of file CSCDQM_Detector.cc.

References relativeConstraints::ring, and relativeConstraints::station.

Referenced by Detector().

                                                                                      {
    float r_max_hv = 0;

    if (station == 1 && ring == 1) {
      if (n_hv == 1)
        r_max_hv = 1500.0;
      if (n_hv == 2)
        r_max_hv = 2565.0;
    }

    if (station == 1 && ring == 2) {
      if (n_hv == 1)
        r_max_hv = 3368.2;
      if (n_hv == 2)
        r_max_hv = 4025.7;
      if (n_hv == 3)
        r_max_hv = 4559.9;
    }

    if (station == 1 && ring == 3) {
      if (n_hv == 1)
        r_max_hv = 5724.1;
      if (n_hv == 2)
        r_max_hv = 6230.2;
      if (n_hv == 3)
        r_max_hv = 6761.5;
    }

    if (station == 2 && ring == 1) {
      if (n_hv == 1)
        r_max_hv = 2152.3;
      if (n_hv == 2)
        r_max_hv = 2763.7;
      if (n_hv == 3)
        r_max_hv = 3365.8;
    }

    if (station == 3 && ring == 1) {
      if (n_hv == 1)
        r_max_hv = 2164.9;
      if (n_hv == 2)
        r_max_hv = 2763.8;
      if (n_hv == 3)
        r_max_hv = 3365.8;
    }

    if (station == 4 && ring == 1) {
      if (n_hv == 1)
        r_max_hv = 2365.9;
      if (n_hv == 2)
        r_max_hv = 2865.0;
      if (n_hv == 3)
        r_max_hv = 3356.3;
    }

    if ((station == 2 || station == 3 || station == 4) && ring == 2) {
      if (n_hv == 1)
        r_max_hv = 4446.3;
      if (n_hv == 2)
        r_max_hv = 5053.2;
      if (n_hv == 3)
        r_max_hv = 5660.1;
      if (n_hv == 4)
        r_max_hv = 6267.0;
      if (n_hv == 5)
        r_max_hv = 6870.8;
    }

    return r_max_hv;
  }

RMinHV()

const float cscdqm::Detector::RMinHV ( const int  station, const int  ring, const int  n_hv  ) const

private

Get R min parameter (used in address eta/phi calculation)

Parameters

station	Station Id
ring	Ring Id
n_hv	HV number

Returns

R min value

Definition at line 464 of file CSCDQM_Detector.cc.

References relativeConstraints::ring, and relativeConstraints::station.

Referenced by Detector().

                                                                                      {
    float r_min_hv = 0;

    if (station == 1 && ring == 1) {
      if (n_hv == 1)
        r_min_hv = 1060.0;
      if (n_hv == 2)
        r_min_hv = 1500.0;
    }

    if (station == 1 && ring == 2) {
      if (n_hv == 1)
        r_min_hv = 2815.0;
      if (n_hv == 2)
        r_min_hv = 3368.2;
      if (n_hv == 3)
        r_min_hv = 4025.7;
    }

    if (station == 1 && ring == 3) {
      if (n_hv == 1)
        r_min_hv = 5120.0;
      if (n_hv == 2)
        r_min_hv = 5724.1;
      if (n_hv == 3)
        r_min_hv = 6230.2;
    }

    if (station == 2 && ring == 1) {
      if (n_hv == 1)
        r_min_hv = 1469.2;
      if (n_hv == 2)
        r_min_hv = 2152.3;
      if (n_hv == 3)
        r_min_hv = 2763.7;
    }

    if (station == 3 && ring == 1) {
      if (n_hv == 1)
        r_min_hv = 1668.9;
      if (n_hv == 2)
        r_min_hv = 2164.9;
      if (n_hv == 3)
        r_min_hv = 2763.8;
    }

    if (station == 4 && ring == 1) {
      if (n_hv == 1)
        r_min_hv = 1876.1;
      if (n_hv == 2)
        r_min_hv = 2365.9;
      if (n_hv == 3)
        r_min_hv = 2865.0;
    }

    if ((station == 2 || station == 3 || station == 4) && ring == 2) {
      if (n_hv == 1)
        r_min_hv = 3640.2;
      if (n_hv == 2)
        r_min_hv = 4446.3;
      if (n_hv == 3)
        r_min_hv = 5053.2;
      if (n_hv == 4)
        r_min_hv = 5660.1;
      if (n_hv == 5)
        r_min_hv = 6267.0;
    }

    return r_min_hv;
  }

Z()

const float cscdqm::Detector::Z ( const int  station, const int  ring  ) const

private

Get Z parameter (used in address eta/phi calculation)

Parameters

station	Station Id
ring	Ring Id

Returns

Z value

Definition at line 438 of file CSCDQM_Detector.cc.

References relativeConstraints::ring, and relativeConstraints::station.

Referenced by Detector().

                                                                 {
    float z_csc = 0;

    if (station == 1 && ring == 1)
      z_csc = (5834.5 + 6101.5) / 2.0;
    if (station == 1 && ring == 2)
      z_csc = (6790.0 + 7064.3) / 2.0;
    if (station == 1 && ring == 3)
      z_csc = 6888.0;
    if (station == 2)
      z_csc = (8098.0 + 8346.0) / 2.0;
    if (station == 3)
      z_csc = (9414.8 + 9166.8) / 2.0;
    if (station == 4)
      z_csc = 10630.0;  // has to be corrected

    return z_csc;
  }

Member Data Documentation

boxes

AddressBox cscdqm::Detector::boxes[9540]

private

Address boxes in epa/phi space

Definition at line 207 of file CSCDQM_Detector.h.

Referenced by Area(), Detector(), NextAddress(), NextAddressBox(), and NextAddressBoxByPartition().

partitions

PartitionMap cscdqm::Detector::partitions

private

Map of partitions and list of it covering addresses indexes

Definition at line 219 of file CSCDQM_Detector.h.

Referenced by Detector(), and NextAddressBoxByPartition().

partitions_x

unsigned int cscdqm::Detector::partitions_x

private

Number of partitions in X axis

Definition at line 213 of file CSCDQM_Detector.h.

Referenced by Detector().

partitions_y

unsigned int cscdqm::Detector::partitions_y

private

Number of partitions in Y axis

Definition at line 216 of file CSCDQM_Detector.h.

Referenced by Detector().

station_area

float cscdqm::Detector::station_area[4]

private

Station areas precalculated

Definition at line 210 of file CSCDQM_Detector.h.

Referenced by Area(), and Detector().