draw_charge_space_coarse_II.cxx File Reference

#include <fstream>

Go to the source code of this file.

Functions
TString	rootfilelist ("rootfile_list.txt")
double	f_real_work ()
void	draw_charge_space_coarse_II ()

Variables
bool	m_debug = false
const double	m_err_standard = 0.005
const double	m_diff_standard = 0.02
int	m_num_in_bin = 100
int	m_num_hists = 0
double	x [1000]
double	xerr [1000]
double	y [1000]
double	yerr [1000]

Function Documentation

draw_charge_space_coarse_II()

void draw_charge_space_coarse_II

(

)

Definition at line 63 of file draw_charge_space_coarse_II.cxx.

                                  {
  double m_error = f_real_work();
  if(m_error>m_err_standard){
    m_num_hists = 0;
    m_num_in_bin = (int)pow(m_error/m_err_standard,2)*m_num_in_bin;
    cout << "new num of bin: " << m_num_in_bin << endl;
    m_error = f_real_work();
  }
  cout << "average error: " << m_error << endl;
 
  TGraphErrors *g = new TGraphErrors(m_num_hists, x, y ,xerr, yerr);
  TCanvas *c1 = new TCanvas("c1", "c1", 10, 10, 800, 642);
  g->Draw("AP");
  for(int k=0; k<m_num_hists; k++){
    if(fabs(y[k+1]-y[k])/y[k]>m_diff_standard){
      TLine *l = new TLine(x[k+1]-xerr[k+1], y[k+1]-5*(yerr[k+1]), x[k+1]-xerr[k+1], y[k+1]+5*(yerr[k+1]));
      l->SetLineColor(2);
      l->Draw();
    }
  }
  c1->SaveAs("figs/charge_space_via_run_coarse_II.eps");
  c1->SaveAs("figs/charge_space_via_run_coarse_II.root");
 
}

f_real_work()

double f_real_work

(

)

Definition at line 10 of file draw_charge_space_coarse_II.cxx.

                    {
  // read input file
  ifstream input_file(rootfilelist);
  if(!input_file){
    cout << "cannot read input file correctly!" << endl;
    exit();
  }
  TString rootfilename("begin");
  TH1F *hx[1000];
  TH1F *h[1000];
  for(int j=0; j<1000; j++){
    h[j] = new TH1F(Form("h%d",j), "h", 3000, 0, 1000);
    hx[j] = new TH1F(Form("hx%d",j), "hx", 100000, 0, 100000);
  }
 
  while(rootfilename!=""){
    input_file >> rootfilename;
    if(m_debug) cout << rootfilename << endl;
    else if(m_num_hists%50==0){
      cout << rootfilename << endl;
      cout << "num of hists " << m_num_hists+1 << endl;
    }
    if(rootfilename=="") continue;
    TFile *f = new TFile(rootfilename);
    if(!f) continue;
    TTree *t = f->Get("n103");
    float runno(0);
    float dedx(0), costheta(0);
    t->SetBranchAddress("runNO", &runno);
    t->SetBranchAddress("dEdx_meas", &dedx);
    t->SetBranchAddress("costheta", &costheta);
    for(int i=0; i<t->GetEntries(); i++){
      t->GetEntry(i);
      if(i==0) hx[m_num_hists]->Fill(runno);
      if(fabs(costheta)<0.15) h[m_num_hists]->Fill(dedx);
    }
    if(h[m_num_hists]->GetEntries()>m_num_in_bin) m_num_hists ++;
    delete t;
    delete f;
  }
  double m_ave_err(0);
  for(int k=0; k<m_num_hists; k++){
    x[k] = hx[k]->GetMean();
    xerr[k] = hx[k]->GetRMS();
    y[k] = h[k]->GetMean();
    yerr[k] = h[k]->GetMeanError();
    m_ave_err += yerr[k]/y[k];
  }
  m_ave_err = m_ave_err/m_num_hists;
  return m_ave_err;
}

Referenced by draw_charge_space_coarse_II().

rootfilelist()

TString rootfilelist

(

"rootfile_list.txt"

)

Referenced by f_real_work().

Variable Documentation

m_debug

bool m_debug = false

Definition at line 3 of file draw_charge_space_coarse_II.cxx.

Referenced by MdcSeg::addHits(), ChargedDReconstruction::ChargedDReconstruction(), DedxCorrecSvc::DedxCorrecSvc(), DQADtag::DQADtag(), DsReconstruction::DsReconstruction(), EsTimeAlg::EsTimeAlg(), f_cals(), f_opt(), f_real_work(), LambdaCReconstruction::LambdaCReconstruction(), main(), main(), MdcHoughFinder::MdcHoughFinder(), MdcMergeDups::MdcMergeDups(), MdcNavigation::MdcNavigation(), MdcTrkRecon::MdcTrkRecon(), MdcUtilitySvc::MdcUtilitySvc(), MdcxCosmicSewer::MdcxCosmicSewer(), MdcxTrackFinder::MdcxTrackFinder(), NeutralDReconstruction::NeutralDReconstruction(), QCMCFilter::QCMCFilter(), and t_to_h().

