99 LogTrace(
"DTDQM|DTMonitorModule|DTRunConditionVar") <<
100 "--- [DTRunConditionVar] Event analysed #Run: " <<
101 event.id().run() <<
" #Event: " <<
event.id().event() << endl;
117 segment != all4DSegments->end(); ++segment){
121 uint32_t indexCh = DTid.
rawId();
124 if( (*segment).hasPhi() ) {
126 int nHitsPhi = (*segment).phiSegment()->degreesOfFreedom()+2;
127 double xdir = (*segment).phiSegment()->localDirection().x();
128 double zdir = (*segment).phiSegment()->localDirection().z();
130 double anglePhiSegm = fabs(atan(xdir/zdir))*180./
TMath::Pi();
134 double segmentVDrift = segment->phiSegment()->vDrift();
140 float vDriftPhi1(0.), vDriftPhi2(0.);
141 float ResPhi1(0.), ResPhi2(0.);
145 if(status1 != 0 || status2 != 0) {
147 throw cms::Exception(
"DTRunConditionVarClient") <<
"Could not find vDrift entry in DB for" 151 float vDriftMed = (vDriftPhi1 + vDriftPhi2) / 2.;
153 segmentVDrift = vDriftMed*(1. - segmentVDrift);
155 double segmentT0 = segment->phiSegment()->t0();
157 if( segment->phiSegment()->ist0Valid() ) (
chamberHistos[indexCh])[
"T0_FromSegm"]->Fill(segmentT0);
158 if( segmentVDrift != vDriftMed ) (
chamberHistos[indexCh])[
"VDrift_FromSegm"]->
Fill(segmentVDrift);
edm::ESHandle< DTMtime > mTime
edm::EDGetTokenT< DTRecSegment4DCollection > dt4DSegmentsToken_
constexpr uint32_t rawId() const
get the raw id
const DTMtime * mTimeMap_
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
std::map< uint32_t, std::map< std::string, MonitorElement * > > chamberHistos
int get(int wheelId, int stationId, int sectorId, int slId, float &mTime, float &mTrms, DTTimeUnits::type unit) const
edm::ESHandle< DTGeometry > dtGeom