Ignore:
Timestamp:
08/08/02 14:01:38 (23 years ago)
Author:
tbretz
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars/macros
Files:
3 edited

Legend:

Unmodified
Added
Removed

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

    r1281 r1490  
    1616!
    1717!
    18 !   Author(s): Thomas Bretz  12/2000 (tbretz@uni-sw.gwdg.de)
     18!   Author(s): Thomas Bretz et al,  12/2000 <mailto:tbretz@astro.uni-wuerzburg.de>
    1919!
    20 !   Copyright: MAGIC Software Development, 2000-2001
     20!   Copyright: MAGIC Software Development, 2000-2002
    2121!
    2222!
    … …  
    2424
    2525
    26 void CT1Hillas()
     26void CT1Hillas(const char *filename)
    2727{
    2828    //
    2929    // This is a demonstration program which calculates the Hillas
    30     // parameter out of a CT1 dataset.
     30    // parameter out of a Magic root file (raw data file).
     31    //
    3132
    3233    //
    … …  
    3435    // The tasklist is identified in the eventloop by its name
    3536    //
    36     MParList  plist;
     37    MParList plist;
     38
    3739
    3840    MTaskList tlist;
    … …  
    4749
    4850    //
    49     // Setup binning for your histograms.
     51    // Use this if you want to change the binning of one of
     52    // the histograms. You can use:
     53    // BinningConc, BinningConc1, BinningAsym, BinningM3Long,
     54    // BinningM3Trans, BinningWidth, BinningLength, BinningDist,
     55    // BinningHeadTail, BinningAlpha, BinningSize, BinningDelta,
     56    // BinningPixels and BinningCamera
    5057    //
    51     MBinning binswidth("BinningWidth");
    52     binswidth.SetEdges(100, 0, 1);   // 100 bins from 0 to 1 deg
     58    // For more information see MBinning and the corresponding
     59    // histograms
     60    //
     61    // MBinning binsalpha("BinningAlpha");
     62    // binsalpha.SetEdges(90, 0, 90);       // 90 bins from 0 to 90 deg
     63    // plist.AddToList(&binsalpha);
    5364
    54     MBinning binslength("BinningLength");
    55     binslength.SetEdges(100, 0, 1);  // 100 bins from 0 to 1 deg
    56 
    57     MBinning binsalpha("BinningAlpha");
    58     binsalpha.SetEdges(90, 0, 90);   // 90 bins from 0 to 90 deg
    59 
    60     MBinning binsdist("BinningDist");
    61     binsdist.SetEdges(100, 0, 2);    // 100 bins from 0 to 2 deg
    62 
    63     plist.AddToList(&binswidth);
    64     plist.AddToList(&binslength);
    65     plist.AddToList(&binsalpha);
    66     plist.AddToList(&binsdist);
     65    // MBinning binssize("BinningSize");
     66    // binssize.SetEdgesLog(50, 1, 1e7);
     67    // plist.AddToList(&binssize);
    6768
    6869    //
    … …  
    7374    MSrcPosCam source("Source")
    7475    source.SetXY(0, 0);
    75 
    7676    plist.AddToList(&source);
    7777
    … …  
    8787    //read.AddFile("../data/CT1_97_off1.dat");
    8888
    89     MImgCleanStd   clean;
     89    MMcPedestalCopy pcopy;
     90    MMcPedestalNSBAdd pnsb;
     91
     92    MCerPhotCalc ncalc;
     93    //
     94    //  Alternative photon calculation:
     95    //  Example: use only 2nd to 6th FADC slices for photon calculation:
     96    //
     97    //    MCerPhotCalc2 ncalc;
     98    //    const Float_t x[15]={0, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0};
     99    //    TArrayF w(15,x);
     100    //    ncalc.SetWeights(w);
     101    //
     102
     103    MCameraSmooth   smooth;
     104    MImgCleanStd    clean;
     105    MBlindPixelCalc blind;
     106    //
     107    // Instead of unmapping the pixels you can also
     108    //
     109    // blind.SetUseInterpolation();
     110    // blind.SetUseCetralPixel();
     111    //
    90112    MHillasCalc    hcalc;
    91113    MHillasSrcCalc csrc1("Source",  "HillasSource");
    92114
    93     MFillH hfill("MHHillas",   "MHillas");
    94     MFillH sfill("MHStarMap",  "MHillas");
    95     MFillH hfill2s("HistSource  [MHHillasSrc]", "HillasSource");
     115    //
     116    // Uncomment this two line if you want to use MHillasExt instead
     117    // of MHillas
     118    //
     119    //MHillasExt hext;
     120    //plist.AddToList(&hext);
     121
     122    MFillH hfill1("MHHillas", "MHillas");
     123    MFillH hfill2("MHHillasExt");
     124    MFillH hfill3("MHStarMap", "MHillas");
     125    MFillH hfill4("HistExtSource [MHHillasExt]", "HillasSource");
     126    MFillH hfill5("HistSource [MHHillasSrc]", "HillasSource");
     127
     128    MWriteRootFile write("hillas.root");
     129    write.AddContainer("HillasSource",  "Hillas");
     130    write.AddContainer("MHStarMap");
     131    write.AddContainer("MMcEvt","Hillas");
     132    write.AddContainer("Source","RunHeaders");
    96133
    97134    tlist.AddToList(&read);
     135    tlist.AddToList(&pcopy);
     136    tlist.AddToList(&pnsb);
     137    tlist.AddToList(&ncalc);
    98138    tlist.AddToList(&clean);
     139    tlist.AddToList(&blind);
     140
    99141    tlist.AddToList(&hcalc);
    100142    tlist.AddToList(&csrc1);
    101     tlist.AddToList(&hfill);
    102     tlist.AddToList(&sfill);
    103     tlist.AddToList(&hfill2s);
     143    tlist.AddToList(&csrc2);
     144
     145    tlist.AddToList(&hfill1);
     146    tlist.AddToList(&hfill2);
     147    tlist.AddToList(&hfill3);
     148    tlist.AddToList(&hfill4);
     149    tlist.AddToList(&hfill5);
     150    tlist.AddToList(&write);
    104151
    105152    //
    … …  
    115162        return;
    116163
     164    tlist.PrintStatistics();
     165
    117166    //
    118167    // After the analysis is finished we can display the histograms
    119168    //
    120169    plist.FindObject("MHHillas")->DrawClone();
     170    plist.FindObject("MHHillasExt")->DrawClone();
     171    plist.FindObject("MHStarMap")->DrawClone();
    121172    plist.FindObject("HistSource")->DrawClone();
    122     plist.FindObject("MHStarMap")->DrawClone();
     173    plist.FindObject("HistExtSource")->DrawClone();
    123174}
     175

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

    r1463 r1490  
    2424
    2525
    26 void MagicHillas(const char *filename="~/data/Gamma*.root")
     26void MagicHillas(const char *filename="~/data/Gamma_20_N*.root")
    2727{
    2828    //
    … …  
    3535    // The tasklist is identified in the eventloop by its name
    3636    //
    37     MParList  plist;
     37    MParList plist;
     38
    3839
    3940    MTaskList tlist;
    4041    plist.AddToList(&tlist);
    41 
    42     //
    43     // Uncomment this two line if you want to use MHillasExt instead
    44     // of MHillas
    45     //
    46     /*
    47      MHillasExt hext;
    48      plist.AddToList(&hext);
    49      */
    5042
    5143    //
    … …  
    5850    //
    5951    // Use this if you want to change the binning of one of
    60     // the histograms
     52    // the histograms. You can use:
     53    // BinningConc, BinningConc1, BinningAsym, BinningM3Long,
     54    // BinningM3Trans, BinningWidth, BinningLength, BinningDist,
     55    // BinningHeadTail, BinningAlpha, BinningSize, BinningDelta,
     56    // BinningPixels and BinningCamera
    6157    //
    62     /*
    63      MBinning binswidth("BinningWidth");
    64      binswidth.SetEdges(100, 0, 0.1);     // 100 bins from 0 to 1 deg
    65      plist.AddToList(&binswidth);
     58    // For more information see MBinning and the corresponding
     59    // histograms
     60    //
     61    // MBinning binsalpha("BinningAlpha");
     62    // binsalpha.SetEdges(90, 0, 90);       // 90 bins from 0 to 90 deg
     63    // plist.AddToList(&binsalpha);
    6664
    67      MBinning binslength("BinningLength");
    68      binslength.SetEdges(100, 0, 1);      // 100 bins from 0 to 1 deg
    69      plist.AddToList(&binslength);
    70 
    71      MBinning binsdist("BinningDist");
    72      binsdist.SetEdges(100, 0, 1.5);      // 100 bins from 0 to 1 deg
    73      plist.AddToList(&binsdist);
    74 
    75      MBinning binsht("BinningHeadTail");
    76      binsht.SetEdges(100, -1.5, 1.5);     // 100 bins from 0 to 1 deg
    77      plist.AddToList(&binsht);
    78 
    79      MBinning binsalpha("BinningAlpha");
    80      binsalpha.SetEdges(90, 0, 90);       // 90 bins from 0 to 90 deg
    81      plist.AddToList(&binsalpha);
    82 
    83      MBinning binscam("BinningCamera");
    84      binscam.SetEdges(50, -1.5, 1.5);     //  20 bins from -2 to 2 deg
    85      plist.AddToList(&binscam);
    86     */
     65    // MBinning binssize("BinningSize");
     66    // binssize.SetEdgesLog(50, 1, 1e7);
     67    // plist.AddToList(&binssize);
    8768
    8869    //
    … …  
    11293    read.DisableAutoScheme();
    11394
    114     MMcPedestalCopy   pcopy;
     95    MMcPedestalCopy pcopy;
    11596    MMcPedestalNSBAdd pnsb;
    11697
    117     MCerPhotCalc      ncalc;
     98    MCerPhotCalc ncalc;
    11899    //
    119100    //  Alternative photon calculation:
    … …  
    126107    //
    127108
    128     MImgCleanStd      clean;
    129     MBlindPixelCalc   blind;
    130     MHillasCalc       hcalc;
    131     MHillasSrcCalc    csrc1("Source",  "HillasSource");
    132     MHillasSrcCalc    csrc2("AntiSrc", "HillasAntiSrc");
     109    MCameraSmooth   smooth;
     110    MImgCleanStd    clean;
     111    MBlindPixelCalc blind;
     112    //
     113    // Instead of unmapping the pixels you can also
     114    //
     115    // blind.SetUseInterpolation();
     116    // blind.SetUseCetralPixel();
     117    //
     118    MHillasCalc    hcalc;
     119    MHillasSrcCalc csrc1("Source",  "HillasSource");
     120    MHillasSrcCalc csrc2("AntiSrc", "HillasAntiSrc");
    133121
    134     MFillH hfill("MHHillas",   "MHillas");
    135     MFillH sfill("MHStarMap",  "MHillas");
    136     MFillH hfill2s("HistSource  [MHHillasSrc]", "HillasSource");
    137     MFillH hfill2a("HistAntiSrc [MHHillasSrc]", "HillasAntiSrc");
     122    //
     123    // Uncomment this two line if you want to use MHillasExt instead
     124    // of MHillas
     125    //
     126    //MHillasExt hext;
     127    //plist.AddToList(&hext);
     128
     129    MFillH hfill1("MHHillas", "MHillas");
     130    MFillH hfill2("MHHillasExt");
     131    MFillH hfill3("MHStarMap", "MHillas");
     132    MFillH hfill4("HistExtSource [MHHillasExt]", "HillasSource");
     133    MFillH hfill5("HistExtAntiSource [MHHillasExt]", "HillasAntiSource");
     134    MFillH hfill6("HistSource [MHHillasSrc]", "HillasSource");
     135    MFillH hfill7("HistAntiSrc [MHHillasSrc]", "HillasAntiSrc");
    138136
    139137    MWriteRootFile write("hillas.root");
    140     write.AddContainer("MHillas",       "Hillas");
    141138    write.AddContainer("HillasSource",  "Hillas");
    142139    write.AddContainer("HillasAntiSrc", "Hillas");
    143140    write.AddContainer("MHStarMap");
    144141    write.AddContainer("MMcEvt","Hillas");
    145 
    146     /*
    147      MWriteAsciiFile write("hillas.txt");
    148      write.AddColumn("MHillas.fLength*MGeomCam.fConvMm2Deg");
    149      write.AddColumn("MHillas", "fConc");
    150      write.AddColumns("MHillas");
    151      */
     142    write.AddContainer("Source","RunHeaders");
     143    write.AddContainer("AntiSource","RunHeaders");
    152144
    153145    tlist.AddToList(&read);
    … …  
    157149    tlist.AddToList(&clean);
    158150    tlist.AddToList(&blind);
     151
    159152    tlist.AddToList(&hcalc);
    160153    tlist.AddToList(&csrc1);
     154    tlist.AddToList(&csrc2);
    161155
    162     tlist.AddToList(&csrc2);
    163     tlist.AddToList(&hfill);
    164     tlist.AddToList(&sfill);
    165     tlist.AddToList(&hfill2s);
    166     tlist.AddToList(&hfill2a);
     156    tlist.AddToList(&hfill1);
     157    tlist.AddToList(&hfill2);
     158    tlist.AddToList(&hfill3);
     159    tlist.AddToList(&hfill4);
     160    tlist.AddToList(&hfill5);
     161    tlist.AddToList(&hfill6);
     162    tlist.AddToList(&hfill7);
    167163    tlist.AddToList(&write);
    168164
    … …  
    184180    // After the analysis is finished we can display the histograms
    185181    //
     182    plist.FindObject("MHHillas")->DrawClone();
     183    plist.FindObject("MHHillasExt")->DrawClone();
     184    plist.FindObject("MHStarMap")->DrawClone();
    186185    plist.FindObject("HistSource")->DrawClone();
    187     plist.FindObject("MHHillas")->DrawClone();
    188     plist.FindObject("HistAntiSrc")->DrawClone();
    189     plist.FindObject("MHStarMap")->DrawClone();
     186    plist.FindObject("HistAntiSource")->DrawClone();
     187    plist.FindObject("HistExtSource")->DrawClone();
     188    plist.FindObject("HistExtAntiSource")->DrawClone();
    190189}
    191190

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

    r1489 r1490  
    7979    read.DisableAutoScheme();
    8080
    81     MFillH hfill("MHHillasExt", "MHillas");
    82     MFillH sfill("MHStarMap",   "MHillas");
    83     MFillH hfill2s("HistSource  [MHHillasSrc]", "HillasSource");
    84     MFillH hfill2a("HistAntiSrc [MHHillasSrc]", "HillasAntiSrc");
     81    MFillH hfill1("MHHillas", "MHillas");
     82    MFillH hfill2("MHHillasExt");
     83    MFillH hfill3("MHStarMap", "MHillas");
     84    MFillH hfill4("HistExtSource [MHHillasExt]", "HillasSource");
     85    MFillH hfill5("HistExtAntiSource [MHHillasExt]", "HillasAntiSource");
     86    MFillH hfill6("HistSource [MHHillasSrc]", "HillasSource");
     87    MFillH hfill7("HistAntiSrc [MHHillasSrc]", "HillasAntiSrc");
    8588
    8689    tlist.AddToList(&read);
    87     tlist.AddToList(&hfill);
    88     tlist.AddToList(&sfill);
    89     tlist.AddToList(&hfill2s);
    90     tlist.AddToList(&hfill2a);
     90    tlist.AddToList(&hfill1);
     91    tlist.AddToList(&hfill2);
     92    tlist.AddToList(&hfill3);
     93    tlist.AddToList(&hfill4);
     94    tlist.AddToList(&hfill5);
     95    tlist.AddToList(&hfill6);
     96    tlist.AddToList(&hfill7);
    9197
    9298    //
    … …  
    107113    // After the analysis is finished we can display the histograms
    108114    //
     115    plist.FindObject("MHHillas")->DrawClone();
     116    plist.FindObject("MHHillasExt")->DrawClone();
     117    plist.FindObject("MHStarMap")->DrawClone();
    109118    plist.FindObject("HistSource")->DrawClone();
    110     plist.FindObject("MHHillasExt")->DrawClone();
    111     plist.FindObject("HistAntiSrc")->DrawClone();
    112     plist.FindObject("MHStarMap")->DrawClone();
     119    plist.FindObject("HistAntiSource")->DrawClone();
     120    plist.FindObject("HistExtSource")->DrawClone();
     121    plist.FindObject("HistExtAntiSource")->DrawClone();
    113122}
    114123
Note: See TracChangeset for help on using the changeset viewer.