m_diff_standard

const double m_diff_standard = 0.02

Definition at line 5 of file draw_charge_space_coarse_II.cxx.

Referenced by draw_charge_space_coarse_II().

m_err_standard

const double m_err_standard = 0.005

Definition at line 4 of file draw_charge_space_coarse_II.cxx.

Referenced by draw_charge_space_coarse_II().

m_num_hists

int m_num_hists = 0

Definition at line 7 of file draw_charge_space_coarse_II.cxx.

Referenced by draw_charge_space_coarse_II(), and f_real_work().

m_num_in_bin

int m_num_in_bin = 100

Definition at line 6 of file draw_charge_space_coarse_II.cxx.

Referenced by draw_charge_space_coarse_II(), and f_real_work().

x

double x[1000]

Definition at line 8 of file draw_charge_space_coarse_II.cxx.

Referenced by draw_charge_space_coarse_II(), and f_real_work().

xerr

double xerr[1000]

Definition at line 8 of file draw_charge_space_coarse_II.cxx.

Referenced by XtCalib::calib(), draw_charge_space_coarse_II(), f_real_work(), and XtMdcCalib::updateConst().

y

double y[1000]

Definition at line 8 of file draw_charge_space_coarse_II.cxx.

Referenced by TrackPool::AddMissTrack(), TrackPool::AddMissTrack(), BesGeoTrack::AddPoint(), TrackPool::AddTrack(), TrackPool::AddTrackVertex(), Dalitz::Amplitude(), EvtVector4C::applyRotateEuler(), EvtVector4R::applyRotateEuler(), KalmanFit::Helix::approach(), TRunge::approach_line(), TRunge::approach_point(), Dalitz::Babar_Amplitude(), BesMarker2D::BesMarker2D(), BesMarker2D::BesMarker2D(), EvtBtoXsgammaFermiUtil::BesselI1(), EvtBtoXsgammaFermiUtil::BesselK1(), G4HepMCInterface::Boost(), DifFourVector::boostFrom(), DifFourVector::boostFromMe(), DifFourVector::boostTo(), DifFourVector::boostToMe(), Bfield::bx(), Bfield::bx(), Bfield::bx(), Bfield::by(), Bfield::by(), Bfield::by(), Bfield::bz(), Bfield::bz(), Bfield::bz(), MdcUtilitySvc::cellTrackPassed(), HTrackParameter::center(), CFCir::CFCir(), TLine0::chi2(), TMLine::chi2(), Dalitz::CLEO_Amplitude(), EvtCubicSpline::Complex_Derivative(), HoughHit::conformalTrans(), BesGeoTrack::Construct3DLine(), BesTofConstruction::ConstructBr1Tof(), BesTofConstruction::ConstructBr2Tof(), BesEvent::ConstructEmcTrackFromRec(), BesEvent::ConstructMdcTrackFromRec(), BesEvent::ConstructMucTrackFromRec(), BesEvent::ConstructTofTrackFromRec(), BesPointErr::covRTPMatrix(), BesVectorErr::covRTPMatrix(), BesPointErr::covRZPMatrix(), BesVectorErr::covRZPMatrix(), Lpar::d(), EmcTrackCnv::DataObjectToTObject(), RecEmcShowerCnv::DataObjectToTObject(), RecExtTrackCnv::DataObjectToTObject(), EvtEtaDalitz::decay(), EvtVectorIsr::decay(), EvtVub::decay(), EvtVubHybrid::decay(), Lpav::delta_chisq(), BesPointErr::determineChisq(), BesVectorErr::determineChisq(), HTrackParameter::dHdp(), HTrackParameter::dHdx(), BesCircle2D::DistancetoPrimitive(), doLineFit(), EvtPHOTOS::doRadCorr(), BesGMenuTitle::DoRedraw(), Lpar::dr(), draw_charge_space_coarse_II(), draw_check_curve(), draw_plots(), EFEventLoopMgr::EF_Version(), TRunge::eloss(), EvtDecayAngleChi(), EvtDecayAnglePhi(), EvtVector3R::EvtVector3R(), CosmicGenerator::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), DQASelHadron::execute(), EventAssemblyAlg::execute(), QCMCFilter::execute(), Bes2DView::ExecuteEvent(), BesView::ExecuteEvent(), BesTView::ExecuteRotateView(), RecMucTrack::Extend(), Lpav::extrapolate(), f_real_work(), EvtBtoXsgammaFermiUtil::FermiExpFunc(), EvtBtoXsgammaFermiUtil::FermiGaussFunc(), EvtBtoXsgammaFermiUtil::FermiGaussRootFcnA(), EvtBtoXsgammaFermiUtil::FermiGaussRootFcnB(), EvtBtoXsgammaFermiUtil::FermiRomanFunc(), EvtBtoXsgammaFermiUtil::FermiRomanRootFcnA(), Bfield::fieldMap(), Bfield::fieldMap(), Bfield::fieldMap(), Bfield::fieldMap(), Bfield::fieldMap(), KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), MucGeoGeneral::FindIntersection(), BesTView::FindNormal(), MucRec2DRoad::Fit(), Lpav::fit(), TLine2D::fit(), TLineFitter::fit(), TLine0::fit2(), TMLine::fit2(), TLine0::fit2p(), TMLine::fit2p(), TLine0::fit2s(), TMLine::fit2s(), TLine0::fit2sp(), TMLine::fit2sp(), HadronCalibration::fitBGCurve(), TRunge::Fly(), TRunge::Fly_SC(), EvtCubicSpline::fprime(), EvtCubicSpline::fprime(), TRunge::Function(), EvtBtoXsgammaFermiUtil::Gamma(), NumRecipes::gammln(), EvtGen::generateEvent(), TRadGlobal::Get_Var_b(), TRadGlobal::Get_VarPar(), Emc2DCrystal::GetAngle(), Muc2DStrip::GetAngle(), Tof2DScin::GetAngle(), MucGeoStrip::GetCenterPos(), TRadSpline3::GetCoeff(), TRadSpline5::GetCoeff(), HoughHit::getConformal_r(), HoughHit::getConformal_u(), HoughHit::getConformal_v(), HelixTraj::getDFInfo(), NeutTraj::getDFInfo(), TrkCircleTraj::getDFInfo(), TrkDifLineTraj::getDFInfo(), HelixTraj::getDFInfo2(), TrkCircleTraj::getDFInfo2(), BesTView::GetDistancetoAxis(), KalFitTrack::getDriftTime(), BesMagneticField::GetFieldValue(), ExtBesMagneticField::GetFieldValue(), TRadSpline3::GetKnot(), TRadSpline5::GetKnot(), TRadSplinePoly::GetKnot(), MucMagneticField::getMucField(), MucGeoGap::GetRotationMatrix(), MucRec2DRoad::GetVertexPos(), TRunge::GetXP(), MucGeoGap::GuessStrip(), MdcUtilitySvc::Hel(), EvtHelSys::Helrotate(), G4HepMCInterface::HepMC2G4(), HTrackParameter::HTrackParameter(), EvtCPUtil::incoherentMix(), EvtVSSBMixCPT::init(), MucROOTGeo::Init2DGeometry(), BesTofDigitizerEcV4::HitStruct::initial(), QCMCFilter::initialize(), InitV(), TRungeFitter::innerwall(), ParticleIDBase::interpolation(), KalFitTrack::intersect_zx_plane(), TRunge::intersect_zx_plane(), MucGeoGap::IsInGap(), MucGeoStrip::IsInStrip(), EvtConExc::islgr(), HoughTrack::Leastfit(), HoughZsFit::leastFit(), HoughZsFit::leastLine(), PartProduce::lg_bareChamber(), EvtConExc::lgr(), MucRecLineFit::LineFit(), MucRecLineFit::LineFit(), EvtConExc::LLr(), main(), HoughHit::makeCir(), TTrackManager::maskCurlHits(), TTrackManager::merge(), EvtConExc::Mhad_sampling(), HTrackParameter::minDistanceTwoHelix(), ExthIonisation::MinPrimaryEnergy(), Bes2DView::Move(), BesView::Move(), MucRec3DRoad::MucRec3DRoad(), MucRecHit::MucRecHit(), MucRecHit::MucRecHit(), DedxPID::mypol3(), TofCPID::mypol3(), DedxPID::mypol5(), TofCPID::mypol5(), BesTView::NormalWCtoNDC(), BesTView::NormalWCtoNDC(), operator+(), operator+(), operator-(), operator==(), EvtTwoBodyKine::p(), BesTView::PadRange(), BesGeoTrack::Paint(), Dalitz::Phase(), Lpar::phi(), ZHelix::Phi2XYZ(), ZHelix::Phi2ZR(), MdcTrackList::pickHits(), KalFitTrack::pivot_numf(), KalFitTrack::pivot_numf(), Dalitz::Point_on_DP(), Dalitz::Point_on_DP2(), MdcUtilitySvc::pointOnHelixPatPar(), ParticleIDBase::pol2(), ParticleIDBase::pol3(), ParticleIDBase::pol4(), EvtPolInt::polynomial(), Gen_HEPEVT::position(), HelixTraj::position(), TrkDifLineTraj::position(), G4HepMCInterface::Print(), BesTofDigitizerEcV4::HitStruct::print(), MucRec3DRoad::PrintHitsInfo(), MucRec2DRoad::Project(), MucRec3DRoad::Project(), RecMucTrack::Project(), MucRec3DRoad::ProjectNoCurrentGap(), MucRec3DRoad::ProjectWithSigma(), TRunge::Propagate(), TRunge::Propagate1(), TRunge::Propagate_QC(), MucRecQuadFit::QuadFit(), FTTrack::r_phi2Fit(), FTTrack::r_phi3Fit(), FTTrack::r_phi4Fit(), FTTrack::r_phiFit(), FTTrack::r_phiReFit(), KalFitTrack::radius_numf(), EvtStdlibRandomEngine::random(), EvtPolInt::ratint(), BesTofDigitizerEcV3::ReadEffTree(), readTree(), TLine0::reducedChi2(), TMLine::reducedChi2(), Lpar::s(), Lpar::sd(), EvtVector3C::set(), EvtVector3C::set(), EvtVector3R::set(), MucBakelite::SetBarrelRpcInBes(), MucGas::SetBarrelRpcInBes(), MucRpc::SetBarrelRpcInLoc(), MucStrip::SetBarrelStripInLoc(), BesCircle2D::SetCenter(), RecMucTrack::SetCurrentDir(), RecMucTrack::SetCurrentInsct(), RecMucTrack::SetCurrentPos(), MucBakelite::SetEndcapRpcInBes(), MucGas::SetEndcapRpcInBes(), MucRpc::SetEndcapRpcInBes(), MucStrip::SetEndcapStripInBes(), RecMucTrack::SetExtMucPos(), TRunge::SetFirst(), MucRecHit::SetIntersectY(), CDPhoton::setIP(), MucAbsorber::SetLocOrgInBes(), MucBox::SetLocOrgInBes(), MucBoxCover::SetLocOrgInBes(), MucGap::SetLocOrgInBes(), MucRpc::SetLocOrgInBes(), MucStrip::SetLocOrgInBes(), MucStripPlane::SetLocOrgInBes(), BesGeoTrack::SetMarker(), RecMucTrack::SetMdcPos(), RecMucTrack::SetMucPos(), RecMucTrack::SetMucPosSigma(), BesPrimaryGeneratorMessenger::SetNewValue(), ZHelix::SetNextPoint(), MucAbsorber::SetObjOrgInBes(), MucBox::SetObjOrgInBes(), MucBoxCover::SetObjOrgInBes(), MucGap::SetObjOrgInBes(), MucStripPlane::SetObjOrgInBes(), MucStripPlane::SetObjOrgInLoc(), BesMarker2D::SetParameter(), ZHelix::SetPoint(), ZHelix::SetPoints(), BesPaveText::SetPos(), BesPrimaryGeneratorAction::SetPosY(), TLorentzVectorC::SetPxPyPzE(), MdcxHel::SetRef(), DstMucTrack::setYPos(), TMucTrack::setYPos(), TRecMucTrack::setYPos(), TMucTrack::setYPosSigma(), TRecMucTrack::setYPosSigma(), Lpav::similarity(), MdcDedxTrk::SpaceChargeCorrec(), SpaceChargeCorrec(), EvtPolInt::splint(), BesEmcConstruction::ThreeVectorTrans(), TLorentzVectorC::TLorentzVectorC(), EmcTrackCnv::TObjectToDataObject(), MdcKalTrackCnv::TObjectToDataObject(), MdcTrackCnv::TObjectToDataObject(), RecEmcShowerCnv::TObjectToDataObject(), RecExtTrackCnv::TObjectToDataObject(), RecMdcKalTrackCnv::TObjectToDataObject(), RecMdcTrackCnv::TObjectToDataObject(), TrackInformation(), TrackKinematics(), BesEmcConstruction::TransformToArb8(), FTSegment::update(), FTSegment::update3D(), RkFitCylinder::updateTrack(), ExtSteppingAction::UserSteppingAction(), wireFit(), writeTree(), Dedx_Helix::x(), Dedx_Helix::x(), Ext_Helix::x(), Ext_Helix::x(), Helix::x(), Helix::x(), KalmanFit::Helix::x(), KalmanFit::Helix::x(), VFHelix::x(), VFHelix::x(), TMDCWire::xyPosition(), and ZHelix::Y2Phi().

yerr

double yerr[1000]

Definition at line 8 of file draw_charge_space_coarse_II.cxx.

Referenced by draw_charge_space_coarse_II(), and f_real_work().