d1/d40/ZDCTask_8cc_source.html

 #include "DQM/HcalTasks/interface/ZDCTask.h"
 #include <map>

 using namespace hcaldqm;
 using namespace hcaldqm::constants;
 ZDCTask::ZDCTask(edm::ParameterSet const& ps)
 {

     //  tags
     _tagQIE10 = ps.getUntrackedParameter<edm::InputTag>("tagQIE10",
         edm::InputTag("hcalDigis"));
     _tokQIE10 = consumes<ZDCDigiCollection>(_tagQIE10);

     //  cuts
     _cut = ps.getUntrackedParameter<double>("cut", 50.0);
     _ped = ps.getUntrackedParameter<int>("ped", 4);
 }

 /* virtual */ void ZDCTask::bookHistograms(DQMStore::IBooker &ib,
     edm::Run const& r, edm::EventSetup const& es)
 {

     //############################## hardcode manually the zdc mapping #############################
     //############################# this follows from https://github.com/cms-sw/cmssw/blob/CMSSW_8_0_X/EventFilter/CastorRawToDigi/src/ZdcUnpacker.cc#L118
     //##############################################################################################
     //PZDC
     std::map<HcalElectronicsId,DetId> myEMap;
     HcalElectronicsId eid = HcalElectronicsId(0, 1, 0, 3);
     eid.setHTR(18, 8, 1);
     myEMap[eid]=DetId(0x54000051);//PZDC EM1

     eid = HcalElectronicsId(1, 1, 0, 3);
     eid.setHTR(18, 8, 1);
     myEMap[eid]=DetId(0x54000052);//PZDC EM2

     eid = HcalElectronicsId(2, 1, 0, 3);
     eid.setHTR(18, 8, 1);
     myEMap[eid]=DetId(0x54000053);//PZDC EM3

     eid = HcalElectronicsId(0, 2, 0, 3);
     eid.setHTR(18, 8, 1);
     myEMap[eid]=DetId(0x54000061);//PZDC HAD1

     eid = HcalElectronicsId(1, 2, 0, 3);
     eid.setHTR(18, 8, 1);
     myEMap[eid]=DetId(0x54000054);//PZDC EM4

     eid = HcalElectronicsId(2, 2, 0, 3);
     eid.setHTR(18, 8, 1);
     myEMap[eid]=DetId(0x54000055);//PZDC EM5

     eid = HcalElectronicsId(0, 3, 0, 3);
     eid.setHTR(18, 8, 1);
     myEMap[eid]=DetId(0x54000062);//PZDC HAD2

     eid = HcalElectronicsId(1, 3, 0, 3);
     eid.setHTR(18, 8, 1);
     myEMap[eid]=DetId(0x54000063);//PZDC HAD3

     eid = HcalElectronicsId(2, 3, 0, 3);
     eid.setHTR(18, 8, 1);
     myEMap[eid]=DetId(0x54000064);//PZDC HAD4

     //NZDC
     eid = HcalElectronicsId(0,1,1,3);
     eid.setHTR(18, 8, 0);
     myEMap[eid]=DetId(0x54000011);//NZDC EM1

     eid = HcalElectronicsId(1,1,1,3);
     eid.setHTR(18, 8, 0);
     myEMap[eid]=DetId(0x54000012);//NZDC EM2

     eid = HcalElectronicsId(2,1,1,3);
     eid.setHTR(18, 8, 0);
     myEMap[eid]=DetId(0x54000013);//NZDC EM3

     eid = HcalElectronicsId(0,2,1,3);
     eid.setHTR(18, 8, 0);
     myEMap[eid]=DetId(0x54000015);//NZDC EM5

     eid = HcalElectronicsId(1,2,1,3);
     eid.setHTR(18, 8, 0);
     myEMap[eid]=DetId(0x54000021);//NZDC HAD1

     eid = HcalElectronicsId(2,2,1,3);
     eid.setHTR(18, 8, 0);
     myEMap[eid]=DetId(0x54000014);//NZDC EM4

     eid = HcalElectronicsId(0,3,1,3);
     eid.setHTR(18, 8, 0);
     myEMap[eid]=DetId(0x54000022);//NZDC HAD2

     eid = HcalElectronicsId(1,3,1,3);
     eid.setHTR(18, 8, 0);
     myEMap[eid]=DetId(0x54000023);//NZDC HAD3

     eid = HcalElectronicsId(2,3,1,3);
     eid.setHTR(18, 8, 0);
     myEMap[eid]=DetId(0x54000024);//NZDC HAD4
     //##################################### end hardcoding ###################################

     ib.cd();

     //quantities for axis
     hcaldqm::quantity::ValueQuantity xAxisShape(hcaldqm::quantity::fTiming_TS);
     hcaldqm::quantity::ValueQuantity yAxisShape(hcaldqm::quantity::ffC_10000);

     hcaldqm::quantity::ValueQuantity xAxisADC(hcaldqm::quantity::fADC_128);

     //book histos per channel
     for (std::map<HcalElectronicsId,DetId>::const_iterator itr=myEMap.begin(); itr!=myEMap.end(); ++itr)
       {
         char histoname[300];

         sprintf(histoname,"%d_%d_%d_%d",itr->first.fiberChanId(),itr->first.fiberIndex(),itr->first.spigot(),itr->first.dccid());

         ib.setCurrentFolder("Hcal/ZDCTask/Shape_perChannel");
         _cShape_EChannel[histoname] = ib.bookProfile(histoname,histoname,xAxisShape.nbins(),xAxisShape.min(),xAxisShape.max(),yAxisShape.nbins(),yAxisShape.min(),yAxisShape.max());
         _cShape_EChannel[histoname]->setAxisTitle("Timing",1);
         _cShape_EChannel[histoname]->setAxisTitle("fC QIE8",2);

         ib.setCurrentFolder("Hcal/ZDCTask/ADC_perChannel");
         _cADC_EChannel[histoname] = ib.book1D(histoname,histoname,xAxisADC.nbins(),xAxisADC.min(),xAxisADC.max());
         _cADC_EChannel[histoname]->getRootObject()->SetBit(BIT(hcaldqm::constants::BIT_OFFSET+hcaldqm::quantity::AxisType::fYAxis));
         _cADC_EChannel[histoname]->setAxisTitle("ADC QIE8",1);

         ib.setCurrentFolder("Hcal/ZDCTask/ADC_vs_TS_perChannel");
         _cADC_vs_TS_EChannel[histoname] = ib.book2D(histoname,histoname,xAxisShape.nbins(),xAxisShape.min(),xAxisShape.max(),xAxisADC.nbins(),xAxisADC.min(),xAxisADC.max());
         _cADC_vs_TS_EChannel[histoname]->getRootObject()->SetBit(BIT(hcaldqm::constants::BIT_OFFSET+hcaldqm::quantity::AxisType::fYAxis));
         _cADC_vs_TS_EChannel[histoname]->setAxisTitle("Timing",1);
         _cADC_vs_TS_EChannel[histoname]->setAxisTitle("ADC QIE8",2);


       }

     //book global histos
     ib.setCurrentFolder("Hcal/ZDCTask");

     _cShape = ib.bookProfile("Shape","Shape",xAxisShape.nbins(),xAxisShape.min(),xAxisShape.max(),yAxisShape.nbins(),yAxisShape.min(),yAxisShape.max());
     _cShape->setAxisTitle("Timing",1);
     _cShape->setAxisTitle("fC QIE8",2);

     _cADC = ib.book1D("ADC","ADC",xAxisADC.nbins(),xAxisADC.min(),xAxisADC.max());
     _cADC->getRootObject()->SetBit(BIT(hcaldqm::constants::BIT_OFFSET+hcaldqm::quantity::AxisType::fYAxis));
     _cADC->setAxisTitle("ADC QIE8",1);

     _cADC_vs_TS = ib.book2D("ADC_vs_TS","ADC_vs_TS",xAxisShape.nbins(),xAxisShape.min(),xAxisShape.max(),xAxisADC.nbins(),xAxisADC.min(),xAxisADC.max());
     _cADC_vs_TS->getRootObject()->SetBit(BIT(hcaldqm::constants::BIT_OFFSET+hcaldqm::quantity::AxisType::fYAxis));
     _cADC_vs_TS->setAxisTitle("Timing",1);
     _cADC_vs_TS->setAxisTitle("ADC QIE8",2);

 }


 /* virtual */ void ZDCTask::analyze(edm::Event const& e, edm::EventSetup const&)
 {
     edm::Handle<ZDCDigiCollection> cqie10;
     if (!e.getByToken(_tokQIE10, cqie10))
       edm::LogError("Collection ZDCDigiCollection isn't available"
             + _tagQIE10.label() + " " + _tagQIE10.instance());


     for (uint32_t i=0; i<cqie10->size(); i++)
       {
         ZDCDataFrame frame = static_cast<ZDCDataFrame>((*cqie10)[i]);
         HcalElectronicsId eid = frame.elecId();

         char histoname[300];
             sprintf(histoname,"%d_%d_%d_%d",eid.fiberChanId(),eid.fiberIndex(),eid.spigot(),eid.dccid());

         //  compute the signal, ped subracted
         //double q = hcaldqm::utilities::sumQ_v10<ZDCDataFrame>(frame, constants::adc2fC[_ped], 0, frame.size()-1);

         //  iterate thru all TS and fill
         for (int j=0; j<frame.size(); j++)
           {
         _cShape_EChannel[histoname]->Fill(j, frame[j].nominal_fC());
         _cShape->Fill(j, frame[j].nominal_fC());

         _cADC_EChannel[histoname]->Fill(frame[j].adc());
         _cADC->Fill(frame[j].adc());

         _cADC_vs_TS_EChannel[histoname]->Fill(j,frame[j].adc());
         _cADC_vs_TS->Fill(j,frame[j].adc());


           }
       }
 }


 DEFINE_FWK_MODULE(ZDCTask);
