CMS 3D CMS Logo
List of all members | Public Member Functions | Private Attributes
DTPreCalibrationTask Class Reference

#include <DTPreCalibrationTask.h>

Inheritance diagram for DTPreCalibrationTask:
one::DQMEDAnalyzer< T > one::dqmimplementation::DQMBaseClass< T... >

Public Member Functions

void analyze (const edm::Event &, const edm::EventSetup &) override
 
void bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
 
void bookOccupancyPlot (DQMStore::IBooker &, int wheel, int sector)
 
void bookTimeBoxes (DQMStore::IBooker &, int wheel, int sector)
 Book histos. More...
 
 DTPreCalibrationTask (const edm::ParameterSet &ps)
 Constructor. More...
 
 ~DTPreCalibrationTask () override
 Destructor. More...
 
- Public Member Functions inherited from one::DQMEDAnalyzer< T >
 DQMEDAnalyzer ()=default
 
 DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
 
 DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
 
 ~DQMEDAnalyzer () override=default
 

Private Attributes

edm::EDGetTokenT< DTDigiCollection > digiLabel
 
std::string folderName
 
int maxTriggerWidth
 
int minTriggerWidth
 
std::map< std::pair< int, int >, MonitorElement * > OccupancyHistos
 
std::map< std::pair< int, int >, MonitorElement * > TimeBoxes
 

Detailed Description

Analysis on DT digis (TB + occupancy) before the calibration step

Author
G. Mila - INFN Torino

Definition at line 30 of file DTPreCalibrationTask.h.

Constructor & Destructor Documentation

DTPreCalibrationTask::DTPreCalibrationTask ( const edm::ParameterSet &  ps)

Constructor.

Definition at line 34 of file DTPreCalibrationTask.cc.

References CastorSimpleReconstructor_cfi::digiLabel, edm::ParameterSet::getUntrackedParameter(), and LogTrace.

34  {
35 
36  LogTrace("DTPreCalibSummary") <<"[DTPrecalibrationTask]: Constructor"<<endl;
37 
38  // Label to retrieve DT digis from the event
39  digiLabel = consumes<DTDigiCollection>(ps.getUntrackedParameter<string>("digiLabel"));
40 
41  // parameter for Time Boxes booking
42  minTriggerWidth = ps.getUntrackedParameter<int>("minTriggerWidth",2000);
43  maxTriggerWidth = ps.getUntrackedParameter<int>("maxTriggerWidth",6000);
44 
45  // get the histo folder name
46  folderName = ps.getUntrackedParameter<string>("folderName");
47 
48 }
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
DTPreCalibrationTask::digiLabel
edm::EDGetTokenT< DTDigiCollection > digiLabel
Definition: DTPreCalibrationTask.h:48
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:602
DTPreCalibrationTask::~DTPreCalibrationTask ( )
override

Destructor.

Definition at line 51 of file DTPreCalibrationTask.cc.

51 {}

Member Function Documentation

void DTPreCalibrationTask::analyze ( const edm::Event &  event,
const edm::EventSetup &  setup 
)
override

Definition at line 70 of file DTPreCalibrationTask.cc.

References CastorSimpleReconstructor_cfi::digiLabel, HcalObjRepresent::Fill(), DTLayerId::layer(), edm::second(), DTChamberId::station(), DTSuperLayerId::superlayer(), and cuy::yBin.

70  {
71 
72  // Get the digis from the event
73  edm::Handle<DTDigiCollection> dtdigis;
74  event.getByToken(digiLabel, dtdigis);
75 
76  // LOOP OVER ALL THE DIGIS OF THE EVENT
77  DTDigiCollection::DigiRangeIterator dtLayerId_It;
78  for (dtLayerId_It=dtdigis->begin(); dtLayerId_It!=dtdigis->end(); ++dtLayerId_It){
79  for (DTDigiCollection::const_iterator digiIt = ((*dtLayerId_It).second).first;
80  digiIt!=((*dtLayerId_It).second).second; ++digiIt){
81 
82  //Fill the Time Boxes
83  int tdcTime = (*digiIt).countsTDC();
84  TimeBoxes[make_pair((*dtLayerId_It).first.superlayerId().chamberId().wheel(),
85  (*dtLayerId_It).first.superlayerId().chamberId().sector())]->Fill(tdcTime);
86 
87  //Fill the occupancy plot
88  const DTLayerId dtLId = (*dtLayerId_It).first;
89  int yBin = (dtLId.station()-1)*12+dtLId.layer()+4*(dtLId.superlayer()-1);
90  if(dtLId.station()==4 && dtLId.superlayer()==3)
91  yBin = (dtLId.station()-1)*12+dtLId.layer()+4*(dtLId.superlayer()-2);
92  if((*dtLayerId_It).first.superlayerId().chamberId().sector()<13)
93  OccupancyHistos[make_pair((*dtLayerId_It).first.superlayerId().chamberId().wheel(),
94  (*dtLayerId_It).first.superlayerId().chamberId().sector())]->Fill((*digiIt).wire(),yBin);
95  else{
96  if(dtLId.superlayer()!=3) yBin = 44 + dtLId.layer();
97  else yBin = 48 + dtLId.layer();
98  if((*dtLayerId_It).first.superlayerId().chamberId().sector()==13)
99  OccupancyHistos[make_pair((*dtLayerId_It).first.superlayerId().chamberId().wheel(),4)]->Fill((*digiIt).wire(),yBin);
100  if((*dtLayerId_It).first.superlayerId().chamberId().sector()==14)
101  OccupancyHistos[make_pair((*dtLayerId_It).first.superlayerId().chamberId().wheel(),10)]->Fill((*digiIt).wire(),yBin);
102  }
103 
104  }
105  }
106 
107 }
DTPreCalibrationTask::digiLabel
edm::EDGetTokenT< DTDigiCollection > digiLabel
Definition: DTPreCalibrationTask.h:48
DTLayerId::layer
int layer() const
Return the layer number.
Definition: DTLayerId.h:53
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:235
HcalObjRepresent::Fill
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:1113
DTSuperLayerId::superlayer
int superlayer() const
Return the superlayer number (deprecated method name)
Definition: DTSuperLayerId.h:55
MuonDigiCollection< DTLayerId, DTDigi >::const_iterator
std::vector< DTDigi >::const_iterator const_iterator
Definition: MuonDigiCollection.h:110
DTChamberId::station
int station() const
Return the station number.
Definition: DTChamberId.h:51
DTPreCalibrationTask::OccupancyHistos
std::map< std::pair< int, int >, MonitorElement * > OccupancyHistos
Definition: DTPreCalibrationTask.h:57
cuy.yBin
yBin
Definition: cuy.py:893
DTPreCalibrationTask::TimeBoxes
std::map< std::pair< int, int >, MonitorElement * > TimeBoxes
Definition: DTPreCalibrationTask.h:54
void DTPreCalibrationTask::bookHistograms ( DQMStore::IBooker &  iBooker,
edm::Run const &  ,
edm::EventSetup const &   
)
override

Definition at line 54 of file DTPreCalibrationTask.cc.

References LogTrace, DQMStore::IBooker::setCurrentFolder(), and makeMuonMisalignmentScenario::wheel.

55  {
56 
57  for(int wheel=-2; wheel<=2; wheel++){
58  for(int sector=1; sector<=14; sector++){
59  LogTrace("DTPreCalibSummary") <<"[DTPrecalibrationTask]: Book histos for wheel "<<wheel<<", sector "<<sector<<endl;
60  iBooker.setCurrentFolder(folderName+"/TimeBoxes");
61  bookTimeBoxes(iBooker, wheel, sector);
62  iBooker.setCurrentFolder(folderName+"/OccupancyHistos");
63  if(sector<13) bookOccupancyPlot(iBooker, wheel, sector);
64  }
65  }
66 
67 }
DTPreCalibrationTask::bookTimeBoxes
void bookTimeBoxes(DQMStore::IBooker &, int wheel, int sector)
Book histos.
Definition: DTPreCalibrationTask.cc:109
DTPreCalibrationTask::bookOccupancyPlot
void bookOccupancyPlot(DQMStore::IBooker &, int wheel, int sector)
Definition: DTPreCalibrationTask.cc:122
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:268
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:602
void DTPreCalibrationTask::bookOccupancyPlot ( DQMStore::IBooker &  iBooker,
int  wheel,
int  sector 
)

Definition at line 122 of file DTPreCalibrationTask.cc.

References DQMStore::IBooker::book2D(), MonitorElement::setAxisTitle(), and makeMuonMisalignmentScenario::wheel.

