View source for Solenoid effect in 2 GeV and up range

<pre>
void Effects_GeV()
{
        struct evt_t
                {
                Int_t event;
                Float_t KE;
                Float_t shift;
                };
        ifstream in;
        evt_t evt;

        double E,delta_phi;
        Int_t nlines=0;

        float Base_Number1,Base_Number2,Base_Number3,Base_Number4,Base_Number5=0;
        float normal1,normal2,normal3,normal4,normal5;
        float neg5T_number1,neg5T_number2,neg5T_number3,neg5T_number4,neg5T_number5;
        float pos5T_number1,pos5T_number2,pos5T_number3,pos5T_number4,pos5T_number5;
        float neg25T_number1,neg25T_number2,neg25T_number3,neg25T_number4,neg25T_number5;
        float pos25T_number1,pos25T_number2,pos25T_number3,pos25T_number4,pos25T_number5;
        float pos6T_number1,pos6T_number2,pos6T_number3,pos6T_number4,pos6T_number5;
        float pos4T_number1,pos4T_number2,pos4T_number3,pos4T_number4,pos4T_number5;
        float avg_0T,avg_pos5T,avg_neg5T,avg_neg25T,avg_pos25T,avg_pos4T,avg_pos6T;
        float total_0T,total_pos5T,total_neg5T,total_pos25T,total_neg25T,total_pos4T,total_pos6T;
        float stderror_0T,stderror_pos5T,stderror_neg5T,stderror_pos25T,stderror_neg25T,stderror_pos4T,stderror_pos6T;

        TFile *f = new TFile("Phi_Shift_counts_GeV.root","RECREATE");
        TTree *tree = new TTree("Phi_Shift_counts_GeV","Phi Shift counts wrt 0T");


//
//Count entries from 0T
//
        in.open("Energy_Phi_Shift_GeV_0T_1.dat");
        tree->Branch("evt",&evt.event,"event/I:KE/F:shift");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
                }
        Base_Number1=nlines;
        normal1=Base_Number1/Base_Number1;
        nlines=0;
        in.close();

        in.open("Energy_Phi_Shift_GeV_0T_2.dat");
        tree->Branch("evt",&evt.event,"event/I:KE/F:shift");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
                }
        Base_Number2=nlines;
        normal2=Base_Number2/Base_Number2;
        nlines=0;
        in.close();

        in.open("Energy_Phi_Shift_GeV_0T_3.dat");

        tree->Branch("evt",&evt.event,"event/I:KE/F:shift");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
                }   
        Base_Number3=nlines;
        normal3=Base_Number3/Base_Number3;
        nlines=0;
        in.close();

        in.open("Energy_Phi_Shift_GeV_0T_4.dat");
        tree->Branch("evt",&evt.event,"event/I:KE/F:shift");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
                }
        Base_Number4=nlines;
        normal4=Base_Number4/Base_Number4;
        nlines=0;
        in.close();
        
        in.open("Energy_Phi_Shift_GeV_0T_5.dat");
        tree->Branch("evt",&evt.event,"event/I:KE/F:shift");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
                }
        Base_Number5=nlines;
        normal5=Base_Number5/Base_Number5;
        nlines=0;
        in.close();
                
                 
//
// Normalize the -5T file
//
        in.open("Energy_Phi_Shift_GeV_neg5T_1.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
                
                }
        neg5T_number1=nlines/Base_Number1;
        nlines=0;
        in.close();
        
        in.open("Energy_Phi_Shift_GeV_neg5T_2.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
                 
                }
        neg5T_number2=nlines/Base_Number2;
        nlines=0;
        in.close(); 
        
        in.open("Energy_Phi_Shift_GeV_neg5T_3.dat");
        while(in.good())
                {  
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
        
                }
        neg5T_number3=nlines/Base_Number3;
        nlines=0;
        in.close();
                 
        in.open("Energy_Phi_Shift_GeV_neg5T_4.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
        
                }
        neg5T_number4=nlines/Base_Number4;
        nlines=0;
        in.close();
                
        in.open("Energy_Phi_Shift_GeV_neg5T_5.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
                 
                }
        neg5T_number5=nlines/Base_Number5;
        nlines=0;
        in.close();
        
                 
