hcaltb::HcalTBTDCUnpacker Class Reference

#include <HcalTBTDCUnpacker.h>

Classes
struct	Hit
struct	WireChamberRecoData

Public Member Functions
	HcalTBTDCUnpacker (bool)
void	setCalib (const std::vector< std::vector< std::string > > &calibLines_)
void	unpack (const FEDRawData &raw, HcalTBEventPosition &pos, HcalTBTiming &timing) const

Private Member Functions
void	reconstructTiming (const std::vector< Hit > &hits, HcalTBTiming &timing) const
void	reconstructWC (const std::vector< Hit > &hits, HcalTBEventPosition &pos) const
void	unpackHits (const FEDRawData &raw, std::vector< Hit > &hits, HcalTBTiming &timing) const

Private Attributes
FILE *	dumpObs_
bool	includeUnmatchedHits_
double	tdc_convers [161]
double	tdc_ped [161]
struct hcaltb::HcalTBTDCUnpacker::WireChamberRecoData	wc_ [PLANECOUNT]

Static Private Attributes
static const int	PLANECOUNT = 16
static const int	WC_CHANNELIDS [PLANECOUNT *3]

Detailed Description

Author

J. Mans, P. Dudero - Minnesota

Definition at line 12 of file HcalTBTDCUnpacker.h.

Constructor & Destructor Documentation

HcalTBTDCUnpacker()

hcaltb::HcalTBTDCUnpacker::HcalTBTDCUnpacker

(

bool

include_unmatched_hits

)

Definition at line 30 of file HcalTBTDCUnpacker.cc.

      : includeUnmatchedHits_(include_unmatched_hits), dumpObs_(nullptr) {
    //  setupWC(); reads it from configuration file
    //  dumpObs_=fopen("dump_obs.csv","w");
  }

Member Function Documentation

reconstructTiming()

void hcaltb::HcalTBTDCUnpacker::reconstructTiming ( const std::vector< Hit > &  hits, HcalTBTiming &  timing  ) const

private

Definition at line 157 of file HcalTBTDCUnpacker.cc.

                                                                                                  {
    std::vector<Hit>::const_iterator j;
    double trigger_time = 0;
    double ttc_l1a_time = 0;
    double laser_flash = 0;
    double qie_phase = 0;
    double TOF1S_time = 0;
    double TOF1J_time = 0;
    double TOF2S_time = 0;
    double TOF2J_time = 0;
 
    std::vector<double> m1hits, m2hits, m3hits, s1hits, s2hits, s3hits, s4hits, bh1hits, bh2hits, bh3hits, bh4hits,
        beam_coinc;
 
    for (j = hits.begin(); j != hits.end(); j++) {
      switch (j->channel) {
        case lcTriggerTime:
          trigger_time = j->time - tdc_ped[lcTriggerTime];
          break;
        case lcTTCL1ATime:
          ttc_l1a_time = j->time - tdc_ped[lcTTCL1ATime];
          break;
        case lcBeamCoincidence:
          beam_coinc.push_back(j->time - tdc_ped[lcBeamCoincidence]);
          break;
        case lcLaserFlash:
          laser_flash = j->time - tdc_ped[lcLaserFlash];
          break;
        case lcQIEPhase:
          qie_phase = j->time - tdc_ped[lcQIEPhase];
          break;
        case lcMuon1:
          m1hits.push_back(j->time - tdc_ped[lcMuon1]);
          break;
        case lcMuon2:
          m2hits.push_back(j->time - tdc_ped[lcMuon2]);
          break;
        case lcMuon3:
          m3hits.push_back(j->time - tdc_ped[lcMuon3]);
          break;
        case lcScint1:
          s1hits.push_back(j->time - tdc_ped[lcScint1]);
          break;
        case lcScint2:
          s2hits.push_back(j->time - tdc_ped[lcScint2]);
          break;
        case lcScint3:
          s3hits.push_back(j->time - tdc_ped[lcScint3]);
          break;
        case lcScint4:
          s4hits.push_back(j->time - tdc_ped[lcScint4]);
          break;
        case lcTOF1S:
          TOF1S_time = j->time - tdc_ped[lcTOF1S];
          break;
        case lcTOF1J:
          TOF1J_time = j->time - tdc_ped[lcTOF1J];
          break;
        case lcTOF2S:
          TOF2S_time = j->time - tdc_ped[lcTOF2S];
          break;
        case lcTOF2J:
          TOF2J_time = j->time - tdc_ped[lcTOF2J];
          break;
        case lcBeamHalo1:
          bh1hits.push_back(j->time - tdc_ped[lcBeamHalo1]);
          break;
        case lcBeamHalo2:
          bh2hits.push_back(j->time - tdc_ped[lcBeamHalo2]);
          break;
        case lcBeamHalo3:
          bh3hits.push_back(j->time - tdc_ped[lcBeamHalo3]);
          break;
        case lcBeamHalo4:
          bh4hits.push_back(j->time - tdc_ped[lcBeamHalo4]);
          break;
        default:
          break;
      }
    }
 
    timing.setTimes(trigger_time, ttc_l1a_time, laser_flash, qie_phase, TOF1S_time, TOF1J_time, TOF2S_time, TOF2J_time);
    timing.setHits(
        m1hits, m2hits, m3hits, s1hits, s2hits, s3hits, s4hits, bh1hits, bh2hits, bh3hits, bh4hits, beam_coinc);
  }

