Ignore:
Timestamp:
05/26/14 16:29:04 (10 years ago)
Author:
Daniela Dorner
Message:
using new image cleaning developed by Thomas
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/Mars/fact/analysis/star.C

    r17885 r17895  
    11#include "MLogManip.h"
    22
    3 int star(const char *seqfile="sequences/20111205_013.seq", Double_t lvl1=4.0, Double_t lvl2=2.5, const char *inpath = "output", const char *outpath = "output")
     3int star(const char *seqfile="sequences/20111205_013.seq",
     4         const char *inpath = "callisto_new", const char *outpath = "callisto_new")
    45{
    56    double deltat = 17.5;
    6     //    lvl1 = 2.5;
    7     //    lvl2 = 0.5;
    87
    98    // The sequence file which defines the files for the analysis
     
    109108
    110109    MContinue cont("MRawEvtHeader.GetTriggerID!=4", "SelectData");
     110    //MContinue cont("(MRawEvtHeader.GetTriggerID&0xff00)!=0",     "SelectDat");
     111    //MContinue cont("(MRawEvtHeader.GetTriggerID&0xff00)!=0x400", "SelectPed");
    111112
    112113    MGeomApply apply;
    113114
    114     MImgCleanStd clean(lvl1, lvl2);
    115     clean.SetMethod(MImgCleanStd::kAbsolute);
    116     clean.SetTimeLvl1(deltat);//17.5/*1.3*/);
    117     clean.SetTimeLvl2(deltat);//17.5/*1.3*/);
    118     clean.SetPostCleanType(3);
     115    MImgCleanTime clean;
     116    clean.SetMinCount(0);
     117    clean.SetMinSize(25);
     118    clean.SetDeltaT(1.25*17.5*0.1111);
    119119
    120120    MHillasCalc hcalc;
     
    296296    return 0;
    297297}
    298 
    299 int star(const ULong64_t seqnum, Double_t lvl1=4.0, Double_t lvl2=2.5, const char *inpath = "output", const char *outpath = "output")
    300 {
    301     UInt_t night = seqnum/1000;
    302     UInt_t num   = seqnum%1000;
    303 
    304     TString file = Form("/scratch/fact/sequences/%04d/%02d/%02d/%06d_%03d.seq",
    305                         night/10000, (night/100)%100, night%100, num);
    306 
    307     return star(file.Data(), lvl1, lvl2, inpath, outpath);
    308 }
Note: See TracChangeset for help on using the changeset viewer.