{

gROOT.Reset();

TFile f1("/Users/Ormes/data/OrmesHiLatDiffuse1day_merit.root");
TTree *t1=(TTree*)f1.Get("MeritTuple");

Long64_t nentries = (Long64_t)MeritTuple->GetEntries();    //Get the number of entries in the "MeritTuple"

Double_t PtTime;											//declare variables from tuple
Float_t PtMagLat; 
Float_t PtPos[3];

MeritTuple->SetBranchAddress("PtTime",&PtTime);
MeritTuple->SetBranchAddress("PtMagLat",&PtMagLat);

for (Int_t j=0; j<3; j++){
MeritTuple->SetBranchAddress(Form("PtPos[%d]",j), &PtPos[j]); 
//MeritTuple->SetBranchAddress("PtPos[j]", &PtPos[j]); 
}


MeritTuple->SetBranchStatus("*",0);						//switch all variables off

MeritTuple->SetBranchStatus("PtTime",1);				//turn on only those we want
MeritTuple->SetBranchStatus("PtMagLat",1);

for (Int_t j=0; j<3; j++){
MeritTuple->SetBranchStatus(Form("PtPos[%d]",j),1); 
}


Double_t Time0=0;
Double_t TimeOut=0;
Double_t PreviousTime;
Double_t TimeInterval;
Double_t TimeSum=0;
Long64_t EventsIn=0;
Long64_t EventsOut=0;

Long64_t Counter0=0;

//Float_t PtDecEarthZenith = "asin(PtPos[2]/sqrt(PtPos[0]^2 + PtPos[1]^2 + PtPos[2]^2))";

TCut Tend = "PtTime> 243700000 && PtTime < 243800000";

for(Long64_t  i=1; i<nentries+1; i++){ //for every event

		MeritTuple->GetEntry(i-1); //grab the previous event
		PreviousTime=PtTime;
		MeritTuple->GetEntry(i);   //Grab current event
		//remember the previous time
		TimeInterval = PtTime-PreviousTime;
		
if (TimeInterval<75 && PtTime> 239700000 && PtTime < 245700000 )
	{
		TimeSum = TimeSum+TimeInterval;
		EventsIn = EventsIn + 1;
	 }
	 
}

//Counter,time displayed in the Console
cout<<endl<<" Sort PtMcIlwainB"<<endl; 
cout<<"Counts  PtLat>0   PtLat<0"<<endl;  

cout<<endl<<"Time (sec)  PtLat>0   PtLat<0"<<endl;  

cout<<endl<<"Number of events input = "<<nentries;
cout<<endl<<"Elapsed time = "<<TimeSum<<" seconds";
cout<<endl<<"Number of events timed = "<<EventsIn;
cout<<endl<<"Number of events not timed = "<<EventsOut;
cout<<endl<<"Time of events not timed = "<<TimeOut<<" seconds"<<endl;

cout<<endl<<"Events at PtMcIlwainB <0  or >=6"<<endl;
cout<<"Time0 = "<<Time0<<" s"<<endl;
cout<<"Counter0 = "<<Counter0<<" entries"<<endl<<endl;

}