References hfClusterShapes_cfi::hits, dqmiolumiharvest::j, lcBeamCoincidence, lcBeamHalo1, lcBeamHalo2, lcBeamHalo3, lcBeamHalo4, lcLaserFlash, lcMuon1, lcMuon2, lcMuon3, lcQIEPhase, lcScint1, lcScint2, lcScint3, lcScint4, lcTOF1J, lcTOF1S, lcTOF2J, lcTOF2S, lcTriggerTime, lcTTCL1ATime, tdc_ped, and ecalMatacq_cfi::timing.

Referenced by unpack().

reconstructWC()

void hcaltb::HcalTBTDCUnpacker::reconstructWC ( const std::vector< Hit > &  hits, HcalTBEventPosition &  pos  ) const

private

Definition at line 297 of file HcalTBTDCUnpacker.cc.

                                                                                                 {
    // process all planes, looping over all hits...
    const int MAX_HITS = 100;
    float hits1[MAX_HITS], hits2[MAX_HITS], hitsA[MAX_HITS];
    int n1, n2, nA, chan1, chan2, chanA;
 
    std::vector<double> x;
 
    for (int plane = 0; plane < PLANECOUNT; plane++) {
      n1 = 0;
      n2 = 0;
      nA = 0;
 
      std::vector<double> plane_hits;
      double hit_time;
      hits1[0] = -1000;
      hits2[0] = -1000;
      hitsA[0] = -1000;
      hitsA[1] = -1000;
 
      chan1 = WC_CHANNELIDS[plane * 3];
      chan2 = WC_CHANNELIDS[plane * 3 + 1];
      chanA = WC_CHANNELIDS[plane * 3 + 2];
 
      for (std::vector<Hit>::const_iterator j = hits.begin(); j != hits.end(); j++) {
        if (j->channel == chan1 && n1 < MAX_HITS) {
          hits1[n1] = j->time - TDC_OFFSET_CONSTANT;
          n1++;
        }
        if (j->channel == chan2 && n2 < MAX_HITS) {
          hits2[n2] = j->time - TDC_OFFSET_CONSTANT;
          n2++;
        }
        if (j->channel == chanA && nA < MAX_HITS) {
          hitsA[nA] = j->time - TDC_OFFSET_CONSTANT;
          nA++;
        }
      }
 
      if (n1 != 0 && n2 != 0 && dumpObs_ != nullptr) {
        fprintf(dumpObs_, "%d,%f,%f,%f,%f\n", plane, hits1[0], hits2[0], hitsA[0], hitsA[1]);
        fflush(dumpObs_);
      }
 
      // anode-matched hits
      for (int ii = 0; ii < n1; ii++) {
        int jmin = -1, lmin = -1;
        float dsumMin = 99999;
        for (int jj = 0; jj < n2; jj++) {
          for (int ll = 0; ll < nA; ll++) {
            float dsum = fabs(wc_[plane].mean - hits1[ii] - hits2[jj] + 2.0 * hitsA[ll]);
            if (dsum < (N_SIGMA * wc_[plane].sigma) && dsum < dsumMin) {
              jmin = jj;
              lmin = ll;
              dsumMin = dsum;
            }
          }
        }
        if (jmin >= 0) {
          hit_time = wc_[plane].b0 + wc_[plane].b1 * (hits1[ii] - hits2[jmin]);
          if ((plane % 2) == 0) {
            plane_hits.push_back(-hit_time);
          } else {
            plane_hits.push_back(hit_time);
          }
          hits1[ii] = -99999;
          hits2[jmin] = -99999;
          hitsA[lmin] = 99999;
        }
      }
 
      if (includeUnmatchedHits_ || plane > 9)  // unmatched hits (all pairs get in here)
        for (int ii = 0; ii < n1; ii++) {
          if (hits1[ii] < -99990)
            continue;
          for (int jj = 0; jj < n2; jj++) {
            if (hits2[jj] < -99990)
              continue;
            hit_time = wc_[plane].b0 + wc_[plane].b1 * (hits1[ii] - hits2[jj]);
            if ((plane % 2) == 0) {
              plane_hits.push_back(-hit_time);
            } else {
              plane_hits.push_back(hit_time);
            }
          }
        }
 
      if ((plane % 2) == 0)
        x = plane_hits;
      else {
        char chamber = 'A' + plane / 2;
        ep.setChamberHits(chamber, x, plane_hits);
      }
    }
  }

