//--------------------------------------------- // Go4 Tasca analysis // Author: Hans G. Essel // H.Essel@gsi.de // GSI, Experiment Electronics, Data Processing // modified by baatar, r.maendl and et al., // July 2021 - cleanup by JAM (j.adamczewski@gsi.de) //--------------------------------------------- #include "TascaAnlProc.h" #include "Riostream.h" #include "TH1.h" #include "TH2.h" #include "TGraph.h" #include "TMultiGraph.h" #include "TGo4WinCond.h" #include "TascaAnlEvent.h" #include "TascaCaliEvent.h" #include "TGo4WinCond.h" #include "TascaAnalysis.h" #include "TascaUnpackEvent.h" //*********************************************************** TascaAnlProc::TascaAnlProc() : TGo4EventProcessor(), fInput(0) { } //*********************************************************** // this one is used in TascaAnlFact TascaAnlProc::TascaAnlProc(const char* name) : TGo4EventProcessor(name) { cout << "Tasca> TascaAnlProc: Create" << endl; //// init user analysis objects: TString obname; TString obtitle; Char_t comments[128]; // this is the correlation input data: setcorr = fopen("set_correlation.dat", "r+"); if (setcorr == 0) { cout << "Tasca> TascaAnlProc: Could not open required input set_correlation.dat. Exiting go4 analysis!" << endl; exit(0); } fscanf(setcorr, "%s \n", comments); printf("%s \n", comments); fscanf(setcorr, "%f %f %d %d %d %d %d \n\n", &ER_STOPXH_Low, &ER_STOPXH_High, &ER_TOFEL_Low, &ER_TOFEL_High, &ER_TOFEH_Low, &ER_MACRO, &ER_MICRO); printf("%6.2f\t %6.2f\t %5d\t %5d\t %5d\t %1d\t %1d \n", ER_STOPXH_Low, ER_STOPXH_High, ER_TOFEL_Low, ER_TOFEL_High, ER_TOFEH_Low, ER_MACRO, ER_MICRO); fscanf(setcorr, "%s \n", comments); printf("%s \n", comments); fscanf(setcorr, "%f %f %d %d %f \n\n", &AL1_STOPXL_Low, &AL1_STOPXL_High, &AL1_TOFEH, &AL1_MICRO, &AL1_Time); printf("%6.2f %6.2f %d %d %6.3f \n", AL1_STOPXL_Low, AL1_STOPXL_High, AL1_TOFEH, AL1_MICRO, AL1_Time); fscanf(setcorr, "%s \n", comments); printf("%s \n", comments); fscanf(setcorr, "%f %f %d %d %f \n\n", &AL2_STOPXL_Low, &AL2_STOPXL_High, &AL2_TOFEH, &AL2_MICRO, &AL2_Time); printf("%6.2f %6.2f %d %d %6.3f \n", AL2_STOPXL_Low, AL2_STOPXL_High, AL2_TOFEH, AL2_MICRO, AL2_Time); fscanf(setcorr, "%s \n", comments); printf("%s \n", comments); fscanf(setcorr, "%f %f %d %d %f \n\n", &AL3_STOPXL_Low, &AL3_STOPXL_High, &AL3_TOFEH, &AL3_MICRO, &AL3_Time); printf("%6.2f %6.2f %d %d %6.3f \n", AL3_STOPXL_Low, AL3_STOPXL_High, AL3_TOFEH, AL3_MICRO, AL3_Time); fscanf(setcorr, "%s \n", comments); printf("%s \n", comments); fscanf(setcorr, "%f %f %d %d %f \n\n", &SF_STOPXL_Low, &SF_STOPXL_High, &SF_TOFEH, &SF_MICRO, &SF_Time); printf("%6.2f %6.2f %d %d %6.3f \n", SF_STOPXL_Low, SF_STOPXL_High, SF_TOFEH, SF_MICRO, SF_Time); //here starts the correlation activation fscanf(setcorr, "%s \n", comments); fscanf(setcorr, "%s \n", comments); fscanf(setcorr, "%d \n", &ER_AL1); printf("%d \n", ER_AL1); fscanf(setcorr, "%s \n\n", comments); printf("%s \n", comments); fscanf(setcorr, "%s \n", comments); fscanf(setcorr, "%d \n", &ER_AL1_AL2); printf("%d \n", ER_AL1_AL2); fscanf(setcorr, "%s \n\n", comments); printf("%s \n", comments); fscanf(setcorr, "%s \n", comments); fscanf(setcorr, "%d \n\n", &AL1_AL2); printf("%d \n", AL1_AL2); fscanf(setcorr, "%s \n", comments); printf("%s \n", comments); fscanf(setcorr, "%s \n", comments); fscanf(setcorr, "%d \n\n", &ER_AL1_AL2_AL3); printf("%d \n", ER_AL1_AL2_AL3); fscanf(setcorr, "%s \n", comments); printf("%s \n", comments); fscanf(setcorr, "%s \n", comments); fscanf(setcorr, "%d \n\n", &AL1_AL2_AL3); printf("%d \n", AL1_AL2_AL3); fscanf(setcorr, "%s \n", comments); printf("%s \n", comments); fscanf(setcorr, "%s \n", comments); fscanf(setcorr, "%d \n\n", &ER_SF); printf("%d \n", ER_SF); fscanf(setcorr, "%s \n", comments); printf("%s \n", comments); fscanf(setcorr, "%s \n", comments); fscanf(setcorr, "%d \n\n", &ER_AL1_SF); printf("%d \n", ER_AL1_SF); fscanf(setcorr, "%s \n", comments); printf("%s \n", comments); fscanf(setcorr, "%s \n", comments); fscanf(setcorr, "%d \n\n", &ER_AL1_AL2_SF); printf("%d \n", ER_AL1_AL2_SF); fscanf(setcorr, "%s \n", comments); printf("%s \n", comments); fscanf(setcorr, "%s \n", comments); fscanf(setcorr, "%d \n\n", &AL1_AL2_SF); printf("%d \n", AL1_AL2_SF); Fissions = fopen("Fissions.dat", "w+"); if (ER_AL1 == 1) { ER_AL_stop = fopen("ER_AL_STOP.dat", "w+"); ER_AL_sb = fopen("ER_AL_SB.dat", "w+"); } if (ER_AL1_AL2 == 1) ER_AL_AL = fopen("ER_AL_AL.dat", "w+"); if (AL1_AL2 == 1) { AL_AL = fopen("AL_AL.dat", "w+"); AL_AL_sb = fopen("AL_AL_SBs.dat", "w+"); } if (ER_AL1_AL2_AL3 == 1) ER_AL_AL_AL = fopen("ER_AL_AL_AL.dat", "w+"); if (AL1_AL2_AL3 == 1) AL_AL_AL = fopen("AL_AL_AL.dat", "w+"); if (ER_SF == 1) ER_FI = fopen("ER_SF.dat", "w+"); if (ER_AL1_SF == 1) ER_AL_SF = fopen("ER_AL_SF.dat", "w+"); if (ER_AL1_AL2_SF == 1) ER_AL_AL_SF = fopen("ER_AL_AL_SF.dat", "w+"); if (AL1_AL2_SF == 1) AL_AL_SF = fopen("AL_AL_SF.dat", "w+"); ERcount = 0; AL1count = 0; AL2count = 0; AL3count = 0; deltaRe_sf_time = 0.; deltaal_sf_time = 30000.; Answers_TX_TX = 0; Answers_TX_BX = 0; Answers_BX_TX = 0; Answers_BX_BX = 0; Com2_up_up = 0; Com2_up_do = 0; Com2_do_up = 0; Com2_do_do = 0; // Com3_up_up=0; // Com3_up_do=0; // Com3_do_up=0; // Com3_do_do=0; 7 // fprintf(Output,"Time1 (s) \t"); // fprintf(Output,"Time2 (us) \t"); // fprintf(Output,"TMP (us) \t"); // fprintf(Output,"Pulse(0) \t"); // fprintf(Output,"Low energy XL \t"); // fprintf(Output,"Front strip (XL) \t"); // fprintf(Output,"High energy XH \t"); // fprintf(Output,"Front strip (XH) \t"); // fprintf(Output,"Back Strip(Y) \t"); // fprintf(Output,"Low energy YL \t"); // fprintf(Output,"Back strip (YL) \t"); // fprintf(Output,"High energy YH \t"); // fprintf(Output,"Back strip (YH) \t"); // fprintf(Output,"TOF_EL \t"); // fprintf(Output,"TOF_EH \t"); // fprintf(Output,"Target \n"); // fhdStopYL_NoSplit=MakeTH1('D',"Anl/NoSplit", "StopX sum low",4000,300,40300); // fhdStopYL_Single=MakeTH1('D',"Anl/Single", "StopX sum low",4000,300,40300); // fhdCorrelation_Time=MakeTH1('D',"Anl/CorrTime", "Correlation Time",100,0,0.1); fhEnergyAnswers_TX = MakeTH1('D', "Anl/Energy_Answers_TX_all", "Energy of compact1", 1600, 1, 16001); fhEnergyAnswers_TX->GetXaxis()->SetTitle("Energy [keV]"); fhEnergyAnswers_TX_pa = MakeTH1('D', "Anl/Energy_Answers_TX_all_pause", "Energy of compact1 pause", 1600, 1, 16001); fhEnergyAnswers_TX_pa->GetXaxis()->SetTitle("Energy [keV]"); fhEnergyAnswers_TX_single = MakeTH1('D', "Anl/Energy_Answers_TX_all_single", "Energy of compact1 single", 1600, 1, 16001); fhEnergyAnswers_TX_single->GetXaxis()->SetTitle("Energy [keV]"); fhEnergyAnswers_BX = MakeTH1('D', "Anl/Energy_Answers_BX_all", "Energy of compact1", 1600, 1, 16001); fhEnergyAnswers_BX->GetXaxis()->SetTitle("Energy [keV]"); fhEnergyAnswers_BX_pa = MakeTH1('D', "Anl/Energy_Answers_BX_all_pause", "Energy of compact1 pause", 1600, 1, 16001); fhEnergyAnswers_BX_pa->GetXaxis()->SetTitle("Energy [keV]"); fhEnergyAnswers_BX_single = MakeTH1('D', "Anl/Energy_Answers_BX_all_single", "Energy of compact1 single", 1600, 1, 16001); fhEnergyAnswers_BX_single->GetXaxis()->SetTitle("Energy [keV]"); fhGamma0 = MakeTH1('D', "Anl/Gamma1_Energy", "Energy of gamma1", 6000, 1, 3001); fhGamma0->GetXaxis()->SetTitle("Energy [keV]"); fhGamma1 = MakeTH1('D', "Anl/Gamma2_Energy", "Energy of gamma2", 6000, 1, 3001); fhGamma1->GetXaxis()->SetTitle("Energy [keV]"); fhGammaMysec0 = MakeTH1('D', "Anl/Gamma1_Time", "hit time of gamma1", 3000, -1000, 2000); fhGammaMysec0->GetXaxis()->SetTitle("Time [10 ns]"); fhGammaMysec1 = MakeTH1('D', "Anl/Gamma2_Time", "hit time of gamma2", 3000, -1000, 2000); fhGammaMysec1->GetXaxis()->SetTitle("Time [10 ns]"); fhTime_Com2_up = MakeTH1('D', "Anl/Distr/Time_Com2_up", "Distibution Time", 2000, 0, 49); fhTime_Com2_up->GetXaxis()->SetTitle("Time [s]"); fhTime_Com2_down = MakeTH1('D', "Anl/Distr/Time_Com2_down", "Distibution Time", 2000, 0, 49); fhTime_Com2_down->GetXaxis()->SetTitle("Time [s]"); fhTime_FPD_fr = MakeTH1('D', "Anl/Distr/Time_FPD_fr", "Distibution Time", 100, 0, 49); fhTime_FPD_fr->GetXaxis()->SetTitle("Time [s]"); fhTime_FPD_ba = MakeTH1('D', "Anl/Distr/Time_FPD_ba", "Distibution Time", 100, 0, 49); fhTime_FPD_ba->GetXaxis()->SetTitle("Time [s]"); fhAfterFi_Si = MakeTH1('D', "Anl/Events_after_al_or_Fi_Up", "Energy of events after ", 6000, 1, 3001); fhAfterFi_Si->GetXaxis()->SetTitle("Energy [keV]"); fhAfterFi_Ga = MakeTH1('D', "Anl/Events_after_al_or_Fi_Gamma", "Energy of events after ", 6000, 1, 3001); fhAfterFi_Ga->GetXaxis()->SetTitle("Energy [keV]"); fhdUP1fpga = MakeTH1('D', "Anl/Distr/UP1fpga", "Up 1 selected", 32, -0.5, 31.5); fhdDO1fpga = MakeTH1('D', "Anl/Distr/DO1fpga", "Down 1 selected", 32, -0.5, 31.5); fhdUP2fpga = MakeTH1('D', "Anl/Distr/UP2fpga", "Up 2 selected", 32, -0.5, 31.5); fhdDO2fpga = MakeTH1('D', "Anl/Distr/DO2fpga", "Down 2 selected", 32, -0.5, 31.5); fhdFPD_fr_fpga = MakeTH1('D', "Anl/Distr/FPD_fr_fpga", "FPD fr selected", 32, -0.5, 31.5); fhdFPD_ba_fpga = MakeTH1('D', "Anl/Distr/FPD_ba_fpga", "FPD ba selected", 32, -0.5, 31.5); fhFPD_fr_ba_fpga = MakeTH2('D', "Anl/Distr/FPD_fr_ba_fpga", "FPD fr ba selected", 25, -1, 24, 49, -1, 48); fhFPD_fr_ba_fpga->GetYaxis()->SetTitle("Front YL"); fhFPD_fr_ba_fpga->GetXaxis()->SetTitle("Back XL"); // Creation of conditions (check if restored from auto save file): // ffpgaKev = MakeWinCond("Anl/fpgaKevL",AL3_STOPXL_Low ,AL3_STOPXL_High); // ffpgaKev->Enable(); // ffpgaKev->ResetCounts(); // /////////////////////////// // fBeamSt = MakeWinCond("Anl/BeamSt",-1 ,2); // fBeamSt->Enable(); // fBeamSt->ResetCounts(); //limits -1, 1 beam OFF; limits 0,2 beam ON //////////////////////// // fStopXY_Low=MakeTH2('D',"Anl/StopXYhits_Low","Hit counters",144,0,144,50,-1,49); // fStopXY_Low->GetXaxis()->SetTitle("X position [stripe]"); // fStopXY_Low->GetYaxis()->SetTitle("Y position [stripe]"); // fStopXY_Low->GetZaxis()->SetTitle("Hits"); // fStopXY_High=MakeTH2('D',"Anl/StopXYhits_High","Hit counters",144,0,144,50,-1,49); // fStopXY_High->GetXaxis()->SetTitle("X position [stripe]"); // fStopXY_High->GetYaxis()->SetTitle("Y position [stripe]"); // fStopXY_High->GetZaxis()->SetTitle("Hits"); // for(i=0;i<48;i++){ // obname.Form("Anl/StopHE/XH_%03d",i); // obtitle.Form("Stop X High %03d",i); // fStopHE[i]=MakeTH2('D',obname.Data(),obtitle.Data(),144,0,144,200,0,300000); // fStopHE[i]->GetXaxis()->SetTitle("X position [stripe]"); // fStopHE[i]->GetYaxis()->SetTitle("Energy [Kev]"); // obname.Form("Anl/StopLE/XL_%03d",i); // obtitle.Form("Stop X Low %03d",i); // fStopLE[i]=MakeTH2('D',obname.Data(),obtitle.Data(),144,0,144,200,0,30000); // fStopLE[i]->GetXaxis()->SetTitle("X position [stripe]"); // fStopLE[i]->GetYaxis()->SetTitle("Energy [Kev]"); // } // // fStopXHsumvsTAC=MakeTH2('D',"Anl/StopXHvsTAC","STOPXH with TAC",2000,0,2000,200,0,300000); // fStopXHsumvsTAC->GetXaxis()->SetTitle("TAC [ns]"); // fStopXHsumvsTAC->GetYaxis()->SetTitle("Energy [keV]"); // fCom1_Com2 = MakeTH2('D', "Anl/Fissions", "FF energies ", 150, 0, 150000, 150, 0, 150000); fCom1_Com2->GetXaxis()->SetTitle("Energy_FF2_Com1 [keV]"); fCom1_Com2->GetYaxis()->SetTitle("Energy_FF2_Com2 [keV]"); fAL_Ga = MakeTH2('D', "Anl/AL_vs_Gamma", "E_both_AL_Ga_Com", 1000, 0, 10000, 2000, 0, 2000); fAL_Ga->GetXaxis()->SetTitle("Energy_AL [keV]"); fAL_Ga->GetYaxis()->SetTitle("Energy_Ga [keV]"); fAL_Ga_Answers_TX = MakeTH2('D', "Anl/AL_vs_Gamma_Answers_TX", "E_AL_Ga_Answers_TX", 4000, 0, 4000, 2000, 0, 2000); fAL_Ga_Answers_TX->GetXaxis()->SetTitle("Energy_AL [keV]"); fAL_Ga_Answers_TX->GetYaxis()->SetTitle("Energy_Ga [keV]"); fAL_Ga_Answers_BX = MakeTH2('D', "Anl/AL_vs_Gamma_Answers_BX", "E_AL_Ga_Answers_BX", 1000, 0, 10000, 2000, 0, 2000); fAL_Ga_Answers_BX->GetXaxis()->SetTitle("Energy_AL [keV]"); fAL_Ga_Answers_BX->GetYaxis()->SetTitle("Energy_Ga [keV]"); fAL_Ga_Answers_TX_1mm = MakeTH2('D', "Anl/AL_vs_Gamma_Answers_TX_1mm", "E_AL_Ga_Answers_TX_1mm", 1000, 0, 10000, 2000, 0, 2000); fAL_Ga_Answers_TX_1mm->GetXaxis()->SetTitle("Energy_AL [keV]"); fAL_Ga_Answers_TX_1mm->GetYaxis()->SetTitle("Energy_Ga [keV]"); fAL_Ga_Answers_TX_2mm = MakeTH2('D', "Anl/AL_vs_Gamma_Answers_TX_2mm", "E_AL_Ga_Answers_TX_2mm", 1000, 0, 10000, 2000, 0, 2000); fAL_Ga_Answers_TX_2mm->GetXaxis()->SetTitle("Energy_AL [keV]"); fAL_Ga_Answers_TX_2mm->GetYaxis()->SetTitle("Energy_Ga [keV]"); fAL_Ga_Answers_TX_5mm = MakeTH2('D', "Anl/AL_vs_Gamma_Answers_TX_5mm", "E_AL_Ga_Answers_TX_5mm", 1000, 0, 10000, 2000, 0, 2000); fAL_Ga_Answers_TX_5mm->GetXaxis()->SetTitle("Energy_AL [keV]"); fAL_Ga_Answers_TX_5mm->GetYaxis()->SetTitle("Energy_Ga [keV]"); fAL_Ga_Answers_TX_10mm = MakeTH2('D', "Anl/AL_vs_Gamma_Answers_TX_10mm", "E_AL_Ga_Answers_TX_10mm", 1000, 0, 10000, 2000, 0, 2000); fAL_Ga_Answers_TX_10mm->GetXaxis()->SetTitle("Energy_AL [keV]"); fAL_Ga_Answers_TX_10mm->GetYaxis()->SetTitle("Energy_Ga [keV]"); fBe_Ga = MakeTH2('D', "Anl/Gam_vs_Gam_prompt", "Ga_Ga_(<100ns)", 1000, 0, 2000, 1000, 0, 2000); fBe_Ga->GetXaxis()->SetTitle("Energy_Gamma1 [keV]"); fBe_Ga->GetYaxis()->SetTitle("Energy_Gamma2 [keV]"); fBe_AL = MakeTH2('D', "Anl/Beta_vs_Alph", "E_both_Bet_AL_Com", 500, 0, 4000, 1000, 0, 10000); fBe_AL->GetXaxis()->SetTitle("Energy_Beta [keV]"); fBe_AL->GetYaxis()->SetTitle("Energy_Alpha [keV]"); fBe_AL_Ga = MakeTH2('D', "Anl/Gam_vs_Alph_Alp", "Gamma_with_Bet_AL_Coms", 500, 0, 4000, 2000, 0, 2000); fBe_AL_Ga->GetXaxis()->SetTitle("Energy_Beta [keV]"); fBe_AL_Ga->GetYaxis()->SetTitle("Energy_Gamma [keV]"); fGa_Ga = MakeTH2('D', "Anl/Gam_vs_Gam_delayed", "Ga_Ga_(0.1-1us)", 1000, 0, 2000, 1000, 0, 2000); fGa_Ga->GetXaxis()->SetTitle("Energy_Gamma1 [keV]"); fGa_Ga->GetYaxis()->SetTitle("Energy_Gamma2 [keV]"); // ALBEGA SETTING BEGIN // fAL_Ga=MakeTH2('D',"Anl/AL_vs_Gamma","E_both_AL_Ga_Com",1000,0,10000,2000,0,2000); // fAL_Ga->GetXaxis()->SetTitle("Energy_AL [keV]"); // fAL_Ga->GetYaxis()->SetTitle("Energy_Ga [keV]"); // fBe_Ga=MakeTH2('D',"Anl/Beta_vs_Gamma","E_both_Bet_Ga_Com",500,0,4000,2000,0,2000); // fBe_Ga->GetXaxis()->SetTitle("Energy_Beta [keV]"); // fBe_Ga->GetYaxis()->SetTitle("Energy_Ga [keV]"); // fBe_AL=MakeTH2('D',"Anl/Beta_vs_Alph","E_both_Bet_AL_Com",500,0,4000,1000,0,10000); // fBe_AL->GetXaxis()->SetTitle("Energy_Beta [keV]"); // fBe_AL->GetYaxis()->SetTitle("Energy_Alpha [keV]"); // fBe_AL_Ga=MakeTH2('D',"Anl/Gam_vs_Alph_Alp","Gamma_with_Bet_AL_Coms",500,0,4000,2000,0,2000); // fBe_AL_Ga->GetXaxis()->SetTitle("Energy_Beta [keV]"); // fBe_AL_Ga->GetYaxis()->SetTitle("Energy_Gamma [keV]"); // fGa_Ga=MakeTH2('D',"Anl/Gam_vs_Gam","Gamma_1_2",2000,0,2000,2000,0,2000); // fGa_Ga->GetXaxis()->SetTitle("Energy_Gamma1 [keV]"); // fGa_Ga->GetYaxis()->SetTitle("Energy_Gamma2 [keV]"); // END ALBEGA SETTING // fStopXHvsTOF=MakeTH2('D',"Anl/StopXHantiTOF","STOPXH with TOF",200,0,300000,1000,0,5000); // fStopXHvsTOF->GetXaxis()->SetTitle("Energy [MEV]"); // fStopXHvsTOF->GetYaxis()->SetTitle("TOF [chan]"); // // fStopXLvsGamma=MakeTH2('D',"Anl/StopXLvsGamma","STOPXL with gamma",2000,400,20000,1300,0,1300); // fStopXLvsGamma->GetXaxis()->SetTitle("Energy alpha [keV]"); // fStopXLvsGamma->GetYaxis()->SetTitle("Energy gamma [kev]"); // // for(i=0;i<144;i++){ // obname.Form("Anl/StopXLvsTACf /XLvsTAC_%03d",i); // obtitle.Form("StopXL_TAC %03d",i); // fStopXLvsTAC[i]=MakeTH2('D',"obname.Data(),obtitle.Data(),2000,0,2000,200,0,30000); // fStopXLvsTAC[i]->GetXaxis()->SetTitle("TAC [ns]"); // fStopXLvsTAC[i]->GetYaxis()->SetTitle("Energy [keV]"); // obname.Form("Anl/StopXHvsTACf/XHvsTAC_%03d",i); // obtitle.Form("StopXH_TAC %03d",i); // fStopXHvsTAC[i]=MakeTH2('D',obname.Data(),obtitle.Data(),2000,0,2000,200,0,300000); // fStopXHvsTAC[i]->GetXaxis()->SetTitle("TAC [ns]"); // fStopXHvsTAC[i]->GetYaxis()->SetTitle("Energy [keV]"); // } /* for(i=0;i<8;i++){ obname.Form("Anl/TACVSEH/TACvsEH_%03d",i); obtitle.Form("TAC vs Energy High %03d",i); fTACvsEH[i]=MakeTH2('D',obname.Data(),obtitle.Data(),4000,0,4000,5000,0,5000); obname.Form("Anl/TACVSEL/TACvsEL_%03d",i); obtitle.Form("TAC vs Energy Low %03d",i); fTACvsEL[i]=MakeTH2('D',"obname.Data(),obtitle.Data(),4000,0,4000,5000,0, 5000); } */ } //*********************************************************** TascaAnlProc::~TascaAnlProc() { cout << "Tasca> TascaAnlProc: Delete" << endl; Int_t end = 10000000; fclose(setcorr); if (ER_AL1_AL2 == 1) { fprintf(ER_AL_AL, "%d \n", end); fclose(ER_AL_AL); } if (ER_AL1 == 1) { fprintf(ER_AL_stop, "%d \n", end); fclose(ER_AL_stop); fprintf(ER_AL_sb, "%d \n", end); fclose(ER_AL_sb); } if (AL1_AL2 == 1) { fprintf(AL_AL, "%d \n", end); fclose(AL_AL); fprintf(AL_AL_sb, "%d \n", end); fclose(AL_AL_sb); } if (ER_AL1_AL2_AL3 == 1) { fprintf(ER_AL_AL_AL, "%d \n", end); fclose(ER_AL_AL_AL); } if (AL1_AL2_AL3 == 1) { fprintf(AL_AL_AL, "%d \n", end); fclose(AL_AL_AL); } if (AL1_AL2_SF == 1) { fprintf(AL_AL_SF, "%d \n", end); fclose(AL_AL_SF); } if (ER_AL1_SF == 1) { fprintf(ER_AL_SF, "%d \n", end); fclose(ER_AL_SF); } if (ER_SF == 1) { fprintf(ER_FI, "%d \n", end); fclose(ER_FI); } fprintf(Fissions, "%d \n", end); fclose(Fissions); } //cout<<"ER's="<fiSystemSec * 0.00000001)); } ffTY_SF_sum = 0.; ffBY_SF_sum = 0.; // Bool_t Beam_Status=fBeamSt->Test(fInput->fiMicro); // prompt gamma_gamma if (fInput->fiMicro == 0 && (-30 < fInput->fiGammaMysec[0]) && (fInput->fiGammaMysec[0] < -20) && fInput->ffGammaKev[0] > 0. && fInput->ffGammaKev[1] > 0.) { if ((-30 < fInput->fiGammaMysec[1]) && (fInput->fiGammaMysec[1] < -20)) fBe_Ga->Fill(fInput->ffGammaKev[0], fInput->ffGammaKev[1]); } if (fInput->fiMicro == 0 && (-30 < fInput->fiGammaMysec[1]) && (fInput->fiGammaMysec[1] < -20) && fInput->ffGammaKev[1] > 0. && fInput->ffGammaKev[0] > 0.) { if ((-30 < fInput->fiGammaMysec[0]) && (fInput->fiGammaMysec[0] < -20)) fBe_Ga->Fill(fInput->ffGammaKev[1], fInput->ffGammaKev[0]); } if (fInput->fiMicro == 0 && fInput->ffGammaKev[0] > 0.) { fhGamma0->Fill(fInput->ffGammaKev[0]); } if (fInput->fiMicro == 0 && fInput->ffGammaKev[1] > 0.) { fhGamma1->Fill(fInput->ffGammaKev[1]); } // if(fInput->fiMicro==1 && fInput->ffGammaKev[0]>0.) {fhGamma0->Fill(fInput->ffGammaKev[0]);} // if(fInput->fiMicro==1 && fInput->ffGammaKev[1]>0.) {fhGamma1->Fill(fInput->ffGammaKev[1]);} //if(fInput->fiMicro==0) fhGamma0->Fill(fInput->ffGammaKev[0]); //if(fInput->fiMicro==0) fhGamma1->Fill(fInput->ffGammaKev[1]); // if(fInput->ffGammaKev[0]>0. || fInput->ffGammaKev[1]>0.){ // fprintf(Fissions,"Gammas \t %6.6f\t %6.2f\t %d\t %6.2f\t %d\n", // (Float_t) (fInput->fiSystemSec*0.00000001), // fInput->ffGammaKev[0], // fInput->fiGammaMysec[0], // fInput->ffGammaKev[1], // fInput->fiGammaMysec[1]); // } fhGammaMysec0->Fill(fInput->fiGammaMysec[0]); fhGammaMysec1->Fill(fInput->fiGammaMysec[1]); //delayed gamma_gamma if ((-30 < fInput->fiGammaMysec[0]) && (fInput->fiGammaMysec[0] < -20) && fInput->ffGammaKev[0] > 0. && fInput->ffGammaKev[1] > 0.) { if ((-20 < fInput->fiGammaMysec[1]) && (fInput->fiGammaMysec[1] < 100)) fGa_Ga->Fill(fInput->ffGammaKev[0], fInput->ffGammaKev[1]); } if ((-30 < fInput->fiGammaMysec[1]) && (fInput->fiGammaMysec[1] < -20) && fInput->ffGammaKev[1] > 0. && fInput->ffGammaKev[0] > 0.) { if ((-20 < fInput->fiGammaMysec[0]) && (fInput->fiGammaMysec[0] < 100)) fGa_Ga->Fill(fInput->ffGammaKev[1], fInput->ffGammaKev[0]); } // if(fInput->fiChopper==1){ // poutevt->fiaChopperSystemSec=fInput->fiChopper_SystemSec; // printf(" chopper on at: %10.2f\n ",(Float_t)(fInput->fiChopper_SystemSec*0.00000001)); // } for (i = 0; i < 24; i++) { // if(fInput->ffTX[i]>AL1_STOPXL_Low && fInput->ffBX[i]>20000.){ // fprintf(Fissions," %10.8f\t %6.2f\t %6.2f\t %6.2f\t %6.2f\t %d\t %d\t %d\t %6.2f\t %d\t %6.2f\t %d\n", // fInput->fiSystemSec*0.00000001, // fInput->ffTX[i], // fInput->ffTY[i], // fInput->ffBX[i], // fInput->ffBY[i], // i, // fInput->fiMicro, // fInput->fievcount, // fInput->ffGammaKev[0], // fInput->fiGammaMysec[0], // fInput->ffGammaKev[1], // fInput->fiGammaMysec[1]); // } // if (fInput->ffTX[i] > 0.) fhEnergyAnswers_TX->Fill(fInput->ffTX[i]); if (fInput->ffBX[i] > 0.) fhEnergyAnswers_BX->Fill(fInput->ffBX[i]); if (fInput->fiMicro == 0) { if (fInput->ffTX[i] > 0.) { if ((-40 < fInput->fiGammaMysec[0]) && (fInput->fiGammaMysec[0] < 100) && fInput->ffGammaKev[0] > 0.) fAL_Ga->Fill(fInput->ffTX[i], fInput->ffGammaKev[0]); if ((-40 < fInput->fiGammaMysec[1]) && (fInput->fiGammaMysec[1] < 100) && fInput->ffGammaKev[1] > 0.) fAL_Ga->Fill(fInput->ffTX[i], fInput->ffGammaKev[1]); fhEnergyAnswers_TX_pa->Fill(fInput->ffTX[i]); if (fInput->fiMicro == 0 && fInput->ffTX[i] > AL1_STOPXL_Low && fInput->ffTX[i] < AL1_STOPXL_High) { fprintf(Fissions, "pUnpackEvent->fievent_number==%d ||\n", fInput->fievcount); fprintf(ER_AL_stop, "((l_feb_id==0 || l_feb_id==2) && l_cha_id==%d && pUnpackEvent->fievent_number==%d) ||\n", i, fInput->fievcount); } Int_t one1 = 2; for (k = 0; k < 48; k++) { if (one1 != 1 && fInput->ffTY[k] > 5000. && fInput->ffTY[k] < (AL1_STOPXL_High + (500.)) && fInput->ffTX[i] > AL1_STOPXL_Low && fInput->ffTX[i] < AL1_STOPXL_High) { //fhFPD_fr_ba_fpga->Fill(i,k); { fprintf(ER_AL_sb, "Top %10.8f\t %d\t %6.2f\t %d\t %6.2f\t %d\t %d\t %6.2f\t %d\t %6.2f\t %d\n", fInput->fiSystemSec * 0.00000001, i, fInput->ffTX[i], k, fInput->ffTY[k], fInput->fiMicro, fInput->fievcount, fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1]); one1 = 1; } } if (fInput->ffTY[k] > SF_STOPXL_Low) { ffTY_SF_sum = ffTY_SF_sum + fInput->ffTY[k]; } if (fInput->ffBY[k] > SF_STOPXL_Low) { ffBY_SF_sum = ffBY_SF_sum + fInput->ffBY[k]; } } } if (fInput->ffBX[i] > 0.) { if ((-40 < fInput->fiGammaMysec[0]) && (fInput->fiGammaMysec[0] < 100) && fInput->ffGammaKev[0] > 0.) fAL_Ga->Fill(fInput->ffBX[i], fInput->ffGammaKev[0]); if ((-40 < fInput->fiGammaMysec[1]) && (fInput->fiGammaMysec[1] < 100) && fInput->ffGammaKev[1] > 0.) fAL_Ga->Fill(fInput->ffBX[i], fInput->ffGammaKev[1]); fhEnergyAnswers_BX_pa->Fill(fInput->ffBX[i]); if (fInput->fiMicro == 0 && fInput->ffBX[i] > AL1_STOPXL_Low && fInput->ffBX[i] < AL1_STOPXL_High) { fprintf(Fissions, "pUnpackEvent->fievent_number==%d ||\n", fInput->fievcount); fprintf(ER_AL_stop, "((l_feb_id==0 || l_feb_id==2) && l_cha_id==%d && pUnpackEvent->fievent_number==%d) ||\n", i, fInput->fievcount); } Int_t one = 2; for (k = 0; k < 48; k++) { if (one != 1 && fInput->ffBY[k] > 5000. && fInput->ffBY[k] < (AL1_STOPXL_High + (500.)) && fInput->ffBX[i] > AL1_STOPXL_Low && fInput->ffBX[i] < AL1_STOPXL_High) { fhFPD_fr_ba_fpga->Fill(i, k); { fprintf(ER_AL_sb, "Bot %10.8f\t %d\t %6.2f\t %d\t %6.2f\t %d\t %d\t %6.2f\t %d\t %6.2f\t %d\n", fInput->fiSystemSec * 0.00000001, i, fInput->ffBX[i], k, fInput->ffBY[k], fInput->fiMicro, fInput->fievcount, fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1]); one = 1; } } } } } if (fInput->ffTX[i] > 0. && fInput->fiAnswers_TX_Sig == 1) fhEnergyAnswers_TX_single->Fill(fInput->ffTX[i]); if (fInput->ffBX[i] > 0. && fInput->fiAnswers_BX_Sig == 1) fhEnergyAnswers_BX_single->Fill(fInput->ffBX[i]); if (fInput->ffTX[i] > AL1_STOPXL_Low && fInput->ffTX[i] < AL1_STOPXL_High) { if (fInput->fiMicro == 0 && i < 8) fhdUP1fpga->Fill(1); if (fInput->fiMicro == 0 && i > 7 && i < 12) fhdUP1fpga->Fill(2); if (fInput->fiMicro == 0 && i > 11 && i < 16) fhdUP1fpga->Fill(3); if (fInput->fiMicro == 0 && i > 15 && i < 18) fhdUP1fpga->Fill(4); if (fInput->fiMicro == 0 && i > 17 && i < 20) fhdUP1fpga->Fill(5); if (fInput->fiMicro == 0 && i > 19 && i < 22) fhdUP1fpga->Fill(6); if (fInput->fiMicro == 0 && i > 21 && i < 24) fhdUP1fpga->Fill(7); if (fInput->fiMicro == 0 && i == 24) fhdUP1fpga->Fill(8); if (fInput->fiMicro == 0 && i == 25) fhdUP1fpga->Fill(9); if (fInput->fiMicro == 0 && i == 26) fhdUP1fpga->Fill(10); if (fInput->fiMicro == 0 && i == 27) fhdUP1fpga->Fill(11); if (fInput->fiMicro == 0 && i == 28) fhdUP1fpga->Fill(12); if (fInput->fiMicro == 0 && i == 29) fhdUP1fpga->Fill(13); if (fInput->fiMicro == 0 && i == 30) fhdUP1fpga->Fill(14); if (fInput->fiMicro == 0 && i == 31) fhdUP1fpga->Fill(15); if (fInput->fiMicro == 0) fhdUP2fpga->Fill(i); // cout<<"(pUnpackEvent->fievent_number== "<fievcount<<" && "<< " ((l_feb_id*16+l_cha_id>="<="<fievent_number== "<fievcount<<" ||"<fiSystemSec*0.00000001, // fInput->ffTX[i], // i, // fInput->fiMicro, // fInput->fievcount, // fInput->ffGammaKev[0], // fInput->fiGammaMysec[0], // fInput->ffGammaKev[1], // fInput->fiGammaMysec[1]); // printf("--- Com1_UP \t %6.6f\t %6.2f\t %d\t %d\t %d\t %6.2f\t %d\t %6.2f\t %d\n", // (Float_t) (fInput->fiSystemSec*0.00000001), // fInput->ffTX[i], // i, // fInput->fiMicro, // fInput->fievcount, // fInput->ffGammaKev[0], // fInput->fiGammaMysec[0], // fInput->ffGammaKev[1], // fInput->fiGammaMysec[1]); // } /* if(fInput->fiMicro==0 && fInput->ffCom2_up[i]>AL1_STOPXL_Low && fInput->ffCom2_up[i]fiSystemSec*0.00000001), fInput->ffCom2_up[i], i, fInput->fiMicro, fInput->fievcount, fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1]); printf("--- Com2_UP \t %6.6f\t %6.2f\t %d\t %d\t %d\t %6.2f\t %d\t %6.2f\t %d\n", (Float_t) (fInput->fiSystemSec*0.00000001), fInput->ffCom2_up[i], i, fInput->fiMicro, fInput->fievcount, fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1]); } if(fInput->ffCom3_up[i]>AL1_STOPXL_Low && fInput->ffCom3_up[i]fiSystemSec*0.00000001), fInput->ffCom3_up[i], i, fInput->fiMicro, fInput->fievcount, fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1]); printf("--- Com3_UP \t %6.6f\t %6.2f\t %d\t %d\t %d\t %6.2f\t %d\t %6.2f\t %d\n", (Float_t) (fInput->fiSystemSec*0.00000001), fInput->ffCom3_up[i], i, fInput->fiMicro, fInput->fievcount, fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1]); } */ } if (ffTY_SF_sum > SF_STOPXL_Low && ffBY_SF_sum > SF_STOPXL_Low) { fCom1_Com2->Fill(ffTY_SF_sum, ffBY_SF_sum); fprintf(ER_FI, " Fission : %d\t %lld\t %6f\t %6f\t %d\t %6.1f\t %d\t %6.1f\t %d\n", fInput->fievcount, fInput->fiSystemSec, ffTY_SF_sum, ffBY_SF_sum, fInput->fiMicro, fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1]); } for (i = 0; i < 32; i++) { if (fInput->ffBX[i] > AL1_STOPXL_Low && fInput->ffBX[i] < AL1_STOPXL_High) { if (fInput->fiMicro == 0 && i < 8) fhdDO1fpga->Fill(1); if (fInput->fiMicro == 0 && i > 7 && i < 12) fhdDO1fpga->Fill(2); if (fInput->fiMicro == 0 && i > 11 && i < 16) fhdDO1fpga->Fill(3); if (fInput->fiMicro == 0 && i > 15 && i < 18) fhdDO1fpga->Fill(4); if (fInput->fiMicro == 0 && i > 17 && i < 20) fhdDO1fpga->Fill(5); if (fInput->fiMicro == 0 && i > 19 && i < 22) fhdDO1fpga->Fill(6); if (fInput->fiMicro == 0 && i > 21 && i < 24) fhdDO1fpga->Fill(7); if (fInput->fiMicro == 0 && i == 24) fhdDO1fpga->Fill(8); if (fInput->fiMicro == 0 && i == 25) fhdDO1fpga->Fill(9); if (fInput->fiMicro == 0 && i == 26) fhdDO1fpga->Fill(10); if (fInput->fiMicro == 0 && i == 27) fhdDO1fpga->Fill(11); if (fInput->fiMicro == 0 && i == 28) fhdDO1fpga->Fill(12); if (fInput->fiMicro == 0 && i == 29) fhdDO1fpga->Fill(13); if (fInput->fiMicro == 0 && i == 30) fhdDO1fpga->Fill(14); if (fInput->fiMicro == 0 && i == 31) fhdDO1fpga->Fill(15); if (fInput->fiMicro == 0) fhdDO2fpga->Fill(i); // cout<<"(pUnpackEvent->fievent_number== "<fievcount<<" && "<< " (l_feb_id*16+l_cha_id=="<fievent_number== "<fievcount<<" && "<< " ((l_feb_id*16+l_cha_id>="<="<fievent_number== "<fievcount<<" ||"<fiSystemSec*0.00000001, // fInput->ffBX[i], // i, // fInput->fiMicro, // fInput->fievcount, // fInput->ffGammaKev[0], // fInput->fiGammaMysec[0], // fInput->ffGammaKev[1], // fInput->fiGammaMysec[1]); // // printf("--- Com1_DO \t %6.6f\t %6.2f\t %d\t %d\t %d\t %6.2f\t %d\t %6.2f\t %d\n", // (Float_t) (fInput->fiSystemSec*0.00000001), // fInput->ffBX[i], // i, // fInput->fiMicro, // fInput->fievcount, // fInput->ffGammaKev[0], // fInput->fiGammaMysec[0], // fInput->ffGammaKev[1], // fInput->fiGammaMysec[1]); } /* if(fInput->fiMicro==0 && fInput->ffCom2_do[i]>AL1_STOPXL_Low && fInput->ffCom2_do[i]fiSystemSec*0.00000001), fInput->ffCom2_do[i], i, fInput->fiMicro, fInput->fievcount, fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1]); printf("--- Com2_DO \t %6.6f\t %6.2f\t %d\t %d\t %d\t %6.2f\t %d\t %6.2f\t %d\n", (Float_t) (fInput->fiSystemSec*0.00000001), fInput->ffCom2_do[i], i, fInput->fiMicro, fInput->fievcount, fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1]); } if(fInput->fiMicro==0 && fInput->ffCom3_do[i]>AL1_STOPXL_Low && fInput->ffCom3_do[i]fiSystemSec*0.00000001), fInput->ffCom3_do[i], i, fInput->fiMicro, fInput->fievcount, fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1]); printf("--- Com3_DO \t %6.6f\t %6.2f\t %d\t %d\t %d\t %6.2f\t %d\t %6.2f\t %d\n", (Float_t) (fInput->fiSystemSec*0.00000001), fInput->ffCom3_do[i], i, fInput->fiMicro, fInput->fievcount, fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1]); } */ } for (i = 0; i < 32; i++) { // if(fInput->ffTX[i]>0. && (Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_TX_time[i])*0.00000001fiSystemSec-poutevt->fia1_Answers_TX_time[i])*0.00000001>0.) // fhdAfterFi_Si->Fill(fInput->ffTX[i]); // if(fInput->ffGammaKev[0]>0.) fhdAfterFi_Ga->Fill(fInput->ffGammaKev[0]); // if(fInput->ffGammaKev[1]>0.) fhdAfterFi_Ga->Fill(fInput->ffGammaKev[1]); // // } /*if(){ if(( fInput->ffTX[i]>0. ) || (fInput->ffBX[i]>0. )) { fhAfterFi_Ga->Fill(fInput->ffGammaKev[0]); fhAfterFi_Ga->Fill(fInput->ffGammaKev[1]); fprintf(AL_AL_AL,"EL_Ga \t %6.6f\t %6.2f\t %6.2f\t %d\t %d\t %d\t %6.2f\t %d\t %6.2f\t %d\t %6.6f\t", (Float_t) (fInput->fiSystemSec*0.00000001), fInput->ffTX[i], fInput->ffBX[i], i, fInput->fiMicro, fInput->fievcount, fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1], (Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_TX_time[i])*0.00000001); // fprintf(AL_AL_AL,"Fi \t %6.6f\t %6.2f\t %6.2f\t %d\t %6.2f\t %d\t %6.2f\t %d\n", // (Float_t) (poutevt->fia1_Answers_TX_time[i]*0.00000001), // poutevt->ffa1_Answers_TX[i], // poutevt->ffa1_Answers_BX[i], // i, // poutevt->ffa1_Answers_TX_GammaKev0[i], // poutevt->fia1_Answers_TX_GammaMysec0[i], // poutevt->ffa1_Answers_TX_GammaKev1[i], // poutevt->fia1_Answers_TX_GammaMysec1[i]); } }*/ /* if((Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_TX_time[i])*0.00000001fiSystemSec-poutevt->fia1_Answers_TX_time[i])*0.00000001>10. && ((fInput->ffTX[i]>0. && fInput->ffTX[i]<3000.) ) && fInput->ffGammaKev[0]>0. && fInput->fiMicro==0) { fhAfterFi_Si->Fill(fInput->ffGammaKev[0]); fprintf(AL_AL_AL,"Fi \t %6.6f\t %6.2f\t %6.2f\t %d\t %6.2f\t %d\t %6.2f\t %d\t", (Float_t) (poutevt->fia1_Answers_TX_time[i]*0.00000001), poutevt->ffa1_Answers_TX[i], poutevt->ffa1_Answers_BX[i], i, poutevt->ffa1_Answers_TX_GammaKev0[i], poutevt->fia1_Answers_TX_GammaMysec0[i], poutevt->ffa1_Answers_TX_GammaKev1[i], poutevt->fia1_Answers_TX_GammaMysec1[i]); fprintf(AL_AL_AL,"EL_Ga \t %6.6f\t %6.2f\t %6.2f\t %d\t %d\t %d\t %6.2f\t %d\t %6.2f\t %d\t %6.6f\n", (Float_t) (fInput->fiSystemSec*0.00000001), fInput->ffTX[i], fInput->ffBX[i], i, fInput->fiMicro, fInput->fievcount, fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1], (Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_TX_time[i])*0.00000001); } if((Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_TX_time[i])*0.00000001fiSystemSec-poutevt->fia1_Answers_TX_time[i])*0.00000001>10. && ((fInput->ffTX[i]>0. && fInput->ffTX[i]<3000.) ) && fInput->ffGammaKev[1]>0. && fInput->fiMicro==0) { fhAfterFi_Si->Fill(fInput->ffGammaKev[1]); } if((Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_BX_time[i])*0.00000001fiSystemSec-poutevt->fia1_Answers_BX_time[i])*0.00000001>10. && ((fInput->ffTX[i]>0. && fInput->ffTX[i]<3000.) ) && fInput->ffGammaKev[0]>0. && fInput->fiMicro==0) { fhAfterFi_Si->Fill(fInput->ffGammaKev[0]); } if((Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_BX_time[i])*0.00000001fiSystemSec-poutevt->fia1_Answers_BX_time[i])*0.00000001>10. && ((fInput->ffTX[i]>0. && fInput->ffTX[i]<3000.) ) && fInput->ffGammaKev[1]>0. && fInput->fiMicro==0) { fhAfterFi_Si->Fill(fInput->ffGammaKev[1]); } */ // if((Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_TX_time[i])*0.00000001fiSystemSec-poutevt->fia1_Answers_TX_time[i])*0.00000001>0. && fInput->ffBX[i]>0. && fInput->ffGammaKev[0]>0. && fInput->fiMicro==0) // { // fhAfterFi_Si->Fill(fInput->ffGammaKev[0]); // } // // if((Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_TX_time[i])*0.00000001fiSystemSec-poutevt->fia1_Answers_TX_time[i])*0.00000001>0. && fInput->ffBX[i]>0. && fInput->ffGammaKev[1]>0. && fInput->fiMicro==0) // { // fhAfterFi_Si->Fill(fInput->ffGammaKev[1]); // // } /* //start corr analysis //correlation for first compact detector if(fInput->ffTX[i]>AL2_STOPXL_Low && fInput->ffTX[i]-1. && ((Double_t)(i)-(Double_t)(k)+1.)<32.){ // printf(" qui"); // printf("fiSystemSec: %lld\t fia1_Answers_TX_time: %10.6f \n",fInput->fiSystemSec,(Float_t)(poutevt->fia1_Answers_TX_time[i-k+1])); if((Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_TX_time[i-k+1])*0.00000001fiSystemSec-poutevt->fia1_Answers_TX_time[i-k+1])*0.00000001>0.){ Answers_TX_TX++; printf(" ----- Answers_TX_TX : Num=%d ----- \n",Answers_TX_TX); printf(" AL1 Answers_TX : %d\t %d\t %lld\t %6.1f\t %d\t %d\t %6.1f\t %d\t %6.1f\t %d\n", poutevt->fia1_Answers_TX_evcount[i-k+1], i-k+1, poutevt->fia1_Answers_TX_time[i-k+1], poutevt->ffa1_Answers_TX[i-k+1], poutevt->fia1_Answers_TX_Micro[i-k+1], poutevt->fia1_Answers_TX_multi[i-k+1], poutevt->ffa1_Answers_TX_GammaKev0[i-k+1], poutevt->fia1_Answers_TX_GammaMysec0[i-k+1], poutevt->ffa1_Answers_TX_GammaKev1[i-k+1], poutevt->fia1_Answers_TX_GammaMysec1[i-k+1]); printf(" AL2 Answers_TX : %d\t %d\t %lld\t %6.1f\t %d\t %d\t %6.1f\t %d\t %6.1f\t %d\t %10.6f\n\n", fInput->fievcount, i, fInput->fiSystemSec, fInput->ffTX[i], fInput->fiMicro, fInput->fiTX_hit[i], fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1], (Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_TX_time[i-k+1])*0.00000001); fprintf(ER_AL_stop," AL1 Answers_TX : %d\t %d\t %lld\t %6.1f\t %d\t %d\t %6.1f\t %d\t %6.1f\t %d\t", poutevt->fia1_Answers_TX_evcount[i-k+1], i-k+1, poutevt->fia1_Answers_TX_time[i-k+1], poutevt->ffa1_Answers_TX[i-k+1], poutevt->fia1_Answers_TX_Micro[i-k+1], poutevt->fia1_Answers_TX_multi[i-k+1], poutevt->ffa1_Answers_TX_GammaKev0[i-k+1], poutevt->fia1_Answers_TX_GammaMysec0[i-k+1], poutevt->ffa1_Answers_TX_GammaKev1[i-k+1], poutevt->fia1_Answers_TX_GammaMysec1[i-k+1]); fprintf(ER_AL_stop," AL2 Answers_TX : %d\t %d\t %lld\t %6.1f\t %d\t %d\t %6.1f\t %d\t %6.1f\t %d\t %10.6f\n", fInput->fievcount, i, fInput->fiSystemSec, fInput->ffTX[i], fInput->fiMicro, fInput->fiTX_hit[i], fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1], (Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_TX_time[i-k+1])*0.00000001); } // printf(" TEST 1 b \n"); } } // printf(" TEST 1 c \n"); for(k=1;k<2;k++){ // printf(" i : %d\t k : %d \t index : %d \n",i, k , i-k+1); if(((Double_t)(i)-(Double_t)(k)+1.)>-1. && ((Double_t)(i)-(Double_t)(k)+1.)<32.){ // printf(" qui"); // printf(" fiSystemSec : %lld\t fia1_Answers_BX_time : %10.6f \n",fInput->fiSystemSec,(Float_t) // (poutevt->fia1_Answers_BX_time[i-k+1])); if((Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_BX_time[i-k+1])*0.00000001fiSystemSec-poutevt->fia1_Answers_BX_time[i-k+1])*0.00000001>0.) { // printf(" TEST 1 d\n"); Answers_BX_TX++; printf(" ----- Answers_BX_TX : Num=%d ----- \n",Answers_BX_TX); printf(" AL1 Answers_BX : %d\t %d\t %lld\t %6.1f\t %d\t %d\t %6.1f\t %d\t %6.1f\t %d\n", poutevt->fia1_Answers_BX_evcount[i-k+1], i-k+1, poutevt->fia1_Answers_BX_time[i-k+1], poutevt->ffa1_Answers_BX[i-k+1], poutevt->fia1_Answers_BX_Micro[i-k+1], poutevt->fia1_Answers_BX_multi[i-k+1], poutevt->ffa1_Answers_BX_GammaKev0[i-k+1], poutevt->fia1_Answers_BX_GammaMysec0[i-k+1], poutevt->ffa1_Answers_BX_GammaKev1[i-k+1], poutevt->fia1_Answers_BX_GammaMysec1[i-k+1]); printf(" AL2 Answers_TX : %d\t %d\t %lld\t %6.1f\t %d\t %d\t %6.1f\t %d\t %6.1f\t %d\t %10.6f\n\n", fInput->fievcount, i, fInput->fiSystemSec, fInput->ffTX[i], fInput->fiMicro, fInput->fiTX_hit[i], fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1], (Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_BX_time[i-k+1])*0.00000001); fprintf(ER_AL_stop," AL1 Answers_BX : %d\t %d\t %lld\t %6.1f\t %d\t %d\t %6.1f\t %d\t %6.1f\t %d\t", poutevt->fia1_Answers_BX_evcount[i-k+1], i-k+1, poutevt->fia1_Answers_BX_time[i-k+1], poutevt->ffa1_Answers_BX[i-k+1], poutevt->fia1_Answers_BX_Micro[i-k+1], poutevt->fia1_Answers_BX_multi[i-k+1], poutevt->ffa1_Answers_BX_GammaKev0[i-k+1], poutevt->fia1_Answers_BX_GammaMysec0[i-k+1], poutevt->ffa1_Answers_BX_GammaKev1[i-k+1], poutevt->fia1_Answers_BX_GammaMysec1[i-k+1]); fprintf(ER_AL_stop," AL2 Answers_TX : %d\t %d\t %lld\t %6.1f\t %d\t %d\t %6.1f\t %d\t %6.1f\t %d\t %10.6f\n", fInput->fievcount, i, fInput->fiSystemSec, fInput->ffTX[i], fInput->fiMicro, fInput->fiTX_hit[i], fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1], (Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_BX_time[i-k+1])*0.00000001); } // printf(" TEST 1 e\n"); } } // printf(" TEST 1 f \n"); //printf("%f\n",AL2_Time); } // printf(" TEST 2 \n"); if(fInput->ffBX[i]>AL2_STOPXL_Low && fInput->ffBX[i]-1. && ((Double_t)(i)-(Double_t)(k)+1.)<32.){ // printf(" qui"); // printf(" fiSystemSec : %lld\t fia1_Answers_BX_time : %10.6f \n",fInput->fiSystemSec,(Float_t) // (poutevt->fia1_Answers_BX_time[i-k+1])); if((Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_TX_time[i-k+1])*0.00000001fiSystemSec-poutevt->fia1_Answers_TX_time[i-k+1])*0.00000001>0.) { Answers_TX_BX++; printf(" ----- Answers_TX_BX : Num=%d ----- \n",Answers_TX_BX); printf(" AL1 Answers_TX : %d\t %d\t %lld\t %6.1f\t %d\t %d\t %6.1f\t %d\t %6.1f\t %d\n", poutevt->fia1_Answers_TX_evcount[i-k+1], i-k+1, poutevt->fia1_Answers_TX_time[i-k+1], poutevt->ffa1_Answers_TX[i-k+1], poutevt->fia1_Answers_TX_Micro[i-k+1], poutevt->fia1_Answers_TX_multi[i-k+1], poutevt->ffa1_Answers_TX_GammaKev0[i-k+1], poutevt->fia1_Answers_TX_GammaMysec0[i-k+1], poutevt->ffa1_Answers_TX_GammaKev1[i-k+1], poutevt->fia1_Answers_TX_GammaMysec1[i-k+1]); printf(" AL2 Answers_BX : %d\t %d\t %lld\t %6.1f\t %d\t %d\t %6.1f\t %d\t %6.1f\t %d\t %10.6f\n\n", fInput->fievcount, i, fInput->fiSystemSec, fInput->ffBX[i], fInput->fiMicro, fInput->fiBX_hit[i], fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1], (Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_TX_time[i-k+1])*0.00000001); fprintf(ER_AL_stop," AL1 Answers_TX : %d\t %d\t %lld\t %6.1f\t %d\t %d\t %6.1f\t %d\t %6.1f\t %d\t", poutevt->fia1_Answers_TX_evcount[i-k+1], i-k+1, poutevt->fia1_Answers_TX_time[i-k+1], poutevt->ffa1_Answers_TX[i-k+1], poutevt->fia1_Answers_TX_Micro[i-k+1], poutevt->fia1_Answers_TX_multi[i-k+1], poutevt->ffa1_Answers_TX_GammaKev0[i-k+1], poutevt->fia1_Answers_TX_GammaMysec0[i-k+1], poutevt->ffa1_Answers_TX_GammaKev1[i-k+1], poutevt->fia1_Answers_TX_GammaMysec1[i-k+1]); fprintf(ER_AL_stop," AL2 Answers_BX : %d\t %d\t %lld\t %6.1f\t %d\t %d\t %6.1f\t %d\t %6.1f\t %d\t %10.6f\n", fInput->fievcount, i, fInput->fiSystemSec, fInput->ffBX[i], fInput->fiMicro, fInput->fiBX_hit[i], fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1], (Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_TX_time[i-k+1])*0.00000001); } // printf(" TEST 2 b \n"); } } // printf(" TEST 2 c \n"); for(k=1;k<2;k++){ // printf(" i : %d\t k : %d \t index : %d \n",i, k , i-k+1); if(((Double_t)(i)-(Double_t)(k)+1.)>-1. && ((Double_t)(i)-(Double_t)(k)+1.)<32.){ // printf(" qui"); // printf(" fiSystemSec : %lld\t fia1_Answers_BX_time : %10.6f \n",fInput->fiSystemSec,(Float_t) // (poutevt->fia1_Answers_BX_time[i-k+1])); if((Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_BX_time[i-k+1])*0.00000001fiSystemSec-poutevt->fia1_Answers_BX_time[i-k+1])*0.00000001>0.) { // printf(" TEST 2 d\n"); Answers_BX_BX++; printf(" ----- Answers_BX_BX : Num=%d ----- \n",Answers_BX_BX); printf(" AL1 Answers_BX : %d\t %d\t %lld\t %6.1f\t %d\t %d\t %6.1f\t %d\t %6.1f\t %d\n", poutevt->fia1_Answers_BX_evcount[i-k+1], i-k+1, poutevt->fia1_Answers_BX_time[i-k+1], poutevt->ffa1_Answers_BX[i-k+1], poutevt->fia1_Answers_BX_Micro[i-k+1], poutevt->fia1_Answers_BX_multi[i-k+1], poutevt->ffa1_Answers_BX_GammaKev0[i-k+1], poutevt->fia1_Answers_BX_GammaMysec0[i-k+1], poutevt->ffa1_Answers_BX_GammaKev1[i-k+1], poutevt->fia1_Answers_BX_GammaMysec1[i-k+1]); printf(" AL2 Answers_BX : %d\t %d\t %lld\t %6.1f\t %d\t %d\t %6.1f\t %d\t %6.1f\t %d\t %10.6f\n\n", fInput->fievcount, i, fInput->fiSystemSec, fInput->ffBX[i], fInput->fiMicro, fInput->fiBX_hit[i], fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1], (Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_BX_time[i-k+1])*0.00000001); fprintf(ER_AL_stop," AL1 Answers_BX : %d\t %d\t %lld\t %6.1f\t %d\t %d\t %6.1f\t %d\t %6.1f\t %d\t", poutevt->fia1_Answers_BX_evcount[i-k+1], i-k+1, poutevt->fia1_Answers_BX_time[i-k+1], poutevt->ffa1_Answers_BX[i-k+1], poutevt->fia1_Answers_BX_Micro[i-k+1], poutevt->fia1_Answers_BX_multi[i-k+1], poutevt->ffa1_Answers_BX_GammaKev0[i-k+1], poutevt->fia1_Answers_BX_GammaMysec0[i-k+1], poutevt->ffa1_Answers_BX_GammaKev1[i-k+1], poutevt->fia1_Answers_BX_GammaMysec1[i-k+1]); fprintf(ER_AL_stop," AL2 Answers_BX : %d\t %d\t %lld\t %6.1f\t %d\t %d\t %6.1f\t %d\t %6.1f\t %d\t %10.6f\n", fInput->fievcount, i, fInput->fiSystemSec, fInput->ffBX[i], fInput->fiMicro, fInput->fiBX_hit[i], fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1], (Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_BX_time[i-k+1])*0.00000001); } // printf(" TEST 2 e\n"); } } // printf(" TEST 2 f\n"); //printf("%f\n",AL2_Time); } //end correlation for first compact if(fInput->ffTX[i]>SF_STOPXL_Low) { for(k=1;k<2;k++){ if(((Double_t)(i)-(Double_t)(k)+1.)>-2. && ((Double_t)(i)-(Double_t)(k)+1.)<32.){ if((Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_TX_time[i-k+1])*0.00000001fiSystemSec-poutevt->fia1_Answers_TX_time[i-k+1])*0.00000001>0.) { // printf(" ----- Answers_TX_TX ----- \n"); printf(" AL Answers_TX : %d\t %d\t %lld\t %6f\t %6.1f\t %d\t %6.1f\t %d\n", poutevt->fia1_Answers_TX_evcount[i-k+1], i-k+1, poutevt->fia1_Answers_TX_time[i-k+1], poutevt->ffa1_Answers_TX[i-k+1], poutevt->ffa1_Answers_TX_GammaKev0[i-k+1], poutevt->fia1_Answers_TX_GammaMysec0[i-k+1], poutevt->ffa1_Answers_TX_GammaKev1[i-k+1], poutevt->fia1_Answers_TX_GammaMysec1[i-k+1]); printf(" SF Answers_TX : %d\t %d\t %lld\t %6f\t %6.1f\t %d\t %6.1f\t %d\n\n", fInput->fievcount, i, fInput->fiSystemSec, fInput->ffTX[i], fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1]); fprintf(ER_AL_stop," AL Answers_TX : %d\t %d\t %lld\t %6f\t %6.1f\t %d\t %6.1f\t %d\t", poutevt->fia1_Answers_TX_evcount[i-k+1], i-k+1, poutevt->fia1_Answers_TX_time[i-k+1], poutevt->ffa1_Answers_TX[i-k+1], poutevt->ffa1_Answers_TX_GammaKev0[i-k+1], poutevt->fia1_Answers_TX_GammaMysec0[i-k+1], poutevt->ffa1_Answers_TX_GammaKev1[i-k+1], poutevt->fia1_Answers_TX_GammaMysec1[i-k+1]); fprintf(ER_AL_stop," SF Answers_TX : %d\t %d\t %lld\t %6f\t %6.1f\t %d\t %6.1f\t %d\n", fInput->fievcount, i, fInput->fiSystemSec, fInput->ffTX[i], fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1]); } } } for(k=1;k<2;k++){ if(((Double_t)(i)-(Double_t)(k)+1.)>-2. && ((Double_t)(i)-(Double_t)(k)+1.)<32.){ if((Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_BX_time[i-k+1])*0.00000001fiSystemSec-poutevt->fia1_Answers_BX_time[i-k+1])*0.00000001>0.) { // printf(" ----- Answers_BX_TX ----- \n"); printf(" AL Answers_BX : %d\t %d\t %lld\t %6f\t %6.1f\t %d\t %6.1f\t %d\n", poutevt->fia1_Answers_BX_evcount[i-k+1], i-k+1, poutevt->fia1_Answers_BX_time[i-k+1], poutevt->ffa1_Answers_BX[i-k+1], poutevt->ffa1_Answers_BX_GammaKev0[i-k+1], poutevt->fia1_Answers_BX_GammaMysec0[i-k+1], poutevt->ffa1_Answers_BX_GammaKev1[i-k+1], poutevt->fia1_Answers_BX_GammaMysec1[i-k+1]); printf(" SF Answers_TX : %d\t %d\t %lld\t %6f\t %6.1f\t %d\t %6.1f\t %d\n\n", fInput->fievcount, i, fInput->fiSystemSec, fInput->ffTX[i], fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1]); fprintf(ER_AL_stop," AL Answers_BX : %d\t %d\t %lld\t %6f\t %6.1f\t %d\t %6.1f\t %d\t", poutevt->fia1_Answers_BX_evcount[i-k+1], i-k+1, poutevt->fia1_Answers_BX_time[i-k+1], poutevt->ffa1_Answers_BX[i-k+1], poutevt->ffa1_Answers_BX_GammaKev0[i-k+1], poutevt->fia1_Answers_BX_GammaMysec0[i-k+1], poutevt->ffa1_Answers_BX_GammaKev1[i-k+1], poutevt->fia1_Answers_BX_GammaMysec1[i-k+1]); fprintf(ER_AL_stop," SF Answers_TX : %d\t %d\t %lld\t %6f\t %6.1f\t %d\t %6.1f\t %d\n", fInput->fievcount, i, fInput->fiSystemSec, fInput->ffTX[i], fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1]); } } } //printf("%f\n",AL2_Time); } if(fInput->ffBX[i]>SF_STOPXL_Low){ for(k=1;k<2;k++){ if(((Double_t)(i)-(Double_t)(k)+1.)>-2. && ((Double_t)(i)-(Double_t)(k)+1.)<32.){ if((Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_TX_time[i-k+1])*0.00000001fiSystemSec-poutevt->fia1_Answers_TX_time[i-k+1])*0.00000001==0) { // printf(" ----- Answers_TX_BX ----- \n"); printf(" AL Answers_TX : %d\t %d\t %lld\t %6f\t %6.1f\t %d\t %6.1f\t %d\n", poutevt->fia1_Answers_TX_evcount[i-k+1], i-k+1, poutevt->fia1_Answers_TX_time[i-k+1], poutevt->ffa1_Answers_TX[i-k+1], poutevt->ffa1_Answers_TX_GammaKev0[i-k+1], poutevt->fia1_Answers_TX_GammaMysec0[i-k+1], poutevt->ffa1_Answers_TX_GammaKev1[i-k+1], poutevt->fia1_Answers_TX_GammaMysec1[i-k+1]); printf(" SF Answers_BX : %d\t %d\t %lld\t %6f\t %6.1f\t %d\t %6.1f\t %d\n\n", fInput->fievcount, i, fInput->fiSystemSec, fInput->ffBX[i], fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1]); fprintf(ER_AL_stop," AL Answers_TX : %d\t %d\t %lld\t %6f\t %6.1f\t %d\t %6.1f\t %d\t", poutevt->fia1_Answers_TX_evcount[i-k+1], i-k+1, poutevt->fia1_Answers_TX_time[i-k+1], poutevt->ffa1_Answers_TX[i-k+1], poutevt->ffa1_Answers_TX_GammaKev0[i-k+1], poutevt->fia1_Answers_TX_GammaMysec0[i-k+1], poutevt->ffa1_Answers_TX_GammaKev1[i-k+1], poutevt->fia1_Answers_TX_GammaMysec1[i-k+1]); fprintf(ER_AL_stop," SF Answers_BX : %d\t %d\t %lld\t %6f\t %6.1f\t %d\t %6.1f\t %d\n", fInput->fievcount, i, fInput->fiSystemSec, fInput->ffBX[i], fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1]); } } } for(k=1;k<2;k++){ if(((Double_t)(i)-(Double_t)(k)+1.)>-2. && ((Double_t)(i)-(Double_t)(k)+1.)<32.){ if((Float_t)(fInput->fiSystemSec-poutevt->fia1_Answers_BX_time[i-k+1])*0.00000001fiSystemSec-poutevt->fia1_Answers_BX_time[i-k+1])*0.00000001==0) { // printf(" ----- Answers_BX_BX ----- \n"); printf(" AL Answers_BX : %d\t %d\t %lld\t %6f\t %6.1f\t %d\t %6.1f\t %d\n", poutevt->fia1_Answers_BX_evcount[i-k+1], i-k+1, poutevt->fia1_Answers_BX_time[i-k+1], poutevt->ffa1_Answers_BX[i-k+1], poutevt->ffa1_Answers_BX_GammaKev0[i-k+1], poutevt->fia1_Answers_BX_GammaMysec0[i-k+1], poutevt->ffa1_Answers_BX_GammaKev1[i-k+1], poutevt->fia1_Answers_BX_GammaMysec1[i-k+1]); printf(" SF Answers_BX : %d\t %d\t %lld\t %6f\t %6.1f\t %d\t %6.1f\t %d\n\n", fInput->fievcount, i, fInput->fiSystemSec, fInput->ffBX[i], fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1]); fprintf(ER_AL_stop," AL Answers_BX : %d\t %d\t %lld\t %6f\t %6.1f\t %d\t %6.1f\t %d\t", poutevt->fia1_Answers_BX_evcount[i-k+1], i-k+1, poutevt->fia1_Answers_BX_time[i-k+1], poutevt->ffa1_Answers_BX[i-k+1], poutevt->ffa1_Answers_BX_GammaKev0[i-k+1], poutevt->fia1_Answers_BX_GammaMysec0[i-k+1], poutevt->ffa1_Answers_BX_GammaKev1[i-k+1], poutevt->fia1_Answers_BX_GammaMysec1[i-k+1]); fprintf(ER_AL_stop," SF Answers_BX : %d\t %d\t %lld\t %6f\t %6.1f\t %d\t %6.1f\t %d\n", fInput->fievcount, i, fInput->fiSystemSec, fInput->ffBX[i], fInput->ffGammaKev[0], fInput->fiGammaMysec[0], fInput->ffGammaKev[1], fInput->fiGammaMysec[1]); } } } //printf("%f\n",AL2_Time); } //end correlation for first compact //start to sort the first alpha if(fInput->ffTX[i]>AL3_STOPXL_Low && fInput->ffTX[i]fia1_Answers_TX_evcount[i]=fInput->fievcount; poutevt->fia1_Answers_TX_time[i]=fInput->fiSystemSec; poutevt->ffa1_Answers_TX[i]=fInput->ffTX[i]; poutevt->fia1_Answers_TX_multi[i]=fInput->fiTX_hit[i]; poutevt->fia1_Answers_TX_Micro[i]=fInput->fiMicro; poutevt->ffa1_Answers_TX_GammaKev0[i]=fInput->ffGammaKev[0]; poutevt->fia1_Answers_TX_GammaMysec0[i]=fInput->fiGammaMysec[0]; poutevt->ffa1_Answers_TX_GammaKev1[i]=fInput->ffGammaKev[1]; poutevt->fia1_Answers_TX_GammaMysec1[i]=fInput->fiGammaMysec[1]; // printf(" ----- found high energy Alpha or SF -----\n "); // printf(" Answers_TX : %d\t %d\t %lld\t %6f\t %d\t %d\t %d \n\n", // fInput->fievcount, // i, // fInput->fiSystemSec, // fInput->ffTX[i], // fInput->fiSatTimeAnswers_TX[i], // fInput->fiTX_hit[i], // fInput->fiMicro); // cout<<"((("<fiSystemSec<<"- fInput->fiSystemSec)>300000000 && ("<fiSystemSec<<"- fInput->fiSystemSec)ffBX[i]=="<ffTX[i]=="<ffBX[i]>AL3_STOPXL_Low && fInput->ffBX[i]fia1_Answers_BX_evcount[i]=fInput->fievcount; poutevt->fia1_Answers_BX_time[i]=fInput->fiSystemSec; poutevt->ffa1_Answers_BX[i]=fInput->ffBX[i]; poutevt->fia1_Answers_BX_multi[i]=fInput->fiBX_hit[i]; poutevt->fia1_Answers_BX_Micro[i]=fInput->fiMicro; poutevt->ffa1_Answers_BX_GammaKev0[i]=fInput->ffGammaKev[0]; poutevt->fia1_Answers_BX_GammaMysec0[i]=fInput->fiGammaMysec[0]; poutevt->ffa1_Answers_BX_GammaKev1[i]=fInput->ffGammaKev[1]; poutevt->fia1_Answers_BX_GammaMysec1[i]=fInput->fiGammaMysec[1]; // printf(" ----- found high energy Alpha or SF -----\n "); // printf(" Answers_BX : %d\t %d\t %lld\t %6f\t %d\t %d\t %d \n\n", // fInput->fievcount, // i, // fInput->fiSystemSec, // fInput->ffBX[i], // fInput->fiSatTimeAnswers_BX[i], // fInput->fiBX_hit[i], // fInput->fiMicro); // cout<<"((("<fiSystemSec<<"- fInput->fiSystemSec)>300000000 && ("<fiSystemSec<<"- fInput->fiSystemSec)ffBX["<0. || fInput->ffTX["<0."<SetValid(kFALSE); // events aInt_t Re not stoInt_t Red until kTRUE is set } // BuildCalEvent