00001 #include "EmcClusterLocalExtMicrov4.h"
00002 #include "EmcSnglClusterLocalExtMicrov4.h"
00003
00004 #include "TClonesArray.h"
00005 #include <iostream>
00006
00007 ClassImp(EmcClusterLocalExtMicrov4)
00008
00009 using namespace std;
00010
00011 static const unsigned int EMCNCLUSV4 = 200;
00012
00013 EmcClusterLocalExtMicrov4::EmcClusterLocalExtMicrov4()
00014 {
00015 EmcNCluster = 0;
00016 EmcClus = new TClonesArray("EmcSnglClusterLocalExtMicrov4",EMCNCLUSV4);
00017 return;
00018 }
00019
00020 EmcClusterLocalExtMicrov4::~EmcClusterLocalExtMicrov4()
00021 {
00022 if (EmcClus)
00023 {
00024 EmcClus->Clear();
00025 delete EmcClus;
00026 }
00027 return;
00028 }
00029
00030 void EmcClusterLocalExtMicrov4::identify(ostream& os) const
00031 {
00032 os << "identify yourself: EmcClusterLocalExtMicrov4 Object" << std::endl;
00033 os << "No of Clusters: " << EmcNCluster << std::endl;
00034 return;
00035 }
00036
00037 void EmcClusterLocalExtMicrov4::Reset()
00038 {
00039 EmcClus->Clear();
00040 if (EmcNCluster>EMCNCLUSV4)
00041 {
00042 EmcClus->Expand(EMCNCLUSV4);
00043 }
00044 EmcNCluster = 0;
00045 return;
00046 }
00047
00048 int EmcClusterLocalExtMicrov4::isValid() const
00049 {
00050 return((EmcNCluster>0) ? 1 : 0);
00051 }
00052
00053 int EmcClusterLocalExtMicrov4::set_TClonesArraySize(unsigned int nclus)
00054 {
00055 if (nclus > EMCNCLUSV4)
00056 {
00057 EmcClus->Expand(nclus);
00058 }
00059 return nclus;
00060 }
00061
00062 void EmcClusterLocalExtMicrov4::AddEmcCluster(unsigned int iclus)
00063 {
00064 TClonesArray &emcclus = *EmcClus;
00065 new(emcclus[iclus]) EmcSnglClusterLocalExtMicrov4();
00066 return;
00067 }
00068
00069 int EmcClusterLocalExtMicrov4::get_index(unsigned int iclus) const
00070 {
00071
00072 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00073 return((emcclus) ? emcclus->get_index() : -999);
00074 }
00075
00076 void EmcClusterLocalExtMicrov4::set_index(unsigned int iclus, int ival)
00077 {
00078 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00079 if (emcclus)
00080 {
00081 emcclus->set_index(ival);
00082 }
00083 else
00084 {
00085 std::cout << PHWHERE << "ERROR no EmcSnglClusterLocalExtMicrov4 object found" << std::endl;
00086 }
00087 return;
00088 }
00089
00090 int EmcClusterLocalExtMicrov4::get_warnmap(unsigned int iclus) const
00091 {
00092
00093 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00094 return((emcclus) ? emcclus->get_warnmap() : -999);
00095 }
00096
00097 void EmcClusterLocalExtMicrov4::set_warnmap(unsigned int iclus, int ival)
00098 {
00099 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00100 if (emcclus)
00101 {
00102 emcclus->set_warnmap(ival);
00103 }
00104 else
00105 {
00106 std::cout << PHWHERE << "ERROR no EmcSnglClusterLocalExtMicrov4 object found" << std::endl;
00107 }
00108 return;
00109 }
00110
00111 int EmcClusterLocalExtMicrov4::get_deadmap(unsigned int iclus) const
00112 {
00113
00114 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00115 return((emcclus) ? emcclus->get_deadmap() : -999);
00116 }
00117
00118 void EmcClusterLocalExtMicrov4::set_deadmap(unsigned int iclus, int ival)
00119 {
00120 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00121 if (emcclus)
00122 {
00123 emcclus->set_deadmap(ival);
00124 }
00125 else
00126 {
00127 std::cout << PHWHERE << "ERROR no EmcSnglClusterLocalExtMicrov4 object found" << std::endl;
00128 }
00129 return;
00130 }
00131
00132
00133 short EmcClusterLocalExtMicrov4::get_arm(const unsigned int iclus) const
00134 {
00135
00136 int key=get_index(iclus);
00137 return key/100000;
00138
00139 }
00140
00141 short EmcClusterLocalExtMicrov4::get_sector(const unsigned int iclus) const
00142 {
00143
00144 int key=get_index(iclus);
00145 return (key%100000)/10000;
00146
00147 }
00148
00149 short EmcClusterLocalExtMicrov4::get_ind(const unsigned int iclus, const short i) const
00150 {
00151
00152
00153
00154
00155 int key=get_index(iclus);
00156
00157 if(i==0) return key%100;
00158 else if(i==1) return (key%10000)/100;
00159
00160 std::cout<<PHWHERE<<"ERROR index "<<i<<"out of range"<<std::endl;
00161 return -999;
00162
00163 }
00164
00165
00166 short EmcClusterLocalExtMicrov4::get_twrhit(unsigned int iclus) const
00167 {
00168
00169 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00170 return((emcclus) ? emcclus->get_twrhit() : -999);
00171 }
00172
00173 void EmcClusterLocalExtMicrov4::set_twrhit(unsigned int iclus, short ival)
00174 {
00175 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00176 if (emcclus)
00177 {
00178 emcclus->set_twrhit(ival);
00179 }
00180 else
00181 {
00182 std::cout << PHWHERE << "ERROR no EmcSnglClusterLocalExtMicrov4 object found" << std::endl;
00183 }
00184 return;
00185 }
00186
00187 float EmcClusterLocalExtMicrov4::get_chi2(unsigned int iclus) const
00188 {
00189
00190 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00191 return((emcclus) ? emcclus->get_chi2() : -9999.9);
00192 }
00193
00194 void EmcClusterLocalExtMicrov4::set_chi2(unsigned int iclus, float rval)
00195 {
00196 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00197 if (emcclus)
00198 {
00199 emcclus->set_chi2(rval);
00200 }
00201 else
00202 {
00203 std::cout << PHWHERE << "ERROR no EmcSnglClusterLocalExtMicrov4 object found" << std::endl;
00204 }
00205 return;
00206 }
00207
00208 float EmcClusterLocalExtMicrov4::get_qual(unsigned int iclus) const
00209 {
00210
00211 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00212 return((emcclus) ? emcclus->get_qual() : -9999.9);
00213 }
00214
00215 void EmcClusterLocalExtMicrov4::set_qual(unsigned int iclus, float rval)
00216 {
00217 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00218 if (emcclus)
00219 {
00220 emcclus->set_qual(rval);
00221 }
00222 else
00223 {
00224 std::cout << PHWHERE << "ERROR no EmcSnglClusterLocalExtMicrov4 object found" << std::endl;
00225 }
00226 return;
00227 }
00228
00229 float EmcClusterLocalExtMicrov4::get_chi2_sh(unsigned int iclus) const
00230 {
00231
00232 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00233 return((emcclus) ? emcclus->get_chi2_sh() : -9999.9);
00234 }
00235
00236 void EmcClusterLocalExtMicrov4::set_chi2_sh(unsigned int iclus, float rval)
00237 {
00238 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00239 if (emcclus)
00240 {
00241 emcclus->set_chi2_sh(rval);
00242 }
00243 else
00244 {
00245 std::cout << PHWHERE << "ERROR no EmcSnglClusterLocalExtMicrov4 object found" << std::endl;
00246 }
00247 return;
00248 }
00249
00250 float EmcClusterLocalExtMicrov4::get_e(unsigned int iclus) const
00251 {
00252
00253 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00254 return((emcclus) ? emcclus->get_e() : -9999.9);
00255 }
00256
00257 void EmcClusterLocalExtMicrov4::set_e(unsigned int iclus, float rval)
00258 {
00259 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00260 if (emcclus)
00261 {
00262 emcclus->set_e(rval);
00263 }
00264 else
00265 {
00266 std::cout << PHWHERE << "ERROR no EmcSnglClusterLocalExtMicrov4 object found" << std::endl;
00267 }
00268 return;
00269 }
00270
00271 float EmcClusterLocalExtMicrov4::get_ecore(unsigned int iclus) const
00272 {
00273
00274 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00275 return((emcclus) ? emcclus->get_ecore() : -9999.9);
00276 }
00277
00278 void EmcClusterLocalExtMicrov4::set_ecore(unsigned int iclus, float rval)
00279 {
00280 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00281 if (emcclus)
00282 {
00283 emcclus->set_ecore(rval);
00284 }
00285 else
00286 {
00287 std::cout << PHWHERE << "ERROR no EmcSnglClusterLocalExtMicrov4 object found" << std::endl;
00288 }
00289 return;
00290 }
00291
00292 float EmcClusterLocalExtMicrov4::get_ecorr(unsigned int iclus) const
00293 {
00294
00295 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00296 return((emcclus) ? emcclus->get_ecorr() : -9999.9);
00297 }
00298
00299 void EmcClusterLocalExtMicrov4::set_ecorr(unsigned int iclus, float rval)
00300 {
00301 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00302 if (emcclus)
00303 {
00304 emcclus->set_ecorr(rval);
00305 }
00306 else
00307 {
00308 std::cout << PHWHERE << "ERROR no EmcSnglClusterLocalExtMicrov4 object found" << std::endl;
00309 }
00310 return;
00311 }
00312
00313 float EmcClusterLocalExtMicrov4::get_ecent(unsigned int iclus) const
00314 {
00315
00316 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00317 return((emcclus) ? emcclus->get_ecent() : -9999.9);
00318 }
00319
00320 void EmcClusterLocalExtMicrov4::set_ecent(unsigned int iclus, float rval)
00321 {
00322 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00323 if (emcclus)
00324 {
00325 emcclus->set_ecent(rval);
00326 }
00327 else
00328 {
00329 std::cout << PHWHERE << "ERROR no EmcSnglClusterLocalExtMicrov4 object found" << std::endl;
00330 }
00331 return;
00332 }
00333
00334 float EmcClusterLocalExtMicrov4::get_e9(unsigned int iclus) const
00335 {
00336
00337 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00338 return((emcclus) ? emcclus->get_e9() : -9999.9);
00339 }
00340
00341 void EmcClusterLocalExtMicrov4::set_e9(unsigned int iclus, float rval)
00342 {
00343 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00344 if (emcclus)
00345 {
00346 emcclus->set_e9(rval);
00347 }
00348 else
00349 {
00350 std::cout << PHWHERE << "ERROR no EmcSnglClusterLocalExtMicrov4 object found" << std::endl;
00351 }
00352 return;
00353 }
00354
00355 float EmcClusterLocalExtMicrov4::get_prob_photon(unsigned int iclus) const
00356 {
00357
00358 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00359 return((emcclus) ? emcclus->get_prob_photon() : -9999.9);
00360 }
00361
00362 void EmcClusterLocalExtMicrov4::set_prob_photon(unsigned int iclus, float rval)
00363 {
00364 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00365 if (emcclus)
00366 {
00367 emcclus->set_prob_photon(rval);
00368 }
00369 else
00370 {
00371 std::cout << PHWHERE << "ERROR no EmcSnglClusterLocalExtMicrov4 object found" << std::endl;
00372 }
00373 return;
00374 }
00375
00376 float EmcClusterLocalExtMicrov4::get_prob_photon_sh(unsigned int iclus) const
00377 {
00378
00379 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00380 return((emcclus) ? emcclus->get_prob_photon_sh() : -9999.9);
00381 }
00382
00383 void EmcClusterLocalExtMicrov4::set_prob_photon_sh(unsigned int iclus, float rval)
00384 {
00385 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00386 if (emcclus)
00387 {
00388 emcclus->set_prob_photon_sh(rval);
00389 }
00390 else
00391 {
00392 std::cout << PHWHERE << "ERROR no EmcSnglClusterLocalExtMicrov4 object found" << std::endl;
00393 }
00394 return;
00395 }
00396
00397 float EmcClusterLocalExtMicrov4::get_re9(unsigned int iclus) const
00398 {
00399
00400 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00401 return((emcclus) ? emcclus->get_re9() : -9999.9);
00402 }
00403
00404 void EmcClusterLocalExtMicrov4::set_re9(unsigned int iclus, float rval)
00405 {
00406 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00407 if (emcclus)
00408 {
00409 emcclus->set_re9(rval);
00410 }
00411 else
00412 {
00413 std::cout << PHWHERE << "ERROR no EmcSnglClusterLocalExtMicrov4 object found" << std::endl;
00414 }
00415 return;
00416 }
00417
00418 float EmcClusterLocalExtMicrov4::get_tofcorr(unsigned int iclus) const
00419 {
00420
00421 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00422 return((emcclus) ? emcclus->get_tofcorr() : -9999.9);
00423 }
00424
00425 void EmcClusterLocalExtMicrov4::set_tofcorr(unsigned int iclus, float rval)
00426 {
00427 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00428 if (emcclus)
00429 {
00430 emcclus->set_tofcorr(rval);
00431 }
00432 else
00433 {
00434 std::cout << PHWHERE << "ERROR no EmcSnglClusterLocalExtMicrov4 object found" << std::endl;
00435 }
00436 return;
00437 }
00438
00439 float EmcClusterLocalExtMicrov4::get_tofmin(unsigned int iclus) const
00440 {
00441
00442 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00443 return((emcclus) ? emcclus->get_tofmin() : -9999.9);
00444 }
00445
00446 void EmcClusterLocalExtMicrov4::set_tofmin(unsigned int iclus, float rval)
00447 {
00448 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00449 if (emcclus)
00450 {
00451 emcclus->set_tofmin(rval);
00452 }
00453 else
00454 {
00455 std::cout << PHWHERE << "ERROR no EmcSnglClusterLocalExtMicrov4 object found" << std::endl;
00456 }
00457 return;
00458 }
00459
00460 float EmcClusterLocalExtMicrov4::get_tofmax(unsigned int iclus) const
00461 {
00462
00463 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00464 return((emcclus) ? emcclus->get_tofmax() : -9999.9);
00465 }
00466
00467 void EmcClusterLocalExtMicrov4::set_tofmax(unsigned int iclus, float rval)
00468 {
00469 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00470 if (emcclus)
00471 {
00472 emcclus->set_tofmax(rval);
00473 }
00474 else
00475 {
00476 std::cout << PHWHERE << "ERROR no EmcSnglClusterLocalExtMicrov4 object found" << std::endl;
00477 }
00478 return;
00479 }
00480
00481 float EmcClusterLocalExtMicrov4::get_disp(unsigned int iclus, short i) const
00482 {
00483
00484 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00485 return((emcclus) ? emcclus->get_disp(i) : -9999.9);
00486 }
00487
00488 void EmcClusterLocalExtMicrov4::set_disp(unsigned int iclus, short i, float rval)
00489 {
00490 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00491 if (emcclus)
00492 {
00493 emcclus->set_disp(rval,i);
00494 }
00495 else
00496 {
00497 std::cout << PHWHERE << "ERROR no EmcSnglClusterLocalExtMicrov4 object found" << std::endl;
00498 }
00499 return;
00500 }
00501
00502 float EmcClusterLocalExtMicrov4::get_padisp(unsigned int iclus, short i) const
00503 {
00504
00505 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00506 return((emcclus) ? emcclus->get_padisp(i) : -9999.9);
00507 }
00508
00509 void EmcClusterLocalExtMicrov4::set_padisp(unsigned int iclus, short i, float rval)
00510 {
00511 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00512 if (emcclus)
00513 {
00514 emcclus->set_padisp(rval,i);
00515 }
00516 else
00517 {
00518 std::cout << PHWHERE << "ERROR no EmcSnglClusterLocalExtMicrov4 object found" << std::endl;
00519 }
00520 return;
00521 }
00522
00523 float EmcClusterLocalExtMicrov4::get_xyz(unsigned int iclus, short i) const
00524 {
00525 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00526 return((emcclus) ? emcclus->get_xyz(i) : -9999.9);
00527 }
00528
00529 void EmcClusterLocalExtMicrov4::set_xyz(unsigned int iclus, short i, float rval)
00530 {
00531 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00532 if (emcclus)
00533 {
00534 emcclus->set_xyz(rval,i);
00535 }
00536 else
00537 {
00538 std::cout << PHWHERE << "ERROR no EmcSnglClusterLocalExtMicrov4 object found" << std::endl;
00539 }
00540 return;
00541 }
00542
00543 float EmcClusterLocalExtMicrov4::get_yz_cg(unsigned int iclus, short i) const
00544 {
00545
00546 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00547 return((emcclus) ? emcclus->get_yz_cg(i) : -9999.9);
00548 }
00549
00550 void EmcClusterLocalExtMicrov4::set_yz_cg(unsigned int iclus, short i, float rval)
00551 {
00552 EmcSnglClusterLocalExtMicrov4 *emcclus = (EmcSnglClusterLocalExtMicrov4 *) GetEmcClus()->UncheckedAt(iclus);
00553 if (emcclus)
00554 {
00555 emcclus->set_yz_cg(rval,i);
00556 }
00557 else
00558 {
00559 std::cout << PHWHERE << "ERROR no EmcSnglClusterLocalExtMicrov4 object found" << std::endl;
00560 }
00561 return;
00562 }
00563