BOSS 7.0.6
BESIII Offline Software System
Loading...
Searching...
No Matches
Public Member Functions | Static Public Member Functions | Protected Member Functions | Friends | List of all members
TreeDedxCalibDataCnv Class Reference

#include <TreeDedxCalibDataCnv.h>

+ Inheritance diagram for TreeDedxCalibDataCnv:

Public Member Functions

const CLID & objType () const
 
 TreeDedxCalibDataCnv (ISvcLocator *svc)
 
virtual ~TreeDedxCalibDataCnv ()
 
virtual long repSvcType () const
 
- Public Member Functions inherited from TreeCalBaseCnv
virtual ~TreeCalBaseCnv ()
 
virtual StatusCode initialize ()
 
virtual StatusCode finalize ()
 
virtual StatusCode createObj (IOpaqueAddress *addr, DataObject *&refpObject)
 
ICalibTreeSvc * getCalibTreeSvc ()
 
 TreeCalBaseCnv (ISvcLocator *svc, const CLID &clid)
 
virtual StatusCode createRoot (const std::string &fname, CalibData::CalibBase1 *pTDSObj)
 
- Public Member Functions inherited from Converter< Ty1, Ty2 >
destination * operator (const source &) const
 

Static Public Member Functions

static const CLID & classID ()
 
- Static Public Member Functions inherited from TreeCalBaseCnv
static const unsigned char storageType ()
 

Protected Member Functions

virtual StatusCode i_createObj (IOpaqueAddress *address, DataObject *&refpObject)
 
- Protected Member Functions inherited from TreeCalBaseCnv
virtual StatusCode internalCreateObj (DataObject *&refpObject, IOpaqueAddress *address)
 
virtual StatusCode i_createObj (IOpaqueAddress *address, DataObject *&refpObject)
 
virtual StatusCode i_processObj (DataObject *pObject, IOpaqueAddress *address)
 In case there is additional work to do on the created object.
 
virtual StatusCode fillRoot (CalibData::CalibBase *pTDSObj, TObject *pRootObj)
 
void setBaseInfo (CalibData::CalibBase1 *pObj)
 Another utility for derived classes to use.
 
- Protected Member Functions inherited from Converter< Ty1, Ty2 >
virtual destination * convert (const source &) const =0
 

Friends

class CnvFactory< TreeDedxCalibDataCnv >
 

Additional Inherited Members

- Public Types inherited from Converter< Ty1, Ty2 >
typedef Ty1 source
 
typedef Ty2 destination
 
- Protected Attributes inherited from TreeCalBaseCnv
ICalibTreeSvc * m_treeSvc
 
ICalibMetaCnvSvc * m_metaSvc
 
IInstrumentName * m_instrSvc
 
int m_serNo
 
int m_runfrm
 
int m_runto
 
TFile * m_outFile
 
TTree * m_ttree
 
TFile * m_inFile
 
TDirectory * m_saveDir
 

Detailed Description

Base class for CAL calibration converters from Ttrees to TCDS. All such converters need to do certain things, which are handled here. Methods common to all calibrations are in the base class TreeCalBaseCnv

Author
huang bin

Definition at line 20 of file TreeDedxCalibDataCnv.h.

Constructor & Destructor Documentation

◆ TreeDedxCalibDataCnv()

TreeDedxCalibDataCnv::TreeDedxCalibDataCnv ( ISvcLocator *  svc)

Definition at line 32 of file TreeDedxCalibDataCnv.cxx.

32 :
33 TreeCalBaseCnv(svc, CLID_Calib_DedxCal) {
34
35}
CLID_Calib_DedxCal
const CLID CLID_Calib_DedxCal
Definition: CalibModel.h:45

◆ ~TreeDedxCalibDataCnv()

virtual TreeDedxCalibDataCnv::~TreeDedxCalibDataCnv ( )
inlinevirtual

Definition at line 29 of file TreeDedxCalibDataCnv.h.

29{};

Member Function Documentation

