LIAnalysis Class Reference

#include <LIAnalysis.h>

Inheritance diagram for LIAnalysis:

List of all members.

Public Member Functions
	LIAnalysis (Int_t analyseChainFlag)
	~LIAnalysis ()
void	AdcVsChannel (Int_t plane, Int_t crate)
void	AdcDistribution ()
void	AdcVsLed ()
void	AdcVsPin ()
void	AdcVsPixel (Int_t plane, Int_t AshtrayMin=-1, Int_t AshtrayMax=-1, Int_t led=-1)
void	AshtrayVsPlane (Int_t nearPb, Int_t farPb)
void	CalcAdcAv ()
void	ChannelAdc (Int_t chan, Int_t chan2)
void	CheckData (Int_t rt=-1)
void	ClearFibres ()
void	GainCurves (Int_t pulserBoxToPlot=0, Int_t maxLedToPlot=6)
void	GainCurves_Led ()
void	GainVsTime ()
void	DebugPins ()
void	HighGainSearch (Int_t maxAdcThreshold, Float_t gainThreshold)
void	IndividualChannels ()
void	LedTuning (Int_t pNum, Int_t pWidth, Int_t idealAdc, const Int_t numGcPoints, Double_t firstGcPoint, Double_t lastGcPoint, Double_t dADCdPHAtSat, Int_t numCalibPoints=-1, Int_t numSecondFile=0)
void	LISnarlProblem ()
void	MiswiringSearch (Int_t adcThreshold)
void	NumPulses ()
void	PinDiode ()
void	PinDiodeChips ()
void	PinMap ()
void	PixelVsPlane (Int_t nearPb, Int_t farPb)
void	PmtGain (Int_t plane, Int_t cr, Int_t printGain)
void	PrintElec ()
void	PrintPmt ()
void	Reflectors (LILookup::ELINearOrFarSide nearOrFarSide)
void	ReflectorsGc (LILookup::ELINearOrFarSide nearOrFarSide)
void	SearchForBadLeds (Int_t minNumHits=200)
void	StripVsPlane (Int_t nearPb=0, Int_t farPb=0, Int_t ledCut=-1)
void	StripVsPlaneWholeDet ()
void	Template ()
void	Test ()
void	TriggerPmt ()
void	WriteGainCurveTextFile ()
void	WriteGainsTextFile (Int_t task=1)
void	WriteOldGainsTextFile ()
void	NoisyChips (Float_t)
Private Member Functions
void	AnalyseChain ()
void	AppendVect (std::vector< Double_t > &trunk, std::vector< Double_t > appendix)
std::string	GetElecString ()
void	InitialiseLoopVariables ()
void	MakeChain ()
void	PrintBigMessage ()
void	PrintBlockInfo (std::string preString)
void	PrintGainTableRow (Int_t pl, Int_t st, Int_t end, std::string sPrefix, std::string sAdcN, std::string sGainN, std::string sGainErrN, std::string sNumN, std::string sAdcF, std::string sGainF, std::string sGainErrF, std::string sNumF, std::string sLogLevel="Debug")
void	SetChainBranches ()
void	SetLoopVariables (Int_t entry, Int_t printOnNewLed=0, Bool_t doInfoPrint=true)
TGraph *	TGraphVect (std::vector< Double_t > &vX, std::vector< Double_t > &vY)
Private Attributes
Int_t	ashtray
Int_t	calibPoint
Int_t	calibType
Int_t	chAdd
TChain *	chain
Int_t	channel
Int_t	chip
Int_t	correlatedHit
Int_t	crate
Int_t	detectorType
Int_t	eastWest
Int_t	elecType
Int_t	farLed
Int_t	farPulserBox
Int_t	firstRunNumber
Int_t	geoAdd
char	histname [80]
Int_t	highRunNumber
Int_t	inRack
Int_t	lastLed
Int_t	lastCalibPoint
Int_t	lastPulserBox
Int_t	lastRunNumber
Int_t	led
Int_t	liEvent
Int_t	liRunNum
LILookup	lookup
Int_t	lowRunNumber
Int_t	masterCh
Int_t	maxLedNum
Int_t	maxPbNum
Int_t	maxCalibPoint
Float_t	mean
Int_t	minderCh
Int_t	nearLed
Int_t	nearPulserBox
Int_t	numCalibPoints
Int_t	numEntries
Int_t	numericMuxBox
Int_t	numEvents
Int_t	numLeds
Int_t	numLiEvents
std::map< Int_t, Int_t >	numLiEventsL
std::map< Int_t, Int_t >	numLiEventsP
Int_t	numLiRuns
Int_t	period
Int_t	pixel
Int_t	pinGain
Int_t	pinInBox
Int_t	plane
Int_t	previousRunNumber
Int_t	pulseHeight
Int_t	pulserBox
Int_t	pulses
Int_t	pulseWidth
Int_t	rackBay
Int_t	rackLevel
Int_t	readoutType
Float_t	rms
Int_t	run
Int_t	runNumber
Int_t	runNumberSub
Int_t	runType
std::string	s
std::string	fS
Int_t	strip
Int_t	stripEnd
Int_t	summaryCounter
Int_t	timestamp
Int_t	timestampNanoSec
Int_t	varc
Int_t	vfb
Int_t	vmm

---

Constructor & Destructor Documentation

LIAnalysis::LIAnalysis (  Int_t  analyseChainFlag  )

