diff --git a/src/THcHodoEff.cxx b/src/THcHodoEff.cxx index 2cbc877574afeb02a21016c2e5964a12fd13324f..c2617624fc42d6fd7f0a18a32ec981da2fd64e9a 100644 --- a/src/THcHodoEff.cxx +++ b/src/THcHodoEff.cxx @@ -183,12 +183,12 @@ Int_t THcHodoEff::ReadDatabase( const TDatime& date ) DBRequest list[]={ {"stat_slop", &fStatSlop, kDouble}, {"stat_maxchisq",&fMaxChisq, kDouble}, + {"HodoEff_CalEnergy_Cut",&fHodoEff_CalEnergy_Cut, kDouble,0,1}, {"hodo_slop", fHodoSlop, kDouble, (UInt_t)fNPlanes}, {0} }; - // fMaxShTrk = 0.05; // For cut on fraction of momentum seen in shower + fHodoEff_CalEnergy_Cut=0.050; // set default value gHcParms->LoadParmValues((DBRequest*)&list,prefix); - cout << "\n\nTHcHodoEff::ReadDatabase nplanes=" << fHod->GetNPlanes() << endl; // Setup statistics arrays // Better method to put this in? @@ -353,7 +353,7 @@ Int_t THcHodoEff::Process( const THaEvData& evdata ) Double_t dist = hitDistance[ip]; if(TMath::Abs(dist) <= fStatSlop && theTrack->GetChi2()/theTrack->GetNDoF() <= fMaxChisq && - theTrack->GetEnergy() >= 0.05 ) + theTrack->GetEnergy() >= fHodoEff_CalEnergy_Cut ) { fStatTrk[fHod->GetScinIndex(ip,hitCounter[ip]-1)]++; // Double_t delta = theTrack->GetDp(); @@ -393,7 +393,7 @@ Int_t THcHodoEff::Process( const THaEvData& evdata ) TMath::Abs(hitcounter-counter) <= checkHit[ip] && fHitPlane[ip] == 0 && theTrack->GetChi2()/theTrack->GetNDoF() <= fMaxChisq && - theTrack->GetEnergy() >= 0.05 ) { + theTrack->GetEnergy() >= fHodoEff_CalEnergy_Cut) { fHitPlane[ip]++; // Need to find out hgood_tdc_pos(igoldentrack,ihit) and neg diff --git a/src/THcHodoEff.h b/src/THcHodoEff.h index 47ab1eee57032bceb3d646e5de37d1afb8237f43..d2369c06389bcbac8f30898481e30b6bc85c8224 100644 --- a/src/THcHodoEff.h +++ b/src/THcHodoEff.h @@ -71,6 +71,7 @@ protected: Int_t* fStatAndSum; Double_t* fStatAndEff; Double_t fStatSlop; + Double_t fHodoEff_CalEnergy_Cut; Double_t fMaxChisq; Double_t* fHodoSlop; Double_t fHodoEff_s1,fHodoEff_s2,fHodoEff_tof,fHodoEff_3_of_4,fHodoEff_4_of_4;