◆ classID()

const CLID & TreeDedxCalibDataCnv::classID ( )
static

Definition at line 42 of file TreeDedxCalibDataCnv.cxx.

42 {
43 return CLID_Calib_DedxCal;
44}

◆ i_createObj()

StatusCode TreeDedxCalibDataCnv::i_createObj ( IOpaqueAddress *  address,
DataObject *&  refpObject 
)
protectedvirtual

This creates the transient representation of an object from the corresponding ROOT object. This actually does the "new" operation and deals with the attributes of the node. This base class implementation does nothing; it should not normally be called because it doesn't correspond to any TCDS class. Instead, i_createObj of some derived class will be called.

Parameters
fnameThe ROOT file to be read in to be used to builds the object
refpObjectthe object to be built
Returns
status depending on the completion of the call

Reimplemented from TreeCalBaseCnv.

Definition at line 46 of file TreeDedxCalibDataCnv.cxx.

47 {
48
49 MsgStream log(msgSvc(), "TreeDedxCalibDataCnv");
50 log<<MSG::DEBUG<<"SetProperty"<<endreq;
51
52 CalibData::DedxCalibData *tmpObject = new CalibData::DedxCalibData;
53 TreeAddress* add = dynamic_cast<TreeAddress*>(addr);
54 DatabaseRecord *records=add->pp();
55
56 /*TBuffer *buf1 = new TBuffer(TBuffer::kRead);
57 TBuffer *buf2 = new TBuffer(TBuffer::kRead);
58 TBuffer *buf3 = new TBuffer(TBuffer::kRead);
59 TBuffer *buf4 = new TBuffer(TBuffer::kRead);
60 TBuffer *buf5 = new TBuffer(TBuffer::kRead);
61 TBuffer *buf6 = new TBuffer(TBuffer::kRead);
62 TBuffer *buf7 = new TBuffer(TBuffer::kRead);
63 TBuffer *buf8 = new TBuffer(TBuffer::kRead);
64 TBuffer *buf9 = new TBuffer(TBuffer::kRead);*/
65
66 TBufferFile *buf1 = new TBufferFile(TBuffer::kRead);
67 TBufferFile *buf2 = new TBufferFile(TBuffer::kRead);
68 TBufferFile *buf3 = new TBufferFile(TBuffer::kRead);
69 TBufferFile *buf4 = new TBufferFile(TBuffer::kRead);
70 TBufferFile *buf5 = new TBufferFile(TBuffer::kRead);
71 TBufferFile *buf6 = new TBufferFile(TBuffer::kRead);
72 TBufferFile *buf7 = new TBufferFile(TBuffer::kRead);
73 TBufferFile *buf8 = new TBufferFile(TBuffer::kRead);
74 TBufferFile *buf9 = new TBufferFile(TBuffer::kRead);
75
76 buf1->SetBuffer((*records)["DriftDist"],512000,kFALSE);
77 buf2->SetBuffer((*records)["EntranceAng"],512000,kFALSE);
78 buf3->SetBuffer((*records)["MeanGain"],512000,kFALSE);
79 buf4->SetBuffer((*records)["GasGain"],512000,kFALSE);
80 buf5->SetBuffer((*records)["LayerGain"],512000,kFALSE);
81 buf6->SetBuffer((*records)["Resolution"],512000,kFALSE);
82 buf7->SetBuffer((*records)["WireGain"],512000,kFALSE);
83 buf8->SetBuffer((*records)["ZDep"],512000,kFALSE);
84 buf9->SetBuffer((*records)["RunGain"],512000,kFALSE);
85
86
87
88 std::cout<<" SftVer is "<<(*records)["SftVer"];
89 std::cout<<" CalVerSft is "<<(*records)["CalParVer"];
90 std::cout<<" File name is "<<(*records)["FileName"]<<std::endl;
91
92
93
94 TTree* ddgtree = new TTree();
95 ddgtree->Streamer(*buf1);
96
97 TTree* entratree = new TTree();
98 entratree->Streamer(*buf2);
99
100 TTree* gaintree= new TTree();
101 gaintree->Streamer(*buf3);
102
103 TTree* ggstree = new TTree();
104 ggstree->Streamer(*buf4);
105
106 TTree* layergtree = new TTree();
107 layergtree->Streamer(*buf5);
108
109 TTree* resoltree = new TTree();
110 resoltree->Streamer(*buf6);
111
112 TTree* wiregtree = new TTree();
113 wiregtree->Streamer(*buf7);
114
115 TTree* zdeptree = new TTree();
116 zdeptree->Streamer(*buf8);
117
118 TTree *rungtree = new TTree();
119 rungtree->Streamer(*buf9);
120
121 double rungain;
122 double runmean;
123 int runno, evtfrom, evtto;
124 double runresol;
125
126 rungtree -> SetBranchAddress("rungain", &rungain);
127 rungtree -> SetBranchAddress("runmean", &runmean);
128 rungtree -> SetBranchAddress("runno", &runno);
129 rungtree -> SetBranchAddress("runresol", &runresol);
130 if(rungtree->GetBranchStatus("evtfrom")){
131 rungtree -> SetBranchAddress("evtfrom", &evtfrom);
132 rungtree -> SetBranchAddress("evtto", &evtto);
133 }
134 else{
135 evtfrom = 0;
136 evtto = 1000000000;
137 }
138 int N = rungtree -> GetEntries();
139 tmpObject -> setrunNO(N);
140 for(int i=0; i<N; i++){
141 rungtree -> GetEntry(i);
142 tmpObject -> setrung(rungain,0);
143 tmpObject -> setrung(runmean,1);
144 tmpObject -> setrung(runno,2);
145 tmpObject -> setrung(runresol,3);
146 tmpObject -> setrung(evtfrom,4);
147 tmpObject -> setrung(evtto,5);
148 }
149
150
151 double ddg0[43];
152 double ddg1[43];
153 double ddg2[43];
154 double ddg3[43];
155 double id_doca[1600];
156 double iner_chi[1600];
157 double iner_gain[1600];
158 double iner_hits[1600];
159 double ip_eangle[1600];
160 double out_chi[1600];
161 double out_gain[1600];
162 double out_hits[1600];
163
164 ddgtree -> SetBranchAddress("ddg0", ddg0);
165 ddgtree -> SetBranchAddress("ddg1", ddg1);
166 ddgtree -> SetBranchAddress("ddg2", ddg2);
167 ddgtree -> SetBranchAddress("ddg3", ddg3);
168 TBranch *bbb = ddgtree->FindBranch("Id_doca");
169 if(bbb){
170 ddgtree -> SetBranchAddress("Id_doca", id_doca);
171 ddgtree -> SetBranchAddress("Iner_chi", iner_chi);
172 ddgtree -> SetBranchAddress("Iner_gain", iner_gain);
173 ddgtree -> SetBranchAddress("Iner_hits", iner_hits);
174 ddgtree -> SetBranchAddress("Ip_eangle", ip_eangle);
175 ddgtree -> SetBranchAddress("Out_chi", out_chi);
176 ddgtree -> SetBranchAddress("Out_gain", out_gain);
177 ddgtree -> SetBranchAddress("Out_hits", out_hits);
178 }
179
180 ddgtree -> GetEntry(0);
181 for(int i=0; i<43; i++){
182
183 tmpObject -> setddg(ddg0[i],0,i);
184 tmpObject -> setddg(ddg1[i],1,i);
185 tmpObject -> setddg(ddg2[i],2,i);
186 tmpObject -> setddg(ddg3[i],3,i);
187 }
188
189 for(int i=0; i<1600; i++){
190 if(!bbb){
191 id_doca[i]=0;
192 iner_chi[i]=0;
193 iner_gain[i]=0;
194 iner_hits[i]=0;
195 ip_eangle[i]=0;
196 out_chi[i]=0;
197 out_gain[i]=0;
198 out_hits[i]=0;
199 }
200 tmpObject -> set_id_doca(id_doca[i],i);
201 tmpObject -> set_iner_chi(iner_chi[i],i);
202 tmpObject -> set_iner_gain(iner_gain[i],i);
203 tmpObject -> set_iner_hits(iner_hits[i],i);
204 tmpObject -> set_ip_eangle(ip_eangle[i],i);
205 tmpObject -> set_out_chi(out_chi[i],i);
206 tmpObject -> set_out_gain(out_gain[i],i);
207 tmpObject -> set_out_hits(out_hits[i],i);
208 }
209
210
211 double entra0[43];
212 double entra1[43];
213 double entra2[43];
214 double entra3[43];
215 double engle[100];
216 int engle_no;
217 entratree -> SetBranchAddress("entra0", entra0);
218 entratree -> SetBranchAddress("entra1", entra1);
219 entratree -> SetBranchAddress("entra2", entra2);
220 entratree -> SetBranchAddress("entra3", entra3);
221 entratree -> SetBranchAddress("1denangle", engle);
222 entratree -> SetBranchAddress("1denangle_entry", &engle_no);
223 entratree -> GetEntry(0);
224 for(int i=0; i<43; i++){
225 tmpObject -> setenta(entra0[i],0,i);
226 tmpObject -> setenta(entra1[i],1,i);
227 tmpObject -> setenta(entra2[i],2,i);
228 tmpObject -> setenta(entra3[i],3,i);
229 }
230
231 tmpObject -> set_enanglesize(engle_no);
232 for(int i=0; i<engle_no; i++){
233 tmpObject -> set_enangle(engle[i],i);
234 }
235
236
237
238
239// read ggscalib ------------------------------------------------------------
240 double ggs0[43];
241 double ggs1[43];
242 double ggs2[43];
243 double ggs3[43];
244 double gcostheta[80];
245 int hadron_entry;
246 double hadron[20];
247 ggstree -> SetBranchAddress("ggs0", ggs0);
248 ggstree -> SetBranchAddress("ggs1", ggs1);
249 ggstree -> SetBranchAddress("ggs2", ggs2);
250 ggstree -> SetBranchAddress("ggs3", ggs3);
251 ggstree -> SetBranchAddress("hadron", hadron);
252 ggstree -> SetBranchAddress("hadronNo", &hadron_entry);
253 if(bbb){
254 ggstree -> SetBranchAddress("costheta", gcostheta);}
255 ggstree -> GetEntry(0);
256 for(int i=0; i<43;i++){
257 tmpObject -> setggs(ggs0[i],0,i);
258 tmpObject -> setggs(ggs1[i],1,i);
259 tmpObject -> setggs(ggs2[i],2,i);
260 tmpObject -> setggs(ggs3[i],3,i);
261 }
262 for(int i=0; i<80;i++){
263 if(!bbb) gcostheta[i]=0;
264 tmpObject ->set_costheta(gcostheta[i],i);
265 }
266 if(hadron_entry>20){
267 log<<MSG::FATAL<<"hadron entry is larger than 20, larger than designed"<<endreq;
268 return StatusCode::FAILURE;
269 }
270 tmpObject->set_hadronNo(hadron_entry);
271 for(int i=0;i<hadron_entry;i++){
272 tmpObject->set_hadron(hadron[i],i);
273 }
274
275
276// read zdepcalib ------------------------------------------------------------
277 double zdep0[43];
278 double zdep1[43];
279 double zdep2[43];
280 double zdep3[43];
281 zdeptree -> SetBranchAddress("zdep0", zdep0);
282 zdeptree -> SetBranchAddress("zdep1", zdep1);
283 zdeptree -> SetBranchAddress("zdep2", zdep2);
284 zdeptree -> SetBranchAddress("zdep3", zdep3);
285 zdeptree -> GetEntry(0);
286
287 for(int i=0; i<43;i++){
288 tmpObject -> setzdep(zdep0[i],0,i);
289 tmpObject -> setzdep(zdep1[i],1,i);
290 tmpObject -> setzdep(zdep2[i],2,i);
291 tmpObject -> setzdep(zdep3[i],3,i);
292 }
293// read gain ------------------------------------------------------------
294 double gain;
295 double gt0[35],gdedx[35];
296 gaintree -> SetBranchAddress("gain", &gain);
297 if(bbb){
298 gaintree -> SetBranchAddress("t0", gt0);
299 gaintree -> SetBranchAddress("dedx", gdedx);}
300 gaintree -> GetEntry(0);
301 tmpObject -> setgain(gain);
302 for(int i=0; i<35;i++){
303 if(!bbb){
304 gt0[i]=0;
305 gdedx[i]=0;
306 }
307 tmpObject->set_t0(gt0[i],i);
308 tmpObject->set_dedx(gdedx[i],i);
309 }
310
311// read resol----------------------------------------------------------------
312 double resol;
313 resoltree -> SetBranchAddress("resol", &resol);
314 resoltree -> GetEntry(0);
315 tmpObject -> setresol(resol);
316
317// read wireg----------------------------------------------------------------
318 double wireg[6796];
319 wiregtree -> SetBranchAddress("wireg",wireg);
320 wiregtree -> GetEntry(0);
321 for(int i=0;i<6796;i++){
322 if(wireg[i]>5||wireg[i]<0.2) wireg[i]=0;
323 tmpObject -> setwireg(wireg[i],i);
324 }
325
326// read layerg----------------------------------------------------------------
327 double layerg[43];
328 layergtree -> SetBranchAddress("layerg",layerg);
329 layergtree -> GetEntry(0);
330
331 for(int i=0;i<43;i++){
332 tmpObject -> setlayerg(layerg[i],i);
333 }
334
335
336
337 // Read in the object
338
339
340 refpObject=tmpObject;
341
342 delete ddgtree;
343 delete entratree;
344 delete gaintree;
345 delete ggstree;
346 delete layergtree;
347 delete resoltree;
348 delete wiregtree;
349 delete zdeptree;
350 delete rungtree;
351
352 return StatusCode::SUCCESS;
353
354
355}
msgSvc
IMessageSvc * msgSvc()
Definition: ServiceAccessor.h:13
CalibData::DedxCalibData::set_costheta
void set_costheta(const double aa, int i)
Definition: DedxCalibData.h:107
CalibData::DedxCalibData::set_hadron
void set_hadron(const double aa, int i)
Definition: DedxCalibData.h:110
CalibData::DedxCalibData::set_t0
void set_t0(const double aa, int i)
Definition: DedxCalibData.h:108
CalibData::DedxCalibData::set_dedx
void set_dedx(const double aa, int i)
Definition: DedxCalibData.h:109
TreeAddress::pp
virtual DatabaseRecord * pp()
Definition: TreeAddress.h:130
GetEntry
curve GetEntry(0)
SetBranchAddress
curve SetBranchAddress("CurveSize",&CurveSize)

◆ objType()

const CLID & TreeDedxCalibDataCnv::objType ( ) const

Definition at line 38 of file TreeDedxCalibDataCnv.cxx.

38 {
39 return CLID_Calib_DedxCal;
40}

◆ repSvcType()

virtual long TreeDedxCalibDataCnv::repSvcType ( ) const
inlinevirtual

Definition at line 34 of file TreeDedxCalibDataCnv.h.

34 {
35 return CALIBTREE_StorageType;
36 }
CALIBTREE_StorageType
unsigned const char CALIBTREE_StorageType
Definition: ICalibTreeSvc.h:20

Friends And Related Function Documentation

◆ CnvFactory< TreeDedxCalibDataCnv >

friend class CnvFactory< TreeDedxCalibDataCnv >
friend

Definition at line 1 of file TreeDedxCalibDataCnv.h.

The documentation for this class was generated from the following files: