ill/OfflineWatchers.C

Go to the documentation of this file.

/***************************************************************************
 *   Copyright (C) 2010 by Olivier Stezowski                               *
 *   stezow(AT)ipnl.in2p3.fr                                               *
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 *   This program is distributed in the hope that it will be useful,       *
 *   but WITHOUT ANY WARRANTY; without even the implied warranty of        *
 *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         *
 *   GNU General Public License for more details.                          *
 *                                                                         *
 *   You should have received a copy of the GNU General Public License     *
 *   along with this program; if not, write to the                         *
 *   Free Software Foundation, Inc.,                                       *
 *   59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.             *
 ***************************************************************************/

// ROOT includes
#include "TROOT.h"
#include "TTreeWatchers.h"

// Gw and Watchers
#include "NOMADFile.h"
#include "NOMADFilter.h"
#include "NOMADBuilder.h"

#include "FrameDispatcher.h"
#include "BlackHole.h"

#include "DEmulator.h"

// watchers
#include "ILLWatchers.h" 

using namespace Gw;
using namespace ADF;


void CheckRawData(TFile *file = 0x0, TFile *file_tag = 0x0, Bool_t load_only = false)
{
// to avoid too much warning on output screen
        Gw::LogMessage::SetDebug(1);

// a root file to save specta in case
        TDirectory *MainDir = file, *MainTagDir = file_tag;
        if ( file == 0 ) {
                MainDir = new TFile("RawData.root","UPDATE");
        }
        if ( file_tag == 0 ) {
                MainTagDir = new TFile("RawDataTag.root","UPDATE");
        }       
        
// Emulator is a chain of a NOMADFile, a NOMADFilter and a FrameDispacher to which Watchers are attached 
        DEmulator *emulator = new DEmulator();

        emulator->SetGlobal ( NamedItem<Short_t>("Generation",0) );     

// customization of the emulator
    emulator->Add<NOMADFile>("-","Producer");
    //
    emulator->SetGlobal ( NamedItem<UInt_t>("SizeDataBlockC",aMByte) );
    emulator->Add<NOMADFilter>("Producer","Calib");
    //
    emulator->SetGlobal ( NamedItem<UInt_t>("SizeDataBlockC",5*aMByte) );
 //   emulator.Add<BlackHole>("Calib","Watchers");  
    emulator->Add<FrameDispatcher>("Calib","Watchers");  

// Create all actors in the topology
        emulator->Init() ;

// Customize
        FrameDispatcher *d = 
                emulator->GetActorByName<FrameDispatcher>("Watchers");
        d->SetName("RawData");
        d->SetDirectory(MainDir, MainTagDir);

        FrameTrigger  *trig = new FrameTrigger("T");
        trig->Add(FactoryItem("Default","UID",Version(1,1)),FactoryItem("Default","data",Version(0,0)));
        trig->Universal(true);
        d->SetTrigger(trig); 

        Watcher *w = d->Add<DoEXOILLCalib>("ILLSpecta","Spectra");
        gROOT->GetListOfTasks()->Add(w);

// Run the emulator  
        if ( !load_only ) {
                emulator->Run();
                d->Snapshot("* !tag");
        }       
}


