checkmomentum.C File Reference

Go to the source code of this file.

Functions
timer	Start ()
gROOT	LoadMacro ("$VMCWORKDIR/gconfig/rootlogon.C")
	rootlogon ()
gSystem	Load ("libSttMvdTracking")
TFile	filereco ("MvdStt_Test_reco.root")
treereco	SetBranchAddress ("SttMvdTrack",&prefit)
treereco	SetBranchAddress ("SttMvdKalmanTrack",&kalfit)
	for (int evt=0;evt< treereco->GetEntriesFast();evt++)
hkalmom	SetLineColor (4)
hpremom	SetLineColor (2)
hkalmom	Draw ()
hpremom	Draw ("SAME")

Variables
TStopwatch	timer
TTree *	treereco = (TTree*) filereco.Get("pndsim")
TClonesArray *	prefit = new TClonesArray("PndTrack")
TClonesArray *	kalfit = new TClonesArray("PndTrack")
TH1F *	hkalmom = new TH1F("hkalmom", "genfit", 100, 0.5, 1.5)
TH1F *	hpremom = new TH1F("hpremom", "prefit", 100, 0.5, 1.5)

Function Documentation

hkalmom 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");
}

hpremom Draw

(

"SAME"

)

TFile filereco

(

"MvdStt_Test_reco.root"

)

Referenced by QAmacro_stt_4(), runDetectorInfo(), runEventInfo(), runOnline(), runOnlineDisplay(), runOnlineDisplayGeomTest(), runOnlineDisplayMCCheck(), runOnlineDisplayMCCheckFaster(), runOnlineDisplayMCCheckFaster2(), runOnlineDisplayMCCheckFaster3(), runOnlineDisplayProduction(), runOnlineDisplayRelease(), runStreamDisplay(), runTimestampCheck(), runTripletAna(), and runTripletFinderMini().

for

(

int

evt = 0; evt < treereco->GetEntriesFast(); evt++

)

Definition at line 20 of file checkmomentum.C.

References evt, PndTrack::GetFlag(), PndTrack::GetParamFirst(), and PndTrack::GetRefIndex().

    {
      treereco->GetEntry(evt);
      if(evt%100 == 0) cout << "EVT " << evt << endl;

      for(int ikal = 0; ikal < kalfit->GetEntriesFast(); ikal++)
        {
          PndTrack* kaltrk = (PndTrack*) kalfit->At(ikal);
          if(!kaltrk) continue;

          Int_t refIndex = kaltrk->GetRefIndex();
          if(refIndex == -1) continue;
          PndTrack *pretrk = (PndTrack*) prefit->At(refIndex);
          if(!pretrk) continue;

          if(kaltrk->GetFlag() < 0) continue;

          // FILL histos only when both kalman track
          // and prefit track are present and kalman
          // flag is not negative
          TVector3 kalmom = kaltrk->GetParamFirst().GetMomentum();
          hkalmom->Fill(kalmom.Mag());
          TVector3 premom = pretrk->GetParamFirst().GetMomentum();
          hpremom->Fill(premom.Mag());
        }
    }

gSystem Load

(

"libSttMvdTracking"

)

gROOT LoadMacro

(

"$VMCWORKDIR/gconfig/rootlogon.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);
}

treereco SetBranchAddress	(	"SttMvdTrack"	,
		&	prefit
	)

treereco SetBranchAddress	(	"SttMvdKalmanTrack"	,
		&	kalfit
	)

hkalmom SetLineColor

(

4

)

hpremom SetLineColor

(

2

)

timer Start

(

)

Variable Documentation

TH1F* hkalmom = new TH1F("hkalmom", "genfit", 100, 0.5, 1.5)

Definition at line 17 of file checkmomentum.C.

TH1F* hpremom = new TH1F("hpremom", "prefit", 100, 0.5, 1.5)

Definition at line 18 of file checkmomentum.C.

TClonesArray* kalfit = new TClonesArray("PndTrack")

Definition at line 14 of file checkmomentum.C.

TClonesArray* prefit = new TClonesArray("PndTrack")

Definition at line 12 of file checkmomentum.C.

TStopwatch timer

Initial value:

{

  gROOT->Reset()

Definition at line 1 of file checkmomentum.C.

TTree* treereco = (TTree*) filereco.Get("pndsim")

Definition at line 11 of file checkmomentum.C.