Definition at line 65 of file LIAnalysis.cxx. References AnalyseChain(), ashtray, calibPoint, calibType, chAdd, chain, channel, chip, correlatedHit, crate, detectorType, eastWest, elecType, farLed, farPulserBox, firstRunNumber, fS, geoAdd, highRunNumber, histname, inRack, lastCalibPoint, lastLed, lastPulserBox, lastRunNumber, led, liEvent, liRunNum, lowRunNumber, MakeChain(), masterCh, maxCalibPoint, maxLedNum, maxPbNum, mean, minderCh, MSG, nearLed, nearPulserBox, numCalibPoints, numEntries, numericMuxBox, numEvents, numLeds, numLiEvents, numLiEventsL, numLiEventsP, numLiRuns, period, pinGain, pinInBox, pixel, plane, previousRunNumber, pulseHeight, pulserBox, pulses, pulseWidth, rackBay, rackLevel, readoutType, rms, run, runNumber, runNumberSub, runType, s, SetChainBranches(), LILookup::SetDetector(), strip, stripEnd, summaryCounter, timestamp, timestampNanoSec, varc, vfb, and vmm. { MSG("LIAnalysis", Msg::kDebug) <<"Running LIAnalysis constructor..."<<endl; //use the pretty palette gStyle->SetPalette(1); //include the under and overflow counts gStyle->SetOptStat(1111111); gStyle->SetOptFit(1111); //data members ashtray=-1; calibPoint=-1; calibType=-1; chAdd=-1; chain=0; channel=-1; chip=-1; correlatedHit=-1; crate=-1; detectorType=-1; eastWest=-1;//0 is east, 1 is west elecType=-1; farLed=-1; farPulserBox=-1; firstRunNumber=-1; geoAdd=-1; highRunNumber=-1; for (Int_t i=0;i<80;i++){ histname[i]='?'; } inRack=-1;//determines mux box in rack 0-7 lastLed=-1; lastCalibPoint=-1; lastPulserBox=-1; lastRunNumber=-1; led=-1; liEvent=-1; liRunNum=-1; lowRunNumber=-1; masterCh=-1; maxLedNum=-1; maxPbNum=-1; maxCalibPoint=-1; mean=-1.; minderCh=-1; nearLed=-1; nearPulserBox=-1; numCalibPoints=0; numEntries=-1; numericMuxBox=-1; numEvents=-1; numLeds=-1; numLiEvents=-1; for (Int_t i=0;i<NUMLEDS;i++){ numLiEventsL[i]=0; } for (Int_t i=0;i<NUMPULSERBOXES;i++){ numLiEventsP[i]=0; } numLiRuns=-1; period=-1; pixel=-1; pinGain=-1;//0 is low, 1 is high gain pinInBox=-1;//0 is low gain (VA=0), 1 is high gain (VA=1) plane=-1; previousRunNumber=-1; pulseHeight=-1; pulserBox=-1; pulses=-1; pulseWidth=-1; rackBay=-1;//determines which rack 0-15 rackLevel=-1;//0 is top, 1 is bottom readoutType=-1; rms=-1.; run=-1; runNumber=-1; runNumberSub=-1; runType=-1; s=""; fS=""; strip=-1; stripEnd=-1; summaryCounter=-1; timestamp=-1; timestampNanoSec=-1; varc=-1; vfb=-1; vmm=-1; //set up the chain and perform some checks this->MakeChain(); this->SetChainBranches(); numEvents=static_cast<Int_t>(chain->GetEntries()); MSG("LIAnalysis",Msg::kInfo) <<"Number of events in chain="<<numEvents<<endl; if (analyseChainFlag==1){ this->AnalyseChain(); } //get the last event chain->GetEvent(numEvents-1); lastRunNumber=runNumber; //get the first event chain->GetEvent(0); firstRunNumber=runNumber;//never need to touch this again //set which detector to use this->SetDetector(detectorType); //set the high run number to be the highest of the first and last //it will get changed if there is a higher run number //between the first and the last highRunNumber=firstRunNumber; if (lastRunNumber>firstRunNumber) highRunNumber=lastRunNumber; //similarly for low run number lowRunNumber=firstRunNumber; if (lastRunNumber<firstRunNumber) lowRunNumber=lastRunNumber; if (firstRunNumber!=lastRunNumber){ MSG("LIAnalysis",Msg::kInfo) <<"First run number = "<<firstRunNumber <<", last = "<<lastRunNumber<<endl; } MSG("LIAnalysis", Msg::kInfo) <<"Finished LIAnalysis constructor"<<endl; }

LIAnalysis::~LIAnalysis (   )

Definition at line 198 of file LIAnalysis.cxx. References MSG. { MSG("LIAnalysis", Msg::kDebug) <<"Running LIAnalysis destructor..."<<endl; MSG("LIAnalysis", Msg::kDebug) <<"Finished LIAnalysis destructor"<<endl; }

---

Member Function Documentation

void LIAnalysis::AdcDistribution (   )

Definition at line 6350 of file LIAnalysis.cxx. References calibPoint, correlatedHit, crate, detectorType, eastWest, farPulserBox, Form(), histname, InitialiseLoopVariables(), led, lookup, mean, MSG, LILookup::NearOrFar(), nearPulserBox, numEntries, period, plane, pulseHeight, pulserBox, pulseWidth, readoutType, rms, runNumber, s, and SetLoopVariables(). { MSG("LIAnalysis",Msg::kInfo) <<endl<<" ** Running the AdcDistribution method... ** "<<endl; //include the under and overflow counts gStyle->SetOptStat(1111111); TH1F *hMean= new TH1F("hMean","Mean",150,-5,15005); hMean->GetXaxis()->SetTitle("Mean"); hMean->GetXaxis()->CenterTitle(); hMean->GetYaxis()->SetTitle("Number"); hMean->GetYaxis()->CenterTitle(); hMean->SetFillColor(0); hMean->SetBit(TH1::kCanRebin); TH1F *hMeanCorrelated= new TH1F("hMeanCorrelated","MeanCorrelated", 150,-5,15005); hMeanCorrelated->GetXaxis()->SetTitle("Mean"); hMeanCorrelated->GetXaxis()->CenterTitle(); hMeanCorrelated->GetYaxis()->SetTitle("Number"); hMeanCorrelated->GetYaxis()->CenterTitle(); hMeanCorrelated->SetFillColor(0); hMeanCorrelated->SetLineColor(4); hMeanCorrelated->SetBit(TH1::kCanRebin); TH1F *hMeanEast= new TH1F("hMeanEast","MeanEast", 150,-5,15005); hMeanEast->GetXaxis()->SetTitle("Mean"); hMeanEast->GetXaxis()->CenterTitle(); hMeanEast->GetYaxis()->SetTitle("Number"); hMeanEast->GetYaxis()->CenterTitle(); hMeanEast->SetLineColor(4); hMeanEast->SetFillStyle(3015); hMeanEast->SetFillColor(4); hMeanEast->SetLineWidth(6); hMeanEast->SetBit(TH1::kCanRebin); TH1F *hMeanWest= new TH1F("hMeanWest","MeanWest", 150,-5,15005); hMeanWest->GetXaxis()->SetTitle("Mean"); hMeanWest->GetXaxis()->CenterTitle(); hMeanWest->GetYaxis()->SetTitle("Number"); hMeanWest->GetYaxis()->CenterTitle(); hMeanWest->SetLineColor(2); hMeanWest->SetFillStyle(3015); hMeanWest->SetFillColor(2); hMeanWest->SetLineWidth(6); hMeanWest->SetBit(TH1::kCanRebin); TH1F *hMeanNearPb= new TH1F("hMeanNearPb","MeanNearPb", 150,-5,15005); hMeanNearPb->GetXaxis()->SetTitle("Mean"); hMeanNearPb->GetXaxis()->CenterTitle(); hMeanNearPb->GetYaxis()->SetTitle("Number"); hMeanNearPb->GetYaxis()->CenterTitle(); hMeanNearPb->SetFillColor(0); hMeanNearPb->SetLineColor(2); hMeanNearPb->SetBit(TH1::kCanRebin); TH1F *hMeanFarPb= new TH1F("hMeanFarPb","MeanFarPb", 150,-5,15005); hMeanFarPb->GetXaxis()->SetTitle("Mean"); hMeanFarPb->GetXaxis()->CenterTitle(); hMeanFarPb->GetYaxis()->SetTitle("Number"); hMeanFarPb->GetYaxis()->CenterTitle(); hMeanFarPb->SetFillColor(0); hMeanFarPb->SetLineColor(3); hMeanFarPb->SetBit(TH1::kCanRebin); //histos for different calib points / pulse heights const Int_t const_numCalibPoints=10; Int_t pulseHeights[const_numCalibPoints]; TH1F **hAdcPh=0; hAdcPh= new TH1F*[const_numCalibPoints]; for (Int_t i=0;i<const_numCalibPoints;i++){ //initialise heights pulseHeights[i]=-1; //initialse histos sprintf(histname,"Adc Distribution, calibpoint %d",i+1); hAdcPh[i]=new TH1F(histname,histname,200,0,15000); hAdcPh[i]->GetXaxis()->SetTitle("ADC"); hAdcPh[i]->GetXaxis()->CenterTitle(); hAdcPh[i]->GetYaxis()->SetTitle("Number of Entries"); hAdcPh[i]->GetYaxis()->CenterTitle(); hAdcPh[i]->SetFillColor(0); hAdcPh[i]->Fill(1); //hAdcPh[i]->SetBit(TH1::kCanRebin); } TH1F **hAdcPhNearPb=0; hAdcPhNearPb= new TH1F*[const_numCalibPoints]; for (Int_t i=0;i<const_numCalibPoints;i++){ sprintf(histname,"AdcPhNearPb Distribution, calibpoint %d",i+1); hAdcPhNearPb[i]=new TH1F(histname,histname,200,0,15000); hAdcPhNearPb[i]->GetXaxis()->SetTitle("ADC"); hAdcPhNearPb[i]->GetXaxis()->CenterTitle(); hAdcPhNearPb[i]->GetYaxis()->SetTitle("Number of Entries"); hAdcPhNearPb[i]->GetYaxis()->CenterTitle(); hAdcPhNearPb[i]->SetFillColor(0); hAdcPhNearPb[i]->SetLineColor(2); hAdcPhNearPb[i]->Fill(1); //hAdcPhNearPb[i]->SetBit(TH1::kCanRebin); } TH1F **hAdcPhFarPb=0; hAdcPhFarPb= new TH1F*[const_numCalibPoints]; for (Int_t i=0;i<const_numCalibPoints;i++){ sprintf(histname,"AdcPh FarPb Distribution, calibpoint %d",i+1); hAdcPhFarPb[i]=new TH1F(histname,histname,200,0,15000); hAdcPhFarPb[i]->GetXaxis()->SetTitle("ADC"); hAdcPhFarPb[i]->GetXaxis()->CenterTitle(); hAdcPhFarPb[i]->GetYaxis()->SetTitle("Number of Entries"); hAdcPhFarPb[i]->GetYaxis()->CenterTitle(); hAdcPhFarPb[i]->SetFillColor(0); hAdcPhFarPb[i]->SetLineColor(3); hAdcPhFarPb[i]->Fill(1); //hAdcPhFarPb[i]->SetBit(TH1::kCanRebin); } //histos for different pulser boxes TH1F **hAdcPb=0; hAdcPb= new TH1F*[NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ //initialse histos sprintf(histname,"Adc Distribution, pulser box %d",i); hAdcPb[i]=new TH1F(histname,histname,200,0,15000); hAdcPb[i]->GetXaxis()->SetTitle("ADC"); hAdcPb[i]->GetXaxis()->CenterTitle(); hAdcPb[i]->GetYaxis()->SetTitle("Number of Entries"); hAdcPb[i]->GetYaxis()->CenterTitle(); hAdcPb[i]->SetFillColor(0); hAdcPb[i]->Fill(1); //hAdcPb[i]->SetBit(TH1::kCanRebin); } TH1F **hAdcPbNearPb=0; hAdcPbNearPb= new TH1F*[NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ sprintf(histname,"AdcPbNearPb Distribution, pulser box %d",i); hAdcPbNearPb[i]=new TH1F(histname,histname,200,0,15000); hAdcPbNearPb[i]->GetXaxis()->SetTitle("ADC"); hAdcPbNearPb[i]->GetXaxis()->CenterTitle(); hAdcPbNearPb[i]->GetYaxis()->SetTitle("Number of Entries"); hAdcPbNearPb[i]->GetYaxis()->CenterTitle(); hAdcPbNearPb[i]->SetFillColor(0); hAdcPbNearPb[i]->SetLineColor(2); hAdcPbNearPb[i]->Fill(1); //hAdcPbNearPb[i]->SetBit(TH1::kCanRebin); } TH1F **hAdcPbFarPb=0; hAdcPbFarPb= new TH1F*[NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ sprintf(histname,"AdcPbFarPb Distribution, pulser box %d",i); hAdcPbFarPb[i]=new TH1F(histname,histname,200,0,15000); hAdcPbFarPb[i]->GetXaxis()->SetTitle("ADC"); hAdcPbFarPb[i]->GetXaxis()->CenterTitle(); hAdcPbFarPb[i]->GetYaxis()->SetTitle("Number of Entries"); hAdcPbFarPb[i]->GetYaxis()->CenterTitle(); hAdcPbFarPb[i]->SetFillColor(0); hAdcPbFarPb[i]->SetLineColor(3); hAdcPbFarPb[i]->Fill(1); //hAdcPbFarPb[i]->SetBit(TH1::kCanRebin); } //histos for different pulse widths const Int_t const_numPulseWidths=3; TH1F **hAdcPw=0; hAdcPw= new TH1F*[const_numPulseWidths]; for (Int_t i=0;i<const_numPulseWidths;i++){ //initialse histos sprintf(histname,"Adc Distribution, pulse width %d",i); hAdcPw[i]=new TH1F(histname,histname,200,0,15000); hAdcPw[i]->GetXaxis()->SetTitle("ADC"); hAdcPw[i]->GetXaxis()->CenterTitle(); hAdcPw[i]->GetYaxis()->SetTitle("Number of Entries"); hAdcPw[i]->GetYaxis()->CenterTitle(); hAdcPw[i]->SetFillColor(0); hAdcPw[i]->Fill(1); //hAdcPw[i]->SetBit(TH1::kCanRebin); } TH1F **hAdcPwNearPb=0; hAdcPwNearPb= new TH1F*[const_numPulseWidths]; for (Int_t i=0;i<const_numPulseWidths;i++){ sprintf(histname,"AdcPwNearPb Distribution, pulser box %d",i); hAdcPwNearPb[i]=new TH1F(histname,histname,200,0,15000); hAdcPwNearPb[i]->GetXaxis()->SetTitle("ADC"); hAdcPwNearPb[i]->GetXaxis()->CenterTitle(); hAdcPwNearPb[i]->GetYaxis()->SetTitle("Number of Entries"); hAdcPwNearPb[i]->GetYaxis()->CenterTitle(); hAdcPwNearPb[i]->SetFillColor(0); hAdcPwNearPb[i]->SetLineColor(2); hAdcPwNearPb[i]->Fill(1); //hAdcPwNearPb[i]->SetBit(TH1::kCanRebin); } TH1F **hAdcPwFarPb=0; hAdcPwFarPb= new TH1F*[const_numPulseWidths]; for (Int_t i=0;i<const_numPulseWidths;i++){ sprintf(histname,"AdcPwFarPb Distribution, pulser box %d",i); hAdcPwFarPb[i]=new TH1F(histname,histname,200,0,15000); hAdcPwFarPb[i]->GetXaxis()->SetTitle("ADC"); hAdcPwFarPb[i]->GetXaxis()->CenterTitle(); hAdcPwFarPb[i]->GetYaxis()->SetTitle("Number of Entries"); hAdcPwFarPb[i]->GetYaxis()->CenterTitle(); hAdcPwFarPb[i]->SetFillColor(0); hAdcPwFarPb[i]->SetLineColor(3); hAdcPwFarPb[i]->Fill(1); //hAdcPwFarPb[i]->SetBit(TH1::kCanRebin); } this->InitialiseLoopVariables(); for(Int_t entry=0;entry<numEvents;entry++){ this->SetLoopVariables(entry,0); //avoid divide by zero errors and skip irrelevant data if (rms==0. || mean==0. || numEntries==0) continue; //only look at scintillator strips if (readoutType!=ReadoutType::kScintStrip) continue; hMean->Fill(mean); if (correlatedHit==1){ if (detectorType==Detector::kCalDet){ if (plane>0 && plane<=LASTPLANE && lookup.NearOrFar(crate,pulserBox,nearPulserBox, farPulserBox,led,detectorType, plane,runNumber)== LILookup::kNearSide){ if (eastWest==1) hMeanEast->Fill(mean); else if (eastWest==2) hMeanWest->Fill(mean); } } //fill histogram for appropriate pulse height if (calibPoint>0 && calibPoint<=const_numCalibPoints){ pulseHeights[calibPoint-1]=pulseHeight; hAdcPh[calibPoint-1]->Fill(mean); } //fill histo for appropriate pulser box if (pulserBox>=FIRSTPULSERBOX && pulserBox<=LASTPULSERBOX){ hAdcPb[pulserBox]->Fill(mean); } //fill histo for appropriate pulse width if (pulseWidth>=1 && pulseWidth<=const_numPulseWidths){ hAdcPw[pulseWidth-1]->Fill(mean); } //fill general histo hMeanCorrelated->Fill(mean); //fill histos with when pulser box is near pulser box if (nearPulserBox==pulserBox) { //fill histogram for appropriate pulse height if (calibPoint>0 && calibPoint<=const_numCalibPoints){ hAdcPhNearPb[calibPoint-1]->Fill(mean); } //fill histo for appropriate pulser box if (pulserBox>=FIRSTPULSERBOX && pulserBox<=LASTPULSERBOX){ hAdcPbNearPb[pulserBox]->Fill(mean); } //fill histo for appropriate pulse width if (pulseWidth>=1 && pulseWidth<=const_numPulseWidths){ hAdcPwNearPb[pulseWidth-1]->Fill(mean); } //fill general histo hMeanNearPb->Fill(mean); } //fill histos with when pulser box is far pulser box else if (farPulserBox==pulserBox) { //fill histogram for appropriate pulse height if (calibPoint>0 && calibPoint<=const_numCalibPoints){ hAdcPhFarPb[calibPoint-1]->Fill(mean); } //fill histo for appropriate pulser box if (pulserBox>=FIRSTPULSERBOX && pulserBox<=LASTPULSERBOX){ hAdcPbFarPb[pulserBox]->Fill(mean); } //fill histo for appropriate pulse width if (pulseWidth>=1 && pulseWidth<=const_numPulseWidths){ hAdcPwFarPb[pulseWidth-1]->Fill(mean); } //fill general histo hMeanFarPb->Fill(mean); } else { MSG("LIAnalysis",Msg::kWarning) <<" ** Wrong Pb **"<<endl; } } }//end of for MSG("LIAnalysis",Msg::kInfo)<<"Finished main loop"<<endl; string sRunNumber=Form("%d",runNumber); TCanvas *cMean=new TCanvas("cMean","cMean",0,0,1000,600); cMean->SetFillColor(0); cMean->cd(); cMean->SetLogy(); hMean->Draw(); hMeanNearPb->Draw("same"); hMeanFarPb->Draw("same"); //hMeanCorrelated->Draw("same"); TCanvas *cMeanEW=new TCanvas("cMeanEW","cMeanEW",0,0,1000,600); cMeanEW->SetFillColor(0); cMeanEW->cd(); hMeanWest->Draw(); hMeanEast->Draw("sames"); //sorted by calibration point section //get the maximum in order to scale all plots Int_t maxNumEnt=0; for (Int_t i=0;i<const_numCalibPoints;i++){ if (hAdcPh[i]->GetMaximum()>maxNumEnt){ maxNumEnt=static_cast<Int_t>(hAdcPh[i]->GetMaximum()); MSG("LIAnalysis",Msg::kInfo) <<"Calculating max in calib point section, current highest=" <<maxNumEnt<<endl; } } //create the canvas and draw TCanvas *cAdcPh=new TCanvas("cAdcPh","cAdcPh",0,0,1000,600); cAdcPh->SetFillColor(0); cAdcPh->cd(); cAdcPh->SetLogy(); MSG("LIAnalysis",Msg::kInfo) <<endl<<"Number of calibration points = "<<const_numCalibPoints <<endl; for (Int_t i=0;i<const_numCalibPoints;i++){ cAdcPh->Clear(); hAdcPh[i]->SetMaximum(maxNumEnt); string sPulseWidth=Form("%d",pulseWidth); string sPulseHeight=Form("%d",pulseHeights[i]); string sPulseFreq=Form("%d",static_cast<Int_t> (ceil(1.0/(period*1.0e-5)))); MSG("LIAnalysis",Msg::kInfo) <<"LI parameters: " <<", PH="<<sPulseHeight <<", PW="<<sPulseWidth <<", PF="<<sPulseFreq <<", period="<<period <<endl; s="ADC Values, Near&Far side, PH="+ sPulseHeight+", PW="+sPulseWidth+ ", PF="+sPulseFreq+" Hz"; hAdcPh[i]->SetTitle(s.c_str()); hAdcPh[i]->Draw(); hAdcPhNearPb[i]->Draw("same"); hAdcPhFarPb[i]->Draw("same"); if (i==0){ s=sRunNumber+"AdcPhHisto.ps("; cAdcPh->Print(s.c_str()); } else if(i==const_numCalibPoints-1){ s=sRunNumber+"AdcPhHisto.ps)"; cAdcPh->Print(s.c_str()); } else{ s=sRunNumber+"AdcPhHisto.ps"; cAdcPh->Print(s.c_str()); } } cAdcPh->Clear(); hAdcPh[0]->Draw(); hAdcPhNearPb[0]->Draw("same"); hAdcPhFarPb[0]->Draw("same"); //sorted by pulser box section //get the maximum in order to scale all plots maxNumEnt=0; for (Int_t i=0;i<NUMPULSERBOXES;i++){ if (hAdcPb[i]->GetMaximum()>maxNumEnt){ maxNumEnt=static_cast<Int_t>(hAdcPb[i]->GetMaximum()); MSG("LIAnalysis",Msg::kInfo) <<"Calculating max in pulser box section, current highest=" <<maxNumEnt<<endl; } } //create the canvas and draw TCanvas *cAdcPb=new TCanvas("cAdcPb","cAdcPb",0,0,1000,600); cAdcPb->SetFillColor(0); cAdcPb->cd(); cAdcPb->SetLogy(); MSG("LIAnalysis",Msg::kInfo) <<endl<<"Number of pulser boxes = "<<NUMPULSERBOXES<<endl; for (Int_t i=0;i<NUMPULSERBOXES;i++){ cAdcPb->Clear(); hAdcPb[i]->SetMaximum(maxNumEnt); string sPulserBox=Form("%d",i); string sPulseWidth=Form("%d",pulseWidth); string sPulseHeight=Form("%d",pulseHeight); string sPulseFreq=Form("%d",static_cast<Int_t> (ceil(1.0/(period*1.0e-5)))); MSG("LIAnalysis",Msg::kInfo) <<"LI parameters: " <<", PB="<<sPulserBox <<", PH="<<sPulseHeight <<", PW="<<sPulseWidth <<", PF="<<sPulseFreq <<", period="<<period <<endl; s="ADC Values, Near&Far side, PB="+sPulserBox+ ", PH="+sPulseHeight+", PW="+sPulseWidth+ ", PF="+sPulseFreq+" Hz"; hAdcPb[i]->SetTitle(s.c_str()); hAdcPb[i]->Draw(); hAdcPbNearPb[i]->Draw("same"); hAdcPbFarPb[i]->Draw("same"); if (i==0){ s=sRunNumber+"AdcPbHisto.ps("; cAdcPb->Print(s.c_str()); } else if(i==NUMPULSERBOXES-1){ s=sRunNumber+"AdcPbHisto.ps)"; cAdcPb->Print(s.c_str()); } else{ s=sRunNumber+"AdcPbHisto.ps"; cAdcPb->Print(s.c_str()); } } cAdcPb->Clear(); hAdcPb[0]->Draw(); hAdcPbNearPb[0]->Draw("same"); hAdcPbFarPb[0]->Draw("same"); //sorted by pulse width //get the maximum in order to scale all plots maxNumEnt=0; for (Int_t i=0;i<const_numPulseWidths;i++){ if (hAdcPw[i]->GetMaximum()>maxNumEnt){ maxNumEnt=static_cast<Int_t>(hAdcPw[i]->GetMaximum()); MSG("LIAnalysis",Msg::kInfo) <<"Calculating max in pulse width section, current highest=" <<maxNumEnt<<endl; } } //create the canvas and draw TCanvas *cAdcPw=new TCanvas("cAdcPw","cAdcPw",0,0,1000,600); cAdcPw->SetFillColor(0); cAdcPw->cd(); cAdcPw->SetLogy(); MSG("LIAnalysis",Msg::kInfo) <<endl<<"Number of pulse widths = "<<const_numPulseWidths<<endl; for (Int_t i=0;i<const_numPulseWidths;i++){ cAdcPw->Clear(); hAdcPw[i]->SetMaximum(maxNumEnt); string sPulseWidth=Form("%d",i+1); string sPulseHeight=Form("%d",pulseHeight); string sPulseFreq=Form("%d",static_cast<Int_t> (ceil(1.0/(period*1.0e-5)))); MSG("LIAnalysis",Msg::kInfo) <<"LI parameters: " <<", PH="<<sPulseHeight <<", PW="<<sPulseWidth <<", PF="<<sPulseFreq <<", period="<<period <<endl; s="ADC Values, Near&Far side, PH="+sPulseHeight+ ", PW="+sPulseWidth+ ", PF="+sPulseFreq+" Hz"; hAdcPw[i]->SetTitle(s.c_str()); hAdcPw[i]->Draw(); hAdcPwNearPb[i]->Draw("same"); hAdcPwFarPb[i]->Draw("same"); if (i==0){ s=sRunNumber+"AdcPwHisto.ps("; cAdcPw->Print(s.c_str()); } else if(i==const_numPulseWidths-1){ s=sRunNumber+"AdcPwHisto.ps)"; cAdcPw->Print(s.c_str()); } else{ s=sRunNumber+"AdcPwHisto.ps"; cAdcPw->Print(s.c_str()); } } cAdcPw->Clear(); hAdcPw[0]->Draw(); hAdcPwNearPb[0]->Draw("same"); hAdcPwFarPb[0]->Draw("same"); MSG("LIAnalysis",Msg::kInfo) <<endl<<" ** Finished the AdcDistribution method ** "<<endl; }

void LIAnalysis::AdcVsChannel (  Int_t  plane, Int_t  crate )

Definition at line 5708 of file LIAnalysis.cxx. References channel, chip, correlatedHit, crate, histname, InitialiseLoopVariables(), mean, MSG, numEntries, pixel, plane, pulserBox, readoutType, rms, SetLoopVariables(), and strip. { //get rid of the stats info gStyle->SetOptStat(0); MSG("LIAnalysis",Msg::kInfo) <<endl<<" ** Running the AdcVsChannel method... **"<<endl; TH2F **hAdcVsChannel=0; hAdcVsChannel= new TH2F*[NUMCHIPS]; for (Int_t i=0;i<NUMCHIPS;i++){ sprintf(histname,"Adc Vs Channel (VA chip %i)",i); hAdcVsChannel[i]=new TH2F (histname,histname,23,0,23,15000,0,15000); hAdcVsChannel[i]->GetXaxis()->SetTitle("Channel"); hAdcVsChannel[i]->GetXaxis()->CenterTitle(); hAdcVsChannel[i]->GetYaxis()->SetTitle("Adc"); hAdcVsChannel[i]->GetYaxis()->CenterTitle(); hAdcVsChannel[i]->SetFillColor(0); hAdcVsChannel[i]->SetBit(TH1::kCanRebin); } this->InitialiseLoopVariables(); for(Int_t entry=0;entry<numEvents;entry++){ this->SetLoopVariables(entry,0); //avoid divide by zero errors and skip irrelevant data if (rms==0. || mean==0. || numEntries==0) continue; //only look at scintillator strips if (readoutType!=ReadoutType::kScintStrip) continue; if (plane==pl && crate==cr){ if (correlatedHit==1){ if (pulserBox==crate){ hAdcVsChannel[chip]->Fill(channel,mean,20); } else if (pulserBox!=crate){ hAdcVsChannel[chip]->Fill(channel,mean,10); } } else{ hAdcVsChannel[chip]->Fill(channel,mean,1); } if (chip==0 && pulserBox==7){ MSG("LIAnalysis",Msg::kVerbose) <<"plane="<<plane <<", pb="<<pulserBox <<", chip="<<chip <<", pixel="<<pixel <<", channel="<<channel <<", strip="<<strip <<", mean="<<mean <<endl; } else if (chip==0 && pulserBox==6){ MSG("LIAnalysis",Msg::kVerbose) <<"****** plane="<<plane <<", pb="<<pulserBox <<", chip="<<chip <<", pixel="<<pixel <<", channel="<<channel <<", strip="<<strip <<", mean="<<mean <<endl; } } }//end of for MSG("LIAnalysis",Msg::kInfo)<<"Finished main loop"<<endl; TCanvas *cAdcVsChannel=new TCanvas ("cAdcVsChannel","Adc Vs Channel",0,0,1200,800); cAdcVsChannel->SetFillColor(0); cAdcVsChannel->Divide(2,2); for (Int_t i=0;i<NUMCHIPS;i++){ cAdcVsChannel->cd(i+1); hAdcVsChannel[i]->Draw("colz"); } MSG("LIAnalysis",Msg::kInfo) <<" ** Finished the AdcVsChannel method ** "<<endl; }

void LIAnalysis::AdcVsLed (   )

Definition at line 9619 of file LIAnalysis.cxx. References correlatedHit, crate, detectorType, draw(), farPulserBox, Form(), histname, InitialiseLoopVariables(), led, lookup, lowRunNumber, mean, MSG, LILookup::NearOrFar(), nearPulserBox, numEntries, period, plane, pulseHeight, pulserBox, pulses, pulseWidth, rackLevel, readoutType, rms, runNumber, s, SetLoopVariables(), and stripEnd. { MSG("LIAnalysis",Msg::kInfo) <<endl<<" ** Running the AdcVsLed method... ** "<<endl; Int_t numAdcBins=150; //histos for different leds TH1F **hAdcLed=0; hAdcLed= new TH1F*[NUMLEDS*NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ //calculate array number Int_t l=i*NUMLEDS+j; //initialse histos sprintf(histname,"Adc Distribution, PB %d, LED %d",i,j+1); hAdcLed[l]=new TH1F(histname,histname,numAdcBins,0,15000); hAdcLed[l]->GetXaxis()->SetTitle("ADC"); hAdcLed[l]->GetXaxis()->CenterTitle(); hAdcLed[l]->GetYaxis()->SetTitle("Number of Entries"); hAdcLed[l]->GetYaxis()->CenterTitle(); hAdcLed[l]->SetFillColor(0); hAdcLed[l]->Fill(1); hAdcLed[l]->SetBit(TH1::kCanRebin); } } TH1F **hAdcLedNearPb=0; hAdcLedNearPb= new TH1F*[NUMLEDS*NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ //calculate array number Int_t l=i*NUMLEDS+j; sprintf(histname,"Adc NearPb Distribution, PB %d, LED %d",i,j+1); hAdcLedNearPb[l]=new TH1F(histname,histname,numAdcBins,0,15000); hAdcLedNearPb[l]->GetXaxis()->SetTitle("ADC"); hAdcLedNearPb[l]->GetXaxis()->CenterTitle(); hAdcLedNearPb[l]->GetYaxis()->SetTitle("Number of Entries"); hAdcLedNearPb[l]->GetYaxis()->CenterTitle(); hAdcLedNearPb[l]->SetFillColor(0); hAdcLedNearPb[l]->SetLineColor(2); hAdcLedNearPb[l]->Fill(1); hAdcLedNearPb[l]->SetBit(TH1::kCanRebin); } } TH1F **hAdcLedFarPb=0; hAdcLedFarPb= new TH1F*[NUMLEDS*NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ //calculate array number Int_t l=i*NUMLEDS+j; sprintf(histname,"Adc FarPb Distribution, PB %d, LED %d",i,j+1); hAdcLedFarPb[l]=new TH1F(histname,histname,numAdcBins,0,15000); hAdcLedFarPb[l]->GetXaxis()->SetTitle("ADC"); hAdcLedFarPb[l]->GetXaxis()->CenterTitle(); hAdcLedFarPb[l]->GetYaxis()->SetTitle("Number of Entries"); hAdcLedFarPb[l]->GetYaxis()->CenterTitle(); hAdcLedFarPb[l]->SetFillColor(0); hAdcLedFarPb[l]->SetLineColor(3); hAdcLedFarPb[l]->Fill(1); hAdcLedFarPb[l]->SetBit(TH1::kCanRebin); } } TH1F **hGain=0; hGain= new TH1F*[NUMLEDS*NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ //calculate array number Int_t l=i*NUMLEDS+j; sprintf(histname,"Gain Distribution, Near Side, PB %d, LED %d", i,j+1); hGain[l]=new TH1F(histname,histname,numAdcBins,0,200); hGain[l]->GetXaxis()->SetTitle("ADC"); hGain[l]->GetXaxis()->CenterTitle(); hGain[l]->GetYaxis()->SetTitle("Number of Entries"); hGain[l]->GetYaxis()->CenterTitle(); hGain[l]->SetFillColor(0); hGain[l]->SetLineColor(2); hGain[l]->SetBit(TH1::kCanRebin); } } TH1F **hGainUE=0; hGainUE= new TH1F*[NUMLEDS*NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ //calculate array number Int_t l=i*NUMLEDS+j; sprintf(histname,"Gain Distribution UE, PB %d, LED %d",i,j+1); hGainUE[l]=new TH1F(histname,histname,numAdcBins,0,200); hGainUE[l]->GetXaxis()->SetTitle("ADC"); hGainUE[l]->GetXaxis()->CenterTitle(); hGainUE[l]->GetYaxis()->SetTitle("Number of Entries"); hGainUE[l]->GetYaxis()->CenterTitle(); hGainUE[l]->SetFillColor(0); hGainUE[l]->SetLineColor(1); hGainUE[l]->SetBit(TH1::kCanRebin); } } TH1F **hGainUW=0; hGainUW= new TH1F*[NUMLEDS*NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ //calculate array number Int_t l=i*NUMLEDS+j; sprintf(histname,"Gain Distribution UW, PB %d, LED %d", i,j+1); hGainUW[l]=new TH1F(histname,histname,numAdcBins,0,200); hGainUW[l]->GetXaxis()->SetTitle("ADC"); hGainUW[l]->GetXaxis()->CenterTitle(); hGainUW[l]->GetYaxis()->SetTitle("Number of Entries"); hGainUW[l]->GetYaxis()->CenterTitle(); hGainUW[l]->SetFillColor(0); hGainUW[l]->SetLineColor(2); hGainUW[l]->SetBit(TH1::kCanRebin); } } TH1F **hGainLE=0; hGainLE= new TH1F*[NUMLEDS*NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ //calculate array number Int_t l=i*NUMLEDS+j; sprintf(histname,"Gain Distribution LE, PB %d, LED %d", i,j+1); hGainLE[l]=new TH1F(histname,histname,numAdcBins,0,200); hGainLE[l]->GetXaxis()->SetTitle("ADC"); hGainLE[l]->GetXaxis()->CenterTitle(); hGainLE[l]->GetYaxis()->SetTitle("Number of Entries"); hGainLE[l]->GetYaxis()->CenterTitle(); hGainLE[l]->SetFillColor(0); hGainLE[l]->SetLineColor(3); hGainLE[l]->SetBit(TH1::kCanRebin); } } TH1F **hGainLW=0; hGainLW= new TH1F*[NUMLEDS*NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ //calculate array number Int_t l=i*NUMLEDS+j; sprintf(histname,"Gain Distribution LW, Near Side, PB %d, LED %d", i,j+1); hGainLW[l]=new TH1F(histname,histname,numAdcBins,0,200); hGainLW[l]->GetXaxis()->SetTitle("ADC"); hGainLW[l]->GetXaxis()->CenterTitle(); hGainLW[l]->GetYaxis()->SetTitle("Number of Entries"); hGainLW[l]->GetYaxis()->CenterTitle(); hGainLW[l]->SetFillColor(0); hGainLW[l]->SetLineColor(4); hGainLW[l]->SetBit(TH1::kCanRebin); } } TH1F **hGainF=0; hGainF= new TH1F*[NUMLEDS*NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ //calculate array number Int_t l=i*NUMLEDS+j; sprintf(histname,"Gain Distribution, Far Side, PB %d, LED %d", i,j+1); hGainF[l]=new TH1F(histname,histname,numAdcBins,0,200); hGainF[l]->GetXaxis()->SetTitle("ADC"); hGainF[l]->GetXaxis()->CenterTitle(); hGainF[l]->GetYaxis()->SetTitle("Number of Entries"); hGainF[l]->GetYaxis()->CenterTitle(); hGainF[l]->SetFillColor(0); hGainF[l]->SetLineColor(3); hGainF[l]->SetBit(TH1::kCanRebin); } } TH1F **hNpe=0; hNpe= new TH1F*[NUMLEDS*NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ //calculate array number Int_t l=i*NUMLEDS+j; sprintf(histname,"Npe Distribution, Near Side, PB %d, LED %d", i,j+1); hNpe[l]=new TH1F(histname,histname,numAdcBins,0,300); hNpe[l]->GetXaxis()->SetTitle("Npe"); hNpe[l]->GetXaxis()->CenterTitle(); hNpe[l]->GetYaxis()->SetTitle("Number of Entries"); hNpe[l]->GetYaxis()->CenterTitle(); hNpe[l]->SetFillColor(0); hNpe[l]->SetLineColor(2); hNpe[l]->SetBit(TH1::kCanRebin); } } TH1F **hNpeF=0; hNpeF= new TH1F*[NUMLEDS*NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ //calculate array number Int_t l=i*NUMLEDS+j; sprintf(histname,"Npe Distribution, Far Side, PB %d, LED %d", i,j+1); hNpeF[l]=new TH1F(histname,histname,numAdcBins,0,300); hNpeF[l]->GetXaxis()->SetTitle("Npe"); hNpeF[l]->GetXaxis()->CenterTitle(); hNpeF[l]->GetYaxis()->SetTitle("Number of Entries"); hNpeF[l]->GetYaxis()->CenterTitle(); hNpeF[l]->SetFillColor(0); hNpeF[l]->SetLineColor(3); hNpeF[l]->SetBit(TH1::kCanRebin); } } TH1F **hFlash=0; hFlash= new TH1F*[NUMLEDS*NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ //calculate array number Int_t l=i*NUMLEDS+j; sprintf(histname,"Number of Flashes Distribution, Near Side, PB %d, LED %d", i,j+1); hFlash[l]=new TH1F(histname,histname,numAdcBins,0,1500); hFlash[l]->GetXaxis()->SetTitle("Number of Flashes"); hFlash[l]->GetXaxis()->CenterTitle(); hFlash[l]->GetYaxis()->SetTitle("Number of Entries"); hFlash[l]->GetYaxis()->CenterTitle(); hFlash[l]->SetFillColor(0); hFlash[l]->SetLineColor(2); hFlash[l]->SetBit(TH1::kCanRebin); } } TH1F **hRms=0; hRms= new TH1F*[NUMLEDS*NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ //calculate array number Int_t l=i*NUMLEDS+j; sprintf(histname,"RMS Distribution, Near Side, PB %d, LED %d", i,j+1); hRms[l]=new TH1F(histname,histname,numAdcBins,0,15000); hRms[l]->GetXaxis()->SetTitle("RMS"); hRms[l]->GetXaxis()->CenterTitle(); hRms[l]->GetYaxis()->SetTitle("Number of Entries"); hRms[l]->GetYaxis()->CenterTitle(); hRms[l]->SetFillColor(0); hRms[l]->SetLineColor(2); hRms[l]->SetBit(TH1::kCanRebin); } } //array to hold the pulse heights Int_t* ph=new Int_t[NUMPULSERBOXES*NUMLEDS]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ Int_t l=i*NUMLEDS+j; ph[l]=0; } } this->InitialiseLoopVariables(); for(Int_t entry=0;entry<numEvents;entry++){ this->SetLoopVariables(entry,0); //ignore any zeros that slipped through if (mean==0 || rms==0 || numEntries==0) continue; //only look at scint strips if (readoutType!=ReadoutType::kScintStrip) continue; //calculate npe and gain Float_t npe=mean*mean/(rms*rms); Float_t gain=0.8*mean/npe; // PMT fudge factor for M64 (near detector): 0.844 if (detectorType==Detector::kNear) gain = 0.844*rms*rms/mean; //led index Int_t l=pulserBox*NUMLEDS+led-1; if (correlatedHit==1){ //fill histograms if (detectorType==Detector::kFar || detectorType==Detector::kNear){ if (pulserBox>=FIRSTPULSERBOX && pulserBox<=LASTPULSERBOX && led>=FIRSTLED && led<=LASTLED){ //fill array to store pulse heights ph[l]=pulseHeight; //fill histo for appropriate led hAdcLed[l]->Fill(mean); //fill histos when pulser box is near pulser box if (nearPulserBox==pulserBox) { hAdcLedNearPb[l]->Fill(mean); if (mean>500 && mean<8000){ hGain[l]->Fill(gain); hNpe[l]->Fill(npe); hFlash[l]->Fill(numEntries); hRms[l]->Fill(rms); //fill histos if (rackLevel==0 && stripEnd==StripEnd::kEast){ hGainLE[l]->Fill(gain); } if (rackLevel==0 && stripEnd==StripEnd::kWest){ hGainLW[l]->Fill(gain); } if (rackLevel==1 && stripEnd==StripEnd::kEast){ hGainUE[l]->Fill(gain); } if (rackLevel==1 && stripEnd==StripEnd::kWest){ hGainUW[l]->Fill(gain); } } } //fill histos when pulser box is far pulser box else if (farPulserBox==pulserBox) { hAdcLedFarPb[l]->Fill(mean); if (mean>500 && mean<8000){ hGainF[l]->Fill(gain); hNpeF[l]->Fill(npe); } } } } //change the farPb to farLed for calDet else if (detectorType==Detector::kCalDet){ if (pulserBox>=FIRSTPULSERBOX && pulserBox<=LASTPULSERBOX && crate>=FIRSTCRATE && crate<=LASTCRATE && led>=FIRSTLED && led<=LASTLED){ //fill array to store pulse heights ph[l]=pulseHeight; //fill histo for appropriate led hAdcLed[l]->Fill(mean); //fill histos when pulser box is near pulser box if (lookup.NearOrFar(crate,pulserBox,nearPulserBox, farPulserBox,led,detectorType, plane,runNumber)== LILookup::kNearSide) { hAdcLedNearPb[l]->Fill(mean); if (mean>500 && mean<8000){ hGain[l]->Fill(gain); hNpe[l]->Fill(npe); hFlash[l]->Fill(numEntries); hRms[l]->Fill(rms); } } //fill histos when pulser box is far pulser box else if (lookup.NearOrFar(crate,pulserBox,nearPulserBox, farPulserBox,led,detectorType, plane,runNumber)== LILookup::kFarSide) { hAdcLedFarPb[l]->Fill(mean); if (mean>500 && mean<8000){ hGainF[l]->Fill(gain); hNpeF[l]->Fill(npe); } } } } } }//end of for MSG("LIAnalysis",Msg::kInfo)<<"Finished main loop"<<endl; //include the under and overflow counts gStyle->SetOptStat(1111111); //set up useful string string sRunNumber=Form("%d",runNumber); string sLowRunNumber=Form("%d",lowRunNumber); //get the maximum in order to scale all plots Double_t maxNumEntAdc=0; Double_t maxNumEntGain=0; Double_t maxNumEntNpe=0; Double_t maxNumEntFlash=0; Double_t maxNumEntRms=0; for (Int_t i=0;i<NUMPULSERBOXES*NUMLEDS;i++){ if (hAdcLed[i]->GetMean()>500){ //find the maxs for each type of histo if (hAdcLed[i]->GetMaximum()>maxNumEntAdc){ maxNumEntAdc=hAdcLed[i]->GetMaximum(); MSG("LIAnalysis",Msg::kInfo) <<"Calculating max num ADC entries, current highest=" <<maxNumEntAdc<<", l="<<i<<endl; } if (hGain[i]->GetMaximum()>maxNumEntGain){ maxNumEntGain=hGain[i]->GetMaximum(); MSG("LIAnalysis",Msg::kInfo) <<"Calculating max num gain entries, current highest=" <<maxNumEntGain<<", l="<<i<<endl; } if (hNpe[i]->GetMaximum()>maxNumEntNpe){ maxNumEntNpe=hNpe[i]->GetMaximum(); MSG("LIAnalysis",Msg::kInfo) <<"Calculating max num npe entries, current highest=" <<maxNumEntNpe<<", l="<<i<<endl; } if (hFlash[i]->GetMaximum()>maxNumEntFlash){ maxNumEntFlash=hFlash[i]->GetMaximum(); MSG("LIAnalysis",Msg::kInfo) <<"Calculating max num flash entries, current highest=" <<maxNumEntFlash<<", l="<<i<<endl; } if (hRms[i]->GetMaximum()>maxNumEntRms){ maxNumEntRms=hRms[i]->GetMaximum(); MSG("LIAnalysis",Msg::kInfo) <<"Calculating max num RMS entries, current highest=" <<maxNumEntRms<<", l="<<i<<endl; } } } //allocate space for TGraphs then loop and fill them below TGraphAsymmErrors *gErrorsAdcLedNear=new TGraphAsymmErrors(NUMPULSERBOXES*NUMLEDS); TGraphAsymmErrors *gErrorsAdcLedFar=new TGraphAsymmErrors(NUMPULSERBOXES*NUMLEDS); TGraph *gAdcLedNear=new TGraph(NUMPULSERBOXES*NUMLEDS); TGraph *gAdcLedFar=new TGraph(NUMPULSERBOXES*NUMLEDS); TGraph *gGain=new TGraph(NUMPULSERBOXES*NUMLEDS); TGraph *gGainF=new TGraph(NUMPULSERBOXES*NUMLEDS); TGraph *gNpe=new TGraph(NUMPULSERBOXES*NUMLEDS); TGraph *gNpeF=new TGraph(NUMPULSERBOXES*NUMLEDS); TGraph *gFlash=new TGraph(NUMPULSERBOXES*NUMLEDS); TGraph *gRms=new TGraph(NUMPULSERBOXES*NUMLEDS); TGraph *gGainUE=new TGraph(7);//((NUMPULSERBOXES-2)/2);//for fardet TGraph *gGainUW=new TGraph(7);//((NUMPULSERBOXES-2)/2);//for fardet TGraph *gGainLE=new TGraph(7);//((NUMPULSERBOXES-2)/2);//for fardet TGraph *gGainLW=new TGraph(7);//((NUMPULSERBOXES-2)/2);//for fardet gErrorsAdcLedNear->SetMinimum(-1); gErrorsAdcLedFar->SetMinimum(-1); gAdcLedNear->SetMinimum(-1); gAdcLedFar->SetMinimum(-1); gErrorsAdcLedNear->SetMaximum(15000); gErrorsAdcLedFar->SetMaximum(15000); gAdcLedNear->SetMaximum(15000); gAdcLedFar->SetMaximum(15000); gGain->SetMinimum(-1); gGainF->SetMinimum(-1); gNpe->SetMinimum(-1); gNpeF->SetMinimum(-1); Double_t nearAdcMax=0; Double_t farAdcMax=0; //create the canvas and draw TCanvas *cAdcLed=new TCanvas("cAdcLed","cAdcLed",0,0,1000,600); cAdcLed->SetFillColor(0); TCanvas *cGainLed=new TCanvas("cGainLed","cGainLed",0,0,1000,600); cGainLed->SetFillColor(0); TCanvas *cNpeLed=new TCanvas("cNpeLed","cNpeLed",0,0,1000,600); cNpeLed->SetFillColor(0); if (detectorType==Detector::kFar){ Float_t* gainUE=new Float_t[NUMPULSERBOXES]; Float_t* gainUW=new Float_t[NUMPULSERBOXES]; Float_t* gainLE=new Float_t[NUMPULSERBOXES]; Float_t* gainLW=new Float_t[NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ MSG("LIAnalysis",Msg::kInfo) <<"Analysing PB="<<i<<endl; gainUE[i]=0; gainUW[i]=0; gainLE[i]=0; gainLW[i]=0; } Float_t maxGain=0; Float_t minGain=100000; for (Int_t i=0;i<NUMPULSERBOXES-2;i++){ MSG("LIAnalysis",Msg::kInfo) <<"Analysing PB="<<i<<endl; for (Int_t j=0;j<NUMLEDS;j++){ Int_t l=i*NUMLEDS+j; gainUE[i]+=hGainUE[l]->GetMean(); gainUW[i]+=hGainUW[l]->GetMean(); gainLE[i]+=hGainLE[l]->GetMean(); gainLW[i]+=hGainLW[l]->GetMean(); } gainUE[i]/=NUMLEDS; gainUW[i]/=NUMLEDS; gainLE[i]/=NUMLEDS; gainLW[i]/=NUMLEDS; if (gainUE[i]>maxGain) maxGain=gainUE[i]; if (gainUW[i]>maxGain) maxGain=gainUW[i]; if (gainLE[i]>maxGain) maxGain=gainLE[i]; if (gainLW[i]>maxGain) maxGain=gainLW[i]; if (i%2==0){ if (gainUE[i]<minGain) minGain=gainUE[i]; if (gainLE[i]<minGain) minGain=gainLE[i]; } if (i%2==1){ if (gainLW[i]<minGain) minGain=gainLW[i]; if (gainUW[i]<minGain) minGain=gainUW[i]; } MSG("LIAnalysis",Msg::kInfo) <<"Max gain="<<maxGain<<", minGain="<<minGain<<endl; } MSG("LIAnalysis",Msg::kInfo) <<"Max gain="<<maxGain<<", minGain="<<minGain<<endl; Float_t z=3.84; Int_t sm2=0; for (Int_t i=0;i<NUMPULSERBOXES-2;i++){ MSG("LIAnalysis",Msg::kInfo) <<"Filling EW UL graphs, PB="<<i<<endl; if (i>=8) sm2=1; //fill gains if (i%2==0){ gGainUE->SetPoint((i+1)/2,(static_cast<Float_t> ((i+1)/2)+1)*z+sm2-2, gainUE[i]); gGainLE->SetPoint((i+1)/2,(static_cast<Float_t> ((i+1)/2)+1)*z+sm2-2, gainLE[i]); MSG("LIAnalysis",Msg::kInfo) <<"East, i="<<i<<", (i+1)/2="<<(i+1)/2 <<", gainUE[i]="<<gainUE[i]<<", gainLE[i]="<<gainLE[i]<<endl; } else if (i%2==1){ gGainUW->SetPoint(i/2,(static_cast<Float_t>(i/2)+1)*z+sm2-2, gainUW[i]); gGainLW->SetPoint(i/2,(static_cast<Float_t>(i/2)+1)*z+sm2-2, gainLW[i]); MSG("LIAnalysis",Msg::kInfo) <<"West, i="<<i<<", i/2="<<i/2 <<", gainUW[i]="<<gainUW[i]<<", gainLW[i]="<<gainLW[i]<<endl; } } //draw the graphs TCanvas *c=new TCanvas("c","c",0,0,1200,600); c->SetFillColor(0); c->cd(); gGainUE->Draw("APL"); s="Average Gain vs Z position"; gGainUE->SetTitle(s.c_str()); gGainUE->GetXaxis()->SetTitle("Z position (m)"); gGainUE->GetYaxis()->SetTitle("Average Gain (Adcs/p.e.)"); gGainUE->GetXaxis()->CenterTitle(); gGainUE->GetYaxis()->CenterTitle(); gGainUE->SetMarkerStyle(3); gGainUE->SetMarkerColor(1); gGainUE->SetMarkerSize(0.3); gGainUE->SetLineColor(1); gGainUE->SetLineStyle(2); gGainUE->SetMaximum(maxGain+0.07*maxGain); gGainUE->SetMinimum(minGain-0.07*minGain); gGainUW->Draw("PL"); gGainUW->SetMarkerStyle(3); gGainUW->SetMarkerColor(1); gGainUW->SetMarkerSize(0.3); gGainUW->SetLineColor(1); gGainLW->Draw("PL"); gGainLW->SetMarkerStyle(3); gGainLW->SetMarkerColor(2); gGainLW->SetMarkerSize(0.3); gGainLW->SetLineColor(2); gGainLE->Draw("PL"); gGainLE->SetMarkerStyle(3); gGainLE->SetMarkerColor(2); gGainLE->SetMarkerSize(0.3); gGainLE->SetLineColor(2); gGainLE->SetLineStyle(2); } for (Int_t i=0;i<NUMPULSERBOXES;i++){ MSG("LIAnalysis",Msg::kInfo) <<"Analysing PB="<<i<<endl; for (Int_t j=0;j<NUMLEDS;j++){ Int_t l=i*NUMLEDS+j; Double_t nearAdc=hAdcLedNearPb[l]->GetMean(); Double_t farAdc=hAdcLedFarPb[l]->GetMean(); //set points for graphs gErrorsAdcLedNear->SetPoint(l,static_cast<Float_t>(l+1),nearAdc); gErrorsAdcLedFar->SetPoint(l,static_cast<Float_t>(l+1),farAdc); gAdcLedNear->SetPoint(l,static_cast<Float_t>(l+1),nearAdc); gAdcLedFar->SetPoint(l,static_cast<Float_t>(l+1),farAdc); gGain->SetPoint(l,static_cast<Float_t>(l+1),hGain[l]->GetMean()); gGainF->SetPoint(l, static_cast<Float_t>(l+1),hGainF[l]->GetMean()); gNpe->SetPoint(l,static_cast<Float_t>(l+1),hNpe[l]->GetMean()); gNpeF->SetPoint(l,static_cast<Float_t>(l+1),hNpeF[l]->GetMean()); gFlash->SetPoint(l,static_cast<Float_t>(l+1),hFlash[l]->GetMean()); gRms->SetPoint(l,static_cast<Float_t>(l+1),hRms[l]->GetMean()); //calculate the maxs if (nearAdc>nearAdcMax) nearAdcMax=nearAdc; if (farAdc>farAdcMax) farAdcMax=farAdc; //set near PB errors Int_t ent=static_cast<Int_t>(hAdcLedNearPb[l]->GetEntries()); Float_t runningTotal=0; Int_t errorSetL=0; Int_t errorSetH=0; Double_t eyl=0; Double_t eyh=0; for (Int_t k=0;k<hAdcLedNearPb[l]->GetNbinsX();k++){ if (ent==0) break; runningTotal+=hAdcLedNearPb[l]->GetBinContent(k); if (runningTotal/ent>0.05 && errorSetL==0){ eyl=static_cast<Double_t>(hAdcLedNearPb[l]->GetBinLowEdge(k)); if (eyl<hAdcLedNearPb[l]->GetMean()){ eyl=hAdcLedNearPb[l]->GetMean()-eyl; } else{ MSG("LIAnalysis",Msg::kWarning) <<"Near: Lower error greater than mean!" <<"mean="<<hAdcLedNearPb[l]->GetMean()<<", eyl="<<eyl <<endl; eyl=0; } errorSetL=1; MSG("LIAnalysis",Msg::kDebug) <<"Near: Low error set, bin="<<k<<", low edge="<<eyl <<", percentage="<<100*runningTotal/ent<<endl; } else if (runningTotal/ent>0.95 && errorSetH==0){ eyh=static_cast<Double_t>(hAdcLedNearPb[l]->GetBinLowEdge(k)); if (eyh>hAdcLedNearPb[l]->GetMean()){ eyh=eyh-hAdcLedNearPb[l]->GetMean(); } else{ MSG("LIAnalysis",Msg::kWarning) <<"Near: Upper error less than mean!" <<"mean="<<hAdcLedNearPb[l]->GetMean()<<", eyh="<<eyh <<endl; eyh=0; } errorSetH=1; MSG("LIAnalysis",Msg::kDebug) <<"Near: High error set, bin="<<k<<", low edge="<<eyh <<", percentage="<<100*runningTotal/ent<<endl; } } gErrorsAdcLedNear->SetPointError(l,0,0,eyl,eyh); MSG("LIAnalysis",Msg::kDebug) <<"eyl="<<eyl<<", get error="<<*gErrorsAdcLedNear->GetEYlow() <<", eyh="<<eyh<<", get error="<<*gErrorsAdcLedNear->GetEYhigh() <<endl; //set far PB errors ent=static_cast<Int_t>(hAdcLedFarPb[l]->GetEntries()); runningTotal=0; errorSetL=0; errorSetH=0; eyl=0; eyh=0; for (Int_t k=0;k<hAdcLedFarPb[l]->GetNbinsX();k++){ if (ent==1) break; runningTotal+=hAdcLedFarPb[l]->GetBinContent(k); if (runningTotal/ent>0.05 && errorSetL==0){ eyl=static_cast<Double_t>(hAdcLedFarPb[l]->GetBinLowEdge(k)); if (eyl<hAdcLedFarPb[l]->GetMean()){ eyl=hAdcLedFarPb[l]->GetMean()-eyl; } else{ MSG("LIAnalysis",Msg::kWarning) <<"Far: Lower error greater than mean!" <<"mean="<<hAdcLedFarPb[l]->GetMean() <<", eyl="<<eyl<<endl; eyl=0; } errorSetL=1; MSG("LIAnalysis",Msg::kDebug) <<"Far: Low error set, bin="<<k<<", low edge="<<eyl <<", percentage="<<100*runningTotal/ent<<endl; } else if (runningTotal/ent>0.95 && errorSetH==0){ eyh=static_cast<Double_t>(hAdcLedFarPb[l]->GetBinLowEdge(k)); if (eyh>hAdcLedFarPb[l]->GetMean()){ eyh=eyh-hAdcLedFarPb[l]->GetMean(); } else{ MSG("LIAnalysis",Msg::kWarning) <<"Far: Upper error less than mean!" <<"mean="<<hAdcLedFarPb[l]->GetMean() <<", eyh="<<eyh<<endl; eyh=0; } errorSetH=1; MSG("LIAnalysis",Msg::kDebug) <<"Far: High error set, bin="<<k<<", low edge="<<eyh <<", percentage="<<100*runningTotal/ent<<endl; } } gErrorsAdcLedFar->SetPointError(l,0,0,eyl,eyh); Int_t draw=0; if (draw==1){ string sPulserBox=Form("%d",i); string sLed=Form("%d",j+1); string sPulseWidth=Form("%d",pulseWidth); string sPulseHeight=Form("%d",ph[l]); string sPulseFreq=Form("%d",static_cast<Int_t> (ceil(1.0/(period*1.0e-5)))); string sConstantBit=", PB="+sPulserBox+ ", LED="+sLed+", PH="+sPulseHeight+", PW="+sPulseWidth+ ", PF="+sPulseFreq+" Hz"; MSG("LIAnalysis",Msg::kDebug) <<"LI parameters: "<<", PB="<<sPulserBox <<", LED="<<sLed<<", PH="<<sPulseHeight <<", PW="<<sPulseWidth<<", PF="<<sPulseFreq<<endl; //set histo titles s="ADC Values, Near&Far side"+sConstantBit; hAdcLed[l]->SetTitle(s.c_str()); s="Gain Values, Near side"+sConstantBit; hGain[l]->SetTitle(s.c_str()); s="Npe Values, Near side"+sConstantBit; hNpe[l]->SetTitle(s.c_str()); s="Flash Values, Near side"+sConstantBit; hFlash[l]->SetTitle(s.c_str()); s="RMS Values, Near side"+sConstantBit; hRms[l]->SetTitle(s.c_str()); //set maximums hAdcLed[l]->SetMaximum(maxNumEntAdc); hGain[l]->SetMaximum(maxNumEntGain); hNpe[l]->SetMaximum(maxNumEntNpe); hFlash[l]->SetMaximum(maxNumEntFlash); hRms[l]->SetMaximum(maxNumEntRms); cAdcLed->Clear(); cAdcLed->cd(); if (sLowRunNumber==sRunNumber) s=sRunNumber+ "AdcLedHisto.ps"; else s=sLowRunNumber+"-"+sRunNumber+"AdcLedHisto.ps"; if (i*NUMLEDS+j==0){ s+="("; gErrorIgnoreLevel=1; } else if(i*NUMLEDS+j==NUMLEDS*NUMPULSERBOXES-1){ s+=")"; gErrorIgnoreLevel=0; MSG("LIAnalysis",Msg::kInfo)<<"Done last plot"<<endl; } hAdcLed[l]->Draw(); hAdcLedNearPb[l]->Draw("same"); hAdcLedFarPb[l]->Draw("same"); cAdcLed->Print(s.c_str()); cGainLed->Clear(); cGainLed->cd(); if (sLowRunNumber==sRunNumber) s=sRunNumber+"GainLedHisto.ps"; else s=sLowRunNumber+"-"+sRunNumber+"GainLedHisto.ps"; if (i*NUMLEDS+j==0){ s+="("; gErrorIgnoreLevel=1; } else if(i*NUMLEDS+j==NUMLEDS*NUMPULSERBOXES-1){ s+=")"; gErrorIgnoreLevel=0; MSG("LIAnalysis",Msg::kInfo)<<"Done last plot"<<endl; } hGain[l]->Draw(); cGainLed->Print(s.c_str()); cNpeLed->Clear(); cNpeLed->cd(); if (sLowRunNumber==sRunNumber) s=sRunNumber+"NpeLedHisto.ps"; else s=sLowRunNumber+"-"+sRunNumber+"NpeLedHisto.ps"; if (i*NUMLEDS+j==0){ s+="("; gErrorIgnoreLevel=1; } else if(i*NUMLEDS+j==NUMLEDS*NUMPULSERBOXES-1){ s+=")"; gErrorIgnoreLevel=0; MSG("LIAnalysis",Msg::kInfo)<<"Done last plot"<<endl; } hNpe[l]->Draw(); cNpeLed->Print(s.c_str()); } } } //print out a list of the near and far average adcs for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ Int_t l=i*NUMLEDS+j; Double_t yNear=-1; Double_t yFar=-1; Double_t x=-1; gAdcLedNear->GetPoint(l,x,yNear); gAdcLedFar->GetPoint(l,x,yFar); Double_t npeNear=-1, gainNear=-1; gNpe->GetPoint(l,x,npeNear); gGain->GetPoint(l,x,gainNear); //Double_t npeFar=-1, gainFar=-1; gNpeF->GetPoint(l,x,npeFar); gGainF->GetPoint(l,x,gainFar); Double_t flashNear=-1, rmsNear=-1; gFlash->GetPoint(l,x,flashNear); gRms->GetPoint(l,x,rmsNear); if (yNear>1){ string sPulseWidth=Form("%d",pulseWidth); string sPulseHeight=Form("%d",ph[l]); string sPulseFreq=Form("%d",static_cast<Int_t> (ceil(1.0/(period*1.0e-5)))); string sNF=Form("%d",pulses); string sConstantBit="NF= "+sNF+" PH= "+sPulseHeight; //string sConstantBit="NF= "+sNF+" PF= "+sPulseFreq+" PH= "+sPulseHeight+" PW= "+sPulseWidth; } } } //set strings for use in titles string sPulseWidth=Form("%d",pulseWidth); string sPulseHeight=Form("%d",pulseHeight); string sPulseFreq=Form("%d",static_cast<Int_t> (ceil(1.0/(period*1.0e-5)))); string sConstantBit=", PW="+sPulseWidth+", PF="+sPulseFreq+" Hz)"; //draw the graph without errors TCanvas *cAdcLedGraph=new TCanvas("cAdcLedGraph","cAdcLedGraph", 0,0,1200,600); cAdcLedGraph->SetFillColor(0); cAdcLedGraph->cd(); gAdcLedNear->Draw("AP"); s="Average ADC vs LED (Near Side"+sConstantBit; gAdcLedNear->SetTitle(s.c_str()); gAdcLedNear->GetXaxis()->SetTitle("Pulser Box * NUM LEDS + LED"); gAdcLedNear->GetYaxis()->SetTitle("Average ADC"); gAdcLedNear->GetXaxis()->CenterTitle(); gAdcLedNear->GetYaxis()->CenterTitle(); gAdcLedNear->SetMarkerStyle(3); gAdcLedNear->SetMarkerColor(2); gAdcLedNear->SetMarkerSize(0.2); gAdcLedNear->SetLineColor(46); gAdcLedNear->SetMaximum(15000); /* gAdcLedFar->Draw("P"); s="Average ADC vs LED (Near&Far"+sConstantBit; gAdcLedFar->SetTitle(s.c_str()); gAdcLedFar->SetMarkerStyle(3); gAdcLedFar->SetMarkerColor(3); gAdcLedFar->SetMarkerSize(0.3); gAdcLedFar->SetLineColor(30); */ //print graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"AdcVsLed.ps("; else s=sLowRunNumber+"-"+sRunNumber+"AdcVsLed.ps("; cAdcLedGraph->Print(s.c_str()); //draw the graphs with errors TCanvas *cErrorsAdcLedGraph=new TCanvas("cErrorsAdcLedGraph", "cErrorsAdcLedGraph", 0,0,1200,600); cErrorsAdcLedGraph->SetFillColor(0); cErrorsAdcLedGraph->cd(); gErrorsAdcLedNear->Draw("AP"); s="Average ADC vs LED (Near&Far, Error=+/-5% Cut"+sConstantBit; gErrorsAdcLedNear->SetTitle(s.c_str()); gErrorsAdcLedNear->GetXaxis()->SetTitle("Pulser Box * NUM LEDS + LED"); gErrorsAdcLedNear->GetYaxis()->SetTitle("Average ADC"); gErrorsAdcLedNear->GetXaxis()->CenterTitle(); gErrorsAdcLedNear->GetYaxis()->CenterTitle(); gErrorsAdcLedNear->SetMarkerStyle(3); gErrorsAdcLedNear->SetMarkerColor(2); gErrorsAdcLedNear->SetMarkerSize(0.2); gErrorsAdcLedNear->SetLineColor(46); gErrorsAdcLedFar->Draw("P"); gErrorsAdcLedFar->SetTitle(s.c_str()); gErrorsAdcLedFar->SetMarkerStyle(3); gErrorsAdcLedFar->SetMarkerColor(3); gErrorsAdcLedFar->SetMarkerSize(0.3); gErrorsAdcLedFar->SetLineColor(30); //print graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"AdcVsLed.ps"; else s=sLowRunNumber+"-"+sRunNumber+"AdcVsLed.ps"; cErrorsAdcLedGraph->Print(s.c_str()); //draw graphs separately TCanvas *cErrorsAdcLedGraphNear=new TCanvas("cErrorsAdcLedGraphNear", "cErrorsAdcLedGraphNear", 0,0,1200,600); cErrorsAdcLedGraphNear->SetFillColor(0); cErrorsAdcLedGraphNear->cd(); gErrorsAdcLedNear->Draw("AP"); s="Average ADC vs LED (Near Side, Error=+/-5% Cut"+sConstantBit; gErrorsAdcLedNear->SetTitle(s.c_str()); gErrorsAdcLedNear->SetLineColor(46); //print graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"AdcVsLed.ps"; else s=sLowRunNumber+"-"+sRunNumber+"AdcVsLed.ps"; cErrorsAdcLedGraphNear->Print(s.c_str()); TCanvas *cErrorsAdcLedGraphFar=new TCanvas("cErrorsAdcLedGraphFar", "cErrorsAdcLedGraphFar", 0,0,1200,600); cErrorsAdcLedGraphFar->SetFillColor(0); cErrorsAdcLedGraphFar->cd(); gErrorsAdcLedFar->SetMaximum();//set to default to trigger recomputing gErrorsAdcLedFar->SetMinimum();//of the range gErrorsAdcLedFar->Draw("AP"); s="Average ADC vs LED (Far Side, Error=+/-5% Cut"+sConstantBit; gErrorsAdcLedFar->SetTitle(s.c_str()); gErrorsAdcLedFar->GetXaxis()->SetTitle("Pulser Box * NUM LEDS + LED"); gErrorsAdcLedFar->GetYaxis()->SetTitle("Average ADC"); gErrorsAdcLedFar->GetXaxis()->CenterTitle(); gErrorsAdcLedFar->GetYaxis()->CenterTitle(); gErrorsAdcLedFar->SetLineColor(30); //print to graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"AdcVsLed.ps)"; else s=sLowRunNumber+"-"+sRunNumber+"AdcVsLed.ps)"; cErrorsAdcLedGraphFar->Print(s.c_str()); //draw gain plots TCanvas *cGain=new TCanvas("cGain","cGain",0,0,1200,600); cGain->SetFillColor(0); cGain->cd(); gGain->Draw("AP"); s="Average Gain vs LED (Near&Far"+sConstantBit; gGain->SetTitle(s.c_str()); gGain->GetXaxis()->SetTitle("Pulser Box * NUM LEDS + LED"); gGain->GetYaxis()->SetTitle("Average Gain"); gGain->GetXaxis()->CenterTitle(); gGain->GetYaxis()->CenterTitle(); gGain->SetMarkerStyle(3); gGain->SetMarkerColor(2); gGain->SetMarkerSize(0.2); gGain->SetLineColor(46); gGain->SetMinimum(-1); gGainF->Draw("P"); s="Average Gain vs LED (Near&Far"+sConstantBit; gGainF->SetTitle(s.c_str()); gGainF->SetMarkerStyle(3); gGainF->SetMarkerColor(3); gGainF->SetMarkerSize(0.2); gGainF->SetLineColor(30); //print to graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"GainVsLed.ps("; else s=sLowRunNumber+"-"+sRunNumber+"GainVsLed.ps("; cGain->Print(s.c_str()); cGain->Clear(); cGain->cd(); gGain->Draw("AP"); s="Average Gain vs LED (Near Side"+sConstantBit; gGain->SetTitle(s.c_str()); gGain->GetXaxis()->SetTitle("Pulser Box * NUM LEDS + LED"); gGain->GetYaxis()->SetTitle("Average Gain"); gGain->GetXaxis()->CenterTitle(); gGain->GetYaxis()->CenterTitle(); gGain->SetMarkerStyle(3); gGain->SetMarkerColor(2); gGain->SetMarkerSize(0.2); gGain->SetLineColor(46); gGain->SetMinimum(-1); //print to graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"GainVsLed.ps"; else s=sLowRunNumber+"-"+sRunNumber+"GainVsLed.ps"; cGain->Print(s.c_str()); cGain->Clear(); cGain->cd(); gGainF->Draw("AP"); s="Average Gain vs LED (Far Side"+sConstantBit; gGainF->SetTitle(s.c_str()); gGainF->GetXaxis()->SetTitle("Pulser Box * NUM LEDS + LED"); gGainF->GetYaxis()->SetTitle("Average Gain"); gGainF->GetXaxis()->CenterTitle(); gGainF->GetYaxis()->CenterTitle(); gGainF->SetMarkerStyle(3); gGainF->SetMarkerColor(3); gGainF->SetMarkerSize(0.2); gGainF->SetLineColor(30); gGainF->SetMinimum(-1); //print to graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"GainVsLed.ps)"; else s=sLowRunNumber+"-"+sRunNumber+"GainVsLed.ps)"; cGain->Print(s.c_str()); //draw npe plots TCanvas *cNpe=new TCanvas("cNpe","cNpe",0,0,1200,600); cNpe->SetFillColor(0); cNpe->cd(); gNpe->Draw("AP"); s="Average Npe vs LED (Near&Far"+sConstantBit; gNpe->SetTitle(s.c_str()); gNpe->GetXaxis()->SetTitle("Pulser Box * NUM LEDS + LED"); gNpe->GetYaxis()->SetTitle("Average Npe"); gNpe->GetXaxis()->CenterTitle(); gNpe->GetYaxis()->CenterTitle(); gNpe->SetMarkerStyle(3); gNpe->SetMarkerColor(2); gNpe->SetMarkerSize(0.2); gNpe->SetLineColor(46); gNpe->SetMinimum(-1); gNpeF->Draw("P"); s="Average Npe vs LED (Near&Far"+sConstantBit; gNpeF->SetTitle(s.c_str()); gNpeF->SetMarkerStyle(3); gNpeF->SetMarkerColor(3); gNpeF->SetMarkerSize(0.2); gNpeF->SetLineColor(30); //print to graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"NpeVsLed.ps("; else s=sLowRunNumber+"-"+sRunNumber+"NpeVsLed.ps("; cNpe->Print(s.c_str()); cNpe->Clear(); cNpe->cd(); gNpe->Draw("AP"); s="Average Npe vs LED (Near Side"+sConstantBit; gNpe->SetTitle(s.c_str()); gNpe->GetXaxis()->SetTitle("Pulser Box * NUM LEDS + LED"); gNpe->GetYaxis()->SetTitle("Average Npe"); gNpe->GetXaxis()->CenterTitle(); gNpe->GetYaxis()->CenterTitle(); gNpe->SetMarkerStyle(3); gNpe->SetMarkerColor(2); gNpe->SetMarkerSize(0.2); gNpe->SetLineColor(46); gNpe->SetMinimum(-1); //print to graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"NpeVsLed.ps"; else s=sLowRunNumber+"-"+sRunNumber+"NpeVsLed.ps"; cNpe->Print(s.c_str()); cNpe->Clear(); cNpe->cd(); gNpeF->Draw("AP"); s="Average Npe vs LED (Far Side"+sConstantBit; gNpeF->SetTitle(s.c_str()); gNpeF->GetXaxis()->SetTitle("Pulser Box * NUM LEDS + LED"); gNpeF->GetYaxis()->SetTitle("Average Npe"); gNpeF->GetXaxis()->CenterTitle(); gNpeF->GetYaxis()->CenterTitle(); gNpeF->SetMarkerStyle(3); gNpeF->SetMarkerColor(3); gNpeF->SetMarkerSize(0.2); gNpeF->SetLineColor(30); gNpeF->SetMinimum(-1); //print to graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"NpeVsLed.ps)"; else s=sLowRunNumber+"-"+sRunNumber+"NpeVsLed.ps)"; cNpe->Print(s.c_str()); MSG("LIAnalysis",Msg::kInfo) <<endl<<" ** Finished the AdcVsLed method ** "<<endl; }

void LIAnalysis::AdcVsPin (   )

Definition at line 1083 of file LIAnalysis.cxx. References chain, chip, correlatedHit, draw(), Form(), GetElecString(), histname, InitialiseLoopVariables(), led, lowRunNumber, mean, MSG, numEntries, period, pinGain, pulseHeight, pulserBox, pulses, pulseWidth, readoutType, rms, runNumber, s, and SetLoopVariables(). { MSG("LIAnalysis",Msg::kInfo) <<endl<<" ** Running the AdcVsPin method... ** "<<endl; Int_t numAdcBins=150; //histos for different pins, high gain TH1F **hAdcHighPin=0; hAdcHighPin= new TH1F*[NUMLEDS*NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ //calculate array number Int_t l=i*NUMLEDS+j; //initialse histos sprintf(histname,"High Gain PIN Adc Distribution, PB %d, LED %d", i,j+1); hAdcHighPin[l]=new TH1F(histname,histname,numAdcBins,0,15000); hAdcHighPin[l]->GetXaxis()->SetTitle("ADC"); hAdcHighPin[l]->GetXaxis()->CenterTitle(); hAdcHighPin[l]->GetYaxis()->SetTitle("Number of Entries"); hAdcHighPin[l]->GetYaxis()->CenterTitle(); hAdcHighPin[l]->SetFillColor(0); hAdcHighPin[l]->SetLineColor(2); hAdcHighPin[l]->Fill(1); //hAdcHighPin[l]->SetBit(TH1::kCanRebin); } } //histos for different pins, low gain TH1F **hAdcLowPin=0; hAdcLowPin= new TH1F*[NUMLEDS*NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ //calculate array number Int_t l=i*NUMLEDS+j; //initialse histos sprintf(histname,"Low Gain PIN Adc Distribution, PB %d, LED %d", i,j+1); hAdcLowPin[l]=new TH1F(histname,histname,numAdcBins,0,15000); hAdcLowPin[l]->GetXaxis()->SetTitle("ADC"); hAdcLowPin[l]->GetXaxis()->CenterTitle(); hAdcLowPin[l]->GetYaxis()->SetTitle("Number of Entries"); hAdcLowPin[l]->GetYaxis()->CenterTitle(); hAdcLowPin[l]->SetFillColor(0); hAdcLowPin[l]->SetLineColor(3); hAdcLowPin[l]->Fill(1); //hAdcLowPin[l]->SetBit(TH1::kCanRebin); } } Float_t* maxAdcHighPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* maxAdcLowPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* rmsHighPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* rmsLowPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* rms2HighPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* rms2LowPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* numHighPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* numLowPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; //initialise arrays for (Int_t i=0;i<NUMLEDS*NUMPULSERBOXES;i++){ maxAdcHighPin[i]=0; maxAdcLowPin[i]=0; rmsHighPin[i]=0; rmsLowPin[i]=0; rms2HighPin[i]=0; rms2LowPin[i]=0; numHighPin[i]=0; numLowPin[i]=0; } this->InitialiseLoopVariables(); for(Int_t entry=0;entry<numEvents;entry++){ this->SetLoopVariables(entry,0); //ignore any zeros if (mean==0 || rms==0 || numEntries==0) continue; //only look at pins if (readoutType!=ReadoutType::kPinDiode) continue; //cut out strange pins with high mean if (numEntries<0.8*pulses) continue; if (pulserBox>=FIRSTPULSERBOX && pulserBox<=LASTPULSERBOX && led>=FIRSTLED && led<=LASTLED){ //define led number Int_t l=pulserBox*NUMLEDS+led-1; if (numEntries<0.5*pulses && mean>1000){ MSG("LIAnalysis",Msg::kInfo) <<"**** Strange pin: ("<<pulserBox<<":"<<led <<") on" <<" "<<this->GetElecString() <<", plex gain="<<pinGain<<", mean="<<mean<<", rms="<<rms <<", num="<<numEntries<<", cHit="<<correlatedHit<<endl; } //will require a correlated hit when plex is working //fill histo for appropriate pin if (chip==1){//high gain (=0) (chip 1) hAdcHighPin[l]->Fill(mean); //find max value for high gain pin if (mean>maxAdcHighPin[l]) { maxAdcHighPin[l]=mean; rmsHighPin[l]=rms; rms2HighPin[l]=rms*rms; numHighPin[l]=numEntries; MSG("LIAnalysis",Msg::kInfo) <<"HG Pin ("<<pulserBox<<":"<<led <<") on" <<" "<<this->GetElecString() <<", plex gain="<<pinGain<<", mean="<<mean<<", rms="<<rms <<", num="<<numEntries<<", cHit="<<correlatedHit<<endl; } } else if (chip==0){//low gain (=1) (chip 0) hAdcLowPin[l]->Fill(mean); //find max value for low gain pin if (mean>maxAdcLowPin[l]) { maxAdcLowPin[l]=mean; rmsLowPin[l]=rms; rms2LowPin[l]=rms*rms; numLowPin[l]=numEntries; MSG("LIAnalysis",Msg::kInfo) <<"LG Pin ("<<pulserBox<<":"<<led <<") on" <<" "<<this->GetElecString() <<", plex gain="<<pinGain<<", mean="<<mean<<", rms="<<rms <<", num="<<numEntries<<", cHit="<<correlatedHit<<endl; } } //print out warnings if there are inconsistencies if ((pinGain!=1 && chip==0) || (pinGain!=0 && chip==1)) { MSG("LIAnalysis",Msg::kDebug) <<"**** Strange pin, wrong gain in plex," <<" "<<this->GetElecString() <<", plex gain="<<pinGain<<", mean="<<mean<<", rms="<<rms <<", num="<<numEntries<<", cHit="<<correlatedHit<<endl; } } }//end of for MSG("LIAnalysis",Msg::kInfo)<<"Finished main loop"<<endl; //include the under and overflow counts gStyle->SetOptStat(1111111); //set up useful string string sRunNumber=Form("%d",runNumber); string sLowRunNumber=Form("%d",lowRunNumber); MSG("LIAnalysis",Msg::kInfo) <<"List of pins:"<<endl; MSG("LIAnalysis",Msg::kInfo) <<" HG pins:"<<endl; for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ Int_t l=i*NUMLEDS+j; MSG("LIAnalysis",Msg::kInfo) <<" ("<<i<<":"<<j+1<<") mean="<<maxAdcHighPin[l] <<", rms="<<rmsHighPin[l] <<", num="<<numHighPin[l]<<endl; } } MSG("LIAnalysis",Msg::kInfo) <<" LG pins:"<<endl; for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ Int_t l=i*NUMLEDS+j; MSG("LIAnalysis",Msg::kInfo) <<" ("<<i<<":"<<j+1<<") mean="<<maxAdcLowPin[l] <<", rms="<<rmsLowPin[l] <<", num="<<numLowPin[l]<<endl; } } //get the maximum in order to scale all plots Int_t maxNumEnt=0; for (Int_t i=0;i<NUMPULSERBOXES*NUMLEDS;i++){ if (hAdcHighPin[i]->GetMaximum()>maxNumEnt){ maxNumEnt=static_cast<Int_t>(hAdcHighPin[i]->GetMaximum()); MSG("LIAnalysis",Msg::kInfo) <<"Calculating max number of entries, current highest=" <<maxNumEnt<<", l="<<i<<endl; } } for (Int_t i=0;i<NUMPULSERBOXES*NUMLEDS;i++){ if (hAdcLowPin[i]->GetMaximum()>maxNumEnt){ maxNumEnt=static_cast<Int_t>(hAdcLowPin[i]->GetMaximum()); MSG("LIAnalysis",Msg::kInfo) <<"Calculating max number of entries (low), current highest=" <<maxNumEnt<<", l="<<i<<endl; } } //allocate space for TGraphs then loop and fill them below TGraphAsymmErrors *gErrorsAdcHighPin=new TGraphAsymmErrors(NUMPULSERBOXES*NUMLEDS); TGraphAsymmErrors *gErrorsAdcLowPin=new TGraphAsymmErrors(NUMPULSERBOXES*NUMLEDS); TGraph *gAdcHighPin=new TGraph(NUMPULSERBOXES*NUMLEDS); TGraph *gAdcLowPin=new TGraph(NUMPULSERBOXES*NUMLEDS); TGraph *gRmsHighPin=new TGraph(NUMPULSERBOXES*NUMLEDS); TGraph *gRmsLowPin=new TGraph(NUMPULSERBOXES*NUMLEDS); TGraph *gRms2HighPin=new TGraph(NUMPULSERBOXES*NUMLEDS); TGraph *gRms2LowPin=new TGraph(NUMPULSERBOXES*NUMLEDS); TGraph *gNumHighPin=new TGraph(NUMPULSERBOXES*NUMLEDS); TGraph *gNumLowPin=new TGraph(NUMPULSERBOXES*NUMLEDS); TGraph *gResHighPin=new TGraph(NUMPULSERBOXES*NUMLEDS); TGraph *gResLowPin=new TGraph(NUMPULSERBOXES*NUMLEDS); gAdcHighPin->SetMinimum(0); gAdcLowPin->SetMinimum(0); gRmsHighPin->SetMinimum(0); gRmsLowPin->SetMinimum(0); gRms2HighPin->SetMinimum(0); gRms2LowPin->SetMinimum(0); gNumHighPin->SetMinimum(0); gNumLowPin->SetMinimum(0); gResHighPin->SetMinimum(0); gResLowPin->SetMinimum(0); //create the canvas and draw TCanvas *cAdcPin=new TCanvas("cAdcPin","cAdcPin",0,0,1000,600); cAdcPin->SetFillColor(0); cAdcPin->cd(); for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ Int_t l=i*NUMLEDS+j; //set points for graphs gErrorsAdcHighPin->SetPoint(l,static_cast<Float_t>(l+1), maxAdcHighPin[l]); gErrorsAdcLowPin->SetPoint(l,static_cast<Float_t>(l+1), maxAdcLowPin[l]); gAdcHighPin->SetPoint(l,static_cast<Float_t>(l+1), maxAdcHighPin[l]); gAdcLowPin->SetPoint(l,static_cast<Float_t>(l+1),maxAdcLowPin[l]); gRmsHighPin->SetPoint(l,static_cast<Float_t>(l+1),rmsHighPin[l]); gRmsLowPin->SetPoint(l,static_cast<Float_t>(l+1),rmsLowPin[l]); gRms2HighPin->SetPoint(l,static_cast<Float_t>(l+1), rms2HighPin[l]); gRms2LowPin->SetPoint(l,static_cast<Float_t>(l+1),rms2LowPin[l]); gNumHighPin->SetPoint(l,static_cast<Float_t>(l+1),numHighPin[l]); gNumLowPin->SetPoint(l,static_cast<Float_t>(l+1),numLowPin[l]); if (maxAdcHighPin[l]>0){ gResHighPin->SetPoint(l,static_cast<Float_t>(l+1), rmsHighPin[l]/maxAdcHighPin[l]); } else { gResHighPin->SetPoint(l,static_cast<Float_t>(l+1),0); } if (maxAdcLowPin[l]>0){ gResLowPin->SetPoint(l,static_cast<Float_t>(l+1), rmsLowPin[l]/maxAdcLowPin[l]); } else { gResLowPin->SetPoint(l,static_cast<Float_t>(l+1),0); } //plot histograms cAdcPin->Clear(); hAdcHighPin[l]->SetMaximum(maxNumEnt); string sPulserBox=Form("%d",i); string sLed=Form("%d",j+1); string sPulseWidth=Form("%d",pulseWidth); string sPulseHeight=Form("%d",pulseHeight); string sPulseFreq=Form("%d",static_cast<Int_t> (ceil(1.0/(period*1.0e-5)))); MSG("LIAnalysis",Msg::kDebug) <<"LI parameters: " <<", PB="<<sPulserBox <<", LED="<<sLed <<", PH="<<sPulseHeight <<", PW="<<sPulseWidth <<", PF="<<sPulseFreq <<endl; s="Pin ADC Values (High&Low Gain, PB="+sPulserBox+ ", LED="+sLed+", PH="+sPulseHeight+", PW="+sPulseWidth+ ", PF="+sPulseFreq+" Hz)"; hAdcHighPin[l]->SetTitle(s.c_str()); Int_t draw=1; if (draw==1){ hAdcHighPin[l]->Draw(); hAdcLowPin[l]->Draw("same"); if (i*NUMLEDS+j==0){ if (sLowRunNumber==sRunNumber) s=sRunNumber+ "PinAdcHisto.ps("; else s=sLowRunNumber+"-"+sRunNumber+"PinAdcHisto.ps("; cAdcPin->Print(s.c_str()); gErrorIgnoreLevel=1; } else if(i*NUMLEDS+j==NUMLEDS*NUMPULSERBOXES-1){ gErrorIgnoreLevel=0; if (sLowRunNumber==sRunNumber) s=sRunNumber+ "PinAdcHisto.ps)"; else s=sLowRunNumber+"-"+sRunNumber+"PinAdcHisto.ps)"; cAdcPin->Print(s.c_str()); MSG("LIAnalysis",Msg::kInfo)<<"Done last plot"<<endl; } else{ if (sLowRunNumber==sRunNumber) s=sRunNumber+ "PinAdcHisto.ps"; else s=sLowRunNumber+"-"+sRunNumber+"PinAdcHisto.ps"; cAdcPin->Print(s.c_str()); } } } } cAdcPin->Clear(); hAdcHighPin[0]->Draw(); hAdcLowPin[0]->Draw("same"); //set strings for use in titles string sPulseWidth=Form("%d",pulseWidth); string sPulseHeight=Form("%d",pulseHeight); string sPulseFreq=Form("%d",static_cast<Int_t> (ceil(1.0/(period*1.0e-5)))); string sPulses=Form("%d",pulses); string sConstantBit=""; if (runNumber>=13123){//don't use pulse height Int_t fph=pulseHeight; chain->GetEvent(numEvents-1); Int_t lph=-10;//pulseHeight; MSG("LIAnalysis",Msg::kInfo) <<"first ph="<<fph<<", last ph="<<lph<<endl; if (lph==fph){ sConstantBit=+", PH="+sPulseHeight+", PW="+sPulseWidth+ ", PF="+sPulseFreq+" Hz, PN="+sPulses+")"; } else if (lph!=fph){ sConstantBit=+", PW="+sPulseWidth+ ", PF="+sPulseFreq+" Hz, PN="+sPulses+")"; } } else{ sConstantBit=+", PH="+sPulseHeight+", PW="+sPulseWidth+ ", PF="+sPulseFreq+" Hz, PN="+sPulses+")"; } //draw the adc graph TCanvas *cAdcPinGraph=new TCanvas("cAdcPinGraph","cAdcPinGraph", 0,0,1200,600); cAdcPinGraph->SetFillColor(0); cAdcPinGraph->cd(); gAdcHighPin->Draw("AP"); s="Max Pin ADC vs LED (High&Low"+sConstantBit; gAdcHighPin->SetTitle(s.c_str()); gAdcHighPin->GetXaxis()->SetTitle("Pulser Box * NUM PINS + LED"); gAdcHighPin->GetYaxis()->SetTitle("Average ADC"); gAdcHighPin->GetXaxis()->CenterTitle(); gAdcHighPin->GetYaxis()->CenterTitle(); gAdcHighPin->SetMarkerStyle(3); gAdcHighPin->SetMarkerColor(2); gAdcHighPin->SetMarkerSize(0.2); gAdcHighPin->SetLineColor(46); gAdcLowPin->Draw("P"); gAdcLowPin->SetTitle(s.c_str()); gAdcLowPin->SetMarkerStyle(3); gAdcLowPin->SetMarkerColor(3); gAdcLowPin->SetMarkerSize(0.3); gAdcLowPin->SetLineColor(30); //print graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"PinAdcVsLed.ps("; else s=sLowRunNumber+"-"+sRunNumber+"PinAdcVsLed.ps("; //cAdcPinGraph->Modified(); gAdcHighPin->SetMinimum(0); cAdcPinGraph->Print(s.c_str()); //draw the adc high graph cAdcPinGraph->Clear(); cAdcPinGraph->cd(); //gAdcHighPin->SetMaximum(4000); gAdcHighPin->Draw("AP"); s="Max Pin Adc vs LED (High Gain"+sConstantBit; gAdcHighPin->SetTitle(s.c_str()); gAdcHighPin->GetXaxis()->SetTitle("Pulser Box * NUMLEDS + Led"); gAdcHighPin->GetYaxis()->SetTitle("Average Adc"); gAdcHighPin->GetXaxis()->CenterTitle(); gAdcHighPin->GetYaxis()->CenterTitle(); gAdcHighPin->SetMarkerStyle(3); gAdcHighPin->SetMarkerColor(2); gAdcHighPin->SetMarkerSize(0.2); gAdcHighPin->SetLineColor(46); //print graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"PinAdcVsLed.ps"; else s=sLowRunNumber+"-"+sRunNumber+"PinAdcVsLed.ps"; gAdcHighPin->SetMinimum(0); cAdcPinGraph->Print(s.c_str()); //draw the adc low graph cAdcPinGraph->Clear(); cAdcPinGraph->cd(); //gAdcLowPin->SetMaximum(2000); gAdcLowPin->Draw("AP"); s="Max Pin Adc vs LED (Low Gain"+sConstantBit; gAdcLowPin->SetTitle(s.c_str()); gAdcLowPin->GetXaxis()->SetTitle("Pulser Box * NUMLEDS + Led"); gAdcLowPin->GetYaxis()->SetTitle("Average Adc"); gAdcLowPin->GetXaxis()->CenterTitle(); gAdcLowPin->GetYaxis()->CenterTitle(); gAdcLowPin->SetMarkerStyle(3); gAdcLowPin->SetMarkerColor(3); gAdcLowPin->SetMarkerSize(0.2); gAdcLowPin->SetLineColor(30); //print graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"PinAdcVsLed.ps)"; else s=sLowRunNumber+"-"+sRunNumber+"PinAdcVsLed.ps)"; gAdcLowPin->SetMinimum(0); cAdcPinGraph->Print(s.c_str()); //draw the rms graph TCanvas *cRmsPinGraph=new TCanvas("cRmsPinGraph","cRmsPinGraph", 0,0,1200,600); cRmsPinGraph->SetFillColor(0); cRmsPinGraph->cd(); gRmsHighPin->Draw("AP"); s="Pin RMS vs LED (High&Low Gain"+sConstantBit; gRmsHighPin->SetTitle(s.c_str()); gRmsHighPin->GetXaxis()->SetTitle("Pulser Box * NUMLEDS + Led"); gRmsHighPin->GetYaxis()->SetTitle("Average Rms"); gRmsHighPin->GetXaxis()->CenterTitle(); gRmsHighPin->GetYaxis()->CenterTitle(); gRmsHighPin->SetMarkerStyle(3); gRmsHighPin->SetMarkerColor(2); gRmsHighPin->SetMarkerSize(0.2); gRmsHighPin->SetLineColor(46); gRmsLowPin->Draw("P"); gRmsLowPin->SetTitle(s.c_str()); gRmsLowPin->SetMarkerStyle(3); gRmsLowPin->SetMarkerColor(3); gRmsLowPin->SetMarkerSize(0.3); gRmsLowPin->SetLineColor(30); //print graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"PinRmsVsLed.ps("; else s=sLowRunNumber+"-"+sRunNumber+"PinRmsVsLed.ps("; gRmsHighPin->SetMinimum(0); cRmsPinGraph->Print(s.c_str()); //draw the rms high graph cRmsPinGraph->Clear(); cRmsPinGraph->cd(); gRmsHighPin->SetMaximum(100); gRmsHighPin->Draw("AP"); s="Pin RMS vs LED (High Gain"+sConstantBit; gRmsHighPin->SetTitle(s.c_str()); gRmsHighPin->GetXaxis()->SetTitle("Pulser Box * NUMLEDS + Led"); gRmsHighPin->GetYaxis()->SetTitle("Average Rms"); gRmsHighPin->GetXaxis()->CenterTitle(); gRmsHighPin->GetYaxis()->CenterTitle(); gRmsHighPin->SetMarkerStyle(3); gRmsHighPin->SetMarkerColor(2); gRmsHighPin->SetMarkerSize(0.2); gRmsHighPin->SetLineColor(46); //print graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"PinRmsVsLed.ps"; else s=sLowRunNumber+"-"+sRunNumber+"PinRmsVsLed.ps"; gRmsHighPin->SetMinimum(0); cRmsPinGraph->Print(s.c_str()); //draw the rms low graph cRmsPinGraph->Clear(); cRmsPinGraph->cd(); gRmsLowPin->SetMaximum(100); gRmsLowPin->Draw("AP"); s="Pin RMS vs LED (Low Gain"+sConstantBit; gRmsLowPin->SetTitle(s.c_str()); gRmsLowPin->GetXaxis()->SetTitle("Pulser Box * NUMLEDS + Led"); gRmsLowPin->GetYaxis()->SetTitle("Average Rms"); gRmsLowPin->GetXaxis()->CenterTitle(); gRmsLowPin->GetYaxis()->CenterTitle(); gRmsLowPin->SetMarkerStyle(3); gRmsLowPin->SetMarkerColor(3); gRmsLowPin->SetMarkerSize(0.2); gRmsLowPin->SetLineColor(30); //print graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"PinRmsVsLed.ps)"; else s=sLowRunNumber+"-"+sRunNumber+"PinRmsVsLed.ps)"; gRmsLowPin->SetMinimum(0); cRmsPinGraph->Print(s.c_str()); //draw the rms2 graph TCanvas *cRms2PinGraph=new TCanvas("cRms2PinGraph","cRms2PinGraph", 0,0,1200,600); cRms2PinGraph->SetFillColor(0); cRms2PinGraph->cd(); gRms2HighPin->Draw("AP"); s="Pin RMS Sqd vs LED (High&Low Gain"+sConstantBit; gRms2HighPin->SetTitle(s.c_str()); gRms2HighPin->GetXaxis()->SetTitle("Pulser Box * NUMLEDS + Led"); gRms2HighPin->GetYaxis()->SetTitle("Average Rms2"); gRms2HighPin->GetXaxis()->CenterTitle(); gRms2HighPin->GetYaxis()->CenterTitle(); gRms2HighPin->SetMarkerStyle(3); gRms2HighPin->SetMarkerColor(2); gRms2HighPin->SetMarkerSize(0.2); gRms2HighPin->SetLineColor(46); gRms2LowPin->Draw("P"); gRms2LowPin->SetTitle(s.c_str()); gRms2LowPin->SetMarkerStyle(3); gRms2LowPin->SetMarkerColor(3); gRms2LowPin->SetMarkerSize(0.3); gRms2LowPin->SetLineColor(30); //print graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"PinRms2VsLed.ps("; else s=sLowRunNumber+"-"+sRunNumber+"PinRms2VsLed.ps("; cRms2PinGraph->Print(s.c_str()); //draw the rms2 high graph cRms2PinGraph->Clear(); cRms2PinGraph->cd(); //gRms2HighPin->SetMaximum(1500); gRms2HighPin->Draw("AP"); s="Pin RMS Sqd vs LED (High Gain"+sConstantBit; gRms2HighPin->SetTitle(s.c_str()); gRms2HighPin->GetXaxis()->SetTitle("Pulser Box * NUMLEDS + Led"); gRms2HighPin->GetYaxis()->SetTitle("Average Rms2"); gRms2HighPin->GetXaxis()->CenterTitle(); gRms2HighPin->GetYaxis()->CenterTitle(); gRms2HighPin->SetMarkerStyle(3); gRms2HighPin->SetMarkerColor(2); gRms2HighPin->SetMarkerSize(0.2); gRms2HighPin->SetLineColor(46); //print graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"PinRms2VsLed.ps"; else s=sLowRunNumber+"-"+sRunNumber+"PinRms2VsLed.ps"; cRms2PinGraph->Print(s.c_str()); //draw the rms2 low graph cRms2PinGraph->Clear(); cRms2PinGraph->cd(); //gRms2LowPin->SetMaximum(300); gRms2LowPin->Draw("AP"); s="Pin RMS Sqd vs LED (Low Gain"+sConstantBit; gRms2LowPin->SetTitle(s.c_str()); gRms2LowPin->GetXaxis()->SetTitle("Pulser Box * NUMLEDS + Led"); gRms2LowPin->GetYaxis()->SetTitle("Average Rms2"); gRms2LowPin->GetXaxis()->CenterTitle(); gRms2LowPin->GetYaxis()->CenterTitle(); gRms2LowPin->SetMarkerStyle(3); gRms2LowPin->SetMarkerColor(3); gRms2LowPin->SetMarkerSize(0.2); gRms2LowPin->SetLineColor(30); //print graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"PinRms2VsLed.ps)"; else s=sLowRunNumber+"-"+sRunNumber+"PinRms2VsLed.ps)"; cRms2PinGraph->Print(s.c_str()); //draw the num graph TCanvas *cNumPinGraph=new TCanvas("cNumPinGraph","cNumPinGraph", 0,0,1200,600); cNumPinGraph->SetFillColor(0); cNumPinGraph->cd(); gNumHighPin->Draw("AP"); gNumHighPin->SetMinimum(0); gNumHighPin->SetMaximum(pulses+200); s="Pin hits vs LED (High&Low Gain"+sConstantBit; gNumHighPin->SetTitle(s.c_str()); gNumHighPin->GetXaxis()->SetTitle("Pulser Box * NUMLEDS + Led"); gNumHighPin->GetYaxis()->SetTitle("Number of Pulses"); gNumHighPin->GetXaxis()->CenterTitle(); gNumHighPin->GetYaxis()->CenterTitle(); gNumHighPin->SetMarkerStyle(3); gNumHighPin->SetMarkerColor(2); gNumHighPin->SetMarkerSize(0.2); gNumHighPin->SetLineColor(46); gNumLowPin->Draw("P"); gNumLowPin->SetMinimum(0); gNumLowPin->SetMaximum(pulses+200); gNumLowPin->SetTitle(s.c_str()); gNumLowPin->SetMarkerStyle(3); gNumLowPin->SetMarkerColor(3); gNumLowPin->SetMarkerSize(0.3); gNumLowPin->SetLineColor(30); //print graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"PinNumVsLed.ps("; else s=sLowRunNumber+"-"+sRunNumber+"PinNumVsLed.ps("; cNumPinGraph->Print(s.c_str()); //draw the num high graph cNumPinGraph->Clear(); cNumPinGraph->cd(); //gNumHighPin->SetMaximum(50); gNumHighPin->Draw("AP"); s="Pin hits vs LED (High Gain"+sConstantBit; gNumHighPin->SetTitle(s.c_str()); gNumHighPin->GetXaxis()->SetTitle("Pulser Box * NUMLEDS + Led"); gNumHighPin->GetYaxis()->SetTitle("Number of Pulses"); gNumHighPin->GetXaxis()->CenterTitle(); gNumHighPin->GetYaxis()->CenterTitle(); gNumHighPin->SetMarkerStyle(3); gNumHighPin->SetMarkerColor(2); gNumHighPin->SetMarkerSize(0.2); gNumHighPin->SetLineColor(46); //print graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"PinNumVsLed.ps"; else s=sLowRunNumber+"-"+sRunNumber+"PinNumVsLed.ps"; cNumPinGraph->Print(s.c_str()); //draw the num low graph cNumPinGraph->Clear(); cNumPinGraph->cd(); //gNumLowPin->SetMaximum(50); gNumLowPin->Draw("AP"); s="Pin hits vs LED (Low Gain"+sConstantBit; gNumLowPin->SetTitle(s.c_str()); gNumLowPin->GetXaxis()->SetTitle("Pulser Box * NUMLEDS + Led"); gNumLowPin->GetYaxis()->SetTitle("Number of Pulses"); gNumLowPin->GetXaxis()->CenterTitle(); gNumLowPin->GetYaxis()->CenterTitle(); gNumLowPin->SetMarkerStyle(3); gNumLowPin->SetMarkerColor(3); gNumLowPin->SetMarkerSize(0.2); gNumLowPin->SetLineColor(30); //print graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"PinNumVsLed.ps)"; else s=sLowRunNumber+"-"+sRunNumber+"PinNumVsLed.ps)"; cNumPinGraph->Print(s.c_str()); //draw the Res graph TCanvas *cResPinGraph=new TCanvas("cResPinGraph","cResPinGraph", 0,0,1200,600); cResPinGraph->SetFillColor(0); cResPinGraph->cd(); gResHighPin->Draw("AP"); gResHighPin->SetMinimum(0); //gResHighPin->SetMaximum(pulses+200); s="Pin Resolution vs LED (High&Low Gain"+sConstantBit; gResHighPin->SetTitle(s.c_str()); gResHighPin->GetXaxis()->SetTitle("Pulser Box * NUMLEDS + Led"); gResHighPin->GetYaxis()->SetTitle("Resolution (rms/mean)"); gResHighPin->GetXaxis()->CenterTitle(); gResHighPin->GetYaxis()->CenterTitle(); gResHighPin->SetMarkerStyle(3); gResHighPin->SetMarkerColor(2); gResHighPin->SetMarkerSize(0.2); gResHighPin->SetLineColor(46); gResLowPin->Draw("P"); gResLowPin->SetMinimum(0); gResLowPin->SetMaximum(pulses+200); gResLowPin->SetTitle(s.c_str()); gResLowPin->SetMarkerStyle(3); gResLowPin->SetMarkerColor(3); gResLowPin->SetMarkerSize(0.3); gResLowPin->SetLineColor(30); //print graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"PinResVsLed.ps("; else s=sLowRunNumber+"-"+sRunNumber+"PinResVsLed.ps("; cResPinGraph->Print(s.c_str()); //draw the Res high graph cResPinGraph->Clear(); cResPinGraph->cd(); gResHighPin->SetMaximum(0.2); gResHighPin->Draw("AP"); s="Pin Resolution vs LED (High Gain"+sConstantBit; gResHighPin->SetTitle(s.c_str()); gResHighPin->GetXaxis()->SetTitle("Pulser Box * NUMLEDS + Led"); gResHighPin->GetYaxis()->SetTitle("Resolution (rms/mean)"); gResHighPin->GetXaxis()->CenterTitle(); gResHighPin->GetYaxis()->CenterTitle(); gResHighPin->SetMarkerStyle(3); gResHighPin->SetMarkerColor(2); gResHighPin->SetMarkerSize(0.2); gResHighPin->SetLineColor(46); //print graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"PinResVsLed.ps"; else s=sLowRunNumber+"-"+sRunNumber+"PinResVsLed.ps"; cResPinGraph->Print(s.c_str()); //draw the Res low graph cResPinGraph->Clear(); cResPinGraph->cd(); gResLowPin->SetMaximum(0.2); gResLowPin->Draw("AP"); s="Pin Resolution vs LED (Low Gain"+sConstantBit; gResLowPin->SetTitle(s.c_str()); gResLowPin->GetXaxis()->SetTitle("Pulser Box * NUMLEDS + Led"); gResLowPin->GetYaxis()->SetTitle("Resolution (rms/mean)"); gResLowPin->GetXaxis()->CenterTitle(); gResLowPin->GetYaxis()->CenterTitle(); gResLowPin->SetMarkerStyle(3); gResLowPin->SetMarkerColor(3); gResLowPin->SetMarkerSize(0.2); gResLowPin->SetLineColor(30); //print graph to postscript if (sLowRunNumber==sRunNumber) s=sRunNumber+"PinResVsLed.ps)"; else s=sLowRunNumber+"-"+sRunNumber+"PinResVsLed.ps)"; cResPinGraph->Print(s.c_str()); MSG("LIAnalysis",Msg::kInfo) <<" ** Finished the AdcVsPin method ** "<<endl; }

void LIAnalysis::AdcVsPixel (  Int_t  plane, Int_t  AshtrayMin = -1, Int_t  AshtrayMax = -1, Int_t  led = -1 )

Definition at line 5446 of file LIAnalysis.cxx. References ashtray, chain, channel, chip, correlatedHit, crate, detectorType, farLed, farPulserBox, Form(), InitialiseLoopVariables(), led, lookup, mean, MSG, nearLed, LILookup::NearOrFar(), nearPulserBox, numEntries, pixel, plane, LILookup::Plane2CrateEast(), LILookup::Plane2CrateWest(), pulserBox, readoutType, rms, runNumber, s, SetLoopVariables(), and strip. { MSG("LIAnalysis",Msg::kInfo) <<endl<<" ** Running the AdcVsPixel method... ** "<<endl; //get first event chain->GetEvent(0); Int_t cr=lookup.Plane2CrateEast(pl,detectorType); Int_t crOpp=lookup.Plane2CrateWest(pl,detectorType); MSG("LIAnalysis",Msg::kInfo) <<"Plane "<<pl<<", crates "<<cr<<" and "<<crOpp<<" selected"<<endl; if (ashtrayMax!=-1){ MSG("LIAnalysis",Msg::kInfo) <<"Only selecting correlated hits in ashtrays "<<ashtrayMin <<"-"<<ashtrayMax<<" inclusive"<<endl; } else{ MSG("LIAnalysis",Msg::kInfo) <<"Not using ashtray information"<<endl; } string sPlane=Form("%d",pl); string sAshtrayMin=Form("%d",ashtrayMin); string sAshtrayMax=Form("%d",ashtrayMax); string sConstantBit=""; if (ashtrayMax!=-1) sConstantBit=", A "+sAshtrayMin+"-"+ sAshtrayMax+")"; else sConstantBit=")"; //variable to scale all histos to have the same colour Int_t maxColour=400; TH2F ***hAdcVsPixel=0; hAdcVsPixel=new TH2F**[NUMCRATES]; for (Int_t i=0;i<NUMCRATES;i++){ hAdcVsPixel[i]=new TH2F*[NUMCHIPS]; for (Int_t j=0;j<NUMCHIPS;j++){ string sChip=Form("%d",j); string sCrate=Form("%d",i); s="Adc Vs Pixel (Chip "+sChip+", Pl "+sPlane+", Cr "+sCrate+ sConstantBit; (hAdcVsPixel[i])[j]=new TH2F(s.c_str(),s.c_str(), 17,0,17,15000,0,15000); (hAdcVsPixel[i])[j]->GetXaxis()->SetTitle("Pixel"); (hAdcVsPixel[i])[j]->GetXaxis()->CenterTitle(); (hAdcVsPixel[i])[j]->GetYaxis()->SetTitle("Adc"); (hAdcVsPixel[i])[j]->GetYaxis()->CenterTitle(); (hAdcVsPixel[i])[j]->SetFillColor(0); (hAdcVsPixel[i])[j]->SetBit(TH1::kCanRebin); (hAdcVsPixel[i])[j]->SetMaximum(maxColour); } } this->InitialiseLoopVariables(); for(Int_t entry=0;entry<numEvents;entry++){ this->SetLoopVariables(entry,0); //avoid divide by zero errors and skip irrelevant data if (rms==0 || mean==0 || numEntries==0) continue; //only look at scintillator strips if (readoutType!=ReadoutType::kScintStrip) continue; if (plane==pl && (crate==cr || crate==crOpp)){ //Int_t ashtray=lookup.Led2Ashtray(led,plane); if (ashtray>=ashtrayMin && ashtray<=ashtrayMax){ if (correlatedHit==1){ //select the appropriate detectorType if (detectorType==Detector::kFar){ if (pulserBox==nearPulserBox){ (hAdcVsPixel[crate])[chip]->Fill(pixel,mean,maxColour); } else if (pulserBox==farPulserBox){ (hAdcVsPixel[crate][chip])->Fill(pixel,mean,340); } } //CalDet else if (detectorType==Detector::kCalDet){ if (led==l){ if (lookup.NearOrFar(crate,pulserBox,nearPulserBox, farPulserBox,led,detectorType, plane,runNumber)== LILookup::kNearSide){//near (hAdcVsPixel[crate])[chip]->Fill(pixel,mean,400); } else if (lookup.NearOrFar(crate,pulserBox,nearPulserBox, farPulserBox,led,detectorType, plane,runNumber)== LILookup::kFarSide){//far (hAdcVsPixel[crate][chip])->Fill(pixel,mean,340); } else { MSG("LIAnalysis",Msg::kInfo) <<"Missing stuff: cr="<<crate <<", pb="<<pulserBox <<", led="<<led<<", nL="<<nearLed<<", fL="<<farLed <<", st="<<strip<<", chp="<<chip<<", pix="<<pixel <<", ch="<<channel<<", m="<<mean <<endl; } } } //catch exceptions to detector type else { MSG("LIAnalysis",Msg::kWarning) <<"Detector Type = "<<detectorType<<" not supported yet" <<endl; } MSG("LIAnalysis",Msg::kInfo) <<"correlated hit: cr="<<crate <<", pb="<<pulserBox<<", led="<<led <<", ash="<<ashtray<<", str="<<strip <<", chp="<<chip<<", pix="<<pixel<<", ch="<<channel <<", m="<<mean<<endl; } else{ //crosstalk //select the appropriate detectorType if (detectorType==Detector::kFar){ if (pulserBox==nearPulserBox){ (hAdcVsPixel[crate])[chip]->Fill(pixel,mean,1); } else if (pulserBox==farPulserBox){ (hAdcVsPixel[crate])[chip]->Fill(pixel,mean,130); } } //CalDet else if (detectorType==Detector::kCalDet){ if (led==l){ if (lookup.NearOrFar(crate,pulserBox,nearPulserBox, farPulserBox,led,detectorType, plane,runNumber)== LILookup::kNearSide){//near (hAdcVsPixel[crate])[chip]->Fill(pixel,mean,1); } else if (lookup.NearOrFar(crate,pulserBox,nearPulserBox, farPulserBox,led,detectorType, plane,runNumber)== LILookup::kFarSide){//far (hAdcVsPixel[crate])[chip]->Fill(pixel,mean,130); } else { MSG("LIAnalysis",Msg::kInfo) <<"Missing stuff: cr="<<crate <<", pb="<<pulserBox <<", led="<<led <<", nL="<<nearLed <<", fL="<<farLed <<", st="<<strip <<", chp="<<chip <<", pix="<<pixel <<", ch="<<channel<<", m="<<mean <<endl; } } } else { MSG("LIAnalysis",Msg::kWarning) <<"Detector Type = "<<detectorType<<" not supported yet" <<endl; } MSG("LIAnalysis",Msg::kInfo) <<"xtalk: cr="<<crate <<", pb="<<pulserBox <<", led="<<led <<", ash="<<ashtray <<", str="<<strip <<", chp="<<chip <<", pix="<<pixel <<", ch="<<channel<<", m="<<mean <<endl; } } } }//end of for MSG("LIAnalysis",Msg::kInfo)<<"Finished main loop"<<endl; //get rid of the stats info gStyle->SetOptStat(0); //HAS TO BE on the heap to stay on the canvas!!!!!!!!!!!!! TPaveText textLabel=TPaveText(0.6,0.1,0.8,0.6); textLabel.SetBorderSize(0); textLabel.SetTextSize(0.8); textLabel.SetTextColor(1); textLabel.SetTextFont(2); textLabel.SetFillColor(2); TText *label=new TText(0.15,0.9,"COLOUR KEY:"); TText *label1=new TText (0.15,0.81," Dark red = Expected hit: Light from near side"); TText *label2=new TText (0.15,0.71," Orange = Expected hit: Light from far side"); TText *label3=new TText (0.15,0.55," Dark blue = Xtalk: Light from near side"); TText *label4=new TText (0.15,0.45," Light blue = Xtalk: Light from far side"); TCanvas *cAdcVsPixel=new TCanvas("cAdcVsPixel","Adc Vs Pixel", 0,0,1100,850); cAdcVsPixel->SetFillColor(0); cAdcVsPixel->Divide(3,3); for (Int_t i=0;i<NUMCHIPS;i++){ cAdcVsPixel->cd(i+3+1); (hAdcVsPixel[cr])[i]->Draw("colz"); } // TCanvas *cAdcVsPixelOpp=new TCanvas //("cAdcVsPixelOpp","Adc Vs PixelOpp",0,0,1200,700); //cAdcVsPixelOpp->SetFillColor(0); //cAdcVsPixelOpp->Divide(2,2); for (Int_t i=0;i<NUMCHIPS;i++){ cAdcVsPixel->cd(i+6+1); (hAdcVsPixel[crOpp])[i]->Draw("colz"); } //cAdcVsPixelOpp->cd(4); //label->Draw(); //label1->Draw(); //label2->Draw(); //label3->Draw(); //label4->Draw(); //scr=Form("%d",crOpp); //s="adcVsPixPla"+spl+"Cr"+scr+"A"+sAshMin+"-"+sAshMax+".eps"; //cAdcVsPixelOpp->Print(s.c_str()); //cAdcVsPixelOpp->Update(); cAdcVsPixel->cd(2); label->Draw(); label1->Draw(); label2->Draw(); label3->Draw(); label4->Draw(); string scr=Form("%d",cr); string sCrOpp=Form("%d",crOpp); string spl=Form("%d",pl); string sAshMin=Form("%d",ashtrayMin); string sAshMax=Form("%d",ashtrayMax); if (ashtrayMax!=-1) sConstantBit="A "+sAshtrayMin+"-"+sAshtrayMax; else sConstantBit=""; s="adcVsPixPla"+spl+"Cr"+scr+"-"+sCrOpp+sConstantBit+".eps"; cAdcVsPixel->Print(s.c_str()); MSG("LIAnalysis",Msg::kInfo) <<" ** Finished the AdcVsPixel method ** "<<endl; }

void LIAnalysis::AnalyseChain (   )  [private]

Definition at line 331 of file LIAnalysis.cxx. References calibPoint, InitialiseLoopVariables(), lastLed, led, liEvent, liRunNum, maxCalibPoint, maxLedNum, maxPbNum, MSG, numCalibPoints, numLeds, numLiEvents, numLiEventsL, numLiEventsP, numLiRuns, pulseHeight, pulserBox, pulseWidth, and SetLoopVariables(). Referenced by LIAnalysis(). { this->InitialiseLoopVariables(); Int_t eventsPerLedCounter=0; Int_t lastLed2=-1; Int_t lastPulserBox2=-1; Int_t lastCalibPoint2=-1; Int_t numLiEventsC[1000]; for (Int_t i=0;i<1000;i++){ numLiEventsC[i]=0; } for(Int_t entry=0;entry<numEvents;entry++){ this->SetLoopVariables(entry,0); if (led>maxLedNum){ maxLedNum=led; } if (pulserBox>maxPbNum){ maxPbNum=pulserBox; } if (calibPoint>maxCalibPoint){ maxCalibPoint=calibPoint; } eventsPerLedCounter++; if (led!=lastLed2){ //if an led is flashed only a few times print if (eventsPerLedCounter<100){ MSG("LIAnalysis",Msg::kInfo) <<"For LED="<<lastLed <<", numEntries="<<eventsPerLedCounter <<endl; } eventsPerLedCounter=0; numLiEventsL[led-1]++; } if (pulserBox!=lastPulserBox2){ numLiEventsP[pulserBox]++; MSG("LIAnalysis",Msg::kInfo)<<"pulserBox="<<pulserBox<<endl; } if (calibPoint!=lastCalibPoint2){ numLiEventsC[calibPoint-1]++; MSG("LIAnalysis",Msg::kDebug) <<"calibPoint="<<calibPoint <<", pulseHeight="<<pulseHeight <<", pulseWidth="<<pulseWidth <<", led="<<led <<endl; } lastLed2=led; lastPulserBox2=pulserBox; lastCalibPoint2=calibPoint; } //set important variables to determine sizes of //arrays of histograms etc numLeds=maxLedNum; numLiRuns=liRunNum+1; numLiEvents=liEvent+1; numCalibPoints=maxCalibPoint; MSG("LIAnalysis",Msg::kInfo) <<endl <<" ** File Summary Information ** "<<endl <<"maxLedNum="<<maxLedNum<<endl <<"maxPbNum="<<maxPbNum<<endl <<"maxCalibPoint="<<maxCalibPoint<<endl <<"numLiRuns="<<numLiRuns <<", numLiRuns*maxLedNum="<<numLiRuns*maxLedNum<<endl <<"numLiEvents="<<numLiEvents<<endl; MSG("LIAnalysis",Msg::kInfo) <<"Number of times each led was selected (any pulser box):" <<endl; for (Int_t i=0;i<NUMLEDS;i++){ MSG("LIAnalysis",Msg::kInfo) <<" LED "<<i+1<<" = "<<numLiEventsL[i] <<endl; } MSG("LIAnalysis",Msg::kInfo) <<"Number of times each pulser box was selected:"<<endl; for (Int_t i=0;i<NUMPULSERBOXES;i++){ MSG("LIAnalysis",Msg::kInfo) <<" Pulser Box "<<i<<" = "<<numLiEventsP[i] <<endl; } MSG("LIAnalysis",Msg::kInfo) <<"Number of times each calibration point was selected:" <<endl; for (Int_t i=0;i<maxCalibPoint;i++){ MSG("LIAnalysis",Msg::kInfo) <<" Calibration Point "<<i+1<<" = "<<numLiEventsC[i] <<endl; } MSG("LIAnalysis",Msg::kInfo) <<" ** End of File Summary Information ** "<<endl; }

void LIAnalysis::AppendVect (  std::vector< Double_t > &  trunk, std::vector< Double_t >  appendix )  [private]

Definition at line 559 of file LIAnalysis.cxx. References MSG. Referenced by PmtGain(). { MSG("LIAnalysis",Msg::kDebug) <<" ** Running AppendVect method... **"<<endl; //should this work? //trunk.insert(trunk.end(),appendix.begin(),appendix.end()); for (UInt_t i=0;i<appendix.size();i++){ trunk.push_back(appendix[i]); } MSG("LIAnalysis",Msg::kDebug) <<" ** Finished AppendVect method **"<<endl; }

void LIAnalysis::AshtrayVsPlane (  Int_t  nearPb, Int_t  farPb )

Definition at line 4915 of file LIAnalysis.cxx. References ashtray, chain, correlatedHit, crate, detectorType, farPulserBox, Form(), histname, InitialiseLoopVariables(), lookup, mean, MSG, nearLed, nearPulserBox, numEntries, pixel, plane, pow(), pulserBox, readoutType, rms, runNumber, s, SetLoopVariables(), and LILookup::SetPbPlanes(). { //get rid of the stats info gStyle->SetOptStat(0); MSG("LIAnalysis",Msg::kInfo) <<endl<<" ** Running the AshtrayVsPlane method... **"<<endl; chain->GetEvent(0); Int_t* planeMax=new Int_t[NUMPULSERBOXES]; Int_t* planeMin=new Int_t[NUMPULSERBOXES]; lookup.SetPbPlanes(planeMin,planeMax,detectorType); TH2F *hAshtrayVsPlaneLedB=new TH2F(" "," ", LASTSCINTPLANE-FIRSTSCINTPLANE+2, FIRSTSCINTPLANE,LASTSCINTPLANE+1, 22,0,22); TH2F **hAshtrayVsPlaneLed=0; hAshtrayVsPlaneLed= new TH2F*[NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ sprintf(histname,"Led (Pulser Box %i)",i); hAshtrayVsPlaneLed[i]=new TH2F (histname,histname, planeMax[i]-planeMin[i],planeMin[i],planeMax[i] ,22,0,22); hAshtrayVsPlaneLed[i]->GetXaxis()->SetTitle("Plane"); hAshtrayVsPlaneLed[i]->GetXaxis()->CenterTitle(); hAshtrayVsPlaneLed[i]->GetYaxis()->SetTitle("Num Entries"); hAshtrayVsPlaneLed[i]->GetYaxis()->CenterTitle(); hAshtrayVsPlaneLed[i]->SetFillColor(0); hAshtrayVsPlaneLed[i]->SetBit(TH1::kCanRebin); } TH2F **hAshtrayVsPlaneMean=0; hAshtrayVsPlaneMean= new TH2F*[NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ sprintf(histname,"ADC Value (Near Side, Pulser Box %i)",i); hAshtrayVsPlaneMean[i]=new TH2F (histname,histname, planeMax[i]-planeMin[i],planeMin[i],planeMax[i] ,22,0,22); hAshtrayVsPlaneMean[i]->GetXaxis()->SetTitle("Plane"); hAshtrayVsPlaneMean[i]->GetXaxis()->CenterTitle(); hAshtrayVsPlaneMean[i]->GetYaxis()->SetTitle("Ashtray"); hAshtrayVsPlaneMean[i]->GetYaxis()->CenterTitle(); hAshtrayVsPlaneMean[i]->SetFillColor(0); hAshtrayVsPlaneMean[i]->SetBit(TH1::kCanRebin); } TH2F **hAshtrayVsPlaneMeanF=0; hAshtrayVsPlaneMeanF= new TH2F*[NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ sprintf(histname,"ADC Value (Far Side, Pulser Box %i)",i); hAshtrayVsPlaneMeanF[i]=new TH2F (histname,histname, planeMax[i]-planeMin[i],planeMin[i],planeMax[i] ,22,0,22); hAshtrayVsPlaneMeanF[i]->GetXaxis()->SetTitle("Plane"); hAshtrayVsPlaneMeanF[i]->GetXaxis()->CenterTitle(); hAshtrayVsPlaneMeanF[i]->GetYaxis()->SetTitle("Ashtray"); hAshtrayVsPlaneMeanF[i]->GetYaxis()->CenterTitle(); hAshtrayVsPlaneMeanF[i]->SetFillColor(0); hAshtrayVsPlaneMeanF[i]->SetBit(TH1::kCanRebin); } TH2F **hAshtrayVsPlaneNum=0; hAshtrayVsPlaneNum= new TH2F*[NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ sprintf(histname,"Num (Pulser Box %i)",i); hAshtrayVsPlaneNum[i]=new TH2F (histname,histname, planeMax[i]-planeMin[i],planeMin[i],planeMax[i] ,22,0,22); hAshtrayVsPlaneNum[i]->GetXaxis()->SetTitle("Plane"); hAshtrayVsPlaneNum[i]->GetXaxis()->CenterTitle(); hAshtrayVsPlaneNum[i]->GetYaxis()->SetTitle("Ashtray"); hAshtrayVsPlaneNum[i]->GetYaxis()->CenterTitle(); hAshtrayVsPlaneNum[i]->SetFillColor(0); hAshtrayVsPlaneNum[i]->SetBit(TH1::kCanRebin); } TH2F **hAshtrayVsPlanePix=0; hAshtrayVsPlanePix= new TH2F*[NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ sprintf(histname,"Pix (Pulser Box %i)",i); hAshtrayVsPlanePix[i]=new TH2F (histname,histname, planeMax[i]-planeMin[i],planeMin[i],planeMax[i] ,22,0,22); hAshtrayVsPlanePix[i]->GetXaxis()->SetTitle("Plane"); hAshtrayVsPlanePix[i]->GetXaxis()->CenterTitle(); hAshtrayVsPlanePix[i]->GetYaxis()->SetTitle("Ashtray"); hAshtrayVsPlanePix[i]->GetYaxis()->CenterTitle(); hAshtrayVsPlanePix[i]->SetFillColor(0); hAshtrayVsPlanePix[i]->SetBit(TH1::kCanRebin); } TH2F **hAshtrayVsPlanePixF=0; hAshtrayVsPlanePixF= new TH2F*[NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ sprintf(histname,"PixF (Pulser Box %i)",i); hAshtrayVsPlanePixF[i]=new TH2F (histname,histname, planeMax[i]-planeMin[i],planeMin[i],planeMax[i] ,22,0,22); hAshtrayVsPlanePixF[i]->GetXaxis()->SetTitle("Plane"); hAshtrayVsPlanePixF[i]->GetXaxis()->CenterTitle(); hAshtrayVsPlanePixF[i]->GetYaxis()->SetTitle("Ashtray"); hAshtrayVsPlanePixF[i]->GetYaxis()->CenterTitle(); hAshtrayVsPlanePixF[i]->SetFillColor(0); hAshtrayVsPlanePixF[i]->SetBit(TH1::kCanRebin); } TH2F **hAshtrayVsPlaneEnt=0; hAshtrayVsPlaneEnt= new TH2F*[NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ sprintf(histname,"Ent (Pulser Box %i, Near)",i); hAshtrayVsPlaneEnt[i]=new TH2F (histname,histname, planeMax[i]-planeMin[i],planeMin[i],planeMax[i] ,22,0,22); hAshtrayVsPlaneEnt[i]->GetXaxis()->SetTitle("Plane"); hAshtrayVsPlaneEnt[i]->GetXaxis()->CenterTitle(); hAshtrayVsPlaneEnt[i]->GetYaxis()->SetTitle("Ashtray"); hAshtrayVsPlaneEnt[i]->GetYaxis()->CenterTitle(); hAshtrayVsPlaneEnt[i]->SetFillColor(0); hAshtrayVsPlaneEnt[i]->SetBit(TH1::kCanRebin); } TH2F **hAshtrayVsPlaneEntF=0; hAshtrayVsPlaneEntF= new TH2F*[NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ sprintf(histname,"EntF (Pulser Box %i, Far)",i); hAshtrayVsPlaneEntF[i]=new TH2F (histname,histname, planeMax[i]-planeMin[i],planeMin[i],planeMax[i] ,22,0,22); hAshtrayVsPlaneEntF[i]->GetXaxis()->SetTitle("Plane"); hAshtrayVsPlaneEntF[i]->GetXaxis()->CenterTitle(); hAshtrayVsPlaneEntF[i]->GetYaxis()->SetTitle("Ashtray"); hAshtrayVsPlaneEntF[i]->GetYaxis()->CenterTitle(); hAshtrayVsPlaneEntF[i]->SetFillColor(0); hAshtrayVsPlaneEntF[i]->SetBit(TH1::kCanRebin); } this->InitialiseLoopVariables(); for(Int_t entry=0;entry<numEvents;entry++){ this->SetLoopVariables(entry,0); //avoid divide by zero errors and skip irrelevant data if (rms==0. || mean==0. || numEntries==0) continue; //only look at scintillator strips if (readoutType!=ReadoutType::kScintStrip) continue; if (pulserBox>=FIRSTPULSERBOX && pulserBox<NUMPULSERBOXES && correlatedHit!=1 && ashtray!=-1 && pulserBox!=crate && ashtray<=NUMASHTRAYS && ashtray>=FIRSTASHTRAY && plane>=FIRSTSCINTPLANE && plane<=LASTSCINTPLANE && pulserBox==farPulserBox ){ hAshtrayVsPlaneMeanF[pulserBox]->Fill(plane,ashtray,mean); hAshtrayVsPlaneEntF[pulserBox]->Fill(plane,ashtray,1); hAshtrayVsPlanePixF[pulserBox]->Fill(plane,ashtray,pixel+1); } if (pulserBox>=FIRSTPULSERBOX && pulserBox<NUMPULSERBOXES && correlatedHit!=1 && ashtray!=-1 && pulserBox==crate && ashtray<=NUMASHTRAYS && ashtray>=FIRSTASHTRAY && plane>=FIRSTSCINTPLANE && plane<=LASTSCINTPLANE && pulserBox==nearPulserBox ){ hAshtrayVsPlaneLed[pulserBox]->Fill(plane,ashtray,nearLed); hAshtrayVsPlaneMean[pulserBox]->Fill(plane,ashtray,mean); hAshtrayVsPlaneNum[pulserBox]->Fill(plane,ashtray,numEntries); hAshtrayVsPlanePix[pulserBox]->Fill(plane,ashtray,pixel+1); hAshtrayVsPlaneEnt[pulserBox]->Fill(plane,ashtray,1); MSG("LIAnalysis",Msg::kVerbose) <<"mean="<<mean <<", rms="<<rms <<", npe="<<pow(mean/rms,2) <<", gain="<<mean/pow(mean/rms,2) <<endl; } }//end of for MSG("LIAnalysis",Msg::kInfo)<<"Finished main loop"<<endl; string sRunNumber=Form("%d",runNumber); Int_t maxMean=0; Int_t maxMeanF=0; //weight histograms according to number of entries for (Int_t pb=FIRSTPULSERBOX;pb<=LASTPULSERBOX;pb++){ for (Int_t pl=FIRSTSCINTPLANE;pl<=LASTSCINTPLANE;pl++){ for (Int_t ash=0;ash<=NUMASHTRAYS+2;ash++){ //take out known miscablings if ((ash==16 || ash==17) && pl==331) continue; if ((ash==7 || ash==8) && pl==129) continue; if ((ash==9 || ash==10) && pl==130) continue; if (hAshtrayVsPlaneEnt[pb]->GetBinContent(pl,ash)>0){ hAshtrayVsPlaneMean[pb]->SetBinContent (pl,ash,hAshtrayVsPlaneMean[pb]->GetBinContent(pl,ash)/ hAshtrayVsPlaneEnt[pb]->GetBinContent(pl,ash)); if (hAshtrayVsPlaneMean[pb]->GetBinContent(pl,ash)>maxMean){ maxMean=static_cast<Int_t> (hAshtrayVsPlaneMean[pb]->GetBinContent(pl,ash)); } } if (hAshtrayVsPlaneEntF[pb]->GetBinContent(pl,ash)>0){ hAshtrayVsPlaneMeanF[pb]->SetBinContent (pl,ash,hAshtrayVsPlaneMeanF[pb]->GetBinContent(pl,ash)/ hAshtrayVsPlaneEntF[pb]->GetBinContent(pl,ash)); if (hAshtrayVsPlaneMeanF[pb]->GetBinContent(pl,ash)>maxMean){ maxMeanF=static_cast<Int_t> (hAshtrayVsPlaneMeanF[pb]->GetBinContent(pl,ash)); } } } } } //do this after it has been weighted for (Int_t pb=FIRSTPULSERBOX;pb<=LASTPULSERBOX;pb++){ //set a ceiling if (maxMean>4000) maxMean=4000; if (maxMeanF>3000) maxMeanF=3000; hAshtrayVsPlaneMean[pb]->SetMaximum(maxMean); hAshtrayVsPlaneMeanF[pb]->SetMaximum(maxMeanF); } // TCanvas *cAshtrayVsPlaneLed=new TCanvas // ("cAshtrayVsPlaneLed","AshtrayVsPlane: Led",0,0,1000,800); // cAshtrayVsPlaneLed->SetFillColor(0); // cAshtrayVsPlaneLed->cd(); // hAshtrayVsPlaneLedB->Draw(); // for (Int_t i=0;i<NUMPULSERBOXES;i+=2){ // hAshtrayVsPlaneLed[i]->Draw("sametextcolz"); // } TCanvas *cAshtrayVsPlaneMean=new TCanvas ("cAshtrayVsPlaneMean","AshtrayVsPlane: Mean",0,0,1000,800); cAshtrayVsPlaneMean->SetFillColor(0); cAshtrayVsPlaneMean->cd(); for (Int_t i=0;i<NUMPULSERBOXES;i++){ hAshtrayVsPlaneMean[i]->Draw("colz"); if (i==0){ s=sRunNumber+"AdcNearAshtray.ps("; cAshtrayVsPlaneMean->Print(s.c_str()); } else if(i==NUMPULSERBOXES-1){ s=sRunNumber+"AdcNearAshtray.ps)"; cAshtrayVsPlaneMean->Print(s.c_str()); } else{ s=sRunNumber+"AdcNearAshtray.ps"; cAshtrayVsPlaneMean->Print(s.c_str()); } } //draw it so it stays there cAshtrayVsPlaneMean->Clear(); hAshtrayVsPlaneMean[nearPb]->Draw("colz"); //delete cAshtrayVsPlaneMean; TCanvas *cAshtrayVsPlaneMeanF=new TCanvas ("cAshtrayVsPlaneMeanF","AshtrayVsPlane: MeanF",0,0,1000,800); cAshtrayVsPlaneMeanF->SetFillColor(0); cAshtrayVsPlaneMeanF->cd(); for (Int_t i=0;i<NUMPULSERBOXES;i++){ cAshtrayVsPlaneMeanF->Clear(); hAshtrayVsPlaneMeanF[i]->Draw("colz"); if (i==0){ s=sRunNumber+"AdcFarAshtray.ps("; cAshtrayVsPlaneMeanF->Print(s.c_str()); } else if(i==NUMPULSERBOXES-1){ s=sRunNumber+"AdcFarAshtray.ps)"; cAshtrayVsPlaneMeanF->Print(s.c_str()); } else{ s=sRunNumber+"AdcFarAshtray.ps"; cAshtrayVsPlaneMeanF->Print(s.c_str()); } } //draw it so it stays there cAshtrayVsPlaneMeanF->Clear(); hAshtrayVsPlaneMeanF[farPb]->Draw("colz"); //delete cAshtrayVsPlaneMeanF; // TCanvas *cAshtrayVsPlaneNum=new TCanvas // ("cAshtrayVsPlaneNum","AshtrayVsPlane: Number Entries",0,0,1000,800); // cAshtrayVsPlaneNum->SetFillColor(0); // cAshtrayVsPlaneNum->cd(); // hAshtrayVsPlaneNum[pb]->Draw("textcolz"); // TCanvas *cAshtrayVsPlanePix=new TCanvas // ("cAshtrayVsPlanePix","AshtrayVsPlane: Ashtray",0,0,1000,800); // cAshtrayVsPlanePix->SetFillColor(0); // cAshtrayVsPlanePix->cd(); // hAshtrayVsPlanePix[nearPb]->Draw("textcolz"); // TCanvas *cAshtrayVsPlanePixF=new TCanvas // ("cAshtrayVsPlanePixF","AshtrayVsPlane: Ashtray",0,0,1000,800); // cAshtrayVsPlanePixF->SetFillColor(0); // cAshtrayVsPlanePixF->cd(); // hAshtrayVsPlanePixF[farPb]->Draw("textcolz"); TCanvas *cAshtrayVsPlaneEntEast=new TCanvas ("cAshtrayVsPlaneEntEast","AshtrayVsPlane: East Side (Near)" ,0,0,1200,800); cAshtrayVsPlaneEntEast->SetFillColor(0); cAshtrayVsPlaneEntEast->cd(); hAshtrayVsPlaneLedB->SetTitle("East Side (Near)"); hAshtrayVsPlaneLedB->Draw(); for (Int_t i=0;i<NUMPULSERBOXES;i+=2){ hAshtrayVsPlaneEnt[i]->Draw("samecolz"); } s=sRunNumber+"EastNearAshtray.eps"; cAshtrayVsPlaneEntEast->Print(s.c_str()); TCanvas *cAshtrayVsPlaneEntWest=new TCanvas ("cAshtrayVsPlaneEntWest","AshtrayVsPlane: West Side (Near)" ,0,0,1200,800); cAshtrayVsPlaneEntWest->SetFillColor(0); cAshtrayVsPlaneEntWest->cd(); hAshtrayVsPlaneLedB->SetTitle("West Side (Near)"); hAshtrayVsPlaneLedB->Draw(); for (Int_t i=1;i<NUMPULSERBOXES;i+=2){ hAshtrayVsPlaneEnt[i]->Draw("samecolz"); } s=sRunNumber+"WestNearAshtray.eps"; cAshtrayVsPlaneEntWest->Print(s.c_str()); TCanvas *cAshtrayVsPlaneEntWestF=new TCanvas ("cAshtrayVsPlaneEntWestF","AshtrayVsPlane: West Side (Far)" ,0,0,1200,800); cAshtrayVsPlaneEntWestF->SetFillColor(0); cAshtrayVsPlaneEntWestF->cd(); hAshtrayVsPlaneLedB->SetTitle("West Side (Far)"); hAshtrayVsPlaneLedB->Draw(); for (Int_t i=0;i<NUMPULSERBOXES;i+=2){ hAshtrayVsPlaneEntF[i]->Draw("samecolz"); } s=sRunNumber+"WestFarAshtray.eps"; cAshtrayVsPlaneEntWestF->Print(s.c_str()); TCanvas *cAshtrayVsPlaneEntEastF=new TCanvas ("cAshtrayVsPlaneEntEastF","AshtrayVsPlane: East Side (Far)" ,0,0,1200,800); cAshtrayVsPlaneEntEastF->SetFillColor(0); cAshtrayVsPlaneEntEastF->cd(); hAshtrayVsPlaneLedB->SetTitle("East Side (Far)"); hAshtrayVsPlaneLedB->Draw(); for (Int_t i=1;i<NUMPULSERBOXES;i+=2){ hAshtrayVsPlaneEntF[i]->Draw("samecolz"); } s=sRunNumber+"EastFarAshtray.eps"; cAshtrayVsPlaneEntEastF->Print(s.c_str()); MSG("LIAnalysis",Msg::kInfo) <<" ** Finished the AshtrayVsPlane method ** "<<endl; }

void LIAnalysis::CalcAdcAv (   )

Definition at line 6872 of file LIAnalysis.cxx. References calibPoint, histname, InitialiseLoopVariables(), led, mean, MSG, numEntries, pulserBox, readoutType, rms, and SetLoopVariables(). { MSG("LIAnalysis",Msg::kInfo) <<endl<<" ** Running the CalcAdcAv method... ** "<<endl; const Int_t const_numCalibPoints=numCalibPoints; MSG("LIAnalysis",Msg::kInfo) <<"Using "<<const_numCalibPoints<<" calibPoints"<<endl; MSG("LIAnalysis",Msg::kInfo) <<"Number of pulser boxes = "<<NUMPULSERBOXES<<endl; //store the pulse heights associated with each calibPoint Int_t *pbPulseHeight=new Int_t[const_numCalibPoints]; //an array of arrays to hold adc values for different pulse //heights for the different pulser boxes Float_t **pbMean=0; Float_t **pbNumEntries=0; pbMean=new Float_t*[const_numCalibPoints]; pbNumEntries=new Float_t*[const_numCalibPoints]; for (int i=0; i<const_numCalibPoints; i++){ pbMean[i]=new Float_t[NUMPULSERBOXES]; pbNumEntries[i]=new Float_t[NUMPULSERBOXES]; for (int j=0; j<NUMPULSERBOXES; j++){ //initialise (pbMean[i])[j]=0.; (pbNumEntries[i])[j]=0.; } pbPulseHeight[i]=0; } //histograms of the adc values for different pulse heights TH1F **hAdc=0; hAdc= new TH1F*[const_numCalibPoints]; for (Int_t i=0;i<const_numCalibPoints;i++){ sprintf(histname,"Adc Distribution, calibpoint %d",i+1); hAdc[i]=new TH1F(histname,histname,200,0,15000); hAdc[i]->GetXaxis()->SetTitle("Average Adc"); hAdc[i]->GetXaxis()->CenterTitle(); hAdc[i]->GetYaxis()->SetTitle("Num Entries"); hAdc[i]->GetYaxis()->CenterTitle(); hAdc[i]->SetFillColor(0); //hAdc[i]->SetBit(TH1::kCanRebin); } //histograms of the adc values for different pulse heights TH1F **hAdcWithCut=0; hAdcWithCut= new TH1F*[const_numCalibPoints]; for (Int_t i=0;i<const_numCalibPoints;i++){ sprintf(histname,"Adc Distribution, calib point %i",i+1); hAdcWithCut[i]=new TH1F(histname,histname,200,0,15000); hAdcWithCut[i]->GetXaxis()->SetTitle("Average Adc"); hAdcWithCut[i]->GetXaxis()->CenterTitle(); hAdcWithCut[i]->GetYaxis()->SetTitle("Num Entries"); hAdcWithCut[i]->GetYaxis()->CenterTitle(); hAdcWithCut[i]->SetFillColor(0); hAdcWithCut[i]->SetLineColor(2); //hAdcWithCut[i]->SetBit(TH1::kCanRebin); } //histograms of the adc values for different pulse heights TH1F ***hAdcL=0; hAdcL= new TH1F**[const_numCalibPoints]; for (Int_t i=0;i<const_numCalibPoints;i++){ hAdcL[i]= new TH1F*[NUMLEDS]; for (Int_t j=0;j<NUMLEDS;j++){ sprintf(histname,"Adc Distribution (LED), calib point %i",i+1); (hAdcL[i])[j]=new TH1F(histname,histname,60,0,15000); (hAdcL[i])[j]->GetXaxis()->SetTitle("Average Adc"); (hAdcL[i])[j]->GetXaxis()->CenterTitle(); (hAdcL[i])[j]->GetYaxis()->SetTitle("Num Entries"); (hAdcL[i])[j]->GetYaxis()->CenterTitle(); (hAdcL[i])[j]->SetFillColor(0); (hAdcL[i])[j]->SetMarkerSize(20); //(hAdc[i])[j]->SetBit(TH1::kCanRebin); } } this->InitialiseLoopVariables(); for(Int_t entry=0;entry<numEvents;entry++){ this->SetLoopVariables(entry,0); //avoid divide by zero errors and skip irrelevant data if (rms==0. || mean==0. || numEntries==0) continue; //only look at scintillator strips if (readoutType!=ReadoutType::kScintStrip) continue; if (pulserBox>=0 && pulserBox<NUMPULSERBOXES && calibPoint>0 && calibPoint<=const_numCalibPoints && led>0 && led<=NUMLEDS){ if (pulserBox==nearPulserBox){ (pbMean[calibPoint-1])[pulserBox]+=mean; (pbNumEntries[calibPoint-1])[pulserBox]++; hAdcWithCut[calibPoint-1]->Fill(mean); (hAdcL[calibPoint-1])[led-1]->Fill(mean); } hAdc[calibPoint-1]->Fill(mean); pbPulseHeight[calibPoint-1]=pulseHeight; } }//end of for MSG("LIAnalysis",Msg::kInfo)<<"Finished main loop"<<endl; MSG("LIAnalysis",Msg::kInfo) <<endl <<"Average Adcs for each pulser box:"<<endl; for (Int_t calPoint=0;calPoint<const_numCalibPoints;calPoint++){ MSG("LIAnalysis",Msg::kInfo) <<" With pulse height "<<pbPulseHeight[calPoint]<<endl; for (Int_t pb=0; pb<NUMPULSERBOXES; pb++){ //avoid divide by zero Float_t averageMean=0.; if ((pbNumEntries[calPoint])[pb]>0){ //calculte average averageMean=(pbMean[calPoint])[pb]/ (pbNumEntries[calPoint])[pb]; } MSG("LIAnalysis",Msg::kInfo) <<" Average ADC = "<<averageMean <<" for pulser box "<<pb <<" (Num Entries = "<<(pbNumEntries[calPoint])[pb]<<")" <<endl; } } //plot adc distributions for different pulse heights TCanvas *cAdc=new TCanvas("cAdc","ADC Distribution",0,0,800,630); cAdc->SetFillColor(0); cAdc->Divide(2,2); TCanvas *cAdc2=new TCanvas("cAdc2","ADC Distribution",0,0,800,630); cAdc2->SetFillColor(0); cAdc2->Divide(2,2); TCanvas *cAdc3=new TCanvas("cAdc3","ADC Distribution",0,0,800,630); cAdc3->SetFillColor(0); cAdc3->Divide(2,2); TLegend *lAdc= new TLegend(0.85, 0.7, 0.9, 0.9); lAdc->SetBorderSize(0); lAdc->SetFillColor(0); lAdc->SetTextSize(0.035); lAdc->AddEntry(hAdc[0]," All hits","l"); lAdc->AddEntry(hAdcWithCut[0]," Hits from near PBox","l"); for (Int_t i=0;i<const_numCalibPoints;i++){ sprintf(histname,"Adc Distribution, pulse Height %i", pbPulseHeight[i]); hAdc[i]->SetTitle(histname); hAdcWithCut[i]->SetTitle(histname); if (i<4){ cAdc->cd(i+1); hAdc[i]->Draw(); hAdcWithCut[i]->Draw("same"); lAdc->Draw(); } if (i>=4 && i<8){ cAdc2->cd(i-3); hAdc[i]->Draw(); hAdcWithCut[i]->Draw("same"); lAdc->Draw(); } if (i>=8 && i<12){ cAdc3->cd(i-7); hAdc[i]->Draw(); hAdcWithCut[i]->Draw("same"); lAdc->Draw(); } } //plot every LED TCanvas *cAdcL=new TCanvas("cAdcL","ADC Distribution",0,0,800,630); cAdcL->SetFillColor(0); cAdcL->Divide(2,2); TLegend **lAdcL=0; lAdcL= new TLegend*[5]; for (Int_t i=0;i<5;i++){ lAdcL[i]= new TLegend(0.85, 0.7, 0.9, 0.9); lAdcL[i]->SetBorderSize(0); lAdcL[i]->SetFillColor(0); lAdcL[i]->SetTextSize(0.035); } for (Int_t i=0;i<const_numCalibPoints;i++){ if (i<4){ cAdcL->cd(i+1); hAdcWithCut[7]->Draw(); for (Int_t j=i*5;j<(i+1)*5;j++){ (hAdcL[7])[j]->SetLineColor(2+j-(i*5)+1); (hAdcL[7])[j]->Draw("same"); sprintf(histname," LED %i",j+1); lAdcL[i]->AddEntry((hAdcL[7])[j],histname,"l"); } lAdcL[i]->Draw(); } } MSG("LIAnalysis",Msg::kInfo) <<endl<<" ** Finished the CalcAdcAv method ** "<<endl; }

void LIAnalysis::ChannelAdc (  Int_t  chan, Int_t  chan2 )

Definition at line 7090 of file LIAnalysis.cxx. References correlatedHit, crate, InitialiseLoopVariables(), mean, minderCh, MSG, and SetLoopVariables(). { MSG("LIAnalysis",Msg::kInfo) <<endl<<" ** Running the ChannelAdc method... ** "<<endl; gStyle->SetOptStat(1111111); MSG("LIAnalysis",Msg::kInfo) <<"Looking at channels "<<chan<<" and "<<chan2<<endl; TH1F *hMean=new TH1F("hMean","Mean",400,-5,15005); hMean->GetXaxis()->SetTitle("Mean"); hMean->GetXaxis()->CenterTitle(); hMean->GetYaxis()->SetTitle(""); hMean->GetYaxis()->CenterTitle(); hMean->SetFillColor(0); hMean->SetBit(TH1::kCanRebin); TH1F *hMeanCorrelated=new TH1F("hMeanCorrelated","MeanCorrelated", 400,-5,15005); hMeanCorrelated->GetXaxis()->SetTitle("Mean"); hMeanCorrelated->GetXaxis()->CenterTitle(); hMeanCorrelated->GetYaxis()->SetTitle(""); hMeanCorrelated->GetYaxis()->CenterTitle(); hMeanCorrelated->SetFillColor(0); hMeanCorrelated->SetLineColor(2); hMeanCorrelated->SetBit(TH1::kCanRebin); TH1F *hMeanCh=new TH1F("hMeanCh","Mean",400,-5,15005); hMeanCh->GetXaxis()->SetTitle("Mean"); hMeanCh->GetXaxis()->CenterTitle(); hMeanCh->GetYaxis()->SetTitle(""); hMeanCh->GetYaxis()->CenterTitle(); hMeanCh->SetFillColor(0); hMeanCh->SetBit(TH1::kCanRebin); TH1F *hMeanChCorrelated=new TH1F("hMeanChCorrelated","MeanCorrelated", 400,-5,15005); hMeanChCorrelated->GetXaxis()->SetTitle("Mean"); hMeanChCorrelated->GetXaxis()->CenterTitle(); hMeanChCorrelated->GetYaxis()->SetTitle(""); hMeanChCorrelated->GetYaxis()->CenterTitle(); hMeanChCorrelated->SetFillColor(0); hMeanChCorrelated->SetLineColor(2); hMeanChCorrelated->SetBit(TH1::kCanRebin); TH1F *hMeanCh2=new TH1F("hMeanCh2","Mean",400,-5,15005); hMeanCh2->GetXaxis()->SetTitle("Mean"); hMeanCh2->GetXaxis()->CenterTitle(); hMeanCh2->GetYaxis()->SetTitle(""); hMeanCh2->GetYaxis()->CenterTitle(); hMeanCh2->SetFillColor(0); hMeanCh2->SetBit(TH1::kCanRebin); TH1F *hMeanCh2Correlated=new TH1F("hMeanCh2Correlated","MeanCorrelated", 400,-5,15005); hMeanCh2Correlated->GetXaxis()->SetTitle("Mean"); hMeanCh2Correlated->GetXaxis()->CenterTitle(); hMeanCh2Correlated->GetYaxis()->SetTitle(""); hMeanCh2Correlated->GetYaxis()->CenterTitle(); hMeanCh2Correlated->SetFillColor(0); hMeanCh2Correlated->SetLineColor(2); hMeanCh2Correlated->SetBit(TH1::kCanRebin); this->InitialiseLoopVariables(); for(Int_t entry=0;entry<numEvents;entry++){ this->SetLoopVariables(entry); if (crate==0) continue; if (correlatedHit==0){ hMean->Fill(mean); if (minderCh==chan) hMeanCh->Fill(mean); if (minderCh==chan2) hMeanCh2->Fill(mean); } else if (correlatedHit==1){ hMeanCorrelated->Fill(mean); if (minderCh==chan) hMeanChCorrelated->Fill(mean); if (minderCh==chan2) hMeanCh2Correlated->Fill(mean); } }//end of for MSG("LIAnalysis",Msg::kInfo)<<"Finished main loop"<<endl; TCanvas *cMeans=new TCanvas("cMeans","cMeans",0,0,1200,800); cMeans->SetFillColor(0); cMeans->Divide(2,2); cMeans->cd(1); hMean->Draw(); cMeans->cd(2); hMean->Draw(); hMeanCorrelated->Draw("sames"); cMeans->cd(3); hMeanCh->Draw(); hMeanChCorrelated->Draw("sames"); cMeans->cd(4); hMeanCh2->Draw(); hMeanCh2Correlated->Draw("sames"); MSG("LIAnalysis",Msg::kInfo) <<"Looked at channels "<<chan<<" and "<<chan2<<endl; MSG("LIAnalysis",Msg::kInfo) <<endl<<" ** Finished the ChannelAdc method ** "<<endl; }

void LIAnalysis::CheckData (  Int_t  rt = -1  )

Definition at line 7212 of file LIAnalysis.cxx. References ashtray, ReadoutType::AsString(), calibPoint, calibType, chAdd, channel, chip, correlatedHit, crate, eastWest, elecType, farLed, farPulserBox, geoAdd, InitialiseLoopVariables(), inRack, led, masterCh, mean, minderCh, MSG, nearLed, nearPulserBox, numEntries, numericMuxBox, period, pinGain, pinInBox, pixel, plane, print(), PrintBigMessage(), pulseHeight, pulserBox, pulses, pulseWidth, rackBay, rackLevel, readoutType, rms, runNumber, runNumberSub, runType, SetLoopVariables(), strip, stripEnd, summaryCounter, timestamp, varc, vfb, and vmm. { MSG("LIAnalysis",Msg::kInfo) <<endl<<" ** Running the CheckData method... ** "<<endl; gStyle->SetOptStat(1111111); TH1F *hAshtray=new TH1F("hAshtray", "Ashtray distribution",24,-2,22); hAshtray->GetXaxis()->SetTitle("Ashtray"); hAshtray->GetXaxis()->CenterTitle(); hAshtray->GetYaxis()->SetTitle(""); hAshtray->GetYaxis()->CenterTitle(); hAshtray->SetFillColor(0); hAshtray->SetBit(TH1::kCanRebin); TH1F *hCalibPoint=new TH1F("hCalibPoint", "CalibPoint distribution",40,-5,35); hCalibPoint->GetXaxis()->SetTitle("CalibPoint"); hCalibPoint->GetXaxis()->CenterTitle(); hCalibPoint->GetYaxis()->SetTitle(""); hCalibPoint->GetYaxis()->CenterTitle(); hCalibPoint->SetFillColor(0); hCalibPoint->SetBit(TH1::kCanRebin); TH1F *hCalibType=new TH1F("hCalibType","CalibType",100,-5,95); hCalibType->GetXaxis()->SetTitle("CalibType"); hCalibType->GetXaxis()->CenterTitle(); hCalibType->GetYaxis()->SetTitle(""); hCalibType->GetYaxis()->CenterTitle(); hCalibType->SetFillColor(0); hCalibType->SetBit(TH1::kCanRebin); TH1F *hChAdd=new TH1F("hChAdd","ChAdd",8210,-205,8005); hChAdd->GetXaxis()->SetTitle("ChAdd"); hChAdd->GetXaxis()->CenterTitle(); hChAdd->GetYaxis()->SetTitle("Number of times hit"); hChAdd->GetYaxis()->CenterTitle(); hChAdd->SetFillColor(0); hChAdd->SetBit(TH1::kCanRebin); TH1F *hChannel=new TH1F("hChannel","VA Channels hit",27,-2,25); hChannel->GetXaxis()->SetTitle("VA Channel"); hChannel->GetXaxis()->CenterTitle(); hChannel->GetYaxis()->SetTitle("Number of times hit"); hChannel->GetYaxis()->CenterTitle(); hChannel->SetFillColor(0); hChannel->SetBit(TH1::kCanRebin); TH1F *hChip=new TH1F("hChip","Chips hit",7,-2,5); hChip->GetXaxis()->SetTitle("Chip"); hChip->GetXaxis()->CenterTitle(); hChip->GetYaxis()->SetTitle("Number of times hit"); hChip->GetYaxis()->CenterTitle(); hChip->SetFillColor(0); TH1F *hCorrelatedHit=new TH1F("hCorrelatedHit", "Correlated Hit",10,-2,8); hCorrelatedHit->GetXaxis()->SetTitle("CorrelatedHit"); hCorrelatedHit->GetXaxis()->CenterTitle(); hCorrelatedHit->GetYaxis()->SetTitle("num entries"); hCorrelatedHit->GetYaxis()->CenterTitle(); hCorrelatedHit->SetFillColor(0); hCorrelatedHit->SetBit(TH1::kCanRebin); TH1F *hCrate=new TH1F("hCrate","Crates hit",22,-2,20); hCrate->GetXaxis()->SetTitle("Crate"); hCrate->GetXaxis()->CenterTitle(); hCrate->GetYaxis()->SetTitle("Number of times hit"); hCrate->GetYaxis()->CenterTitle(); hCrate->SetFillColor(0); TH1F *hEastWest=new TH1F("hEastWest","EastWest",6,-2,4); hEastWest->GetXaxis()->SetTitle("EastWest"); hEastWest->GetXaxis()->CenterTitle(); hEastWest->GetYaxis()->SetTitle(""); hEastWest->GetYaxis()->CenterTitle(); hEastWest->SetFillColor(0); hEastWest->SetBit(TH1::kCanRebin); TH1F *hEastWestCor=new TH1F("hEastWestCor","EastWestCor",6,-2,4); hEastWestCor->GetXaxis()->SetTitle("EastWest"); hEastWestCor->GetXaxis()->CenterTitle(); hEastWestCor->GetYaxis()->SetTitle(""); hEastWestCor->GetYaxis()->CenterTitle(); hEastWestCor->SetFillColor(0); hEastWestCor->SetLineColor(2); hEastWestCor->SetBit(TH1::kCanRebin); TH1F *hElecType=new TH1F("hElecType","ElecType",12,-2,10); hElecType->GetXaxis()->SetTitle("ElecType"); hElecType->GetXaxis()->CenterTitle(); hElecType->GetYaxis()->SetTitle(""); hElecType->GetYaxis()->CenterTitle(); hElecType->SetFillColor(0); hElecType->SetLineColor(2); hElecType->SetBit(TH1::kCanRebin); TH1F *hFarLed=new TH1F("hFarLed","FarLed",30,-5,25); hFarLed->GetXaxis()->SetTitle("FarLed"); hFarLed->GetXaxis()->CenterTitle(); hFarLed->GetYaxis()->SetTitle("Num Entries"); hFarLed->GetYaxis()->CenterTitle(); hFarLed->SetFillColor(0); hFarLed->SetBit(TH1::kCanRebin); TH1F *hFarPulserBox=new TH1F("hFarPulserBox","FarPulserBox", 30,-5,25); hFarPulserBox->GetXaxis()->SetTitle("FarPulserBox"); hFarPulserBox->GetXaxis()->CenterTitle(); hFarPulserBox->GetYaxis()->SetTitle("Num Entries"); hFarPulserBox->GetYaxis()->CenterTitle(); hFarPulserBox->SetFillColor(0); hFarPulserBox->SetBit(TH1::kCanRebin); TH1F *hGeoAdd=new TH1F("hGeoAdd","GeoAdd", 30,-5,25); hGeoAdd->GetXaxis()->SetTitle("GeoAdd"); hGeoAdd->GetXaxis()->CenterTitle(); hGeoAdd->GetYaxis()->SetTitle("Num Entries"); hGeoAdd->GetYaxis()->CenterTitle(); hGeoAdd->SetFillColor(0); hGeoAdd->SetBit(TH1::kCanRebin); TH1F *hInRack=new TH1F("hInRack","InRack",20,-5,15); hInRack->GetXaxis()->SetTitle("InRack"); hInRack->GetXaxis()->CenterTitle(); hInRack->GetYaxis()->SetTitle(""); hInRack->GetYaxis()->CenterTitle(); hInRack->SetFillColor(0); hInRack->SetBit(TH1::kCanRebin); TH1F *hInRackCor=new TH1F("hInRackCor","InRackCor",20,-5,15); hInRackCor->GetXaxis()->SetTitle("InRack"); hInRackCor->GetXaxis()->CenterTitle(); hInRackCor->GetYaxis()->SetTitle(""); hInRackCor->GetYaxis()->CenterTitle(); hInRackCor->SetFillColor(0); hInRackCor->SetLineColor(2); hInRackCor->SetBit(TH1::kCanRebin); TH1F *hLed=new TH1F("hLed","Leds",30,-5,25); hLed->GetXaxis()->SetTitle("Led"); hLed->GetXaxis()->CenterTitle(); hLed->GetYaxis()->SetTitle(""); hLed->GetYaxis()->CenterTitle(); hLed->SetFillColor(0); hLed->SetBit(TH1::kCanRebin); TH1F *hMasterCh=new TH1F("hMasterCh","MasterCh",30,-5,25); hMasterCh->GetXaxis()->SetTitle("MasterCh"); hMasterCh->GetXaxis()->CenterTitle(); hMasterCh->GetYaxis()->SetTitle(""); hMasterCh->GetYaxis()->CenterTitle(); hMasterCh->SetFillColor(0); hMasterCh->SetBit(TH1::kCanRebin); TH1F *hMean=new TH1F("hMean","Mean",20000,-20,15005); hMean->GetXaxis()->SetTitle("Mean"); hMean->GetXaxis()->CenterTitle(); hMean->GetYaxis()->SetTitle("Num Entries"); hMean->GetYaxis()->CenterTitle(); hMean->SetFillColor(0); hMean->SetBit(TH1::kCanRebin); TH1F *hMeanCorrelated=new TH1F("hMeanCorrelated","MeanCorrelated", 1000,-5,15005); hMeanCorrelated->GetXaxis()->SetTitle("Mean"); hMeanCorrelated->GetXaxis()->CenterTitle(); hMeanCorrelated->GetYaxis()->SetTitle("Num Entries"); hMeanCorrelated->GetYaxis()->CenterTitle(); hMeanCorrelated->SetFillColor(0); hMeanCorrelated->SetLineColor(2); hMeanCorrelated->SetBit(TH1::kCanRebin); TH1F *hMeanFarPB=new TH1F("hMeanFarPB","MeanFarPB", 1000,-5,15005); hMean->GetXaxis()->SetTitle("Mean"); hMeanFarPB->GetXaxis()->CenterTitle(); hMeanFarPB->GetYaxis()->SetTitle("Num Entries"); hMeanFarPB->GetYaxis()->CenterTitle(); hMeanFarPB->SetFillColor(0); hMeanFarPB->SetLineColor(4); hMeanFarPB->SetBit(TH1::kCanRebin); TH1F *hMeanNearPB=new TH1F("hMeanNearPB","MeanNearPB", 1000,-5,15005); hMean->GetXaxis()->SetTitle("Mean"); hMeanNearPB->GetXaxis()->CenterTitle(); hMeanNearPB->GetYaxis()->SetTitle("Num Entries"); hMeanNearPB->GetYaxis()->CenterTitle(); hMeanNearPB->SetFillColor(0); hMeanNearPB->SetLineColor(3); hMeanNearPB->SetBit(TH1::kCanRebin); TH1F *hMinderCh=new TH1F("hMinderCh","MinderCh",30,-5,25); hMinderCh->GetXaxis()->SetTitle("MinderCh"); hMinderCh->GetXaxis()->CenterTitle(); hMinderCh->GetYaxis()->SetTitle(""); hMinderCh->GetYaxis()->CenterTitle(); hMinderCh->SetFillColor(0); hMinderCh->SetBit(TH1::kCanRebin); TH1F *hNearLed=new TH1F("hNearLed","NearLed",30,-5,25); hNearLed->GetXaxis()->SetTitle("NearLed"); hNearLed->GetXaxis()->CenterTitle(); hNearLed->GetYaxis()->SetTitle(""); hNearLed->GetYaxis()->CenterTitle(); hNearLed->SetFillColor(0); hNearLed->SetBit(TH1::kCanRebin); TH1F *hNearPulserBox=new TH1F("hNearPulserBox","NearPulserBox", 30,-5,25); hNearPulserBox->GetXaxis()->SetTitle("NearPulserBox"); hNearPulserBox->GetXaxis()->CenterTitle(); hNearPulserBox->GetYaxis()->SetTitle(""); hNearPulserBox->GetYaxis()->CenterTitle(); hNearPulserBox->SetFillColor(0); hNearPulserBox->SetBit(TH1::kCanRebin); TH1F *hNumEntries=new TH1F("hNumEntries","NumEntries",10000,-5,1000); hNumEntries->GetXaxis()->SetTitle("NumEntries"); hNumEntries->GetXaxis()->CenterTitle(); hNumEntries->GetYaxis()->SetTitle(""); hNumEntries->GetYaxis()->CenterTitle(); hNumEntries->SetFillColor(0); hNumEntries->SetBit(TH1::kCanRebin); TH1F *hNumEntriesCor=new TH1F("hNumEntriesCor","NumEntriesCor", 1000,-5,35); hNumEntriesCor->GetXaxis()->SetTitle("NumEntriesCor"); hNumEntriesCor->GetXaxis()->CenterTitle(); hNumEntriesCor->GetYaxis()->SetTitle(""); hNumEntriesCor->GetYaxis()->CenterTitle(); hNumEntriesCor->SetFillColor(0); hNumEntriesCor->SetLineColor(2); hNumEntriesCor->SetBit(TH1::kCanRebin); TH2F *hNumEntVsPb=new TH2F("hNumEntVsPb", "Number of Pulses Detected Vs Pulser Box (All hits)", NUMPULSERBOXES,FIRSTPULSERBOX, LASTPULSERBOX+1, 10000,-5,1550); hNumEntVsPb->GetXaxis()->SetTitle("Pulser box"); hNumEntVsPb->GetXaxis()->CenterTitle(); hNumEntVsPb->GetYaxis()->SetTitle("Number of Pulses detected"); hNumEntVsPb->GetYaxis()->CenterTitle(); hNumEntVsPb->SetFillColor(0); hNumEntVsPb->SetBit(TH1::kCanRebin); TH1F *hNumericMuxBox=new TH1F("hNumericMuxBox","NumericMuxBox", 6010,-10,6000); hNumericMuxBox->GetXaxis()->SetTitle("NumericMuxBox"); hNumericMuxBox->GetXaxis()->CenterTitle(); hNumericMuxBox->GetYaxis()->SetTitle(""); hNumericMuxBox->GetYaxis()->CenterTitle(); hNumericMuxBox->SetFillColor(0); hNumericMuxBox->SetBit(TH1::kCanRebin); TH1F *hNumericMuxBoxCor=new TH1F("hNumericMuxBoxCor", "NumericMuxBoxCor", 6010,-10,6000); hNumericMuxBoxCor->GetXaxis()->SetTitle("NumericMuxBox"); hNumericMuxBoxCor->GetXaxis()->CenterTitle(); hNumericMuxBoxCor->GetYaxis()->SetTitle(""); hNumericMuxBoxCor->GetYaxis()->CenterTitle(); hNumericMuxBoxCor->SetFillColor(0); hNumericMuxBoxCor->SetLineColor(2); hNumericMuxBoxCor->SetBit(TH1::kCanRebin); TH1F *hPeriod=new TH1F("hPeriod","Period",3010,-5,3005); hPeriod->GetXaxis()->SetTitle("Period"); hPeriod->GetXaxis()->CenterTitle(); hPeriod->GetYaxis()->SetTitle(""); hPeriod->GetYaxis()->CenterTitle(); hPeriod->SetFillColor(0); hPeriod->SetBit(TH1::kCanRebin); TH1F *hPinGain=new TH1F("hPinGain","PinGain",10,-5,5); hPinGain->GetXaxis()->SetTitle("PinGain"); hPinGain->GetXaxis()->CenterTitle(); hPinGain->GetYaxis()->SetTitle(""); hPinGain->GetYaxis()->CenterTitle(); hPinGain->SetFillColor(0); hPinGain->SetBit(TH1::kCanRebin); TH1F *hPinGainCor=new TH1F("hPinGainCor","PinGainCor",10,-5,5); hPinGainCor->GetXaxis()->SetTitle("PinGain"); hPinGainCor->GetXaxis()->CenterTitle(); hPinGainCor->GetYaxis()->SetTitle(""); hPinGainCor->GetYaxis()->CenterTitle(); hPinGainCor->SetFillColor(0); hPinGainCor->SetLineColor(2); hPinGainCor->SetBit(TH1::kCanRebin); TH1F *hPinInBox=new TH1F("hPinInBox","PinInBox",10,-5,5); hPinInBox->GetXaxis()->SetTitle("PinInBox"); hPinInBox->GetXaxis()->CenterTitle(); hPinInBox->GetYaxis()->SetTitle(""); hPinInBox->GetYaxis()->CenterTitle(); hPinInBox->SetFillColor(0); hPinInBox->SetBit(TH1::kCanRebin); TH1F *hPinInBoxCor=new TH1F("hPinInBoxCor","PinInBoxCor",10,-5,5); hPinInBoxCor->GetXaxis()->SetTitle("PinInBox"); hPinInBoxCor->GetXaxis()->CenterTitle(); hPinInBoxCor->GetYaxis()->SetTitle(""); hPinInBoxCor->GetYaxis()->CenterTitle(); hPinInBoxCor->SetFillColor(0); hPinInBoxCor->SetLineColor(2); hPinInBoxCor->SetBit(TH1::kCanRebin); TH1F *hPixel=new TH1F("hPixel","Pixels hit",75,-5,70); hPixel->GetXaxis()->SetTitle("Pixel"); hPixel->GetXaxis()->CenterTitle(); hPixel->GetYaxis()->SetTitle("Number of times hit"); hPixel->GetYaxis()->CenterTitle(); hPixel->SetFillColor(0); hPixel->SetBit(TH1::kCanRebin); TH1F *hPlane=new TH1F("hPlane","Planes hit",805,-5,800); hPlane->GetXaxis()->SetTitle("Plane"); hPlane->GetXaxis()->CenterTitle(); hPlane->GetYaxis()->SetTitle("Number of times hit"); hPlane->GetYaxis()->CenterTitle(); hPlane->SetFillColor(0); TH1F *hPulseHeight=new TH1F("hPulseHeight","PulseHeight", 1055,-5,1050); hPulseHeight->GetXaxis()->SetTitle("PulseHeight"); hPulseHeight->GetXaxis()->CenterTitle(); hPulseHeight->GetYaxis()->SetTitle(""); hPulseHeight->GetYaxis()->CenterTitle(); hPulseHeight->SetFillColor(0); hPulseHeight->SetBit(TH1::kCanRebin); TH1F *hPulserBox=new TH1F("hPulserBox","PulserBox",25,-5,20); hPulserBox->GetXaxis()->SetTitle("PulserBox"); hPulserBox->GetXaxis()->CenterTitle(); hPulserBox->GetYaxis()->SetTitle(""); hPulserBox->GetYaxis()->CenterTitle(); hPulserBox->SetFillColor(0); hPulserBox->SetBit(TH1::kCanRebin); TH1F *hPulses=new TH1F("hPulses","Pulses",11000,0,11000); hPulses->GetXaxis()->SetTitle("Pulses"); hPulses->GetXaxis()->CenterTitle(); hPulses->GetYaxis()->SetTitle(""); hPulses->GetYaxis()->CenterTitle(); hPulses->SetFillColor(0); hPulses->SetBit(TH1::kCanRebin); TH1F *hPulseWidth=new TH1F("hPulseWidth","PulseWidth",10,-2,8); hPulseWidth->GetXaxis()->SetTitle("PulseWidth"); hPulseWidth->GetXaxis()->CenterTitle(); hPulseWidth->GetYaxis()->SetTitle("Num Entries"); hPulseWidth->GetYaxis()->CenterTitle(); hPulseWidth->SetFillColor(0); hPulseWidth->SetBit(TH1::kCanRebin); TH1F *hRackBay=new TH1F("hRackBay","RackBay",22,-2,20); hRackBay->GetXaxis()->SetTitle("rackBay"); hRackBay->GetXaxis()->CenterTitle(); hRackBay->GetYaxis()->SetTitle(""); hRackBay->GetYaxis()->CenterTitle(); hRackBay->SetFillColor(0); hRackBay->SetBit(TH1::kCanRebin); TH1F *hRackBayCor=new TH1F("hRackBayCor","RackBayCor",20,-5,15); hRackBayCor->GetXaxis()->SetTitle("rackBay"); hRackBayCor->GetXaxis()->CenterTitle(); hRackBayCor->GetYaxis()->SetTitle(""); hRackBayCor->GetYaxis()->CenterTitle(); hRackBayCor->SetFillColor(0); hRackBayCor->SetLineColor(2); hRackBayCor->SetBit(TH1::kCanRebin); TH1F *hRackLevel=new TH1F("hRackLevel","RackLevel",5,-2,3); hRackLevel->GetXaxis()->SetTitle("RackLevel"); hRackLevel->GetXaxis()->CenterTitle(); hRackLevel->GetYaxis()->SetTitle(""); hRackLevel->GetYaxis()->CenterTitle(); hRackLevel->SetFillColor(0); hRackLevel->SetBit(TH1::kCanRebin); TH1F *hRackLevelCor=new TH1F("hRackLevelCor","RackLevelCor",5,-2,3); hRackLevelCor->GetXaxis()->SetTitle("RackLevel"); hRackLevelCor->GetXaxis()->CenterTitle(); hRackLevelCor->GetYaxis()->SetTitle(""); hRackLevelCor->GetYaxis()->CenterTitle(); hRackLevelCor->SetFillColor(0); hRackLevelCor->SetLineColor(2); hRackLevelCor->SetBit(TH1::kCanRebin); TH1F *hReadoutType=new TH1F("hReadoutType", "ReadOut Type",100,-2,20); hReadoutType->GetXaxis()->SetTitle("ReadoutType"); hReadoutType->GetXaxis()->CenterTitle(); hReadoutType->GetYaxis()->SetTitle("Num Entries"); hReadoutType->GetYaxis()->CenterTitle(); hReadoutType->SetFillColor(0); hReadoutType->SetBit(TH1::kCanRebin); TH1F *hRms=new TH1F("hRms","Rms",9010,-5,9005); hRms->GetXaxis()->SetTitle("Rms"); hRms->GetXaxis()->CenterTitle(); hRms->GetYaxis()->SetTitle("Entries"); hRms->GetYaxis()->CenterTitle(); hRms->SetFillColor(0); hRms->SetBit(TH1::kCanRebin); TH1F *hRunNumber=new TH1F("hRunNumber","RunNumber", 30010,-5,15000); hRunNumber->GetXaxis()->SetTitle("RunNumber"); hRunNumber->GetXaxis()->CenterTitle(); hRunNumber->GetYaxis()->SetTitle("Entries"); hRunNumber->GetYaxis()->CenterTitle(); hRunNumber->SetFillColor(0); hRunNumber->SetBit(TH1::kCanRebin); TH1F *hRunNumberSub=new TH1F("hRunNumberSub","RunNumberSub", 420,-5,100); hRunNumberSub->GetXaxis()->SetTitle("RunNumberSub"); hRunNumberSub->GetXaxis()->CenterTitle(); hRunNumberSub->GetYaxis()->SetTitle("Entries"); hRunNumberSub->GetYaxis()->CenterTitle(); hRunNumberSub->SetFillColor(0); hRunNumberSub->SetBit(TH1::kCanRebin); TH1F *hRunType=new TH1F("hRunType","RunType", 55,-5,50); hRunType->GetXaxis()->SetTitle("RunType"); hRunType->GetXaxis()->CenterTitle(); hRunType->GetYaxis()->SetTitle("Entries"); hRunType->GetYaxis()->CenterTitle(); hRunType->SetFillColor(0); hRunType->SetBit(TH1::kCanRebin); TH1F *hStrip=new TH1F("hStrip","Strips hit",205,-5,200); hStrip->GetXaxis()->SetTitle("Strip"); hStrip->GetXaxis()->CenterTitle(); hStrip->GetYaxis()->SetTitle("Number of times hit"); hStrip->GetYaxis()->CenterTitle(); hStrip->SetFillColor(0); hStrip->SetBit(TH1::kCanRebin); TH1F *hStripEnd=new TH1F("hStripEnd","StripEnds hit",40,-2,4); hStripEnd->GetXaxis()->SetTitle("StripEnd"); hStripEnd->GetXaxis()->CenterTitle(); hStripEnd->GetYaxis()->SetTitle("Number of times hit"); hStripEnd->GetYaxis()->CenterTitle(); hStripEnd->SetFillColor(0); hStripEnd->SetBit(TH1::kCanRebin); TH1F *hSummaryCounter=new TH1F("hSummaryCounter","SummaryCounter", 1000,-5,100); hSummaryCounter->GetXaxis()->SetTitle("SummaryCounter"); hSummaryCounter->GetXaxis()->CenterTitle(); hSummaryCounter->GetYaxis()->SetTitle("Entries"); hSummaryCounter->GetYaxis()->CenterTitle(); hSummaryCounter->SetFillColor(0); hSummaryCounter->SetBit(TH1::kCanRebin); TH1F *hSummaryCounterCor=new TH1F("hSummaryCounterCor", "SummaryCounterCor", 1010,-5,1005); hSummaryCounterCor->GetXaxis()->SetTitle("SummaryCounter"); hSummaryCounterCor->GetXaxis()->CenterTitle(); hSummaryCounterCor->GetYaxis()->SetTitle("Entries"); hSummaryCounterCor->GetYaxis()->CenterTitle(); hSummaryCounterCor->SetFillColor(0); hSummaryCounterCor->SetLineColor(2); hSummaryCounterCor->SetBit(TH1::kCanRebin); TH1F *hTimestamp=new TH1F("hTimestamp","Timestamp", 1000000,1.0339e9,1.0341e9); hTimestamp->GetXaxis()->SetTitle("Timestamp"); hTimestamp->GetXaxis()->CenterTitle(); hTimestamp->GetYaxis()->SetTitle(""); hTimestamp->GetYaxis()->CenterTitle(); hTimestamp->SetFillColor(0); hTimestamp->SetBit(TH1::kCanRebin); TH1F *hVarc=new TH1F("hVarc","Varcs hit",7,-2,5); hVarc->GetXaxis()->SetTitle("Varc"); hVarc->GetXaxis()->CenterTitle(); hVarc->GetYaxis()->SetTitle("Number of times hit"); hVarc->GetYaxis()->CenterTitle(); hVarc->SetFillColor(0); TH1F *hVfb=new TH1F("hVfb","Vfbs hit",7,-2,5); hVfb->GetXaxis()->SetTitle("Vfb"); hVfb->GetXaxis()->CenterTitle(); hVfb->GetYaxis()->SetTitle("Number of times hit"); hVfb->GetYaxis()->CenterTitle(); hVfb->SetFillColor(0); TH1F *hVmm=new TH1F("hVmm","Vmms hit",10,-2,8); hVmm->GetXaxis()->SetTitle("Vmm"); hVmm->GetXaxis()->CenterTitle(); hVmm->GetYaxis()->SetTitle("Number of times hit"); hVmm->GetYaxis()->CenterTitle(); hVmm->SetFillColor(0); ReadoutType::EReadoutType ert= static_cast<ReadoutType::EReadoutType>(rt); MSG("LIAnalysis",Msg::kInfo) <<endl<<"CheckData for readoutType = " <<ReadoutType::AsString(ert)<<endl<<endl; this->InitialiseLoopVariables(); for(Int_t entry=0;entry<numEvents;entry++){ this->SetLoopVariables(entry,0); //user selection of what to include in histos if (rt==-1);//include everything - this is the default else if (readoutType!=rt) continue; //if (led!=0) continue; //if (crate==0 || crate==2 || geoAdd<=18) continue; if (mean<5) { this->PrintBigMessage(); } if (masterCh>60000) masterCh=-2; if (minderCh>60000) minderCh=-2; if (geoAdd>60000) geoAdd=-2; if (varc>60000) varc=-2; if (vmm>60000) vmm=-2; if (vfb>60000) vfb=-2; if (chip>60000) chip=-2; if (channel>60000) channel=-2; hAshtray->Fill(ashtray); hCalibPoint->Fill(calibPoint); hPulserBox->Fill(pulserBox); hLed->Fill(led); hPulseHeight->Fill(pulseHeight); hPulseWidth->Fill(pulseWidth); hPeriod->Fill(period); hCalibType->Fill(calibType); hTimestamp->Fill(timestamp); hPulses->Fill(pulses); hPlane->Fill(plane); hStrip->Fill(strip); hStripEnd->Fill(stripEnd); hPixel->Fill(pixel); hCrate->Fill(crate); hVarc->Fill(varc); hVmm->Fill(vmm); hVfb->Fill(vfb); hChip->Fill(chip); hChannel->Fill(channel); hMasterCh->Fill(masterCh); hMinderCh->Fill(minderCh); hGeoAdd->Fill(geoAdd); hElecType->Fill(elecType); hChAdd->Fill(chAdd); hNearLed->Fill(nearLed); hFarLed->Fill(farLed); hNearPulserBox->Fill(nearPulserBox); hFarPulserBox->Fill(farPulserBox); hCorrelatedHit->Fill(correlatedHit); hReadoutType->Fill(readoutType); hMean->Fill(mean); if (correlatedHit==1){ hMeanCorrelated->Fill(mean); hSummaryCounterCor->Fill(summaryCounter); hNumEntriesCor->Fill(numEntries); //mux and pin stuff hRackLevelCor->Fill(rackLevel); hEastWestCor->Fill(eastWest); hNumericMuxBoxCor->Fill(numericMuxBox); hInRackCor->Fill(inRack); hRackBayCor->Fill(rackBay); hPinGainCor->Fill(pinGain); hPinInBoxCor->Fill(pinInBox); if (nearPulserBox==pulserBox){ hMeanNearPB->Fill(mean); } else if (farPulserBox==pulserBox){ hMeanFarPB->Fill(mean); } } hRms->Fill(rms); hNumEntries->Fill(numEntries); hNumEntVsPb->Fill(pulserBox,numEntries); hRunNumber->Fill(runNumber); hRunNumberSub->Fill(runNumberSub); hRunType->Fill(runType); hSummaryCounter->Fill(summaryCounter); hRackLevel->Fill(rackLevel); hEastWest->Fill(eastWest); hNumericMuxBox->Fill(numericMuxBox); hInRack->Fill(inRack); hRackBay->Fill(rackBay); hPinGain->Fill(pinGain); hPinInBox->Fill(pinInBox); }//end of for MSG("LIAnalysis",Msg::kInfo)<<"Finished main loop"<<endl; //do plots of all the data in the root file Int_t print=0; TCanvas *cPlane=new TCanvas("cPlane","cPlane",0,0,1200,800); cPlane->SetFillColor(0); cPlane->Divide(2,2); cPlane->cd(1); hPlane->Draw(); cPlane->cd(2); hStrip->Draw(); cPlane->cd(3); hPixel->Draw(); cPlane->cd(4); hStripEnd->Draw(); TCanvas *cElec=new TCanvas("cElec","cElec",0,0,1200,800); cElec->SetFillColor(0); cElec->Divide(1,3); cElec->cd(1); hCrate->Draw(); cElec->cd(2); hVarc->Draw(); cElec->cd(3); hVmm->Draw(); TCanvas *cElec2=new TCanvas("cElec2","cElec2",0,0,1200,800); cElec2->SetFillColor(0); cElec2->Divide(1,3); cElec2->cd(1); hVfb->Draw(); cElec2->cd(2); hChip->Draw(); cElec2->cd(3); hChannel->Draw(); TCanvas *cMean=new TCanvas("cMean","cMean",0,0,1200,800); cMean->SetFillColor(0); cMean->Divide(2,2); cMean->cd(1); hMean->Draw(); hMeanCorrelated->Draw("same"); //hMeanNearPB->Draw("same"); //hMeanFarPB->Draw("same"); cMean->cd(2); hRms->Draw(); cMean->cd(3); hNumEntries->Draw(); hNumEntriesCor->Draw("same"); cMean->cd(4); hSummaryCounter->Draw(); hSummaryCounterCor->Draw("same"); TCanvas *c18=new TCanvas("c18","c18",0,0,1200,800); c18->SetFillColor(0); c18->Divide(1,3); c18->cd(1); hCalibPoint->Draw(); c18->cd(2); hPeriod->Draw(); c18->cd(3); hCalibType->Draw(); if (print){ TCanvas *c19=new TCanvas("c19","c19",0,0,1200,800); c19->SetFillColor(0); c19->Divide(1,2); c19->cd(1); hTimestamp->Draw(); hTimestamp->GetXaxis()->SetTimeFormat("%H:%M-%d/%m"); hTimestamp->GetXaxis()->SetTimeDisplay(1); //hTimestamp->GetXaxis()->SetNdivisions(520); c19->cd(2); hChAdd->Draw(); } TCanvas *c20=new TCanvas("c20","c20",0,0,1200,800); c20->SetFillColor(0); c20->Divide(2,2); c20->cd(1); hPulserBox->Draw(); c20->cd(2); hPulseHeight->Draw(); c20->cd(3); hPulseWidth->Draw(); c20->cd(4); hPulses->Draw(); TCanvas *c21=new TCanvas("c21","c21",0,0,1200,800); c21->SetFillColor(0); c21->Divide(1,2); c21->cd(1); hPinGain->Draw(); hPinGainCor->Draw("same"); c21->cd(2); hPinInBox->Draw(); hPinInBoxCor->Draw("same"); TCanvas *c22=new TCanvas("c22","c22",0,0,1200,800); c22->SetFillColor(0); c22->Divide(3,2); c22->cd(1); hRackLevel->Draw(); hRackLevelCor->Draw("same"); c22->cd(2); hEastWest->Draw(); hEastWestCor->Draw("same"); c22->cd(3); hNumericMuxBox->Draw(); hNumericMuxBoxCor->Draw("same"); c22->cd(4); hInRack->Draw(); hInRackCor->Draw("same"); c22->cd(5); hRackBay->Draw(); hRackBayCor->Draw("same"); TCanvas *c16=new TCanvas("c16","c16",0,0,1200,800); c16->SetFillColor(0); c16->Divide(2,2); c16->cd(1); hLed->Draw(); c16->cd(2); hNearLed->Draw(); c16->cd(3); hFarLed->Draw(); c16->cd(4); hAshtray->Draw(); TCanvas *c23=new TCanvas("c23","c23",0,0,1200,800); c23->SetFillColor(0); c23->Divide(3,3); c23->cd(1); hFarPulserBox->Draw(); c23->cd(2); hNearPulserBox->Draw(); c23->cd(3); hReadoutType->Draw(); c23->cd(4); hRunNumber->Draw(); c23->cd(5); hRunNumberSub->Draw(); c23->cd(6); hRunType->Draw(); c23->cd(7); hCorrelatedHit->Draw(); TCanvas *cElecNear=new TCanvas("cElecNear","cElecNear",0,0,1200,800); cElecNear->SetFillColor(0); cElecNear->Divide(2,2); cElecNear->cd(1); hGeoAdd->Draw(); cElecNear->cd(2); hMasterCh->Draw(); cElecNear->cd(3); hMinderCh->Draw(); cElecNear->cd(4); hElecType->Draw(); if (print==1){ TCanvas *cNumEntVsPb=new TCanvas("cNumEntVsPb","NumEntVsPb", 0,0,1100,600); cNumEntVsPb->cd(); cNumEntVsPb->SetFillColor(0); hNumEntVsPb->Draw("colz"); } MSG("LIAnalysis",Msg::kInfo) <<endl<<" ** Finished the CheckData method ** "<<endl; }

void LIAnalysis::ClearFibres (   )

Definition at line 824 of file LIAnalysis.cxx. References led, MSG, and s. { MSG("LIAnalysis",Msg::kInfo) <<endl<<" ** Running the ClearFibres method... ** "<<endl; //these are the PHs that give you 8000 ADCs minus the //turn on point for each led //rather than the explicit turn on point I chose the //value at which 1000 ADCs because light //level as a function of pulse height would seem //less likely to be linear when the led first turns on //I know that light level is proportional to PH //for most of the PHs from looking at the pin response //(data for PB 11 with the new-new UV leds) Double_t ph[20]={155,175,140,125,145, 100,110,100,120,80, 80,80,60,125,120, 110,110,125,95,135}; Double_t led[20]={1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20}; TCanvas *cPhVsLed=new TCanvas("cPhVsLed","cPhVsLed", 0,0,1200,600); cPhVsLed->SetFillColor(0); cPhVsLed->cd(); TGraph* gPhVsLed=new TGraph(20,led,ph); gPhVsLed->Draw("AP"); s="Corrected Pulse Height Vs Led Number"; gPhVsLed->SetTitle(s.c_str()); gPhVsLed->GetXaxis()->SetTitle("Led Number"); gPhVsLed->GetYaxis()->SetTitle ("Corrected Pulse Height (for 8000 ADCs)"); gPhVsLed->GetXaxis()->CenterTitle(); gPhVsLed->GetYaxis()->CenterTitle(); gPhVsLed->SetMarkerStyle(3); gPhVsLed->SetMarkerColor(2); gPhVsLed->SetMarkerSize(0.5); gPhVsLed->SetMinimum(0); MSG("LIAnalysis",Msg::kInfo) <<" ** Finished the ClearFibres method ** "<<endl; }

void LIAnalysis::DebugPins (   )

Definition at line 1815 of file LIAnalysis.cxx. References chip, correlatedHit, crate, detectorType, eastWest, Form(), geoAdd, GetElecString(), InitialiseLoopVariables(), inRack, MsgService::Instance(), MsgService::IsActive(), led, lowRunNumber, masterCh, mean, MSG, numEntries, pinGain, pinInBox, PrintBigMessage(), pulserBox, pulses, rackBay, rackLevel, readoutType, rms, runNumber, and SetLoopVariables(). { MSG("LIAnalysis",Msg::kInfo) <<endl<<" ** Running the DebugPins method... ** "<<endl; Float_t* maxAdcHighPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* maxAdcLowPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* corrHitHighPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* corrHitLowPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* ledHighPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* nearLedHighPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* ledLowPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* nearLedLowPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* eastWestHighPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* eastWestLowPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* rackLevelHighPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* rackLevelLowPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* rackBayHighPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* rackBayLowPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* inRackHighPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* inRackLowPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* inBoxHighPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* inBoxLowPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* gainHighPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* gainLowPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* crateHighPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* crateLowPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* varcHighPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* varcLowPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* vmmHighPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* vmmLowPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* vfbHighPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* vfbLowPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* chipHighPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* chipLowPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* chAddHighPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* chAddLowPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* planeHighPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; Float_t* planeLowPin=new Float_t[NUMLEDS*NUMPULSERBOXES]; //initialise arrays for (Int_t i=0;i<NUMLEDS*NUMPULSERBOXES;i++){ maxAdcHighPin[i]=0; maxAdcLowPin[i]=0; corrHitHighPin[i]=-1; corrHitLowPin[i]=-1; ledHighPin[i]=0; nearLedHighPin[i]=0; ledLowPin[i]=0; nearLedLowPin[i]=0; eastWestHighPin[i]=0; eastWestLowPin[i]=0; rackLevelHighPin[i]=0; rackLevelLowPin[i]=0; rackBayHighPin[i]=0; rackBayLowPin[i]=0; inRackHighPin[i]=0; inRackLowPin[i]=0; inBoxHighPin[i]=0; inBoxLowPin[i]=0; gainHighPin[i]=0; gainLowPin[i]=0; crateHighPin[i]=0; crateLowPin[i]=0; varcHighPin[i]=0; varcLowPin[i]=0; vmmHighPin[i]=0; vmmLowPin[i]=0; vfbHighPin[i]=0; vfbLowPin[i]=0; chipHighPin[i]=0; chipLowPin[i]=0; chAddHighPin[i]=0; chAddLowPin[i]=0; planeHighPin[i]=0; planeLowPin[i]=0; } this->InitialiseLoopVariables(); for(Int_t entry=0;entry<numEvents;entry++){ this->SetLoopVariables(entry,0); //CalDet ND pins if (crate==1 && geoAdd==19 && masterCh==4 && detectorType==Detector::kCalDet){ this->PrintBigMessage(); } //ignore any zeros that slipped through if (mean==0 || rms==0 || numEntries==0) continue; //only look at scint strips if (readoutType!=ReadoutType::kPinDiode) continue; //printout everything! if (MsgService::Instance()->IsActive("LIAnalysis",Msg::kDebug)){ MSG("LIAnalysis",Msg::kInfo) <<"("<<pulserBox<<":"<<led <<") CH="<<correlatedHit <<" ("<<static_cast<Int_t>(mean) <<","<<static_cast<Int_t>(rms) <<","<<numEntries <<"/"<<pulses <<")," <<" "<<this->GetElecString() <<", mux=("<<eastWest <<","<<rackLevel <<","<<rackBay <<","<<inRack <<","<<pinInBox <<")" <<", gain="<<pinGain <<endl; } //cut out strange pins with high mean if (numEntries<0.8*pulses) continue; if (pulserBox>=FIRSTPULSERBOX && pulserBox<=LASTPULSERBOX && led>=FIRSTLED && led<=LASTLED){ //define led number Int_t l=pulserBox*NUMLEDS+led-1; //will require a correlated hit when plex is working //fill histo for appropriate pin if (chip==1){//high gain (=0) (chip 1) //find max value for high gain pin if (mean>maxAdcHighPin[l]){ maxAdcHighPin[l]=mean; corrHitHighPin[l]=correlatedHit; ledHighPin[l]=led; nearLedHighPin[l]=nearLed; eastWestHighPin[l]=eastWest; rackLevelHighPin[l]=rackLevel; rackBayHighPin[l]=rackBay; inRackHighPin[l]=inRack; inBoxHighPin[l]=pinInBox; gainHighPin[l]=pinGain; crateHighPin[l]=crate; varcHighPin[l]=varc; vmmHighPin[l]=vmm; vfbHighPin[l]=vfb; chipHighPin[l]=chip; chAddHighPin[l]=chAdd; planeHighPin[l]=plane; } } else if (chip==0){//low gain (=1) (chip 0) //find max value for low gain pin if (mean>maxAdcLowPin[l]) { maxAdcLowPin[l]=mean; corrHitLowPin[l]=correlatedHit; ledLowPin[l]=led; nearLedLowPin[l]=nearLed; eastWestLowPin[l]=eastWest; rackLevelLowPin[l]=rackLevel; rackBayLowPin[l]=rackBay; inRackLowPin[l]=inRack; inBoxLowPin[l]=pinInBox; gainLowPin[l]=pinGain; crateLowPin[l]=crate; varcLowPin[l]=varc; vmmLowPin[l]=vmm; vfbLowPin[l]=vfb; chipLowPin[l]=chip; chAddLowPin[l]=chAdd; planeLowPin[l]=plane; } } //print out warnings if there are inconsistencies if ((pinGain!=1 && chip==0) || (pinGain!=0 && chip==1)) { MSG("LIAnalysis",Msg::kVerbose) <<"Strange pin, wrong gain in plex," <<" "<<this->GetElecString() <<", plex gain="<<pinGain<<", mean="<<mean<<", rms="<<rms <<", num="<<numEntries<<endl; } } }//end of for MSG("LIAnalysis",Msg::kInfo)<<"Finished main loop"<<endl; //include the under and overflow counts gStyle->SetOptStat(1111111); //set up useful string string sRunNumber=Form("%d",runNumber); string sLowRunNumber=Form("%d",lowRunNumber); //loop over pins and print out low ones for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ Int_t l=i*NUMLEDS+j; //print out which pins have low pulse height if (maxAdcHighPin[l]<500 && i!=11){ MSG("LIAnalysis",Msg::kInfo) <<"High pin has low max: PB="<<i<<", led="<<j+1 <<", mean="<<maxAdcHighPin[l]<<endl; } if (maxAdcLowPin[l]<200 && i!=11){ MSG("LIAnalysis",Msg::kInfo) <<"Low pin has low max: PB="<<i<<", led="<<j+1 <<", mean="<<maxAdcLowPin[l]<<endl; } } } //loop over pins to compare led and nearLed MSG("LIAnalysis",Msg::kInfo) <<endl<<"Comparing led and nearLed for pins"<<endl; for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ Int_t l=i*NUMLEDS+j; MSG("LIAnalysis",Msg::kInfo) <<"("<<i<<":"<<j+1<<") High gain pin" <<" has nearLed="<<nearLedHighPin[l] <<" (mean="<<maxAdcHighPin[l]<<", plane="<<planeHighPin[l]<<")" <<endl; MSG("LIAnalysis",Msg::kInfo) <<"("<<i<<":"<<j+1<<") High gain pin" <<" has nearLed="<<nearLedLowPin[l] <<" (mean="<<maxAdcLowPin[l]<<", plane="<<planeLowPin[l]<<")" <<endl<<endl; } } TH1F *hVarc=new TH1F("hVarc","Varcs hit",6,-1,5); hVarc->GetXaxis()->SetTitle("Varc"); hVarc->GetXaxis()->CenterTitle(); hVarc->GetYaxis()->SetTitle("Number of times hit"); hVarc->GetYaxis()->CenterTitle(); hVarc->SetFillColor(0); TH1F *hVfb=new TH1F("hVfb","Vfbs hit",6,-1,5); hVfb->GetXaxis()->SetTitle("Vfb"); hVfb->GetXaxis()->CenterTitle(); hVfb->GetYaxis()->SetTitle("Number of times hit"); hVfb->GetYaxis()->CenterTitle(); hVfb->SetFillColor(0); TH1F *hVmm=new TH1F("hVmm","Vmms hit",9,-1,8); hVmm->GetXaxis()->SetTitle("Vmm"); hVmm->GetXaxis()->CenterTitle(); hVmm->GetYaxis()->SetTitle("Number of times hit"); hVmm->GetYaxis()->CenterTitle(); hVmm->SetFillColor(0); TH1F *hVarcW=new TH1F("hVarcW","Varcs hit",6,-1,5); hVarcW->GetXaxis()->SetTitle("Varc"); hVarcW->GetXaxis()->CenterTitle(); hVarcW->GetYaxis()->SetTitle("Number of times hit"); hVarcW->GetYaxis()->CenterTitle(); hVarcW->SetFillColor(0); hVarcW->SetLineColor(2); TH1F *hVfbW=new TH1F("hVfbW","Vfbs hit",6,-1,5); hVfbW->GetXaxis()->SetTitle("Vfb"); hVfbW->GetXaxis()->CenterTitle(); hVfbW->GetYaxis()->SetTitle("Number of times hit"); hVfbW->GetYaxis()->CenterTitle(); hVfbW->SetFillColor(0); hVfbW->SetLineColor(2); TH1F *hVmmW=new TH1F("hVmmW","Vmms hit",9,-1,8); hVmmW->GetXaxis()->SetTitle("Vmm"); hVmmW->GetXaxis()->CenterTitle(); hVmmW->GetYaxis()->SetTitle("Number of times hit"); hVmmW->GetYaxis()->CenterTitle(); hVmmW->SetFillColor(0); hVmmW->SetLineColor(2); MSG("LIAnalysis",Msg::kInfo) <<endl<<"Making pin location histos..."<<endl; for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ Int_t l=i*NUMLEDS+j; Int_t side=i%2; if (side==0){ hVarc->Fill(varcHighPin[l]); hVarc->Fill(varcLowPin[l]); hVmm->Fill(vmmHighPin[l]); hVmm->Fill(vmmLowPin[l]); hVfb->Fill(vfbHighPin[l]); hVfb->Fill(vfbLowPin[l]); } else if (side==1){ hVarcW->Fill(varcHighPin[l]); hVarcW->Fill(varcLowPin[l]); if (varcLowPin[l]==0){ MSG("LIAnalysis",Msg::kInfo) <<"Low pin ("<<i<<":"<<j+1 <<") nL="<<nearLedLowPin[l] <<", corrHit="<<corrHitLowPin[l] <<" VA("<<crateLowPin[l] <<","<<varcLowPin[l] <<","<<vmmLowPin[l] <<","<<vfbLowPin[l] <<","<<chipLowPin[l] <<")="<<chAddLowPin[l] <<", ("<<eastWestLowPin[l] <<","<<rackLevelLowPin[l] <<","<<rackBayLowPin[l] <<","<<inRackLowPin[l] <<","<<inBoxLowPin[l] <<","<<gainLowPin[l] <<")" <<endl<<endl; } if (varcHighPin[l]==0){ MSG("LIAnalysis",Msg::kInfo) <<"High pin ("<<i<<":"<<j+1 <<") nL="<<nearLedHighPin[l] <<", corrHit="<<corrHitHighPin[l] <<" VA("<<crateHighPin[l] <<","<<varcHighPin[l] <<","<<vmmHighPin[l] <<","<<vfbHighPin[l] <<","<<chipHighPin[l] <<")="<<chAddHighPin[l] <<", ("<<eastWestHighPin[l] <<","<<rackLevelHighPin[l] <<","<<rackBayHighPin[l] <<","<<inRackHighPin[l] <<","<<inBoxHighPin[l] <<","<<gainHighPin[l] <<")" <<endl; } hVmmW->Fill(vmmHighPin[l]); hVmmW->Fill(vmmLowPin[l]); hVfbW->Fill(vfbHighPin[l]); hVfbW->Fill(vfbLowPin[l]); } else{ MSG("LIAnalysis",Msg::kInfo) <<"Warning"<<endl; } } } MSG("LIAnalysis",Msg::kInfo) <<endl<<"Comparing mux and elec info..."<<endl; for (Int_t i=0;i<NUMPULSERBOXES;i++){ //leave a line MSG("LIAnalysis",Msg::kInfo)<<endl; for (Int_t j=0;j<NUMLEDS;j++){ //get led index Int_t l=i*NUMLEDS+j; MSG("LIAnalysis",Msg::kInfo) <<"HG ("<<i<<":"<<j+1 <<") nL="<<nearLedHighPin[l] <<", corrHit="<<corrHitHighPin[l] <<", VA("<<crateHighPin[l] <<","<<varcHighPin[l] <<","<<vmmHighPin[l] <<","<<vfbHighPin[l] <<","<<chipHighPin[l] <<")="<<chAddHighPin[l] <<", Position=("<<eastWestHighPin[l] <<","<<rackLevelHighPin[l] <<","<<rackBayHighPin[l] <<","<<inRackHighPin[l] <<","<<inBoxHighPin[l] <<","<<gainHighPin[l] <<")" <<endl; MSG("LIAnalysis",Msg::kInfo) <<"LG ("<<i<<":"<<j+1 <<") nL="<<nearLedLowPin[l] <<", corrHit="<<corrHitLowPin[l] <<", VA("<<crateLowPin[l] <<","<<varcLowPin[l] <<","<<vmmLowPin[l] <<","<<vfbLowPin[l] <<","<<chipLowPin[l] <<")="<<chAddLowPin[l] <<", Position=("<<eastWestLowPin[l] <<","<<rackLevelLowPin[l] <<","<<rackBayLowPin[l] <<","<<inRackLowPin[l] <<","<<inBoxLowPin[l] <<","<<gainLowPin[l] <<")" <<endl<<endl; } } MSG("LIAnalysis",Msg::kInfo) <<endl<<endl<<"No data for these pins:"<<endl; for (Int_t i=0;i<NUMPULSERBOXES;i++){ //leave a line MSG("LIAnalysis",Msg::kInfo)<<endl; for (Int_t j=0;j<NUMLEDS;j++){ //get led index Int_t l=i*NUMLEDS+j; if (corrHitHighPin[l]==-1){ MSG("LIAnalysis",Msg::kInfo) <<"HG ("<<i<<":"<<j+1 <<") nL="<<nearLedHighPin[l] <<", corrHit="<<corrHitHighPin[l] <<", VA("<<crateHighPin[l] <<","<<varcHighPin[l] <<","<<vmmHighPin[l] <<","<<vfbHighPin[l] <<","<<chipHighPin[l] <<")="<<chAddHighPin[l] <<", Position=("<<eastWestHighPin[l] <<","<<rackLevelHighPin[l] <<","<<rackBayHighPin[l] <<","<<inRackHighPin[l] <<","<<inBoxHighPin[l] <<","<<gainHighPin[l] <<")" <<endl; } if (corrHitLowPin[l]==-1){ MSG("LIAnalysis",Msg::kInfo) <<"LG ("<<i<<":"<<j+1 <<") nL="<<nearLedLowPin[l] <<", corrHit="<<corrHitLowPin[l] <<", VA("<<crateLowPin[l] <<","<<varcLowPin[l] <<","<<vmmLowPin[l] <<","<<vfbLowPin[l] <<","<<chipLowPin[l] <<")="<<chAddLowPin[l] <<", Position=("<<eastWestLowPin[l] <<","<<rackLevelLowPin[l] <<","<<rackBayLowPin[l] <<","<<inRackLowPin[l] <<","<<inBoxLowPin[l] <<","<<gainLowPin[l] <<")" <<endl<<endl; } } } //create the canvas and draw TCanvas *cPin=new TCanvas("cAdcPin","cAdcPin",0,0,1000,600); cPin->SetFillColor(0); cPin->Divide(2,2); cPin->cd(1); hVarc->Draw(); hVarcW->Draw("same"); cPin->cd(2); hVmm->Draw(); hVmmW->Draw("same"); cPin->cd(3); hVfb->Draw(); hVfbW->Draw("same"); //The pin diode information in the plex is as follows //The electronics information is in descending hierarchial order //(x:y) is pulser box x and led y //e.g. (crate,varc,vmm,vfb,VA chip,VA channel) //(0:1) Planes=41&43 Pin1=VA(0,2,0,0,0,18) Pin2=VA(0,2,0,0,1,18) //(0:2) Planes=45&47 Pin1=VA(0,2,0,1,0,18) Pin2=VA(0,2,0,1,1,18) //(0:3) Planes=49&51 Pin1=VA(0,2,1,0,0,18) Pin2=VA(0,2,1,0,1,18) //(0:4) Planes=40&42 Pin1=VA(0,2,2,0,0,18) Pin2=VA(0,2,2,1,1,18) //(0:5) Planes=44&46 Pin1=VA(0,2,3,0,0,18) Pin2=VA(0,2,3,0,1,18) //(0:6) Planes=48&50 Pin1=VA(0,2,3,1,0,18) Pin2=VA(0,2,3,1,1,18) MSG("LIAnalysis",Msg::kInfo) <<" ** Finished the DebugPins method ** "<<endl; }

void LIAnalysis::GainCurves (  Int_t  pulserBoxToPlot = 0, Int_t  maxLedToPlot = 6 )

Definition at line 12258 of file LIAnalysis.cxx. References LILookup::CalcGain(), calibPoint, chain, chip, correlatedHit, crate, detectorType, farPulserBox, Form(), GetElecString(), highRunNumber, InitialiseLoopVariables(), MsgService::Instance(), MsgService::IsActive(), led, lookup, lowRunNumber, mean, MSG, LILookup::NearOrFar(), nearPulserBox, numEntries, period, pinGain, plane, print(), LITuning::PrintAll(), pulseHeight, pulserBox, pulses, pulseWidth, readoutType, rms, runNumber, s, SetLoopVariables(), stripEnd, and TGraphVect(). { MSG("LIAnalysis",Msg::kInfo) <<endl<<" ** Running the GainCurves method... ** "<<endl; chain->GetEvent(0); const Int_t const_numCalibPoints=calibType; //just select the first pulser box if not told if (pulserBoxToPlot==-1) pulserBoxToPlot=pulserBox; MSG("LIAnalysis",Msg::kInfo) <<"Will only look at pulser box "<<pulserBoxToPlot<<endl; MSG("LIAnalysis",Msg::kInfo) <<"Using "<<const_numCalibPoints<<" calibPoints"<<endl; MSG("LIAnalysis",Msg::kInfo) <<"Number of pulser boxes = "<<NUMPULSERBOXES <<", leds = "<<NUMLEDS <<endl<<endl; MSG("LIAnalysis",Msg::kInfo) <<"Creating LIRun objects"<<endl; vector<LIRun> gCurveData(NUMSIDES*NUMPLANES*NUMSTRIPS); Float_t **maxAdcHighPin=0; maxAdcHighPin=new Float_t*[NUMLEDS*NUMPULSERBOXES]; Float_t **maxAdcLowPin=0; maxAdcLowPin=new Float_t*[NUMLEDS*NUMPULSERBOXES]; //initialise arrays for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ Int_t l=i*NUMLEDS+j; maxAdcHighPin[l]=new Float_t[const_numCalibPoints]; maxAdcLowPin[l]=new Float_t[const_numCalibPoints]; for (Int_t k=0;k<const_numCalibPoints;k++){ maxAdcHighPin[l][k]=0; maxAdcLowPin[l][k]=0; } } } Bool_t usePinPlex=true; MSG("LIAnalysis",Msg::kInfo) <<endl<<"Looping to find pins"<<endl; this->InitialiseLoopVariables(); for(Int_t entry=0;entry<numEvents;entry++){ this->SetLoopVariables(entry,0); //ignore any zeros that slipped through if (mean==0 || rms==0 || numEntries==0) continue; //only look at one pulser box at a time //especially important at the far detector if (pulserBox!=pulserBoxToPlot) continue; //limit the number of leds that are plotted if (led>maxLedToPlot) continue; //cut out strange pins with high mean if (numEntries<0.5*pulses) continue; Int_t c=calibPoint-1; //look at pins if (readoutType==ReadoutType::kPinDiode){ if (pulserBox>=FIRSTPULSERBOX && pulserBox<=LASTPULSERBOX && led>=FIRSTLED && led<=LASTLED && c<const_numCalibPoints){ //define led number Int_t l=pulserBox*NUMLEDS+led-1; if (usePinPlex){ if (correlatedHit==1){ if (chip==1){//high gain (chip 1) //store high gain pin maxAdcHighPin[l][c]=mean; } else if (chip==0){//low gain (chip 0) //store low gain pin maxAdcLowPin[l][c]=mean; } } } else{ //fill histo for appropriate pin if (chip==1){//high gain (chip 1) //find max value for high gain pin if (mean>maxAdcHighPin[l][c]) { maxAdcHighPin[l][c]=mean; MSG("LIAnalysis",Msg::kDebug) <<"("<<pulserBox<<":"<<led<<") High pin found (" <<", "<<this->GetElecString() <<", plex gain="<<pinGain<<", mean="<<mean <<", rms="<<rms <<", num="<<numEntries<<", cHit="<<correlatedHit<<endl; } } else if (chip==0){//low gain (chip 0) //find max value for low gain pin if (mean>maxAdcLowPin[l][c]) { maxAdcLowPin[l][c]=mean; MSG("LIAnalysis",Msg::kDebug) <<"("<<pulserBox<<":"<<led<<") Low pin found (" <<", "<<this->GetElecString() <<", plex gain="<<pinGain<<", mean="<<mean <<", rms="<<rms <<", num="<<numEntries<<", cHit="<<correlatedHit<<endl; } } } } } }//end of for MSG("LIAnalysis",Msg::kInfo) <<"Finished loop to find pins"<<endl; MSG("LIAnalysis",Msg::kInfo) <<endl<<"Starting main loop"<<endl; this->InitialiseLoopVariables(); for(Int_t entry=0;entry<numEvents;entry++){ this->SetLoopVariables(entry,1); //ignore any zeros that slipped through if (mean==0 || rms==0 || numEntries==0) continue; //only look at scint strips if (readoutType!=ReadoutType::kScintStrip) continue; //only look at one pulser box at a time //especially important at the far detector if (pulserBox!=pulserBoxToPlot) continue; //limit the number of leds that are plotted if (led>maxLedToPlot) continue; //Int_t side=crate%2;//even=0, odd=1 //at CalDet crate 0 is leds 4,5,6 if (correlatedHit==1 && mean<500 && 1.*numEntries/pulses<0.99){ if (lookup.NearOrFar(crate,pulserBox,nearPulserBox, farPulserBox,led,detectorType, plane,runNumber)== LILookup::kNearSide && pulses>0){ MSG("LIAnalysis",Msg::kInfo) <<"mean="<<mean <<" ("<<mean*numEntries/pulses<<")" <<", rms="<<rms <<", num="<<numEntries <<"/"<<pulses <<"="<<1.*numEntries/pulses <<endl; //correct the mean for zeros mean*=(numEntries/pulses); } } //indexes Int_t c=calibPoint-1; Int_t l=pulserBox*NUMLEDS+led-1; Int_t se=(stripEnd-1)*NUMSTRIPS*NUMPLANES+plane*NUMSTRIPS+strip; //create LIRun objects if (correlatedHit==1){ if (detectorType==Detector::kFar){ if (pulserBox>=FIRSTPULSERBOX && pulserBox<=LASTPULSERBOX && led>=FIRSTLED && led<=LASTLED && calibPoint>=1 && calibPoint<=const_numCalibPoints){ //fill histos with when pulser box is near pulser box if (nearPulserBox==pulserBox){ gCurveData[se].SetLIInfo(pulserBox,led,LIRun::kGainCurve); gCurveData[se].AddPoint(pulseHeight,mean, maxAdcHighPin[l][c], maxAdcLowPin[l][c], rms,//a bit of a hack (adcF) numEntries); } //fill histos with when pulser box is far pulser box else if (farPulserBox==pulserBox) { } } } //change the farPb to farLed for calDet else if (detectorType==Detector::kCalDet){ MSG("LIAnalysis",Msg::kVerbose) <<"("<<pulserBox<<":"<<led <<") caldet se="<<se<<endl; if (pulserBox>=FIRSTPULSERBOX && pulserBox<=LASTPULSERBOX && crate>=FIRSTCRATE && crate<=LASTCRATE && led>=FIRSTLED && led<=LASTLED && calibPoint>=1 && calibPoint<=const_numCalibPoints && strip>=FIRSTSTRIP && strip<=LASTSTRIP && plane>=FIRSTPLANE && plane<=LASTPLANE){ //fill object when pulser box is near pulser box if (lookup.NearOrFar(crate,pulserBox,nearPulserBox, farPulserBox,led,detectorType, plane,runNumber)== LILookup::kNearSide){ //add the points to the LIRun object gCurveData[se].SetLIInfo(pulserBox,led,LIRun::kGainCurve); gCurveData[se].AddPoint(pulseHeight,mean, maxAdcHighPin[l][c], maxAdcLowPin[l][c], rms,//a bit of a hack (adcF) numEntries); } //fill object when pulser box is far pulser box else if (lookup.NearOrFar(crate,pulserBox,nearPulserBox, farPulserBox,led,detectorType, plane,runNumber)== LILookup::kFarSide) { //MSG("LIAnalysis",Msg::kInfo) //<<"("<<pulserBox<<":"<<led //<<") filling far se="<<se<<endl; } } } } }//end of for MSG("LIAnalysis",Msg::kInfo) <<"Finished main loop"<<endl; //include the under and overflow counts gStyle->SetOptStat(1111111); //turn off the fit stats box gStyle->SetOptFit(0000); //set up useful string string sRunNumber=Form("%d",runNumber); string sLowRunNumber=Form("%d",lowRunNumber); string sHighRunNumber=Form("%d",highRunNumber); s=""; //set strings for use in titles string sPulseWidth=Form("%d",pulseWidth); string sPulseFreq=Form("%d",static_cast<Int_t> (ceil(1.0/(period*1.0e-5)))); //loop and print all the info stored in the object if (MsgService::Instance()->IsActive("LIAnalysis",Msg::kDebug)){ for (Int_t k=0;k<NUMSIDES;k++){ MSG("LIAnalysis",Msg::kDebug)<<"Side="<<k<<endl; for (Int_t i=0;i<NUMPLANES;i++){ MSG("LIAnalysis",Msg::kDebug)<<"Plane="<<i<<endl; for (Int_t j=0;j<NUMSTRIPS;j++){ Int_t se=k*NUMPLANES*NUMSTRIPS+i*NUMSTRIPS+j; MSG("LIAnalysis",Msg::kVerbose)<<"se="<<se<<endl; gCurveData[se].PrintAll(); } } } } //create canvas TCanvas *cGc=new TCanvas("cGc","cGc",0,0,1200,600); cGc->SetFillColor(0); TCanvas *cResiduals=new TCanvas("cResiduals","cResiduals", 0,0,1200,600); cResiduals->SetFillColor(0); TCanvas *cResidPin=new TCanvas("cResidPin","cResidPin",0,0,1200,600); cResidPin->SetFillColor(0); TCanvas *cNonLin=new TCanvas("cNonLin","cNonLin", 0,0,1200,600); cNonLin->SetFillColor(0); TCanvas *cNonLinPin=new TCanvas("cNonLinPin","cNonLinPin", 0,0,1200,600); cNonLinPin->SetFillColor(0); TCanvas *cGain=new TCanvas("cGain","cGain",0,0,1200,600); cGain->SetFillColor(0); TCanvas *cNumEntries=new TCanvas("cNumEntries","cNumEntries", 0,0,1200,600); cNumEntries->SetFillColor(0); //create the graphs TGraphErrors* gGc=new TGraphErrors(const_numCalibPoints); TGraphErrors* gResiduals=new TGraphErrors(const_numCalibPoints); TGraphErrors* gResidPin=new TGraphErrors(const_numCalibPoints); TGraph* gNonLin=new TGraph(const_numCalibPoints); TGraph* gNonLinPin=new TGraph(const_numCalibPoints); TGraph* gGain=new TGraph(const_numCalibPoints); TGraph* gNumEntries=new TGraph(const_numCalibPoints); vector<TGraph*> vResid(5); vResid[0]=new TGraph(const_numCalibPoints); vResid[1]=new TGraph(const_numCalibPoints); vResid[2]=new TGraph(const_numCalibPoints); vResid[3]=new TGraph(const_numCalibPoints); vResid[4]=new TGraph(const_numCalibPoints); string sResiduals=sRunNumber+"GcResidVsPmt.ps"; string sResidPin=sRunNumber+"GcResidVsPin.ps"; string sGc=sRunNumber+"Gc.ps"; string sNonLin=sRunNumber+"GcNonLinVsPmt.ps"; string sNonLinPin=sRunNumber+"GcNonLinVsPin.ps"; string sGain=sRunNumber+"GcGain.ps"; string sNumEntries=sRunNumber+"GcNumEntries.ps"; Int_t fitPoint=-1; Int_t firstGoodAdcPoint=-1; Int_t nonLinNormPoint=-1; Double_t maxRatio=-1; Double_t pinAdcCut=10; //Int_t firstPlotDone=0; MSG("LIAnalysis",Msg::kInfo) <<"Starting loop over planes..."<<endl; for (Int_t i=0;i<NUMPLANES;i++){ MSG("LIAnalysis",Msg::kDebug)<<"Plane="<<i<<endl; //ignore the pain in the arse two ends in one side on plane zero if (i==0 && detectorType==Detector::kCalDet) continue; //not too many plots either if (i>=10 && detectorType==Detector::kCalDet) continue; //print out occasionally if (i%(NUMPLANES/5)==0){ MSG("LIAnalysis",Msg::kInfo)<<"Plane="<<i<<endl; } for (Int_t j=0;j<NUMSTRIPS;j++){ for (Int_t k=FIRSTEND;k<=SECONDEND;k++){ //calculate the index Int_t se=(k-1)*NUMPLANES*NUMSTRIPS+i*NUMSTRIPS+j; string sPlane=Form("%d",i); string sStrip=Form("%d",j); string sEnd=Form("%d",k); string sPlStEnd=" (Pl,St,End)=("+sPlane+","+sStrip+","+sEnd+")"; MSG("LIAnalysis",Msg::kVerbose)<<sPlStEnd<<endl; vector<Double_t> pin; vector<Double_t> adc; pin=gCurveData[se].GetPin(2); adc=gCurveData[se].GetAdc(); //use built in functionality but with wrong name vector<Double_t> numEntries; numEntries=gCurveData[se].GetAdcLow(); vector<Double_t> rms; rms=gCurveData[se].GetAdcF();//a quick way to store it gGc->Clear(); cGc->Clear(); gResiduals->Clear(); cResiduals->Clear(); gResidPin->Clear(); cResidPin->Clear(); gNonLin->Clear(); cNonLin->Clear(); gNonLinPin->Clear(); cNonLinPin->Clear(); gGain->Clear(); cGain->Clear(); gNumEntries->Clear(); cNumEntries->Clear(); //reset the fit point fitPoint=-1; firstGoodAdcPoint=-1; nonLinNormPoint=-1; maxRatio=-1; //only look at decent gain curves if (adc.size()>0){ //print out occasionally if (i%3==0 && j%6==0){ MSG("LIAnalysis",Msg::kInfo) <<"Processing plane="<<sPlStEnd<<endl; } //find a load of points to help in fitting and plotting for (UInt_t p=0;p<adc.size();p++){ //calcuate error Double_t ey=0; if (numEntries[p]>0){ ey=rms[p]/sqrt(numEntries[p]); } //set the points on the graph gGc->SetPoint(p,pin[p],adc[p]); gGc->SetPointError(p,0.,ey); //find the max point on the graph to which to fit if (adc[p]>5500 && fitPoint==-1){ fitPoint=p; MSG("LIAnalysis",Msg::kDebug) <<sPlStEnd<<" fitPoint="<<fitPoint <<", adc="<<adc[fitPoint] <<", pin="<<pin[fitPoint]<<endl; } //find the first good adc point on the graph to which to fit if (adc[p]>400 && firstGoodAdcPoint==-1){// >5pe to use firstGoodAdcPoint=p; MSG("LIAnalysis",Msg::kDebug) <<sPlStEnd<<" firstGoodAdcPoint="<<firstGoodAdcPoint <<", adc="<<adc[firstGoodAdcPoint] <<", pin="<<pin[firstGoodAdcPoint]<<endl; //set the pin cut to be just below the first good adc pinAdcCut=pin[firstGoodAdcPoint]-1; } Double_t ratio=-1; if (pin[p]!=0){ ratio=adc[p]/pin[p]; } //find the point to normalise to if (ratio>maxRatio && pin[p]>pinAdcCut){ maxRatio=ratio; nonLinNormPoint=p; MSG("LIAnalysis",Msg::kVerbose) <<sPlStEnd <<", nonLinNormPoint="<<nonLinNormPoint <<", adc="<<adc[nonLinNormPoint] <<", pin="<<pin[nonLinNormPoint]<<endl; } } if (firstGoodAdcPoint==fitPoint){ MSG("LIAnalysis",Msg::kWarning) <<sPlStEnd<<" Can't fit GC, firstGoodAdcPoint=" <<firstGoodAdcPoint<<", fitPoint="<<fitPoint<<endl; } //protect against not finding appropriate points if (fitPoint==-1) fitPoint=0; if (firstGoodAdcPoint==-1) firstGoodAdcPoint=0; if (nonLinNormPoint==-1) nonLinNormPoint=0; //fit the gain curve in the appropriate range gGc->Fit("pol1","q","",pinAdcCut,pin[fitPoint]); TF1* func=gGc->GetFunction("pol1"); Double_t funcC=func->GetParameter(0); Double_t funcM=func->GetParameter(1); //draw the gain curve cGc->cd(); gGc->Draw("AP"); s="Gain Curve"+sPlStEnd; gGc->Draw("AP"); gGc->SetTitle(s.c_str()); gGc->GetXaxis()->SetTitle("Pin Adc"); gGc->GetYaxis()->SetTitle("Pmt Adc"); gGc->GetXaxis()->CenterTitle(); gGc->GetYaxis()->CenterTitle(); gGc->SetMarkerStyle(3); gGc->SetMarkerColor(2); gGc->SetMarkerSize(0.55); gGc->SetLineColor(46); //get the axis values Axis_t xMax=gGc->GetXaxis()->GetXmax(); Axis_t xMin=gGc->GetXaxis()->GetXmin(); MSG("LIAnalysis",Msg::kDebug) <<sPlStEnd<<", y="<<funcM<<"x+"<<funcC <<", xMin="<<xMin<<", xMax="<<xMax<<endl; for (Int_t k=0;k<const_numCalibPoints;k++){ Double_t actualX=-1; Double_t actualYval=-1; Double_t ey=0; Double_t fractEy=0; gGc->GetPoint(k,actualX,actualYval); ey=gGc->GetErrorY(k); if (actualYval>0){ fractEy=ey/actualYval; } else{ MSG("LIAnalysis",Msg::kWarning) <<sPlStEnd<<" actualY="<<actualYval<<", ey="<<ey<<endl; } //do y=mx+c to get value on straight line fit Double_t fittedYval=funcM*actualX+funcC; //calculate the residual and normalise Double_t resid=actualYval-fittedYval; if (fittedYval!=0) resid/=fittedYval; else{ resid=0; MSG("LIAnalysis",Msg::kWarning) <<sPlStEnd<<" Fitted y val=0, setting residual to zero" <<endl; } Double_t eResid=fractEy*resid; MSG("LIAnalysis",Msg::kVerbose) <<sPlStEnd<<" actualY="<<actualYval<<", ey="<<ey <<", resid="<<resid<<", eResid="<<eResid <<endl; //set the point on the graph //if (actualX>pinAdcCut){ gResiduals->SetPoint(k,actualYval,resid); gResiduals->SetPointError(k,0,eResid); gResidPin->SetPoint(k,actualX,resid); gResidPin->SetPointError(k,0,eResid); //section for special 5 strips plot Int_t resIndex=se-25; if (resIndex>=0 && resIndex<5){ vResid[resIndex]->SetPoint(k,actualYval,resid); } //} //else { //MSG("LIAnalysis",Msg::kVerbose) //<<"Setting a zero in residuals because pin too low" //<<endl; //gResiduals->SetPoint(k,0,1); //gResidPin->SetPoint(k,0,1); //} if (i<5){ MSG("LIAnalysis",Msg::kDebug) <<sPlStEnd<<", x="<<actualX <<", y="<<actualYval <<", fitY="<<fittedYval<<", resid="<<resid<<i<<endl; } } cResiduals->cd(); s="Gain Curve: Residuals to Fit Vs PMT ADC"+sPlStEnd; gResiduals->Draw("AP"); gResiduals->SetTitle(s.c_str()); gResiduals->GetXaxis()->SetTitle("Pmt Adc"); gResiduals->GetYaxis()->SetTitle("Fractional residual"); gResiduals->GetXaxis()->CenterTitle(); gResiduals->GetYaxis()->CenterTitle(); gResiduals->SetMarkerStyle(3); gResiduals->SetMarkerColor(2); gResiduals->SetMarkerSize(0.55); gResiduals->SetLineColor(46); gResiduals->SetMinimum(-0.2); gResiduals->SetMaximum(0.2); cResidPin->cd(); s="Gain Curve: Residuals to Fit Vs PIN ADC"+sPlStEnd; gResidPin->Draw("AP"); gResidPin->SetTitle(s.c_str()); gResidPin->GetXaxis()->SetTitle("Pin Adc"); gResidPin->GetYaxis()->SetTitle("Fractional residual"); gResidPin->GetXaxis()->CenterTitle(); gResidPin->GetYaxis()->CenterTitle(); gResidPin->SetMarkerStyle(3); gResidPin->SetMarkerColor(2); gResidPin->SetMarkerSize(0.55); gResidPin->SetLineColor(46); gResidPin->SetMinimum(-0.2); gResidPin->SetMaximum(0.2); gNumEntries=TGraphVect(adc,numEntries); cNumEntries->cd(); s="Num Flashes Vs PMT ADC"+sPlStEnd; gNumEntries->Draw("AP"); gNumEntries->SetTitle(s.c_str()); gNumEntries->GetXaxis()->SetTitle("Pmt Adc"); gNumEntries->GetYaxis()->SetTitle("Number of Entries"); gNumEntries->GetXaxis()->CenterTitle(); gNumEntries->GetYaxis()->CenterTitle(); gNumEntries->SetMarkerStyle(3); gNumEntries->SetMarkerColor(2); gNumEntries->SetMarkerSize(0.55); gNumEntries->SetLineColor(46); gNumEntries->SetMaximum(pulses+0.1*pulses); //fill the graphs for (UInt_t p=0;p<adc.size();p++){ //set the points on the graph if (pin[p]!=0 && pin[nonLinNormPoint]!=0 && adc[nonLinNormPoint]!=0){ //calculate variables Double_t ratio=adc[p]/pin[p]; Double_t normRatio=adc[nonLinNormPoint]/ pin[nonLinNormPoint]; //fill the graphs gNonLin->SetPoint(p,adc[p],ratio/normRatio); gNonLinPin->SetPoint(p,pin[p],ratio/normRatio); // PMT fudge factor for M64 (near detector): 0.844 if (detectorType==Detector::kNear) { gGain->SetPoint (p,adc[p],0.844*lookup.CalcGain(adc[p],rms[p],0)); } else { gGain->SetPoint (p,adc[p],0.8*lookup.CalcGain(adc[p],rms[p],0)); } } else { gNonLin->SetPoint(p,0,0); gNonLinPin->SetPoint(p,0,0); gGain->SetPoint(p,0,0); } } cNonLin->cd(); s="Gain Curve: Nonlinearity Vs PMT ADC"+sPlStEnd; gNonLin->Draw("AP"); gNonLin->SetTitle(s.c_str()); gNonLin->GetXaxis()->SetTitle("Pmt Adc"); gNonLin->GetYaxis()->SetTitle ("Fractional Nonlinearity (Norm. Pmt/Pin)"); gNonLin->GetXaxis()->CenterTitle(); gNonLin->GetYaxis()->CenterTitle(); gNonLin->SetMarkerStyle(3); gNonLin->SetMarkerColor(2); gNonLin->SetMarkerSize(0.55); gNonLin->SetLineColor(46); gNonLin->SetMaximum(1.2); cNonLinPin->cd(); s="Gain Curve: Nonlinearity Vs PIN ADC"+sPlStEnd; gNonLinPin->Draw("AP"); gNonLinPin->SetTitle(s.c_str()); gNonLinPin->GetXaxis()->SetTitle("Pin Adc"); gNonLinPin->GetYaxis()->SetTitle ("Fractional Nonlinearity (Norm. Pmt/Pin)"); gNonLinPin->GetXaxis()->CenterTitle(); gNonLinPin->GetYaxis()->CenterTitle(); gNonLinPin->SetMarkerStyle(3); gNonLinPin->SetMarkerColor(2); gNonLinPin->SetMarkerSize(0.55); gNonLinPin->SetLineColor(46); gNonLinPin->SetMaximum(1.2); cGain->cd(); s="Gain (mean/npe) Vs PMT Mean"+sPlStEnd; gGain->Draw("AP"); gGain->SetTitle(s.c_str()); gGain->GetXaxis()->SetTitle("Pmt Adc"); gGain->GetYaxis()->SetTitle("Gain (mean/npe)"); gGain->GetXaxis()->CenterTitle(); gGain->GetYaxis()->CenterTitle(); gGain->SetMarkerStyle(3); gGain->SetMarkerColor(2); gGain->SetMarkerSize(0.55); gGain->SetLineColor(46); gGain->SetMaximum(130); gGain->SetMinimum(0); static Bool_t firstPlotDone=false; if (!firstPlotDone){ //print canvases cResiduals->Print((sResiduals+"(").c_str()); cResidPin->Print((sResidPin+"(").c_str()); cGc->Print((sGc+"(").c_str()); cNonLin->Print((sNonLin+"(").c_str()); cNonLinPin->Print((sNonLinPin+"(").c_str()); cGain->Print((sGain+"(").c_str()); cNumEntries->Print((sNumEntries+"(").c_str()); //set variables gErrorIgnoreLevel=1; firstPlotDone=true; } else { cResiduals->Print(sResiduals.c_str()); cResidPin->Print(sResidPin.c_str()); cGc->Print(sGc.c_str()); cNonLin->Print(sNonLin.c_str()); cNonLinPin->Print(sNonLinPin.c_str()); cGain->Print(sGain.c_str()); cNumEntries->Print(sNumEntries.c_str()); } } } } } //close the files MSG("LIAnalysis",Msg::kInfo)<<"Closing the postscript files..."<<endl; gErrorIgnoreLevel=0; cResiduals->Print((sResiduals+")").c_str()); cResidPin->Print((sResidPin+")").c_str()); cGc->Print((sGc+")").c_str()); cNonLin->Print((sNonLin+")").c_str()); cNonLinPin->Print((sNonLinPin+")").c_str()); cGain->Print((sGain+")").c_str()); cNumEntries->Print((sNumEntries+")").c_str()); MSG("LIAnalysis",Msg::kInfo)<<"... closed"<<endl; //a special plot for 5 on one canvas Int_t print=0; if (print==1){ cResiduals->Clear(); s="Gain Curve Residuals (5 strips)"; vResid[0]->Draw("AP"); vResid[0]->SetTitle(s.c_str()); vResid[0]->GetXaxis()->SetTitle("Pmt Adc"); vResid[0]->GetYaxis()->SetTitle("Normalised residual"); vResid[0]->GetXaxis()->CenterTitle(); vResid[0]->GetYaxis()->CenterTitle(); vResid[0]->SetMarkerStyle(3); vResid[0]->SetMarkerColor(2); vResid[0]->SetMarkerSize(0.55); vResid[0]->SetLineColor(46); vResid[0]->SetMinimum(-0.2); vResid[0]->SetMaximum(0.2); for (Int_t i=1;i<5;i++){ vResid[i]->Draw("P"); vResid[i]->SetTitle(s.c_str()); vResid[i]->SetMarkerStyle(3); vResid[i]->SetMarkerColor(2+i); vResid[i]->SetMarkerSize(0.55); vResid[i]->SetLineColor(46); vResid[i]->SetMinimum(-0.2); vResid[i]->SetMaximum(0.2); } cResiduals->Print("residuals5.eps"); } MSG("LIAnalysis",Msg::kInfo) <<endl<<" ** Finished the GainCurves method ** "<<endl; }

void LIAnalysis::GainCurves_Led (   )

Definition at line 11678 of file LIAnalysis.cxx. References calibPoint, chain, chip, correlatedHit, crate, detectorType, farLed, farPulserBox, Form(), GetElecString(), highRunNumber, InitialiseLoopVariables(), led, lowRunNumber, mean, MSG, nearLed, nearPulserBox, numEntries, period, pinGain, pulserBox, pulseWidth, readoutType, rms, runNumber, s(), s, and SetLoopVariables(). { MSG("LIAnalysis",Msg::kInfo) <<endl<<" ** Running the GainCurves_Led method... ** "<<endl; chain->GetEvent(0); Int_t numAdcBins=20; const Int_t const_numCalibPoints=calibType; MSG("LIAnalysis",Msg::kInfo) <<"Using "<<const_numCalibPoints<<" calibPoints"<<endl; MSG("LIAnalysis",Msg::kInfo) <<"Number of pulser boxes = "<<NUMPULSERBOXES <<", leds = "<<NUMLEDS <<endl<<endl; //store the pulse heights associated with each calibPoint Int_t **pulseHeights=0; pulseHeights=new Int_t*[NUMLEDS*NUMPULSERBOXES]; Float_t **maxAdcHighPin=0; maxAdcHighPin=new Float_t*[NUMLEDS*NUMPULSERBOXES]; Float_t **maxAdcLowPin=0; maxAdcLowPin=new Float_t*[NUMLEDS*NUMPULSERBOXES]; Float_t **maxAdcPin=0; maxAdcPin=new Float_t*[NUMLEDS*NUMPULSERBOXES]; //initialise arrays TH1F ***hAdcPh=0; hAdcPh=new TH1F**[NUMLEDS*NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ Int_t l=i*NUMLEDS+j; hAdcPh[l]=new TH1F*[const_numCalibPoints]; pulseHeights[l]=new Int_t[const_numCalibPoints]; maxAdcHighPin[l]=new Float_t[const_numCalibPoints]; maxAdcLowPin[l]=new Float_t[const_numCalibPoints]; maxAdcPin[l]=new Float_t[const_numCalibPoints]; for (Int_t k=0;k<const_numCalibPoints;k++){ string sCalibPoint=Form("%d",k); string sPulserBox=Form("%d",i); string sLed=Form("%d",j+1); s="ADC Near&Far, PB "+sPulserBox+", LED "+sLed+ ", calibPoint "+sCalibPoint; (hAdcPh[l])[k]=new TH1F(s.c_str(),s.c_str(), numAdcBins,0,15000); (hAdcPh[l])[k]->GetXaxis()->SetTitle("ADC"); (hAdcPh[l])[k]->GetXaxis()->CenterTitle(); (hAdcPh[l])[k]->GetYaxis()->SetTitle("Number of Entries"); (hAdcPh[l])[k]->GetYaxis()->CenterTitle(); (hAdcPh[l])[k]->SetFillColor(0); (hAdcPh[l])[k]->SetLineColor(1); (hAdcPh[l])[k]->Fill(1,0.0001); //(hAdcPh[l])[k]->SetBit(TH1::kCanRebin); pulseHeights[l][k]=0; maxAdcHighPin[l][k]=0; maxAdcLowPin[l][k]=0; maxAdcPin[l][k]=0; } } } TH1F ***hAdcPhNear=0; hAdcPhNear=new TH1F**[NUMLEDS*NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ Int_t l=i*NUMLEDS+j; hAdcPhNear[l]=new TH1F*[const_numCalibPoints]; for (Int_t k=0;k<const_numCalibPoints;k++){ string sCalibPoint=Form("%d",k); string sPulserBox=Form("%d",i); string sLed=Form("%d",j+1); s="ADC Near Side, PB "+sPulserBox+", LED "+sLed+ ", calibPoint "+sCalibPoint; Int_t l=i*NUMLEDS+j; hAdcPhNear[l][k]=new TH1F(s.c_str(),s.c_str(), numAdcBins,0,15000); (hAdcPhNear[l])[k]->GetXaxis()->SetTitle("ADC"); (hAdcPhNear[l])[k]->GetXaxis()->CenterTitle(); (hAdcPhNear[l])[k]->GetYaxis()->SetTitle("Number of Entries"); (hAdcPhNear[l])[k]->GetYaxis()->CenterTitle(); (hAdcPhNear[l])[k]->SetFillColor(0); (hAdcPhNear[l])[k]->SetLineColor(2); (hAdcPhNear[l])[k]->Fill(1,0.0001); //(hAdcPhNear[l])[k]->SetBit(TH1::kCanRebin); } } } TH1F ***hAdcPhFar=0; hAdcPhFar=new TH1F**[NUMLEDS*NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ Int_t l=i*NUMLEDS+j; hAdcPhFar[l]=new TH1F*[const_numCalibPoints]; for (Int_t k=0;k<const_numCalibPoints;k++){ string sCalibPoint=Form("%d",k); string sPulserBox=Form("%d",i); string sLed=Form("%d",j+1); s="ADC Far Side, PB "+sPulserBox+", LED "+sLed+ ", calibPoint "+sCalibPoint; (hAdcPhFar[l])[k]=new TH1F(s.c_str(),s.c_str(), numAdcBins,0,15000); (hAdcPhFar[l])[k]->GetXaxis()->SetTitle("ADC"); (hAdcPhFar[l])[k]->GetXaxis()->CenterTitle(); (hAdcPhFar[l])[k]->GetYaxis()->SetTitle("Number of Entries"); (hAdcPhFar[l])[k]->GetYaxis()->CenterTitle(); (hAdcPhFar[l])[k]->SetFillColor(0); (hAdcPhFar[l])[k]->SetLineColor(3); (hAdcPhFar[l])[k]->Fill(1,0.0001); //(hAdcPhFar[l])[k]->SetBit(TH1::kCanRebin); } } } TH2F **hGainCurvePmt=0; hGainCurvePmt=new TH2F*[NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ string sPb=Form("%d",i); s="Gain Curve Points in PMT ADCs, Pulser Box "+sPb; hGainCurvePmt[i]=new TH2F(s.c_str(),s.c_str(),21,0,21,200,0,15000); hGainCurvePmt[i]->GetXaxis()->SetTitle("LED"); hGainCurvePmt[i]->GetXaxis()->CenterTitle(); hGainCurvePmt[i]->GetYaxis()->SetTitle("Pmt ADCs"); hGainCurvePmt[i]->GetYaxis()->CenterTitle(); hGainCurvePmt[i]->SetFillColor(0); hGainCurvePmt[i]->SetBit(TH1::kCanRebin); } TH2F **hGainCurvePin=0; hGainCurvePin=new TH2F*[NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ string sPb=Form("%d",i); s="Gain Curve Points in Pin diode ADCs, Pulser Box "+sPb; hGainCurvePin[i]=new TH2F(s.c_str(),s.c_str(),21,0,21,200,0,15000); hGainCurvePin[i]->GetXaxis()->SetTitle("LED"); hGainCurvePin[i]->GetXaxis()->CenterTitle(); hGainCurvePin[i]->GetYaxis()->SetTitle("Pin ADCs"); hGainCurvePin[i]->GetYaxis()->CenterTitle(); hGainCurvePin[i]->SetFillColor(0); hGainCurvePin[i]->SetBit(TH1::kCanRebin); } TH2F **hGainCurvePh=0; hGainCurvePh=new TH2F*[NUMPULSERBOXES]; for (Int_t i=0;i<NUMPULSERBOXES;i++){ string sPb=Form("%d",i); s="Gain Curve Points Pulse Height, Pulser Box "+sPb; hGainCurvePh[i]=new TH2F(s.c_str(),s.c_str(),21,0,21,200,0,1023); hGainCurvePh[i]->GetXaxis()->SetTitle("LED"); hGainCurvePh[i]->GetXaxis()->CenterTitle(); hGainCurvePh[i]->GetYaxis()->SetTitle("Pulse Height"); hGainCurvePh[i]->GetYaxis()->CenterTitle(); hGainCurvePh[i]->SetFillColor(0); hGainCurvePh[i]->SetBit(TH1::kCanRebin); } this->InitialiseLoopVariables(); for(Int_t entry=0;entry<numEvents;entry++){ this->SetLoopVariables(entry,1); //ignore any zeros that slipped through if (mean==0 || rms==0 || numEntries==0) continue; Int_t c=calibPoint-1; //look at pins if (readoutType==ReadoutType::kPinDiode){ if (pulserBox>=FIRSTPULSERBOX && pulserBox<=LASTPULSERBOX && led>=FIRSTLED && led<=LASTLED){ //define led number Int_t l=pulserBox*NUMLEDS+led-1; //will require a correlated hit when plex is working //fill histo for appropriate pin if (chip==1){//high gain (chip 1) //find max value for high gain pin if (mean>maxAdcHighPin[l][c]) { maxAdcHighPin[l][c]=mean; } } else if (chip==0){//low gain (chip 0) //find max value for low gain pin if (mean>maxAdcLowPin[l][c]) { maxAdcLowPin[l][c]=mean; } } //print out warnings if there are inconsistencies if ((pinGain!=1 && chip==0) || (pinGain!=0 && chip==1)) { MSG("LIAnalysis",Msg::kInfo) <<"**** Strange pin, wrong gain in plex, (" <<", "<<this->GetElecString() <<", plex gain="<<pinGain<<", mean="<<mean<<", rms="<<rms <<", num="<<numEntries<<", cHit="<<correlatedHit<<endl; } } } //only look at scint strips if (readoutType!=ReadoutType::kScintStrip) continue; //fill histograms if (correlatedHit==1){ if (detectorType==Detector::kFar){ if (pulserBox>=FIRSTPULSERBOX && pulserBox<=LASTPULSERBOX && led>=FIRSTLED && led<=LASTLED && calibPoint>=1 && calibPoint<=const_numCalibPoints){ Int_t l=pulserBox*NUMLEDS+led-1; pulseHeights[l][c]=pulseHeight; //fill histo for appropriate led hAdcPh[l][c]->Fill(mean); //fill histos with when pulser box is near pulser box if (nearPulserBox==pulserBox){ hAdcPhNear[l][c]->Fill(mean); } //fill histos with when pulser box is far pulser box else if (farPulserBox==pulserBox) { hAdcPhFar[l][c]->Fill(mean); } } } //change the farPb to farLed for calDet else if (detectorType==Detector::kCalDet){ if (pulserBox>=FIRSTPULSERBOX && pulserBox<=LASTPULSERBOX && crate>=FIRSTCRATE && crate<=LASTCRATE && led>=FIRSTLED && led<=LASTLED && calibPoint>=1 && calibPoint<=const_numCalibPoints){ Int_t l=pulserBox*NUMLEDS+led-1; pulseHeights[l][c]=pulseHeight; //fill histo for appropriate led hAdcPh[l][c]->Fill(mean); //fill histos with when pulser box is near pulser box if (nearLed==led){ hAdcPhNear[l][c]->Fill(mean); } //fill histos with when pulser box is far pulser box else if (farLed==led) { hAdcPhFar[l][c]->Fill(mean); } } } } }//end of for MSG("LIAnalysis",Msg::kInfo)<<"Finished main loop"<<endl; //include the under and overflow counts gStyle->SetOptStat(1111111); //set up useful string string sRunNumber=Form("%d",runNumber); string sLowRunNumber=Form("%d",lowRunNumber); string sHighRunNumber=Form("%d",highRunNumber); s=""; //set strings for use in titles string sPulseWidth=Form("%d",pulseWidth); string sPulseFreq=Form("%d",static_cast<Int_t> (ceil(1.0/(period*1.0e-5)))); //allocate space for TGraphs then loop and create them all TGraph **gAdcVsPinNear=0; TGraph **gAdcVsPinFar=0; gAdcVsPinNear=new TGraph*[NUMPULSERBOXES*NUMLEDS]; gAdcVsPinFar=new TGraph*[NUMPULSERBOXES*NUMLEDS]; TGraph **gNonlinearity=0; gNonlinearity=new TGraph*[NUMPULSERBOXES*NUMLEDS]; TGraph **gResiduals=0; gResiduals=new TGraph*[NUMPULSERBOXES*NUMLEDS]; for ( Int_t l=0;l<NUMPULSERBOXES*NUMLEDS;l++){ gAdcVsPinNear[l]=new TGraph(const_numCalibPoints); gAdcVsPinFar[l]=new TGraph(const_numCalibPoints); gNonlinearity[l]=new TGraph(const_numCalibPoints); gResiduals[l]=new TGraph(const_numCalibPoints); //set minimum on all graphs gAdcVsPinNear[l]->SetMinimum(-1); gAdcVsPinFar[l]->SetMinimum(-1); //set maximum on all graphs gAdcVsPinNear[l]->SetMaximum(15000); gAdcVsPinFar[l]->SetMaximum(15000); } Int_t maxAdc=-1; Float_t* maxRatio=new Float_t[NUMPULSERBOXES*NUMLEDS]; //find the max ratio for each point for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ Int_t l=i*NUMLEDS+j; maxRatio[l]=-1; for (Int_t k=0;k<const_numCalibPoints;k++){ if (maxAdcLowPin[l][k]>0){ if (hAdcPhNear[l][k]->GetMean()/maxAdcLowPin[l][k]> maxRatio[l]){ maxRatio[l]=hAdcPhNear[l][k]->GetMean()/maxAdcLowPin[l][k]; } } } } } //fill graphs for (Int_t i=0;i<NUMPULSERBOXES;i++){ for (Int_t j=0;j<NUMLEDS;j++){ Int_t l=i*NUMLEDS+j; for (Int_t k=0;k<const_numCalibPoints;k++){ //set points for graphs gAdcVsPinNear[l]->SetPoint(k,maxAdcLowPin[l][k], hAdcPhNear[l][k]->GetMean()); //what is the point of this one? gAdcVsPinFar[l]->SetPoint(k,maxAdcLowPin[l][k], hAdcPhFar[l][k]->GetMean()); if (maxAdcLowPin[l][k]>0){ //do pmt_adc/pin_adc VS pin_adc gNonlinearity[l]->SetPoint(k,maxAdcLowPin[l][k], (hAdcPhNear[l][k]->GetMean()/ maxAdcLowPin[l][k])/maxRatio[l]); } else{ gNonlinearity[l]->SetPoint(k,maxAdcLowPin[l][k],-1); MSG("LIAnalysis",Msg::kInfo) <<"("<<i<<":"<<j+1<<") Zero point"<<endl; } //get max for plot if (hAdcPhNear[l][k]->GetMean()>maxAdc) { maxAdc=static_cast<Int_t>(hAdcPhNear[l][k]->GetMean()); } } } } Axis_t pinMax[6]={1500,1200,1700,900,1600,1000}; //calculate the residuals //and fill the residuals graph for (Int_t i=0;i<NUMPULSERBOXES;i++){ MSG("LIAnalysis",Msg::kInfo)<<"Doing plots for PB "<<i<<endl; for (Int_t j=0;j<NUMLEDS;j++){ Int_t l=i*NUMLEDS+j; //a quick and nasty hack to get results for caldet if (l<6){ gAdcVsPinNear[l]->Fit("pol1","q","",10.0,pinMax[l]); } else{ gAdcVsPinNear[l]->Fit("pol1","q","",10.0,1500); } TF1* func=gAdcVsPinNear[l]->GetFunction("pol1"); Double_t funcC=func->GetParameter(0); Double_t funcM=func->GetParameter(1); gAdcVsPinNear[l]->Draw("AP"); Axis_t xMax=gAdcVsPinNear[l]->GetXaxis()->GetXmax(); Axis_t xMin=gAdcVsPinNear[l]->GetXaxis()->GetXmin(); MSG("LIAnalysis",Msg::kInfo) <<"("<<i<<":"<<j+1<<") y="<<funcM<<"x+"<<funcC <<", xMin="<<xMin<<", xMax="<<xMax<<endl; for (Int_t k=0;k<const_numCalibPoints;k++){ Double_t actualX=-1; Double_t actualYval=-1; gAdcVsPinNear[l]->GetPoint(k,actualX,actualYval); //do y=mx+c to get value on straight line fit Double_t fittedYval=funcM*actualX+funcC; //calculate the residual and normalise Double_t resid=actualYval-fittedYval; resid/=fittedYval; //set the point on the graph if (actualX>0){ gResiduals[l]->SetPoint(k,actualX,resid); } else { MSG("LIAnalysis",Msg::kInfo) <<"Setting a zero in residuals"<<endl; gResiduals[l]->SetPoint(k,actualX,0); } MSG("LIAnalysis",Msg::kInfo) <<"("<<i<<":"<<j+1<<") k="<<k<<", x="<<actualX <<", y="<<actualYval <<", fitY="<<fittedYval<<", resid="<<resid<<i<<endl; } } } //create canvas TCanvas *c=new TCanvas("c","c",0,0,1200,600); c->SetFillColor(0); //turn off the fit stats box gStyle->SetOptFit(0000); //plot the graphs for (Int_t i=0;i<NUMPULSERBOXES;i++){ MSG("LIAnalysis",Msg::kInfo)<<"Doing plots for PB "<<i<<endl; for (Int_t j=0;j<NUMLEDS;j++){ Int_t l=i*NUMLEDS+j; //draw the graphs of the pin adc vs pulse height c->cd(); c->Clear(); string sPulserBox=Form("%d",i); string sLed=Form("%d",j+1); s="Gain Curve (PB="+sPulserBox+ ", LED="+sLed+", PW="+sPulseWidth+", PF="+sPulseFreq+" Hz)"; gAdcVsPinNear[l]->Draw("AP"); gAdcVsPinNear[l]->SetTitle(s.c_str()); gAdcVsPinNear[l]->GetXaxis()->SetTitle("Pin Adc"); gAdcVsPinNear[l]->GetYaxis()->SetTitle("Pmt Adc"); gAdcVsPinNear[l]->GetXaxis()->CenterTitle(); gAdcVsPinNear[l]->GetYaxis()->CenterTitle(); gAdcVsPinNear[l]->SetMarkerStyle(3); gAdcVsPinNear[l]->SetMarkerColor(2); gAdcVsPinNear[l]->SetMarkerSize(0.55); gAdcVsPinNear[l]->SetLineColor(46); gAdcVsPinNear[l]->SetMinimum(-1); gAdcVsPinNear[l]->SetMaximum(maxAdc); s="GainCurves_Led.ps"; //print graph to postscript if(i*NUMLEDS+j==0){ if (sLowRunNumber==sHighRunNumber) s=sHighRunNumber+s+"("; else s=sLowRunNumber+"-"+sHighRunNumber+s+"("; c->Print(s.c_str()); } else if(i*NUMLEDS+j==NUMLEDS*NUMPULSERBOXES-1){ gErrorIgnoreLevel=0; if (sLowRunNumber==sHighRunNumber) s=sHighRunNumber+s+")"; else s=sLowRunNumber+"-"+sHighRunNumber+s+")"; c->Print(s.c_str()); MSG("LIAnalysis",Msg::kInfo)<<"Done last plot"<<endl; } else{ gErrorIgnoreLevel=1; if (sLowRunNumber==sHighRunNumber) s=sHighRunNumber+s; else s=sLowRunNumber+"-"+sHighRunNumber+s; c->Print(s.c_str()); } } } //plot the graphs for (Int_t i=0;i<NUMPULSERBOXES;i++){ MSG("LIAnalysis",Msg::kDebug)<<"Doing plots for PB "<<i<<endl; for (Int_t j=0;j<NUMLEDS;j++){ Int_t l=i*NUMLEDS+j; MSG("LIAnalysis",Msg::kDebug)<<"Doing plots for led "<<j+1<<endl; //draw the graphs of the pin adc vs pulse height c->cd(); c->Clear(); string sPulserBox=Form("%d",i); string sLed=Form("%d",j+1); s="PMT Nonlinearity (PB="+sPulserBox+ ", LED="+sLed+", PW="+sPulseWidth+", PF="+sPulseFreq+" Hz)"; gNonlinearity[l]->Draw("AP"); gNonlinearity[l]->SetTitle(s.c_str()); gNonlinearity[l]->GetXaxis()->SetTitle("Pin Adc"); gNonlinearity[l]->GetYaxis()->SetTitle("Pmt Adc / Pin Adc"); gNonlinearity[l]->GetXaxis()->CenterTitle(); gNonlinearity[l]->GetYaxis()->CenterTitle(); gNonlinearity[l]->SetMarkerStyle(3); gNonlinearity[l]->SetMarkerColor(2); gNonlinearity[l]->SetMarkerSize(0.35); gNonlinearity[l]->SetLineColor(46); //gNonlinearity[l]->SetMinimum(-1); //gNonlinearity[l]->SetMaximum(maxAdc); s="GainCurvesNonLin_Led.ps"; //print graph to postscript if(i*NUMLEDS+j==0){ if (sLowRunNumber==sHighRunNumber) s=sHighRunNumber+s+"("; else s=sLowRunNumber+"-"+sHighRunNumber+s+"("; c->Print(s.c_str()); } else if(i*NUMLEDS+j==NUMLEDS*NUMPULSERBOXES-1){ gErrorIgnoreLevel=0; if (sLowRunNumber==sHighRunNumber) s=sHighRunNumber+s+")"; else s=sLowRunNumber+"-"+sHighRunNumber+s+")"; c->Print(s.c_str()); MSG("LIAnalysis",Msg::kInfo)<<"Done last plot"<<endl; } else{ gErrorIgnoreLevel=1; if (sLowRunNumber==sHighRunNumber) s=sHighRunNumber+s; else s=sLowRunNumber+"-"+sHighRunNumber+s; c->Print(s.c_str()); } } } //plot the graphs for (Int_t i=0;i<NUMPULSERBOXES;i++){ MSG("LIAnalysis",Msg::kDebug)<<"Doing plots for PB "<<i<<endl; for (Int_t j=0;j<NUMLEDS;j++){ Int_t l=i*NUMLEDS+j; MSG("LIAnalysis",Msg::kDebug)<<"Doing plots for led "<<j+1<<endl; //draw the graphs of the pin adc vs pulse height c->cd(); c->Clear(); string sPulserBox=Form("%d",i); string sLed=Form("%d",j+1); s="Residuals (PB="+sPulserBox+ ", LED="+sLed+", PW="+sPulseWidth+", PF="+sPulseFreq+" Hz)"; gResiduals[l]->Draw("AP"); gResiduals[l]->SetTitle(s.c_str()); gResiduals[l]->GetXaxis()->SetTitle("Pin Adc"); gResiduals[l]->GetYaxis()->SetTitle("Residual"); gResiduals[l]->GetXaxis()->CenterTitle(); gResiduals[l]->GetYaxis()->CenterTitle(); gResiduals[l]->SetMarkerStyle(3); gResiduals[l]->SetMarkerColor(2); gResiduals[l]->SetMarkerSize(0.55); gResiduals[l]->SetLineColor(46); gResiduals[l]->SetMinimum(-1); gResiduals[l]->SetMaximum(1); s="GainCurvesResid_Led.ps"; //print graph to postscript if(i*NUMLEDS+j==0){ if (sLowRunNumber==sHighRunNumber) s=sHighRunNumber+s+"("; else s=sLowRunNumber+"-"+sHighRunNumber+s+"("; c->Print(s.c_str()); } else if(i*NUMLEDS+j==NUMLEDS*NUMPULSERBOXES-1){ gErrorIgnoreLevel=0; if (sLowRunNumber==sHighRunNumber) s=sHighRunNumber+s+")"; else s=sLowRunNumber+"-"+sHighRunNumber+s+")"; c->Print(s.c_str()); MSG("LIAnalysis",Msg::kInfo)<<"Done last plot"<<endl; } else{ gErrorIgnoreLevel=1; if (sLowRunNumber==sHighRunNumber) s=sHighRunNumber+s; else s=sLowRunNumber+"-"+sHighRunNumber+s; c->Print(s.c_str());