00001
00002
00003
00004
00005
00007
00008 #include <iostream>
00009 using namespace std;
00010
00011 #include "MessageService/MsgService.h"
00012 #include "CandNtupleEM/NtpEMShower.h"
00013
00014 ClassImp(NtpEMShower)
00015
00016 std::ostream& operator << (std::ostream& os, const NtpEMShower& ntp)
00017 { return ntp.Print(os); }
00018
00019 CVSID("$Id: NtpEMShower.cxx,v 1.2 2004/08/25 05:30:32 cbs Exp $");
00020
00021 void NtpEMShower::AddStripAt(Int_t stripindex, Int_t ind) {
00022
00023 if ( ind < 0 || ind >= nstrip ) {
00024 MSG("NtpEM",Msg::kWarning)
00025 << "Attempt to add strip at index " << ind
00026 << " outside of array size " << nstrip << " ignored." << endl;
00027 return;
00028 }
00029
00030 stp[ind] = stripindex;
00031 return;
00032
00033 }
00034
00035 NtpEMShower::NtpEMShower(): index(0),ndigit(0),
00036 nstrip(0),stp(0),
00037 status(-1) {
00038
00039 evector[0] = 0; evector[2] = 0; evector[3] = 0; evector[3] = 0;
00040 evector[4] = 0; evector[5] = 0; evector[6] = 0; evector[7] = 0;
00041 evalue[0] = 0; evalue[1] = 0; evalue[2] = 0; evalue[3] = 0;
00042 avgdev[0] = 0; avgdev[1] = 0; avgdev[2] = 0; avgdev[3] = 0;
00043 outph[0] = 0; outph[1] = 0; outph[2] = 0; outph[3] = 0; outph[4] = 0;
00044
00045 }
00046
00047 NtpEMShower::NtpEMShower(Int_t ninit):index(0),ndigit(0),
00048 nstrip(ninit),stp(0),
00049 status(-1) {
00050
00051 if ( nstrip ) {
00052 stp = new Int_t[nstrip];
00053 for ( Int_t i = 0; i < nstrip; i++ ) stp[i] = -1;
00054 }
00055
00056 evector[0] = 0; evector[2] = 0; evector[3] = 0; evector[3] = 0;
00057 evector[4] = 0; evector[5] = 0; evector[6] = 0; evector[7] = 0;
00058 evalue[0] = 0; evalue[1] = 0; evalue[2] = 0; evalue[3] = 0;
00059 avgdev[0] = 0; avgdev[1] = 0; avgdev[2] = 0; avgdev[3] = 0;
00060 outph[0] = 0; outph[1] = 0; outph[2] = 0; outph[3] = 0; outph[4] = 0;
00061
00062 }
00063
00064 std::ostream& NtpEMShower::Print(std::ostream& os) const {
00065
00066 os << "NtpEMShower::Print" << endl;
00067 return os;
00068
00069 }
00070
00071 void NtpEMShower::Print(Option_t* ) const {
00072
00073 Print(std::cout);
00074 return;
00075
00076 }