- InitRecon();
-}
-
-void AliRICHRecon::InitRecon()
-{
-
- outputfile = new TFile("Anal.root","RECREATE","My Analysis histos");
- if(kDISPLAY) Display = new TCanvas("Display","RICH Display",0,0,1200,750);
-
- h1_photons = new TH1F("h1_photons","photons",750,0.,0.75);
- h1_photacc = new TH1F("h1_photacc","photons",750,0.,0.75);
- h1_hough = new TH1F("h1_hough","hough",750,0.,0.75);
- h1_houghpos= new TH1F("h1_houghpos","hough",750,0.,0.75);
- h1_houghneg= new TH1F("h1_houghneg","hough",750,0.,0.75);
-
- h2_tvsppos = new TH2F("h2_tvsppos","thetac vs p",100,0.,5.,750,0.,0.75);
- h2_tvspneg = new TH2F("h2_tvspneg","thetac vs p",100,0.,5.,750,0.,0.75);
- h2_func = new TH2F("h2_func"," func ",800,0.,0.8,100,-100.,100.);
- h2_mvsp = new TH2F("h2_mvsp","mass vs p",100,0.,5.,200,0.,2.);
- h2_mvst = new TH2F("h2_mvst","mass vs t",750,0.,0.75,200,0.,2.);
- h2_map = new TH2F("h2_map","h2_map",160,0.,160.,96,0.,96.);
- h2_mapw = new TH2F("h2_mapw","h2_mapw",160,0.,160.,96,0.,96.);
-
- h2_dist_p = new TH2F("h2_dist_p","h2_dist_p",100,0.,5.,100,0.,5.);
- //
-
- h2_disp = new TH2F("h2_disp","STAR-RICH Event Display",165,Xmin,Xmax,100,Ymin,Ymax);
-
- // h2_test1 = new TH2F("h2_test1","test1 map",165,-64.,64.,100,-42.,42.);
- h2_test2 = new TH2F("h2_test2","test2 map",165,-64.,64.,100,-42.,42.);
- // h2_test4 = new TH2F("h2_test4","test4 map",165,-64.,64.,100,-42.,42.);
- h2_testmap= new TH2F("h2_testmap","test map",165,-64.,64.,100,-42.,42.);
-
- //
- h1_photons1 = new TH1F("h1_photons1","photons",750,0.,0.75);
- h1_photons2 = new TH1F("h1_photons2","photons",750,0.,0.75);
- //
- h1_hcs = new TH1F("h1_hcs","hcs",750,0.,750.);
- h1_hcsw = new TH1F("h1_hcsw","hcsw",750,0.,750.);
- //
- h1_nprotons = new TH1F("h1_nprotons","n prot",30,0.,30.);
- //
- hp_1pos = new TProfile("hp_1pos","Nphot vs thetac pos",250,0.,0.75);
- hp_1neg = new TProfile("hp_1neg","Nphot vs thetac neg",250,0.,0.75);
- hp_1posnorm = new TProfile("hp_1posnorm","Nphot vs thetac pos norm",250,0.,0.75);
- hp_1negnorm = new TProfile("hp_1negnorm","Nphot vs thetac neg norm",250,0.,0.75);
- //
- h2_1pos = new TH2F("h2_1pos","Nphot vs p pos",100,0.,5.,30,0.,30.);
- h2_1neg = new TH2F("h2_1neg","Nphot vs p neg",100,0.,5.,30,0.,30.);
- h2_1posnorm = new TH2F("h2_1posnorm","Nphot vs p pos norm",100,0.,5.,30,0.,30.);
- h2_1negnorm = new TH2F("h2_1negnorm","Nphot vs p neg norm",100,0.,5.,30,0.,30.);
-
- h1_deltap = new TH1F("h1_deltap","delta_p",200,-0.5,0.5);
- h1_deltapop = new TH1F("h1_deltapop","deltapop",200,-1.,1.);
- h1_diffTrackTheta = new TH1F("h1_diffTrackTheta","delta theta",200,-0.25,0.25);
- h1_diffTrackPhi = new TH1F("h1_diffTrackPhi","delta phi",200,-0.25,0.25);
-
- h1_photaccspread = new TH1F("h1_photaccspread","photons spread",200,-0.1,0.1);
-
- //
-
- h1_mass = new TH1F("h1_mass","mass",200,0.,2.);
- photris = new TH1F("photris","photris",1000,0.,1.);
- h2_diffneg = new TH2F("h2_diffneg","diff neg",100,-2.5,2.5,100,-2.5,2.5);
- h2_diffpos = new TH2F("h2_diffpos","diff pos",100,-2.5,2.5,100,-2.5,2.5);
- gChargeMipH1 = new TH1F("gChargeMipH1"," Charge Mip ",2000,0.,2000.);
- gMultMipH1 = new TH1F("gMultMipH1"," Cluster Size Mip ",50,0.,50.);
-
-
- hn = new TNtuple("hn","ntuple",