#pragma rtGlobals=3		// Use modern global access method and strict wave access.
Function SumIntensity_traj(TotalPar) // to calculate the total intensity of all particles detected at each frame from SPT to plot total intensity vs. frame
	// setting up variables	
	variable TotalPar	//the number for all N Par
	variable FirstTime	//for specifying where to store the extracted data
	variable TrajLength  //the length of wave for ParN
	Make /o/N=100 Total_int	//the new wave to store the extracted data
	
	// looping all particles
	variable i, j
	
	for(i=0;i<TotalPar;i+=1)
		SetDataFolder "root:data:Traces:Par"+num2str(i) //select folder for certain Par
		TrajLength=DimSize($"Vol",0)
		NVAR FisrtTime=$"StartAndEndOfTrace(1)"	//specify where to store extracted data - is this correct?
		
		// looping all traaces of the i particle
		for (j=0;j<TrajLength;j+=1)
			variable temp=FirstTime+j //element of where to start data storagev
			String timetag="Vol["+num2str(j) //I was trying to figure out how to loop elements here Vol(j)...this is wrong syntax I think..can you help?
			timetag += "]"
			Total_int[temp] += $timetag
		endfor
		
	endfor
	
End