CorrSpaceRead.C

Go to the documentation of this file.

/**************************************************************************
 * Copyright(c) 1998, Agata Experiment, All rights reserved. *
 *                                                                        *
 * Author: The Agata Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/


//
//
#if !defined(__CINT__) || defined(__MAKECINT__)

#include <Riostream.h>
#include "TFile.h"
#include "TROOT.h"
#include "TCanvas.h"
#include "TFolder.h"
#include "TTree.h"
#include "TString.h"
#include "TH1D.h"
#include "TStopwatch.h"
#include "CorrelatedSpaceTree2I.h"
#include "CorrelatedSpaceTH2.h"
#include "CubicFilter.h"

#endif

void CorrSpaceTreeRead(Int_t gateInf = 2815, Int_t gateSup = 2824, TString fileName = "test1.root")
{
  
  
  TStopwatch watch;
  Option_t* axis = "x";
  
  TFile* fileIn = new TFile(fileName.Data(), "READ");
  
  watch.Start();
  
  Gw::CorrelatedSpaceTree2I* coor = (Gw::CorrelatedSpaceTree2I*)fileIn->Get("Correlated_Space");
  Gw::CubicFilter* filter = new Gw::CubicFilter();
  filter->SetName("Gate0");
  
  TFolder *folder = gROOT->GetRootFolder()->AddFolder("Filter","Folder");
  folder->Add(filter);  
  
  filter->AddGate(gateInf, gateSup);
  //filter->AddGate(gateInf+5, gateSup+5);
  //filter->ImportGates("gate.dat");
  
  coor->Project(axis);
  
  watch.Stop();
  
  TH1D* h = coor->GetHisto(axis);
  h->Draw();
  printf("Real time %6.3f\n", watch.RealTime());
}

void CorrSpaceH2Read(Int_t gateInf = 2715, Int_t gateSup = 2824, TString fileName = "test3.root")
{
  
  
  TStopwatch watch;
  Option_t* axis = "x";
  
  TFile* fileIn = new TFile(fileName.Data(), "READ");
  
  watch.Start();
  
  Gw::CorrelatedSpaceTH2* coor = (Gw::CorrelatedSpaceTH2*)fileIn->Get("Correlated_Space");
  Gw::CubicFilter* filter = new Gw::CubicFilter();
  filter->SetName("Gate0");
  
  TFolder *folder = gROOT->GetRootFolder()->AddFolder(Gw::CSFilter::GetFolderName(),"Folder");
  folder->Add(filter);
  
  
  filter->AddGate(gateInf, gateSup);
  //filter->AddGate(gateInf+200, gateSup+200);
  //filter->ImportGates("gate.dat");
  
  coor->Project(axis);
  
  watch.Stop();
  
  TH1D* h = coor->GetHisto(axis);
  h->Draw();
  printf("Real time %6.3f\n", watch.RealTime());
}