Adding new Vendeta macro for RawSpectra checking

Projects/Vendeta/macro/RawSpectra/ShowRawFCSpectra.C

+TChain* chain;
+
+int NumberOfAnodes= 1;
+int NumberOfEvents= 1e6;
+
+//////////////////////////////////////////////////
+void OpenRootFile(){
+		chain = new TChain("RawTree");
+		chain->Add("/home/faster/fastertonptool/data/rootfiles/V4B_SAMPLING_6_0001.root");
+}
+
+//////////////////////////////////////////////////
+void ShowRawFCSpectra(string Nucleus="Cf"){
+		OpenRootFile();
+
+		if(Nucleus=="Cf" || Nucleus=="252Cf")
+				NumberOfAnodes=1;
+		if(Nucleus=="U" || Nucleus=="238U")
+				NumberOfAnodes=11;
+
+		// Canvas Definition for Low Gain //
+		TCanvas* c1   = new TCanvas("Charge Q1","Charge Q1",1800,1800);
+		TCanvas* c2   = new TCanvas("Charge Q2","Charge Q2",1800,1800);
+		TCanvas* c3   = new TCanvas("Charge Qmax","Charge Qmax",1800,1800);
+		TCanvas* c4   = new TCanvas("Q2 vs Q1","Q2 vs Q1",1800,1800);
+		TCanvas* c5   = new TCanvas("Qmax vs Q1","Qmax vs Q1",1800,1800);
+
+		if(Nucleus=="U" || Nucleus=="238U"){
+				c1->Divide(3,4);
+				c2->Divide(3,4);
+				c3->Divide(3,4);
+				c4->Divide(3,4);
+				c5->Divide(3,4);
+		}
+
+		for(int i=0; i<NumberOfAnodes; i++){
+				// Draw //
+				TString draw_Q1 = Form("fFC_Q1>>hQ1_Anode%i(500,0,100e3)",i+1);
+				TString draw_Q2 = Form("fFC_Q2>>hQ2_Anode%i(500,0,20e3)",i+1);
+				TString draw_Qmax = Form("fFC_Qmax>>hQmax_Anode%i(500,0,10e3)",i+1);
+				TString draw_Q2vsQ1 = Form("fFC_Q2:fFC_Q1>>hQ2vsQ1_Anode%i(500,0,100e3,500,0,20e3)",i+1);
+				TString draw_QmaxvsQ1 = Form("fFC_Qmax:fFC_Q1>>hmax2vsQ1_Anode%i(500,0,100e3,500,0,10e3)",i+1);
+
+				TString condition;
+				if(Nucleus=="Cf" || Nucleus=="252Cf")
+						condition = Form("fFC_AnodeNbr==%i",6);
+
+				if(Nucleus=="U" || Nucleus=="238U")
+						condition = Form("fFC_AnodeNbr==%i",i+1);
+
+				if(Nucleus=="U")c1->cd(i+1);
+				else c1->cd();
+				chain->Draw(draw_Q1,condition,"",NumberOfEvents);
+				c1->Update();
+
+				if(Nucleus=="U")c2->cd(i+1);
+				else c2->cd();
+				chain->Draw(draw_Q2,condition,"",NumberOfEvents);
+
+				if(Nucleus=="U")c3->cd(i+1);
+				else c3->cd();
+				chain->Draw(draw_Qmax,condition,"",NumberOfEvents);
+
+				if(Nucleus=="U")c4->cd(i+1)->SetLogz();
+				else c4->cd()->SetLogz();
+				chain->Draw(draw_Q2vsQ1,condition,"colz",NumberOfEvents);
+
+				if(Nucleus=="U") c5->cd(i+1)->SetLogz();
+				else c5->cd()->SetLogz();
+				chain->Draw(draw_QmaxvsQ1,condition,"colz",NumberOfEvents);
+		}
+}

Projects/Vendeta/macro/RawSpectra/ShowRawPSD.C