ecalMGPA::adc
int adc(sample_type sample)
get the ADC sample (12 bits)
Definition: EcalMGPASample.h:11

edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const

HcalElectronicsId::fiberIndex
int fiberIndex() const
get the fiber index. For VME 1-8 (which of eight fibers carried by a spigot), for uTCA fibers are zer...
Definition: HcalElectronicsId.h:61

mps_fire.i
i
Definition: mps_fire.py:156

hcaldqm::quantity::fTiming_TS
Definition: ValueQuantity.h:25

DQMStore::IBooker::bookProfile
MonitorElement * bookProfile(Args &&...args)
Definition: DQMStore.h:157

edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:457

DQMStore::IBooker::cd
void cd(void)
Definition: DQMStore.cc:269

DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17

edm::Handle
Definition: AssociativeIterator.h:47

ZDCDataFrame
Definition: ZDCDataFrame.h:15

MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37

ZDCTask
Definition: ZDCTask.py:1

ZDCTask::analyze
virtual void analyze(edm::Event const &, edm::EventSetup const &)
Definition: ZDCTask.cc:157

edm::LogError
Definition: MessageLogger.h:166

ZDCTask::ZDCTask
ZDCTask(edm::ParameterSet const &)
Definition: ZDCTask.cc:7

HcalElectronicsId::dccid
int dccid() const
get the (Hcal local) DCC id for VME, crate number for uTCA
Definition: HcalElectronicsId.h:76

DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115

edm::EventSetup
Definition: EventSetup.h:45

HcalElectronicsId::fiberChanId
int fiberChanId() const
get the fiber channel id (which of channels on a fiber)
Definition: HcalElectronicsId.h:59

hcaldqm::quantity::ValueQuantity::min
virtual double min()
Definition: ValueQuantity.h:281

HcalElectronicsId::spigot
int spigot() const
get the spigot (input number on DCC, AMC card number for uTCA)
Definition: HcalElectronicsId.h:74

alignCSCRings.r
r
Definition: alignCSCRings.py:92

DetId
Definition: DetId.h:18

HcalElectronicsId::setHTR
void setHTR(int crate, int slot, int tb)
Definition: HcalElectronicsId.cc:52

hcaldqm::constants
Definition: Constants.h:10

MonitorElement::getRootObject
TObject * getRootObject(void) const
Definition: MonitorElement.cc:1451

DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:277

hcaldqm
Definition: Constants.h:8

DQMStore::IBooker::book2D
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:133

ZDCTask.h

hcaldqm::quantity::fADC_128
Definition: ValueQuantity.h:19

hcaldqm::quantity::fYAxis
Definition: Quantity.h:31

ZDCDataFrame::size
int size() const
total number of samples in the digi
Definition: ZDCDataFrame.h:26

runTauDisplay.eid
eid
Definition: runTauDisplay.py:295

edm::InputTag
Definition: InputTag.h:15

edm::SortedCollection::size
size_type size() const
Definition: SortedCollection.h:236

ZDCDataFrame::elecId
const HcalElectronicsId & elecId() const
Definition: ZDCDataFrame.h:23

ZDCTask::bookHistograms
virtual void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)
Definition: ZDCTask.cc:20

hcaldqm::constants::BIT_OFFSET
int const BIT_OFFSET
Definition: Constants.h:271

edm::ParameterSet
Definition: ParameterSet.h:36

hcaldqm::quantity::ValueQuantity::max
virtual double max()
Definition: ValueQuantity.h:282

edm::Event
Definition: Event.h:66

hcaldqm::quantity::ffC_10000
Definition: ValueQuantity.h:22

DQMStore::IBooker
Definition: DQMStore.h:90

MonitorElement::setAxisTitle
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:955

HcalElectronicsId
Readout chain identification for Hcal.
Definition: HcalElectronicsId.h:32

hcaldqm::quantity::ValueQuantity
Definition: ValueQuantity.h:260

edm::Run
Definition: Run.h:42

cuy.ib
ib
Definition: cuy.py:660

hcaldqm::quantity::ValueQuantity::nbins
virtual int nbins()
Definition: ValueQuantity.h:280