DTRunConditionVar.cc

Go to the documentation of this file.

/******* \class DTRunConditionVar *******
 *
 * Description:
 *
 *  detailed description
 *
 * \author : Paolo Bellan, Antonio Branca
 * $date   : 23/09/2011 15:42:04 CET $
 *
 * Modification:
 *
 *********************************/

#include "DTRunConditionVar.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "Geometry/Records/interface/MuonGeometryRecord.h"

#include "DataFormats/MuonDetId/interface/MuonSubdetId.h"
#include "CondFormats/DataRecord/interface/DTMtimeRcd.h"

#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"

#include "MagneticField/Records/interface/IdealMagneticFieldRecord.h"

#include "RecoMuon/Navigation/interface/DirectMuonNavigation.h"

#include "DataFormats/DetId/interface/DetId.h"
#include "TrackingTools/DetLayers/interface/DetLayer.h"
#include "DataFormats/Common/interface/RefToBase.h"

#include <TMath.h>
#include <cmath>

using namespace std;
using namespace edm;

DTRunConditionVar::DTRunConditionVar(const ParameterSet& pSet):
  // Get the debug parameter for verbose output
  debug(pSet.getUntrackedParameter<bool>("debug",false)),
  nMinHitsPhi(pSet.getUntrackedParameter<int>("nMinHitsPhi")),
  maxAnglePhiSegm(pSet.getUntrackedParameter<double>("maxAnglePhiSegm")),
  dt4DSegmentsToken_(consumes<DTRecSegment4DCollection>(
      pSet.getParameter<InputTag>("recoSegments")))
{

}

DTRunConditionVar::~DTRunConditionVar()
{
  LogTrace("DTDQM|DTMonitorModule|DTRunConditionVar")
    << "DTRunConditionVar: destructor called";

  // free memory
}

void DTRunConditionVar::bookHistograms(DQMStore::IBooker & ibooker,
                                             edm::Run const & iRun,
                                             edm::EventSetup const & /* iSetup */) {

  LogTrace("DTDQM|DTMonitorModule|DTRunConditionVar")
    << "DTRunConditionVar: bookHistograms";

  for(int wheel=-2;wheel<=2;wheel++){
    for(int sec=1; sec<=14; sec++) {
      for(int stat=1; stat<=4; stat++) {

        bookChamberHistos(ibooker,DTChamberId(wheel,stat,sec),"VDrift_FromSegm",100,0.0043,0.0065);
        bookChamberHistos(ibooker,DTChamberId(wheel,stat,sec),"T0_FromSegm",100,-25.,25.);

      }
    }
  }


  return;
}

void DTRunConditionVar::dqmBeginRun(const Run& run, const EventSetup& setup)
{
  // Get the DT Geometry
  setup.get<MuonGeometryRecord>().get(dtGeom);

  return;
}

void DTRunConditionVar::analyze(const Event & event,
    const EventSetup& eventSetup)
{

  LogTrace("DTDQM|DTMonitorModule|DTRunConditionVar") <<
    "--- [DTRunConditionVar] Event analysed #Run: " <<
    event.id().run() << " #Event: " << event.id().event() << endl;

  // Get the DT Geometry
  ESHandle<DTGeometry> dtGeom;
  eventSetup.get<MuonGeometryRecord>().get(dtGeom);

  // Get the map of vdrift from the setup
  eventSetup.get<DTMtimeRcd>().get(mTime);
  mTimeMap_ = &*mTime;

  // Get the segment collection from the event
  Handle<DTRecSegment4DCollection> all4DSegments;
  event.getByToken(dt4DSegmentsToken_, all4DSegments);

  // Loop over the segments
  for(DTRecSegment4DCollection::const_iterator segment  = all4DSegments->begin();
      segment != all4DSegments->end(); ++segment){

    // Get the chamber from the setup
    DTChamberId DTid = (DTChamberId) segment->chamberId();
    uint32_t indexCh = DTid.rawId();

    // Fill v-drift values
    if( (*segment).hasPhi() ) {

      int nHitsPhi = (*segment).phiSegment()->degreesOfFreedom()+2;
      double xdir = (*segment).phiSegment()->localDirection().x();
      double zdir = (*segment).phiSegment()->localDirection().z();

      double anglePhiSegm = fabs(atan(xdir/zdir))*180./TMath::Pi();

      if( nHitsPhi >= nMinHitsPhi && anglePhiSegm <= maxAnglePhiSegm ) {

        double segmentVDrift = segment->phiSegment()->vDrift();


        DTSuperLayerId indexSLPhi1(DTid,1);
        DTSuperLayerId indexSLPhi2(DTid,3);

        float vDriftPhi1(0.), vDriftPhi2(0.);
        float ResPhi1(0.), ResPhi2(0.);
        int status1 = mTimeMap_->get(indexSLPhi1,vDriftPhi1,ResPhi1,DTVelocityUnits::cm_per_ns);
        int status2 = mTimeMap_->get(indexSLPhi2,vDriftPhi2,ResPhi2,DTVelocityUnits::cm_per_ns);

        if(status1 != 0 || status2 != 0) {
          DTSuperLayerId sl = (status1 != 0) ? indexSLPhi1 : indexSLPhi2;
          throw cms::Exception("DTRunConditionVarClient") << "Could not find vDrift entry in DB for"
            << sl << endl;
        }

        float vDriftMed = (vDriftPhi1 + vDriftPhi2) / 2.;

        segmentVDrift = vDriftMed*(1. - segmentVDrift);

        double segmentT0 = segment->phiSegment()->t0();

        if( segment->phiSegment()->ist0Valid() ) (chamberHistos[indexCh])["T0_FromSegm"]->Fill(segmentT0);
        if( segmentVDrift != vDriftMed ) (chamberHistos[indexCh])["VDrift_FromSegm"]->Fill(segmentVDrift);

      }
    }

  } //end loop on segment

} //end analyze

void DTRunConditionVar::bookChamberHistos(DQMStore::IBooker & ibooker,const DTChamberId& dtCh, string histoType, int nbins, float min, float max) {

  int wh = dtCh.wheel();
  int sc = dtCh.sector();
  int st = dtCh.station();
  stringstream wheel; wheel << wh;
  stringstream station; station << st;
  stringstream sector; sector << sc;

  string bookingFolder = "DT/02-Segments/Wheel" + wheel.str() + "/Sector" + sector.str() + "/Station" + station.str();
  string histoTag      = "_W" + wheel.str() + "_Sec" + sector.str() + "_St" + station.str();

  ibooker.setCurrentFolder(bookingFolder);

  LogTrace ("DTDQM|DTMonitorModule|DTRunConditionVar")
    << "[DTRunConditionVar]: booking histos in " << bookingFolder << endl;

  string histoName = histoType  +  histoTag;
  string histoLabel = histoType;

  (chamberHistos[dtCh.rawId()])[histoType] =
    ibooker.book1D(histoName,histoLabel,nbins,min,max);

}

// Local Variables:
// show-trailing-whitespace: t
// truncate-lines: t
// End: