nhits1 = HITS1.hit_num
h_Pedestal = TFile("Pedestal","",100,150,250)
for ev in range(nevents1):
fChain1.GetEntry(ev)
for i in range(NSCA):
TS[i] = HITS.rollposition[i]
for hit in range(nhits):
H=HITS1.Hits.at(nhits1)
if(H.CCorNC == 0):
continue
if(H.chip == 0 and H.ch == 20):
if(TS[i]<9):
h_Pedestal.Fill(H.SCA_hg[0])
h_Pedestal.Fit(Full_Gaussian1,"R")
Ped_Chi=Full_Gaussian1.GetChisquare()
Ped_NDF=Full_Gaussian1.GetNDF()
Ped_Degree_part1=Ped_Chi/Ped_NDF
Ped_Mean=Full_Gaussian1.GetParameter(1)
Ped_Error=Full_Gaussian1.GetParameter(2)
#=======================================For the injection channel analysis=============================#
c2 = TCanvas("c1", "c1",0,0,1000,1000)
nevents2= fChain2.GetEntries()
nhits2 = HITS2.hit_num
h_Charge = TFile("Charge","",100,150,250)
for ev in range(nevents2):
fChain2.GetEntry(ev)
for i in range(NSCA):
