Changeset 698 for trunk/MagicSoft


Ignore:
Timestamp:
03/22/01 12:24:47 (24 years ago)
Author:
tbretz
Message:
*** empty log message ***
Location:
trunk/MagicSoft
Files:
17 added
4 deleted
49 edited

Legend:

Unmodified
Added
Removed

  • trunk/MagicSoft/Mars/Changelog

    r695 r698  
    11                                                                  -*-*- END -*-*-
     2
     3 2000/03/21: Thomas Bretz
     4                                                                 
     5   * macros/CT1Hillas.C, manalysis/MImgCleanStd.cc, manalysis/MImgCleanStd.h,
     6     manalysis/MMcPedestalCopy.cc, manalysis/MMcPedestalCopy.h,
     7     manalysis/MCerPhotCalc.h, manalysis/MCerPhotCalc.cc, mdatacheck/MHHillas.h,
     8     mdatacheck/MHHillas.cc, mdatacheck/MFillHFadc.cc, mdatacheck/MFillHFadc.h,
     9     mdatacheck/MHFadcCam.cc, mdatacheck/MHFadcCam.h, mdatacheck/MHFadcPix.cc,
     10     mdatacheck/MHFadcPix.h, mdatacheck/MFillHHillas.cc,
     11     mdatacheck/MFillHHillas.h
     12     - added
     13
     14   * mdatacheck/MHistosAdc.[h,cc]:
     15     - replaced by MHFadc[Cam,Pix]
     16     
     17   * mdatacheck/MFillAdcSpect.[h,cc]:
     18     - replaced by MFillHFadc
     19     
     20   * macros/readCT1.C:
     21     - uses the new classes now
     22     - calculates and displays hillas
     23     
     24   * manalysis/AnalysisLinkDef.h:
     25     - added MCerPhotCalc, MImgCleanStd, MMcPedestalCopy
     26     
     27   * manalysis/MCT1ReadAscii.cc:
     28     - PreProcess uses FindCreateObj now
     29     
     30   * manalysis/MCerPhotEvt.cc, manalysis/MPedestalCam.cc:
     31     - changed the new back to 'new with placement'
     32     
     33   * manalysis/MCerPhotEvt.[h, cc]:
     34     - removed cleaning from this class
     35     
     36   * manalysis/MHillas.cc:
     37     - some small changes in the naming inside Calc()
     38     
     39   * manalysis/MHillas.h:
     40     - added Get-functions
     41     
     42   * manalysis/MPedestalCam.h:
     43     - added InitSize-function
     44     
     45   * manalysis/MPedestalPix.h:
     46     - made Get-functions const
     47     
     48   * manalysis/Makefile:
     49     - added MMcPedestalCopy, MImgCleanStd, MCerPhotCalc
     50     
     51   * mbase/MLogManip.h:
     52     - "l." before line number
     53     
     54   * mbase/MParList.cc:
     55     - added the autodelete feature for automatically created classes
     56     
     57   * mdatacheck/DataCheckLinkDef.h:
     58     - added MHFadc[Cam,Pix], MHHillas
     59     - removed MHistosAdc
     60     - added MFillHFadc, MFillHHillas
     61     
     62   * mdatacheck/MGDisplayAdc.[h,cc], MShowSpect.[h,cc],
     63     mdatacheck/MViewAdcSpectra.[h,cc], :
     64     - changed from MHistaosAdc to MHFadcCam
     65     
     66   * mdatacheck/Makefile:
     67     - added MFillHFadc.cc, MFillHHillas.cc
     68     - removed MFillAdcSpect.cc
     69
     70   * manalysis/MCT1ReadAscii.h, mbase/MEvtLoop.h, mbase/MLog.h,
     71     mbase/MParList.h, mbase/MReadTree.h, mbase/MTask.h,
     72     mbase/MTaskList.h, mdatacheck/MDumpEvtHeader.h,
     73     meventdisp/MGFadcDisp.h, mgui/MCamDisplay.h,
     74     mgui/MCamNeighbor.h, mgui/MGDataCheckMain.h,
     75     mgui/MGEventDispMain.h, mgui/MGMarsMain.h,
     76     mgui/MGMonteCarloMain.h, mgui/MGPrototyp.h, mgui/MHexagon.h,
     77     mmontecarlo/MCollAreaTrigger.h, mraw/MRawEvtPixelIter.h,
     78     mraw/MRawFileRead.h, mraw/MRawFileWrite.h:
     79     - set version number to '0' (no i/o for this class)
     80
     81   * mgui/MGeomCamCT1.cc:
     82     - new claculation algorith for the CT1 camera. It is more
     83       exact in the camera plain
     84       
     85   * mgui/MHexagon.cc:
     86     - removed Draw-function (it is inherited from TObject)
     87     
     88   * mmc/MMcFadcHeader.hxx:
     89     - Added get-functions
     90     
     91
    292
    393 2000/03/20: Thomas Bretz

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

    r695 r698  
    33    MParList plist;
    44
    5     MGeomCamMagic camct1;
    6     MHillas       hillas;
    7     MTaskList     tlist;
     5    MGeomCamCT1 camct1;
     6    MHillas     hillas;
     7    MTaskList   tlist;
    88
    99    plist->AddToList(&camct1);
    … …  
    1212
    1313    MCT1ReadAscii readct1("CT1_99_on1.dat") ;
    14     MHillasCalc hcalc;
     14    MImgCleanStd  clean;
     15    MHillasCalc   hcalc;
    1516
    1617    tlist.AddToList(&readct1);
     18    tlist.AddToList(&clean);
    1719    tlist.AddToList(&hcalc);
    1820
    … …  
    4042            break;
    4143
    42         phevt.CleanLevel1();
    43         phevt.CleanLevel2();
    44         display.DrawPhotNum(&phevt);
    45         gClient->HandleInput();
    46         if(getchar()=='q')
    47             break;
    48 
    49         phevt.CleanLevel3();
    50 
     44        clean.Process();
    5145        hcalc.Process();
    5246

  • trunk/MagicSoft/Mars/manalysis/AnalysisLinkDef.h

    r695 r698  
    77#pragma link C++ class MCerPhotPix;
    88#pragma link C++ class MCerPhotEvt;
     9#pragma link C++ class MCerPhotCalc;
     10
     11#pragma link C++ class MImgCleanStd;
    912
    1013#pragma link C++ class MCT1ReadAscii;
    … …  
    1316#pragma link C++ class MPedestalPix;
    1417#pragma link C++ class MPedestalCam;
     18#pragma link C++ class MMcPedestalCopy;
    1519
    1620#pragma link C++ class MHillas;

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

    r654 r698  
    3434    // skip the next 4 values
    3535    //
    36     fin >> dummyF;
    37     fin >> dummyF;
    38     fin >> dummyF;
    39     fin >> dummyF;
     36    fin >> dummyF; 
     37    fin >> dummyF; 
     38    fin >> dummyF; 
     39    fin >> dummyF; 
    4040
    4141    //

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

    r656 r698  
    1010
    1111#include "MLog.h"
     12#include "MLogManip.h"
     13
    1214#include "MParList.h"
    1315#include "MCerPhotEvt.h"
    … …  
    4143    if (!(*fIn))
    4244    {
    43         *fLog << "Error: MCT1ReadAscii::PreProcess: Cannot open file." << endl;
     45        *fLog << dbginf << "Cannot open file." << endl;
    4446        return kFALSE;
    4547    }
    … …  
    4850    //  look for the MCerPhotEvt class in the plist
    4951    //
    50     fNphot = (MCerPhotEvt*)pList->FindObject("MCerPhotEvt");
     52    fNphot = (MCerPhotEvt*)pList->FindCreateObj("MCerPhotEvt");
    5153    if (!fNphot)
    52     {
    53         *fLog << "MRawCT1Ascii::PreProcess - WARNING: MCerPhotEvt not found... creating." << endl;
    54         fNphot = new MCerPhotEvt;
    55         pList->AddToList(fNphot);
    56     }
     54        return kFALSE;
    5755
    5856    //
    5957    //  look for the pedestal class in the plist
    6058    //
    61     fPedest = (MCT1Pedestals*)pList->FindObject("MCT1Pedestals");
     59    fPedest = (MCT1Pedestals*)pList->FindCreateObj("MCT1Pedestals");
    6260    if (!fPedest)
    63     {
    64         *fLog << "MRawCT1Ascii::PreProcess - WARNING: MCT1Pedestals not found... creating." << endl;
    65         fPedest = new MCT1Pedestals;
    66         pList->AddToList(fPedest);
    67     }
     61        return kFALSE;
    6862
    6963    return kTRUE;
    … …  
    7367{
    7468    //
     69    // FIXME. This function should switch between reading pedestals and
     70    // reading event data by the 'switch entry'.
     71    // After reading it should set the InputStreamID correctly.
     72    // ( should use MPedestalCam )
     73    //
     74 
     75    //
    7576    // read in a dummy number (event number)
    7677    //
    … …  
    7879    *fIn >> dummyI;
    7980
     81    //
     82    // check if we are done
     83    //
    8084    if (fIn->eof())
    81     {
    82         *fLog << "MRawCT1Ascii::Process - Error: EOF reached." << endl;
    8385        return kFALSE;
    84     }
    8586
    8687    //
    … …  
    9495    // five unsused numbers
    9596    //
    96     *fIn >> dummyI;
    97     *fIn >> dummyI;
    98     *fIn >> dummyI;
    99     *fIn >> dummyI;
     97    *fIn >> dummyI;   // ener
     98    *fIn >> dummyI;   // zenang
     99    *fIn >> dummyI;   // sec1
     100    *fIn >> dummyI;   // sec2
    100101
    101102    //
    … …  
    110111    for (Int_t i = 0; i<127; i++ )
    111112    {
    112         Float_t dummyF;
     113        Float_t nphot;
    113114
    114         *fIn >> dummyF;
     115        *fIn >> nphot;
    115116
    116         if (dummyF > 0.0)
    117             fNphot->AddPixel(i, dummyF, (*fPedest)[i]);
     117        if (nphot > 0.0)
     118            fNphot->AddPixel(i, nphot, (*fPedest)[i]);
    118119    }
    119120

  • trunk/MagicSoft/Mars/manalysis/MCT1ReadAscii.h

    r653 r698  
    2626    Bool_t PostProcess();
    2727
    28     ClassDef(MCT1ReadAscii, 1)  // Reads the CT1 data file
     28    ClassDef(MCT1ReadAscii, 0)  // Reads the CT1 data file
    2929};
    3030

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

    r695 r698  
    4747}
    4848
    49 Int_t MCerPhotEvt::GetNumPixels()
    50 {
    51     return fNumPixels;
    52 }
    53 
    5449void MCerPhotEvt::AddPixel(Int_t id, Float_t nph, Float_t err)
    5550{
    … …  
    5853    // of valid pixels in the list by one
    5954    //
    60     (*fPixels)[fNumPixels++] = new MCerPhotPix( id, nph, err);
     55
     56    // TClonesArray -> 'operator new with placement' should be used
     57    new ((*fPixels)[fNumPixels++]) MCerPhotPix( id, nph, err);
    6158}
    6259
    … …  
    8380}
    8481
     82/*
    8583void MCerPhotEvt::CleanLevel1()
    8684{
    … …  
    228226    }
    229227}
    230 
     228*/
    231229
    232230Bool_t MCerPhotEvt::IsPixelExisting(Int_t id)
    … …  
    293291
    294292    Float_t testval;
    295     for (Int_t i=1 ; i<fNumPixels; i++ )
     293    for (UInt_t i=1 ; i<fNumPixels; i++ )
    296294    {
    297295        testval = (*this)[i].GetNumPhotons();
    … …  
    315313
    316314    Float_t testval;
    317     for (Int_t i=1; i<fNumPixels; i++)
     315    for (UInt_t i=1; i<fNumPixels; i++)
    318316    {
    319317        testval = (*this)[i].GetNumPhotons();

  • trunk/MagicSoft/Mars/manalysis/MCerPhotEvt.h

    r695 r698  
    66#endif
    77#ifndef ROOT_TClonesArray
    8 #include "TClonesArray.h"
     8#include <TClonesArray.h>
    99#endif
    1010#ifndef MPARCONTAINER_H
    … …  
    2121private:
    2222
    23     Int_t         fNumPixels;
     23    UInt_t         fNumPixels;
    2424    TClonesArray *fPixels;
    2525
    … …  
    2727
    2828public:
    29   MCerPhotEvt(const char *name=NULL, const char *title=NULL) ;
     29    MCerPhotEvt(const char *name=NULL, const char *title=NULL) ;
     30    ~MCerPhotEvt() { delete fPixels; }
    3031
    3132  void Draw(Option_t* option = "" ) ;
    3233
    33   Int_t    GetNumPixels() ;
     34  UInt_t GetNumPixels() const { return fNumPixels; }
    3435
    3536  void AddPixel(Int_t id, Float_t nph, Float_t err );
    … …  
    3940  void Print(Option_t *opt=NULL) ;
    4041
    41   void CleanLevel1() ;
    42   void CleanLevel2() ;
    43   void CleanLevel3() ;
     42  //  void CleanLevel1() ;
     43  //  void CleanLevel2() ;
     44  //  void CleanLevel3() ;
    4445 
    45   Bool_t  IsPixelExisting( Int_t id ) ;
    46   Bool_t  IsPixelUsed    ( Int_t id ) ;
    47   Bool_t  IsPixelCore    ( Int_t id ) ;
     46  Bool_t  IsPixelExisting( Int_t id );
     47  Bool_t  IsPixelUsed    ( Int_t id );
     48  Bool_t  IsPixelCore    ( Int_t id );
    4849 
    4950  Float_t GetNumPhotonsMin();

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

    r695 r698  
    7979    // calculate mean valu of pixels
    8080    //
    81     float xav =0;
    82     float yav =0;
    83     float spix=0;
     81    float xmean =0;
     82    float ymean =0;
     83
     84    fSize = 0;
    8485
    8586    for (UInt_t i=0; i<nevt; i++)
    … …  
    9495        const float nphot = pix.GetNumPhotons();
    9596
    96         spix += nphot;
    97         xav  += nphot * gpix.GetX();
    98         yav  += nphot * gpix.GetY();
     97        fSize += nphot;
     98        xmean += nphot * gpix.GetX(); // [mm]
     99        ymean += nphot * gpix.GetY(); // [mm]
    99100    }
    100101
    101     xav /= spix;
    102     yav /= spix;
     102    xmean /= fSize; // [mm]
     103    ymean /= fSize; // [mm]
    103104
    104105    //
    105106    // calculate sdev
    106107    //
    107     float dis11=0;
    108     float dis12=0;
    109     float dis22=0;
     108    float sigmaxx=0;
     109    float sigmaxy=0;
     110    float sigmayy=0;
    110111
    111112    for (UInt_t i=0; i<nevt; i++)
    … …  
    118119        const MGeomPix &gpix = geom[pix.GetPixId()];
    119120
    120         const float dx = gpix.GetX() - xav;
    121         const float dy = gpix.GetY() - yav;
     121        const float dx = gpix.GetX() - xmean;
     122        const float dy = gpix.GetY() - ymean;
    122123
    123124        const float nphot = pix.GetNumPhotons();
    124125
    125         dis11 += nphot * dx*dx;
    126         dis12 += nphot * dx*dy;
    127         dis22 += nphot * dy*dy;
     126        sigmaxx += nphot * dx*dx; // [mm^2]
     127        sigmaxy += nphot * dx*dy; // [mm^2]
     128        sigmayy += nphot * dy*dy; // [mm^2]
    128129    }
    129130
    … …  
    131132    // check for orientation
    132133    //
    133     const float theta = atan(dis12/(dis11-dis22)*2)/2;
     134    const float theta = atan(sigmaxy/(sigmaxx-sigmayy)*2)/2;
    134135
    135     float c = cos(theta);
    136     float s = sin(theta);
     136    float c = cos(theta); // [1]
     137    float s = sin(theta); // [1]
    137138
    138     const float direction = c*xav+s*yav;
     139    const float direction = c*xmean + s*ymean;
    139140
    140141    if (direction<0)
    … …  
    144145    }
    145146
    146     float rot1 =  2.0*c*s*dis12 + c*c*dis11 + s*s*dis22;
    147     float rot2 = -2.0*c*s*dis12 + s*s*dis11 + c*c*dis22;
     147    float axis1 =  2.0*c*s*sigmaxy + c*c*sigmaxx + s*s*sigmayy; // [mm^2]
     148    float axis2 = -2.0*c*s*sigmaxy + s*s*sigmaxx + c*c*sigmayy; // [mm^2]
    148149
    149     rot1 /= spix;
    150     rot2 /= spix;
     150    axis1 /= fSize; // [mm^2]
     151    axis2 /= fSize; // [mm^2]
    151152
    152153    //
    153154    // check for numerical negatives
    154155    //
    155     if (rot1 < 0) rot1=0;
    156     if (rot2 < 0) rot2=0;
     156    if (axis1 < 0) axis1=0;
     157    if (axis2 < 0) axis2=0;
    157158
    158159    //
    159160    // check the rotation of the axis
    160161    //
    161     const int rotation = rot1<rot2;
     162    const int rotation = axis1<axis2;
    162163
    163     fLength = rotation ? sqrt(rot2) : sqrt(rot1);
    164     fWidth  = rotation ? sqrt(rot1) : sqrt(rot2);
     164    fLength = rotation ? sqrt(axis2) : sqrt(axis1);      // [mm]
     165    fWidth  = rotation ? sqrt(axis1) : sqrt(axis2);      // [mm]
    165166
    166     fAlpha = 180/kPI*atan((-xav*s+yav*c)/direction);
     167    fAlpha = 180/kPI*atan((-xmean*s+ymean*c)/direction); // [deg]
    167168
    168     fDist  = sqrt(xav*xav + yav*yav);
     169    fDist  = sqrt(xmean*xmean + ymean*ymean);            // [mm]
    169170
    170     fTheta = atan(yav/xav);
    171     if (xav<0) fTheta += kPI;
     171    fTheta = atan(ymean/xmean);                          // [rad]
     172    if (xmean<0) fTheta += kPI;                          // [deg]
    172173}

  • trunk/MagicSoft/Mars/manalysis/MHillas.h

    r695 r698  
    3434    void Clear(Option_t *opt=NULL);
    3535
     36    Float_t GetAlpha() const  { return fAlpha; }
     37    Float_t GetWidth() const  { return fWidth; }
     38    Float_t GetLength() const { return fLength; }
     39    Float_t GetDist() const   { return fDist; }
     40
    3641    ClassDef(MHillas, 1) // Storage Container for Hillas Parameter
    3742};

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

    r695 r698  
    1818
    1919    fArray = new TClonesArray("MPedestalPix", 577);
     20
     21    //
     22    // TClonesArray: The 'new operator with placement' must be used
     23    //
     24    for (int i=0; i<577; i++)
     25        new ((*fArray)[i]) MPedestalPix;
    2026}
    2127

  • trunk/MagicSoft/Mars/manalysis/MPedestalCam.h

    r695 r698  
    2323    ~MPedestalCam();
    2424
     25    void InitSize(const UInt_t i) { fArray->ExpandCreateFast(i); }
     26
    2527    MPedestalPix &operator[](Int_t i) { return *(MPedestalPix*)fArray->At(i); }
    2628

  • trunk/MagicSoft/Mars/manalysis/MPedestalPix.h

    r695 r698  
    1717    MPedestalPix();
    1818
    19     Float_t GetMean()     { return fMean;     }
    20     Float_t GetSigma()    { return fSigma;    }
    21     Float_t GetMeanRms() { return fMeanRms;  }
    22     Float_t GetSigmaRms() { return fSigmaRms; }
     19    Float_t GetMean() const     { return fMean;     }
     20    Float_t GetSigma() const    { return fSigma;    }
     21    Float_t GetMeanRms() const { return fMeanRms;  }
     22    Float_t GetSigmaRms() const { return fSigmaRms; }
    2323
    2424    void SetMean(Float_t f)     { fMean = f;     }
    … …  
    2727    void SetSigmaRms(Float_t f) { fSigmaRms = f; }
    2828
     29    void SetPedestal(Float_t m, Float_t s) { fMean = m; fSigma = s; }
     30    void SetPedestalRms(Float_t m, Float_t s) { fMeanRms = m; fSigmaRms = s; }
     31
    2932    ClassDef(MPedestalPix, 1)   // Storage Container for Pedestal information of one pixel
    3033};

  • trunk/MagicSoft/Mars/manalysis/Makefile

    r695 r698  
    2222#  connect the include files defined in the config.mk file
    2323#
    24 INCLUDES = -I. -I../mbase -I../mgui
     24INCLUDES = -I. -I../mbase -I../mgui -I../mmc -I../mraw
    2525
    2626#------------------------------------------------------------------------------
    … …  
    3232           MPedestalCam.cc \
    3333           MPedestalPix.cc \
     34           MMcPedestalCopy.cc \
     35           MImgCleanStd.cc \
    3436           MHillas.cc \
    3537           MHillasCalc.cc \
     38           MCerPhotCalc.cc \
    3639           MCerPhotEvt.cc \
    3740           MCerPhotPix.cc

  • trunk/MagicSoft/Mars/mbase/MEvtLoop.h

    r604 r698  
    3939    void Eventloop(Int_t maxcnt=-1, const char *tlist="MTaskList");
    4040
    41     ClassDef(MEvtLoop, 1) // Class to execute the tasks in a tasklist
     41    ClassDef(MEvtLoop, 0) // Class to execute the tasks in a tasklist
    4242};
    4343

  • trunk/MagicSoft/Mars/mbase/MLog.h

    r609 r698  
    116116    }
    117117
    118     ClassDef(MLog, 1)
     118    ClassDef(MLog, 0)
    119119};
    120120

  • trunk/MagicSoft/Mars/mbase/MLogManip.h

    r654 r698  
    5353//
    5454#ifndef __CINT__
    55 #define dbginf        __FILE__ << " " << __LINE__ << ": "
     55#define dbginf        __FILE__ << " l." << __LINE__ << ": "
    5656#endif
    5757//

  • trunk/MagicSoft/Mars/mbase/MParList.cc

    r667 r698  
    2424ClassImp(MParList)
    2525
     26    //
     27    // FIXME: The Automatic created classes are NEVER deleted!!!!!
     28    //
     29
    2630MParList::MParList(const char *name, const char *title)
    2731{
    … …  
    3236    *fName  = name  ? name  : "MParList";
    3337    *fTitle = title ? title : "List of Parameter Containers";
     38
     39    //
     40    // This sets a flag that the list is the owner, which means
     41    // that the destructor of the list deletes all it's objects
     42    //
     43    fAutodelete.SetOwner();
    3444}
    3545
    … …  
    151161    AddToList(pcont);
    152162
     163    fAutodelete.Add(pcont);
     164
    153165    //
    154166    //  Find an object in the list.

  • trunk/MagicSoft/Mars/mbase/MParList.h

    r667 r698  
    2727private:
    2828    TOrdCollection fContainer;  // Collection of Parameter and Data Containers
     29    TOrdCollection fAutodelete; // All what this list contains is deleted in the destructor
    2930
    3031public:
    … …  
    4546    void Print(Option_t *t = NULL);
    4647
    47     ClassDef(MParList, 1)       //list of parameter containers (MParContainer)
     48    ClassDef(MParList, 0) // list of parameter containers (MParContainer)
    4849};
    4950

  • trunk/MagicSoft/Mars/mbase/MReadTree.h

    r585 r698  
    3636    UInt_t GetEntries() {return fNumEntries;}
    3737
    38     ClassDef(MReadTree, 1)      // Reads one tree
     38    ClassDef(MReadTree, 0)      // Reads one tree
    3939};
    4040

  • trunk/MagicSoft/Mars/mbase/MTask.h

    r609 r698  
    2727    virtual Bool_t PostProcess();
    2828
    29     ClassDef(MTask, 1)          //Abstract base class for a task
     29    ClassDef(MTask, 0)          //Abstract base class for a task
    3030};
    3131

  • trunk/MagicSoft/Mars/mbase/MTaskList.h

    r609 r698  
    4141    void Print(Option_t *t = NULL);
    4242
    43     ClassDef(MTaskList, 1)      //collection of tasks to be performed in the eventloop
     43    ClassDef(MTaskList, 0)      //collection of tasks to be performed in the eventloop
    4444};
    4545

  • trunk/MagicSoft/Mars/mdatacheck/DataCheckLinkDef.h

    r666 r698  
    55#pragma link off all functions;
    66
     7#pragma link C++ class MHFadcCam;
     8#pragma link C++ class MHFadcPix;
     9#pragma link C++ class MHHillas;
     10
    711#pragma link C++ class MShowSpect;
    8 #pragma link C++ class MHistosAdc;
    9 #pragma link C++ class MFillAdcSpect;
     12
     13#pragma link C++ class MFillHFadc;
     14#pragma link C++ class MFillHHillas;
    1015
    1116#pragma link C++ class MDumpEvtHeader;

  • trunk/MagicSoft/Mars/mdatacheck/MDumpEvtHeader.h

    r666 r698  
    2626  Bool_t Process() ;
    2727 
    28   ClassDef(MDumpEvtHeader, 1)   // Class to dump the pixel ids of a raw evt to the screen
     28  ClassDef(MDumpEvtHeader, 0)   // Class to dump the pixel ids of a raw evt to the screen
    2929
    3030};

  • trunk/MagicSoft/Mars/mdatacheck/MGDisplayAdc.cc

    r667 r698  
    2323} ;
    2424
    25 MGDisplayAdc::MGDisplayAdc ( MHistosAdc *histos,
     25MGDisplayAdc::MGDisplayAdc ( MHFadcCam *histos,
    2626                            const TGWindow *p, const TGWindow *main,
    2727                            UInt_t w, UInt_t h,
    … …  
    218218    //   Histogramms in there.
    219219    //
    220     //   In the class MHistosAdc are in fact two lists. One for the high and
     220    //   In the class MHFadcCam are in fact two lists. One for the high and
    221221    //   one for the low gain. Here we will use only the high gain list!!!
    222222    //   With some special options (settings in the gui) we will also be able
    … …  
    224224    //
    225225
    226     const Int_t nhi = fHists->GetHiEntries();
     226    const Int_t nhi = fHists->GetEntries();
    227227    for ( Int_t i=0 ; i < nhi; i++ ) {
    228228        fHistoList->AddEntry(fHists->GetHistHi(i)->GetName(), i+1) ;

  • trunk/MagicSoft/Mars/mdatacheck/MGDisplayAdc.h

    r665 r698  
    1212#endif
    1313
    14 #ifndef MHISTOSADC_H
    15 #include "MHistosAdc.h"
     14#ifndef MHFADCCAM_H
     15#include "MHFadcCam.h"
    1616#endif
    1717
    … …  
    2424    private:
    2525       
    26         MHistosAdc *fHists;             // Pointer to Container with the histograms
     26        MHFadcCam *fHists;              // Pointer to Container with the histograms
    2727       
    2828        // Create a main frame with a number of different buttons.
    … …  
    7171        TCanvas *fCanvas;
    7272       
    73         MGDisplayAdc(MHistosAdc *fHists ,
     73        MGDisplayAdc(MHFadcCam *fHists ,
    7474                     const TGWindow *p, const TGWindow *main,
    7575                     UInt_t w, UInt_t h,

  • trunk/MagicSoft/Mars/mdatacheck/MShowSpect.cc

    r667 r698  
    3434  //   Do the preprocessing for MShowSpect
    3535  //
    36   //   Connects Histogramms in the MHistosAdc container as the input
     36  //   Connects Histogramms in the MHFadcCam container as the input
    3737  //
    3838
    39   fHists = (MHistosAdc*) pList->FindObject( fHistName );
     39  fHists = (MHFadcCam*) pList->FindObject( fHistName );
    4040 
    4141  if (!fHists)
    4242    {
    43         *fLog << dbginf << " Error: MHistosAdc '" << fHistName << "' not found!" << endl;
     43        *fLog << dbginf << " Error: MHFadcCam '" << fHistName << "' not found!" << endl;
    4444        return kFALSE;
    4545    }

  • trunk/MagicSoft/Mars/mdatacheck/MShowSpect.h

    r667 r698  
    1111
    1212class MGDisplayAdc;
    13 class MHistosAdc;
     13class MHFadcCam;
    1414class MParList ;
    1515
    … …  
    1919    char       fHistName[256] ;
    2020
    21     MHistosAdc *fHists;         // Pointer to Container with the histograms
     21    MHFadcCam *fHists;          // Pointer to Container with the histograms
    2222 
    2323 public:
    … …  
    2929     Bool_t PostProcess();
    3030
    31      ClassDef(MShowSpect, 1)    // Fill the raw ADC in the histograms
     31     ClassDef(MShowSpect, 0)    // Fill the raw ADC in the histograms
    3232};
    3333

  • trunk/MagicSoft/Mars/mdatacheck/MViewAdcSpectra.cc

    r667 r698  
    1515#include "MReadTree.h"
    1616#include "MDumpEvtHeader.h"
    17 #include "MFillAdcSpect.h"
     17#include "MFillHFadc.h"
    1818#include "MShowSpect.h"
    19 #include "MHistosAdc.h"
     19#include "MHFadcCam.h"
    2020
    2121MViewAdcSpectra::MViewAdcSpectra()
    … …  
    6565  plist.AddToList(&evttime);
    6666
    67   fHistosAdc= new MHistosAdc;
     67  fHistosAdc= new MHFadcCam;
    6868  plist.AddToList( fHistosAdc ) ;
    6969
    … …  
    8080  //  tasks->AddToList( dumpheader ) ;
    8181 
    82   MFillAdcSpect fillspect;
     82  MFillHFadc fillspect;
    8383  tasks.AddToList( &fillspect ) ;
    8484
    85   MShowSpect showspect( "MHistosAdc" ) ;
     85  MShowSpect showspect( "MHFadcCam" ) ;
    8686  tasks.AddToList( &showspect ) ;
    8787

  • trunk/MagicSoft/Mars/mdatacheck/MViewAdcSpectra.h

    r667 r698  
    66#endif
    77
    8 class MHistosAdc;
     8class MHFadcCam;
    99
    1010class MViewAdcSpectra {
    … …  
    1212 private:
    1313
    14      MHistosAdc *fHistosAdc;
     14     MHFadcCam *fHistosAdc;
    1515
    1616 public:

  • trunk/MagicSoft/Mars/mdatacheck/Makefile

    r665 r698  
    2222#  connect the include files defined in the config.mk file
    2323#
    24 INCLUDES = -I. -I../mbase -I../mraw
     24INCLUDES = -I. -I../mbase -I../mraw -I../manalysis
    2525
    2626#------------------------------------------------------------------------------
    … …  
    2929
    3030SRCFILES = MDumpEvtHeader.cc \
    31            MFillAdcSpect.cc \
     31           MFillHFadc.cc \
     32           MFillHHillas.cc \
    3233           MGDisplayAdc.cc \
    33            MHistosAdc.cc \
     34           MHFadcPix.cc \
     35           MHFadcCam.cc \
     36           MHHillas.cc \
    3437           MShowSpect.cc \
    3538           MViewAdcSpectra.cc

  • trunk/MagicSoft/Mars/meventdisp/MGFadcDisp.cc

    r695 r698  
    3636
    3737  //   first connect the file with this Object
    38  
    39   pList = new MParList();
     38
     39    pList = new MParList();
    4040
    4141  fRunHeader = new MRawRunHeader();

  • trunk/MagicSoft/Mars/meventdisp/MGFadcDisp.h

    r695 r698  
    118118  virtual Bool_t ProcessMessage(Long_t msg, Long_t parm1, Long_t parm2);
    119119 
    120   ClassDef(MGFadcDisp, 1)
     120  ClassDef(MGFadcDisp, 0)
    121121};
    122122

  • trunk/MagicSoft/Mars/mgui/MCamDisplay.h

    r695 r698  
    7373    }
    7474
    75   ClassDef(MCamDisplay, 1) // Display the magic camera
     75  ClassDef(MCamDisplay, 0) // Display the magic camera
    7676};
    7777

  • trunk/MagicSoft/Mars/mgui/MCamNeighbor.h

    r653 r698  
    2020  void Print(Option_t *opt=NULL) ;
    2121
    22   ClassDef(MCamNeighbor, 1) // Conatins the neighbors in the magic camera
     22  ClassDef(MCamNeighbor, 0) // Conatins the neighbors in the magic camera
    2323} ;
    2424

  • trunk/MagicSoft/Mars/mgui/MGDataCheckMain.cc

    r571 r698  
    252252  if ( strcmp ( fInputFile, "\n") == 0 )
    253253    {
    254       return ( kFALSE ) ;
     254      return kFALSE ;
    255255    }
    256256 
    257   return (kTRUE) ;
     257  return kTRUE ;
    258258}
    259259

  • trunk/MagicSoft/Mars/mgui/MGDataCheckMain.h

    r669 r698  
    8686  Bool_t ProcessMessage(Long_t msg, Long_t parm1, Long_t parm2);
    8787
    88   ClassDef(MGDataCheckMain, 1)
     88  ClassDef(MGDataCheckMain, 0)
    8989
    9090} ;

  • trunk/MagicSoft/Mars/mgui/MGEventDispMain.h

    r587 r698  
    9797  Bool_t ProcessMessage(Long_t msg, Long_t parm1, Long_t parm2);
    9898
    99   ClassDef(MGEventDispMain, 1)
     99  ClassDef(MGEventDispMain, 0)
    100100
    101101} ;

  • trunk/MagicSoft/Mars/mgui/MGMarsMain.h

    r466 r698  
    5858  Bool_t ProcessMessage(Long_t msg, Long_t parm1, Long_t parm2);
    5959
    60   ClassDef(MGMarsMain, 1)
     60  ClassDef(MGMarsMain, 0)
    6161} ;
    6262

  • trunk/MagicSoft/Mars/mgui/MGMonteCarloMain.h

    r571 r698  
    8181  Bool_t ProcessMessage(Long_t msg, Long_t parm1, Long_t parm2);
    8282
    83   ClassDef(MGMonteCarloMain, 1)
     83  ClassDef(MGMonteCarloMain, 0)
    8484
    8585} ;

  • trunk/MagicSoft/Mars/mgui/MGPrototyp.h

    r609 r698  
    9393  Bool_t ProcessMessage(Long_t msg, Long_t parm1, Long_t parm2);
    9494
    95   ClassDef(MGPrototyp, 1)
     95  ClassDef(MGPrototyp, 0)
    9696
    9797} ;

  • trunk/MagicSoft/Mars/mgui/MGeomCamCT1.cc

    r695 r698  
    3434    // it was punt into a root/C++ context by Harald Kornmayer and Thomas Bretz
    3535   
    36     const Float_t pixdiameter = 21 ;    // units are cm
     36    const Float_t diameter = 21 ;    // units are mm
     37    const Float_t kS32  = sqrt(3)/2;
    3738
    3839    //
    3940    //  add the first pixel to the list
    4041    //
    41     (*this)[0].Set(0, 0, pixdiameter);
    42     //fPixels->Add( new MGeomPix( ) ;
    43 
    44     const Float_t kS32  = sqrt(3)/2;
    45     const Float_t kPI23 = kPI*2/3;
    46 
    4742    Int_t pixnum = 1;
    4843
    49     for (Int_t ringcounter=1; ringcounter<7; ringcounter++) {
    50         //
    51         // calc. numofpix in ring number i first
    52         //
    53         const Int_t numpixthisring = ringcounter*6;
     44    (*this)[pixnum++].Set(0, 0, diameter);
    5445
     46    for (Int_t ring=1; ring<7; ring++)
     47    {
    5548        //
    5649        // calc. coords for this ring counting from the
    5750        // starting number to the ending number
    5851        //
    59         for (Int_t ipix = 0; ipix < numpixthisring; ipix++) {
     52        for (int i=0; i<ring; i++)
     53            (*this)[pixnum++].Set((-ring+i*0.5)*diameter,
     54                                  i*kS32*diameter,
     55                                  diameter);
    6056
    61             Float_t ang     = 60.0/ringcounter * ipix;
    62             Float_t fracang = ang - 60*(int)floor(ang/60);
     57        for (int i=0; i<ring; i++)
     58            (*this)[pixnum++].Set((-ring*0.5+i)*diameter,
     59                                  ring*kS32 * diameter,
     60                                  diameter);
    6361
    64             ang     /= kRad2Deg;
    65             fracang /= kRad2Deg;
     62        for (int i=0; i<ring; i++)
     63            (*this)[pixnum++].Set((ring+i)*0.5*diameter,
     64                                  (ring-i)*kS32*diameter,
     65                                  diameter);
    6666
    67             Float_t rad  = pixdiameter * ringcounter;
    68                     rad *= sin(kPI23-fracang) * kS32;
     67        for (int i=0; i<ring; i++)
     68            (*this)[pixnum++].Set((ring-0.5*i)*diameter,
     69                                  -i*kS32*diameter,
     70                                  diameter);
    6971
    70             //
    71             //   fill the Data into the array
    72             //
     72        for (int i=0; i<ring; i++)
     73            (*this)[pixnum++].Set((ring*0.5-i)*diameter,
     74                                  -ring*kS32 * diameter,
     75                                  diameter);
    7376
    74             (*this)[pixnum++].Set(rad * cos(ang),
    75                                   rad * sin(ang),
    76                                   pixdiameter);
    77             // fPixels->Add( new MGeomPix( ) ;
    78         }
     77        for (int i=0; i<ring; i++)
     78            (*this)[pixnum++].Set((-ring-i)*0.5*diameter,
     79                                  (-ring+i)*kS32*diameter,
     80                                  diameter);
    7981    }
    8082}

  • trunk/MagicSoft/Mars/mgui/MHexagon.cc

    r653 r698  
    9292}
    9393
    94 
    95 void MHexagon::Draw(Option_t *option )
    96 {
    97   //   Draw this ellipse with its current attributes
    98   AppendPad(option) ;
    99 }
    100 
    101 void MHexagon::DrawHexagon( Float_t x, Float_t y, Float_t d )
     94void MHexagon::DrawHexagon( Float_t x, Float_t y, Float_t d )
    10295{
    10396  //   Draw this ellipse with new coordinate

  • trunk/MagicSoft/Mars/mgui/MHexagon.h

    r665 r698  
    4343
    4444  virtual Int_t  DistancetoPrimitive(Int_t px, Int_t py ) ;
    45   virtual void   Draw( Option_t *Option="") ;
    4645  virtual void   DrawHexagon( Float_t x, Float_t y, Float_t d ) ;
    4746 
    … …  
    5352  virtual void   Print(Option_t *Option="") ;
    5453
    55   ClassDef ( MHexagon, 1 )    //  a hexagon for MAGIC
     54  ClassDef ( MHexagon, 0 )    //  a hexagon for MAGIC
    5655
    5756} ;

  • trunk/MagicSoft/Mars/mmontecarlo/MCollAreaTrigger.h

    r686 r698  
    2424  Bool_t PostProcess() ;
    2525 
    26   ClassDef(MCollAreaTrigger, 1) // Task to fill the collection area histograms
     26  ClassDef(MCollAreaTrigger, 0) // Task to fill the collection area histograms
    2727
    2828};

  • trunk/MagicSoft/Mars/mraw/MRawEvtPixelIter.h

    r664 r698  
    8989    void Draw(Option_t *t="GRAPH");
    9090
    91     ClassDef(MRawEvtPixelIter, 1) // iterates over all pixels of one MRawEvtData object
     91    ClassDef(MRawEvtPixelIter, 0) // iterates over all pixels of one MRawEvtData object
    9292};
    9393

  • trunk/MagicSoft/Mars/mraw/MRawFileRead.h

    r690 r698  
    3333    Bool_t Process();
    3434
    35     ClassDef(MRawFileRead, 1)   // Task to read the raw data binary file
     35    ClassDef(MRawFileRead, 0)   // Task to read the raw data binary file
    3636};
    3737

  • trunk/MagicSoft/Mars/mraw/MRawFileWrite.h

    r454 r698  
    4444    Bool_t PostProcess();
    4545
    46     ClassDef(MRawFileWrite, 1)  // Task to write the raw data containers to a root file
     46    ClassDef(MRawFileWrite, 0)  // Task to write the raw data containers to a root file
    4747};
    4848

  • trunk/MagicSoft/include-Classes/MMcFormat/MMcFadcHeader.hxx

    r693 r698  
    5656      fElecNoise[i]=sigma[i];
    5757  }
     58 
     59  Float_t GetPedestal(UInt_t i) const    { return fPedesMean[i]; }
     60  Float_t GetPedestalRms(UInt_t i) const { return fElecNoise[i]; }
     61
     62  const UInt_t GetNumPixel() const { return MFADC_CHANNELS; }
    5863
    5964  ClassDef(MMcFadcHeader, 1)  //Stores Montecarlo Information
Note: See TracChangeset for help on using the changeset viewer.