checkdedx_helixhit.C File Reference

#include <vector>

Go to the source code of this file.

Functions
timer	Start ()
TFile	filerun ("testrun.root")
treemc	SetBranchAddress ("STTPoint",&pnt)
treemc	SetBranchAddress ("MCTrack",&mctrack)
TFile	filedigi ("testdigi.root")
treedigi	SetBranchAddress ("STTHit",&digi)
TFile	filehelix ("testreco.root")
treereco	SetBranchAddress ("SttHelixHit",&hh)
treereco	SetBranchAddress ("STTTrack",&track)
treereco	SetBranchAddress ("STTTrackMatch",&matchtrack)
c	Divide (1, 2)
c1	Divide (2, 2)
cout<< treereco->	GetEntriesFast ()
c	cd (1)
hdedxvsp	Draw ()
c	cd (2)
hdedxvsp_e	SetMarkerColor (2)
hdedxvsp_k	SetMarkerColor (3)
c1	cd (3)
hdedxvsp_pi	SetMarkerColor (4)
hdedxvsp_mu	SetMarkerColor (7)
hdedxvsp_mu	Draw ("SAME")
c1	cd (4)
hdedxvsp_p	SetMarkerColor (5)
c2	cd ()
hdedxvsp	Write ()
out_dedx	Close ()

Variables
TStopwatch	timer
TTree *	treemc = (TTree*) filerun.Get("pndsim")
TClonesArray *	pnt = new TClonesArray("PndSttPoint")
TClonesArray *	mctrack = new TClonesArray("PndMCTrack")
TTree *	treedigi = (TTree*) filedigi.Get("pndsim")
TClonesArray *	digi = new TClonesArray("PndSttHit")
TTree *	treereco = (TTree*) filehelix.Get("pndsim")
TClonesArray *	hh = new TClonesArray("PndSttHelixHit")
TClonesArray *	track = new TClonesArray("PndSttTrack")
TClonesArray *	matchtrack = new TClonesArray("PndSttTrackMatch")
TCanvas *	c = new TCanvas("c", "c", 0, 0, 600, 600)
TH1F *	samplenum = new TH1F("samplenum","number of sampling",50,0.,50.)
TH2F *	hdedxvsp = new TH2F("hdedxvsp","dedx vs p",100,0.,1.5, 100,0.,40.)
TCanvas *	c1 = new TCanvas("c1", "c1", 0, 0, 600, 600)
TH2F *	hdedxvsp_p = new TH2F("hdedxvsp_p","dedx vs p for p",100,0.,1.5, 100,0.,40.)
TH2F *	hdedxvsp_pi = new TH2F("hdedxvsp_pi","dedx vs p for #pi",100,0.,1.5, 100,0.,40.)
TH2F *	hdedxvsp_e = new TH2F("hdedxvsp_e","dedx vs p for e",100,0.,1.5, 100,0.,40.)
TH2F *	hdedxvsp_k = new TH2F("hdedxvsp_k","dedx vs p for K",100,0.,1.5, 100,0.,40.)
TH2F *	hdedxvsp_mu = new TH2F("hdedxvsp_mu","dedx vs p for #mu",100,0.,1.5, 100,0.,40.)
TCanvas *	c2 = new TCanvas("c2", "c2", 0, 0, 600, 600)
TH2F *	hdedxvsp_reco = new TH2F("hdedxvsp_reco","dedx vs reco p",100,0.,1.5, 100,0.,40.)
TFile *	out_dedx = new TFile("out_dedx_helixhit.root", "RECREATE")

Function Documentation

c cd

(

1

)

c cd

(

2

)

c1 cd

(

3

)

c1 cd

(

4

)

c2 cd

(

)

out_dedx Close

(

)

c Divide	(	1	,
		2
	)

c1 Divide	(	2	,
		2
	)

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

fdedx5 Draw

(

"SAME"

)

TFile filedigi

(

"testdigi.root"

)

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

TFile filehelix

(

"testreco.root"

)

TFile filerun

(

"testrun.root"

)

Referenced by QAmacro_stt_4().

cout<< treereco-> GetEntriesFast

(

)

Definition at line 51 of file checkdedx_helixhit.C.

References cos(), d0, Double_t, evt, PndSttTrack::GetCharge(), PndSttHelixHit::GetdEdx(), PndSttTrack::GetDist(), PndSttTrack::GetHelixHitIndex(), PndSttHelixHit::GetHitIndex(), PndSttTrackMatch::GetMCTrackID(), PndMCTrack::GetMomentum(), PndSttTrack::GetNofHelixHits(), PndMCTrack::GetPdgCode(), PndSttTrack::GetPhi(), PndSttTrack::GetRad(), PndSttTrack::GetTanL(), PndSttTrack::GetZ(), h, m, phi0, plong, ptot, ptran, R, sin(), sqrt(), tanl, x0, y0, and z0.

                                                             {
    if(evt%100 == 0) cout << evt << endl;

    treemc->GetEntry(evt);
    treedigi->GetEntry(evt);
    treereco->GetEntry(evt);

    // tracks loop
    for (Int_t k = 0; k < track->GetEntriesFast(); k++) {

      PndSttTrack *stttrack = (PndSttTrack*) track->At(k);
      if(!stttrack) continue;
      PndSttTrackMatch *mtrack = (PndSttTrackMatch *) matchtrack->At(k);
      if(!mtrack) continue;
      PndMCTrack *MCtrack = (PndMCTrack*) mctrack->At(mtrack->GetMCTrackID());
      if(!MCtrack) continue;
      Int_t PDGcode = MCtrack->GetPdgCode();

      Double_t d0  = stttrack->GetDist();
      Double_t phi0 = stttrack->GetPhi();
      Double_t R = stttrack->GetRad();
      Double_t z0 = stttrack->GetZ();
      Double_t tanl = stttrack->GetTanL();
      Double_t h = -(Int_t) stttrack->GetCharge();
      Double_t ptran = 0.003 * 2 * R;

      Double_t plong = ptran * tanl;
      Double_t ptot = sqrt(plong*plong + ptran*ptran);

      Double_t x0 = (d0 + R) * cos(phi0);
      Double_t y0 = (d0 + R) * sin(phi0);

      Int_t hitcounter = stttrack->GetNofHelixHits();
      std::vector<double> dedxvec;
      dedxvec.clear();
      TVector3 momentum = MCtrack->GetMomentum();
      Int_t losthit = 0;
      for (Int_t j = 0; j < hitcounter; j++) {
        Int_t iHHit = stttrack->GetHelixHitIndex(j);
        PndSttHelixHit *helixhit = (PndSttHelixHit*) hh->At(iHHit);
        if(!helixhit) continue;

        Int_t hitindex = helixhit->GetHitIndex();
//      PndSttHit* hit = (PndSttHit*) digi->At(hitindex);
//      PndSttPoint *point = (PndSttPoint*) pnt->At(hit->GetRefIndex());

        if(helixhit->GetdEdx() != 0) dedxvec.push_back(helixhit->GetdEdx());
        else losthit++;
      }

      hitcounter -= losthit;

      if(hitcounter > 0) {
        // truncated mean
        Double_t perc = 0.60;
        // sort
        std::sort(dedxvec.begin(), dedxvec.end());

        //truncated mean
        Double_t sum = 0;
        Int_t endnum = floor(hitcounter * perc);


        for(Int_t m = 0; m < endnum; m++) sum += dedxvec[m];
        Double_t tmean;
        if(endnum > 0) {
          tmean = sum/(Double_t) endnum;
          hdedxvsp->Fill(momentum->Mag(), tmean);
          hdedxvsp_reco->Fill(ptot, tmean);

          samplenum->Fill(hitcounter);

          if(abs(PDGcode) == 11)       { hdedxvsp_e->Fill(momentum->Mag(), tmean);}
          else if(abs(PDGcode) == 13)  { hdedxvsp_mu->Fill(momentum->Mag(), tmean);}
          else if(abs(PDGcode) == 211) { hdedxvsp_pi->Fill(momentum->Mag(), tmean);}
          else if(abs(PDGcode) == 321) { hdedxvsp_k->Fill(momentum->Mag(), tmean); }
          else if(abs(PDGcode) == 2212){ hdedxvsp_p->Fill(momentum->Mag(), tmean);}

        }
      }
    }
  }

