#include "TCanvas.h"

Functions
gROOT	Macro ("$VMCWORKDIR/gconfig/rootlogon.C")

gROOT	LoadMacro ("Tools.C")

	rootlogon ()

timer	Start ()

tree2	SetBranchAddress ("LMDPoint",&point)

tree	SetBranchAddress ("LMDStripDigis",&digiStrip_array)

hisxy	SetTitle ("LMD MC Points, X vs Y view;X (cm);X (cm)")

hisrz	SetTitle ("LMD MC Point, rz view;z / cm;r/ cm")

hisde	SetTitle ("Lmd MC Energy Loss; X / keV ; Y / count")

hisStripTop	SetTitle ("FE numbers vs Channels counts; X / FE; Y / channel")

hisStripStrip	SetTitle ("Strips on Front/Back side; X / channels; Y / count")

hisStripChargeTop	SetLineColor (kBlue)

hisStripChargeBot	SetLineColor (kRed)

	if (verbose)

can1	Divide (a, b)

can1	cd (1)

	DrawNice2DHisto (hisxy)

can1	cd (2)

gPad	SetLogy ()

hisde	DrawCopy ()

can1	cd (3)

	DrawNice2DHisto (hisStripTop)

can1	cd (4)

	DrawNice2DHisto (hisStripBot)

can1	cd (5)

can1	cd (6)

hisStripChargeTop	DrawCopy ("sames")

can1	Update ()

	BetterStatBox (mypad)

timer	Stop ()

Variables
int	nEvents = 1e4

bool	verbose = false

TStopwatch	timer

std::string	inFile = "/private/huagen/simdata/Lmd_DPM_elastic_6.2_1.9mrad_5M_1.root"

std::string	digiFile = "/private/huagen/simdata/Lmd_Digi_DPM_elastic_6.2_1.9mrad_5M_1.root"

std::string	parfile = "/private/huagen/simdata/LmdParams.root"

TFile *	Digis = new TFile(digiFile.c_str(),"READ")

TTree *	tree =(TTree *) Digis->Get("pndsim")

TFile *	MCPoint = new TFile(inFile.c_str(),"READ")

TTree *	tree2 = (TTree*)MCPoint->Get("pndsim")

TClonesArray *	point =new TClonesArray("PndSdsMCPoint")

TClonesArray *	digiStrip_array =new TClonesArray("PndSdsDigiStrip")

TH2D *	hisxy = new TH2D("hisxy","",100,-10.,10.,100,-10.,10.)

TH2D *	hisrz = new TH2D("hisrz","",100,1070.,1150.,100,1070.,1150.)

TH1D *	hisde = new TH1D("hisde","LMD MC Points, Energyloss",1000,0.,1000)

TH2I *	hisStripTop = new TH2I("StripFront","Strip Front channel & front end numbers",20,0,40,130,0,130)

TH2I *	hisStripBot = new TH2I("StripBack","Strip Back channel & front end numbers",20,0,40,130,0,130)

TH1I *	hisStripStrip = new TH1I("hisstripstrip","Strip numbers",2330,-10,4650)

TH1D *	hisStripCharge = new TH1D("StripChargeTotal","Strip Charge content",1000,0.,2.5e5)

TH1D *	hisStripChargeTop = new TH1D("StripChargeFront","Strip Charge content",1000,0.,2.5e5)

TH1D *	hisStripChargeBot = new TH1D("StripChargeBack","Strip Charge content",1000,0.,2.5e5)

TVector3	vecmc

TVector3	mommc

Double_t	tmpx

Double_t	tmpy

Double_t	tmpz

TVector2	locals

TVector2	localmc

TVector2	localdiff

int	col

int	row

int	fe

double	x

double	y

double	X =0

double	Y =0

double	Z =0

double	de =0

TFile *	parDB = new TFile(parfile.c_str())

PndSdsStripDigiPar *	par = (PndSdsStripDigiPar*)gROOT->FindObject("LmdStripDigiParTrap")

int	nrFeChannels = 128

int	nrStrips = 1280

Int_t	a = 2

Int_t	b = 3

TCanvas *	can1 = new TCanvas("Lmd MC point & Digis Plot","MCHit view in LMD",0,0,2400,2400)

TPad *	mypad =0

Double_t	rtime = timer.RealTime()

Double_t	ctime = timer.CpuTime()

Function Documentation

BetterStatBox ( mypad )

Referenced by anaclust().

can1 cd ( 1 )

can1 cd ( 2 )

can1 cd ( 3 )

can1 cd ( 4 )

can1 cd ( 5 )

can1 cd ( 6 )

can1 Divide	(	a	,
		b
	)

hisde DrawCopy ( )

hisStripChargeBot DrawCopy ( "sames" )

DrawNice2DHisto ( hisxy )

Referenced by anaclust(), anaGemPointrate(), anaGemSmearing(), and plothistos().

DrawNice2DHisto ( hisStripTop )

DrawNice2DHisto ( hisStripBot )

if ( verbose )

Definition at line 77 of file anaLmdDigi.C.

References col, de, PndSdsDigiStrip::GetChannel(), PndSdsDigi::GetCharge(), PndSdsDigi::GetFE(), PndSdsMCPoint::GetPosition(), i, Y, and Z.

                                                      :"<< nrFeChannels<<","<<nrStrips<<endl;
 
   for (Int_t j=0; j<nEvents&&j<tree->GetEntriesFast(); j++)
   {
         if(verbose) std::cout<< "the number of "<< j <<"event is been processed"<<std::endl;
      tree->GetEntry(j);
      tree2->GetEntry(j);
     if(verbose)         cout<<"Event No "<<j<<endl;
 
     // ----- MC Points -----
     for (Int_t i=0; i<point->GetEntriesFast(); i++)
     {
         if(verbose)    cout<<"Point No "<<i<<endl;
       PndSdsMCPoint *Mypoint=(PndSdsMCPoint*)point->At(i);
           X = Mypoint->GetPosition().X();
           Y = Mypoint->GetPosition().Y();
           Z = Mypoint->GetPosition().Z();
           hisxy->Fill(X,Y);
       de = Mypoint->GetEnergyLoss()*1000000;
      // de = Mypoint->GetEnergyLoss();
       hisde->Fill(de);
       if(verbose) cout<<"the energy loss is "<<de<<endl;
     }
 
     // ----- STRIP DIGIS -----
     for (Int_t i=0; i<digiStrip_array->GetEntriesFast(); i++)
     {
       PndSdsDigiStrip *stripdigi = digiStrip_array->At(i);
       fe = stripdigi->GetFE();
       col = stripdigi->GetChannel();
   //    std::cout<<"the fe and channel are:"<<fe<<","<<col<<std::endl;
 
       int strip = fe * nrFeChannels + col;
       hisStripStrip->Fill(strip);
       hisStripCharge->Fill(stripdigi->GetCharge());
   //    std::cout<<"the strip charge is :"<<stripdigi->GetCharge()<<std::endl;
 
       if (strip <= nrStrips)
       {
         hisStripTop->Fill(fe,col);
         hisStripChargeTop->Fill(stripdigi->GetCharge());
       } else {
         hisStripBot->Fill(fe,col);
         hisStripChargeBot->Fill(stripdigi->GetCharge());
       }
     }
 
   }// end for j (events)

gROOT LoadMacro ( "Tools.C" )

gROOT Macro ( "$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);
 }