33{
34 MsgStream log(
msgSvc(), name());
35 log<<MSG::INFO<<"DedxSaveWireGain::FillHists()"<<endreq;
36 int run, evt;
37 double gain;
38 m_wiregain->Branch("run", &run, "run/I");
39 m_wiregain->Branch("evt", &evt, "evt/I");
40 m_wiregain->Branch("gain", &gain, "gain/D");
41
43 TTree* n102;
44 string runlist;
45
46 double dedx=0;
47 float runNO=0,evtNO=0,runFlag=0,pathlength=0,wid=0,layid=0,dd_in=0,driftdist=0,eangle=0,zhit=0,
costheta=0,tes=0,
ptrk=0;
50 {
52 cout<<"runlist: "<<runlist.c_str()<<endl;
53 f =
new TFile(runlist.c_str());
54 n102 = (TTree*)
f->Get(
"n102");
68 n102->SetBranchAddress(
"ptrk1",&
ptrk);
69
70 cout <<"entries in this file" << n102->GetEntries() << endl;
71 for(int j=0;j<n102->GetEntries();j++)
72 {
73 n102->GetEntry(j);
74 if(tes>1400) continue;
76 if(layid <8)
78 else
80 dedx =
exsvc->
StandardHitCorrec(0,runFlag,2,runNO,evtNO,pathlength,wid,layid,dedx,dd_in,eangle,
costheta);
82 run = (int)runNO;
83 evt = (int)evtNO;
84 gain = dedx;
85
86 m_wiregain->Fill();
87 }
88 }
89}
TFile f("ana_bhabha660a_dqa_mcPat_zy_old.root")
data SetBranchAddress("time",&time)
#define Iner_DriftDistCut
vector< string > m_recFileVector
virtual double StandardHitCorrec(int calib_rec_Flag, int runFlag, int ntpFlag, int runNO, int evtNO, double pathl, int wid, int layid, double adc, double dd, double eangle, double costheta) const =0
virtual double StandardTrackCorrec(int calib_rec_Flag, int typFlag, int ntpFlag, int runNO, int evtNO, double ex, double costheta, double t0) const =0