#include <DTTTrigResidualCorrection.h>

Inheritance diagram for DTTTrigResidualCorrection:

Public Member Functions
virtual DTTTrigData	correction (const DTSuperLayerId &)
	DTTTrigResidualCorrection (const edm::ParameterSet &)
virtual void	setES (const edm::EventSetup &setup)
virtual	~DTTTrigResidualCorrection ()
Private Member Functions
const TH1F *	getHisto (const DTSuperLayerId &)
std::string	getHistoName (const DTSuperLayerId &slID)
Private Attributes
std::string	dbLabel_
const DTMtime *	mTimeMap_
std::string	rootBaseDir_
TFile *	rootFile_
const DTTtrig *	tTrigMap_
bool	useFit_
bool	useSlopesCalib_
double	vDriftEff_ [5][14][4][3]

Detailed Description

Concrete implementation of a DTTTrigBaseCorrection. Computes residual correction for tTrig

Revision:: 1.4

Author:: A. Vilela Pereira

Definition at line 26 of file DTTTrigResidualCorrection.h.

Constructor & Destructor Documentation

DTTTrigResidualCorrection::DTTTrigResidualCorrection ( const edm::ParameterSet & pset )

Definition at line 40 of file DTTTrigResidualCorrection.cc.

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

                                                                             {
  string residualsRootFile = pset.getParameter<string>("residualsRootFile");
  rootFile_ = new TFile(residualsRootFile.c_str(),"READ");
  rootBaseDir_ = pset.getUntrackedParameter<string>("rootBaseDir","/DQMData/DT/DTCalibValidation"); 
  useFit_ = pset.getParameter<bool>("useFitToResiduals");
  //useConstantvDrift_ = pset.getParameter<bool>("useConstantDriftVelocity");
  dbLabel_  = pset.getUntrackedParameter<string>("dbLabel", "");
  useSlopesCalib_ = pset.getUntrackedParameter<bool>("useSlopesCalib",false);

  // Load external slopes
  if(useSlopesCalib_){ 
     ifstream fileSlopes; 
     fileSlopes.open( pset.getParameter<FileInPath>("slopesFileName").fullPath().c_str() );
     //fin.open("slopesDB.txt");

     int tmp_wheel = 0;
     int tmp_sector = 0;
     int tmp_station = 0;
     int tmp_SL = 0;
     double tmp_ttrig = 0.;
     double tmp_t0 = 0.;
     double tmp_kfact = 0.;
     int tmp_a = 0;
     int tmp_b = 0;
     int tmp_c = 0;
     int tmp_d = 0;
     double tmp_v_eff = 0.;

     while(!fileSlopes.eof()){

        fileSlopes >> tmp_wheel >> tmp_sector >> tmp_station >> tmp_SL  >> tmp_a >> tmp_b >>
                      tmp_ttrig >> tmp_t0 >> tmp_kfact >> tmp_c >> tmp_d >> tmp_v_eff;

        vDriftEff_[tmp_wheel+2][tmp_sector-1][tmp_station-1][tmp_SL-1] = -tmp_v_eff;

     }

     fileSlopes.close();
  } 

}

DTTTrigResidualCorrection::~DTTTrigResidualCorrection ( ) [virtual]

Definition at line 82 of file DTTTrigResidualCorrection.cc.

                                                      {
  delete rootFile_;
}

Member Function Documentation

DTTTrigData DTTTrigResidualCorrection::correction ( const DTSuperLayerId & slId ) [virtual]

Implements DTTTrigBaseCorrection.

Definition at line 100 of file DTTTrigResidualCorrection.cc.

References DTVelocityUnits::cm_per_ns, Exception, cropTnPTrees::frac, getHisto(), LogTrace, plotscripts::mean(), DTTimeUnits::ns, DTChamberId::sector(), relativeConstraints::station, DTChamberId::station(), ntuplemaker::status, DTSuperLayerId::superLayer(), DTChamberId::wheel(), and ExpressReco_HICollisions_FallBack::x.

                                                                            {

  float tTrigMean,tTrigSigma,kFactor;
  int status = tTrigMap_->get(slId,tTrigMean,tTrigSigma,kFactor,DTTimeUnits::ns);
  if(status != 0) throw cms::Exception("[DTTTrigResidualCorrection]") << "Could not find tTrig entry in DB for"
                                                                      << slId << endl;

  float vDrift,hitResolution;
  status = mTimeMap_->get(slId,vDrift,hitResolution,DTVelocityUnits::cm_per_ns);
  if(status != 0) throw cms::Exception("[DTTTrigResidualCorrection]") << "Could not find vDrift entry in DB for"
                                                                      << slId << endl;
  TH1F residualHisto = *(getHisto(slId));
  double fitMean = -1.;

  if(useFit_){

    RooRealVar x("x","residual",-1.,1.);
    RooRealVar mean("mean","mean",residualHisto.GetMean(),-0.3,0.3);
    RooRealVar sigma1("sigma1","sigma1",0.,0.5);
    RooRealVar sigma2("sigma2","sigma2",0.,0.5);

    RooRealVar frac("frac","frac",0.,1.);

    RooGaussian myg1("myg1","Gaussian distribution",x,mean,sigma1);
    RooGaussian myg2("myg2","Gaussian distribution",x,mean,sigma2);

    RooAddPdf myg("myg","myg",RooArgList(myg1,myg2),RooArgList(frac));

    RooDataHist hdata("hdata","Binned data",RooArgList(x),&residualHisto);
    myg.fitTo(hdata,RooFit::Minos(0),RooFit::Range(-0.2,0.2));

    fitMean = mean.getVal();
    LogTrace("Calibration") << "[DTTTrigResidualCorrection]: \n"
                            << "                             Mean, Fit Mean    = " << residualHisto.GetMean()
                            << ", " << fitMean << "\n"
                            << "                             RMS, Fit RMS      = " << residualHisto.GetRMS()
                            << ", " << sigma1.getVal();

    //static int count = 0;
    /*
    if(count == 0){
      RooPlot *xframe = x.frame();
      hdata.plotOn(xframe);
      myg.plotOn(xframe);
      TCanvas c1;
      c1.cd();
      xframe->Draw();
      c1.SaveAs("prova.eps");
      count++;
    }
    */

  }

  int wheel = slId.wheel();
  int sector = slId.sector();
  int station = slId.station();
  int superLayer = slId.superLayer();

  double resTime = 0.;
  if(useSlopesCalib_){
     double vdrift_eff = vDriftEff_[wheel+2][sector-1][station-1][superLayer-1];
     if(vdrift_eff == 0) vdrift_eff = vDrift;

     if(vdrift_eff) resTime = ( (useFit_) ? fitMean : residualHisto.GetMean() )/vdrift_eff;

     LogTrace("Calibration") << "[DTTTrigResidualCorrection]: Effective vDrift, correction to tTrig = "
                             << vdrift_eff << ", " << resTime;
  } else{
     if(vDrift) resTime = ( (useFit_) ? fitMean : residualHisto.GetMean() )/vDrift;

     LogTrace("Calibration") << "[DTTTrigResidualCorrection]: vDrift from DB, correction to tTrig = "
                             << vDrift << ", " << resTime;
  }
 
  double corrMean = tTrigMean;
  double corrSigma = (tTrigSigma != 0.) ? tTrigSigma : 1.;
  double corrKFact = (tTrigSigma != 0.) ? (kFactor + resTime/tTrigSigma) : resTime;

  return DTTTrigData(corrMean,corrSigma,corrKFact);  
}

const TH1F * DTTTrigResidualCorrection::getHisto ( const DTSuperLayerId & slId ) [private]

Definition at line 182 of file DTTTrigResidualCorrection.cc.

References Exception, trackerHits::histo, and LogTrace.

                                                                          {
  string histoName = getHistoName(slId);
  LogTrace("Calibration") << "[DTTTrigResidualCorrection]: Accessing histogram " << histoName.c_str();
  TH1F* histo = static_cast<TH1F*>(rootFile_->Get(histoName.c_str()));
  if(!histo) throw cms::Exception("[DTTTrigResidualCorrection]") << "residual histogram not found:"
                                                                 << histoName << endl; 
  return histo;
}

string DTTTrigResidualCorrection::getHistoName ( const DTSuperLayerId & slID ) [private]

Definition at line 191 of file DTTTrigResidualCorrection.cc.

References DTChamberId::sector(), relativeConstraints::station, DTChamberId::station(), ExpressReco_HICollisions_FallBack::step, cmsPerfCommons::Step, DTSuperLayerId::superlayer(), and DTChamberId::wheel().

                                                                         {

  int step = 3;
  stringstream wheel; wheel << slId.wheel();    
  stringstream station; station << slId.station();      
  stringstream sector; sector << slId.sector(); 
  stringstream superLayer; superLayer << slId.superlayer();
  stringstream Step; Step << step;

  string histoName =
    rootBaseDir_ + "/Wheel" + wheel.str() + 
    "/Station" + station.str() +
    "/Sector" + sector.str() +
    "/hResDist_STEP" + Step.str() +
    "_W" + wheel.str() +
    "_St" + station.str() +
    "_Sec" + sector.str() +
    "_SL" + superLayer.str();

  return histoName;
}

void DTTTrigResidualCorrection::setES ( const edm::EventSetup & setup ) [virtual]

Implements DTTTrigBaseCorrection.

Definition at line 86 of file DTTTrigResidualCorrection.cc.

References edm::EventSetup::get().

                                                             {
  // Get tTrig record from DB
  ESHandle<DTTtrig> tTrig;
  //setup.get<DTTtrigRcd>().get(tTrig);
  setup.get<DTTtrigRcd>().get(dbLabel_,tTrig);
  tTrigMap_ = &*tTrig;

  // Get vDrift record
  ESHandle<DTMtime> mTimeHandle;
  setup.get<DTMtimeRcd>().get(mTimeHandle);
  mTimeMap_ = &*mTimeHandle;

}