+TChain* chain;
+
+int NumberOfDetectors= 2;
+int NumberOfAnodes= 1;
+int NumberOfEvents= 1e6;
+
+//////////////////////////////////////////////////
+void OpenRootFile(){
+		chain = new TChain("RawTree");
+		chain->Add("/home/faster/fastertonptool/data/rootfiles/V4B_SAMPLING_6_0001.root");
+}
+
+//////////////////////////////////////////////////
+void ShowRawPSD(){
+		OpenRootFile();
+
+		// Canvas Definition for Low Gain //
+		TCanvas* cLG_RI   = new TCanvas("Det-LG RI","Det-LG RI",1800,1800);
+		TCanvas* cLG_RII  = new TCanvas("Det-LG RII","Det-LG RII",1800,1800);
+		TCanvas* cLG_RIII = new TCanvas("Det-LG RIII","Det-LG RIII",1800,1800);
+		TCanvas* cLG_LI   = new TCanvas("Det-LG LI","Det-LG LI",1800,1800);
+		TCanvas* cLG_LII  = new TCanvas("Det-LG LII","Det-LG LII",1800,1800);
+		TCanvas* cLG_LIII = new TCanvas("Det-LG LIII","Det-LG LIII",1800,1800);
+
+		cLG_RI->Divide(3,4);
+		cLG_RII->Divide(3,4);
+		cLG_RIII->Divide(3,4);
+		cLG_LI->Divide(3,4);
+		cLG_LII->Divide(3,4);
+		cLG_LIII->Divide(3,4);
+
+		// Canvas Definition for High Gain //
+		TCanvas* cHG_RI   = new TCanvas("Det-HG RI","Det-HG RI",1800,1800);
+		TCanvas* cHG_RII  = new TCanvas("Det-HG RII","Det-HG RII",1800,1800);
+		TCanvas* cHG_RIII = new TCanvas("Det-HG RIII","Det-HG RIII",1800,1800);
+		TCanvas* cHG_LI   = new TCanvas("Det-HG LI","Det-HG LI",1800,1800);
+		TCanvas* cHG_LII  = new TCanvas("Det-HG LII","Det-HG LII",1800,1800);
+		TCanvas* cHG_LIII = new TCanvas("Det-HG LIII","Det-HG LIII",1800,1800);
+
+		cHG_RI->Divide(3,4);
+		cHG_RII->Divide(3,4);
+		cHG_RIII->Divide(3,4);
+		cHG_LI->Divide(3,4);
+		cHG_LII->Divide(3,4);
+		cHG_LIII->Divide(3,4);
+
+		for(int i=0; i<NumberOfDetectors; i++){
+				// LG //	
+				TString to_draw_LG = Form("fVendeta_LG_Q2/fVendeta_LG_Q1:fVendeta_LG_Q1>>hLG_det%i(500,0,500e3,500,0,1)",i+1);
+				TString condition_LG = Form("fVendeta_LG_DetectorNbr==%i",i+1);
+
+				// HG //
+				TString to_draw_HG = Form("fVendeta_HG_Q2/fVendeta_HG_Q1:fVendeta_HG_Q1>>hHG_det%i(500,0,800e3,500,0,1)",i+1);
+				TString condition_HG = Form("fVendeta_HG_DetectorNbr==%i",i+1);
+
+				if(i+1<12){
+						cLG_RI->cd(i+1);
+						chain->Draw(to_draw_LG,condition_LG,"colz",NumberOfEvents);
+
+						cHG_RI->cd(i+1);
+						chain->Draw(to_draw_HG,condition_HG,"colz",NumberOfEvents);
+				}
+				if(i+1>12 && i+1<24){
+						cLG_RII->cd(i+1);
+						chain->Draw(to_draw_LG,condition_LG,"colz",NumberOfEvents);
+
+						cHG_RII->cd(i+1);
+						chain->Draw(to_draw_HG,condition_HG,"colz",NumberOfEvents);
+				}
+				if(i+1>24 && i+1<36){
+						cLG_RIII->cd(i+1);
+						chain->Draw(to_draw_LG,condition_LG,"colz",NumberOfEvents);
+
+						cHG_RIII->cd(i+1);
+						chain->Draw(to_draw_HG,condition_HG,"colz",NumberOfEvents);
+				}
+				if(i+1>46 && i+1<48){
+						cLG_LI->cd(i+1);
+						chain->Draw(to_draw_LG,condition_LG,"colz",NumberOfEvents);
+
+						cHG_LI->cd(i+1);
+						chain->Draw(to_draw_HG,condition_HG,"colz",NumberOfEvents);
+				}
+				if(i+1>48 && i+1<60){
+						cLG_LII->cd(i+1);
+						chain->Draw(to_draw_LG,condition_LG,"colz",NumberOfEvents);
+
+						cHG_LII->cd(i+1);
+						chain->Draw(to_draw_HG,condition_HG,"colz",NumberOfEvents);
+				}
+				if(i+1>60){
+						cLG_LIII->cd(i+1);
+						chain->Draw(to_draw_LG,condition_LG,"colz",NumberOfEvents);
+				
+						cHG_LIII->cd(i+1);
+						chain->Draw(to_draw_HG,condition_HG,"colz",NumberOfEvents);
+				}
+
+		}
+
+}