src/PtFlowHe4View.cc

#include "PtFlowHe4View.hh"
#include "LogViewer.hh"
#include <TROOT.h>
#include <TGraph.h>

#include <iostream>
#include <iomanip>

using namespace std;

PtFlowHe4View::PtFlowHe4View() :
  PtView( LogViewer::tgFlwHe4 ),
  sep_( NULL ), eva_( NULL ), sepDo_( NULL ), precool_( NULL )
{
  
  this->title() = "Flow (He4)";
  this->add( "separator",               "Sep.",      kBlue );
  this->add( "evaporator",              "Eva.",      kMagenta );
  this->add( "separator_do",            "Sep. (DO)", kCyan );
  this->add( "precooling",              "Precool.",  kRed );
  
  sep_     = new TGraph;
  eva_     = new TGraph;
  sepDo_   = new TGraph;
  precool_ = new TGraph;
  
}

PtFlowHe4View::~PtFlowHe4View(){
}

void PtFlowHe4View::analysis(){

  
}

void PtFlowHe4View::average( TGraph* org, TGraph* ave, 
                             const double& interval ){
  
  ave->SetLineWidth( org->GetLineWidth() );
  ave->SetLineStyle( org->GetLineStyle() );
  ave->SetLineColor( org->GetLineColor() );
  
  double start( -100.0 );
  double volume = 0.0;
  for( int i = 0; i < org->GetN(); i++ ){
    double x, y;
    org->GetPoint( i, x, y );
    volume += y;
    if( start == -100.0 ) start = x;
    if( x - start < interval ) continue;
    
    volume /= ( x - start );
    int n = ave->GetN();
    ave->Set( n + 1 );
    ave->SetPoint( n, 0.5 * ( x + start ), volume );
    start  = -100.0;
    volume = 0.0;
  }
}

void PtFlowHe4View::setRange(){
  
  TGraph& se = (*this)[ 2 ].graph();
  TGraph& ev = (*this)[ 3 ].graph();
  TGraph& sd = (*this)[ 4 ].graph();
  TGraph& pr = (*this)[ 5 ].graph();

  this->average( &se, sep_ );
  this->average( &ev, eva_ );
  this->average( &sd, sepDo_ );
  this->average( &pr, precool_ );
  
  vector< TGraph* > gs;
  gs.push_back( sep_ );
  gs.push_back( eva_ );
  gs.push_back( sepDo_ );
  gs.push_back( precool_ );

  double x,y( 0.0 );
  for( int i = 0; i < gs.size(); i++ ){
    for( int j = 0; j < gs[ i ]->GetN(); j++ ){
      gs[ i ]->GetPoint( j, x, y );
      ymin_ = ( y < ymin_ ? y : ymin_ );
      ymax_ = ( y > ymax_ ? y : ymax_ );
    }
  }
  
}

void PtFlowHe4View::drawGraph(){

  sep_->Draw( "L" );
  eva_->Draw( "L" );
  sepDo_->Draw( "L" );
  precool_->Draw( "L" );

}