slope.cxx File Reference

Go to the source code of this file.

Functions
void	slope ()

Variables
const double	k_high = 0.0
const double	k_low = 0.0
const double	b = 1.0
const bool	b_trunc = true

Function Documentation

slope()

void slope

(

)

Definition at line 12 of file slope.cxx.

{
    TFile *f=new TFile("eangle.root");
    ostringstream strout;
    strout<<"ddgcalib";
    TTree *tree = f->Get(strout.str().c_str()); 
    double dedx[100],eangle[100], fitmean[100], fitmeanerr[100];
    double pos_fitmean[100], pos_fitmeanerr[100], neg_fitmean[100], neg_fitmeanerr[100];
    //    tree->SetBranchAddress("gain",dedx);
    if(b_trunc) tree->SetBranchAddress("truncgain",dedx);
    else tree->SetBranchAddress("gain",dedx);
    tree->SetBranchAddress("eangle",eangle);
    tree->SetBranchAddress("fitmean",fitmean);
    tree->SetBranchAddress("fitmeanerr",fitmeanerr);
    tree->SetBranchAddress("pos_fitmean",pos_fitmean);
    tree->SetBranchAddress("pos_fitmeanerr",pos_fitmeanerr);
    tree->SetBranchAddress("neg_fitmean",neg_fitmean);
    tree->SetBranchAddress("neg_fitmeanerr",neg_fitmeanerr);
 
    cout<<"Entries: "<<tree->GetEntries()<<endl; 
    tree->GetEntry(0);
 
    int entry[1] = {100};
    double denangle[100];
    // for(int i=0; i<100; i++) denangle[i] = (dedx[i])*(k*(eangle[i])+b); // correct by a slope
    // next we try sub-region correction with a slope 
    for(int i=0; i<100; i++){
      if(fabs(eangle[i])<=0.04) denangle[i] = (dedx[i])*(k_high*(eangle[i])+b);
      //      cout << "eangle: " << eangle[i] << endl;
      else denangle[i] = (dedx[i])*(k_low*(eangle[i])+b); // correct by a slope
    }
    TFile *f_save = new TFile("eangle_new.root","recreate");
    TTree* entra = new TTree("entracalib", "entracalib");
    entra->Branch("1denangle_entry", entry, "entry[1]/I");
    entra->Branch("1denangle", denangle, "denangle[100]/D");
    entra->Fill();
    entra->Write();
    f_save->Close();
    f->Close();
}

Referenced by MdcSegInfoAxialO::calcFromOrigin(), MdcSegInfoAxialO::calcFromOrigin(), MdcSegInfoSterO::calcStereo(), MdcSegInfoCsmc::calcStraight(), MdcSegInfoCsmc::calcStraight(), XtInteCalib::calib(), dEdxVsCosNP(), DedxCorrecSvc::EntaCorrec(), EsTimeAlg::execute(), MucRec2DRoad::Fit(), MucRec2DRoad::GetSimpleFitParams(), MucRec2DRoad::SimpleFit(), and MucRec2DRoad::SimpleFitNoCurrentGap().

Variable Documentation

b

const double b = 1.0

Definition at line 9 of file slope.cxx.

