Ignore:
Timestamp:
04/01/03 17:49:19 (21 years ago)
Author:
moralejo
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/manalysis/MEnergyEstParam.cc

    r1861 r1886  
    5858    fA[3] =  0.006;
    5959    fA[4] =   58.3;
     60
    6061    /*
    61      fA[0] =   39.667402; // [cm]
    62      fA[1] =  143.081912; // [cm]
    63      fA[2] = -395.501677; // [cm]
    64      fA[3] =    0.006193;
    65      fA[4] =    0;
    66      */
     62      fA[0] =   39.667402; // [cm]
     63      fA[1] =  143.081912; // [cm]
     64      fA[2] = -395.501677; // [cm]
     65      fA[3] =    0.006193;
     66      fA[4] =    0;
     67    */
    6768
    6869    fB.Set(7);
     
    7475    fB[5] =  4.13e-8;   // [GeV]
    7576    fB[6] =     0.05;
     77
    7678    /*
    77      fB[0] =   45.037867; // [GeV]
    78      fB[1] =    0.087222; // [GeV]
    79      fB[2] =   -0.208065; // [GeV/cm]
    80      fB[3] =   78.164942; // [GeV]
    81      fB[4] = -159.361283; // [GeV]
    82      fB[5] =   -0.130455; // [GeV/cm]
    83      fB[6] =    0.051139;
    84      */
     79      fB[0] =   45.037867; // [GeV]
     80      fB[1] =    0.087222; // [GeV]
     81      fB[2] =   -0.208065; // [GeV/cm]
     82      fB[3] =   78.164942; // [GeV]
     83      fB[4] = -159.361283; // [GeV]
     84      fB[5] =   -0.130455; // [GeV/cm]
     85      fB[6] =    0.051139;
     86    */
    8587}
    8688
     
    340342        const Double_t dist = fMatrix ? GetVal(col++) : ((MHillasSrc*)NextH())->GetDist();
    341343
    342 
    343344        /* DANIEL - MARCOS */
    344345        const Double_t ir = i0 * (i1 + fA[1]*dist); // [cm]
    345346        Double_t er = e0 * (e1 + e2*ir);       // [GeV]
    346347
    347 
    348         /* MY PARAM */
     348        /* THOMAS BRETZ */
    349349        // if (width==0) return kCONTINUE;
    350350        // const Double_t ir = i0 * (i1 + dist*(fA[1]/width + fA[4]/log10(size))); // [cm]
Note: See TracChangeset for help on using the changeset viewer.