Collected all scripts in macros/exp1016, all raw files in clb_raw_files, run...

Collected all scripts in macros/exp1016, all raw files in clb_raw_files, run folders with calibration files are also added

calibration1.cxx → macros/exp1016/calibration1.cxx

 {
-//	gSystem->Load("/home/dariak/AculUtils/libAculData.so");
-	gSystem->Load("/home/vratik/workspace/AculUtils/libAculData.so");
+	gSystem->Load("/home/dariak/AculUtils/libAculData.so");
+//	gSystem->Load("/home/vratik/workspace/AculUtils/libAculData.so");
 
 	AculCalibration cal;
 	cal.SetInputParameters(); //from .par
 //	cal.CalculateCalibParameters("clb08_0001.root", "SQ22", "AnalysisxTree", 100, 4095);
-	cal.CalculateCalibParameters("macros/myMacros/clb01_0001.root", "SQ22", "AnalysisxTree", 100, 4095);
+	cal.CalculateCalibParameters("clb08_0001.root", "SQ22", "AnalysisxTree", 100, 4095);
 
 //	    CalculateCalibParameters(const char* inputfile, const char* block, const Int_t address, const char* treename, Int_t lowerchannel = 0, Int_t upperchannel = 4095, Int_t nEBins = 1000, Int_t lowersubaddress = 0, Int_t uppersubaddress = ADDRESSNUMBER-1);
 

calibration2.cxx → macros/exp1016/calibration2.cxx

@@ -9,8 +9,9 @@ void calibration2()
 	cal->PrintCalibrationParameters(22, 22);
 
 	TCanvas *c1 = new TCanvas();
+	TCanvas *c2 = new TCanvas();
 	cal->CalibrateRawSpectra("clb09_0001.root", "SQ22", "AnalysisxTree", 100, 800, 500); //takes data from raw file and calibrates it with obtained calibration parameters
-	cal->CalibrateRawSpectra("clb09_0001.root", "SQ22", "AnalysisxTree", 100, 800, 500); //takes data from raw file and calibrates it with obtained calibration parameters
+	//cal->CalibrateRawSpectra("clb09_0001.root", "SQ22", "AnalysisxTree", 100, 800, 500); //takes data from raw file and calibrates it with obtained calibration parameters
 //return;
 	TFile fr("SQ22[]E.root");
 	cout << fr.GetListOfKeys()->GetEntries() << " histograms" << endl;
@@ -28,6 +29,8 @@ void calibration2()
 		hWork->Draw();
 	}
 
+	cal.FindEnergyPeaks(c2,"SQ22[]E.root","output.txt");
+
 /*	c1->Divide(2, 2);
 
 	c1->cd(1);

macros/exp1016/deadlayers.cxx

+{
+	gSystem->Load("../libTELoss.so");
+
+	TELoss	mSi;
+//set parameters for Si and alpha
+	mSi.SetEL(1, 2.321); // density in g/cm3
+	mSi.AddEL(14., 28.086, 1);  //Z, mass
+//	mSi.SetZP(1., 1.);		//protons
+	mSi.SetZP(2., 4.);		//alphas, Z, A
+	mSi.SetEtab(100000, 200.);	// ?, MeV calculate ranges
+	mSi.SetDeltaEtab(300);        // calculate delta E
+
+//	TCanvas *c1 = new TCanvas("name", "title which can consist more than one word.", 617, 0, 1058, 972);
+//	cout << mSi.GetE(50, 1000) << endl;	//(input energy E0 in MeV, microns)  
+//	cout << mSi.GetE0(9.04, 1000) << endl;	//(output energy E in MeV, microns) 
+
+//calculated energies for run1 - run3 
+/*
+	Double_t deltaL1 = mSi.GetR(4.789, 4.623); //(doule E0, double E in MeV) calculates layer in mcm for averages 
+	Double_t deltaL2 = mSi.GetR(5.5,5.355);  
+	Double_t deltaL3 = mSi.GetR(6.017,5.880); 
+	Double_t deltaL4 = mSi.GetR(7.693, 7.577); 
+*/
+//calculated energies for run4 - run5
+/*
+	Double_t deltaL1 = mSi.GetR(4.782,4.613); //(doule E0, double E in MeV) calculates layer in mcm for averages 
+	Double_t deltaL2 = mSi.GetR(5.503,5.348);  
+	Double_t deltaL3 = mSi.GetR(6.016,5.875); 
+	Double_t deltaL4 = mSi.GetR(7.692,7.571); 
+*/
+
+//calculated energies for run6 - run7
+/*
+	Double_t deltaL1 = mSi.GetR(4.788,4.675); //(doule E0, double E in MeV) calculates layer in mcm for averages 
+	Double_t deltaL2 = mSi.GetR(5.502,5.340);  
+	Double_t deltaL3 = mSi.GetR(6.017,5.923); 
+	Double_t deltaL4 = mSi.GetR(7.693,7.612); 
+*/
+//calculated energies for run8 - run9
+
+	Double_t deltaL1 = mSi.GetR(4.787,4.648); //(doule E0, double E in MeV) calculates layer in mcm for averages 
+	Double_t deltaL2 = mSi.GetR(5.501,5.380);  
+	Double_t deltaL3 = mSi.GetR(6.018,5.903); 
+	Double_t deltaL4 = mSi.GetR(7.692,7.594); 
+
+	cout << deltaL1 << " mcm" << endl;	//MeV, microns -  delta layer
+	cout << deltaL2 << " mcm" << endl;	
+	cout << deltaL3 << " mcm" << endl;	
+	cout << deltaL4 << " mcm" << endl << endl;	
+
+	cout << "dead layer 1 is: " << deltaL1/(TMath::Sqrt(2) - 1) << " mcm" << endl;	//MeV, microns
+	cout << "dead layer 2 is: " <<deltaL2/(TMath::Sqrt(2) - 1) << " mcm" << endl;	//MeV, microns
+	cout << "dead layer 3 is: " <<deltaL3/(TMath::Sqrt(2) - 1) << " mcm" << endl;	//MeV, microns
+	cout << "dead layer 4 is: " <<deltaL4/(TMath::Sqrt(2) - 1) << " mcm" << endl << endl;	//MeV, microns
+
+	cout << "Back dead layer for detector 1-11 is : " << ( deltaL1/(TMath::Sqrt(2) - 1) + deltaL2/(TMath::Sqrt(2) - 1)+ deltaL3/(TMath::Sqrt(2) - 1) + deltaL4/(TMath::Sqrt(2) - 1) )/4. << " mcm" << endl; //dead layer width
+
+}

macros/exp1016/one_strip_raw.cxx

+{
+	gSystem->Load("/home/dariak/AculUtils/libAculData.so");
+
+	AculCalibration cal;
+	cal.SetInputParameters(); //from .par
+	cal.CalculateCalibParameters("clb09_0001.root", "SQ22", "AnalysisxTree", 100, 4095);
+
+//	    CalculateCalibParameters(const char* inputfile, const char* block, const Int_t address, const char* treename, Int_t lowerchannel = 0, Int_t upperchannel = 4095, Int_t nEBins = 1000, Int_t lowersubaddress = 0, Int_t uppersubaddress = ADDRESSNUMBER-1);
+
+	TFile fr("SQ22[].root");
+	cout << fr.GetListOfKeys()->GetEntries() << " histograms" << endl;
+	TList *histList = fr.GetListOfKeys();
+	TH1 *hWork = 0;
+	for (Int_t i = 0; i < 32; i++) 
+	{
+		fr.GetObject(histList->At(i)->GetName(), hWork);		
+		cal->PeaksFitting(hWork);
+	}
+}