Changeset 3140 for trunk/MagicSoft


Ignore:
Timestamp:
02/13/04 15:28:02 (21 years ago)
Author:
wittek
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars
Files:
2 added
10 edited

Legend:

Unmodified
Added
Removed

  • trunk/MagicSoft/Mars/Changelog

    r3139 r3140  
    44
    55                                                 -*-*- END OF LINE -*-*-
     6
     7
     8 2004/02/13: Wolfgang Wittek
     9
     10  * mcalib/MCalibrateData.[h,cc]
     11    - new class;
     12      like MCalibrate, but in ReInit it fills the MPedPhot container
     13      using informations from MPedestal, MExtracteSignalCam and
     14      MCalibrationCam
     15
     16  * mcalib/Makefile
     17           CalibLinkDef.h
     18    - added MCalibrateData
     19
     20  * manalysis/MPedestalWorkaround.[h,cc]
     21    - create MMcEvt container and fill with theta = 10.0
     22    - set pedestalRMS = 0. if it is > 10.0
     23
     24  * manalysis/MSupercuts.[h,cc]
     25    - change default values of parameters
     26
     27  * manalysis/MSupercutsCalc.[h,cc]
     28    - change dNOMLOGSIZE from 4.1 to 5.0
     29
     30
     31  * mhist/MHBlindPixels.[h,cc]
     32    - change MPedestalCam to MPedPhotCam
     33
     34  * mhist/MHSigmaTheta.[h,cc]
     35    - automatic binning for pixel number doesn't work
     36      use default binning
     37
     38
    639 2004/02/13: Markus Gaug
    740

  • trunk/MagicSoft/Mars/macros/ONOFFAnalysis.C

    r2798 r3140  
    8080
    8181        MBinning *binsdiff = new MBinning("BinningDiffsigma2");
    82         binsdiff->SetEdges(100, -5.0, 20.0);
     82        binsdiff->SetEdges(100, -10.0, 15.0);
    8383        plist->AddToList(binsdiff);
    8484
    … …  
    176176      //const char *onfile  = "~magican/ct1test/wittek/mkn421_on.preproc";
    177177      //      const char *onfile  = "~magican/ct1test/wittek/mkn421_00-01";
    178       const char *onfile  = "/data/MAGIC/rootdata/2003_11_29/20031128_031*_D_Crab-Nebula_E.root";
     178      const char *onfile  = "20040126_Crab_";
    179179
    180180     const char *mcfile  = "/data/MAGIC/mc_eth/magLQE_3/gh/0/0/G_M0_00_0_550*.root";
    … …  
    184184
    185185      // path for input for Mars
    186       TString inPath = "/data/MAGIC/scratch/wittek/";
     186      TString inPath = "/.magic/magicserv01/scratch/";
    187187
    188188      // path for output from Mars
    189       TString outPath = "/data/MAGIC/scratch/wittek/";
     189      TString outPath = "~wittek/datacrab/";
    190190
    191191      //-----------------------------------------------
    … …  
    203203    //  - write root file for ON (or OFF or MC) data (ON1.root, ...);
    204204
    205     Bool_t JobA    = kTRUE; 
     205    Bool_t JobA    = kFALSE; 
    206206    Bool_t GPad    = kFALSE;    // generate padding histograms?
    207207    Bool_t WPad    = kFALSE;   // write out padding histograms ?
    … …  
    234234
    235235    Bool_t JobB_SC_UP  = kTRUE;
    236     Bool_t CMatrix     = kTRUE;  // create training and test matrices
    237     Bool_t RMatrix     = kFALSE;  // read training and test matrices from file
    238     Bool_t WOptimize   = kFALSE;  // do optimization using the training sample
     236    Bool_t CMatrix     = kFALSE;  // create training and test matrices
     237    Bool_t RMatrix     = kTRUE;  // read training and test matrices from file
     238    Bool_t WOptimize   = kTRUE;  // do optimization using the training sample
    239239                                  // and write supercuts parameter values
    240240                                  // onto the file parSCfile
    241     Bool_t RTest       = kFALSE;  // test the supercuts using the test matrix
     241    Bool_t RTest       = kTRUE;  // test the supercuts using the test matrix
    242242    Bool_t WSC         = kFALSE;  // update input root file ?
    243243
    … …  
    312312    // names of ON and OFF files to be read
    313313    // for generating the histograms to be used in the padding
    314     TString fileON  = onfile;
     314
     315    TString fileON  = inPath;
     316    fileON += onfile;
     317    fileON += "CalibratedEvts";
     318    fileON += ".root";
    315319    TString fileOFF = offfile;
    316320    TString fileMC  = mcfile;
    … …  
    333337    // name of input root file
    334338    if (typeInput == "ON")
    335       TString filenamein(onfile);
     339      TString filenamein(fileON);
    336340    else if (typeInput == "OFF")
    337       TString filenamein(offfile);
     341      TString filenamein(fileOFF);
    338342    else if (typeInput == "MC")
    339       TString filenamein(mcfile);
     343      TString filenamein(fileMC);
    340344    gLog << "data to be padded : " << filenamein << endl;
    341345
    342346    // name of output root file
    343347    TString outNameImage = outPath;
     348    outNameImage += onfile;
     349    outNameImage += "Hillas";
    344350    outNameImage += typeInput;
    345351    outNameImage += "1.root";
    … …  
    375381
    376382    MReadMarsFile  readON("Events", fileON);
    377     read.DisableAutoScheme();
    378       MCT1ReadPreProc readON(fileON);
     383    readON.DisableAutoScheme();
     384    //MCT1ReadPreProc readON(fileON);
    379385
    380386      //MFSelBasic selthetaon;
    … …  
    454460
    455461    MReadMarsFile  readOFF("Events", fileOFF);
    456     read.DisableAutoScheme();
     462    readOFF.DisableAutoScheme();
    457463    //      MCT1ReadPreProc readOFF(fileOFF);
    458464
    … …  
    534540
    535541    MReadMarsFile  readMC("Events", fileMC);
    536     read.DisableAutoScheme();
     542    readMC.DisableAutoScheme();
    537543    //      MCT1ReadPreProc readMC(fileMC);
    538544
    … …  
    656662
    657663    MGeomApply        apply;
    658     MMcPedestalCopy   pcopy;
    659     MMcPedestalNSBAdd pnsb;
    660664
    661665    MPedestalWorkaround waround;
    … …  
    667671    f2.SetName("Select Data");
    668672
    669     MCerPhotCalc      ncalc;
    670     ncalc.SetFilter(&f1);
    671     MCerPhotAnal2     nanal;
    672     nanal.SetFilter(&f2);
    673 
    674673    //if (typeInput ==  "ON")
    675674    //{
    … …  
    690689    contbasic.SetName("SelBasic");
    691690
    692     MFillH fillblind("BlindPixels[MHBlindPixels]", "MBlindPixels");
    693     fillblind.SetName("HBlind");
     691    //MFillH fillblind("BlindPixels[MHBlindPixels]", "MBlindPixels");
     692    //fillblind.SetName("HBlind");
    694693
    695694    MSigmabarCalc sigbarcalc;
    … …  
    734733    selstandard.SetHillasName(fHilName);
    735734    selstandard.SetImgParName(fImgParName);
    736     selstandard.SetCuts(92, 4, 60, 0.4, 1.05, 0.0, 0.0);
     735    selstandard.SetCuts(100, 4, 20, 0.0, 1.0, 0.0, 0.0);
    737736    MContinue contstandard(&selstandard);
    738737    contstandard.SetName("SelStandard");
    739738
    740 
     739   
    741740      MWriteRootFile write(outNameImage);
    742741
    743742      write.AddContainer("MRawRunHeader", "RunHeaders");
    744       write.AddContainer("MMcRunHeader",  "RunHeaders", kFALSE);
     743      //write.AddContainer("MMcRunHeader",  "RunHeaders", kFALSE);
    745744      //write.AddContainer("MTime",         "Events");
    746745      write.AddContainer("MMcEvt",        "Events");
    … …  
    752751      write.AddContainer("MHillasSrc",    "Events");
    753752      write.AddContainer("MNewImagePar",  "Events");
    754 
     753   
    755754
    756755    //*****************************
    … …  
    767766    tliston.AddToList(&read);
    768767
    769     tliston.AddToList(&f1);
    770     tliston.AddToList(&f2);
     768    //tliston.AddToList(&f1);
     769    //tliston.AddToList(&f2);
    771770    tliston.AddToList(&apply);
    772     tliston.AddToList(&pcopy);
    773     //tliston.AddToList(&waround);
    774 
    775     tliston.AddToList(&pnsb);
    776     tliston.AddToList(&ncalc);
    777     tliston.AddToList(&nanal);
    778 
    779     tliston.AddToList(&blindbeforepad);
     771    tliston.AddToList(&waround);
     772
     773    //tliston.AddToList(&blindbeforepad);
    780774    //  tliston.AddToList(&pad);
    781775    //  if (typeInput ==  "ON")
    782776    //  tliston.AddToList(&pointcorr);
    783     tliston.AddToList(&blind);
     777
     778    //tliston.AddToList(&blind);
    784779    tliston.AddToList(&contbasic);
    785780
    786     tliston.AddToList(&fillblind);
     781    //tliston.AddToList(&fillblind);
    787782    tliston.AddToList(&sigbarcalc);
    788783    tliston.AddToList(&fillsigtheta);
    … …  
    825820
    826821    pliston.FindObject("SigmaTheta", "MHSigmaTheta")->DrawClone();
    827     pliston.FindObject("BlindPixels", "MHBlindPixels")->DrawClone();
     822    //pliston.FindObject("BlindPixels", "MHBlindPixels")->DrawClone();
    828823
    829824    pliston.FindObject("MHHillas")->DrawClone();
    … …  
    17701765
    17711766    TString parSCinit = outPath;
    1772     //parSCinit += "parSC_1709d";
    1773     parSCinit = "";
    1774 
    1775     gLog << "parSCinit = " << parSCinit << endl;
     1767    parSCinit += "parSC_060204a";
     1768    //parSCinit = "";
     1769
     1770    if (parSCinit != "")
     1771    {
     1772      gLog << "Initial values of parameters are taken from file '"
     1773           << parSCinit << "'" << endl;
     1774    }
     1775    else
     1776    {
     1777      gLog << "Initial values of parameters are taken from the constructor of MSupercuts"
     1778           << endl;
     1779    }
     1780
    17761781
    17771782    //---------------
    … …  
    17801785
    17811786    TString parSCfile = outPath;
    1782     parSCfile += "parSC_2310a";
     1787    parSCfile += "parSC_060204b";
    17831788
    17841789    gLog << "parSCfile = " << parSCfile << endl;
    … …  
    17871792    // file to be updated (either ON or MC)
    17881793
    1789     //TString typeInput = "ON";
     1794    TString typeInput = "ON";
    17901795    //TString typeInput = "OFF";
    1791     TString typeInput = "MC";
     1796    //TString typeInput = "MC";
    17921797    gLog << "typeInput = " << typeInput << endl;
    17931798
    17941799    // name of input root file
    17951800    TString filenameData = outPath;
     1801    filenameData += onfile;
     1802    filenameData += "Hillas";
    17961803    filenameData += typeInput;
    1797     filenameData += "2.root";
     1804    filenameData += "1.root";
    17981805    gLog << "filenameData = " << filenameData << endl;
    17991806
    18001807    // name of output root file
    18011808    TString outNameImage = outPath;
     1809    outNameImage += onfile;
     1810    outNameImage += "Hillas";
    18021811    outNameImage += typeInput;
    1803     outNameImage += "3.root";
     1812    outNameImage += "2.root";
    18041813   
    18051814
    … …  
    18131822    // for the training
    18141823    TString filenameTrain = outPath;
    1815     filenameTrain += "ON";
     1824    filenameTrain += onfile;
     1825    filenameTrain += "Hillas";
     1826    filenameTrain += typeInput;
    18161827    filenameTrain += "1.root";
    1817     Int_t howManyTrain = 800000;
     1828    Int_t howManyTrain = 7000;
    18181829    gLog << "filenameTrain = " << filenameTrain << ",   howManyTrain = "
    18191830         << howManyTrain  << endl;
    … …  
    18211832    // for testing
    18221833    TString filenameTest = outPath;
    1823     filenameTest += "ON";
     1834    filenameTest += onfile;
     1835    filenameTest += "Hillas";
     1836    filenameTest += typeInput;
    18241837    filenameTest += "1.root";
    1825     Int_t howManyTest = 800000;
     1838    Int_t howManyTest = 7000;
    18261839
    18271840    gLog << "filenameTest = " << filenameTest << ",   howManyTest = "
    … …  
    19301943    TArrayD steps(0);
    19311944 
     1945
    19321946    if (parSCinit == "")
    19331947    {
    19341948      Double_t vparams[104] = {
    19351949      // LengthUp
    1936         0.315585,  0.001455, 0.203198, 0.005532, -0.001670, -0.020362,
    1937         0.007388, -0.013463,
     1950        0.2,       0.,       0.,       0.,        0.,        0.0,
     1951        0.,        0.,
    19381952      // LengthLo
    1939         0.151530,  0.028323, 0.510707, 0.053089,  0.013708,  2.357993,
    1940         0.000080, -0.007157,
     1953        0.,        0.,       0.,       0.,        0.,        0.0,
     1954        0.,        0.,
    19411955      // WidthUp
    1942         0.145412, -0.001771, 0.054462, 0.022280, -0.009893,  0.056353,
    1943         0.020711, -0.016703,
     1956        0.1,       0.,       0.,       0.,        0.,        0.0,
     1957        0.,        0.,
    19441958      // WidthLo
    1945         0.089187, -0.006430, 0.074442, 0.003738, -0.004256, -0.014101,
    1946         0.006126, -0.002849,
     1959        0.,        0.,       0.,       0.,        0.,        0.0,
     1960        0.,        0.,
    19471961      // DistUp
    1948         1.787943,  0.0,      2.942310, 0.199815,  0.0,       0.249909,
    1949         0.189697,  0.0,
     1962        0.,        0.,       0.,       0.,        0.,        0.0,
     1963        0.,        0.,
    19501964      // DistLo
    1951         0.589406,  0.0,     -0.083964,-0.007975,  0.0,       0.045374,
    1952        -0.001750,  0.0,
     1965        0.,        0.,       0.,       0.,        0.,        0.0,
     1966        0.,        0.,
    19531967      // AsymUp
    19541968        1.e10,     0.0,      0.0,      0.0,       0.0,       0.0,
    … …  
    19751989      Double_t vsteps[104] = {
    19761990      // LengthUp
    1977         0.03,      0.0002,   0.02,     0.0006,    0.0002,    0.002,
    1978         0.0008,    0.002,
     1991        0.02,      0.,       0.,       0.,        0.,        0.0,
     1992        0.,        0.,
    19791993      // LengthLo
    1980         0.02,      0.003,    0.05,     0.006,     0.002,     0.3,
    1981         0.0001,    0.0008,
     1994        0.,        0.,       0.,       0.,        0.,        0.0,
     1995        0.,        0.,
    19821996      // WidthUp
    1983         0.02,      0.0002,   0.006,    0.003,     0.002,     0.006,
    1984         0.002,     0.002,
     1997        0.01,      0.,       0.,       0.,        0.,        0.0,
     1998        0.,        0.,
    19851999      // WidthLo
    1986         0.009,     0.0007,   0.008,    0.0004,    0.0005,    0.002,
    1987         0.0007,    0.003,
     2000        0.,        0.,       0.,       0.,        0.,        0.0,
     2001        0.,        0.,
    19882002      // DistUp
    1989         0.2,       0.0,      0.3,      0.02,      0.0,       0.03,
    1990         0.02,      0.0
     2003        0.,        0.,       0.,       0.,        0.,        0.0,
     2004        0.,        0.,
    19912005      // DistLo
    1992         0.06,      0.0,      0.009,    0.0008,    0.0,       0.005,
    1993         0.0002,    0.0
    1994       // AsymUp 
     2006        0.,        0.,       0.,       0.,        0.,        0.0,
     2007        0.,        0.,
     2008      // AsymUp
    19952009        0.0,       0.0,      0.0,      0.0,       0.0,       0.0,
    19962010        0.0,       0.0,
    1997       // AsymLo 
     2011      // AsymLo
    19982012        0.0,       0.0,      0.0,      0.0,       0.0,       0.0,
    19992013        0.0,       0.0,
    2000       // ConcUp 
     2014      // ConcUp
    20012015        0.0,       0.0,      0.0,      0.0,       0.0,       0.0,
    20022016        0.0,       0.0,
    2003       // ConcLo 
     2017      // ConcLo
    20042018        0.0,       0.0,      0.0,      0.0,       0.0,       0.0,
    20052019        0.0,       0.0,
    2006       // Leakage1Up 
     2020      // Leakage1Up
    20072021        0.0,       0.0,      0.0,      0.0,       0.0,       0.0,
    20082022        0.0,       0.0,
    2009       // Leakage1Lo 
     2023      // Leakage1Lo
    20102024        0.0,       0.0,      0.0,      0.0,       0.0,       0.0,
    20112025        0.0,       0.0,
    2012       // AlphaUp 
    2013         0.0,       0.0,      0.0,      0.0,       0.0,       0.0,
     2026      // AlphaUp
     2027        0.0,       0.0,      0.0,      0.0,       0.0,       0.0,
    20142028        0.0,       0.0                                                 };
     2029
    20152030
    20162031      params.Set(104, vparams);
    … …  
    21422157
    21432158      write.AddContainer("MRawRunHeader", "RunHeaders");
    2144       write.AddContainer("MTime",         "Events");
     2159      //write.AddContainer("MTime",         "Events");
    21452160      write.AddContainer("MMcEvt",        "Events");
    2146       write.AddContainer("ThetaOrig",     "Events");
     2161      //write.AddContainer("ThetaOrig",     "Events");
    21472162      write.AddContainer("MSrcPosCam",    "Events");
    21482163      write.AddContainer("MSigmabar",     "Events");
    … …  
    21522167      write.AddContainer("MNewImagePar",  "Events");
    21532168
    2154       write.AddContainer("HadRF",         "Events");
     2169      //write.AddContainer("HadRF",         "Events");
    21552170      write.AddContainer(hadSCName,       "Events");
    21562171   

  • trunk/MagicSoft/Mars/manalysis/MPedestalWorkaround.cc

    r2746 r3140  
    5151#include "MGeomCam.h"
    5252
    53 #include "MPedestalCam.h"
    54 #include "MPedestalPix.h"
    5553#include "MPedPhotCam.h"
     54#include "MPedPhotPix.h"
    5655#include "MMcEvt.hxx"
    5756#include "MRawRunHeader.h"
    … …  
    7372Int_t MPedestalWorkaround::PreProcess(MParList *pList)
    7473{
     74  /*
    7575   fPed = (MPedestalCam*)pList->FindObject("MPedestalCam");
    7676   if (!fPed)
    … …  
    7979       return kFALSE;
    8080     }
     81  */
    8182
    8283   fPedPhot = (MPedPhotCam*)pList->FindObject("MPedPhotCam");
    … …  
    101102   }
    102103
    103    fMcEvt = (MMcEvt*)pList->FindObject("MMcEvt");
     104   fMcEvt = (MMcEvt*)pList->FindCreateObj("MMcEvt");
    104105   if (!fMcEvt)
    105106   {
    … …  
    119120  // from MPedPhotCam into MPedestalCam   
    120121
     122  // set pedestalRMS = 0 if it is too high
    121123  UInt_t imaxnumpix = fCam->GetNumPixels();
    122 
    123124  for (UInt_t i=0; i<imaxnumpix; i++)
    124125  {
    125     Int_t type = 0;
    126     Double_t val;
    127     Float_t valout;
    128     fPedPhot->GetPixelContent( val, i, *fCam, type);
    129     valout = (*fPed)[i].GetPedestal();
    130     (*fPed)[i].SetPedestal(val);
    131     *fLog << "i, val, valout : " << i <<",  "<<  val<<",  " << valout << endl;
    132     type = 1;
    133     fPedPhot->GetPixelContent( val, i, *fCam, type);
    134     valout = (*fPed)[i].GetPedestalRms();
    135     (*fPed)[i].SetPedestalRms(val);
    136     *fLog << "RMS : i, val, valout : " << i <<",  "<<  val<<",  " << valout
    137           << endl;
     126    Double_t val = (*fPedPhot)[i].GetRms();
     127    if (val < 0.0 || val > 10.0)
     128      (*fPedPhot)[i].SetRms(0.0);
    138129  }
    139130
    … …  
    142133  // put the zenith angle into MMcEvt
    143134
    144   Int_t run = fRun->GetRunNumber();
    145135  Double_t thetadeg;
    146136  Double_t thetarad;
     137
     138
     139  /*
     140  Int_t run = fRun->GetRunNumber();
    147141
    148142  if      (run == 3127)  thetadeg = 27.2;  // Crab
    … …  
    537531
    538532  thetarad = thetadeg / kRad2Deg;
     533  */
     534
     535  thetadeg = 10.0;
     536  thetarad = thetadeg / kRad2Deg;
    539537  fMcEvt->SetTelescopeTheta(thetarad);
    540538 

  • trunk/MagicSoft/Mars/manalysis/MSupercuts.cc

    r2663 r3140  
    7474{
    7575    //---------------------------------------------------
    76     //  these are Daniel's original values for Mkn 421
    77 
    78     fLengthUp[0] =  0.315585;
    79     fLengthUp[1] =  0.001455;
    80     fLengthUp[2] =  0.203198;
    81     fLengthUp[3] =  0.005532;
    82     fLengthUp[4] = -0.001670;
    83     fLengthUp[5] = -0.020362;
    84     fLengthUp[6] =  0.007388;
    85     fLengthUp[7] = -0.013463;
    86 
    87     fLengthLo[0] =  0.151530;
    88     fLengthLo[1] =  0.028323;
    89     fLengthLo[2] =  0.510707;
    90     fLengthLo[3] =  0.053089;
    91     fLengthLo[4] =  0.013708;
    92     fLengthLo[5] =  2.357993;
    93     fLengthLo[6] =  0.000080;
    94     fLengthLo[7] = -0.007157;
    95 
    96     fWidthUp[0] =  0.145412;
    97     fWidthUp[1] = -0.001771;
    98     fWidthUp[2] =  0.054462;
    99     fWidthUp[3] =  0.022280;
    100     fWidthUp[4] = -0.009893;
    101     fWidthUp[5] =  0.056353;
    102     fWidthUp[6] =  0.020711;
    103     fWidthUp[7] = -0.016703;
    104 
    105     fWidthLo[0] =  0.089187;
    106     fWidthLo[1] = -0.006430;
    107     fWidthLo[2] =  0.074442;
    108     fWidthLo[3] =  0.003738;
    109     fWidthLo[4] = -0.004256;
    110     fWidthLo[5] = -0.014101;
    111     fWidthLo[6] =  0.006126;
    112     fWidthLo[7] = -0.002849;
    113 
    114     fDistUp[0] =  1.787943;
    115     fDistUp[1] =  0;
    116     fDistUp[2] =  2.942310;
    117     fDistUp[3] =  0.199815;
    118     fDistUp[4] =  0;
    119     fDistUp[5] =  0.249909;
    120     fDistUp[6] =  0.189697;
    121     fDistUp[7] =  0;
    122 
    123     fDistLo[0] =  0.589406;
    124     fDistLo[1] =  0;
    125     fDistLo[2] = -0.083964;
    126     fDistLo[3] = -0.007975;
    127     fDistLo[4] =  0;
    128     fDistLo[5] =  0.045374;
    129     fDistLo[6] = -0.001750;
    130     fDistLo[7] =  0;
     76    //  these are the default values
     77
     78    fLengthUp[0] =  0.2;
     79    fLengthUp[1] =  0.0;
     80    fLengthUp[2] =  0.0;
     81    fLengthUp[3] =  0.0;
     82    fLengthUp[4] =  0.0;
     83    fLengthUp[5] =  0.0;
     84    fLengthUp[6] =  0.0;
     85    fLengthUp[7] =  0.0;
     86
     87    fLengthLo[0] =  0.;
     88    fLengthLo[1] =  0.;
     89    fLengthLo[2] =  0.;
     90    fLengthLo[3] =  0.;
     91    fLengthLo[4] =  0.;
     92    fLengthLo[5] =  0.;
     93    fLengthLo[6] =  0.;
     94    fLengthLo[7] =  0.;
     95
     96    fWidthUp[0] =  0.1;
     97    fWidthUp[1] =  0.0;
     98    fWidthUp[2] =  0.0;
     99    fWidthUp[3] =  0.0;
     100    fWidthUp[4] =  0.0;
     101    fWidthUp[5] =  0.0;
     102    fWidthUp[6] =  0.0;
     103    fWidthUp[7] =  0.0;
     104
     105    fWidthLo[0] =  0.;
     106    fWidthLo[1] =  0.;
     107    fWidthLo[2] =  0.;
     108    fWidthLo[3] =  0.;
     109    fWidthLo[4] =  0.;
     110    fWidthLo[5] =  0.;
     111    fWidthLo[6] =  0.;
     112    fWidthLo[7] =  0.;
     113
     114    fDistUp[0] =  1.e10;
     115    fDistUp[1] =  0.0;
     116    fDistUp[2] =  0.0;
     117    fDistUp[3] =  0.0;
     118    fDistUp[4] =  0.0;
     119    fDistUp[5] =  0.0;
     120    fDistUp[6] =  0.0;
     121    fDistUp[7] =  0.0;
     122
     123    fDistLo[0] =  0.0;
     124    fDistLo[1] =  0.0;
     125    fDistLo[2] =  0.0;
     126    fDistLo[3] =  0.0;
     127    fDistLo[4] =  0.0;
     128    fDistLo[5] =  0.0;
     129    fDistLo[6] =  0.0;
     130    fDistLo[7] =  0.0;
    131131   
    132132
    … …  
    202202
    203203    // LengthUp
    204     fStepsizes[0] = 0.03;
    205     fStepsizes[1] = 0.0002;
    206     fStepsizes[2] = 0.02;
    207     fStepsizes[3] = 0.0006;
    208     fStepsizes[4] = 0.0002;
    209     fStepsizes[5] = 0.002;
    210     fStepsizes[6] = 0.0008;
    211     fStepsizes[7] = 0.002;
     204    fStepsizes[0] = 0.02;
     205    fStepsizes[1] = 0.0;
     206    fStepsizes[2] = 0.0;
     207    fStepsizes[3] = 0.0;
     208    fStepsizes[4] = 0.0;
     209    fStepsizes[5] = 0.0;
     210    fStepsizes[6] = 0.0;
     211    fStepsizes[7] = 0.0;
    212212
    213213    // LengthLo
    214     fStepsizes[8]  = 0.02;
    215     fStepsizes[9]  = 0.003;
    216     fStepsizes[10] = 0.05;
    217     fStepsizes[11] = 0.006;
    218     fStepsizes[12] = 0.002;
    219     fStepsizes[13] = 0.3;
    220     fStepsizes[14] = 0.0001;
    221     fStepsizes[15] = 0.0008;
     214    fStepsizes[8]  = 0.0;
     215    fStepsizes[9]  = 0.0;
     216    fStepsizes[10] = 0.0;
     217    fStepsizes[11] = 0.0;
     218    fStepsizes[12] = 0.0;
     219    fStepsizes[13] = 0.0;
     220    fStepsizes[14] = 0.0;
     221    fStepsizes[15] = 0.0;
    222222
    223223    // WidthUp
    224     fStepsizes[16] = 0.02;
    225     fStepsizes[17] = 0.0002;
    226     fStepsizes[18] = 0.006;
    227     fStepsizes[19] = 0.003;
    228     fStepsizes[20] = 0.002;
    229     fStepsizes[21] = 0.006;
    230     fStepsizes[22] = 0.002;
    231     fStepsizes[23] = 0.002;
     224    fStepsizes[16] = 0.01;
     225    fStepsizes[17] = 0.0;
     226    fStepsizes[18] = 0.0;
     227    fStepsizes[19] = 0.0;
     228    fStepsizes[20] = 0.0;
     229    fStepsizes[21] = 0.0;
     230    fStepsizes[22] = 0.0;
     231    fStepsizes[23] = 0.0;
    232232
    233233    // WidthLo
    234     fStepsizes[24] = 0.009;
    235     fStepsizes[25] = 0.0007;
    236     fStepsizes[26] = 0.008;
    237     fStepsizes[27] = 0.0004;
    238     fStepsizes[28] = 0.0005;
    239     fStepsizes[29] = 0.002;
    240     fStepsizes[30] = 0.0007;
    241     fStepsizes[31] = 0.003;
     234    fStepsizes[24] = 0.0;
     235    fStepsizes[25] = 0.0;
     236    fStepsizes[26] = 0.0;
     237    fStepsizes[27] = 0.0;
     238    fStepsizes[28] = 0.0;
     239    fStepsizes[29] = 0.0;
     240    fStepsizes[30] = 0.0;
     241    fStepsizes[31] = 0.0;
    242242
    243243    // DistUp
    244     fStepsizes[32] = 0.2;
     244    fStepsizes[32] = 0.0;
    245245    fStepsizes[33] = 0.0;
    246     fStepsizes[34] = 0.3;
    247     fStepsizes[35] = 0.02;
     246    fStepsizes[34] = 0.0;
     247    fStepsizes[35] = 0.0;
    248248    fStepsizes[36] = 0.0;
    249     fStepsizes[37] = 0.03;
    250     fStepsizes[38] = 0.02;
     249    fStepsizes[37] = 0.0;
     250    fStepsizes[38] = 0.0;
    251251    fStepsizes[39] = 0.0;
    252252
    253253    // DistLo
    254     fStepsizes[40] = 0.06;
     254    fStepsizes[40] = 0.0;
    255255    fStepsizes[41] = 0.0;
    256     fStepsizes[42] = 0.009;
    257     fStepsizes[43] = 0.0008;
     256    fStepsizes[42] = 0.0;
     257    fStepsizes[43] = 0.0;
    258258    fStepsizes[44] = 0.0;
    259     fStepsizes[45] = 0.005;
    260     fStepsizes[46] = 0.0002;
     259    fStepsizes[45] = 0.0;
     260    fStepsizes[46] = 0.0;
    261261    fStepsizes[47] = 0.0;
    262262

  • trunk/MagicSoft/Mars/manalysis/MSupercutsCalc.cc

    r2663 r3140  
    157157    // define cut-function
    158158    //
    159     //    dNOMLOGSIZE = 4.1 (=log(60.0)
     159    //    dNOMLOGSIZE = 5.0 (=log(150.0)
    160160    //    dNOMCOSZA   = 1.0
    161161    //

  • trunk/MagicSoft/Mars/mcalib/CalibLinkDef.h

    r3108 r3140  
    66
    77#pragma link C++ class MCalibrate+;
     8#pragma link C++ class MCalibrateData+;
    89
    910#pragma link C++ class MCalibrationPix+;

  • trunk/MagicSoft/Mars/mcalib/Makefile

    r3108 r3140  
    3737
    3838SRCFILES = MCalibrate.cc \
     39           MCalibrateData.cc \
    3940           MCalibrationCalc.cc \
    4041           MCalibrationPix.cc  \

  • trunk/MagicSoft/Mars/mhist/MHBlindPixels.cc

    r2173 r3140  
    3535#include "MMcEvt.hxx"
    3636#include "MBlindPixels.h"
    37 #include "MPedestalCam.h"
     37#include "MPedPhotCam.h"
    3838#include "MParList.h"
    3939#include "MBinning.h"
    … …  
    8585    }
    8686
    87     fPed = (MPedestalCam*)plist->FindObject("MPedestalCam");
    88     if (!fPed)
     87   
     88    fPedPhot = (MPedPhotCam*)plist->FindObject("MPedPhotCam");
     89    if (!fPedPhot)
    8990    {
    90         *fLog << err << "MPedestalCam not found... aborting." << endl;
     91        *fLog << err << "MPedPhotCam not found... aborting." << endl;
    9192        return kFALSE;
    9293    }
    93 
     94   
    9495
    9596    // Get Theta Binning
    … …  
    102103
    103104    // Get binning for pixel number
    104     const UInt_t npix1 = fPed->GetSize()+1;
     105    const UInt_t npix1 = fPedPhot->GetSize()+1;
    105106
    106107    MBinning binspix("BinningPixel");
    … …  
    166167
    167168    // FIXME: Slow.
    168     const UInt_t npix = fPed->GetSize();
     169    const UInt_t npix = fPedPhot->GetSize();
    169170
    170171    UInt_t nb = 0;

  • trunk/MagicSoft/Mars/mhist/MHBlindPixels.h

    r2128 r3140  
    99#endif
    1010
    11 class MPedestalCam;
     11class MPedPhotCam;
    1212class MMcEvt;
    1313class MParList;
    … …  
    1717{
    1818private:
    19     MPedestalCam  *fPed;      //!
     19    MPedPhotCam   *fPedPhot;      //!
    2020    MMcEvt        *fMcEvt;    //!
    2121

  • trunk/MagicSoft/Mars/mhist/MHSigmaTheta.cc

    r2798 r3140  
    9595    MBinning binst;
    9696    MBinning binsd;
    97     binsd.SetEdges(100, -5, 20);
    98     binsb.SetEdges(100, 0, 5);
     97    binsd.SetEdges(100, -10, 20);
     98    binsb.SetEdges(100, 0, 10);
    9999    binst.SetEdgesCos(10, 0, 90);
    100100
    … …  
    166166
    167167    // Get binning for pixel number
    168     const UInt_t npix1 = fPed->GetSize()+1;
    169 
     168    //const UInt_t npix1 = fPed->GetSize()+1;
     169    //*fLog << "npix1 = " << npix1 << endl;
     170    //MBinning binspix("BinningPixel");
     171    //binspix.SetEdges(npix1, -0.5, npix1-0.5);
     172    // the above gives npix1 = 1; therefore :
    170173    MBinning binspix("BinningPixel");
    171     binspix.SetEdges(npix1, -0.5, npix1-0.5);
     174    binspix.SetEdges(578, -0.5, 577.5);
    172175
    173176    // Set binnings in histograms
    … …  
    191194{
    192195    Double_t theta = fMcEvt ? fMcEvt->GetTelescopeTheta()*kRad2Deg : 0;
    193     Double_t mysig = fSigmabar->Calc(*fCam, *fPed, *fEvt);
     196    fSigmabar->Calc(*fCam, *fPed, *fEvt);
     197    Double_t mysig = fSigmabar->GetSigmabarInner();
     198
     199    //*fLog << "theta, mysig = " << theta << ",  " << mysig << endl;
    194200
    195201    fSigmaTheta.Fill(theta, mysig);
    … …  
    211217        const Double_t sigma = pix.GetRms();
    212218
    213         fSigmaPixTheta.Fill(theta, (Double_t)id, sigma);
     219        fSigmaPixTheta.Fill(theta, (Double_t)id, sigma*sqrt(ratio));
    214220
    215221        const Double_t diff = sigma*sigma*ratio - mysig*mysig;
Note: See TracChangeset for help on using the changeset viewer.