//
// Normalize the 5T file
//
        in.open("Energy_Phi_Shift_GeV_pos5T_1.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
        
                }
        pos5T_number1=nlines/Base_Number1;
        nlines=0;
        in.close();
                
        in.open("Energy_Phi_Shift_GeV_pos5T_2.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
        
                }
        pos5T_number2=nlines/Base_Number2;
        nlines=0;
        in.close();
                
        
        in.open("Energy_Phi_Shift_GeV_pos5T_3.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
        
                }
        pos5T_number3=nlines/Base_Number3;
        nlines=0;
        in.close();
        
        in.open("Energy_Phi_Shift_GeV_pos5T_4.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
        
                }
        pos5T_number4=nlines/Base_Number4;
        nlines=0;
        in.close();
                 
        in.open("Energy_Phi_Shift_GeV_pos5T_5.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
  
                }
        pos5T_number5=nlines/Base_Number5;
        nlines=0;
        in.close();
                 
                
//
// Normalize the -2.5T file
//      
        in.open("Energy_Phi_Shift_GeV_neg25T_1.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
        
                }
        neg25T_number1=nlines/Base_Number1;
        nlines=0;
        in.close();
        
        in.open("Energy_Phi_Shift_GeV_neg25T_2.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
        
                }
        neg25T_number2=nlines/Base_Number2;
        nlines=0;
        in.close();
                
        in.open("Energy_Phi_Shift_GeV_neg25T_3.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
        
                }
        neg25T_number3=nlines/Base_Number3;
        nlines=0;
        in.close();
                
        in.open("Energy_Phi_Shift_GeV_neg25T_4.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
        
                }
        neg25T_number4=nlines/Base_Number4;
        nlines=0;
        in.close();
                
        in.open("Energy_Phi_Shift_GeV_neg25T_5.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
                
                }
        neg25T_number5=nlines/Base_Number5;
        nlines=0;
        in.close();
        
                 
//
// Normalize the 2.5T file
//
        in.open("Energy_Phi_Shift_GeV_pos25T_1.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
        
                }
        pos25T_number1=nlines/Base_Number1;
        nlines=0;
                
        in.close();
        
        in.open("Energy_Phi_Shift_GeV_pos25T_2.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
        
                }
        pos25T_number2=nlines/Base_Number2;
        nlines=0;
                
        in.close();
                 
        in.open("Energy_Phi_Shift_GeV_pos25T_3.dat");
        while(in.good())
                {  
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
                 
                }
        pos25T_number3=nlines/Base_Number3;
        nlines=0;
        
        in.close();
        
        in.open("Energy_Phi_Shift_GeV_pos25T_4.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
                
                }
        pos25T_number4=nlines/Base_Number4;
        nlines=0;
                 
        in.close();
        
        in.open("Energy_Phi_Shift_GeV_pos25T_5.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
        
                }
        pos25T_number5=nlines/Base_Number5;
        nlines=0;
                
        in.close();
        
                 
//
// Normalize the 6T file
//              
        in.open("Energy_Phi_Shift_GeV_pos6T_1.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
                
                }
        pos6T_number1=nlines/Base_Number1;
        nlines=0;
        in.close();
        
        in.open("Energy_Phi_Shift_GeV_pos6T_2.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
                
                }
        pos6T_number2=nlines/Base_Number2;
        nlines=0;
        in.close();
        
        in.open("Energy_Phi_Shift_GeV_pos6T_3.dat");
        while(in.good())
                {  
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
                 
                }
        pos6T_number3=nlines/Base_Number3;
        nlines=0;
        in.close();
                 
        in.open("Energy_Phi_Shift_GeV_pos6T_4.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
        
                }
        pos6T_number4=nlines/Base_Number4;
        nlines=0;
        in.close();
        
        in.open("Energy_Phi_Shift_GeV_pos6T_5.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
                 
                }
        pos6T_number5=nlines/Base_Number5;
        nlines=0;
        in.close();
        
                 
//
// Normalize the 4T file
//
        in.open("Energy_Phi_Shift_GeV_pos4T_1.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
        
                }
        pos4T_number1=nlines/Base_Number1;
        nlines=0;
        in.close();
                
        in.open("Energy_Phi_Shift_GeV_pos4T_2.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
        
                }
        pos4T_number2=nlines/Base_Number2;
        nlines=0;
        in.close();
                
        in.open("Energy_Phi_Shift_GeV_pos4T_3.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
        
                }
        pos4T_number3=nlines/Base_Number3;
        nlines=0;
        in.close();
                
        in.open("Energy_Phi_Shift_GeV_pos4T_4.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
        
                }
        pos4T_number4=nlines/Base_Number4;
        nlines=0;
        in.close();
                
        in.open("Energy_Phi_Shift_GeV_pos4T_5.dat");
        while(in.good())
                {
                evt.event=nlines;
                in >> evt.KE >> evt.shift;
                nlines++;
                 
                }
        pos4T_number5=nlines/Base_Number5;
        nlines=0;
        in.close();
        
                 
//
// Calculate the error bars
//
//      B=0T
        total_0T=normal1+normal2+normal3+normal4+normal5;
        avg_0T=total_0T/5;
        
//      B=5T
        total_pos5T=pos5T_number1+pos5T_number2+pos5T_number3+pos5T_number4+pos5T_number5;
        avg_pos5T=total_pos5T/5;
        
//      B=-5T
        total_neg5T=neg5T_number1+neg5T_number2+neg5T_number3+neg5T_number4+neg5T_number5;
        avg_neg5T=total_neg5T/5;
                
//      B=2.5T
        total_pos25T=pos25T_number1+pos25T_number2+pos25T_number3+pos25T_number4+pos25T_number5;
        avg_pos25T=total_pos25T/5;
        
//      B=-2.5T
        total_neg25T=neg25T_number1+neg25T_number2+neg25T_number3+neg25T_number4+neg25T_number5;
        avg_neg25T=total_neg25T/5;
        
//      B=4T
        total_pos4T=pos4T_number1+pos4T_number2+pos4T_number3+pos4T_number4+pos4T_number5;
        avg_pos4T=total_pos4T/5;
                
