Changeset 19160


Ignore:
Timestamp:
08/14/18 10:20:02 (6 years ago)
Author:
Daniela Dorner
Message:
new disp parametrization for analysis + cut that had not been committed yet
File:
1 edited

Legend:

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

    r18753 r19160  
    4343    // Crab (05h34'32" 22d00'52")
    4444    MPointingPos source("MSourcePos");
    45     //source.SetSkyPosition(5.5755, 22.0144); // ra[h], dec[deg]
     45    //source.SetSkyPosition(5.5755, 22.0144); // ra[h], dec[deg] // Crab
    4646    //source.SetSkyPosition(12.356089, 30.176897); // ra[h], dec[deg]//1218
    4747    //source.SetSkyPosition(16.897867,39.760201);// mrk501
     
    5454
    5555    // FIXME: If fPathIn read cuts and energy estimator from file!
    56     MContinue contq("MImagePar.fNumIslands>3.5"
    57                     "|| MNewImagePar.fNumUsedPixels<5.5"
     56    MContinue contq("MNewImagePar.fNumUsedPixels<5.5"
     57                 // "|| MImagePar.fNumIslands>3.5"
    5858                    "|| MNewImagePar.fLeakage1>0.1",
    5959                    "CutQ");
    6060
    6161    MContinue cont0("MHillas.fSize<0"
    62                     "|| MHillasExt.fSlopeSpreadWeighted>0.68"
    63                     "|| MHillasExt.fSlopeSpreadWeighted<0.18"
    64                     "|| log10(MHillas.GetArea*MGeomCam.fConvMm2Deg^2)>(log10(MHillas.fSize)-2)*1.1-1.55"
    65                     "|| MNewImagePar.fConcCore < 0.13"
    66                     "|| MNewImagePar.fConcCOG  < 0.15"
     62                    "|| MHillas.GetArea > log10(MHillas.fSize)*898-1535 "
    6763                    , "Cut0");
    6864
     
    8278    // Parametrization of Disp
    8379    // for new cleaning
    84     param[0]  =  1.39252;  // constant
    85     param[8]  =  0.154247; // slope
     80    //param[0]  =  1.39252;  // constant
     81    param[0]  =  1.34;  // constant
     82    //param[8]  =  0.154247; // slope
     83    param[8]  =  0.0755; // slope
    8684    param[9]  =  1.67972;  // leak
    8785    param[10] =  4.86232;  // leak
     
    9694    param[7]  =  0.5;
    9795    // ThetaSq-Cut
    98     param[1]  =  0.11;  // 0.215
     96    param[1]  =  0.156;  // 0.215
    9997    // Area-Cut
    10098    param[2]  =  0.215468;
     
    114112    // --------------------------------------------------------------------------------
    115113
     114    //MBinning bins3( 67, -0.005, 0.665, "BinningTheta",  "asin");
    116115    MBinning bins3(100, -0.005, 0.995, "BinningTheta",  "asin");
    117116    MBinning binsD( 25,   0,      2.5, "BinningDist");
     
    466465        path += fname1;
    467466
    468         d->SaveAs(path);
     467        //d->SaveAs(path);
     468        d->Write();
    469469    }
    470470
Note: See TracChangeset for help on using the changeset viewer.