Logo ROOT  
Reference Guide
 
Loading...
Searching...
No Matches
seism.C
Go to the documentation of this file.
1/// \file
2/// \ingroup tutorial_graphs
3/// \notebook -js
4/// Strip chart example.
5///
6/// \macro_image
7/// \macro_code
8///
9/// \author Rene Brun
10
11#include "TStopwatch.h"
12#include "TDatime.h"
13#include "TStyle.h"
14#include "TH1F.h"
15#include "TCanvas.h"
16#include "TSystem.h"
17#include "TRandom.h"
18#include <cstdio>
19
20void seism() {
21
22 TStopwatch sw; sw.Start();
23 //set time offset
24 TDatime dtime;
26
27 TCanvas *c1 = new TCanvas("c1","Time on axis",10,10,1000,500);
28 c1->SetGrid();
29
30 Float_t bintime = 1; //one bin = 1 second. change it to set the time scale
31 TH1F *ht = new TH1F("ht","The ROOT seism",10,0,10*bintime);
32 Float_t signalval = 1000;
33 ht->SetMaximum( signalval);
34 ht->SetMinimum(-signalval);
35 ht->SetStats(false);
36 ht->SetLineColor(2);
37 ht->GetXaxis()->SetTimeDisplay(1);
38 ht->GetYaxis()->SetNdivisions(520);
39 ht->Draw();
40
41 for (Int_t i=1;i<2300;i++) {
42 //======= Build a signal : noisy damped sine ======
43 Float_t noise = gRandom->Gaus(0,120);
44 if (i > 700) noise += signalval*sin((i-700.)*6.28/30)*exp((700.-i)/300.);
45 ht->SetBinContent(i,noise);
46 c1->Modified();
47 c1->Update();
48 gSystem->ProcessEvents(); //canvas can be edited during the loop
49 }
50 printf("Real Time = %8.3fs, Cpu Time = %8.3fs\n",sw.RealTime(),sw.CpuTime());
51}
int Int_t
Definition RtypesCore.h:45
float Float_t
Definition RtypesCore.h:57
R__EXTERN TRandom * gRandom
Definition TRandom.h:62
R__EXTERN TStyle * gStyle
Definition TStyle.h:433
R__EXTERN TSystem * gSystem
Definition TSystem.h:560
virtual void SetNdivisions(Int_t n=510, Bool_t optim=kTRUE)
Set the number of divisions for this axis.
Definition TAttAxis.cxx:233
virtual void SetLineColor(Color_t lcolor)
Set the line color.
Definition TAttLine.h:40
virtual void SetTimeDisplay(Int_t value)
Definition TAxis.h:171
The Canvas class.
Definition TCanvas.h:23
This class stores the date and time with a precision of one second in an unsigned 32 bit word (950130...
Definition TDatime.h:37
UInt_t Convert(Bool_t toGMT=kFALSE) const
Convert fDatime from TDatime format to the standard time_t format.
Definition TDatime.cxx:182
1-D histogram with a float per channel (see TH1 documentation)}
Definition TH1.h:577
TAxis * GetXaxis()
Definition TH1.h:322
virtual void SetMaximum(Double_t maximum=-1111)
Definition TH1.h:400
TAxis * GetYaxis()
Definition TH1.h:323
void Draw(Option_t *option="") override
Draw this histogram with options.
Definition TH1.cxx:3067
virtual void SetMinimum(Double_t minimum=-1111)
Definition TH1.h:401
virtual void SetBinContent(Int_t bin, Double_t content)
Set bin content see convention for numbering bins in TH1::GetBin In case the bin number is greater th...
Definition TH1.cxx:9139
virtual void SetStats(Bool_t stats=kTRUE)
Set statistics option on/off.
Definition TH1.cxx:8907
virtual Double_t Gaus(Double_t mean=0, Double_t sigma=1)
Samples a random number from the standard Normal (Gaussian) Distribution with the given mean and sigm...
Definition TRandom.cxx:274
Stopwatch class.
Definition TStopwatch.h:28
Double_t RealTime()
Stop the stopwatch (if it is running) and return the realtime (in seconds) passed between the start a...
void Start(Bool_t reset=kTRUE)
Start the stopwatch.
Double_t CpuTime()
Stop the stopwatch (if it is running) and return the cputime (in seconds) passed between the start an...
void SetTimeOffset(Double_t toffset)
Change the time offset for time plotting.
Definition TStyle.cxx:1909
virtual Bool_t ProcessEvents()
Process pending events (GUI, timers, sockets).
Definition TSystem.cxx:403
RVec< PromoteType< T > > exp(const RVec< T > &v)
Definition RVec.hxx:1800
RVec< PromoteType< T > > sin(const RVec< T > &v)
Definition RVec.hxx:1814
return c1
Definition legend1.C:41