Skip to content
Snippets Groups Projects
Commit ff106a07 authored by hallc-online's avatar hallc-online Committed by Mark K Jones
Browse files

Updated THcHodoscope.cxx

Modify so default starttime is -1000.

Modify so if no hits in scintillator then EstimateFocalPlane is not called.
parent 04242df6
No related branches found
No related tags found
No related merge requests found
...@@ -163,6 +163,7 @@ THaAnalysisObject::EStatus THcHodoscope::Init( const TDatime& date ) ...@@ -163,6 +163,7 @@ THaAnalysisObject::EStatus THcHodoscope::Init( const TDatime& date )
// maximum number of hits after setting up the detector map // maximum number of hits after setting up the detector map
// But it needs to happen before the sub detectors are initialized // But it needs to happen before the sub detectors are initialized
// so that they can get the pointer to the hitlist. // so that they can get the pointer to the hitlist.
cout << " Hodo tdc ref time cut = " << fTDC_RefTimeCut << " " << fADC_RefTimeCut << endl;
InitHitList(fDetMap, "THcRawHodoHit", fDetMap->GetTotNumChan()+1, InitHitList(fDetMap, "THcRawHodoHit", fDetMap->GetTotNumChan()+1,
fTDC_RefTimeCut, fADC_RefTimeCut); fTDC_RefTimeCut, fADC_RefTimeCut);
...@@ -645,6 +646,7 @@ Int_t THcHodoscope::Decode( const THaEvData& evdata ) ...@@ -645,6 +646,7 @@ Int_t THcHodoscope::Decode( const THaEvData& evdata )
Int_t nexthit = 0; Int_t nexthit = 0;
fNfptimes=0; fNfptimes=0;
Int_t thits=0;
for(Int_t ip=0;ip<fNPlanes;ip++) { for(Int_t ip=0;ip<fNPlanes;ip++) {
fPlaneCenter[ip] = fPlanes[ip]->GetPosCenter(0) + fPlanes[ip]->GetPosOffset(); fPlaneCenter[ip] = fPlanes[ip]->GetPosCenter(0) + fPlanes[ip]->GetPosOffset();
...@@ -654,11 +656,12 @@ Int_t THcHodoscope::Decode( const THaEvData& evdata ) ...@@ -654,11 +656,12 @@ Int_t THcHodoscope::Decode( const THaEvData& evdata )
// GN: select only events that have reasonable TDC values to start with // GN: select only events that have reasonable TDC values to start with
// as per the Engine h_strip_scin.f // as per the Engine h_strip_scin.f
nexthit = fPlanes[ip]->ProcessHits(fRawHitList,nexthit); nexthit = fPlanes[ip]->ProcessHits(fRawHitList,nexthit);
thits+=fPlanes[ip]->GetNScinHits();
} }
EstimateFocalPlaneTime(); fStartTime=-1000;
if (thits>0 ) EstimateFocalPlaneTime();
if (fdebugprintscinraw == 1) { if (fdebugprintscinraw == 1) {
cout << " Event number = " << evdata.GetEvNum()<<endl;
for(UInt_t ihit = 0; ihit < fNRawHits ; ihit++) { for(UInt_t ihit = 0; ihit < fNRawHits ; ihit++) {
// THcRawHodoHit* hit = (THcRawHodoHit *) fRawHitList->At(ihit); // THcRawHodoHit* hit = (THcRawHodoHit *) fRawHitList->At(ihit);
// cout << ihit << " : " << hit->fPlane << ":" << hit->fCounter << " : " // cout << ihit << " : " << hit->fPlane << ":" << hit->fCounter << " : "
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment