digi_analys.C File Reference

Go to the source code of this file.

Functions
gROOT	LoadMacro ("$VMCWORKDIR/gconfig/basiclibs.C")
	rootlogon ()
	basiclibs ()
t	SetBranchAddress ("EmcDigi",&digi_array)
	for (Int_t j=0;j< t->GetEntriesFast();j++)
h1	Draw ()
c1	Update ()
	if (0)

Variables
TFile *	f = new TFile("digi_emc.root")
TTree *	t =(TTree *) f->Get("pndsim")
TClonesArray *	digi_array =new TClonesArray("EmcDigi")
double	digi_energy
double	hit_energy
TH1F *	h1 = new TH1F("h1","Energy ratio",100,0.,10.)
TCanvas *	c1 = new TCanvas("c1", "", 100, 100, 800, 800)

Function Documentation

basiclibs

(

)

h1 Draw

(

)

Definition at line 64 of file NHitsPerEvent.C.

           {

        std::cout << "mvdHitArray: " << mvdHitArray->GetEntries() << " projectedHitArray " << mvdProjectedArray->GetEntries() << std::endl;

//      gHitPoints->Clear();
        ExtractData(mvdHitArray, gHitPoints);
        ExtractData(mvdProjectedArray, gProjections);
        //if (gHitPoints->GetN() != x.size()) gHitPoints->Expand(x.size());

        std::cout << "NHits: " << gHitPoints->GetN() << std::endl;
        gHitPoints->Draw("apl");
        gProjections->Draw("pl");
}

for

(

Int_t

j = 0; j< t->GetEntriesFast(); j++

)

Definition at line 19 of file digi_analys.C.

References digi, and i.

        {
                t->GetEntry(j);
                for (Int_t i=0; i<digi_array->GetEntriesFast(); i++)
                {
                        EmcDigi *digi=(EmcDigi*)digi_array->At(i);
                        digi_energy=digi->GetEnergy();
                        std::vector<CbmEmcHit*> hitList=digi->itsHitList();
                        if (hitList.size()!=0)
                        {
                                hit_energy=hitList[0]->GetEnergy();
                                if (hit_energy>1e-3) h1->Fill(digi_energy/hit_energy);
                        }
                }
        }

if

(

0

)

Definition at line 57 of file digi_analys.C.

References c2, h2, and i.

        {
        TClonesArray *wfa=new TClonesArray("EmcWaveform");
        t->SetBranchAddress("EmcWaveform",&wfa);

        t->GetEntry(0);

        char ch;

        TH1F *h2= new TH1F("h2","Waveform",64,1.,64.);
        TCanvas* c2 = new TCanvas("c2", "", 100, 100, 800, 800);
        for (Int_t i=0; i<wfa->GetEntriesFast(); i++)
        {
                h2->Reset();
                EmcWaveform *wf=(EmcWaveform*)wfa->At(i);
                vector<double> signal = wf->GetSignal();

                for (Int_t j=0;j<signal.size();j++)
                {
                        h2->SetBinContent(j,signal[j]);
                        j++;
                }

                h2->Draw();
                c2->Update();
                gPad->WaitPrimitive();
                cin.get(ch);
        }
        }

gROOT LoadMacro

(

"$VMCWORKDIR/gconfig/basiclibs.C"

)

rootlogon

(

)

Definition at line 1 of file outdated/mpiTools/macros/emc/rootlogon.C.

{
  gSystem->AddIncludePath("-I$VMCWORKDIR/emc");
  gSystem->AddIncludePath("-I$VMCWORKDIR/emc/EmcDigi");
  gSystem->AddIncludePath("-I$VMCWORKDIR/emc/EmcMC");
  gSystem->AddIncludePath("-I$VMCWORKDIR/emc/EmcTools");
  gSystem->AddIncludePath("-I$VMCWORKDIR/emc/EmcReco");
  gSystem->AddIncludePath("-I$VMCWORKDIR/emc/EmcData");
  gSystem->AddIncludePath("-I$VMCWORKDIR/pnddata");
  gSystem->AddIncludePath("-I$VMCWORKDIR/base");
  gROOT->LoadMacro("$VMCWORKDIR/gconfig/basiclibs.C");
  basiclibs();

  gSystem->Load("libDpmEvtGen.so");
  gSystem->Load("libGeoBase");
  gSystem->Load("libParBase");
  gSystem->Load("libBase");
  gSystem->Load("libTrkBase");
  gSystem->Load("libPndData");
  gSystem->Load("libField");
  gSystem->Load("libPassive");
  gSystem->Load("libGen");
  gSystem->Load("libPGen");
  gSystem->Load("libEmc"); 

  delete gRandom;
  gRandom=new TRandom3(0);
}

t SetBranchAddress	(	"EmcDigi"	,
		&	digi_array
	)

c1 Update

(

)

Variable Documentation

TCanvas* c1 = new TCanvas("c1", "", 100, 100, 800, 800)

Definition at line 35 of file digi_analys.C.

TClonesArray* digi_array =new TClonesArray("EmcDigi")

Definition at line 13 of file digi_analys.C.

double digi_energy

Definition at line 17 of file digi_analys.C.

Referenced by reco_analys().

TFile* f = new TFile("digi_emc.root")

Definition at line 10 of file digi_analys.C.

TH1F* h1 = new TH1F("h1","Energy ratio",100,0.,10.)

Definition at line 18 of file digi_analys.C.

double hit_energy

Definition at line 17 of file digi_analys.C.

Referenced by for(), and reco_analys().

TTree* t =(TTree *) f->Get("pndsim")

Definition at line 11 of file digi_analys.C.