References hcaltb::HcalTBTDCUnpacker::WireChamberRecoData::b0, hcaltb::HcalTBTDCUnpacker::WireChamberRecoData::b1, relativeConstraints::chamber, dumpObs_, SiStripBadComponentsDQMServiceTemplate_cfg::ep, hfClusterShapes_cfi::hits, cuy::ii, includeUnmatchedHits_, dqmiolumiharvest::j, findQualityFiles::jj, SiStripPI::mean, hcaltb::N_SIGMA, PLANECOUNT, hcaltb::TDC_OFFSET_CONSTANT, wc_, and WC_CHANNELIDS.

Referenced by unpack().

setCalib()

void hcaltb::HcalTBTDCUnpacker::setCalib

(

const std::vector< std::vector< std::string > > &

calibLines_

)

Definition at line 35 of file HcalTBTDCUnpacker.cc.

                                                                                      {
    for (int i = 0; i < 161; i++) {
      tdc_ped[i] = 0.;
      tdc_convers[i] = 1.;
    }
    for (unsigned int ii = 0; ii < calibLines_.size(); ii++) {
      //   TDC configuration
      if (calibLines_[ii][0] == "TDC") {
        if (calibLines_[ii].size() == 4) {
          int channel = atoi(calibLines_[ii][1].c_str());
          tdc_ped[channel] = atof(calibLines_[ii][2].c_str());
          tdc_convers[channel] = atof(calibLines_[ii][3].c_str());
          //        printf("Got TDC %i ped %f , conversion %f\n",channel, tdc_ped[channel],tdc_convers[channel]);
        } else {
          throw cms::Exception("Incomplete configuration")
              << "Wrong TDC configuration format : expected 3 parameters, got " << calibLines_[ii].size() - 1;
        }
      }  // End of the TDCs
 
      //   Wire chambers calibration
      if (calibLines_[ii][0] == "WC") {
        if (calibLines_[ii].size() == 6) {
          int plane = atoi(calibLines_[ii][1].c_str());
          wc_[plane].b0 = atof(calibLines_[ii][2].c_str());
          wc_[plane].b1 = atof(calibLines_[ii][3].c_str());
          wc_[plane].mean = atof(calibLines_[ii][4].c_str());
          wc_[plane].sigma = atof(calibLines_[ii][5].c_str());
          //         printf("Got WC plane %i b0 %f, b1 %f, mean %f, sigma %f\n",plane,
          //         wc_[plane].b0,wc_[plane].b1,wc_[plane].mean,wc_[plane].sigma);
        } else {
          throw cms::Exception("Incomplete configuration")
              << "Wrong Wire Chamber configuration format : expected 5 parameters, got " << calibLines_[ii].size() - 1;
        }
      }  // End of the Wire Chambers
 
    }  // End of the CalibLines
  }

