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)
	if (evt%100==0) cout<< "EVT "<< evt<< endl
	if (!kaltrk) continue

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)
Int_t	refIndex = kaltrk->GetRefIndex()

Function Documentation

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().

if

(

evt%

100 = =0

)

if

(

!

kaltrk

)

Definition at line 4 of file checktubes.C.

References c, exit(), PndSttMapCreator::FillTubeArray(), geoMan, PndSttTube::GetHalfLength(), PndSttTube::GetNeighborings(), PndSttTube::GetPosition(), PndSttTube::GetWireDirection(), gGeoManager, h2, and PndSttTube::IsParallel().

                    {
  file->Get("FairBaseParSet"); 
  TGeoManager *geoMan = gGeoManager; 
  if(!geoMan) { 
    file->Get("FairGeoParSet"); 
    geoMan = gGeoManager; 
    if(!geoMan) { 
      cout << "Could not find valid GeoManager. Abort now!" << endl; 
      exit(1); 
    } 
  }

  PndGeoSttPar *parameters = file.Get("PndGeoSttPar");
  PndSttMapCreator *mapper = new PndSttMapCreator(parameters);
  TClonesArray *fTubeArray = mapper->FillTubeArray();

  TCanvas *c = new TCanvas("c", "", 400, 0, 900, 900);
  TH2F *h2 = new TH2F("h2", "", 100, -42, 42, 100, -42, 42);
  h2->SetStats(kFALSE);
  h2->Draw();

  TArc *arc = NULL;
  TMarker *mrk = NULL;
  PndSttTube *tube = NULL;
  int itube = 0;
  while(itube != -1) {
    c->Clear();
    c->Draw();
    h2->Draw();


    // draw all the tubes
    for(itube = 1; itube < fTubeArray->GetEntriesFast(); itube++) {
      tube = (PndSttTube*) fTubeArray->At(itube);
      if(tube->IsParallel()) {
        arc = new TArc(tube->GetPosition().X(), tube->GetPosition().Y(), 0.5);
        arc->Draw("SAME");
      }
      else {
        mrk = new TMarker(tube->GetPosition().X(), tube->GetPosition().Y(), 6);
        mrk->Draw("SAME");
      }
    }
    // ............................

    
    
    cout << "WHICH tube do you want to test? [to exit type -1]" << endl;
    cin >> itube;
    if(itube == -1) exit(0);
    cout << "TESTING tube " << itube << endl;
    tube = (PndSttTube*) fTubeArray.At(itube);
    if(!tube) {
      cout << "this tube does not exist" << endl; 
      continue; 
    }

    TArrayI neighborings = tube->GetNeighborings();
    cout << "nof neighboring tubes: " << neighborings->GetSize() << endl;
    arc = new TArc(tube->GetPosition().X(), tube->GetPosition().Y(), 0.5);
    arc->SetFillColor(kYellow);
    arc->Draw("SAME");

    for(int inei = 0; inei < neighborings.GetSize(); inei++) { 
      tube = (PndSttTube*) fTubeArray.At(neighborings.At(inei));
      TVector3 position = tube->GetPosition(); 
      cout << " " << neighborings.At(inei);
      if(tube->IsParallel()) {
        arc = new TArc(position.X(), position.Y(), 0.5); 
        arc->SetFillColor(3);
        arc->Draw("SAME");
      }
      else {
        mrk = new TMarker(position.X(), position.Y(), 6);

        TVector3 extr1 = position - tube->GetWireDirection() * tube->GetHalfLength();
        TVector3 extr2 = position + tube->GetWireDirection() * tube->GetHalfLength();
        //      mrk = new TMarker(extr1.X(), extr1.Y(), 6);
        //      mrk = new TMarker(extr2.X(), extr2.Y(), 6);

        mrk->SetMarkerColor(3);
        mrk->Draw("SAME");
      }
    }

    cout << endl;
    c->Update();
    c->Modified();
  }
 
}

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
	)

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.

Int_t refIndex = kaltrk->GetRefIndex()

Definition at line 30 of file checkmomentum.C.

Referenced by PndSttMvdGemTracking::EvaluatePerformances(), PndGemFindHits::FindHits(), PndGemFindHits::FindHits2(), PndSttMvdGemTracking::UpdateMCTrackId(), and PndGemFindClustersTB::WriteClusters().

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.