reco_analys2.C File Reference

Go to the source code of this file.

Functions
gROOT	LoadMacro ("$VMCWORKDIR/gconfig/rootlogon.C")
gROOT	LoadMacro ("$VMCWORKDIR/gconfig/basiclibs.C")
	rootlogon ()
	basiclibs ()
t	SetBranchAddress ("EmcCluster",&cluster_array)
	for (Int_t j=0;j< t->GetEntriesFast();j++)
h1	SetTitle ("#pi^{0} Invariant mass reconstruction")
h1	SetStats (kTRUE)
h1	GetXaxis () -> SetTitle("Energy, GeV")
h1	Draw ()
f1	SetParameters (par1)
f1	SetLineColor (2)
h1	Fit ("f1")

Variables
TFile *	f = new TFile("full_emc.root")
TTree *	t =(TTree *) f->Get("pndsim")
TClonesArray *	cluster_array =new TClonesArray("PndEmcCluster")
double	cluster_energy
TH1F *	h1 = new TH1F("h1","Invariant mass",100,0.05,0.18)
TVector3	v1
TVector3	v2
TCanvas *	c1 = new TCanvas("c1", "", 100, 100, 800, 800)
TF1 *	f1 = new TF1("f1","gaus(0)",50.,150.)
Double_t	par1 [3] ={60,0.130,0.01}
double	mu1 =f1->GetParameter(1)
double	sigma1 =f1->GetParameter(2)

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

h1 Fit

(

"f1"

)

Referenced by PndCAGBTracker::CATrackFinder().

for

(

Int_t

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

)

Definition at line 22 of file reco_analys2.C.

References alpha, cos(), PndEmcCluster::energy(), i, sqrt(), threshold, and PndEmcCluster::where().

        {
                t->GetEntry(j);
                for (Int_t i=0; i<cluster_array->GetEntriesFast(); i++)
                {
                        PndEmcCluster *cluster1=(PndEmcCluster*)cluster_array->At(i);
                        if (cluster1->energy()<threshold) continue;
                        v1=cluster1->where();

                        for (Int_t k=i; k<cluster_array->GetEntriesFast(); k++)
                        {
                                PndEmcCluster *cluster2=(PndEmcCluster*)cluster_array->At(k);
                                if (cluster2->energy()<threshold) continue;
                                v2=cluster2->where();
                                double alpha=v1.Angle(v2);
                                double invMass=sqrt(2*cluster1->energy()*cluster2->energy()*(1-cos(alpha)));
                                h1->Fill(invMass);
                        }

                }
        }

h1 GetXaxis

(

)

-> SetTitle("Energy, GeV")

Definition at line 653 of file anal_hit_digi_cluster_fwendcap.C.

{Reconstructed}-#theta_{MC} (#circ)");

gROOT LoadMacro

(

"$VMCWORKDIR/gconfig/rootlogon.C"

)

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	(	"EmcCluster"	,
		&	cluster_array
	)

f1 SetLineColor

(

2

)

f1 SetParameters

(

par1

)

h1 SetStats

(

kTRUE

)

h1 SetTitle

(

"#pi^{0} Invariant mass reconstruction"

)

Variable Documentation

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

Definition at line 44 of file reco_analys2.C.

TClonesArray* cluster_array =new TClonesArray("PndEmcCluster")

Definition at line 14 of file reco_analys2.C.

double cluster_energy

Definition at line 17 of file reco_analys2.C.

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

Definition at line 10 of file reco_analys2.C.

TF1 * f1 = new TF1("f1","gaus(0)",50.,150.)

Definition at line 50 of file reco_analys2.C.

Referenced by analyse_etac1(), analyse_etac1_slc(), analyse_J(), analyse_J_slc(), analyse_phi(), analyse_phi2(), analyse_phi_base(), breitWigner(), compareDalitzPlot(), confgraph(), drawMvdHistos(), emc_correction_hist(), emc_correction_parametrization(), emc_correction_QA(), error_matrix_fit(), F32vec4::F32vec4(), PndRichRecoTask::FinishTask(), fit_resolution(), getpars(), KFParticle::GetPseudoProperDecayTime(), KFParticleSIMD::GetPseudoProperDecayTime(), HypHitAnasigBg_AR(), PndTrkCTGeometryCalculations::IsInsideArc(), PndSttTrackFinderReal::IsInsideArc(), PndMvdAdvancedPidAlgo::LandauGaus(), PndHypAdvancedPidAlgo::LandauGaus(), PndFsmMvdPid::LandauGaus(), PndFsmMvd2::LandauGaus(), LimitsAlign(), merge(), ModulesStat(), ModulesStatFast(), nastyFunction(), GFMaterialEffects::noiseBetheBloch(), pandaResponse(), pid_plot(), qadalitz(), qadalitz2(), PndFTSCAGBTracker::Refit_1(), SamplesSummaryAlign(), PndSdsFE::SaveInterpolatorList(), SenseAlign(), signalResponse(), sim_Gamma(), PndSttSingleStraw::STUrban(), PndFtsSingleStraw::STUrban(), SummaryAlign(), TOTALSummaryAlign(), and trickySummaryAlign().

TH1F* h1 = new TH1F("h1","Invariant mass",100,0.05,0.18)

Definition at line 18 of file reco_analys2.C.

double mu1 =f1->GetParameter(1)

Definition at line 56 of file reco_analys2.C.

Referenced by DecayTreeFitter::InternalParticle::initPar1().

Double_t par1[3] ={60,0.130,0.01}

Definition at line 51 of file reco_analys2.C.

Referenced by createdirc(), createdirc_prism(), createdircPix(), and PndDrcTimeDigiTask::FuncD1().

double sigma1 =f1->GetParameter(2)

Definition at line 57 of file reco_analys2.C.

Referenced by QAmacro_mvd_ana(), and QAmacro_stt_4().

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

Definition at line 12 of file reco_analys2.C.

TVector3 v1

Definition at line 20 of file reco_analys2.C.

TVector3 v2

Definition at line 20 of file reco_analys2.C.