References hcaltb::HcalTBTDCUnpacker::WireChamberRecoData::b0, hcaltb::HcalTBTDCUnpacker::WireChamberRecoData::b1, Exception, mps_fire::i, cuy::ii, hcaltb::HcalTBTDCUnpacker::WireChamberRecoData::mean, hcaltb::HcalTBTDCUnpacker::WireChamberRecoData::sigma, findQualityFiles::size, tdc_convers, tdc_ped, and wc_.

Referenced by HcalTBObjectUnpacker::HcalTBObjectUnpacker().

unpack()

void hcaltb::HcalTBTDCUnpacker::unpack	(	const FEDRawData &	raw,
		HcalTBEventPosition &	pos,
		HcalTBTiming &	timing
	)		const

Definition at line 73 of file HcalTBTDCUnpacker.cc.

                                                                                                            {
    std::vector<Hit> hits;
 
    unpackHits(raw, hits, timing);
 
    reconstructWC(hits, pos);
    reconstructTiming(hits, timing);
  }

References hfClusterShapes_cfi::hits, reconstructTiming(), reconstructWC(), ecalMatacq_cfi::timing, and unpackHits().

Referenced by HcalTBObjectUnpacker::produce().

unpackHits()

void hcaltb::HcalTBTDCUnpacker::unpackHits ( const FEDRawData &  raw, std::vector< Hit > &  hits, HcalTBTiming &  timing  ) const

private

Definition at line 98 of file HcalTBTDCUnpacker.cc.

                                                                                                            {
    const ClassicTDCDataFormat* tdc = (const ClassicTDCDataFormat*)raw.data();
 
    if (raw.size() < 3 * 8) {
      throw cms::Exception("Missing Data") << "No data in the TDC block";
    }
 
    const unsigned int* hitbase = nullptr;
    unsigned int totalhits = 0;
 
    // old TDC (767)
    if (tdc->n_max_hits != 192) {
      const CombinedTDCQDCDataFormat* qdctdc = (const CombinedTDCQDCDataFormat*)raw.data();
      hitbase = (const unsigned int*)(qdctdc);
      hitbase += 6;                             // header
      hitbase += qdctdc->n_qdc_hits / 2;        // two unsigned short per unsigned long
      totalhits = qdctdc->n_tdc_hits & 0xFFFF;  // mask off high bits
 
      //    for (unsigned int i=0; i<qdctdc->n_qdc_hits; i++)
      //      printf("QADC: %02d %d\n",i,qdctdc->qdc_values[i]&0xFFF);
 
    } else {
      hitbase = &(tdc->hits[0]);
      totalhits = tdc->n_hits;
    }
 
    for (unsigned int i = 0; i < totalhits; i++) {
      Hit h;
      h.channel = (hitbase[i] & 0x7FC00000) >> 22;  // hardcode channel assignment
      h.time = (hitbase[i] & 0xFFFFF) * tdc_convers[h.channel];
      hits.push_back(h);
      //        printf("V767: %d %f\n",h.channel,h.time);
    }
 
    // new TDC (V775)
    int v775[32];
    for (int i = 0; i < 32; i++)
      v775[i] = -1;
    if (tdc->n_max_hits != 192) {
      const CombinedTDCQDCDataFormat* qdctdc = (const CombinedTDCQDCDataFormat*)raw.data();
      hitbase = (const unsigned int*)(qdctdc);
      hitbase += 6;                                         // header
      hitbase += qdctdc->n_qdc_hits / 2;                    // two unsigned short per unsigned long
      hitbase += (qdctdc->n_tdc_hits & 0xFFFF);             // same length
      totalhits = (qdctdc->n_tdc_hits & 0xFFFF0000) >> 16;  // mask off high bits
      for (unsigned int i = 0; i < totalhits; i++) {
        Hit h;
        //      h.channel=129+i;
        h.channel = 129 + ((hitbase[i] & 0x3F0000) >> 16);
        h.time = (hitbase[i] & 0xFFF) * tdc_convers[h.channel];
        hits.push_back(h);
        if ((h.channel - 129) < 32)
          v775[(h.channel - 129)] = (hitbase[i] & 0xFFF);
        //      printf("V775: %d %f\n",h.channel,h.time);
      }
    }
    timing.setV775(v775);
  }

References FEDRawData::data(), Exception, h, hfClusterShapes_cfi::hits, hcaltb::ClassicTDCDataFormat::hits, mps_fire::i, hcaltb::ClassicTDCDataFormat::n_hits, hcaltb::ClassicTDCDataFormat::n_max_hits, hcaltb::CombinedTDCQDCDataFormat::n_qdc_hits, hcaltb::CombinedTDCQDCDataFormat::n_tdc_hits, FEDRawData::size(), tdc_convers, and ecalMatacq_cfi::timing.

Referenced by unpack().

Member Data Documentation

dumpObs_

FILE* hcaltb::HcalTBTDCUnpacker::dumpObs_

private

Definition at line 38 of file HcalTBTDCUnpacker.h.

Referenced by reconstructWC().

includeUnmatchedHits_

bool hcaltb::HcalTBTDCUnpacker::includeUnmatchedHits_

private

Definition at line 35 of file HcalTBTDCUnpacker.h.

Referenced by reconstructWC().

PLANECOUNT

const int hcaltb::HcalTBTDCUnpacker::PLANECOUNT = 16

staticprivate

Definition at line 29 of file HcalTBTDCUnpacker.h.

Referenced by reconstructWC().

tdc_convers

double hcaltb::HcalTBTDCUnpacker::tdc_convers[161]

private

Definition at line 37 of file HcalTBTDCUnpacker.h.

Referenced by setCalib(), and unpackHits().

tdc_ped

double hcaltb::HcalTBTDCUnpacker::tdc_ped[161]

private

Definition at line 36 of file HcalTBTDCUnpacker.h.

Referenced by reconstructTiming(), and setCalib().

wc_

struct hcaltb::HcalTBTDCUnpacker::WireChamberRecoData hcaltb::HcalTBTDCUnpacker::wc_[PLANECOUNT]

private

Referenced by reconstructWC(), and setCalib().

WC_CHANNELIDS

const int hcaltb::HcalTBTDCUnpacker::WC_CHANNELIDS

staticprivate

Initial value:

= {
      12,  13,  14,   
      10,  11,  14,   
      22,  23,  24,   
      20,  21,  24,   
      32,  33,  34,   
      30,  31,  34,   
      101, 102, 104,  
      107, 108, 110,  
      113, 114, 116,  
      97,  98,  99,   
      42,  43,  -1,   
      44,  60,  -1,   
      40,  41,  -1,   
      45,  61,  -1,   
      52,  53,  -1,   
      54,  62,  -1    
  }

Definition at line 30 of file HcalTBTDCUnpacker.h.

Referenced by reconstructWC().