EmcClusterLocalExtMicrov4.C

Go to the documentation of this file.
00001 #include "EmcClusterLocalExtMicrov4.h"
00002 #include "EmcSnglClusterLocalExtMicrov4.h"
00003 //INCLUDECHECKER: Removed this line: #include "phool.h"
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   // i==0 corresponds to the X direction in the local EMC sector frame
00152   // (which coincides with the Z direction in the global PHENIX frame)
00153   // i==1 corresponds to the Y direction in the local EMC sector fram
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