CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
HcalTBQADCUnpacker.cc
Go to the documentation of this file.
3 
4 // QADC channels
5 static const int N_QADCS_ALLOWED = 6;
6 
7 // Channel to logical unit
8 // TB04
9 static const int aScint1 = 52; // 14x14 cm
10 static const int aScint2 = 53; // 4x4 cm
11 static const int aScint3 = 54; // 2x2 cm
12 static const int aScint4 = 66; // 14x14 cm
13 static const int aMuonV = 56; // Behind the table
14 static const int aMuonV3 = 64; // on HB1
15 static const int aMuonV6 = 62; // on HB2
16 static const int aMuonVH1 = 67; // Oct. muon veto wall
17 static const int aMuonVH2 = 68; // Oct. muon veto wall
18 static const int aMuonVH3 = 69; // Oct. muon veto wall
19 static const int aMuonVH4 = 70; // Oct. muon veto wall
20 static const int aCerenkov2 = 49; // el id
21 static const int aCerenkov3 = 59; // pi/proton separation
22 static const int aSCI_VLE = 65; // VLE line
23 static const int aSCI_521 = 57;
24 static const int aSCI_528 = 58;
25 // TB06
26 static const int bMuonV1 = 0; // Muon veto wall
27 static const int bMuonV2 = 1; // Muon veto wall
28 static const int bMuonV3 = 2; // Muon veto wall
29 static const int bMuonV4 = 3; // Muon veto wall
30 static const int bMuonV5 = 4; // Muon veto wall
31 static const int bMuonV6 = 5; // Muon veto wall
32 static const int bMuonV7 = 6; // Muon veto wall
33 static const int bMuonV8 = 7; // Muon veto wall
34 static const int bMuonVF = 10; // Behind the table
35 static const int bMuonVB = 11; // Behind a beam dump
36 static const int bScint1 = 12; // 14x14 cm
37 static const int bScint2 = 13; // 4x4 cm
38 static const int bScint3 = 14; // 2x2 cm
39 static const int bScint4 = 15; // 14x14 cm
40 static const int bCerenkov1 = 16; // in HE line
41 static const int bCerenkov2 = 17; // el id in VLE
42 static const int bCerenkov3 = 18; // pi/proton id
43 static const int bTOF1S = 20; // TOF1S (S=Saleve side, TOF1 upstream)
44 static const int bTOF1J = 21; // TOF1J (J=Jura side)
45 static const int bTOF2S = 22; // TOF2S (TOF2 downstream)
46 static const int bTOF2J = 23; // TOF2J
47 static const int bSCI_521 = 24; // In HE beam line for now
48 static const int bSCI_528 = 25; // NC
49 static const int bVH2 = 27; // beam halo left from particle view
50 static const int bVH4 = 28; // beam halo down
51 static const int bVH3 = 29; // beam halo up
52 static const int bVH1 = 30; // beam halo right from particle view
53 
54 namespace hcaltb {
55 
57 
60  unsigned short data[N_QADCS_ALLOWED*32];
61  unsigned int cdfTrailer0,cdfTrailer1;
62 };
63 
66  unsigned int n_qdc_hits; // Count of QDC channels
67  unsigned int n_tdc_hits; // upper/lower TDC counts
68  unsigned short qdc_values[4];
69  };
70 
71 // Sets the pedestal and gain
72 void HcalTBQADCUnpacker::setCalib(const std::vector<std::vector<std::string> >& calibLines_) {
73 // The default pedestal and gain
74  for(int i=0;i<N_QADCS_ALLOWED*32;i++)
75  {
76  qdc_ped[i]=0.;qdc_gain[i]=1.;
77  }
78 // Pedestal and gains from configuration file.
79  for(unsigned int ii=0;ii<calibLines_.size();ii++)
80  {
81  if(calibLines_[ii][0]=="QDC")
82  {
83  if(calibLines_[ii].size()==4)
84  {
85  int channel=atoi(calibLines_[ii][1].c_str());
86  qdc_ped[channel]=atof(calibLines_[ii][2].c_str());
87  qdc_gain[channel]=atof(calibLines_[ii][3].c_str());
88  // printf("Got QDC %i ped %f , gain %f\n",channel, qdc_ped[channel],qdc_gain[channel]);
89  }
90  else
91  {
92  throw cms::Exception("Incomplete configuration") <<
93  "Wrong QADC configuration format: expected 3 parameters, got "<<calibLines_[ii].size()-1;
94  }
95  }
96  } // End of calibLines.
97  }
98 
100  HcalTBBeamCounters& beamadc, bool is04) const {
101 
102  if (raw.size()<3*8) {
103  throw cms::Exception("Missing Data") << "No data in the QDC block";
104  }
105 
106 
107  if(is04){
108  const ClassicQADCDataFormat* qadc=(const ClassicQADCDataFormat*)raw.data();
109  double qdc_calib_hits[N_QADCS_ALLOWED*32];
110  // Applying mask, pedestal subtraction and gain.
111  for (unsigned int i=0;i<N_QADCS_ALLOWED*32;i++)
112  qdc_calib_hits[i]=((qadc->data[i]&0xFFF)-qdc_ped[i])/qdc_gain[i];
113 
114  // Ecal energy sum should go here.
115  double Ecal7x7=0.;
116  for(int i=0;i<49;i++)Ecal7x7+=qdc_calib_hits[i];
117 
118 
119  beamadc.setADCs04(qdc_calib_hits[aMuonV],qdc_calib_hits[aMuonV3],qdc_calib_hits[aMuonV6],
120  qdc_calib_hits[aMuonVH1],qdc_calib_hits[aMuonVH2],qdc_calib_hits[aMuonVH3],
121  qdc_calib_hits[aMuonVH4],qdc_calib_hits[aCerenkov2],qdc_calib_hits[aCerenkov3],
122  qdc_calib_hits[aSCI_VLE],qdc_calib_hits[aSCI_521],qdc_calib_hits[aSCI_528],
123  qdc_calib_hits[aScint1],qdc_calib_hits[aScint2],qdc_calib_hits[aScint3],
124  qdc_calib_hits[aScint4],Ecal7x7);
125  }
126  else{
127  const CombinedTDCQDCDataFormat* qdctdc=(const CombinedTDCQDCDataFormat*)raw.data();
128  double qdc_calib_hits[32];
129  for (unsigned int i=0;i<qdctdc->n_qdc_hits;i++)
130  qdc_calib_hits[i]=((qdctdc->qdc_values[i]&0xFFF)-qdc_ped[i])/qdc_gain[i];
131 
132  beamadc.setADCs06( qdc_calib_hits[bMuonVF], qdc_calib_hits[bMuonVB],
133  qdc_calib_hits[bMuonV1],qdc_calib_hits[bMuonV2],qdc_calib_hits[bMuonV3],
134  qdc_calib_hits[bMuonV4],qdc_calib_hits[bMuonV5],qdc_calib_hits[bMuonV6],
135  qdc_calib_hits[bMuonV7], qdc_calib_hits[bMuonV8],
136  qdc_calib_hits[bCerenkov1], qdc_calib_hits[bCerenkov2], qdc_calib_hits[bCerenkov3],
137  qdc_calib_hits[bScint1], qdc_calib_hits[bScint2],
138  qdc_calib_hits[bScint3], qdc_calib_hits[bScint4],
139  qdc_calib_hits[bTOF1S], qdc_calib_hits[bTOF1J],
140  qdc_calib_hits[bTOF2S], qdc_calib_hits[bTOF2J],
141  qdc_calib_hits[bSCI_521],qdc_calib_hits[bSCI_528],
142  qdc_calib_hits[bVH1],qdc_calib_hits[bVH2],
143  qdc_calib_hits[bVH3],qdc_calib_hits[bVH4]);
144 
145  }
146  }
147 
148 }
static const int aSCI_521
static const int aScint4
static const int bMuonV1
int i
Definition: DBlmapReader.cc:9
static const int aMuonVH2
static const int bTOF1S
static const int aCerenkov2
static const int bVH2
static const int bTOF1J
static const int bCerenkov1
void unpack(const FEDRawData &raw, HcalTBBeamCounters &beamadc, bool is04_=true) const
static const int bScint1
int ii
Definition: cuy.py:588
static const int bMuonV3
static const int bMuonV2
size_t size() const
Lenght of the data buffer in bytes.
Definition: FEDRawData.h:47
static const int bSCI_528
static const int bTOF2J
static const int bVH1
static const int bCerenkov3
static const int aScint1
static const int aMuonV3
static const int bMuonV7
static const int bMuonVB
static const int bScint3
static const int bMuonV4
static const int aMuonVH1
static const int aMuonV6
static const int aCerenkov3
static const int aScint2
static const int bScint2
static const int aScint3
static const int bMuonV6
static const int bMuonV5
static const int bSCI_521
static const int bVH3
static const int bMuonV8
unsigned short data[N_QADCS_ALLOWED *32]
static const int bVH4
static const int bCerenkov2
static const int aMuonVH3
static const int aMuonV
void setADCs04(double VMadc, double V3adc, double V6adc, double VH1adc, double VH2adc, double VH3adc, double VH4adc, double CK2adc, double CK3adc, double SciVLEadc, double Sci521adc, double Sci528adc, double S1adc, double S2adc, double S3adc, double S4adc, double Ecal7x7)
void setCalib(const std::vector< std::vector< std::string > > &calibLines_)
const unsigned char * data() const
Return a const pointer to the beginning of the data buffer.
Definition: FEDRawData.cc:28
static const int aMuonVH4
static const int aSCI_528
void setADCs06(double VMFadc, double VMBadc, double VM1adc, double VM2adc, double VM3adc, double VM4adc, double VM5adc, double VM6adc, double VM7adc, double VM8adc, double CK1adc, double CK2adc, double CK3adc, double S1adc, double S2adc, double S3adc, double S4adc, double TOF1Sadc, double TOF1Jadc, double TOF2Sadc, double TOF2Jadc, double Sci521adc, double Sci528adc, double BH1adc, double BH2adc, double BH3adc, double BH4adc)
static const int aSCI_VLE
static const int bScint4
static const int bTOF2S
static const int N_QADCS_ALLOWED
tuple size
Write out results.
static const int bMuonVF