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

i
int i
Definition: DBlmapReader.cc:9

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

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

cuy.ib
int ib
Definition: cuy.py:660

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:462

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

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

edm::Handle
Definition: AssociativeIterator.h:47

ZDCDataFrame
Definition: ZDCDataFrame.h:15

ZDCTask
Definition: ZDCTask.h:22

HLT_25ns10e33_v2_cff.InputTag
tuple InputTag
Definition: HLT_25ns10e33_v2_cff.py:58467

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

edm::LogError
Definition: MessageLogger.h:164

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

j
int j
Definition: DBlmapReader.cc:9

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:109

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

DetId
Definition: DetId.h:18

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

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

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

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

alignCSCRings.e
list e
Definition: alignCSCRings.py:90

edm::InputTag
Definition: InputTag.h:15

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

alignCSCRings.r
list r
Definition: alignCSCRings.py:92

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

edm::ParameterSet
Definition: ParameterSet.h:36

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

edm::Event
Definition: Event.h:65

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:947

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

hcaldqm::quantity::ValueQuantity
Definition: ValueQuantity.h:88

edm::Run
Definition: Run.h:43

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