Changeset 14913

Timestamp:

02/13/13 10:48:11 (12 years ago)

Author:

tbretz

Message:

Reverting to last revision.

File:

• trunk/Mars/mmuon/MMuonSearchPar.cc (modified) (1 diff)

trunk/Mars/mmuon/MMuonSearchPar.cc

@@ -295 +295 @@
 
     //SetReadyToSave();
-}
-
-Bool_t MMuonSearchPar::CalcFact(const MGeomCam &geom, const MSignalCam &evt)
-{
-    // ===================== Reset cleaning ======================
-
-    for (UInt_t idx=0; idx<evt.GetNumPixels(); idx++)
-    {
-        MSignalPix &pix = evt[idx];
-        if (pix.IsPixelUnmapped())
-            continue;
-
-        pix.SetPixelUnused();
-        pix.SetPixelCore(kFALSE);
-    }
-
-    // ============ Do muon cleaning / calculate COG =============
-
-    // The window must be large enough that the earliest and latest
-    // events do not get a biased timerms
-    double time_min   =  -19.5;
-    double time_max   =  -4.5;
-    double signal_min =  2.30;
-    double delta_t    =  1.75;
-
-    Double_t sumx = 0.;
-    Double_t sumy = 0.;
-    Double_t sumw = 0.;
-    for (UInt_t i=0; i<evt.GetNumPixels(); i++)
-    {
-        MSignalPix &pix = evt[i];
-
-        if (pix.GetNumPhotons()<signal_min)
-            continue;
-        if (pix.GetArrivalTime()>=time_max || pix.GetArrivalTime()<time_min)
-            continue;
-        if (pix.IsPixelUnmapped())
-            continue;
-
-        const MGeom &gpix = geom[i];
-
-        const double x = gpix.GetX();
-        const double y = gpix.GetY();
-
-        int counter = 0;
-        for (int j=0; j<gpix.GetNumNeighbors(); j++)
-        {
-            const int idx = gpix.GetNeighbor(j);
-
-            const MSignalPix &spix = evt[idx];
-
-            if (spix.GetNumPhotons()<signal_min)
-                continue;
-            if (spix.GetArrivalTime()>=time_max || spix.GetArrivalTime()<time_min)
-                continue;
-            if (spix.IsPixelUnmapped())
-                continue;
-
-            if (TMath::Abs(pix.GetArrivalTime()-spix.GetArrivalTime())>delta_t)
-                continue;
-
-            counter++;
-        }
-
-        if (counter==0)
-            continue;
-
-        sumx += pix.GetNumPhotons()*x;
-        sumy += pix.GetNumPhotons()*y;
-        sumw += pix.GetNumPhotons();
-
-        pix.SetPixelUsed();
-        pix.SetPixelCore();
-    }
-
-    if (sumw==0)
-        return kFALSE;
-
-    sumx /= sumw;
-    sumy /= sumw;
-
-    // ==============  Fit circle to resulting data ===============
-
-    Double_t sigma, rad;
-    CalcMinimumDeviation(geom, evt, sumx, sumy, sigma, rad);
-
-    fCenterX   = sumx;
-    fCenterY   = sumy;
-    fRadius    = rad;
-    fDeviation = sigma;
-
-    return kTRUE;
-}
+} 
 
 void MMuonSearchPar::Print(Option_t *) const