CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
DTPreCalibrationTask.cc
Go to the documentation of this file.
1 #include "DQMOffline/CalibMuon/interface/DTPreCalibrationTask.h"
2 
3 
4 // Framework
5 #include "FWCore/Framework/interface/Event.h"
6 #include "FWCore/Framework/interface/ESHandle.h"
7 #include <FWCore/Framework/interface/EventSetup.h>
8 #include "FWCore/ParameterSet/interface/ParameterSet.h"
9 #include "DQMServices/Core/interface/DQMStore.h"
10 #include "DQMServices/Core/interface/MonitorElement.h"
11 #include "FWCore/ServiceRegistry/interface/Service.h"
12 #include "FWCore/MessageLogger/interface/MessageLogger.h"
13 
14 // Geometry
15 #include "Geometry/DTGeometry/interface/DTLayer.h"
16 #include "Geometry/DTGeometry/interface/DTGeometry.h"
17 #include "Geometry/DTGeometry/interface/DTTopology.h"
18 #include "Geometry/Records/interface/MuonGeometryRecord.h"
19 
20 // Digis
21 #include <DataFormats/DTDigi/interface/DTDigi.h>
22 #include "DataFormats/DTDigi/interface/DTDigiCollection.h"
23 #include "CondFormats/DataRecord/interface/DTStatusFlagRcd.h"
24 #include "CondFormats/DTObjects/interface/DTStatusFlag.h"
25 #include "CondFormats/DTObjects/interface/DTReadOutMapping.h"
26 
27 
28 #include "TH2F.h"
29 #include "TFile.h"
30 
31 using namespace edm;
32 using namespace std;
33 
34 
35 
36 DTPreCalibrationTask::DTPreCalibrationTask(const edm::ParameterSet& ps){
37 
38  LogTrace("DTPreCalibSummary") <<"[DTPrecalibrationTask]: Constructor"<<endl;
39 
40  dbe = Service<DQMStore>().operator->();
41 
42  // Label to retrieve DT digis from the event
43  digiLabel = ps.getUntrackedParameter<string>("digiLabel");
44 
45  // parameter for Time Boxes booking
46  minTriggerWidth = ps.getUntrackedParameter<int>("minTriggerWidth",2000);
47  maxTriggerWidth = ps.getUntrackedParameter<int>("maxTriggerWidth",6000);
48 
49  // histo saving on file
50  saveFile = ps.getUntrackedParameter<bool>("SaveFile",false);
51  // output file name
52  outputFileName = ps.getUntrackedParameter<string>("outputFileName");
53  // get the histo folder name
54  folderName = ps.getUntrackedParameter<string>("folderName");
55 
56 }
57 
58 
59 DTPreCalibrationTask::~DTPreCalibrationTask(){}
60 
61 
62 void DTPreCalibrationTask::beginJob(){
63 
64  for(int wheel=-2; wheel<=2; wheel++){
65  for(int sector=1; sector<=14; sector++){
66  LogTrace("DTPreCalibSummary") <<"[DTPrecalibrationTask]: Book histos for wheel "<<wheel<<", sector "<<sector<<endl;
67  dbe->setCurrentFolder(folderName+"/TimeBoxes");
68  bookTimeBoxes(wheel, sector);
69  dbe->setCurrentFolder(folderName+"/OccupancyHistos");
70  if(sector<13) bookOccupancyPlot(wheel, sector);
71  }
72  }
73 
74 }
75 
76 
77 void DTPreCalibrationTask::analyze(const edm::Event& event, const edm::EventSetup& setup) {
78 
79  // Get the digis from the event
80  edm::Handle<DTDigiCollection> dtdigis;
81  event.getByLabel(digiLabel, dtdigis);
82 
83  // LOOP OVER ALL THE DIGIS OF THE EVENT
84  DTDigiCollection::DigiRangeIterator dtLayerId_It;
85  for (dtLayerId_It=dtdigis->begin(); dtLayerId_It!=dtdigis->end(); ++dtLayerId_It){
86  for (DTDigiCollection::const_iterator digiIt = ((*dtLayerId_It).second).first;
87  digiIt!=((*dtLayerId_It).second).second; ++digiIt){
88 
89  //Fill the Time Boxes
90  int tdcTime = (*digiIt).countsTDC();
91  TimeBoxes[make_pair((*dtLayerId_It).first.superlayerId().chamberId().wheel(),
92  (*dtLayerId_It).first.superlayerId().chamberId().sector())]->Fill(tdcTime);
93 
94  //Fill the occupancy plot
95  const DTLayerId dtLId = (*dtLayerId_It).first;
96  int yBin = (dtLId.station()-1)*12+dtLId.layer()+4*(dtLId.superlayer()-1);
97  if(dtLId.station()==4 && dtLId.superlayer()==3)
98  yBin = (dtLId.station()-1)*12+dtLId.layer()+4*(dtLId.superlayer()-2);
99  if((*dtLayerId_It).first.superlayerId().chamberId().sector()<13)
100  OccupancyHistos[make_pair((*dtLayerId_It).first.superlayerId().chamberId().wheel(),
101  (*dtLayerId_It).first.superlayerId().chamberId().sector())]->Fill((*digiIt).wire(),yBin);
102  else{
103  if(dtLId.superlayer()!=3) yBin = 44 + dtLId.layer();
104  else yBin = 48 + dtLId.layer();
105  if((*dtLayerId_It).first.superlayerId().chamberId().sector()==13)
106  OccupancyHistos[make_pair((*dtLayerId_It).first.superlayerId().chamberId().wheel(),4)]->Fill((*digiIt).wire(),yBin);
107  if((*dtLayerId_It).first.superlayerId().chamberId().sector()==14)
108  OccupancyHistos[make_pair((*dtLayerId_It).first.superlayerId().chamberId().wheel(),10)]->Fill((*digiIt).wire(),yBin);
109  }
110 
111  }
112  }
113 
114 }
115 
116 
117 void DTPreCalibrationTask::endJob(){
118 
119  // save file for offLine analysis
120  if(saveFile)
121  dbe->save(outputFileName);
122 
123 }
124 
125 
126 void DTPreCalibrationTask::bookTimeBoxes(int wheel, int sector) {
127 
128  stringstream wh; wh << wheel;
129  stringstream sec; sec << sector;
130 
131  // book the time boxes
132  TimeBoxes[make_pair(wheel, sector)]= dbe->book1D("TimeBox_W"+wh.str()+"_Sec"+sec.str(), "Time Box W"+wh.str()+"_Sec"+sec.str(),(maxTriggerWidth-minTriggerWidth)/50, minTriggerWidth, maxTriggerWidth);
133  TimeBoxes[make_pair(wheel, sector)]->setAxisTitle("TDC counts");
134 
135 }
136 
137 
138 
139 void DTPreCalibrationTask::bookOccupancyPlot(int wheel, int sector) {
140 
141  stringstream wh; wh << wheel;
142  stringstream sec; sec << sector;
143 
144  // book the occpancy plot
145  if(sector==4 || sector==10)
146  OccupancyHistos[make_pair(wheel, sector)]= dbe->book2D("Occupancy_W"+wh.str()+"_Sec"+sec.str(), "Occupancy W"+wh.str()+"_Sec"+sec.str(),100,1,100,52,1,53);
147  else
148  OccupancyHistos[make_pair(wheel, sector)]= dbe->book2D("Occupancy_W"+wh.str()+"_Sec"+sec.str(), "Occupancy W"+wh.str()+"_Sec"+sec.str(),100,1,100,44,1,45);
149  OccupancyHistos[make_pair(wheel, sector)]->setAxisTitle("wire number", 1);
150  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(1,"M1L1",2);
151  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(2,"M1L2",2);
152  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(3,"M1L3",2);
153  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(4,"M1L4",2);
154  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(5,"M1L5",2);
155  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(6,"M1L6",2);
156  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(7,"M1L7",2);
157  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(8,"M1L8",2);
158  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(9,"M1L9",2);
159  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(10,"M1L10",2);
160  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(11,"M1L11",2);
161  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(12,"M1L12",2);
162  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(13,"M2L1",2);
163  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(14,"M2L2",2);
164  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(15,"M2L3",2);
165  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(16,"M2L4",2);
166  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(17,"M2L5",2);
167  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(18,"M2L6",2);
168  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(19,"M2L7",2);
169  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(20,"M2L8",2);
170  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(21,"M2L9",2);
171  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(22,"M2L10",2);
172  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(23,"M2L11",2);
173  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(24,"M2L12",2);
174  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(25,"M3L1",2);
175  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(26,"M3L2",2);
176  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(27,"M3L3",2);
177  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(28,"M3L4",2);
178  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(29,"M3L5",2);
179  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(30,"M3L6",2);
180  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(31,"M3L7",2);
181  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(32,"M3L8",2);
182  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(33,"M3L9",2);
183  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(34,"M3L10",2);
184  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(35,"M3L11",2);
185  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(36,"M3L12",2);
186  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(37,"M4L1",2);
187  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(38,"M4L2",2);
188  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(39,"M4L3",2);
189  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(40,"M4L4",2);
190  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(41,"M4L5",2);
191  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(42,"M4L6",2);
192  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(43,"M4L7",2);
193  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(44,"M4L8",2);
194  if(sector==4){
195  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(45,"M4Sec13L1",2);
196  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(46,"M4Sec13L2",2);
197  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(47,"M4Sec13L3",2);
198  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(48,"M4Sec13L4",2);
199  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(49,"M4Sec13L5",2);
200  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(50,"M4Sec13L6",2);
201  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(51,"M4Sec13L7",2);
202  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(52,"M4Sec13L8",2);
203  }
204  if(sector==10){
205  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(45,"M4Sec14L1",2);
206  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(46,"M4Sec14L2",2);
207  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(47,"M4Sec14L3",2);
208  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(48,"M4Sec14L4",2);
209  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(49,"M4Sec14L5",2);
210  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(50,"M4Sec14L6",2);
211  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(51,"M4Sec14L7",2);
212  OccupancyHistos[make_pair(wheel, sector)]->setBinLabel(52,"M4Sec14L8",2);
213  }
214 
215 }
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
DTPreCalibrationTask::DTPreCalibrationTask
DTPreCalibrationTask(const edm::ParameterSet &ps)
Constructor.
Definition: DTPreCalibrationTask.cc:36
cuy.yBin
tuple yBin
Definition: cuy.py:891
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:250
DTPreCalibrationTask::~DTPreCalibrationTask
virtual ~DTPreCalibrationTask()
Destructor.
Definition: DTPreCalibrationTask.cc:59
dumpDBToFile_GT_ttrig_cfg.outputFileName
tuple outputFileName
Definition: dumpDBToFile_GT_ttrig_cfg.py:31
HcalObjRepresent::Fill
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:611
DTPreCalibrationTask::bookTimeBoxes
void bookTimeBoxes(int wheel, int sector)
Book histos.
Definition: DTPreCalibrationTask.cc:126
event
How EventSelector::AcceptEvent() decides whether to accept an event for output otherwise it is excluding the probing of A single or multiple positive and the trigger will pass if any such matching triggers are PASS or EXCEPTION[A criterion thatmatches no triggers at all is detected and causes a throw.] A single negative with an expectation of appropriate bit checking in the decision and the trigger will pass if any such matching triggers are FAIL or EXCEPTION A wildcarded negative criterion that matches more than one trigger in the trigger but the state exists so we define the behavior If all triggers are the negative crieriion will lead to accepting the event(this again matches the behavior of"!*"before the partial wildcard feature was incorporated).The per-event"cost"of each negative criterion with multiple relevant triggers is about the same as!*was in the past
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:502
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:105
DTChamberId::station
int station() const
Return the station number.
Definition: DTChamberId.h:51
DTPreCalibrationTask::bookOccupancyPlot
void bookOccupancyPlot(int wheel, int sector)
Definition: DTPreCalibrationTask.cc:139
DTPreCalibrationTask::analyze
void analyze(const edm::Event &e, const edm::EventSetup &c)
Analyze.
Definition: DTPreCalibrationTask.cc:77
HcalObjRepresent::setup
void setup(std::vector< TH2F > &depth, std::string name, std::string units="")
Definition: HcalObjRepresent.h:628
dtTPAnalyzer_cfg.digiLabel
tuple digiLabel
Definition: dtTPAnalyzer_cfg.py:47
Generated for CMSSW Reference Manual by doxygen 1.8.5