source: branches/MarsISDCBranchBasedOn17887/macros/tutorials/readraw.C@ 19085

/* ======================================================================== *\
!
! *
! * This file is part of MARS, the MAGIC Analysis and Reconstruction
! * Software. It is distributed to you in the hope that it can be a useful
! * and timesaving tool in analysing Data of imaging Cerenkov telescopes.
! * It is distributed WITHOUT ANY WARRANTY.
! *
! * Permission to use, copy, modify and distribute this software and its
! * documentation for any purpose is hereby granted without fee,
! * provided that the above copyright notice appear in all copies and
! * that both that copyright notice and this permission notice appear
! * in supporting documentation. It is provided "as is" without express
! * or implied warranty.
! *
!
!
!   Author(s): Thomas Bretz  12/2000 <mailto:tbretz@astro.uni-wuerzburg.de>
!
!   Copyright: MAGIC Software Development, 2000-2004
!
!
\* ======================================================================== */


void readraw(const char *fname="/data/MAGIC/Period016/mcdata/spot_1cm/standard/gamma/Gamma_zbin9_90_7_1740to1749_w0.root")
{
    //
    // open the file
    //
    TFile input(fname, "READ");

    //
    // open the Run Header tree
    //
    TTree *runtree = (TTree*) input.Get("RunHeaders") ;

    if (!runtree)
        return;

    cout << " Entries in Tree RunHeaders: " << runtree->GetEntries() << endl ;

    //
    // create an instance of MRawRunHeader, enable the branch and
    // read the header, print it
    //
    MRawRunHeader *runheader = new MRawRunHeader();
    runtree->GetBranch("MRawRunHeader.")->SetAddress(&runheader);
    runtree->GetEvent(0);
    runheader->Print();

    //
    // open the Data Tree
    //
    TTree *evttree = (TTree*) input.Get("Events");

    //
    // create the instances of the data to read in
    //
    MRawEvtHeader  *evtheader = 0;
    MTime          *evttime   = 0;
    MRawEvtData    *evtdata   = 0;
    MRawCrateArray *evtcrate  = 0;

    //
    // enable the corresponding branches
    //
    evttree->GetBranch("MRawEvtHeader.")->SetAddress(&evtheader);
    evttree->GetBranch("MRawEvtData.")->SetAddress(&evtdata);

    // Use this for real data only
    //evttree->GetBranch("MTime.")->SetAddress(&evttime);
    //evttree->GetBranch("MRawCrateArray.")->SetAddress(&evtcrate);

    //
    // loop over all events and print the information
    //
    Int_t iEnt = evttree->GetEntries();
    cout << " Entries in Tree Data: " << iEnt << endl;

    for (Int_t i = 0; i<iEnt; i++)
    {
      // readin event with the selected branches
      evttree->GetEvent(i);

      evtheader->Print();
      evtdata->Print();

      // Use this for real data only
      //evttime->Print();
      //evtcrate->Print();
    } 
}