//      B=6T
        total_pos6T=pos6T_number1+pos6T_number2+pos6T_number3+pos6T_number4+pos6T_number5;
        avg_pos6T=total_pos6T/5;
        
        stderror_0T=0;
        stderror_0T=(normal1-avg_0T)*(normal1-avg_0T)+stderror_0T;
        stderror_0T=(normal2-avg_0T)*(normal2-avg_0T)+stderror_0T;
        stderror_0T=(normal3-avg_0T)*(normal3-avg_0T)+stderror_0T;
        stderror_0T=(normal4-avg_0T)*(normal4-avg_0T)+stderror_0T;
        stderror_0T=(normal5-avg_0T)*(normal5-avg_0T)+stderror_0T;
        stderror_0T=sqrt(stderror_0T);
        stderror_0T=stderror_0T/4;
        
        stderror_pos5T=0;
        stderror_pos5T=(pos5T_number1-avg_pos5T)*(pos5T_number1-avg_pos5T)+stderror_pos5T;
        stderror_pos5T=(pos5T_number2-avg_pos5T)*(pos5T_number2-avg_pos5T)+stderror_pos5T;
        stderror_pos5T=(pos5T_number3-avg_pos5T)*(pos5T_number3-avg_pos5T)+stderror_pos5T;
        stderror_pos5T=(pos5T_number4-avg_pos5T)*(pos5T_number4-avg_pos5T)+stderror_pos5T;
        stderror_pos5T=(pos5T_number5-avg_pos5T)*(pos5T_number5-avg_pos5T)+stderror_pos5T;
        stderror_pos5T=sqrt(stderror_pos5T);
        stderror_pos5T=stderror_pos5T/4;
                
        stderror_neg5T=0;
        stderror_neg5T=(neg5T_number1-avg_neg5T)*(neg5T_number1-avg_neg5T)+stderror_neg25T;
        stderror_neg5T=(neg5T_number1-avg_neg5T)*(neg5T_number1-avg_neg5T)+stderror_neg25T;
        stderror_neg5T=(neg5T_number1-avg_neg5T)*(neg5T_number1-avg_neg5T)+stderror_neg25T;
        stderror_neg5T=(neg5T_number1-avg_neg5T)*(neg5T_number1-avg_neg5T)+stderror_neg25T;
        stderror_neg5T=(neg5T_number1-avg_neg5T)*(neg5T_number1-avg_neg5T)+stderror_neg25T;
        stderror_neg5T=sqrt(stderror_neg5T);
        stderror_neg5T=stderror_neg5T/4;
                 
        stderror_neg25T=0;
        stderror_neg25T=(neg25T_number1-avg_neg25T)*(neg25T_number1-avg_neg25T)+stderror_neg25T;
        stderror_neg25T=(neg25T_number2-avg_neg25T)*(neg25T_number2-avg_neg25T)+stderror_neg25T;
        stderror_neg25T=(neg25T_number3-avg_neg25T)*(neg25T_number3-avg_neg25T)+stderror_neg25T;
        stderror_neg25T=(neg25T_number4-avg_neg25T)*(neg25T_number4-avg_neg25T)+stderror_neg25T;
        stderror_neg25T=(neg25T_number5-avg_neg25T)*(neg25T_number5-avg_neg25T)+stderror_neg25T;
        stderror_neg25T=sqrt(stderror_neg25T);
        stderror_neg25T=stderror_neg25T/4;
        
        stderror_pos25T=0;
        stderror_pos25T=(pos25T_number1-avg_pos25T)*(pos25T_number1-avg_pos25T)+stderror_pos25T;
        stderror_pos25T=(pos25T_number2-avg_pos25T)*(pos25T_number2-avg_pos25T)+stderror_pos25T;
        stderror_pos25T=(pos25T_number3-avg_pos25T)*(pos25T_number3-avg_pos25T)+stderror_pos25T;
        stderror_pos25T=(pos25T_number4-avg_pos25T)*(pos25T_number4-avg_pos25T)+stderror_pos25T;
        stderror_pos25T=(pos25T_number5-avg_pos25T)*(pos25T_number5-avg_pos25T)+stderror_pos25T;
        stderror_pos25T=sqrt(stderror_pos25T);
        stderror_pos25T=stderror_pos25T/4;
        
        stderror_pos4T=0;
        stderror_pos4T=(pos4T_number1-avg_pos4T)*(pos4T_number1-avg_pos4T)+stderror_pos4T;
        stderror_pos4T=(pos4T_number2-avg_pos4T)*(pos4T_number2-avg_pos4T)+stderror_pos4T;
        stderror_pos4T=(pos4T_number3-avg_pos4T)*(pos4T_number3-avg_pos4T)+stderror_pos4T;
        stderror_pos4T=(pos4T_number4-avg_pos4T)*(pos4T_number4-avg_pos4T)+stderror_pos4T;
        stderror_pos4T=(pos4T_number5-avg_pos4T)*(pos4T_number5-avg_pos4T)+stderror_pos4T;
        stderror_pos4T=sqrt(stderror_pos4T);
        stderror_pos4T=stderror_pos4T/4;
                
        stderror_pos6T=0;
        stderror_pos6T=(pos6T_number1-avg_pos6T)*(pos6T_number1-avg_pos6T)+stderror_pos6T;
        stderror_pos6T=(pos6T_number2-avg_pos6T)*(pos6T_number2-avg_pos6T)+stderror_pos6T;
        stderror_pos6T=(pos6T_number3-avg_pos6T)*(pos6T_number3-avg_pos6T)+stderror_pos6T;
        stderror_pos6T=(pos6T_number4-avg_pos6T)*(pos6T_number4-avg_pos6T)+stderror_pos6T;
        stderror_pos6T=(pos6T_number5-avg_pos6T)*(pos6T_number5-avg_pos6T)+stderror_pos6T;
        stderror_pos6T=sqrt(stderror_pos6T);
        stderror_pos6T=stderror_pos6T/4;
       
// Make an xmgrace file
        
        fprintf(stdout,"%d %g %g\n",-5,avg_neg5T,stderror_neg5T);
        fprintf(stdout,"%d %g %g\n",-2.5,avg_neg25T,stderror_neg25T);
        fprintf(stdout,"%d %g %g\n",0,avg_0T,stderror_0T);
        fprintf(stdout,"%d %g %g\n",2.5,avg_pos25T,stderror_pos25T);
        fprintf(stdout,"%d %g %g\n",4,avg_pos4T,stderror_pos4T);
        fprintf(stdout,"%d %g %g\n",5,avg_pos5T,stderror_pos5T);
        fprintf(stdout,"%d %g %g\n",6,avg_pos6T,stderror_pos6T);
        
        
        tree->Print();
        tree->Write();   
        f->Write();
        
        delete tree;
        delete f;
        
}
</pre>


----

=Links=

[[Looking_at_effects_of_Solenoid_on_Phi_Shifts|Back]]