Vanwdani: Created page with "
 #include "Riostream.h" #include "TFile.h" #include "TTree.h" #include "TCanvas.h" #include  #include  #include  #include "TLatex.h" #include 
2018-01-11T17:11:49Z
Created page with "<pre> #include "Riostream.h" #include "TFile.h" #include "TTree.h" #include "TCanvas.h" #include <TH1.h> #include <TH2.h> #include <TString.h> #include "TLatex.h" #include <TM..."
New page
<pre>

#include "Riostream.h"

#include "TFile.h"

#include "TTree.h"

#include "TCanvas.h"

#include <TH1.h>

#include <TH2.h>

#include <TString.h>

#include "TLatex.h"

#include <TMath.h>

#include <math.h> 

// In order to run this code as a "macro" in Root, the  main function must

// have the same name as the source file



void plotter_OUT()

{

        gROOT->Reset();

        gStyle->SetStatY(0.75);

        gStyle->SetStatX(0.75);

        gStyle->SetStatW(0.1);

        gStyle->SetStatH(0.12);

        //

        // Initialize input filestream and output file pointers 

        //

        int SolenoidField=0;

        double TorusField;

        char ShieldStatus[10]="Out";

        int Step_Number;

        int i;

        for(i=0;i<2;i++)                //loop over Torus Field Variables

        {

                Step_Number=i+1;

                TorusField=0.5*Step_Number-1.5;

                printf("ShieldStatus=%s\n",ShieldStatus);

                printf("Step_Number=%d\n",Step_Number);

                char root_name[200];

                char conditions[50];

                double temp;

                printf("TorusField=%.1f\n",TorusField);

                temp=TorusField;

                        printf("temp=%f\n",temp);

                temp=fmod(temp,1.0);

                        printf("Modulus=%f\n",temp);



                if(temp==0.0)

                {

                        sprintf(root_name,"/home/lds/src/CLAS/GEMC/experiments/eg12/MolrBckGrd/DV/Isotropic_study/LH2_%dSol_%dTor_11GeV_IsotropicPhi_Shield%s_v2_6/LH2_%dSol_%dTor_11GeV_IsotropicPhi_Shield%s.root",SolenoidField,Tor$

                        sprintf(conditions,"%dxSol\t %dxTor\t Shield%s",SolenoidField,TorusField,ShieldStatus);

                        

                }

                if(temp!=0.0)

                {

                        sprintf(root_name,"/home/lds/src/CLAS/GEMC/experiments/eg12/MolrBckGrd/DV/Isotropic_study/LH2_%dSol_%.1fTor_11GeV_IsotropicPhi_Shield%s_v2_6/LH2_%dSol_%.1fTor_11GeV_IsotropicPhi_Shield%s.root",SolenoidField$

                        sprintf(conditions,"%dxSol\t %.1fxTor\t Shield%s",SolenoidField,TorusField,ShieldStatus);

                }

                printf("root name=%s\n",root_name);

                TFile *file = new TFile(root_name); //load the input file

                TCanvas *c = new TCanvas; // make a canvas

                c->SetRightMargin(0.18);

                TText *xlabel = new TText();

                        xlabel-> SetNDC();

                        xlabel -> SetTextAngle(0);

                        xlabel -> SetTextSize(0.03);









        //get the histograms, plot to canvas and save to file

                TH2F *h = (TH2F*)file->Get("MolEThetaLabLUND");

                h->GetZaxis()->SetTitleOffset(1.2);

                h->Draw();

                char png_filename[200];

                sprintf(png_filename,"/home/lds/src/CLAS/GEMC/experiments/eg12/MolrBckGrd/DV/Isotropic_study/ShieldOUT_gifs/%s_MolEThetaLabLUND/MolEThetaLabLUND_%d.png",ShieldStatus,Step_Number);

                printf("PNG File=%s\n",png_filename);

                        xlabel -> DrawText(0.6, 0.8,conditions);

                c->Print(png_filename);

                c->Clear();  //clear the canvas for a new histogram



                TH2F *h1 = (TH2F*)file->Get("MolEThetaLabHITS");

                h1->GetZaxis()->SetTitleOffset(1.2);

                h1->Draw();

                char png_filename1[200];

                sprintf(png_filename1,"/home/lds/src/CLAS/GEMC/experiments/eg12/MolrBckGrd/DV/Isotropic_study/ShieldOUT_gifs/%s_MolEThetaLabHITS/MolEThetaLabHITS_%d.png",ShieldStatus,Step_Number);

                        xlabel -> DrawText(0.6, 0.8,conditions);

                c->Print(png_filename1);

                c->Clear();



                TH2F *h2 = (TH2F*)file->Get("MolEThetaLabLUNDWeighted");

                h2->GetZaxis()->SetTitleOffset(1.4);

                h2->Draw();

                char png_filename2[200];

                sprintf(png_filename2,"/home/lds/src/CLAS/GEMC/experiments/eg12/MolrBckGrd/DV/Isotropic_study/ShieldOUT_gifs/%s_MolEThetaLabLUNDWeighted/MolEThetaLabLUNDWeighted_%d.png",ShieldStatus,Step_Number);

                        xlabel -> DrawText(0.6, 0.8,conditions);

                c->Print(png_filename2);

                c->Clear();



                TH2F *h3 = (TH2F*)file->Get("MolEThetaLabHITSWeighted");

                h3->GetZaxis()->SetTitleOffset(1.2);

                h3->Draw();

                char png_filename3[200];

                sprintf(png_filename3,"/home/lds/src/CLAS/GEMC/experiments/eg12/MolrBckGrd/DV/Isotropic_study/ShieldOUT_gifs/%s_MolEThetaLabHITSWeighted/MolEThetaLabHITSWeighted_%d.png",ShieldStatus,Step_Number);

                        xlabel -> DrawText(0.6, 0.8,conditions);

                c->Print(png_filename3);

                c->Clear();



                TH1F *h4 = (TH1F*)file->Get("MolEPhiLabLUND");

                h4->GetZaxis()->SetTitleOffset(1.2);

                h4->Draw();

                char png_filename4[200];

                sprintf(png_filename4,"/home/lds/src/CLAS/GEMC/experiments/eg12/MolrBckGrd/DV/Isotropic_study/ShieldOUT_gifs/%s_MolEPhiLabLUND/MolEPhiLabLUND_%d.png",ShieldStatus,Step_Number);

                        xlabel -> DrawText(0.6, 0.8,conditions);

                c->Print(png_filename4);

                c->Clear();



                TH1F *h5 = (TH1F*)file->Get("MolEPhiLabLUNDWeighted");

                h5->GetZaxis()->SetTitleOffset(1.2);

                h5->Draw();

                char png_filename5[200];

                sprintf(png_filename5,"/home/lds/src/CLAS/GEMC/experiments/eg12/MolrBckGrd/DV/Isotropic_study/ShieldOUT_gifs/%s_MolEPhiLabLUNDWeighted/MolEPhiLabLUNDWeighted_%d.png",ShieldStatus,Step_Number);

                        xlabel -> DrawText(0.6, 0.8,conditions);

                c->Print(png_filename5);

                c->Clear();



                TH1F *h6 = (TH1F*)file->Get("MolEPhiLabHITS");

                h6->GetZaxis()->SetTitleOffset(1.2);

                h6->Draw();

                char png_filename6[200];

                sprintf(png_filename6,"/home/lds/src/CLAS/GEMC/experiments/eg12/MolrBckGrd/DV/Isotropic_study/ShieldOUT_gifs/%s_MolEPhiLabHITS/MolEPhiLabHITS_%d.png",ShieldStatus,Step_Number);

                        xlabel -> DrawText(0.6, 0.8,conditions);

                c->Print(png_filename6);

                c->Clear();



                TH1F *h7 = (TH1F*)file->Get("MolEPhiLabHITSWeighted");

                h7->GetZaxis()->SetTitleOffset(1.2);

                h7->Draw();

                char png_filename7[200];

                sprintf(png_filename7,"/home/lds/src/CLAS/GEMC/experiments/eg12/MolrBckGrd/DV/Isotropic_study/ShieldOUT_gifs/%s_MolEPhiLabHITSWeighted/MolEPhiLabHITSWeighted_%d.png",ShieldStatus,Step_Number);

                        xlabel -> DrawText(0.6, 0.8,conditions);

                c->Print(png_filename7);

                c->Clear();



                file->Close();



        }//End of loop

        return 0;

}



</pre>
Plotter OUT.C - Revision history