treemc SetBranchAddress	(	"STTPoint"	,
		&	pnt
	)

treemc SetBranchAddress	(	"MCTrack"	,
		&	mctrack
	)

treedigi SetBranchAddress	(	"STTHit"	,
		&	digi
	)

treereco SetBranchAddress	(	"SttHelixHit"	,
		&	hh
	)

treereco SetBranchAddress	(	"STTTrack"	,
		&	track
	)

treereco SetBranchAddress	(	"STTTrackMatch"	,
		&	matchtrack
	)

hdedxvsp_e SetMarkerColor

(

2

)

hdedxvsp_k SetMarkerColor

(

3

)

hdedxvsp_pi SetMarkerColor

(

4

)

hdedxvsp_mu SetMarkerColor

(

7

)

hdedxvsp_p SetMarkerColor

(

5

)

timer Start

(

)

hdedxvsp Write

(

)

Variable Documentation

TCanvas* c = new TCanvas("c", "c", 0, 0, 600, 600)

Definition at line 35 of file checkdedx_helixhit.C.

TCanvas* c1 = new TCanvas("c1", "c1", 0, 0, 600, 600)

Definition at line 40 of file checkdedx_helixhit.C.

TCanvas* c2 = new TCanvas("c2", "c2", 0, 0, 600, 600)

Definition at line 48 of file checkdedx_helixhit.C.

TClonesArray* digi = new TClonesArray("PndSttHit")

Definition at line 20 of file checkdedx_helixhit.C.

Referenced by PndDrcDigiTask::ActivatePixel(), PndDrcTimeDigiTask::ActivatePixel(), PndEmcCluster::addDigi(), PndMasterRunAna::AddDigiOnlyTasks(), PndMasterRunAna::AddDigiTasks(), PndGemFindClustersTB::AddDigiToCluster(), PndGemFindClusters::AddDigiToCluster(), PndEmcCluster::addLocalMax(), PndGemFindClustersTB::AnalyzeClusters(), PndEmcClusterMoments::Centre1(), PndEmcClusterMoments::Centre2(), PndGemFindClustersTB::CheckCluster(), PndGemFindClustersTB::CompareDigiToClusters(), PndGemFindClusters::CompareDigiToClustersDigis(), PndGemFindClusters::CreateClusters(), PndGemFindClustersTB::CreateClusters(), PndSdsChargeWeightingAlgorithms::DigiCharge(), PndSdsIdealChargeConversion::DigiValueToCharge(), PndEmcClusterEnergySums::E9(), PndEmcClusterProperties::Energy(), PndEmcClusterEnergySums::energy(), PndMQSdsHitLoader::Exec(), for(), PndEmcDigiCaloDraw::GetEnergy(), PndEmcDigiCaloDraw::GetVector(), PndEmcClusterProperties::GravWhere(), PndEmcXClMoments::Init(), PndEmc2DLocMaxFinder::isALocalMax(), PndEmcClusterMoments::Major1(), PndEmcClusterMoments::Major2(), PndEmcClusterProperties::Major_axis(), PndEmcClusterProperties::Mass(), PndEmcCluster::Maxima(), PndEmcClusterMoments::Minor1(), PndEmcClusterMoments::Minor2(), PndEmcClusterMoments::Phi1(), PndEmcClusterMoments::Phi2(), PndEmcXClMoments::Print(), PndMvdPastaTorino::ReadTorinoDaqFrame(), reco_analys(), PndEmcCluster::removeDigi(), PndEmcClusterMoments::SecondMoment(), PndEmcClusterMoments::SecondMomentTP(), PndGemFindHits::SortDigis(), PndEmcClusterMoments::Theta1(), PndEmcClusterMoments::Theta2(), timeStructurePasta(), PndGemFindClustersTB::WriteClusters(), and PndCATracking::WriteMVDHits().

