diff --git a/src/THcCherenkov.cxx b/src/THcCherenkov.cxx index 179e4cffd6df985d3715a251e48226c7d3c16a92..52eaa7b7f58855336f4f0c95fa1de554babc87bc 100644 --- a/src/THcCherenkov.cxx +++ b/src/THcCherenkov.cxx @@ -272,7 +272,7 @@ void THcCherenkov::Clear(Option_t* opt) // Clear Cherenkov variables from h_trans_cer.f fNhits = 0; // Don't really need to do this. (Be sure this called before Decode) - fNPEsum = 0; + fNPEsum = 0.0; fNCherHit = 0; for(Int_t itube = 0;itube < fNelem;itube++) { @@ -412,15 +412,7 @@ Int_t THcCherenkov::FineProcess( TClonesArray& tracks ) } } - - // if ( fCerEvent > 5880 ) { - // cout << "Event = " << fCerEvent - // << " region = " << ir + 1 - // << " track counter = " << fCerTrackCounter[ir] - // << " fired coutner = " << fCerFiredCounter[ir] - // << endl; - // } - + } // loop over regions // cout << endl; @@ -512,6 +504,7 @@ Int_t THcCherenkov::GetCerIndex( Int_t nRegion, Int_t nValue ) { return fCerNRegions * nValue + nRegion; } + //_____________________________________________________________________________ void THcCherenkov::Print( const Option_t* opt) const { THaNonTrackingDetector::Print(opt); @@ -530,6 +523,12 @@ void THcCherenkov::Print( const Option_t* opt) const { cout << endl; } +//_____________________________________________________________________________ +Double_t THcCherenkov::GetCerNPE() { + + return fNPEsum; +} + ClassImp(THcCherenkov) //////////////////////////////////////////////////////////////////////////////// diff --git a/src/THcCherenkov.h b/src/THcCherenkov.h index 4a0a90b75297151b2e156822bcda9e3f26258077..9d6eee02693bde6c49ffedf1b5baa0f86c6bb76a 100644 --- a/src/THcCherenkov.h +++ b/src/THcCherenkov.h @@ -38,6 +38,9 @@ class THcCherenkov : public THaNonTrackingDetector, public THcHitList { Int_t GetCerIndex(Int_t nRegion, Int_t nValue); + // Double_t GetCerNPE() { return fNPEsum;} + Double_t GetCerNPE(); + THcCherenkov(); // for ROOT I/O protected: Int_t fAnalyzePedestals; @@ -52,8 +55,8 @@ class THcCherenkov : public THaNonTrackingDetector, public THcHitList { Double_t* fADC; // [fNelem] Array of ADC amplitudes Double_t* fADC_P; // [fNelem] Array of ADC amplitudes Double_t* fNPE; // [fNelem] Array of ADC amplitudes - Double_t fNPEsum; // [fNelem] Array of ADC amplitudes - Double_t fNCherHit; // [fNelem] Array of ADC amplitudes + Double_t fNPEsum; + Int_t fNCherHit; Double_t* fCerRegionValue; Double_t fCerChi2Max;