DTTTrigOffsetCalibration Class Reference

#include <DTTTrigOffsetCalibration.h>

Inheritance diagram for DTTTrigOffsetCalibration:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &eventSetup)
void	beginRun (const edm::Run &run, const edm::EventSetup &setup)
void	bookHistos (DTChamberId)
	DTTTrigOffsetCalibration (const edm::ParameterSet &pset)
void	endJob ()
virtual	~DTTTrigOffsetCalibration ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

Private Types
typedef std::map< DTChamberId, std::vector< TH1F * > >	ChamberHistosMap

Private Attributes
bool	checkNoisyChannels_
std::string	dbLabel
bool	doTTrigCorrection_
std::string	theCalibChamber_
TFile *	theFile_
double	theMaxChi2_
double	theMaxPhiAngle_
double	theMaxZAngle_
edm::InputTag	theRecHits4DLabel_
ChamberHistosMap	theT0SegHistoMap_
const DTTtrig *	tTrigMap

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

No description available.

Date:

2010/02/16 10:03:23

Revision:

1.2

Author

A. Vilela Pereira

Definition at line 28 of file DTTTrigOffsetCalibration.h.

Member Typedef Documentation

typedef std::map<DTChamberId, std::vector<TH1F*> > DTTTrigOffsetCalibration::ChamberHistosMap

private

Definition at line 45 of file DTTTrigOffsetCalibration.h.

Constructor & Destructor Documentation

DTTTrigOffsetCalibration::DTTTrigOffsetCalibration

(

const edm::ParameterSet &

pset

)

Definition at line 41 of file DTTTrigOffsetCalibration.cc.

References edm::ParameterSet::getParameter(), and edm::ParameterSet::getUntrackedParameter().

                                                                           {

  LogVerbatim("Calibration") << "[DTTTrigOffsetCalibration] Constructor called!";

  // the root file which will contain the histos
  string rootFileName = pset.getUntrackedParameter<string>("rootFileName","DTT0SegHistos.root");
  theFile_ = new TFile(rootFileName.c_str(), "RECREATE");
  theFile_->cd();

  dbLabel  = pset.getUntrackedParameter<string>("dbLabel", "");

  // Do t0-seg correction to ttrig
  doTTrigCorrection_ = pset.getUntrackedParameter<bool>("doT0SegCorrection", false);

  // Chamber/s to calibrate
  theCalibChamber_ =  pset.getUntrackedParameter<string>("calibChamber", "All");

  // the name of the 4D rec hits collection
  theRecHits4DLabel_ = pset.getParameter<InputTag>("recHits4DLabel");

  //get the switch to check the noisy channels
  checkNoisyChannels_ = pset.getParameter<bool>("checkNoisyChannels");

  // get maximum chi2 value 
  theMaxChi2_ =  pset.getParameter<double>("maxChi2");

  // Maximum incidence angle for Phi SL 
  theMaxPhiAngle_ =  pset.getParameter<double>("maxAnglePhi");

  // Maximum incidence angle for Theta SL 
  theMaxZAngle_ =  pset.getParameter<double>("maxAngleZ");
}

DTTTrigOffsetCalibration::~DTTTrigOffsetCalibration ( )

virtual

Definition at line 83 of file DTTTrigOffsetCalibration.cc.

                                                   {
  theFile_->Close();
  LogVerbatim("Calibration") << "[DTTTrigOffsetCalibration] Destructor called!";
}

Member Function Documentation

void DTTTrigOffsetCalibration::analyze ( const edm::Event &  event, const edm::EventSetup &  eventSetup  )

virtual

Implements edm::EDAnalyzer.

Definition at line 88 of file DTTTrigOffsetCalibration.cc.

References bookHistos(), DTChamberId, edm::EventSetup::get(), DTRecSegment2D::localDirection(), DTRecSegment2D::localPosition(), LogTrace, Geom::pi(), DTRecSegment2D::specificRecHits(), crabStatusFromReport::statusMap, DTChamber::superLayer(), DTSLRecSegment2D::superLayerId(), DTLayerId::superlayerId(), GeomDet::toGlobal(), GeomDet::toLocal(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                        {
  theFile_->cd();
  DTChamberId chosenChamberId;

  if(theCalibChamber_ != "All") {
    stringstream linestr;
    int selWheel, selStation, selSector;
    linestr << theCalibChamber_;
    linestr >> selWheel >> selStation >> selSector;
    chosenChamberId = DTChamberId(selWheel, selStation, selSector);
    LogVerbatim("Calibration") << " chosen chamber " << chosenChamberId << endl;
  }

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

  // Get the rechit collection from the event
  Handle<DTRecSegment4DCollection> all4DSegments;
  event.getByLabel(theRecHits4DLabel_, all4DSegments); 

  // Get the map of noisy channels
  ESHandle<DTStatusFlag> statusMap;
  if(checkNoisyChannels_) eventSetup.get<DTStatusFlagRcd>().get(statusMap);

  // Loop over segments by chamber
  DTRecSegment4DCollection::id_iterator chamberIdIt;
  for (chamberIdIt = all4DSegments->id_begin();
       chamberIdIt != all4DSegments->id_end();
       ++chamberIdIt){

    // Get the chamber from the setup
    const DTChamber* chamber = dtGeom->chamber(*chamberIdIt);
    LogTrace("Calibration") << "Chamber Id: " << *chamberIdIt;

    // Book histos
    if(theT0SegHistoMap_.find(*chamberIdIt) == theT0SegHistoMap_.end()){
      bookHistos(*chamberIdIt);
    }
   
    // Calibrate just the chosen chamber/s    
    if((theCalibChamber_ != "All") && ((*chamberIdIt) != chosenChamberId)) continue;

    // Get the range for the corresponding ChamberId
    DTRecSegment4DCollection::range range = all4DSegments->get((*chamberIdIt));

    // Loop over the rechits of this DetUnit
    for (DTRecSegment4DCollection::const_iterator segment = range.first;
         segment!=range.second; ++segment){

      LogTrace("Calibration") << "Segment local pos (in chamber RF): " << (*segment).localPosition()
                              << "\nSegment global pos: " << chamber->toGlobal((*segment).localPosition());
      
      //get the segment chi2
      double chiSquare = ((*segment).chi2()/(*segment).degreesOfFreedom());
      // cut on the segment chi2 
      if(chiSquare > theMaxChi2_) continue;

      // get the Phi 2D segment and plot the angle in the chamber RF
      if(!((*segment).phiSegment())){
        LogTrace("Calibration") << "No phi segment";
      }
      LocalPoint phiSeg2DPosInCham;  
      LocalVector phiSeg2DDirInCham;

      bool segmNoisy = false;
      map<DTSuperLayerId,vector<DTRecHit1D> > hitsBySLMap; 

      if((*segment).hasPhi()){
        const DTChamberRecSegment2D* phiSeg = (*segment).phiSegment();  // phiSeg lives in the chamber RF
        phiSeg2DPosInCham = phiSeg->localPosition();  
        phiSeg2DDirInCham = phiSeg->localDirection();

        vector<DTRecHit1D> phiHits = phiSeg->specificRecHits();
        for(vector<DTRecHit1D>::const_iterator hit = phiHits.begin();
            hit != phiHits.end(); ++hit) {
          DTWireId wireId = (*hit).wireId();
          DTSuperLayerId slId =  wireId.superlayerId();
          hitsBySLMap[slId].push_back(*hit); 

          // Check for noisy channels to skip them
          if(checkNoisyChannels_) {
            bool isNoisy = false;
            bool isFEMasked = false;
            bool isTDCMasked = false;
            bool isTrigMask = false;
            bool isDead = false;
            bool isNohv = false;
            statusMap->cellStatus(wireId, isNoisy, isFEMasked, isTDCMasked, isTrigMask, isDead, isNohv);
            if(isNoisy) {
              LogTrace("Calibration") << "Wire: " << wireId << " is noisy, skipping!";
              segmNoisy = true;
            }      
          }
        }
      }

      // get the Theta 2D segment and plot the angle in the chamber RF
      LocalVector zSeg2DDirInCham;
      LocalPoint zSeg2DPosInCham;
      if((*segment).hasZed()) {
        const DTSLRecSegment2D* zSeg = (*segment).zSegment();  // zSeg lives in the SL RF
        const DTSuperLayer* sl = chamber->superLayer(zSeg->superLayerId());
        zSeg2DPosInCham = chamber->toLocal(sl->toGlobal((*zSeg).localPosition())); 
        zSeg2DDirInCham = chamber->toLocal(sl->toGlobal((*zSeg).localDirection()));
        hitsBySLMap[zSeg->superLayerId()] = zSeg->specificRecHits();

        // Check for noisy channels to skip them
        vector<DTRecHit1D> zHits = zSeg->specificRecHits();
        for(vector<DTRecHit1D>::const_iterator hit = zHits.begin();
            hit != zHits.end(); ++hit) {
          DTWireId wireId = (*hit).wireId();
          if(checkNoisyChannels_) {
            bool isNoisy = false;
            bool isFEMasked = false;
            bool isTDCMasked = false;
            bool isTrigMask = false;
            bool isDead = false;
            bool isNohv = false;
            statusMap->cellStatus(wireId, isNoisy, isFEMasked, isTDCMasked, isTrigMask, isDead, isNohv);
            if(isNoisy) {
              LogTrace("Calibration") << "Wire: " << wireId << " is noisy, skipping!";
              segmNoisy = true;
            }      
          }
        }
      } 

      if (segmNoisy) continue;

      LocalPoint segment4DLocalPos = (*segment).localPosition();
      LocalVector segment4DLocalDir = (*segment).localDirection();
      if(fabs(atan(segment4DLocalDir.y()/segment4DLocalDir.z())* 180./Geom::pi()) > theMaxZAngle_) continue; // cut on the angle
      if(fabs(atan(segment4DLocalDir.x()/segment4DLocalDir.z())* 180./Geom::pi()) > theMaxPhiAngle_) continue; // cut on the angle
      // Fill t0-seg values
      if((*segment).hasPhi()) {
    //if((segment->phiSegment()->t0()) != 0.00){
        if(segment->phiSegment()->ist0Valid()){
      (theT0SegHistoMap_[*chamberIdIt])[0]->Fill(segment->phiSegment()->t0());
    }
      }
      if((*segment).hasZed()){
        //if((segment->zSegment()->t0()) != 0.00){
        if(segment->zSegment()->ist0Valid()){
      (theT0SegHistoMap_[*chamberIdIt])[1]->Fill(segment->zSegment()->t0());
    }
      }
      
      // Fill t0-seg values
    //      if((*segment).hasPhi()) (theT0SegHistoMap_[*chamberIdIt])[0]->Fill(segment->phiSegment()->t0());
    //  if((*segment).hasZed()) (theT0SegHistoMap_[*chamberIdIt])[1]->Fill(segment->zSegment()->t0());
      //if((*segment).hasZed() && (*segment).hasPhi()) {}

      /*//loop over the segments 
      for(map<DTSuperLayerId,vector<DTRecHit1D> >::const_iterator slIdAndHits = hitsBySLMap.begin(); slIdAndHits != hitsBySLMap.end();  ++slIdAndHits) {
        if (slIdAndHits->second.size() < 3) continue;
        DTSuperLayerId slId =  slIdAndHits->first;

      }*/
    }
  }
}

void DTTTrigOffsetCalibration::beginRun ( const edm::Run &  run, const edm::EventSetup &  setup  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 74 of file DTTTrigOffsetCalibration.cc.

References edm::EventSetup::get().

                                                                                     {
  if(doTTrigCorrection_){
    ESHandle<DTTtrig> tTrig;
    setup.get<DTTtrigRcd>().get(dbLabel,tTrig);
    tTrigMap = &*tTrig;
    LogVerbatim("Calibration") << "[DTTTrigOffsetCalibration]: TTrig version: " << tTrig->version() << endl; 
  }
}

void DTTTrigOffsetCalibration::bookHistos

(

DTChamberId

chId

)

Definition at line 298 of file DTTTrigOffsetCalibration.cc.

References combine::histos, LogTrace, DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, and DTChamberId::wheel().

                                                          {

  LogTrace("Calibration") << "   Booking histos for Chamber: " << chId;

  // Compose the chamber name
  stringstream wheel; wheel << chId.wheel();
  stringstream station; station << chId.station();
  stringstream sector; sector << chId.sector();

  string chHistoName =
    "_W" + wheel.str() +
    "_St" + station.str() +
    "_Sec" + sector.str();

  /*// Define the step
  stringstream Step; Step << step;

  string chHistoName =
    "_STEP" + Step.str() +
    "_W" + wheel.str() +
    "_St" + station.str() +
    "_Sec" + sector.str();

  theDbe->setCurrentFolder("DT/DTCalibValidation/Wheel" + wheel.str() +
                           "/Station" + station.str() +
                           "/Sector" + sector.str());
  // Create the monitor elements
  vector<MonitorElement *> histos;
  // Note hte order matters
  histos.push_back(theDbe->book1D("hRPhiSegT0"+chHistoName, "t0 from Phi segments", 200, -25., 25.));
  histos.push_back(theDbe->book1D("hRZSegT0"+chHistoName, "t0 from Z segments", 200, -25., 25.));*/

  vector<TH1F*> histos;
  // Note the order matters
  histos.push_back(new TH1F(("hRPhiSegT0"+chHistoName).c_str(), "t0 from Phi segments", 250, -60., 60.));
  if(chId.station() != 4) histos.push_back(new TH1F(("hRZSegT0"+chHistoName).c_str(), "t0 from Z segments", 250, -60., 60.));

  theT0SegHistoMap_[chId] = histos;
}

void DTTTrigOffsetCalibration::endJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 251 of file DTTTrigOffsetCalibration.cc.

References DTSuperLayerId, DTTimeUnits::ns, DTTtrig::set(), DTChamberId::station(), and DTCalibDBUtils::writeToDB().

                                      {
  theFile_->cd();
  
  LogVerbatim("Calibration") << "[DTTTrigOffsetCalibration]Writing histos to file!" << endl;

  for(ChamberHistosMap::const_iterator itChHistos = theT0SegHistoMap_.begin(); itChHistos != theT0SegHistoMap_.end(); ++itChHistos){
    for(vector<TH1F*>::const_iterator itHist = (*itChHistos).second.begin();
                                      itHist != (*itChHistos).second.end(); ++itHist) (*itHist)->Write();
  }

  if(doTTrigCorrection_){
    // Create the object to be written to DB
    DTTtrig* tTrig = new DTTtrig();

    for(ChamberHistosMap::const_iterator itChHistos = theT0SegHistoMap_.begin(); itChHistos != theT0SegHistoMap_.end(); ++itChHistos){
      DTChamberId chId = itChHistos->first;
      // Get SuperLayerId's for each ChamberId
      vector<DTSuperLayerId> slIds;
      slIds.push_back(DTSuperLayerId(chId,1));
      slIds.push_back(DTSuperLayerId(chId,3));
      if(chId.station() != 4) slIds.push_back(DTSuperLayerId(chId,2));

      for(vector<DTSuperLayerId>::const_iterator itSl = slIds.begin(); itSl != slIds.end(); ++itSl){      
        // Get old values from DB
        float ttrigMean = 0;
        float ttrigSigma = 0;
    float kFactor = 0;
        tTrigMap->get(*itSl,ttrigMean,ttrigSigma,kFactor,DTTimeUnits::ns);
        //FIXME: verify if values make sense
        // Set new values
        float ttrigMeanNew = ttrigMean;
        float ttrigSigmaNew = ttrigSigma;
    float t0SegMean = (itSl->superLayer() != 2)?itChHistos->second[0]->GetMean():itChHistos->second[1]->GetMean();

    float kFactorNew = (kFactor*ttrigSigma+t0SegMean)/ttrigSigma;

        tTrig->set(*itSl,ttrigMeanNew,ttrigSigmaNew,kFactorNew,DTTimeUnits::ns);
      }
    }
    LogVerbatim("Calibration")<< "[DTTTrigOffsetCalibration]Writing ttrig object to DB!" << endl;
    // Write the object to DB
    string tTrigRecord = "DTTtrigRcd";
    DTCalibDBUtils::writeToDB(tTrigRecord, tTrig);
  } 
}

Member Data Documentation

bool DTTTrigOffsetCalibration::checkNoisyChannels_

private

Definition at line 60 of file DTTTrigOffsetCalibration.h.

std::string DTTTrigOffsetCalibration::dbLabel

private

Definition at line 77 of file DTTTrigOffsetCalibration.h.

bool DTTTrigOffsetCalibration::doTTrigCorrection_

private

Definition at line 54 of file DTTTrigOffsetCalibration.h.

std::string DTTTrigOffsetCalibration::theCalibChamber_

private

Definition at line 75 of file DTTTrigOffsetCalibration.h.

TFile* DTTTrigOffsetCalibration::theFile_

private

Definition at line 51 of file DTTTrigOffsetCalibration.h.

double DTTTrigOffsetCalibration::theMaxChi2_

private

Definition at line 66 of file DTTTrigOffsetCalibration.h.

double DTTTrigOffsetCalibration::theMaxPhiAngle_

private

Definition at line 69 of file DTTTrigOffsetCalibration.h.

double DTTTrigOffsetCalibration::theMaxZAngle_

private

Definition at line 72 of file DTTTrigOffsetCalibration.h.

edm::InputTag DTTTrigOffsetCalibration::theRecHits4DLabel_

private

Definition at line 48 of file DTTTrigOffsetCalibration.h.

ChamberHistosMap DTTTrigOffsetCalibration::theT0SegHistoMap_

private

Definition at line 63 of file DTTTrigOffsetCalibration.h.

const DTTtrig* DTTTrigOffsetCalibration::tTrigMap

private

Definition at line 57 of file DTTTrigOffsetCalibration.h.