Referenced by AsymGauss(), MucRec2DRoad::AttachHit(), TofQCorrSvc::BQRaw1(), TofQCorrSvc::BQRaw2(), EmcRecGeoPlane::Build(), EvtbTosllVectorAmp::CalcAmp(), EvtbTosllAmp::CalcMaxProb(), EvtSemiLeptonicAmp::CalcMaxProb(), HoughStereo::cald(), HoughTrack::cald_layer(), TMLink::cDrift(), CFCir::CFCir(), Lpar::circle(), compare_zsfit(), MucGeoGap::CompareIntersection(), BesEmcEndGeometry::ComputeDimAndPos(), ExtBesEmcEndGeometry::ComputeDimAndPos(), db_mysql::Connect(), BesTofConstruction::ConstructBr1Tof(), BesTofConstruction::ConstructBr2Tof(), BesTofConstruction::ConstructEcTof(), DifNumber::correlation(), DifNumber::correlation(), TMLink::dDrift(), EvtBsquark::decay(), BesPolygon2D::DistancetoPrimitive(), TRunge::DisToWire(), TTrackManager::divide(), TTrackManager::divideByIp(), TMLink::drift(), EmcRecGeoPlane::EmcRecGeoPlane(), TrkHotListFull::endFoundRange(), Kinem< R >::eq(), TofQCorrSvc::EQRaw(), EvtPawt(), MagFieldReader::execute(), EvtSpinAmp::extcont(), MagneticFieldSvc::fieldVector(), MucGeoGeneral::FindIntersection(), MucGeoGeneral::FindIntersectionQuadLocal(), MucGeoGeneral::FindIntersectionSurface(), TConformalFinder0::findSegments(), MdcxFindTracks::findz_cyl(), MucRec2DRoad::Fit(), TCosmicFitter::fit(), TCosmicFitter::fit(), RealDBUtil::ConnectionProvider::getcal(), TRadSpline3::GetCoeff(), TRadSpline5::GetCoeff(), BesTView::GetDistancetoAxis(), BesMagneticField::GetFieldValue(), ExtBesMagneticField::GetFieldValue(), MucGeometron::GetIntersectionQuadPlane(), MucGeometron::GetIntersectionQuadPlaneLocal(), MucMagneticField::getMucField(), MucRec2DRoad::GetSimpleFitParams(), BesTofConstruction::getXYZ(), DedxCorrecSvc::I2D(), HadronSaturation::I2D(), HadronSaturation::I2D(), EvtSpinAmp::intcont(), ParticleIDBase::interpolation(), HoughTrack::Leastfit(), HoughZsFit::leastFit(), HoughZsFit::leastLine(), less_layer(), less_layer_3D(), MucRecLineFit::LineFit(), MucRecLineFit::LineFit(), main(), main(), HTrackParameter::minDistanceTwoHelix(), mwCalloc(), mwCalloc_(), Kinem< R >::neq(), rdbModel::ColCompare::operator()(), rdbModel::FieldValCompare::operator()(), MdcUtilitySvc::pointOnHelixPatPar(), HTrackParameter::positionTwoHelix(), MucRec2DRoad::Project(), MucRec3DRoad::Project(), MucRec3DRoad::ProjectNoCurrentGap(), MucGeoGap::ProjectToGap(), MucGeoGap::ProjectToGapQuadLocal(), MucGeoGap::ProjectToGapSurface(), MucRec3DRoad::ProjectWithSigma(), MucRecQuadFit::QuadFit(), FTTrack::r_phi2Fit(), FTTrack::r_phi3Fit(), FTTrack::r_phi4Fit(), FTTrack::r_phiReFit(), KalFitTrack::radius_numf(), MdcTSF::ReadLut(), MdcTF::ReadTFLut(), BesEmcConstruction::RotAngleFromNewZ(), DifVector::rotate(), Lpar::s(), TMSelector::select(), eformat::write::set(), RootInterface::setBranchAddress(), rdbModel::Manager::setBuilder(), BesEventHeader::SetEventHeaderGeneral(), XmlRpc::XmlRpcSource::setKeepOpen(), EvtGammaMatrix::sigmaLower(), EvtGammaMatrix::sigmaUpper(), slope(), small_layer(), EvtConExc::SoftPhoton_xs(), solveQuad(), TTrackManager::sortBanksByPt(), sortByArcLength(), SortByB(), SortByPt(), sortBySequentialLength(), SortByWireId(), sortByWireSerialNumber(), SortByX(), EvtPolInt::splint(), EvtSVVHelAmp::SVVHel(), ers::Core::to_string(), TrackStatus(), and MagneticFieldSvc::uniFieldVector().

b_trunc

const bool b_trunc = true

Definition at line 10 of file slope.cxx.

Referenced by slope().

k_high

const double k_high = 0.0

Definition at line 7 of file slope.cxx.

Referenced by slope().

k_low

const double k_low = 0.0

Definition at line 8 of file slope.cxx.

Referenced by slope().