122  {
123 
124  stringstream wh; wh << wheel;
125  stringstream sec; sec << sector;
126 
127  // book the occpancy plot
128  if(sector==4 || sector==10)
129  OccupancyHistos[make_pair(wheel, sector)]= iBooker.book2D("Occupancy_W"+wh.str()+"_Sec"+sec.str(), "Occupancy W"+wh.str()+"_Sec"+sec.str(),100,1,100,52,1,53);
130  else
131  OccupancyHistos[make_pair(wheel, sector)]= iBooker.book2D("Occupancy_W"+wh.str()+"_Sec"+sec.str(), "Occupancy W"+wh.str()+"_Sec"+sec.str(),100,1,100,44,1,45);
132  OccupancyHistos[make_pair(wheel, sector)]->setAxisTitle("wire number", 1);
133  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(1,"M1L1",2);
134  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(2,"M1L2",2);
135  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(3,"M1L3",2);
136  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(4,"M1L4",2);
137  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(5,"M1L5",2);
138  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(6,"M1L6",2);
139  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(7,"M1L7",2);
140  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(8,"M1L8",2);
141  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(9,"M1L9",2);
142  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(10,"M1L10",2);
143  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(11,"M1L11",2);
144  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(12,"M1L12",2);
145  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(13,"M2L1",2);
146  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(14,"M2L2",2);
147  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(15,"M2L3",2);
148  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(16,"M2L4",2);
149  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(17,"M2L5",2);
150  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(18,"M2L6",2);
151  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(19,"M2L7",2);
152  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(20,"M2L8",2);
153  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(21,"M2L9",2);
154  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(22,"M2L10",2);
155  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(23,"M2L11",2);
156  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(24,"M2L12",2);
157  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(25,"M3L1",2);
158  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(26,"M3L2",2);
159  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(27,"M3L3",2);
160  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(28,"M3L4",2);
161  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(29,"M3L5",2);
162  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(30,"M3L6",2);
163  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(31,"M3L7",2);
164  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(32,"M3L8",2);
165  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(33,"M3L9",2);
166  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(34,"M3L10",2);
167  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(35,"M3L11",2);
168  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(36,"M3L12",2);
169  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(37,"M4L1",2);
170  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(38,"M4L2",2);
171  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(39,"M4L3",2);
172  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(40,"M4L4",2);
173  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(41,"M4L5",2);
174  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(42,"M4L6",2);
175  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(43,"M4L7",2);
176  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(44,"M4L8",2);
177  if(sector==4){
178  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(45,"M4Sec13L1",2);
179  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(46,"M4Sec13L2",2);
180  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(47,"M4Sec13L3",2);
181  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(48,"M4Sec13L4",2);
182  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(49,"M4Sec13L5",2);
183  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(50,"M4Sec13L6",2);
184  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(51,"M4Sec13L7",2);
185  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(52,"M4Sec13L8",2);
186  }
187  if(sector==10){
188  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(45,"M4Sec14L1",2);
189  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(46,"M4Sec14L2",2);
190  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(47,"M4Sec14L3",2);
191  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(48,"M4Sec14L4",2);
192  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(49,"M4Sec14L5",2);
193  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(50,"M4Sec14L6",2);
194  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(51,"M4Sec14L7",2);
195  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(52,"M4Sec14L8",2);
196  }
197 
198 }
DQMStore::IBooker::book2D
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:109
DTPreCalibrationTask::OccupancyHistos
std::map< std::pair< int, int >, MonitorElement * > OccupancyHistos
Definition: DTPreCalibrationTask.h:57
void DTPreCalibrationTask::bookTimeBoxes ( DQMStore::IBooker &  iBooker,
int  wheel,
int  sector 
)

Book histos.

Definition at line 109 of file DTPreCalibrationTask.cc.

References DQMStore::IBooker::book1D(), MonitorElement::setAxisTitle(), and makeMuonMisalignmentScenario::wheel.

109  {
110 
111  stringstream wh; wh << wheel;
112  stringstream sec; sec << sector;
113 
114  // book the time boxes
115  TimeBoxes[make_pair(wheel, sector)]= iBooker.book1D("TimeBox_W"+wh.str()+"_Sec"+sec.str(), "Time Box W"+wh.str()+"_Sec"+sec.str(),(maxTriggerWidth-minTriggerWidth)/50, minTriggerWidth, maxTriggerWidth);
116  TimeBoxes[make_pair(wheel, sector)]->setAxisTitle("TDC counts");
117 
118 }
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:106
DTPreCalibrationTask::TimeBoxes
std::map< std::pair< int, int >, MonitorElement * > TimeBoxes
Definition: DTPreCalibrationTask.h:54

Member Data Documentation

edm::EDGetTokenT<DTDigiCollection> DTPreCalibrationTask::digiLabel
private

Definition at line 48 of file DTPreCalibrationTask.h.

std::string DTPreCalibrationTask::folderName
private

Definition at line 51 of file DTPreCalibrationTask.h.

int DTPreCalibrationTask::maxTriggerWidth
private

Definition at line 50 of file DTPreCalibrationTask.h.

int DTPreCalibrationTask::minTriggerWidth
private

Definition at line 49 of file DTPreCalibrationTask.h.

std::map<std::pair<int,int>, MonitorElement* > DTPreCalibrationTask::OccupancyHistos
private

Definition at line 57 of file DTPreCalibrationTask.h.

std::map<std::pair<int,int>, MonitorElement* > DTPreCalibrationTask::TimeBoxes
private

Definition at line 54 of file DTPreCalibrationTask.h.

Generated for CMSSW Reference Manual by doxygen 1.8.11