void AnalyseData(TFile *file = 0x0, TFile *file_tag = 0x0, Bool_t load_only = false)
{
// to avoid too much warning on output screen
        Gw::LogMessage::SetDebug(1);
        
// a root file to save specta in case
        TDirectory *MainDir = file, *MainTagDir = file_tag;
        if ( file == 0 ) {
                MainDir = new TFile("Analysis.root","UPDATE");
        }
        if ( file_tag == 0 ) {
                MainTagDir = new TFile("AnalysisTag.root","UPDATE");
        }       
        
// Emulator is a chain of a NOMADFile, a NOMADFilter and a FrameDispacher to which Watchers are attached 
        DEmulator *emulator = new DEmulator();
        
        emulator->SetGlobal ( NamedItem<Short_t>("Generation",0) );     
        
        // customization of the emulator
        emulator->Add<NOMADFile>("-","Producer");
        // emulator->Add<NOMADGene>("-","Producer");    
        //
        emulator->SetGlobal ( NamedItem<UInt_t>("SizeDataBlockC",aMByte) );
        emulator->Add<NOMADBuilder>("Producer","Calib");
        //
        emulator->SetGlobal ( NamedItem<UInt_t>("SizeDataBlockC",5*aMByte) );
        //   emulator.Add<BlackHole>("Calib","Watchers");  
        emulator->Add<FrameDispatcher>("Calib","Watchers");
        // Create all actors in the topology
        emulator->Init() ;
        
// Customize the final consumer
        FrameDispatcher *d = 
        emulator->GetActorByName<FrameDispatcher>("Watchers");
        d->SetName("RawData");
        d->SetDirectory(MainDir, MainTagDir);
        
        FrameTrigger  *trig = new FrameTrigger("T");
        trig->Add(FactoryItem("Default","UID",Version(1,1)),FactoryItem("Default","data",Version(0,1)));
        trig->Universal(true);
        d->SetTrigger(trig); 
        
        // list of watchers added
    // singles
        Watcher *w = d->Add<DoEXOILLCalib>("ILLSpecta","Spectra");
        gROOT->GetListOfTasks()->Add(w);
    // raw matrices
        w = d->Add<RawMatrices>("RawMatrices","Coincidences");
        gROOT->GetListOfTasks()->Add(w);
    // gated matrix
    std::ifstream in("GatedGG.conf"); std::string tmp, file_with_gates;
        if ( in.is_open() == true ) {
        getline(in,tmp);
        tmp += " ";
        while ( in.good() && !in.eof() ) { // read input stream line by line
            
            if ( tmp[0] == '#' ) {
                getline(in,tmp);
                continue;
            } // this line is a comment
            
            std::istringstream decode(tmp); decode >> file_with_gates ;
            w = d->Add<GatedGG>(file_with_gates.data(),"Gated Gamma Gamma Matrix");
            gROOT->GetListOfTasks()->Add(w);
            
            getline(in,tmp);
            tmp += " ";
        }
        in.close();
        }
    
// Run the emulator  
        if ( !load_only )
                emulator->Run();

// save spectra produced
        d->Snapshot();
}

void DoMyTree()
{
        // to avoid too much warning on output screen
        Gw::LogMessage::SetDebug(1);
        
        // Emulator is a chain of a NOMADFile, a NOMADFilter and a FrameDispacher to which Watchers are attached 
        DEmulator *emulator = 
                new DEmulator();
        emulator->SetGlobal ( NamedItem<Short_t>("Generation",0) );     
        
        // customization of the emulator
        emulator->Add<NOMADFile>("-","Producer");       
        emulator->SetGlobal ( NamedItem<UInt_t>("SizeDataBlockC",aMByte) );
        emulator->Add<NOMADBuilder>("Producer","Calib");
        emulator->SetGlobal ( NamedItem<UInt_t>("SizeDataBlockC",5*aMByte) );
        emulator->Add<FrameDispatcher>("Calib","Watchers");  
        
        // Create all actors in the topology
        emulator->Init() ;
        
        // Customize
        FrameDispatcher *d = 
                emulator->GetActorByName<FrameDispatcher>("Watchers");
        d->SetName("RawData");
        
        FrameTrigger *trig = new FrameTrigger("T");
        trig->Add(FactoryItem("Default","UID",Version(1,1)),FactoryItem("Default","data",Version(0,1)));
        trig->Universal(true);
        d->SetTrigger(trig); 
        
        
        TTreeMaster *rootmaster = d->Add<TTreeMaster>("ILLTree","Own the tree");
        rootmaster->Add<MyILLTree>("EXO","Branch on EXOGAM DATA");
        
        // Run the emulator  
        emulator->Run();
        
        delete emulator;        
}


void CheckReading(TFile *file = 0x0)
{
        // to avoid too much warning on output screen
        Gw::LogMessage::SetDebug(1);
        
        // a root file to save specta in case
        TDirectory *MainDir = file;
        if ( file == 0 ) 
                MainDir = new TFile("Offline.root","UPDATE");
        
        // Emulator is a chain of a NOMADFile, a NOMADFilter and a FrameDispacher to which Watchers are attached 
        DEmulator *emulator = new DEmulator();
        
        emulator->SetGlobal ( NamedItem<Short_t>("Generation",0) );     
        
        // customization of the emulator
        emulator->Add<NOMADFile>("-","Producer");
        emulator->SetGlobal ( NamedItem<UInt_t>("SizeDataBlockC",aMByte) );
        emulator->Add<NOMADFilter>("Producer","Calib");
        emulator->SetGlobal ( NamedItem<UInt_t>("SizeDataBlockC",5*aMByte) );
        //   emulator.Add<BlackHole>("Calib","Watchers");  
//      emulator->Add<FrameDispatcher>("Calib","Watchers");  
        emulator->Add<BlackHole>("Calib","Watchers");
        
        // Run the emulator        
        emulator->Init() ;
        emulator->Run();
        
        
}

#if !defined(__CINT__) 
int main(int /*argc*/, char ** /*argv*/)
{               
        // a new root file is created
        TFile f("Offline.root","RECREATE");
        // run your analysis on it
        CheckRawData(&f);
        
        return 1;
}
#endif