TH2F* hdedxvsp = new TH2F("hdedxvsp","dedx vs p",100,0.,1.5, 100,0.,40.)

Definition at line 38 of file checkdedx_helixhit.C.

TH2F* hdedxvsp_e = new TH2F("hdedxvsp_e","dedx vs p for e",100,0.,1.5, 100,0.,40.)

Definition at line 44 of file checkdedx_helixhit.C.

TH2F* hdedxvsp_k = new TH2F("hdedxvsp_k","dedx vs p for K",100,0.,1.5, 100,0.,40.)

Definition at line 45 of file checkdedx_helixhit.C.

TH2F* hdedxvsp_mu = new TH2F("hdedxvsp_mu","dedx vs p for #mu",100,0.,1.5, 100,0.,40.)

Definition at line 46 of file checkdedx_helixhit.C.

TH2F* hdedxvsp_p = new TH2F("hdedxvsp_p","dedx vs p for p",100,0.,1.5, 100,0.,40.)

Definition at line 42 of file checkdedx_helixhit.C.

TH2F* hdedxvsp_pi = new TH2F("hdedxvsp_pi","dedx vs p for #pi",100,0.,1.5, 100,0.,40.)

Definition at line 43 of file checkdedx_helixhit.C.

TH2F* hdedxvsp_reco = new TH2F("hdedxvsp_reco","dedx vs reco p",100,0.,1.5, 100,0.,40.)

Definition at line 49 of file checkdedx_helixhit.C.

TClonesArray* hh = new TClonesArray("PndSttHelixHit")

Definition at line 26 of file checkdedx_helixhit.C.

Referenced by PndDrcReco::CanvasSave(), and PndFTSCAGBTracker::FitTracks().

TClonesArray* matchtrack = new TClonesArray("PndSttTrackMatch")

Definition at line 32 of file checkdedx_helixhit.C.

TClonesArray* mctrack = new TClonesArray("PndMCTrack")

Definition at line 14 of file checkdedx_helixhit.C.

TFile* out_dedx = new TFile("out_dedx_helixhit.root", "RECREATE")

Definition at line 160 of file checkdedx_helixhit.C.

TClonesArray* pnt = new TClonesArray("PndSttPoint")

Definition at line 12 of file checkdedx_helixhit.C.

Referenced by RhoSimpleVertexSelector::Accept(), PndSttMvdGemTracking::EvaluatePerformances(), PndRichRecoTask::Exec(), PndRichHitProducer::Exec(), PndFtsTrackerIdeal::Exec(), PndLmdBPtestTask::Exec(), PndTrkCombiLegendreTask::Exec(), RhoCalculationTools::GetBz(), PndPidCorrelator::GetTrackInfo(), PndFsmTrack::HelixRep(), PndRichReco::Init(), PndFsmTrack::Propagate(), and QAmacro_stt_4().

TH1F* samplenum = new TH1F("samplenum","number of sampling",50,0.,50.)

Definition at line 37 of file checkdedx_helixhit.C.

TStopwatch timer

Initial value:

{
  gROOT->Reset()

Definition at line 1 of file checkdedx_helixhit.C.

TClonesArray* track = new TClonesArray("PndSttTrack")

Definition at line 30 of file checkdedx_helixhit.C.

TTree* treedigi = (TTree*) filedigi.Get("pndsim")

Definition at line 19 of file checkdedx_helixhit.C.

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

TTree* treemc = (TTree*) filerun.Get("pndsim")

Definition at line 11 of file checkdedx_helixhit.C.

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

Definition at line 25 of file checkdedx_helixhit.C.

Referenced by QAmacro_stt_4().